Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Subsurface”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A simulation study of carbon storage with active reservoir management

As part of the Integrated Midcontinent Stacked Carbon Storage Hub (IMSCS-Hub) project led by Battelle Memorial Institute, a study was conducted to determine the feasibility of storing carbon dioxide (CO 2 ) in the stacked saline rock formations of the Sleepy Hollow Field (SHF), located in Red Willow County, southern Nebraska. A series of CO 2 injection simulation scenarios, with and without active reservoir management (ARM; brine extraction), were evaluated to investigate the feasibility of storing 50+ million tonnes (Mt) of CO 2 . The results indicated CO 2 injection combined with ARM may enable permanent storage of 50+ Mt of CO 2 . The area of review (AOR), the area in which underground sources of drinking water (USDWs) might be endangered during CO 2 injection, was assessed for the simulation scenarios. In comparison to a case without ARM, brine extraction resulted in a much smaller AOR, covering an area of 42 square miles (108.8 km 2 ), roughly one-fourth the size of an AOR resulting from CO 2 injection without ARM. Our findings, presented in this paper, indicate that CO 2 injection with ARM can improve the CO 2 storage capacity of a geologic storage complex up to 100% while also reduce the rate of pressure buildup in the subsurface, resulting in a 75% reduction in AOR. This may help in lowering carbon capture and storage project costs, risks, and effort needed to meet monitoring requirements for a storage project.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Quantitative interpretation of time-lapse seismic data at Farnsworth field unit: Rock physics modeling, and calibration of simulated time-lapse velocity responses

Here, this study investigates the contribution of fluid saturation variation to the time-lapse velocity response by performing fluid substitution modeling. The methodology is exemplified by the time-lapse seismic monitoring of carbon dioxide at Farnsworth field unit (FWU). In order to evaluate the fluid distribution in a matured oil reservoir, the Southwest Regional Partnership (SWP) acquired multiple vertical seismic profile (VSP) surveys at different times during the CO 2 –water alternatinggas (WAG) injection period. In this work, we present a thorough methodology for computing the elastic response of the saturated rock for different fluid saturations using a site-specific petro-elastic model (PEM). The output from the PEM was combined with results from a fluid compositional model to compute the seismic velocities at times corresponding to each VSP survey. To produce a calibrated simulated response, the measured time-lapse seismic velocities were integrated into the numerical simulation model. The mismatches between the predicted and measured time-lapse velocities were minimized through an iterative calibration process using a trained artificial neural network proxy (ANN) coupled with a particle swarm optimizer (PSO). Our study indicates that the hybrid optimization workflow can effectively perform the history matching. With an accurate prediction of the hydrodynamic properties, the migration of CO 2 within the subsurface was modeled by predicting the spatial velocity distribution for a radius of 305 m around the injection well. The technology demonstrated and the expertise gained from this study can guide similar CO 2 -WAG projects.

58 GEOSCIENCES↗
DOE OSTI2020

Dynamic domain kinematic modelling for predicting interflow over leaky impeding layers

Traditional Boussinesq or kinematic simulations of interflow (i.e., lateral subsurface flow) assume no leakage through the impeding layer and require a no-flow boundary condition at the ridge top. However, recent analyses of many interflow-producing landscapes indicate that leaky impeding layers are common, that most interflow percolates well before reaching the toe slope, and therefore, the downslope contributing length is shorter than the hillslope length. In watersheds characterised by perched interflow over a low conductivity layer through permeable topsoil, interflow with percolation may be modelled with a kinematic wave model using a mobile upslope boundary condition defining the hillslope portion contributing interflow to valleys. Here, we developed and applied a dynamic interflow model to simulate interflow using a downslope travel distance concept such that only the active contributing length is modelled at any time. The model defines a variable active area based on the depth of the perched layer, the topographic slope and the ratio of the hydraulic conductivity of topsoil to that of the impeding layer. It incorporates a two-layer soil moisture accounting water balance analysis, a pedo-transfer function, and percolation and evaporation routines to predict interflow rates in continuous and event-based scenarios. We tested the modelling concept on two sets of data (2-year dataset of rainfall observations for the continuous simulation and a multi-day irrigation experiment for the event simulation) from a 121-m-long open interflow collection trench on an experimental hillslope at the Savannah River Site, South Carolina. The continuous model simulation partially represented the observed interflow hydrograph and perched water depth in the experimental hillslope with correlation coefficients of 0.85 and 0.35, respectively. Model performance improved significantly at event-scale analysis. The modelling approach realistically represents interflow dynamics in hillslopes with leaky impeding layers and can be integrated into catchment-scale hydrology models for more detailed hillslope process modelling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Catchment storage and residence time in a periodically irrigated watershed

Understanding anthropogenic impacts on water storage and water flow pathways in catchments is an ongoing challenge in hydrology. Here, we study the dynamics of subsurface storage and residence time of water in a catchment in Berkeley, California, that is within a regional park but contains diverse land use within its perimeter, including a periodically irrigated golf course. Our study combines several isotopic tracers with water budget data to examine sources of water in a stream draining the site. Irrigation water, applied to a small area of the watershed, is a minor component of the water budget. However, geochemical tracers reveal that irrigation water is a significant fraction of stream flow downstream of the golf course during baseflow and during precipitation events. Isotopic tracers indicate that the watershed has a preference to release young water for stream flow generation, resulting in contrasting tritium ages for stream water and groundwater of 1.3 ± 0.5 year and 8.2 ± 1.7 year, respectively. We determined that the older water is a very small component (0.7%) of the stream water in the tail of an assumed exponential distribution. We used the seasonal variation of stable water isotopes in precipitation and stream water over two water years to explain the damping of the isotopic signature of stream water, which yields information about the catchment's response to the input signal. In conclusion, the methods described here may be applicable to other urban or suburban headwater catchments in areas with a component of non-natural recharge from, for example, leaky infrastructure, storm water routing or dry season irrigation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Vegetation water sources in California's Sierra Nevada (USA) are young and change over time, a multi-isotope ( δ 18 O, δ 2 H, 3 H) tracer approach.

Sierra Nevada forests transpire a significant amount of California's water resources, sparking interest in applying forest management to improve California's water supply. Determining the source water of evapotranspiration enables forest managers to make informed decisions. To this end, a significant interest in critical zone science is to develop new methods to work across time scales to predict subsurface water storage and use. Here, forest vegetation accessed young water and switched sources depending on availability, suggesting that forest drought vulnerability may depend on the range of water sources available (rain, snowmelt and deeply stored water). This finding also suggests that changes in transpiration rates may have immediate effects on water sources in close proximity to vegetation, and delayed effects on storage and runoff. New δ 18 O, δ 2 H and 3 H data were used to track precipitation, runoff, evapotranspiration and storage through the critical zone seasonally, including seasons where evapotranspiration and snowmelt were in phase (winter snowmelt) and out of phase (seasonally dry summer). The main source of this headwater catchment's runoff is derived from its meadow saturated zone water, which was dominated by snowmelt. Water that originated as snowmelt contributed to transpiration, unless other sources, such as recent rain, became available. In cases where xylem δ 18 O and δ 2 H signatures matched those of deeper saturated zone water, 3 H data showed that xylem water was distinctly younger than the deep saturated zone water. During 2016, which experienced relatively normal snowpack in winter and seasonally dry summer conditions, mean summer saturated zone water and vegetation water were similar in δ 18 O, –12.4 ± 0.04 ‰ and – 12.5 ± 0.3 ‰, respectively, but were distinctly different in 3 H, 5.5 ± 0.2 pCi/L and 13.7 ± 1.1 pCi/L, respectively. While δ 18 O shows that vegetation and meadow saturated zone water have similar origins, 3 H shows they have dissimilar ages.

58 GEOSCIENCES↗
DOE OSTI2021

The stable isotopes of natural waters at the Marcell Experimental Forest

Here, we present a new data set from the Marcell Experimental Forest (MEF) that compiles water isotope measurements from multiple research catchments, some of which have been studied since the 1960s. The MEF is located in northern Minnesota, USA, and is home to heavily studied and monitored forests, streams, bogs, and fens. Peat-forming systems (bogs and fens) are an important component of the MEF landscape and have a profound impact on the water cycle in these catchments. Within the last decade, analysis of stable isotopes of water (expressed as δD and δ 18 O) has been implemented to characterize the different components of the water budget, and to allow researchers to look at catchment and peatland-specific hydrologic effects in the watershed. This δD and δ 18 O data set of natural waters from MEF catchments is primarily composed of measurements from three peatlands (S1, S2, S6) during an 11-year period. More recently collection and analysis were expanded to also include samples from the Spruce and Peatland Responses Under Changing Environments (SPRUCE) project in the S1 bog, peatlands S3, S4, S5, as well as nearby lakes. We establish a local meteoric water line by analyzing the isotopic composition of precipitation, which fills a void in regional meteoric water lines for Minnesota. Furthermore, we establish baseline isotopic composition for bog outlet streams, bog porewater, aquifer groundwater, overland flow, subsurface stormflow, and snowpack, as well as runoff from the SPRUCE experimental chambers. These data are publicly available and will be expanded upon in the future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The Wasatch Environmental Observatory: A mountain to urban research network in the semi‐arid western US

Abstract The 2085 km 2 Jordan River Basin, and its seven sub‐catchments draining the Central Wasatch Range immediately east of Salt Lake City, UT, are home to an array of hydrologic, atmospheric, climatic and chemical research infrastructure that collectively forms the Wasatch Environmental Observatory (WEO). WEO is geographically nested within a wildland to urban land‐use gradient and built upon a strong foundation of over a century of discharge and climate records. A 2200 m gradient in elevation results in variable precipitation, temperature and vegetation patterns. Soil and subsurface structure reflect systematic variation in geology from granitic, intrusive to mixed sedimentary clastic across headwater catchments, all draining to the alluvial or colluvial sediments of the former Lake Bonneville. Winter snowfall and spring snowmelt control annual hydroclimate, rapid population growth dominates geographic change in lower elevations and urban gas and particle emissions contribute to episodes of severe air pollution in this closed‐basin. Long‐term hydroclimate observations across this diverse landscape provide the foundation for an expanding network of infrastructure in both montane and urban landscapes. Current infrastructure supports both basic and applied research in atmospheric chemistry, biogeochemistry, climate, ecology, hydrology, meteorology, resource management and urban redesign that is augmented through strong partnerships with cooperating agencies. These features allow WEO to serve as a unique natural laboratory for addressing research questions facing seasonally snow‐covered, semi‐arid regions in a rapidly changing world and an excellent facility for providing student education and research training.

Follstad Shah, Jennifer J.↗
DOE OSTI2022

Conservative solute transport processes and associated transient storage mechanisms: Comparing streams with contrasting channel morphologies, land use and land cover

Abstract Land use within a watershed impacts stream channel morphology and hydrology and, therefore, in‐stream solute transport processes and associated transient storage mechanisms. This study evaluated transport processes in two contrasting stream sites where channel morphology was influenced by the surrounding land use, land cover, climate and geologic controls: Como Creek, CO, a relatively undisturbed, high gradient, forested stream with a gravel bed and complex channel morphology, and Clear Creek, IA, an incised, low‐gradient stream with low‐permeability substrate draining an agricultural landscape. We performed conservative stream tracer injections at these sites to address the following questions: (1) How does solute transport vary between streams with differing morphologies? and (2) How does solute transport at each stream site change as a function of discharge? We analysed in‐stream tracer time series data and compared results quantifying solute attenuation in surface and subsurface transient storage zones. Significant trends were observed in these metrics with varying discharge conditions at the forested site but not at the agricultural site. There was a broad range of transport mechanisms and evidence of substantial exchange with both surface and hyporheic transient storage in the relatively undisturbed, forested stream. Changing discharge conditions activated or deactivated different solute transport mechanisms in the forested site and greatly impacted advective travel time. Conversely in the simplified agricultural stream, there was a narrow range of solute transport behaviour across flows and predominantly surface transient storage at all measured discharge conditions. These results demonstrate how channel simplification inhibits available solute transport mechanisms across varying discharge conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Temporal Patterns in Soil Redox Potential Vary Across a Freshwater Coastal Delta

Widespread and persistent flooding is submerging coastal ecosystems, particularly along the Louisiana Gulf coast, where flooding results from the compounding effects of ground subsidence and rising sea level. Restoration projects aim to mitigate land loss by diverting sediment loads from rivers into degraded areas to increase ground elevation. To predict how coastal ecosystems will change over time in response to projected changes in relative sea level and restoration, it is necessary to understand how subsurface biogeochemical processes respond to dynamic hydrologic forcings. Here, this study evaluates how environmental parameters that integrate biogeochemical processes vary with water table fluctuations in the freshwater Wax Lake Delta (WLD) in Louisiana, USA, where water diversions have formed one of the only active deltas along the coast. High‐frequency observations of water level, soil redox potential, specific conductance and pH were made for 1 year along elevation transects located on the older, proximal and younger, distal ends of a deltaic island. Redox responded rapidly to changing water tables, with fluctuations occurring primarily in shallow soils (< 20 cm) and at higher elevations. Deeper soils and those at lower elevation remained inundated and reduced. Semi‐diurnal tidal fluctuations were pronounced in younger, distal soils, presumably due to rapid groundwater exchange with the river channel. Tidal signals were muted in older soils that instead exhibited seasonal variability associated with river discharge and evapotranspiration. Although much of the delta sediments are persistently reducing and anoxic, redox fluctuations in the natural levees that border the deltaic islands likely drive high rates of biogeochemical activity. Evaluating how hydrology drives the frequency and duration of redox fluctuations provides a basis for understanding how biogeochemical processes might vary with complex hydrological interactions in coastal systems.

Herndon, Elizabeth [Oak Ridge National Laboratory ↗
DOE OSTI2025

Phosphorus and cover crop management practices affect phosphorus speciation in soils and eroded sediments

Abstract Agricultural runoff often contains P in dissolved and sediment‐bound forms, decreasing surface water quality. No‐till and cover cropping conservation practices have been recommended for reducing erosion and nutrient loss from cropping systems. The overall aims of this study were to characterize and evaluate the effects of fertilizer (placement and source) and cover crop management on P speciation in surface runoff sediments and source soil. In 2014, a field‐scale experiment was established in a no‐till, corn (Zea maysL.)–soybean (Glycine maxL.) cropping system with two cover crop treatments (with and without a winter crop; winter wheat [Triticum aestivumL.], rapeseed [Brassica napusL.], hairy vetch [Vicia villosaRoth], winter triticale [×Triticosecale Wittm.], and cereal rye [Secale cerealeL.]) and three P fertilizer management treatments (no P, fall broadcast diammonium phosphate, and spring subsurface injected ammonium polyphosphate). Phosphorus fractionation in the source soil collected in the fall of 2019 and sediment samples collected throughout 2020 were analyzed using a modified sequential P extraction method to evaluate the cumulative effects of imposing the treatment factors over 5 years. The direct P speciation was done using X‐ray absorption near edge structure spectroscopy. The indirect P speciation (fractionation) results showed that the management practices affected the exchangeable, organic matter‐associated, and Fe‐bound P fractions in sediments and the exchangeable and residual fractions in source soil. Direct P speciation results showed a depletion of Fe‐associated P in soil and sediment from cover crop treatment, suggesting that Fe‐associated P species were affected by cover crops. Changes in soil and runoff sediment P speciation would change the proportions and forms of soluble and particulate P in runoff sediments and may influence P bioavailability in aquatic ecosystems. Developing P fertilizer and cropping system management options with an understanding of soil P transformations helps maintain environmental sustainability.

Environmental Sciences & Ecology↗
DOE OSTI2022

Influence of copper and aluminum substitution on high-temperature oxidation of the FeCoCrNiMn “Cantor” alloy

In this study, the oxidation behavior of FeCoCrNiMn (HEA + Mn) is compared to three modified HEAs manufactured by substituting Mn with Al, Cu, or Al + Cu. Oxidation tests were conducted between 600°C and 800°C for up to 500 h in synthetic air. Substitution of Mn leads to a significant improvement in the oxidation resistance for the three modified HEAs. For FeCoCrNiCu (HEA + Cu), a local attack of a Cu-rich phase was observed, leading to the formation of CuO blisters on the surface. The FeCoCrNiAl (HEA + Al) alloy was characterized by the formation of a thin Al 2 O 3 surface layer for all temperatures. However, for the HEA + Al alloy the formation of AlN was observed after 300 h at 800°C, leading to a partial breakdown of the protective scale. FeCoCrNiCuAl (HEA + Cu + Al) by far showed the best oxidation resistance, characterized by the formation of a highly protective Al 2 O 3 scale that effectively inhibited nitrogen penetration into the metal subsurface and local attack of the Cu-rich phase.

36 MATERIALS SCIENCE↗
DOE OSTI2020

An anisotropic viscoplasticity model for shale based on layered microstructure homogenization

Viscoplastic deformation of shale is frequently observed in many subsurface applications. Many studies have suggested that this viscoplastic behavior is anisotropic—specifically, transversely isotropic—and closely linked to the layered composite structure at the microscale. In this work, we develop a two-scale constitutive model for shale in which anisotropic viscoplastic behavior naturally emerges from semianalytical homogenization of a bilayer microstructure. The microstructure is modeled as a composite of soft layers, representing a ductile matrix formed by clay and organics, and hard layers, corresponding to a brittle matrix composed of stiff minerals. This layered microstructure renders the macroscopic behavior anisotropic, even when the individual layers are modeled with isotropic constitutive laws. Using a common correlation between clay and organic content and magnitude of creep, we apply a viscoplastic modified Cam-Clay plasticity model to the soft layers, while treating the hard layers as a linear elastic material to minimize the number of calibration parameters. We then describe the implementation of the proposed model in a standard material update subroutine. The model is validated with laboratory creep data on samples from three gas shale formations. We also demonstrate the computational behavior of the proposed model through simulation of time-dependent borehole closure in a shale formation with different bedding plane directions.

58 GEOSCIENCES↗
DOE OSTI2021

Randomized algorithms for generalized singular value decomposition with application to sensitivity analysis

The generalized singular value decomposition (GSVD) is a valuable tool that has many applications in computational science. However, computing the GSVD for large-scale problems is challenging. Motivated by applications in hyper-differential sensitivity analysis (HDSA), in this work we propose new randomized algorithms for computing the GSVD which use randomized subspace iteration and weighted QR factorization. Detailed error analysis is given which provides insight into the accuracy of the algorithms and the choice of the algorithmic parameters. We demonstrate the performance of our algorithms on test matrices and a large-scale model problem where HDSA is used to study subsurface flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Convergence analysis of single rate and multirate fixed stress split iterative coupling schemes in heterogeneous poroelastic media

Recently, the accurate modeling of flow–structure interactions has gained more attention and importance for both petroleum and environmental engineering applications. Of particular interest is the coupling between subsurface flow and reservoir geomechanics. Different single rate and multirate iterative and explicit coupling schemes have been proposed and analyzed in the past. In addition, Banach fixed point contraction results were obtained for iterative coupling schemes, and conditionally stable results were obtained for explicit coupling schemes. In this work, we will consider the mathematical analysis of the single rate and multirate fixed stress split iterative coupling schemes for spatially heterogeneous poroelastic media. We will re–establish the contractivity for both schemes in the localized case, and we will show that heterogeneities come at the expense of imposing more restricted conditions on the number of fine flow time steps that can be taken within one coarse mechanics time step in the multirate case. Our mathematical analysis is supplemented by numerical simulations validating our derived upper bounds. Finally, to the best of our knowledge, this is the first rigorous mathematical analysis of the multirate fixed–stress split iterative coupling scheme in heterogeneous poroelastic media.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Atomic Dynamics of Multi‐Interfacial Migration and Transformations

Redox-induced interconversions of metal oxidation states typically result in multiple phase boundaries that separate chemically and structurally distinct oxides and suboxides. Directly probing such multi-interfacial reactions is challenging because of the difficulty in simultaneously resolving the multiple reaction fronts at the atomic scale. Using the example of CuO reduction in H 2 gas, a reaction pathway of CuO → monoclinic m-Cu 4 O 3 → Cu 2 O is demonstrated and identifies interfacial reaction fronts at the atomic scale, where the Cu 2 O/m-Cu 4 O 3 interface shows a diffuse-type interfacial transformation; while the lateral flow of interfacial ledges appears to control the m-Cu 4 O 3 /CuO transformation. Together with atomistic modeling, it is shown that such a multi-interface transformation results from the surface-reaction-induced formation of oxygen vacancies that diffuse into deeper atomic layers, thereby resulting in the formation of the lower oxides of Cu 2 O and m-Cu 4 O 3 , and activate the interfacial transformations. In conclusion, these results demonstrate the lively dynamics at the reaction fronts of the multiple interfaces and have substantial implications for controlling the microstructure and interphase boundaries by coupling the interplay between the surface reaction dynamics and the resulting mass transport and phase evolution in the subsurface and bulk.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Exploring Charged Defects in Ferroelectrics by the Switching Spectroscopy Piezoresponse Force Microscopy

Monitoring the charged defect concentration at the nanoscale is of critical importance for both the fundamental science and applications of ferroelectrics. However, up-to-date, high-resolution study methods for the investigation of structural defects, such as transmission electron microscopy, X-ray tomography, etc., are expensive and demand complicated sample preparation. Here, with an example of the lanthanum-doped bismuth ferrite ceramics, a novel method is proposed based on the switching spectroscopy piezoresponse force microscopy (SSPFM) that allows probing the electric potential from buried subsurface charged defects in the ferroelectric materials with a nanometer-scale spatial resolution. When compared with the composition-sensitive methods, such as neutron diffraction, X-ray photoelectron spectroscopy, and local time-of-flight secondary ion mass spectrometry, the SSPFM sensitivity to the variation of the electric potential from the charged defects is shown to be equivalent to less than 0.3 at% of the defect concentration. Additionally, the possibility to locally evaluate dynamics of the polarization screening caused by the charged defects is demonstrated, which is of significant interest for further understanding defect-mediated processes in ferroelectrics.

42 ENGINEERING↗
DOE OSTI2021

Long-Term Degradation of Passivated Emitter and Rear Contact Silicon Solar Cell under Light and Heat

Advanced designs enable high-efficiency solar cells; however, more complex structures create new long-term stability concerns. Herein, the long-term degradation processes affecting advanced silicon solar cells using laboratory-based illumination and heating over hundreds of hours are investigated. The activation energy for the degradation of voltage is estimated and the degradation rates to normal solar cell operating temperature ranges are extrapolated. The cell degradation observed at high temperatures in the lab is kinetically similar to the process affecting field-deployed modules contributing to 0.37% year-1 of annualized degradation. Electroluminescence and photoluminescence mapping show that the degradation is dominated by minority carrier lifetime reduction. Suns-open-circuit voltage and light beam-induced current results indicate that the degradation could result from passivation degradation at the surface or defect formation in the near-subsurface region, leading to increased minority carrier recombination. This work highlights a long-term degradation process under elevated temperature and illumination that may continue to affect cells in an irreversible manner that is separate from recoverable light-induced degradation and light- and elevated temperature-induced degradation.

14 SOLAR ENERGY↗
DOE OSTI2020

Algebraic expressions for estimating the impact depths of a surface barrier over a homogeneous soil

Engineered surface barriers are used to isolate subsurface contaminants for effective long-term containment of municipal solid waste, other nonhazardous solid and liquid waste, hazardous and toxic wastes, and radioactive waste. The impact depths of a surface barrier are affected by the pre-barrier recharge rate and the properties of the soil beneath the barrier. In this paper, the pore-size-specific (PSS) water velocity is defined and an algebraic expression of PSS velocity is derived based on the stream tube concept and the Brooks and Corey hydraulic retention model. Algebraic expressions are developed to estimate drainage velocities and barrier impact depths after the emplacement of a surface barrier. Four impact depth terms are used to convey the protective effect: drainage front, average drainage, the location with 50% impact, and drainage tail. The drainage front depth is the deepest point at which the barrier has a detectable impact at a specific time (also called the near zero-impact depth). The average-impact depth is the depth at which average drainage occurs. At the 50% impact depth, the water flux rate is reduced by half because of the surface barrier. Lastly, the drainage tail depth (also called the full-impact depth) is the deepest depth at which the water conditions above it are in equilibrium with the barrier. The algebraic expressions show that the average-impact depth is no more than 1/3 of near zero-impact depth, while the 50% impact depth is slightly larger than 1/2 of the near zero-impact depth. The full-impact depth, depending on the final recharge rate from the surface barrier, is usually much smaller than the other impact depths. These differences lead to a very large transition zone beneath a surface barrier. Numerical simulations were conducted to replicate the same conditions. The numerical results corroborated the analytical models by predicting very similar water content profiles and near zero-, average-, 50%, and full-impact depths. The algebraic expressions provided in this paper are useful for quickly identifying sites where the depth of the existing contaminants could be beyond the protection of a surface barrier.

54 ENVIRONMENTAL SCIENCES↗