Search

Peatland drainage may degrade system resilience to high intensity, soil-consuming fires. Peat soil fires are unique in that they can smoulder vertically through the soil column, with a multitude of consequences including large carbon emissions, altered hydrology, and dramatic shifts in vegetation communities. In this work we developed and verified a new method to model peat burn depths with readily available water level and peat hydraulic property data at the Great Dismal Swamp National Wildlife Refuge (VA and NC, USA). To model peat burn depths across 11 sites in the Great Dismal Swamp National Wildlife Refuge we combined water table time series data and soil moisture release curves, developed at multiple depths, with moisture-to-ignition thresholds. A subset of the results from this empirical modelling approach of peat burn depth severity were compared against those made using a mechanistic model of soil moisture, HYDRUS 1-D. By comparing modelled burn depth potentials between these two approaches for a range of peats, we confirmed that our simpler, water table-based approach had similar performance to HYDRUS 1-D in drained and degraded peats, like those found in the Great Dismal Swamp National Wildlife Refuge. A comparative analysis of modelled burn depths across our study site found that water table position and peat water holding capacity were the key governing controls on burn depth potential. Our findings suggest that drainage weakens both short- and long-term controls on peat burn depths by reducing soil moisture and by decreasing peat water holding capacity. This new approach offers land managers with an additional tool for assessing risk while offering insight into the drivers of peatland wildfire severity. Published version Public domain – authored by a U.S. government employee

Conservation management practices - including agroforestry, cover cropping, no-till, reduced tillage, and residue return - have been applied for decades to control surface runoff and soil erosion, yet results have not been integrated and evaluated across cropping systems. In this study we collected data comparing agricultural production with and without conservation management strategies. We used a bootstrap resampling analysis to explore interactions between practice type, soil texture, surface runoff, and soil erosion. We then used a correlation analysis to relate changes in surface runoff and soil erosion to 13 other soil health and agronomic indicators, including soil organic carbon, soil aggregation, infiltration, porosity, subsurface leaching, and cash crop yield. Across all conservation management practices, surface runoff and erosion had respective mean decreases of 67% and 80% compared with controls. Use of cover cropping provided the largest decreases in erosion and surface runoff, thus emphasizing the importance of maintaining continuous vegetative cover on soils. Coarse- and medium-textured soils had greater decreases in both erosion and runoff than fine-textured soils. Changes in surface runoff and soil erosion under conservation management were highly correlated with soil organic carbon, aggregation, porosity, infiltration, leaching, and yield, showing that conservation practices help drive important interactions between these different facets of soil health. This study offers the first large-scale comparison of how different conservation agriculture practices reduce surface runoff and soil erosion, and at the same time provides new insight into how these interactions influence the improvement or loss of soil health. (C) 2021 International Research and Training Center on Erosion and Sedimentation, China Water and Power Press, and China Institute of Water Resources and Hydropower Research. Publishing services by Elsevier B.V. on behalf of KeAi Communications Co. Ltd. Research, Terrestrial Ecosystem Sciences Program [DE-AC05-76RL01830]; Yangling Vocational & Technical College 2019 Natural Science Research Fund Project [A2019048]; U.S. Department of Agriculture NRCS Conservation Innovation Grant [69-3A75-14-260]; Virginia Agricultural Experiment Station; Hatch Program of the National Institute of Food and Agriculture, U.S. Department of Agriculture Published version Research as part of the Terrestrial Ecosystem Sciences Program, under contract DE-AC05-76RL01830. Can Du was supported by the Yangling Vocational & Technical College 2019 Natural Science Research Fund Project, under grant number: A2019048. Ryan Stewart was supported by the U.S. Department of Agriculture NRCS Conservation Innovation Grant (No. 69-3A75-14-260) and the Vir-ginia Agricultural Experiment Station and the Hatch Program of the National Institute of Food and Agriculture, U.S. Department of Agriculture.

Subsurface injection of liquid manure is a no-till compatible application method that may fundamentally change nutrient cycling dynamics. This study investigated subsurface injection and surface broadcast application methods on soil nitrate-nitrogen (NO3–-N) at 0, 1, 2, 4, 6, 8, 10, 14, 18, and 22 weeks after first-year applications of liquid dairy slurry in 0 to 15 cm and 15 to 30 cm sampling depths. The effect on corn (Zea mays L.) grain yield and quality parameters was considered. Soil health was assessed by proxy on weeks 0, 2, 8, 14, and 22 by two microbial indicators sensitive to change: carbon (C) mineralization and microbial biomass C. Soil NO3–-N and the microbial indicators were also measured at distances from an injected band of manure: In-band, 10 cm, 20 cm, and 36 cm. The injection application resulted in significantly (p 0.05) different between treatments. Carbon mineralization was significantly (p 0.05) different from the no-manure control. Microbial biomass C was significantly (p

Soil organic matter (SOM) retains and attenuates many contaminants; however, its interactions with neonicotinoid insecticides under field conditions remain poorly understood. The goal of this study was to determine if SOM influences the persistence or leaching of two neonicotinoid insecticides: thiamethoxam (TMX) and its transformation-product clothianidin (CLO). Thiamethoxam-coated soybean [Glycine max (L.) Merr.] was planted into a clay soil containing different soil organic carbon (SOC) concentrations. Leachate and soil samples were collected for 10 wk after planting and were analyzed for insecticide concentrations using liquid chromatography–tandem mass spectrometry. Single and multiple linear regressions were performed between SOC, leaching volumes, and measured insecticide concentrations, focusing on rainfall events near the beginning, middle, and end of the study. Correlations were also tested between SOC and cumulative mass of leached insecticides. Neither SOC nor per-event leachate volumes explained variability in TMX leaching or residual CLO concentrations in soils; however, by the conclusion of the study residual thiamethoxam concentrations in soil were negatively correlated with cumulative volume of leached water. Initially, the concentration and total mass of leached CLO were significantly and negatively correlated with SOC content; however, this effect faded with time. Leachate dynamics also affected CLO transport, with positive correlations between leachate volume and CLO concentration during the latter events. This analysis demonstrates that SOM can reduce peak loading of neonicotinoids but may not alter cumulative leaching over the entire growing season. Funding for this work was provided by the USDA National Institute of Food and Agriculture (Grant 2018-67019-27851) and by the Virginia Agricultural Experiment Station and the Hatch Program of the USDA National Institute of Food and Agriculture (1026126). Published version

Peatland functions (e.g., carbon sequestration and flora diversity) are largely driven by soil moisture dynamics and thus dependent on interactions between hydrologic regimes and organic soil properties. Understanding these interactions is particularly important in drained peatlands, where drier conditions may alter soil properties with feedbacks to soil water retention and associated ecosystem functions. In this work, we focused on the Great Dismal Swamp (GDS) in Virginia, USA, a historically drained, temperate peatland with ongoing hydrologic restoration efforts. Two distinct soil layers varying in thickness exist at GDS: an upper layer with subangular blocky structure thought to be a result of past drainage, and a sapric lower layer with a massive structure more representative of an undisturbed state. To understand the occurrence and consequences of these distinct layers, we used continuous water table data and analysed soil physical and hydraulic properties to characterize soil profiles at 16 locations. We found significant differences between layer properties, where upper layers had lower fibre and organic matter contents and higher bulk densities. Further, moisture release curves demonstrated lower water retention in upper layers compared with lower layers and key differences in pore structure, with upper layers having higher macroporosity. Upper layers varied in thickness across sampling locations (similar to 0.30 to 1.0 m) with a transition to lower soil layers typically occurring at depths below contemporary water level observations, suggesting that the upper layer may be a result of historical drainage and deeper water table conditions. Yet, upper layers with more frequent saturation exhibited higher water retention and lower macroporosity compared with drier upper layers, thus indicating potential recovery following re-wetting efforts. These findings highlight how past drainage influences soil properties and water retention, with important implications for current management objectives at GDS and other drained peatland systems. U.S. Fish and Wildlife ServiceUS Fish & Wildlife Service Published version This work was funded by the U.S. Fish and Wildlife Service. Special thanks to Brooke Bass, Gavin Washburn, Karen Balentine, and Daniel Franz for their data collection support and to David Mitchem and Brandon Lester for their laboratory assistance.

Liquid manures are typically applied via surface broadcasting; however, subsurface injection is an alternative characterized by greater nutrient retention and a spatially distinct application pattern, altering management strategies and nutrient cycling dynamics. Thus, a field study was conducted from spring of 2016 through fall of 2018 on seven sites to assess pre-sidedress nitrate test (PSNT) methodology, seasonal soil nitrate nitrogen (NO3−-N) trends, corn (Zea mays L.) silage and grain yield, estimated milk production via MILK 2006, and biological soil health among surface broadcast and subsurface injection applications of dairy slurry. A weighted sampling method had a coefficient of variation of 37%, ~8% higher relative to random (28%) and equispaced (30%) sampling methods. Soil NO3−-N was greater in 7 out of 25 measurements under subsurface injection and 30% higher under injection on average during the corn PSNT. There were no significant differences in crop yield or milk production between surface and injected slurry applications, but means were always higher for injection. Biological soil health tests were highly variable, and analyzing carbon mineralization (C-min) took considerably more time than other tests. There were no significant differences in C-min between manure application methods; although, mineralization values increased with soil organic matter. Estimated microbial biomass was on average 46% lower under subsurface injection relative to surface broadcast in 2017, but results were inconsistent in 2016 and 2018. Overall, the biological indicators of soil health were not productive in showing differences between application methods. Nevertheless, it is apparent that injection can decrease chemical sidedress N applications, and either the standard method of PSNT soil sampling or an equispaced method can be used in injected fields.

Land application of manure, while beneficial to soil health and plant growth, can lead to an overabundance of nutrients and introduction of emerging contaminants into agricultural fields. Compared with surface application of manure, subsurface injection has been shown to reduce nutrients and antibiotics in surface runoff. However, less is known about the influence of subsurface injection on the transport and persistence of antibiotic-resistant microorganisms. We simulated rainfall to field plots at two sites (one in Virginia and one in Pennsylvania) 1 or 7 d after liquid dairy manure surface and subsurface application (56 Mg ha–1) and monitored the abundance of culturable antibiotic-resistant fecal coliform bacteria (ARFCB) in surface runoff and soils for 45 d. We performed these tests at both sites in spring 2018 and repeated the test at the Virginia site in fall 2019. Manure subsurface injection, compared with surface application, resulted in less ARFCB in surface runoff, and this reduction was greater at Day 1 after application compared with Day 7. The reductions of ARFCB in surface runoff because of manure subsurface injection were 2.5–593 times at the Virginia site in spring 2018 and fall 2019 and 4–5 times at the Pennsylvania site in spring 2018. The ARFCB were only detectable in the 0-to-5-cm soil depth within 14 d of manure surface application but remained detectable in the injection slits of manure subsurface-injected plots even at Day 45. This study demonstrated that subsurface injection can significantly reduce surface runoff of ARFCB from manure-applied fields. National Institute of Food and Agriculture,Grant/Award Number: 2017-67019-26401; USDA REISS, Grant/Award Number:PEN04574/1004448 Published version

Hydraulic conductivity (K) is a key hydrologic parameter widely recognized to be difficult to estimate and constrain, with little consistent assessment in disturbed, urbanized soils. To estimate K, it is either measured, or simulated by pedotransfer functions, which relate K to easily measured soil properties. We measured K in urbanized soils by double-ring infiltrometer (K(dring)), near-saturated tension infiltrometry (K(minidisk)), and constant head borehole permeametry (K(borehole)), along with other soil properties across the major soil orders in 12 United States cities. We compared measured K with that predicted from the pedotransfer function, ROSETTA. We found that regardless of soil texture, K(dring) was consistently larger than K(minidisk); with the latter having slightly less sample variance. K(borehole) was dependent upon specific subsurface conditions, and contrary to common expectations, did not always decrease with depth. Based on either soil textural class, or percent textural separates (sand, silt clay), ROSETTA did not accurately predict measured K for surface nor subsurface soils. We go on to discuss how K varies in urban landscapes, the role of measurement methods and artifacts in the perception of this metric, and implications for hydrologic modeling. Overall, we aim to inspire consistency and coherence when addressing K-related challenges in sustainable urban water management.

Soil water repellency induced by wildfires can alter hydraulic properties and hydrologic processes; however, the persistence and vertical position (i.e., depth) of water‐repellent layers can vary between systems and fires, with limited understanding of how those variations affect infiltration processes. This study occurred in two forested locations in the south‐central Appalachian Mountains that experienced wildfires in late 2016: Mount Pleasant Wildfire Refuge, Virginia, and Chimney Rock State Park, North Carolina. In each location, sites were selected to represent unburned conditions and low to moderate burn intensities. At each site, we measured the soil water repellency at the surface (ash layer or O horizon) and ~2 cm below the surface (A horizon) using the water drop penetration time method (n = 10–14). Soil water content was also measured over the upper 10 cm of the soil (n = 10), and infiltration tests were conducted using a tension infiltrometer (n = 6–8). The results showed that soil repellency was highest in the surface layer at the Mount Pleasant location and was highest in the subsurface layer at the Chimney Rock location. Soil water content was lower in unburned soil than in burned soil, especially for measurements taken immediately postfire, with soil water content negatively correlated with water repellency. Water repellency in the surface layer significantly reduced relative infiltration rates (estimated as differences between initial and steady‐state rates), whereas subsurface water repellency did not affect relative infiltration. As a result, water repellency persisted longer in sites with surface as opposed to subsurface water repellency. Finally, differences between burned and unburned sites showed that although the wildfires increased the occurrence of water repellency, they did not alter the underlying relationship between relative infiltration and water repellency of the surface soil.

Vertisols shrink and swell with changes in soil moisture, influencing hydraulic properties. Vertisols are often in floodplains, yet the importance of flooding as a source of soil moisture remains poorly understood. We used blue dye and deuterated water as tracers to determine the role of the macropore network in matrix recharge under artificial flood durations of 3 and 31 days in large soil monoliths extracted from a forested soil. Gravimetric soil moisture content increased by 47% in the first three days, then increased only 3.5% from day 3-31. Post-flood moisture content was greatest in the organic-rich, top 10 cm and was lower at 10-75 cm where organic matter was less. Deuterium concentration revealed that soil moisture in the top 10 cm was quickly dominated by artificial flood water, but at depth remained

Novel in situ sensor technologies can measure water chemistry at high temporal frequencies, yet few studies have evaluated how installation affects measurements. In this study, we assessed the effects of commonly used protective housings on in situ sensor readings. Working in two mountain streams, we colocated specific conductance sensors in four different housing types that varied in openings for water exchange (mesh, screen, holes, and open). We compared measured conductance values through time and performed repeated salt tracer additions to evaluate the influence of housing type on calculated discharge. Sensors readings in mesh and, to a lesser extent, screen housings frequently diverged from housings with larger openings (i.e., holes and open), indicating reduced water exchange between stream water and housed sensors. Further, mesh and screen housings recorded more damped and delayed response to salt tracer additions compared to the other two housings, resulting in markedly different discharge values. From these findings, we recommend that water chemistry sensors should be deployed in a protective housing with large openings for sufficient water exchange. Accepted version

International audience; Soil water repellency can substantially alter hydrologic processes, particularly the abilityof soils to infiltrate water. Water repellency often changes through time, making it difficult to simulate infiltration behaviors of water-repellent soils using standard models. Here, we propose a simple rate-based correction term that starts with a value of zero at the beginning of the infiltration process (t = 0) and asymptotically approaches 1 as time increases, thus simulating decreasing soil water repellency through time. The correction term can be used with any infiltration model. For this study, we selected a simple two-term infiltration equation and then, using two data sets of infiltration measurements conductedin soils with varying water repellency, compared model error with versus without the added term. The correction substantially reduced model error, particularly in more repellent soils. At the same time, the rate constant parameter introduced in the new model may be useful to better understand dynamics of soil water repellency and to provide more consistent interpretations of hydraulic properties in water-repellent soils.

International audience; Comprehensive infiltration models can simultaneously describe transient and steady-state infiltration behaviors, and therefore can be applied to a range of experimental conditions. However, satisfactory model accuracy re- quires proper parameterization, including estimating the transition time from transient to steady-state flow conditions (τcrit). This study focused on improving the estimation of two parameters – τcrit and a second constant called a – used in a comprehensive, explicit, two-term model for single ring infiltration (hereafter referred to as the SA model). Different studies have recommended that a should be as low as 0.45 to as high as 0.91. Furthermore, τcrit is often obtained a-priori by assuming that steady-state conditions are reached before the end of an infiltration run. However, there has not been a systematic analysis of those terms for different soils and infiltration conditions. To investigate these open issues related to the use of the SA model, here we introduce a novel, iterative method for estimating τcrit and the parameter a. We then applied this method to both analytical and experimental infiltration data, and compared it with two existing empirical methods. The analytical infil- tration experiments showed that τcrit was approximately 1.5 times larger than the maximum validity time of a similar two-term transient infiltration model. Further, the iterative method for obtaining τcrit had minimal effects on the a term, which varied between 0.706 and 0.904 and was larger for finer soils and when small water sources were used. Application of the proposed method was less efficient with experimental data. Only ~ 33% of the experiments yielding plausible estimates of a (i.e., a < 1), indicating that these infiltration model parameters often have high uncertainty. The successful runs indicated that a depended on the rate at which the initial infiltration rate approached the final infiltration rate. Depending on the fitting algorithm used, a had mean values of 0.74–0.78, which were intermediate between those suggested by previous studies. Altogether, these findings expand the applicability of the SA model by providing new methods for estimating τcrit and by showing that a does not need to be fixed a-priori. We expect that these advances will result in more reliable estimations of soil hydrodynamic parameters, including hydraulic conductivity.

Uncontrolled overland flow drives flooding, erosion, and contaminant transport, with the severity of these outcomes often amplified in urban areas. In pervious media such as urban soils, overland flow is initiated via either infiltration-excess (where precipitation rate exceeds infiltration capacity) or saturation-excess (when precipitation volume exceeds soil profile storage) mechanisms. These processes call for different management strategies, making it important for municipalities to discern between them. In this study, we derived a generalized one-dimensional model that distinguishes between infiltration-excess overland flow (IEOF) and saturation-excess overland flow (SEOF) using Green-Ampt infiltration concepts. Next, we applied this model to estimate overland flow generation from pervious areas in 11 U.S. cities. We used rainfall forcing that represented low- and high-intensity events and compared responses among measured urban versus predevelopment reference soil hydraulic properties. The derivation showed that the propensity for IEOF versus SEOF is related to the equivalence between two nondimensional ratios: (a) precipitation rate to depth-weighted hydraulic conductivity and (b) depth of soil profile restrictive layer to soil capillary potential. Across all cities, reference soil profiles were associated with greater IEOF for the high-intensity set of storms, and urbanized soil profiles tended towards production of SEOF during the lower intensity set of storms. Urban soils produced more cumulative overland flow as a fraction of cumulative precipitation than did reference soils, particularly under conditions associated with SEOF. These results will assist cities in identifying the type and extent of interventions needed to manage storm water produced from pervious areas.

International audience; We present a new open-source and modular instrumentation package composed of up to ten automatic in- filtrometers connected to data acquisition systems for automatic recording of multiple infiltration experiments. The infiltrometers are equipped with differential transducers to monitor water level changes in a Mariotte re- servoir, and, in turn, to quantify water infiltration rates. The data acquisition systems consist of low-cost components and operate on the open-source microcontroller platform Arduino. The devices were tested both in the laboratory and on different urban and agricultural soils in France and India. More specifically, we tested three procedures to treat the transducers readings, including a filtering algorithm that substantially improved the ability to determine cumulative infiltration from raw data. We combined these three procedures with four methods for estimating the soil parameters from infiltrometer data, showing pros and cons of each scenario. We also demonstrated advantages in using the automatic infiltrometers when infiltration measurements were hin- dered by: i) linearity in cumulative infiltration curves owing to gravity-driven flow, ii) an imprecise description of the transient state of infiltration, and iii) the occurrence of soil water repellency. The use of the automatic infiltrometers allows the user to obtain more accurate estimates of soil hydraulic parameters, while also reducing the amount of effort needed to run multiple experiments.

International audience; The macroscopic capillary length, λc, is a fundamental soil parameter expressing the relative importance of the capillary over gravity forces during water movement in unsaturated soil. In this investigation, we propose a simple field method for estimating λc using only a single-ring infiltration experiment of the Beerkan type and measurements of initial and saturated soil water contents. We assumed that the intercept of the linear regression fitted to the steady-state portion of the experimental infiltration curve could be used as a reliable predictor of λc. This hypothesis was validated by assessing the proposed calculation approach using both analytical and field data. The analytical validation demonstrated that the proposed method was able to provide reliable λc estimates over a wide range of soil textural characteristics and initial soil water contents. The field testing was performed on a large database including 433 Beerkan infiltration experiments, with the 99% of the experiments yielding realistic λc values. The generated λc values were then used in conjunction with four different methods for estimating saturated soil hydraulic conductivity, Ks. Estimated Ks values were close to those generated by a reference method, with relative error < 25% in nearly all cases. By comparison, assuming constant or soil-dependent λc values caused relative errors in Ks of up to 600%. Altogether, the proposed method constitutes an easy solution for estimating λc, which can improve our ability to estimate Ks in the field.

Widespread observations of ecohydrological separation are interpreted by suggesting that water flowing through highly conductive soil pores resists mixing with matrix storage over periods of days to months (i.e., two 'water worlds' exist). These interpretations imply that heterogeneous flow can produce ecohydrological separation in soils, yet little mechanistic evidence exists to explain this phenomenon. We quantified the separation between mobile water moving through preferential flow paths versus less mobile water remaining in the soil matrix after free-drainage to identify the amount of preferential flow necessary to maintain a two water world's scenario. Soil columns of varying macropore structure were subjected to simulated rainfall of increasing rainfall intensity (26 mm h(-1), 60 mm h(-1), and 110 mm h(-1)) whose stable isotope signatures oscillated around known baseline values. Prior to rainfall, soil matrix water delta H-2 nearly matched the known value used to initially wet the pore space whereas soil delta O-18 deviated from this value by up to 3.4 parts per thousand, suggesting that soils may strongly fractionate O-18. All treatments had up to 100% mixing between rain and matrix water under the lowest (26 mm h(-1)) and medium (60 mm h(-1)) rainfall intensities. The highest rainfall intensity (110 mm h(-1)), however, reduced mixing of rain and matrix water for all treatments and produced significantly different preferential flow estimates between columns with intact soil structure compared to columns with reduced soil structure. Further, artificially limiting exchange between preferential flow paths and matrix water reduced bypass flow under the most intense rainfall. We show that (1) precipitation offset metrics such as lc-excess and d-excess may yield questionable interpretations when used to identify ecohydrological separation, (2) distinct domain separation may require extreme rainfall intensities and (3) domain exchange is an important component of macropore flow. Virginia Agricultural Experiment Station; Hatch Program of the National Institute of Food and Agriculture, USDA [1007839] Published version Funding for this work was provided in part by the Virginia Agricultural Experiment Station and the Hatch Program of the National Institute of Food and Agriculture, USDA (1007839). We would like to thank Aaron Cleveland for his help and hard work conducting this experiment. Thanks to Matthias Sprenger for initial criticism and discussion which helped improve the manuscript. Thanks to Kevin McGuire for the motivational discussions regarding mixing and preferential flow calculations with stable isotopes, all of which influenced the structure of this study. Thank you, Kelly Peeler and Durelle Scott, for help with liquid isotope analysis. We thank the two anonymous reviewers for their valuable and detailed comments.

Made available in DSpace on 2022-04-28T19:41:56Z (GMT). No. of bitstreams: 0 Previous issue date: 2021-12-01 European Regional Development Fund Water-repellent soils usually experience water flow impedance during the early stage of a wetting process followed by progressive increase of infiltration rate. Current infiltration models are not formulated to describe this peculiar process. Similarly, simplified methods of soil hydraulic characterization (e.g., BEST) are not equipped to handle water-repellent soils. Here, we present an adaptation of the BEST method, named BEST-WR, for the hydraulic characterization of soils at any stage of water-repellency. We modified the Haverkamp explicit transient infiltration model, included in BEST for modeling infiltration data, by embedding a scaling factor describing the rate of attenuation of infiltration rate due to water repellency. The new model was validated using analytically generated data, involving soils with different texture and a dataset that included data from 60 single-ring infiltration tests. The scaling factor was used as a new index to assess soil water repellency in a Mediterranean wooded grassland, where the scattered evergreen oak trees induced more noticeable water repellency under the canopies as compared to the open spaces. The new index produced results in line with those obtained using the water drop penetration time test, which is one of the most widely test applied for quantifying soil water repellency persistence. Finally, we used BEST-WR to determine the hydraulic characteristic curves under both hydrophilic and hydrophobic conditions. Department of Agricultural Sciences University of Sassari, Viale Italia, 39A Université de Lyon UMR5023 Ecologie des Hydrosystèmes Naturels et Anthropisés CNRS ENTPE Université Lyon 1 School of Plant and Environmental Sciences Virginia Polytechnic Institute and State University Department of Land Air and Water Resources University of California Department of Architecture Design and Urban Planning University of Sassari, Viale Piandanna, 4 School of Agriculture São Paulo State University (UNESP) Fazenda Experimental Lageado Civil Engineering Department Engineering Faculty Munzur University School of Agriculture São Paulo State University (UNESP) Fazenda Experimental Lageado

Municipalities and agencies use green infrastructure to combat pollution and hydrological impacts (e.g., flooding) related to excess stormwater. Bioretention cells are one type of infiltration green infrastructure (GI) intervention that infiltrate and redistribute otherwise uncontrolled stormwater volume. However, the effects of these installations on the rest of the local water cycle is understudied; in particular, impacts on stormwater return flows and groundwater levels are not fully understood. In this study, full water cycle monitoring data was used to construct and calibrate a two-dimensional Richards equation model (HYDRUS-2D/3D) detailing hydrological implications of an unlined bioretention cell (Cleveland, Ohio) that accepts direct runoff from surrounding impervious surfaces. Using both pre- and post-installation data, the model was used to: 1) establish a mass balance to determine reduction in stormwater return flow, 2) evaluate GI effects on subsurface water dynamics, and 3) determine model sensitivity to measured soil properties. Comparisons of modeled versus observed data indicated that the model captured many hydrological aspects of the bioretention cell, including subsurface storage and transient groundwater mounding. Model outputs suggested that the bioretention cell reduced stormwater return flows into the local sewer collection system, though the extent of this benefit was attenuated during high inflow events that may have exhausted detention capacity. The model also demonstrated how, prior to bioretention cell installation, surface and subsurface hydrology were largely decoupled, whereas after installation, exfiltration from the bioretention cell activated a new groundwater dynamic. Still, the extent of groundwater mounding from the cell was limited in spatial extent, and did not threaten other subsurface infrastructure. Finally, the sensitivity analysis demonstrated that the overall hydrological response was regulated by the hydraulics of the bioretention cell fill material, which controlled water entry into the system, and by the water retention parameters of the native soil, which controlled connectivity between the surface and groundwater.

Subsurface flow dominates water movement from hillslopes to streams in most forested headwater catchments. Hewlett and Hibbert (1963, https://doi.org/10.1029/JZ068i004p01081) constructed an idealized hillslope model (0.91 × 0.91 × 15.0 m; 21.8°) using reconstituted C horizon soil to investigate importance of interflow, a type of subsurface flow. They saturated the model, covered it to prevent evaporation, and allowed free drainage for 145 days. The resulting recession drainage curve suggested two phases: fast drainage of saturated soil in the first 1.5 days and then slow drainage of unsaturated soil. Hydrologists interpreted the latter as evidence interflow could sustain baseflow, even during extended drought. Since that experiment, typical forest vegetation grew in the model, providing root and litter inputs for 55 years. We removed all aboveground live biomass and repeated the experiment physically and numerically (HYDRUS‐2D), hypothesizing that pedogenesis would change the drainage curve and further elucidate the role of unsaturated flow from hillslopes. Contrary to this hypothesis, drainage curves in our twice‐repeated physical experiments and numerical simulation were unchanged for the first ~10 days, indicating pedogenesis and biological processes had not largely altered bulk hydraulic conductivities or soil moisture release characteristics. However, drainage unexpectedly ceased after about 2 weeks (14.3 ± 2.5 days), an order of magnitude sooner than in the original experiment, due to an apparent leak in the hillslope analogous to commonly observed bedrock fractures in natural systems. Thus, our results are a more natural recession behavior that highlight how incorporation of alternative hydrologic outputs can reduce drainage duration and volume from soils to baseflow.

© 2019 Elsevier B.V. The objective of this paper was to evaluate a recently proposed comprehensive model for three-dimensional single-ring infiltration and its suitability for estimating soil hydraulic properties. Infiltration data from four different soils with contrasting characteristics were inverted to estimate field-saturated soil hydraulic conductivity, K fs , values using a total of fourteen different scenarios. Those scenarios differed by: i) the way they constrained the macroscopic capillary length, λ, and the initial and saturated soil water contents, θ i and θ s , ii) the use of transient or steady-state data, and iii) the fitting methods applied to transient data. For comparative purposes, the SSBI method (Steady version of the Simplified method based on a Beerkan Infiltration run) was also applied. For validation purposes K fs data estimated from the different scenarios were compared with those values obtained by numerical inverse modeling with HYDRUS-2D/3D. This comparison identified Approaches 1 and 3, which respectively estimate K fs via optimization and using analytical expressions, as the most accurate methods. The steady-state scenario of Approach 4 and the SSBI method, both of which use a λ value of first approximation, appeared preferable for field campaigns aimed to sample remote or large areas, given that they do not need additional data and still provide acceptable estimates. The reliability of K fs data was also checked through a comparison with unsaturated hydraulic conductivity, K h , values measured in laboratory on extracted soil cores, in order to discriminate between theoretically possible (K fs > K h ) and impossible (K fs ≤ K h ) situations. Physically possible K fs values were always obtained with the exception of the crusted soil, where K fs < K h situations suggested that the crust layer reduced water flow during ponding experiments in the field. The new comprehensive model tested in this study represents a valuable tool for analyzing both transient and steady-state infiltration data, as well as experiments carried out with different depths of ponded water, ring sizes and ring insertion depths.

Pedotransfer functions (PTFs) provide point predictions of soil hydraulic properties from more readily measured soil characteristics, yet uncertainties and biases in measurement methods, sampling distributions, and boundary conditions can limit accuracy when estimating near-saturated hydraulic conductivity (Kn). These limitations may be particularly problematic in understudied urban landscapes that often contain altered hydraulic properties. To better treat deficiencies in PTF performance, we addressed three objectives, which were to: 1) develop PTFs to predict urban Kn, 2) assess bulk density and coarse fragments as explanatory variables; and 3) evaluate the predictive capability of these PTFs by comparing their output to measured hydraulic conductivity values from three other studies of urban soil hydraulics. We used artificial neural networks (ANN) and random forest (RF) approaches to predict urban Kn, with the training dataset including 307 tension infiltrometer tests and other measurements drawn from urban soil assessments in 11 U.S. cities. The PTFs utilized a hierarchy of inputs, starting with percentage sand, silt, clay, and then adding percentage coarse fragments and bulk density. The ANN models performed similar to the RF models, and all models exhibited similar or better predictive performance as models results collected from published articles. The inclusion of bulk density or coarse fragments did not improve accuracy over soil texture alone. Possible reasons for this result include low correlation between Kn and bulk density and the exclusion of large voids during flow measurements with tension infiltrometers. The models have been made available as an open-source software package to encourage adoption by users working in urban systems.

Shrink-swell soils crack and become compacted as they dry, changing properties such as bulk density and hydraulic conductivity. Multidomain models divide soil into independent realms that allow soil cracks to be incorporated into classical flow and transport models. Incongruously, most applications of multidomain models assume that the porosity distributions, bulk density, and effective saturated hydraulic conductivity of the soil are constant. This study builds on a recently derived soil shrinkage model to develop a new multidomain, dual-permeability model that can accurately predict variations in soil hydraulic properties due to dynamic changes in crack size and connectivity. The model only requires estimates of soil gravimetric water content and a minimal set of parameters, all of which can be determined using laboratory and/or field measurements. We apply the model to eight clayey soils, and demonstrate its ability to quantify variations in volumetric water content (as can be determined during measurement of a soil water characteristic curve) and transient saturated hydraulic conductivity, Ks (as can be measured using infiltration tests). The proposed model is able to capture observed variations in Ks of one to more than two orders of magnitude. In contrast, other dual-permeability models assume that Ks is constant, resulting in the potential for large error when predicting water movement through shrink-swell soils. Overall, the multidomain model presented here successfully quantifies fluctuations in the hydraulic properties of shrink-swell soil matrices, and are suitable for use in physical flow and transport models based on Darcy's Law, the Richards Equation, and the advection-dispersion equation.

Recognizing that elevated stream temperatures can harm aquatic organisms, regulatory agencies have begun to enforce thermal total daily maximum load (TMDL) limits on wastewater discharges. Infiltration discharge, where treated wastewater is allowed to infiltrate into the soil and then percolate toward the stream or river, represents a potential method for meeting temperature requirements and for supplementing flow. This study combines observed and simulated temperature data to assess the efficacy of using infiltration discharge to meet effluent temperature limits. Observational data collected during operation of a 0.15-ha pilot-scale infiltration discharge system revealed the pattern of groundwater heating beneath the site by the wastewater, with the largest temperature increases occurring near the infiltration point. Numerical simulations used to examine the long-term groundwater response to operation of a 5.5-ha large-scale infiltration wetland system confirmed this trend and demonstrated how mu...

Irrigation experiments on 12 instrumented field plots were used to assess the impact of dynamic soil crack networks on infiltration and run-off. During applications of intensity similar to a heavy rainstorm, water was seen being preferentially delivered within the soil profile. However, run-off was not observed until soil water content of the profile reached field capacity, and the apertures of surface-connected cracks had closed >60%. Electrical resistivity measurements suggested that subsurface cracks persisted and enhanced lateral transport, even in wet conditions. Likewise, single-ring infiltration measurements taken before and after irrigation indicated that infiltration remained an important component of the water budget at high soil water content values, despite apparent surface sealing. Overall, although the wetting and sealing of the soil profile showed considerable complexity, an emergent property at the hillslope scale was observed: all of the plots demonstrated a strikingly similar threshold run-off response to the cumulative precipitation amount. Copyright © 2014 John Wiley & Sons, Ltd.

[1] A soil's capillarity, associated with the parameter sorptivity, is a dominant control on infiltration, particularly at the onset of rainfall or irrigation. Many mathematical models used to estimate sorptivity are only valid for dry soils. This paper examines how sorptivity and its capillary component (as wetting front potential) change with initial degree of saturation. We capture these effects with a simple modification to the classic Green-Ampt model of sorptivity. The modified model has practical applications, including (1) accurately describing the relative sorptivity of a soil at various water contents and (2) allowing for quantification of a soil's saturated hydraulic conductivity from sorptivity measurements, given estimates of the soil's characteristic curve and initial water content. The latter application is particularly useful in soils of low permeability, where the time required to estimate hydraulic conductivity through steady-state methods can be impractical.

Efforts to prevent temperature increases in streams and rivers have led to regulation of the maximum permissible temperature for wastewater discharges. Hyporheic discharge treatment, where warm wastewater is allowed to infiltrate into the soil and percolate toward the stream, has been proposed as an alternative management practice, potentially providing increased streamflow at decreased temperature. To assess the feasibility of such a system, a pilot study was performed in which groundwater level, temperature, and chemistry were monitored to characterize percolation from a 0.15-ha infiltration wetland. The temperature patterns suggested that water seepage rates exceeded 1.6 md−1 near the wetland, while the chemistry data indicated that infiltrated water had laterally percolated more than 100 m during the study period. The water level and temperature measurements were also used to calibrate a software model. Observations and simulations showed that the local water table became directly connected t...

Recharge rates evaluated from pilot projects of surface managed aquifer recharge are extrapolated to develop design criteria for full-scale projects. Field experiments at surface-managed artificial aquifer recharge facilities in the Walla Walla River basin, Oregon, United States, were used to estimate recharges rates in relation to groundwater mounding and to the expansion in surface area of infiltration basins. Analysis of the results shows that in the case scenarios where the water table mounding does not reach the infiltration basin floor (thereby maintaining an unsaturated zone between the water table and the infiltration basin), the recharge rates from pilot tests scale linearly with the basin’s surface area expansion. However, in the case where the groundwater mound reaches the bottom of the basin floor (thereby providing a full hydraulic connection between the infiltration basins and the aquifer), recharge rates should be extrapolated using the perimeter of the infiltration basin. The expla...

We propose a novel, accurate quantification of precipitation and evaporation, as needed to understand fundamental hydrologic processes. Our system uses a collection vessel placed on top of a slender rod that is securely fixed at its base. As the vessel is deflected, either by manual perturbation or ambient forcing (for example, wind), its oscillatory response is measured, here by a miniature accelerometer. This response can be modeled as a damped mass-spring system. As the mass of water within the collection vessel changes, either through the addition of precipitation or by evaporative loss, the resonant frequency experiences an inverse shift. This shift can be measured and used to estimate the change in the mass of water. We tested this concept by creating a simple prototype which was used in field conditions for a period of 1 month. The instrument was able to detect changes in mass due to precipitation with an accuracy of approximately 1 mm.