Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Aquifer”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Evaluation of Zero-Net-Rate Pumping Tests

Accurately estimating the distribution of aquifer properties is key to understanding contaminant movement in the subsurface. The distribution of aquifer properties is typically addressed using slug or constant-rate well tests, and the pros and cons of these tests are well known. Slug tests are appealing because they avoid removing contaminated water, but their results are affected by well skin and the small volume of displaced water limits the volume of aquifer that can be evaluated. Constant-rate well tests have the disadvantage of requiring disposal of potentially contaminated water, but they can generate properties that are more representative than slug tests, and they can be used to estimate well efficiency and storativity, which are difficult to characterize using slug tests. Periodic pumping tests are appealing because they have many of the advantages and few of the disadvantages of slug and constant-rate well tests. Periodic pumping tests involve cycling the pumping rate with a regular period and measuring the resulting response in monitoring wells. Some of these tests involve imposing a periodic rate on a constant mean rate. However, other tests involve moving water out and back into the well at rates that are balanced so the net rate after each period is zero. Some zero-net-rate (ZNR) well tests use rates that follow a sinusoidal pattern, whereas others use a square wave pattern that switches between a constant rate of pumping and a constant rate of injection to achieve a zero-net rate. ZNR well tests appear to be a useful compromise between slug and constant rate tests, but methods for conducting and analyzing the results from these tests have received limited evaluation. Periodic pumping tests are typically evaluated through recorded pressure responses to a disturbance source. Recent work has demonstrated that the strain field in the vadose zone shows a response to disturbances in the underlying aquifer. Various instrument designs are available that can record vertical and horizontal strain to very precise resolutions (10-9 ε). Measuring the strain in the vadose zone can present lower costs than measuring pressure in a monitoring well, so using strain could improve the resolution of ZNR tests, as well as make these tests potentially cheaper. A ZNR periodic pumping system was constructed and used to perform tests in the Clemson, South Carolina area. The system was designed to pump and inject with a periodic square-wave, so pumping occurs at a constant rate for half the period and is followed by injection at the same rate for half the period resulting in no net gain or loss of water from the aquifer. The system is designed to generate flow rates from 1 to 3.5 gallons per minute (gpm) with a capacity of 875-gallon per half period. Pressure and strain monitoring points are located in the vicinity of the pumping well. A series of 10 ZNR periodic pumping tests were conducted at the field site, with periods ranging from 2 to 540 minutes. Pressure data were collected and recorded in four monitoring wells in various locations around the pumping well. In addition, vertical and horizontal strain and tilt data were recorded in various locations around the pumping well. Traditional aquifer tests, including slug tests and constant-rate pumping tests were conducted at the field site to provide a baseline of aquifer parameter estimates that are compared the results of ZNR tests. Pressure data measured at monitoring wells during ZNR tests were analyzed utilizing an analytical solution (Streltsova, 1988) that assumes a confined aquifer and estimates hydraulic diffusivity and transmissivity using the time lag and amplitude of the pressure signal. The storativity can be separated out from the hydraulic diffusivity using the transmissivity estimate. The time-lag and amplitude of the pressure signals were estimated using a Fourier transform. The time-lag of the primary period of the pressure increases as a linear function of distance from the pumping well, and it increases as the square root of the period of the ZNR test over a range of periods spanning two orders of magnitudes. The first few harmonics follow similar trends, but the higher frequency harmonics diverge from this trend. Strain and tilt data from various instruments in the vadose and saturated zones are periodic during ZNR tests. The time lags of the strain components increase roughly linearly with distance from the well. Many of the time lags of the strain data from the vadose zone are similar to the time lags observed in the pressure data from a similar distance and pumping period. Theoretical experiments were performed to understand how pressure and strain responds to ZNR periodic pumping. The time-lags of the simulated pressure responses match the field data trends. The results were used to validate the Streltsova (1988) solution for periodic pumping to estimate hydraulic diffusivity. Hydraulic diffusivity was estimated to be 1.8 x 10-2 m2s-1 < Dh < 9.0 x 10-2 m2s-1 using pressure data from all the wells during ZNR tests. Transmissivity was estimated to be 0.8 x 10-4 m2s-1 < T < 4.3 x 10-4 m2s-1. Assuming Storativity = T/Dh gives 1.1 x 10-3 < S < 19 x10-3. Transmissivity is used to estimate the hydraulic conductivity, K, by dividing by the assumed aquifer thickness. Two constant-rate pumping tests were conducted and analyzed using two conceptual models: The Hantush (1961) solution was used to analyze data assuming confined conditions and the Neuman (1974) solution was used to analyze data assuming unconfined conditions. The Hantush (1961) solution estimated transmissivity to be 1.6 x 10-4 m2s-1 < T < 1.8 x 10-4 m2s-1 when using data from all the wells. Storativity was estimated to be 4.4 x 10-3 < S < 8.2 x 10-3. Assuming Dh = T/S gives 0.22 x 10-2 m2s-1 < Dh < 3.6 x 10-2 m2s-1. The Neuman (1974) solution estimated transmissivity to be 0.80 x 10-4 m2s-1 < T < 0.84 x 10-4 m2s-1. Total storativity (S + Sy) was estimated to be 22 x 10-3 < S < 110 x 10-3. Storativity excluding the Sy term was estimated to be 2.2 x 10-3 < S < 2.5 x 10-3. Assuming Dh = T/S gives 0.073 x 10-2 m2s-1 < S < 0.36 x 10-2 m2s-1. Two slug tests were performed on the pumping well, PW-2. The Bouwer and Rice (1976) solution for slug tests was used to analyze the data. This solution works for both confined and unconfined conditions, so both conceptual models were covered through the analysis. Hydraulic conductivity was estimated to be 1.95 x 10-6 ms-1 < K < 2.49 x 10-6 ms-1. Statistical comparison show generally no significant difference in the estimates of Dh, K, or S made using ZNR and conventional pumping tests.. However, the total storage estimated by the Neuman (1974) solution for unconfined settings is larger than that estimate using the Streltsove analysis applied to the ZNR tests. This is likely because the Neuman solution considers delayed yield from storage that the water table, whereas the Streltsova (1988) solution assumes confined conditions. The analysis also indicates that there is a statistically significant difference between parameters estimated with conventional slug tests and those measured with either constant-rate pumping tests or ZNR tests. Hydraulic diffusivity is estimated from the time-lag and the distance from the well. The time-lags of the strain data from the vadose zone are similar to the lags of the pressure data from a similar distance. This indicates that hydraulic diffusivities estimated from the strain data measured in the vadose zone would be similar to the estimates from the pressure data measured in the aquifer. There appear to be sign reversals in some of the strain data that were corrected to estimate the time lag. The dissolved oxygen (CO) concentration was measured during several ZNR tests to evaluate the feasibility of using the procedure to increase contaminant degradation kinetics that are related to DO. The results indicate that DO can be increased during the cycling of ZNR tests in some cases.

Smith-Jones, Austin↗
DOE OSTI2020

Interplay of changing irrigation technologies and water reuse: example from the upper Snake River basin, Idaho, USA

Abstract. Careful allotment of water resources for irrigation is critical for ensuring the resiliency of agriculture in semiarid regions, and modernizing irrigation technology to minimize inefficient water losses is an important tool for farmers and agricultural economies. While modernizing irrigation technology can achieve reductions in the nonbeneficial use of water, such as bare soil evaporation and nonconsumptive losses, water returned to the landscape is also reduced, often eliminating flow paths that other users rely on. In basins using a combination of surface and groundwater, replenishing aquifer storage by the managed aquifer recharge (MAR) of seasonally available water can mitigate the aquifer drawdown that results from reduced recharge when irrigation efficiency is improved. We examine the effects of MAR on the system-scale efficiency of modernizing irrigation technology and the resulting changes in the reuse of nonconsumptive losses, using a macroscale hydrologic model applied to the semiarid upper Snake River basin (USRB) of western Wyoming and southern Idaho, USA. Irrigation technologies were represented explicitly in the model, and available data informed baseline parameterizations of the irrigation technology. A suite of parameterizations were simulated that updated the existing technologies to be more efficient, both with and without sufficient MAR to cause a stabilization of the aquifer at the present-day head. As expected, simulated changes in irrigation technology resulted in greater downstream export of pristine water and a higher rate of aquifer drawdown when MAR was not simulated. Under current water use and cropping patterns, we were not able to simulate aquifer stabilization and maintain discharge downstream at any level of irrigation efficiency. We found support for the hypothesis that, as efficiency improves, less MAR is required to maintain a stable aquifer than when return flows are reduced due to increased efficiency. To evaluate the hypothesis, we defined the management benefit as a metric that compared the difference between the change in irrigation's net recharge and the change in MAR required as irrigation technology became more efficient. The metric generally indicated that less MAR was needed than net recharge was lost, but only for the most efficient case did the management benefit exceed the MAR needed at the baseline to stabilize the aquifer. Increasing efficiency of irrigation technology reduced the reuse of the gross irrigation derived from prior nonconsumptive losses, but simulating MAR increased reuse for a given parameterization, leading to higher effective irrigation efficiency. We find that local groundwater storage that users depend on is generally more sensitive to management decisions than downstream flows, and the drawdown of the aquifer without MAR always exceeded any decrease in discharge induced by MAR. Improving resource sufficiency in semiarid systems like the USRB will require an array of solutions that will need to balance benefits to local and downstream users.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Post-remediation geophysical assessment: Investigating long-term electrical geophysical signatures resulting from bioremediation at a chlorinated solvent contaminated site

There is a growing need to assess long-term impacts of active remediation strategies on treated aquifers. A variety of biogeochemical alterations can result from interactions of the amendment with the aquifer, conceivably leading to a geophysical footprint of the long-term alteration of an aquifer. This concept of post-remediation geophysical assessment was investigated in a shallow, chlorinated solvent-contaminated aquifer six to eight years after amendment delivery. Surface resistivity imaging and cross-borehole resistivity and induced polarization (IP) imaging were performed on a transect that spanned treated and untreated zones of the aquifer. Established relationships between IP parameters and surface electrical conductivity were used to predict vertical profiles of electrolytic conductivity and surface conductivity from the inverted cross-borehole images. Aqueous geochemistry data, along with natural gamma and magnetic susceptibility logs, were used to constrain the interpretation. The electrical conductivity structure determined from surface and borehole imaging was foremost controlled by the electrolytic conductivity of the interconnected pore space, being linearly related to fluid specific conductance. The electrolytic conductivity (and thus the conductivity images alone) did not discriminate between treated and untreated zones of the aquifer. In contrast, inverted phase angles and surface conductivities did discriminate between treated and untreated zones of the aquifer, with the treated zone being up to an order of magnitude more polarizable in places. Supporting aqueous chemistry and borehole logging datasets indicate that this geophysical footprint of the long-term impact of the remediation on the aquifer is most likely associated with the formation of polarizable, dispersed iron sulfide minerals. The study opens the door to the possibility of employing time-lapse electrical geophysical measurements to assist with long-term environmental stewardship of legacy sites undergoing active remediation.

bioremediation↗
DOE OSTI2023

Using 81 Kr and isotopic tracers to characterise old groundwater in the Bangkok metropolitan and vicinity areas

Krypton-81 was applied to investigate the age of groundwater in the aquifer system in the Bangkok metropolitan and vicinity areas. Stable ( 2 H, 18 O and 13 C) and radioactive ( 3 H, 85 Kr and 14 C) isotopes and noble gases were applied in parallel. Low levels of 14 C and significant radiogenic 4 He confirm that groundwater in the deep aquifers is older than 30 ka. 81 Kr analysis identified groundwater with ages ranging from 17 to 300 ka. At some sites, large age discrepancies between 81 Kr and 14 C indicated that inter-aquifer mixing is likely occurring. The interpretation of the noble gases suggests that groundwaters in the deeper aquifers, with apparent ages of 300 to 10 ka, have recharged in slightly colder and wetter climates than those found in the upper aquifers with apparent ages < 10 ka. Degradation of water quality from seawater intrusion was identified in the upper four aquifers. This was also evidenced by higher δ 18 O and δ 2 H values, typical of seawater. The four deeper aquifers contain high quality water characterised by less enriched 18 O and 2 H. In conclusion, this work presents new findings of very old groundwater in the Bangkok aquifer system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Recent Developments in Hydrogeologic Applications for Strain Tensor Analyses

Changes in fluid pressure deform porous media and this effect occurs in a variety of hydrogeologic processes, from the change in storage during pumping or injection to fluctuations in water levels caused by barometric pressure. We have developed instruments for measuring small strains in porous media, and we have used the resulting strain data to evaluate well testing, hydraulic fracturing, manual loading at the ground surface, and ambient hydrologic processes, like rainfall and evaporation. A particularly important application is the use of strain tensor data measured at shallow depths to analyze well tests or hydraulic fractures conducted at much greater depths. An early demonstration of this technique was conducted at the North Avant Field north of Tulsa, Oklahoma, where Pennsylvanian sandstone creates a confined aquifer and oil reservoir at a depth of 530m. We have showed that the strains caused by injecting into the aquifer could be measured at a depth of 30m and used to evaluate the properties of the aquifer. We recently expanded the array of strainmeters at the North Avant Field by deploying three more instruments at shallow depth (30m) along with a deep instrument at 520m depth in the winter, 2021. The deep instrument is deployed in shale caprock slightly above the aquifer. To our knowledge, the deep strainmeter at the North Avant Field is the deepest strainmeter ever deployed and this required refining methods originally developed for shallow deployments. The instrument was lowered to depth on oil field tubing and cemented in place using techniques and materials developed for use in oil wells. Optical fiber used to communicate with the instrument was cemented in the annulus on the outside of tubing. This is significant because the techniques we used could readily be extended to greater depth, suggesting that strainmeters can be deployed over a wide range of depths for monitoring critical subsurface processes. For example, it suggests that strainmeters could be deployed through the caprock to monitor for leaks from underlying CO2 storage reservoirs. The strainmeter array at the North Avant Field has been used to characterize deformation during a series of injection tests in the spring and summer, 2021. All the new strainmeters respond to pumping, and the strainmeters we installed earlier also responded. To our knowledge, this is the first application of well testing in a deep aquifer that was monitored by an array of strainmeters—our earlier work used strainmeters at a single location. We are currently analyzing the strain data using an analytical solution, a proxy-based Bayesian inversion algorithm, and other methods. Strainmeter data has also been used to characterize periodic pumping tests by us and Riley Blais. A periodically varying pumping rate causes both hydraulic head and strain signals that vary with the same period as the pumping. The peaks and troughs of the head in monitoring wells lag behind the peaks and troughs of the head in the pumping well, and this lag time increases with distance from the pumping well. The lag time of the pressure and the distance to the monitoring well can be used in a simple analysis to estimate the hydraulic diffusivity of the aquifer. The lag time determined from strain data can be used to estimate aquifer properties using the same analysis that works for the pressure only for strain data measured at particular locations. That is because the strain field in a confining unit advances upward, laterally and then downward even though the pressure in the underlying aquifer only advances laterally, according to our recent simulations. We have field data showing that a small periodic signal superimposed on an injection rate at the North Avant Field will create a periodic strain signal at shallow strainmeters. The field data and the recent simulations suggest that including a periodic component to injection or pumping and then monitoring the resulting strain signal could be a way to monitor the subsurface.

Murdoch, Larry↗
DOE OSTI2023

Recent Developments in Hydrogeologic Applications for Strain Tensor Analyses

Changes in fluid pressure deform porous media and this effect occurs in a variety of hydrogeologic processes, from the change in storage during pumping or injection to fluctuations in water levels caused by barometric pressure. We have developed instruments for measuring small strains in porous media, and we have used the resulting strain data to evaluate well testing, hydraulic fracturing, manual loading at the ground surface, and ambient hydrologic processes, like rainfall and evaporation. A particularly important application is the use of strain tensor data measured at shallow depths to analyze well tests or hydraulic fractures conducted at much greater depths. An early demonstration of this technique was conducted at the North Avant Field north of Tulsa, Oklahoma, where Pennsylvanian sandstone creates a confined aquifer and oil reservoir at a depth of 530m. We have showed that the strains caused by injecting into the aquifer could be measured at a depth of 30m and used to evaluate the properties of the aquifer. We recently expanded the array of strainmeters at the North Avant Field by deploying three more instruments at shallow depth (30m) along with a deep instrument at 520m depth in the winter, 2021. The deep instrument is deployed in shale caprock slightly above the aquifer. To our knowledge, the deep strainmeter at the North Avant Field is the deepest strainmeter ever deployed and this required refining methods originally developed for shallow deployments. The instrument was lowered to depth on oil field tubing and cemented in place using techniques and materials developed for use in oil wells. Optical fiber used to communicate with the instrument was cemented in the annulus on the outside of tubing. This is significant because the techniques we used could readily be extended to greater depth, suggesting that strainmeters can be deployed over a wide range of depths for monitoring critical subsurface processes. For example, it suggests that strainmeters could be deployed through the caprock to monitor for leaks from underlying CO2 storage reservoirs. The strainmeter array at the North Avant Field has been used to characterize deformation during a series of injection tests in the spring and summer, 2021. All the new strainmeters respond to pumping, and the strainmeters we installed earlier also responded. To our knowledge, this is the first application of well testing in a deep aquifer that was monitored by an array of strainmeters—our earlier work used strainmeters at a single location. We are currently analyzing the strain data using an analytical solution, a proxy-based Bayesian inversion algorithm, and other methods. Strainmeter data has also been used to characterize periodic pumping tests by us and Riley Blais. A periodically varying pumping rate causes both hydraulic head and strain signals that vary with the same period as the pumping. The peaks and troughs of the head in monitoring wells lag behind the peaks and troughs of the head in the pumping well, and this lag time increases with distance from the pumping well. The lag time of the pressure and the distance to the monitoring well can be used in a simple analysis to estimate the hydraulic diffusivity of the aquifer. The lag time determined from strain data can be used to estimate aquifer properties using the same analysis that works for the pressure only for strain data measured at particular locations. That is because the strain field in a confining unit advances upward, laterally and then downward even though the pressure in the underlying aquifer only advances laterally, according to our recent simulations. We have field data showing that a small periodic signal superimposed on an injection rate at the North Avant Field will create a periodic strain signal at shallow strainmeters. The field data and the recent simulations suggest that including a periodic component to injection or pumping and then monitoring the resulting strain signal could be a way to monitor the subsurface.

Murdoch, Larry↗
DOE OSTI2024

Globally Gridded Groundwater Extraction Volumes and Costs under Six Depletion and Ponded Depth Targets

This repository contains simulated outputs from superwell – a hydro-economic tool for long-term assessment of groundwater cost and supply – providing globally gridded groundwater extractable volumes and associated unit costs ($/km³) for accessible groundwater production, based on a variety of user-defined depletion and ponded depth scenarios. Key model documentation: Niazi, H., Ferencz, S. B., Graham, N. T., Yoon, J., Wild, T. B., Hejazi, M., Watson, D. J., & Vernon, C. R. (2025). Long-term hydro-economic analysis tool for evaluating global groundwater cost and supply: Superwell v1.1. Geoscientific Model Development, 18(5), 1737-1767. https://doi.org/10.5194/gmd-18-1737-2025 Find the source code of the superwell model on GitHub: https://github.com/JGCRI/superwell Repository Overview Main output: superwell_outputs.7z contains 6 files (4.5 GB) named as superwell_py_deep_all_0.*PD_0.*DL.csv. These files present superwell outputs of global groundwater extraction volumes and cost estimates on a 0.5° scale for six scenarios with different Ponded Depth (PD; 0.3 and 0.6 m) and Depletion Limit (DL; 5%, 25%, and 40% of available volume) targets over the entire pumping lifetime of a grid cell superwell_py_deep_all_0.3PD_0.25DL_sample_100.csv contains superwell outputs for 100 data points sampled to match the global inputs' distribution superwell_py_deep_all_0.3PD_0.25DL_Grid_72548.csv contains superwell output for a single grid cell concept_v5.png provides an overview of the superwell workflow Outputs Description year_number: year of pumping depletion_limit: set depletion limit (DL) as a volume fraction of total available groundwater Mappings: continent, country, gcam_basin_id, Basin_long_name, grid_id: geographic identifiers and basin information Inputs: grid_area (km²): area of the grid cell whyclass: hydrogeological classification of the aquifer permeability (m/day), porosity (%), total_thickness (m), depth_to_water (m): aquifer properties. The geo-processed input data has been published separately: https://doi.org/10.57931/2307831 Model outputs: orig_aqfr_sat_thickness (m), aqfr_sat_thickness (m): original and remaining/instantaneous saturated thickness of the aquifer hydraulic_conductivity (m/day), transmissivity (m²/day): hydraulic properties of the aquifer radius_of_influence (m), areal_extent (km²): well radius and area of influence from the center of the well number_of_wells (-): number of wells in a grid cell determined by a ratio of well area and grid area max_drawdown (m), drawdown (m), drawdown_interference (m): well and aquifer drawdown during extraction total_head (m): total lift for the groundwater (depth to water plus drawdown) total_well_length (m): total depth of wells drilled well_yield (m³/day): pumping rate or well yield power (kW), energy (kWh): power and energy required for pumping groundwater Volume Outputs: volume_produced_perwell (m³), cumulative_vol_produced_perwell (m³): production volume metrics per well volume_produced_allwells (m³), cumulative_vol_produced_allwells (m³): aggregate extraction volumes for all wells in a grid cell available_volume (m³): available groundwater in storage for the grid cell as determined by aquifer properties depleted_vol_fraction: fraction of total volume pumped over available volumes in a grid cell (same as depletion limit) Cost Outputs: well_installation_cost ($): well installation cost based on the hydrogeological complexity of the aquifer annual_capital_cost, maintenance_cost, nonenergy_cost ($): nonenergy costs energy_cost_rate ($/kWh): electricity rate energy_cost ($): energy cost of pumping groundwater total_cost_perwell ($), total_cost_allwells ($): total annual energy and non-energy cost for each and all wells in a grid cell a unit_cost ($/m³), unit_cost_per_km3 ($/km³), unit_cost_per_acreft ($/acre-ft): total cost of pumping a unit of groundwater, indicated for different spatial units Key Resources Model documentation: Niazi, H., Ferencz, S., Graham, N., Yoon, J., Wild, T., Hejazi, M., Watson, D., & Vernon, C. (2024; In-prep). Long-term Hydro-economic Assessment Tool for Evaluating Global Groundwater Cost and Supply: Superwell v1. Geoscientific Model Development. Input data: Niazi, H., Watson, D., Hejazi, M., Yonkofski, C., Ferencz, S., Vernon, C., Graham, N., Wild, T., & Yoon, J. (2024). Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins. MSD-LIVE Data repository. https://doi.org/10.57931/2307831 superwell source code: https://github.com/JGCRI/superwell Cite as Niazi, H., Ferencz, S., Yoon, J., Graham, N., Wild, T., Hejazi, M., Watson, D., & Vernon, C. (2024). Globally Gridded Groundwater Extraction Volumes and Costs under Six Depletion and Ponded Depth Targets. MSD-LIVE Data repository. https://doi.org/10.57931/2307832 Contact Reach out to Hassan Niazi or Stephen Ferencz or open an issue in the superwell repository for questions or suggestions.

Earth Systems↗
NASA NTRS2016

Implications of Projected Climate Change for Groundwater Recharge in the Western United States

Existing studies on the impacts of climate change on groundwater recharge are either global or basin/ location-specific. The global studies lack the specificity to inform decision making, while the local studies do little to clarify potential changes over large regions (major river basins, states, or groups of states), a scale often important in the development of water policy. An analysis of the potential impact of climate change on groundwater recharge across the western United States (west of 100 degrees longitude) is presented synthesizing existing studies and applying current knowledge of recharge processes and amounts. Eight representative aquifers located across the region were evaluated. For each aquifer published recharge budget components were converted into four standard recharge mechanisms: diffuse, focused, irrigation, and mountain-systems recharge. Future changes in individual recharge mechanisms and total recharge were then estimated for each aquifer. Model-based studies of projected climate-change effects on recharge were available and utilized for half of the aquifers. For the remainder, forecasted changes in temperature and precipitation were logically propagated through each recharge mechanism producing qualitative estimates of direction of changes in recharge only (not magnitude). Several key patterns emerge from the analysis. First, the available estimates indicate average declines of 10-20% in total recharge across the southern aquifers, but with a wide range of uncertainty that includes no change. Second, the northern set of aquifers will likely incur little change to slight increases in total recharge. Third, mountain system recharge is expected to decline across much of the region due to decreased snowpack, with that impact lessening with higher elevation and latitude. Factors contributing the greatest uncertainty in the estimates include: (1) limited studies quantitatively coupling climate projections to recharge estimation methods using detailed, process-based numerical models; (2) a generally poor understanding of hydrologic flowpaths and processes in mountain systems; (3) difficulty predicting the response of focused recharge to potential changes in the frequency and intensity of extreme precipitation events; and (4) unconstrained feedbacks between climate, irrigation practices, and recharge in highly developed aquifer systems.

longitude↗
DOE OSTI2020

Analysis of Hydrology and Interim Measure Performance in the Chromium Plume at Los Alamos National Laboratory - 20498

Hexavalent chromium, Cr(VI), is present in the regional aquifer beneath Los Alamos National Laboratory (LANL) at concentrations greater than the New Mexico groundwater standard (50 ppb). The source of the Cr is blowdown from a power plant where potassium dichromate was used as a corrosion inhibitor from 1956 to 1972. Blowdown comingled with effluent from another source and migrated approximately 3 miles down a canyon, primarily as surface flow, until reaching an alluvial infiltration zone. The infiltrating water moved vertically through a thick (approximately 900 ft), complex stratigraphic sequence of unsaturated tuffs, basalt, and unconsolidated alluvial deposits. Horizontal transport also took place within spatially limited perched groundwater horizons in the vadose zone (VZ). Breakthrough at the regional aquifer water table is thought to have occurred in multiple hydraulic windows or 'drip points' located approximately 3 miles from the initial release site at the power plant. The regional aquifer is highly heterogeneous and anisotropic, and displays unconfined behavior in an upper zone and confined behavior at depth. The Cr(VI) contamination has so far remained confined to the upper approximately 25 m of the aquifer. The ambient hydraulic gradient in the regional aquifer is toward the east/southeast, but the water table is extremely flat in the Cr plume area compared to surrounding areas. There is some local variability in flow direction, potentially related to effects of seasonal pumping from nearby Los Alamos county water-supply wells, which draw from the regional aquifer several hundred feet below the water table. An interim measure (IM) is currently operating at the site to achieve hydraulic control of plume migration using a series of extraction wells, ion-exchange treatment, and injection along the downgradient portion of the plume. The migration of Cr contamination in the aquifer has been previously evaluated using a calibrated probabilistic numerical model, which has also been used to design the spatial configuration of extraction and injection wells for the IM, and to identify locations for monitoring wells. The model is built using the Finite Element Heat and Mass (FEHM) transfer code, and the Model Analysis and Decision Support (MADS) software is used for analyses. To date, the portion of the IM system that has operated in a continuous manner has been very effective at hydraulic control. Modifications to the flow field associated with the influence of IM pumping and injection are evaluated here, with modeling constraints to be provided by recent observations of tracers injected into both monitoring wells and IM injection wells, as well as opportunistic observations of geochemical signatures of treated ion-exchange effluent injected into IM injection wells. The tracer and geochemical data provide insights into modified flow directions within the Cr plume area that cannot be otherwise obtained. Results of the modeling work are being used to inform adaptive management of IM operations and will be a valuable tool for development and refinement of long-term remediation strategies. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2022

Analyzing Recharge Dynamics and Storage in a Thick, Karstic Vadose Zone

Karst vadose zone heterogeneity creates complex transmission and storage dynamics that affect the timing and magnitude of aquifer recharge. Young, high-matrix permeability eogenetic karst aquifers may have significantly higher matrix storage than older, lower matrix permeability counterparts. In vulnerable and water-limited karst regions, the timescales of storage may be important for seasonal and sub-seasonal water resource management. We create a framework to quantify storage dynamics in karst aquifers using high-resolution precipitation and groundwater levels from the Northern Guam Lens Aquifer in the US territory of Guam. We estimate recharge using the Water Table Fluctuation method, and then develop transfer functions between precipitation and recharge to quantify storage and release of water from the vadose zone. The transfer functions are partitioned into different flow pathways including conduit, conduit/matrix, and slowly draining matrix. Probability distributions are fit to each pathway to determine the average travel times of infiltrated waters. The results show that aquifer recharge through secondary porosity features typically occurs within a few hours of a rainfall event, but this rapid recharge accounts for only 12%–28% of total recharge. The majority of aquifer recharge (>70%) occurs within a month of a contributing storm event. An additional 10% of recharge, on average, took longer than a month to reach the water table. The framework established can help improve the hydrological modeling and freshwater management for karst aquifers.

Patricia Spellman↗
DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗
NASA NTRS2023

Large-Scale Groundwater Monitoring in Brazil Assisted With Satellite-Based Artificial Intelligence Techniques

Here, we develop and test an artificial intelligence (AI)-based approach to monitor major Brazilian aquifers. The approach combines Gravity Recovery and Climate Experiment (GRACE) data and ground-based hydrogeological measurements from Brazil’s Integrated Groundwater Monitoring Network at hundreds of wells distributed in twelve aquifers across the country. We tested model ensembles based on three AI approaches: Extreme Gradient Boost, Light Gradient Boosting Model and CatBoost, followed by a Linear Regression (LR) step. The approach is further boosted with wavelet and seasonal decomposition processes applied to GRACE data. To determine the AI-based model’s sensitivity to data availability, we propose four experiments combining hydrogeological measurements from different aquifers. Groundwater storage estimates from the Global Land Data Assimilation System (GLDAS) are used as benchmark. A sensitivity analysis shows that the LR-based model ensemble is the best suited and to reproduce groundwater storage change in all studied Brazilian aquifers. Results show that the proposed approach outperforms GLDAS in all experiments, with an RMSE value of 2.68cm for the experiment that covers all monitored wells in Brazil. GLDAS resulted in RMSE=6.76cm. Using our AI model outputs, we quantified the groundwater storage change of two major aquifers, Urucuia and Bauru-Caiuá, over the past two decades: -31km 3 and -6km 3 , respectively. Water loss is driven by a prolonged drought across most of the country and intensification of groundwater pumping for irrigation. This study demonstrates that combining satellite data and AI can be a cost-effective alternative to monitor poorly equipped aquifers at the continental scale, with possible global replicability.

GRACE↗
DOE OSTI2023

SALSA_python (SALSython) v1

SALSA_python (Semi-Analytical Leakage Solutions for Aquifers) is a software that computed semi-analytical solutions for hydraulic head and brine leakage in multilayered aquifer–aquitard systems with geologic pressure forcing. It can simulate brine leakage into aquifers in a multi-aquifer-aquitard system with multiple injection and leaky wells. This situation is encountered in underground CO2 storage wherein brine leakage from pressurized reservoirs into aquifers is of concern. SALSA_python calls the original SALSA[1] subroutines in python by using the salsa2.so library. This enables the incorporation and coupling of SALSA computations into existing python-based codes and tools. SALSA_python is available for Linux and Mac operating systems. [1] Cihan, A., Oldenburg, C. M., & Birkholzer, J. T. (2022). Leakage from coexisting geologic forcing and injection-induced pressurization: A semi-analytical solution for multilayered aquifers with multiple wells. Water Resources Research, 58, e2022WR032343.

Bhuvankar, Pramod↗
DOE OSTI2023

Assessment of Baseline Monitoring Data for the Remote-Handled Low-Level Waste Disposal Facility at Idaho National Laboratory

The purpose of this document is to summarize environmental monitoring data collected at the Remote-Handled Low-Level Waste (RHLLW) Disposal Facility during the first four years of facility operations (FY 2019 through FY 2022). This summary provides a “baseline” condition for the facility, which can be used to distinguish contaminant releases from the RHLLW Disposal Facility from pre existing contamination as well as potential future releases from other sources (i.e., upgradient aquifer sources), as specified in the facility monitoring plan. Sufficient data has been collected over the first four years of operation of the RHLLW Disposal Facility to establish baseline conditions for future compliance monitoring of aquifer wells and performance monitoring of vadose zone lysimeters. Except for seven elevated tritium measurements from lysimeter HFEF-South believed to have been impacted by waste disposals, all data was deemed appropriate for establishing baseline concentrations. Monitoring of the aquifer detected indicator analytes gross alpha and gross beta and target analytes tritium and C-14. There were a few C-14 detections prior to an increase in the required detection level (RDL) after the first round of sampling. I-129 and Tc-99 were not detected above RDLs in aquifer samples. Lysimeter samples detected indicator analytes gross alpha and gross beta and target analyte tritium. Target analytes C-14, I-129, and Tc-99 were not detected above RDLs in lysimeter samples. Except for tritium, baseline concentrations of indicator and target analytes measured during the first four years of RHLLW facility operations are expected to represent conditions through the expected 20 year operating period of the facility. Elevated tritium in the aquifer, a result of past discharges of tritium at upgradient facilities, showed a decline over the first four years of operations, and levels are expected to continue to decline with time as a result of dilution and decay. Current tritium levels in lysimeter samples may be elevated due to tritium in the water applied during construction and infiltration testing of the facility. As a result, tritium concentrations in lysimeters may decline with time as a result of dilution and decay. Statistical measures (i.e., mean, standard deviation, and 99% upper confidence level [UCL]) were calculated for all detected analytes using the four years of concentration data. These measures will be used to evaluate future monitoring results and to demonstrate the facility is performing, or not performing, as established in the facility performance assessment (PA). Future monitoring results will also be combined with the baseline data in this report to further establish temporal trends and natural variability in the measurements. Based on the data collected during the baseline period, it is recommended the gross alpha action level for performance monitoring increase from 10 to 20 pCi/L due to several baseline measurements exceeding the preliminary action level of 10. An action level of 20 pCi/L is slightly greater than the 99% UCL, is protective of the aquifer, and would reduce unnecessary sampling. It is also recommended that tritium be added to the lysimeter analyte list for post-baseline period monitoring. Tritium, while not a dose concern, is a good tracer that can provide valuable information on water flow in and around the RHLLW Disposal Facility. Recommendations for a tritium action level are provided if considered necessary.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

baseflow: a MATLAB and GNU Octave package for baseflow recession analysis

baseflow is a MATLAB® toolbox designed for baseflow recession analysis, a technique used in hydrologic science to infer aquifer properties from streamflow. By leveraging widely available streamflow data, baseflow can be used to estimate aquifer properties such as hydraulic conductivity and drainable porosity over the modern instrumental stream gage record. The toolbox is intended for analysis of measured streamflow values recorded on a daily timestep, and is tailored for shallow, unconfined riparian aquifers that discharge groundwater laterally into adjacent streams. Additionally, baseflow can analyze the collective behavior of individual hillslope aquifers constituting hydrologic catchments, known as “watersheds”, from a nonlinear dynamical systems perspective. The toolbox incorporates recent advances in baseflow recession analysis to enable objective estimations of aquifer properties, and their sensitivity to methodological decisions, at both hillslope and catchment scales.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

CO 2 zonal injection rate allocation and plume extent evaluation through wellbore temperature analysis

Temperature analysis during a pause in injection operations, known as warmback analysis, has been used in the petroleum industry for evaluating the injection conformance and estimating the location of the flooded front in applications, such as waterflooding oil reservoirs. Here in this work, methods are introduced to extend the application of temperature warmback analysis to estimate the zonal CO 2 injection rate and zonal CO 2 plume extent during geologic CO 2 storage in a saline aquifer. First, novel analytical solutions are developed to model transient temperature in the aquifer during the injection and subsequent shut-in periods considering two-phase flow (gaseous CO 2 and aqueous brine) conditions in the aquifer. The solution involves a discretization of the aquifer into regions; the energy and mass conservation equations for the regions are solved simultaneously considering appropriate boundary conditions at the interfaces. Two solutions techniques are presented: multi-region and three-region solutions. Inverse models are developed accordingly to evaluate the injection profile and estimate the extent of the plume front in the reservoir during the injection period. The multi-region solution results in an inversion approach that requires regression analysis. However, the three-region formulation results in a simple graphical technique for inverse modeling. The analytical solutions are validated against a thermally coupled reservoir simulation tool using different synthetic cases for CO 2 injection in deep saline aquifers. The results of the developed solutions provide a good match with numerical results during forward and inverse modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Redox Heterogeneities Promote Thioarsenate Formation and Release into Groundwater from Low Arsenic Sediments

Groundwater contamination by As from natural and anthropogenic sources is a worldwide concern. Redox heterogeneities over space and time are common and can influence the molecular-level speciation of As, and thus, As release/retention but are largely unexplored. Here, we present results from a dual-domain column experiment, with natural organic-rich, fine-grained, and sulfidic sediments embedded as lenses (referred to as “reducing lenses”) within natural aquifer sand. We show that redox interfaces in sulfur-rich, alkaline aquifers may release concerning levels of As, even when sediment As concentration is low (<2 mg/kg), due to the formation of mobile thioarsenates at aqueous sulfide/Fe molar ratios <1. In our experiments, this behavior occurred in the aquifer sand between reducing lenses and was attributed to the spreading of sulfidic conditions and subsequent Fe reductive dissolution. In contrast, inside reducing lenses (and some locations in the aquifer) the aqueous sulfide/Fe molar ratios exceeded 1 and aqueous sulfide/As molar ratios exceeded 100, which partitioned As(III)–S to the solid phase (associated with organics or as realgar (As 4 S 4 )). These results highlight the importance of thioarsenates in natural sediments and indicate that redox interfaces and sediment heterogeneities could locally degrade groundwater quality, even in aquifers with unconcerning solid-phase As concentrations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

New Insights from Legacy Seismic Data regarding Basalt Elevations and Variability on the Hanford Site

Migration of groundwater contaminants in the Gable Gap area of the Hanford Site in southeastern Washington State is strongly influenced by the distribution and permeability of basalts that lie beneath an unconfined aquifer. Locally, folding and faulting of the Columbia River Basalt associated with the Yakima fold and thrust belt followed by erosion due to the Lake Missoula floods resulted in a complex basalt surface that represents either an impermeable lower boundary to the unconfined aquifer system or localized regions of increased permeability that potentially promote communication between the unconfined aquifer system and deeper, confined aquifer systems. Paleo-channels carved into the basalt by floodwaters are thought to provide preferential flow paths for groundwater contaminants. In 2011, a seismic landstreamer campaign was carried out to image the basalt surface and produced pre-stack depth migrated p-wave reflection images. The reflection images identified two large troughs that may represent paleo-channels and several areas of possible faulting. Here, the streamer data are re-analyzed using refraction travel-time and Rayleigh wave dispersion analyses to obtain images of compressional and shear wave velocities within the suprabasalt sediment sections and the upper basalt surface. The combined interpretation of reflection and seismic velocity images shows complexity in the basalt velocity and elevation, which varies by 50 m or more within the study area. These results, along with other ongoing geophysical investigations, will be used to inform the site geologic model and potentially guide placement of future boreholes needed to quantify vertical flow between the confined and unconfined aquifers.

58 GEOSCIENCES↗