Search for: All records

Abstract Hyporheic zones are commonly regarded as resilient and enduring interfaces between groundwater and surface water in river corridors. In particular, bedform‐induced advective pumping hyporheic exchange (bedform‐induced exchange) is often perceived as a relatively persistent mechanism in natural river systems driving water, solutes, and energy exchanges between the channel and its surrounding streambed sediments. Numerous studies have been based on this presumption. To evaluate the persistence of hyporheic zones under varying hydrologic conditions, we use a multi‐physics framework to model advective pumping bedform‐induced hyporheic exchange in response to a series of seasonal‐ and event‐scale groundwater table fluctuation scenarios, which lead to episodic river‐aquifer disconnections and reconnections. Our results suggest that hyporheic exchange is not as ubiquitous as generally assumed. Instead, the bedform‐induced hyporheic exchange is restricted to a narrow range of conditions characterized by minor river‐groundwater head differences, is intermittent, and can be easily obliterated by minor losing groundwater conditions. These findings shed light on the fragility of bedform‐induced hyporheic exchange and have important implications for biogeochemical transformations along river corridors. 

Abstract Hyporheic exchange in streams is critical to ecosystem functions such as nutrient cycling along river corridors, especially for slowly moving or small stream systems. The transient storage model (TSM) has been widely used for modeling of hyporheic exchange. TSM calibration, for hyporheic exchange, is typically used to estimate four parameters, including the mass exchange rate coefficient, the dispersion coefficient, stream cross‐sectional area, and hyporheic zone cross‐sectional area. Prior studies have raised concerns regarding the non‐uniqueness of the inverse problem for the TSM, that is, the occurrence of different parameter vectors resulting in TSM solution that reproduces the observed in‐stream tracer break through curve (BTC) with the same error. This leads to practical non‐identifiability in determining the unknown parameter vector values even when global‐optimal values exist, and the parameter optimization becomes practically non‐unique. To address this problem, we applied the simulated annealing method to calibrate the TSM to BTCs, because it is less susceptible to local minima‐induced non‐identifiability. A hypothetical (or synthetic) tracer test data set with known parameters was developed to demonstrate the capability of the simulated annealing method to find the global minimum parameter vector, and it identified the “hypothetically‐true” global minimum parameter vector even with input data that were modified with up to 10% noise without increasing the number of iterations required for convergence. The simulated annealing TSM was then calibrated using two in‐stream tracer tests conducted in East Fork Poplar Creek, Tennessee. Simulated annealing was determined to be appropriate for quantifying the TSM parameter vector because of its search capability for the global minimum parameter vector. 

Heterogeneity of soil hydraulic (e.g., hydraulic conductivity (KS), porosity (θS)) and chemical (e.g., solid-phase adsorption (Kd)) properties complicates contaminant transport by creating spatial variability in sources of contaminant leaching. There is a knowledge gap on the effect of the interplay between these properties on the retardation and transport of per- and polyfluoroalkyl substances (PFAS) with different properties including carbon–fluorine chain-length and functional groups even in water-saturated conditions. Breakthrough curves have been used to evaluate PFAS transport behavior through heterogeneous media, including arrival time, maximum concentration, and tailing behavior. Contaminant mass flux reduction and mass removal correlations are also compared using numerical modeling to characterize PFAS transport through different source zones within a two-domain, heterogeneous system with comparison to homogeneous scenarios under water-saturated conditions. With heterogeneous properties, model sensitivity to KS was the highest among the other parameters and was controlled by the KS ratio between the different soils. The PFAS models in the homogeneous and heterogeneous scenarios were both sensitive to θS, depending on PFAS chain length. However, long-chain PFAS were less sensitive to θS variability compared to short-chain PFAS due to their higher Kd. The homogeneous and heterogeneous scenarios were equally sensitive to Kd variability, which was dependent on PFAS chain length. 

This study evaluates the performance of multiple machine learning (ML) algorithms and electrical resistivity (ER) arrays for inversion with comparison to a conventional Gauss-Newton numerical inversion method. Four different ML models and four arrays were used for the estimation of only six variables for locating and characterizing hypothetical subsurface targets. The combination of dipole-dipole with Multilayer Perceptron Neural Network (MLP-NN) had the highest accuracy. Evaluation showed that both MLP-NN and Gauss-Newton methods performed well for estimating the matrix resistivity while target resistivity accuracy was lower, and MLP-NN produced sharper contrast at target boundaries for the field and hypothetical data. Both methods exhibited comparable target characterization performance, whereas MLP-NN had increased accuracy compared to Gauss-Newton in prediction of target width and height, which was attributed to numerical smoothing present in the Gauss-Newton approach. MLP-NN was also applied to a field dataset acquired at U.S. DOE Hanford site. 

Physically based hydrologic models require significant effort and extensive information for development, calibration, and validation. The study explored the use of the random forest regression (RFR), a supervised machine learning (ML) model, as an alternative to the physically based Soil and Water Assessment Tool (SWAT) for predicting streamflow in the Rio Grande Headwaters near Del Norte, a snowmelt-dominated mountainous watershed of the Upper Rio Grande Basin. Remotely sensed data were used for the random forest machine learning analysis (RFML) and RStudio for data processing and synthesizing. The RFML model outperformed the SWAT model in accuracy and demonstrated its capability in predicting streamflow in this region. We implemented a customized approach to the RFR model to assess the model’s performance for three training periods, across 1991–2010, 1996–2010, and 2001–2010; the results indicated that the model’s accuracy improved with longer training periods, implying that the model trained on a more extended period is better able to capture the parameters’ variability and reproduce streamflow data more accurately. The variable importance (i.e., IncNodePurity) measure of the RFML model revealed that the snow depth and the minimum temperature were consistently the top two predictors across all training periods. The paper also evaluated how well the SWAT model performs in reproducing streamflow data of the watershed with a conventional approach. The SWAT model needed more time and data to set up and calibrate, delivering acceptable performance in annual mean streamflow simulation, with satisfactory index of agreement (d), coefficient of determination (R2), and percent bias (PBIAS) values, but monthly simulation warrants further exploration and model adjustments. The study recommends exploring snowmelt runoff hydrologic processes, dust-driven sublimation effects, and more detailed topographic input parameters to update the SWAT snowmelt routine for better monthly flow estimation. The results provide a critical analysis for enhancing streamflow prediction, which is valuable for further research and water resource management, including snowmelt-driven semi-arid regions.