Geostatistical features of streambed vertical hydraulic conductivities in Frenchman Creek Watershed in Western Nebraska

This study evaluated the spatial variability of streambed vertical hydraulic conductivity (K_v) in different stream morphologies in the Frenchman Creek Watershed, Western Nebraska, using different variogram models. Streambed K_v values were determined in situ using permeameter tests at 10 sites in Frenchman, Stinking Water and Spring Creeks during the dry season at baseflow conditions. Measurements were taken both in straight and meandering stream channels during a 5 day period at similar flow conditions. Each test site comprised of at least three transects and each transect comprised of at least three K_v measurements. Linear, Gaussian, exponential and spherical variogram models were used with Kriging gridding method for the 10 sites. As a goodness-of-fit statistic for the variogram models, cross-validation results showed differences in the median absolute deviation and the standard deviation of the cross-validation residuals. Results show that using the geometric means of the 10 sites for gridding performs better than using either all the K_v values from the 93 permeameter tests or 10 K_v values from the middle transects and centre permeameters. Incorporating both the spatial variability and the uncertainty involved in the measurement at a reach segment can yield more accurate grid results that can be useful in calibrating K_v at watershed or sub-watershed scales in distributed hydrological models.     

河床水力传导度及其各向异性的测定

采用直接测定法观测了黑河中游段河床水力传导度及其各向异性,结果表明:河床水力传导度不仅存在较强的各向异性,而且存在空间尺度上的变异性.河床中心位置在垂直、水平和θ=30°方向的平均水力传导度分别为0.45、22.49和1.71 m/d,河床边分别为5.95、29.69、16.80 m/d.在同一测点,水力传导度随着与水平方向的夹角增大呈幂函数曲线下降.试验结果表明河床边是河水的主要渗漏区,并且以侧渗为主.     

河床沉积物非均质性影响下的潜流交换数值模拟

为了揭示河床沉积物非均质性对潜流交换的影响,构建了沙波地形作用下的地表水-地下水耦合模型,通过生成不同的渗透系数随机场,讨论了不同地表水动力过程和河床沉积物非均质性对潜流交换通量、交换空间及平均停留时间的影响规律。结果表明,所构建的地表水-地下水耦合模型能够准确刻画水沙界面附近流场,模型具有良好的适用性;均质或非均质河床沉积物情景下,水沙界面上的平均交换通量、停留时间与雷诺数之间均呈现幂函数关系,潜流交换深度则在地表水进入完全紊流之后趋于稳定。结果还表明,较强的河床沉积物非均质性能够有效增强水沙界面上潜流交换通量和空间交换频率,但会制约潜流交换空间并缩短水流在潜流带中的停留时间。     

Investigation of the thermal regime and subsurface properties of a tidally affected,variably saturated streambed

Temperature and moisture content in the variably saturated subsurface are two of the most important physical parameters that govern a wide variety of geochemical and ecological processes. An understanding of thermal and hydraulic processes and properties of transient vadose zones is therefore fundamental in the evaluation of such processes. Here, an investigation of the thermal regime and subsurface properties of a tidally affected, variably saturated streambed is presented. Field and laboratory measurements, as well as a forward numerical model, are jointly employed in the investigation. Temperature, soil moisture, surface level, and water level data were recorded in a transect perpendicular to a tidally driven stream. Frequency‐domain analysis of the subsurface temperature measurements revealed the rapid decay of the tidal temperature driver within the top ～30 cm of sediment. Several techniques were used to evaluate subsurface thermal and hydraulic properties, including thermal conductivity and the soil water retention curve. These properties were used to constrain a forward numerical model that included coupled treatment of relevant variable saturation thermal and hydraulic physics. Even though the investigated vadose zone is intermittent and relatively shallow ( 20 cm), the results illustrate how error can be introduced into heat‐transport calculations if unsaturated conditions are not taken into account.     

Heterogeneous characteristics of streambed saturated hydraulic conductivity of the Touchet River,south eastern Washington,USA

Traditionally a streambed is treated as a layer of uniform thickness and low saturated hydraulic conductivity (K) in surface‐ and ground‐water studies. Recent findings have shown a high level of spatial heterogeneity within a streambed and such heterogeneity directly affects surface‐ and ground‐water exchange and can have ecological implications for biogeochemical transformations, nutrient cycling, organic matter decomposition, and reproduction of gravel spawning fish. In this study a detailed field investigation of K was conducted in two selected sites in Touchet River, a typical salmon spawning stream in arid south eastern Washington, USA. In‐stream slug tests were conducted to determine K following the Bouwer and Rice method. For the upper and lower sites, each 50 m long and 9 m wide and roughly 20 m apart, a sampling grid of 5 m longitudinally and 3 m transversely was used. The slug tests were performed for each horizontal coordinate at 0·3–0·45, 0·6–0·75, 0·9–1·05 and 1·2–1·35 m depth intervals unless a shallower impenetrable obstruction was encountered. Additionally, water levels were measured to obtain vertical hydraulic gradient (VHG) between each two adjacent depth intervals. Results indicated that K ranged over three orders of magnitude at both the upper and lower sites and differed between the two sites. At the upper site, K did not differ significantly among different depth intervals based on nonparametric statistical tests for mean, median, and empirical cumulative distribution, but the spatial pattern of K varied among different depth intervals. At the lower site, K for the 0·3–0·45 m depth interval differed statistically from those at other depth intervals, and no similar spatial pattern was found among different depth intervals. Zones of upward and downward water flow based on VHG also varied among different depth intervals, reflecting the complexities of the water flow regime. Detailed characterization of the streambed as attempted in this study should be helpful in providing information on spatial variations of streambed hydraulic properties as well as surface‐ and ground‐water interaction. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical modelling of stream–aquifer interaction: Quantifying the impact of transient streambed permeability and aquifer heterogeneity

Stream–aquifer interaction plays a vital role in the water cycle, and a proper study of this interaction is needed for understanding groundwater recharge, contaminants migration, and for managing surface water and groundwater resources. A model‐based investigation of a field experiment in a riparian zone of the Schwarzbach river, a tributary of the Rhine River in Germany, was conducted to understand stream–aquifer interaction under alternative gaining and losing streamflow conditions. An equivalent streambed permeability, estimated by inverting aquifer responses to flood waves, shows that streambed permeability increased during infiltration of stream water to aquifer and decreased during exfiltration. Aquifer permeability realizations generated by multiple‐point geostatistics exhibit a high degree of heterogeneity and anisotropy. A coupled surface water groundwater flow model was developed incorporating the time‐varying streambed permeability and heterogeneous aquifer permeability realizations. The model was able to reproduce varying pressure heads at two observation wells near the stream over a period of 55 days. A Monte Carlo analysis was also carried out to simulate groundwater flow, its age distribution, and the release of a hypothetical wastewater plume into the aquifer from the stream. Results of this uncertainty analysis suggest (a) stream–aquifer exchange flux during the infiltration periods was constrained by aquifer permeability; (b) during exfiltration, this flux was constrained by the reduced streambed permeability; (c) the effect of temporally variable streambed permeability and aquifer heterogeneity were found important to improve the accurate capture of the uncertainty; and (d) probabilistic infiltration paths in the aquifer reveal that such pathways and the associated prediction of the extent of the contaminant plume are highly dependent on aquifer heterogeneity.     

Influence of water flux and redox conditions on chlorobenzene concentrations in a contaminated streambed

Significant natural attenuation may occur on the passage of groundwater plumes through streambed sediments because of the transition from anaerobic to aerobic conditions and an increased microbial activity. Varying directions and magnitudes of water flow in the streambed may enhance or inhibit the supply of oxygen to the streambed and thus influence the redox zoning. In a field study at a small stream in the industrial area of Bitterfeld‐Wolfen, we observed the variability of hydraulic gradients, streambed temperatures, redox conditions and monochlorobenzene (MCB) concentrations in the streambed over the course of 5 months. During the observation period, the hydrologic conditions changed from losing to gaining. Accordingly, the temperature‐derived water fluxes changed from recharge to discharge. Redox conditions were highly variable between ? 170 and 368 mV in the shallow streambed at a depth of 0·1 m below the streambed surface. Deeper in the streambed, at depths of 0·3 m and 0·5 m, the redox conditions were more stable between ? 198 and ? 81 mV and comparable to those typically found in the aquifer. MCB concentrations in the streambed at 0·3 and 0·5 m depth increased with increasing upward water flux. The MCB concentrations in the shallow streambed at 0·1 m depth appeared to be independent of the hydrologic conditions suggesting that degradation of MCB may have occured. Copyright © 2010 John Wiley & Sons, Ltd.     

Mapping river bathymetries: Evaluating topobathymetric LiDAR survey

Advances in topobathymetric LiDARs could enable rapid surveys at sub-meter resolution over entire stream networks. This is the first step to improving our knowledge of riverine systems, both their morphology and role in ecosystems. The Experimental Advanced Airborne Research LiDAR B (EAARL-B) system is one such topobathymetric sensor, capable of mapping both terrestrial and aquatic systems. Whereas the original EAARL was developed to survey littoral areas, the new version, EAARL-B, was also designed for riverine systems but has yet to be tested. Thus, we evaluated the ability of EAARL-B to map bathymetry and floodplain topography at sub-meter resolution in a mid-size gravel-bed river. We coupled the EAARL-B survey with highly accurate field surveys (0.03 m vertical accuracy and approximately 0.6 by 0.6 m resolution) of three morphologically distinct reaches, approximately 200 m long 15 m wide, of the Lemhi River (Idaho, USA). Both point-to-point and raster-to-raster comparisons between ground and EAARL-B surveyed elevations show that differences (ground minus EAARL-B surveyed elevations) over the entire submerged topography are small (root mean square error, RMSE, and median absolute error, M, of 0.11 m), and large differences (RMSE, between 0.15 and 0.38 m and similar M) are mainly present in areas with abrupt elevation changes and covered by dense overhanging vegetation. RMSEs are as low as 0.03 m over paved smooth surfaces, 0.07 m in submerged, gradually varying topography, and as large as 0.24 m along banks with and without dense, tall vegetation. EAARL-B performance is chiefly limited by point density in areas with strong elevation gradients and by LiDAR footprint size (0.2 m) in areas with topographic features of similar size as the LiDAR footprint. © 2018 John Wiley & Sons, Ltd.     

Model of particulate organic carbon transport in an agriculturally impacted stream

The present contribution focuses on modeling the total particulate organic carbon (POC) and benthic POC transport from a lowland stream impacted by agricultural land‐use. A mass balance, reach scale model is verified that accounts for water, sediment and POC transport, sediment and POC temporary storage and exchange with the streambed, and production and degradation of carbon pools in the benthos. We found that the POC load is highly variable during individual hydrologic events and is influenced by transport of mixed carbon sources including upland, streambank and benthic POC sources. Benthic POC stocks and transport were found to vary seasonally and annually but are in a state of long‐term equilibrium. Equilibrium is governed by negative feedback mechanisms whereby high POC export due to extreme hydrologic events and high‐frequency hydrologic events reduces benthic POC stocks and inhibits benthic POC growth. Benthic POC accounted for 4 tC y^?1 or 22% of the total annual POC loading in the stream's main stem and 8.9 tC y^?1 or 48% of the POC yield for the entire watershed. These results suggest that further attention should be given to benthic‐derived POC when budgeting stream ecosystem carbon for low‐order stream systems. Copyright © 2012 John Wiley & Sons, Ltd.     

Estimating streambed parameters for a disconnected river

Evaluation of stream–aquifer interaction and water balance for a catchment often requires specific information on streambed parameters, such as streambed hydraulic conductivity, seepage flux across the streambed and so on. This paper describes a simple, inexpensive instrument that is used to measure these streambed parameters under the condition of a stream disconnected from groundwater. Our method includes a seepage cylinder for simulation of river water depth. The proposed method was applied to estimate the vertical hydraulic conductivity of a streambed and the changes in vertical seepage rate from stream to groundwater with varied stream water depth in the Manasi River of Xinjiang Uygur Autonomous Region, China. The vertical hydraulic conductivities of the streambed determined from 12 sites along the Manasi River vary from 1.01 to 29.m/day where the stream disconnects from the groundwater. The experimental results suggest that there are two kinds of relations between the vertical seepage rate and the simulated stream water depth. One is a linear relation between the two variables with low Reynolds numbers (less than 10); the other is a nonlinear relation (exponential relation) between the two variables with larger Reynolds numbers (greater than 10). This second relationship is quite different from the traditional model that usually calculates the vertical seepage rate from stream to groundwater under the condition of disconnection using a linear relation (Darcy's Law). Our results suggest that a linear relation can only be used for a limited range of river water depth. This method gives a convenient tool for rapidly estimating the streambed hydraulic conductivity and the changes in the vertical seepage rate across streambed with varied stream water depths for the case of a stream disconnected from groundwater. Copyright © 2013 John Wiley & Sons, Ltd.