Search NASA⌕ Search

SEARCH · Search NASA

Results for “Groundwater Modeling”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Data-Driven Method for Groundwater-Level Mapping and Monitoring-Well Network Optimization at Hanford

This report summarizes the initial results and outcomes of a physics-informed, data-driven groundwater level (GWL) mapping capability for the Hanford Site. GWL mapping at Hanford is typically conducted annually and requires a significant amount of computational and expert resources, and it does not allow assessment of the informational value of specific monitoring wells. The proposed method produces spatially and temporally resolved fields consistent with sparse, irregularly sampled, and nonuniformly distributed well measurements. Implemented successfully, this capability will allow rapid mapping of groundwater levels and provide an opportunity to optimize monitoring activities (both location and sampling frequency) based on data information value evaluation. The approach integrates a diffusion-based generative model – trained on MODFLOW simulation data from the Plateau-to-River (P2R) model – with score-based data assimilation (SDA), allowing observation-conditioned mapping without retraining for each monitoring-network layout.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2001

RFI to CMS: An Approach to Regulatory Acceptance of Site Remediation Technologies

Lockheed Martin made a smooth transition from RCRA Facility Investigation (RFI) at the National Aeronautics and Space Administrations'(NASA) Michoud Assembly Facility (MA-F) to its Corrective Measures Study (CMS) phase within the RCRA Corrective Action Process. We located trichloroethylene (TCE) contamination that resulted from the manufacture of the Apollo Program Saturn V rocket and the Space Shuttle External Tank, began the cleanup, and identified appropriate technologies for final remedies. This was accomplished by establishing a close working relationship with the state environmental regulatory agency through each step of the process, and resulted in receiving approvals for each of those steps. The agency has designated Lockheed Martin's management of the TCE-contamination at the MAF site as a model for other manufacturing sites in a similar situation. In February 1984, the Louisiana Department of Environmental Quality (LDEQ) issued a compliance order to begin the clean up of groundwater contaminated with TCE. In April 1984 Lockheed Martin began operating a groundwater recovery well to capture the TCE plume. The well not only removes contaminants, but also sustains an inward groundwater hydraulic gradient so that the potential offsite migration of the TCE plume is greatly diminished. This effort was successful, and for the agency to give orders and for a regulated industry to follow them is standard procedure, but this is a passive approach to solving environmental problems. The goal of the company thereafter was to take a leadership, proactive role and guide the MAF contamination clean up to its best conclusion at minimum time and lowest cost to NASA. To accomplish this goal, we have established a positive working relationship with LDEQ, involving them interactively in the implementation of advanced remedial activities at MAF as outlined in the following paragraphs.

Rowland, Martin A.↗
DOE OSTI2024

Gradient-based surface nuclear magnetic resonance for groundwater investigation

In medical magnetic resonance imaging, spatial localization (imaging) is based upon the application of controlled magnetic field gradients on top of the main magnetic field to spatially modulate the frequency and/or phase of the nuclear magnetic resonance (NMR) signal across the volume of investigation. In this work, we have applied similar physical principles to produce controlled magnetic field gradients during surface NMR-based groundwater investigations. In this approach, a gradient pulse of variable amplitude or duration is applied immediately after the excitation pulse to cause predictable phase encoding of the NMR signal as a function of depth. This approach is also applicable to emerging surface NMR detection methods that use a prepolarization field with fast nonadiabatic turn-off to generate detectable NMR signals from the shallow subsurface. In this case, the gradient pulse is applied after terminating the prepolarization field and provides a heretofore unavailable means of localizing the NMR response as a function of depth. The application of gradients can also be combined with tip-angle-based modulation to yield higher imaging resolution than can be achieved through either gradient- or tip-angle-based imaging alone. We implemented this new gradient-based capability into a surface NMR gradient generation accessory that is compatible with the GMR-Flex instrument and developed surface NMR-specific forward modeling and linear inverse models. We validated the accuracy of this novel gradient-based sNMR technology using computer simulations, experiments using a small pool filled with a discrete layer of bulk water, and field experiments at well-characterized groundwater test sites along Ebey Island, WA, and Larned, KS. The gradient-based sNMR imaging observations were compared with high-resolution direct push NMR results observed at these sites. The results of computer simulations and field experiments indicate improvements in both the detection (signal-to-noise ratio) and spatial resolution of shallow subsurface water content using gradient-based surface NMR, compared with traditional surface NMR imaging methods.

Geochemistry & Geophysics↗
DOE OSTI2020

Achieving the End State for the Pahute Mesa Corrective Action Units at the Nevada National Security Site - 20355

Underground nuclear testing at the Nevada National Security Site (NNSS) ended in 1992. To address the impact of radionuclide contamination from underground nuclear testing on the groundwater resources of Nevada, the U.S. Department of Energy (DOE) Environmental Management (EM) Nevada Program created the Underground Test Area (UGTA) activity to characterize the radionuclides in groundwater, understand the nature and extent of radionuclide migration, and forecast the distribution of radionuclides in groundwater for 1,000 years. The UGTA activity is regulated by the Federal Facility Agreement and Consent Order (FFACO), an agreement negotiated between the Nevada Division of Environmental Protection (NDEP), DOE, and the U.S. Department of Defense (DoD). DOE is close to achieving the end state (closure in place with monitoring and institutional control) for three of the five UGTA corrective action units (CAUs) on the NNSS. A number of challenges remain to achieve the end state of the other two CAUs, both located on Pahute Mesa. Pahute Mesa was the location of 82 underground nuclear tests, less than 10% of the total number of tests on the NNSS. Yet, Pahute Mesa contains slightly more than 60% of the total radionuclides (in curies). Based on groundwater sampling in monitoring wells, two radionuclide plumes have migrated several kilometers (km) in groundwater from selected test cavities and have crossed the boundary of the NNSS (yet remain within the boundaries of Federally controlled land). The path to achieve the end state continues to evolve as more data are collected and a better understanding of radionuclide migration in groundwater is developed. DOE has invested in drilling and sampling more than 50 characterization and monitoring groundwater wells over the past 25 years. The data indicate that the primary radionuclide of concern is tritium as it is about 89% of the total radionuclide inventory (in curies). As well, only tritium has been measured in groundwater outside of cavities at concentrations exceeding the Safe Drinking Water Act (SDWA) standard. For one of the tritium plumes that has migrated across the NNSS boundary, the average rate of migration has been measured as about 45 meters (m) per year with the rate of migration at the leading edge of the plume of about 85 meters per year. At that rate of migration, the tritium plume will decay to safe levels and not reach the publicly accessible environment at concentrations above the SDWA standard. Other radionuclides, at concentrations below the SDWA standards, will be monitored to ensure they remain at safe levels. Consequently, the end state path forward for Pahute Mesa has evolved to take full advantage of the data from the monitoring network to constrain uncertainty in model forecasts and to reduce reliance on probabilistic simulations. The focus of the end state approach relies on developing a monitoring well network that is protective of human health by increasing confidence that no radionuclide plume in the groundwater will migrate undetected to the accessible environment, located about 22 km from the nearest up-gradient underground nuclear test. Using the measured data to remove uncertainty, the evaluation of radionuclide migration from Pahute Mesa is directed toward meeting the goals of the end state. (authors)

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

Quantifying Subsurface Flow and Solute Transport in a Snowmelt‐Recharged Hillslope With Multiyear Water Balance

Abstract Quantifying flow and transport from hillslopes is vital for understanding water quantity and quality in rivers, but remains obscure because of limited subsurface measurements. Using measured hydraulic conductivity K profiles and water balance over a single year to calibrate a transmissivity feedback model for a hillslope in the East River watershed (Colorado) proved unsatisfactory for predicting flow over the subsequent years. Well‐constrained field‐scale K were obtained by optimizing subsurface flux predictions over years having large differences in recharge, and by including estimates of interannual transfer of excess snowmelt recharge. Water and solute exports during high snowmelt recharge occur predominantly via shallow groundwater flow through weathered rock and soil because of their enlarged transmissivities under saturated conditions. Conversely, these shallow pathways are less active in snow drought years when the water table remains deeper within the weathering zone. Hillslope soil water monitoring showed that rainfall does not infiltrate deeply during summer and fall months, and revealed water losses consistent with model ET predictions. By combining water table‐dependent fluxes with pore water chemistry in different zones, time‐dependent rates of solute exports become predictable. As an example, calibrated K were combined with dissolved nitrogen concentrations in pore waters to show the snowmelt‐dependence of reactive nitrogen exported from the hillslope, further supporting the recent finding that the weathering zone is the dominant source of reactive nitrogen at this site. Subsurface export predictions can now be obtained for wide ranges of recharge based on measurements of water table elevation and profiles of pore water chemistry.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Sensitivity of geophysical techniques for monitoring secondary CO 2 storage plumes

For geologic carbon storage, the ability to detect secondary CO 2 plumes—defined as those CO 2 plumes accumulating outside the intended storage reservoir—is fundamental to preventing unexpected CO 2 migration into groundwater resources and for risk and liability management. Understanding the sensitivity of various geophysical methods to secondary plumes is crucial for designing cost-effective monitoring schemes. We use several modeling scenarios to demonstrate the process of assessing sensitivities and detection thresholds of three primary geophysical techniques—surface seismic, borehole-to-surface electromagnetic (EM), and surface and borehole gravity—for early detection of secondary CO 2 plumes in the post-injection phase. While seismic reflection methods are often considered in monitoring strategies to track the evolution of CO 2 plumes, they are also the most expensive. Due to cost considerations, especially for long-term post-injection monitoring, other techniques complement seismic monitoring when designing an adaptive monitoring network. Borehole-to-surface EM or surface gravity surveys are feasible for time-lapse monitoring of deep secondary CO 2 plumes. Furthermore, these surveys could be carried at intervals defined by site-specific conditions. If time-lapse EM and/or gravity surveys detect any signal responses beyond the expected change, it would trigger a need for the higher resolution seismic survey.

58 GEOSCIENCES↗
DOE OSTI2020

FeS colloids – formation and mobilization pathways in natural waters

For this study, we have used synchrotron-based X-ray absorption spectroscopy (structure of Fe–S clusters), transmission electron microscopy (solid-phase crystallinity), Fourier-transform ion-cyclotron-resonance mass spectrometry (identity and composition of natural organic carbon compounds), inductively coupled plasma optical emission spectrometry (total aqueous Fe), and the revised ferrozine method (aqueous Fe( II ) and Fe( III ) concentrations) to determine the stability and nature of colloids generated by sulfidation of ferrihydrite nanoparticles in the absence and presence of organic compounds. We observed that reductive dissolution of ferrihydrite by aqueous sulfide generates nm-scale FeS clusters. Their subsequent aggregation, which promotes settling of FeS aggregates into the solid fraction, was directly correlated with sulfide/Fe ratio. At sulfide/Fe ratios ≤0.5, FeS clusters and larger colloids remained in suspension for at least 14 days (and up to several months). At sulfide/Fe ratios >0.5, sulfidation reaction rates were rapid and FeS cluster aggregation was accelerated. Moreover, the presence of organic compounds increased the time of suspension of FeS colloids, whereas increased ionic strength inhibited the generation of FeS colloids. We present a general conceptual model to predict when and where FeS colloids can form and enhance or inhibit the mobility of contaminants and nutrients associated with them. Our study indicates that in lowsalinity fresh groundwater systems poor in sulfate (i.e. low sulfidation; lakes, floodplains, peatlands etc.), the ferrihydrite sulfidation reaction generates aqueous FeS clusters and larger colloids that remain suspended over long time periods, thus mobilizing a substantial fraction of the total aqueous Fe and S.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2016

National Climate Assessment - Land Data Assimilation System (NCA-LDAS) Data and Services at NASA GES DISC

The National Climate Assessment-Land Data Assimilation System (NCA-LDAS) is an Integrated Terrestrial Water Analysis, and is one of NASAs contributions to the NCA of the United States. The NCA-LDAS has undergone extensive development, including multi-variate assimilation of remotely-sensed water states and anomalies as well as evaluation and verification studies, led by the Goddard Space Flight Centers Hydrological Sciences Laboratory (HSL). The resulting NCA-LDAS data have recently been released to the general public and include those from the Noah land-surface model (LSM) version 3.3 (Noah-3.3) and the Catchment LSM version Fortuna-2.5 (CLSM-F2.5). Standard LSM output variables including soil moistures temperatures, surface fluxes, snow cover depth, groundwater, and runoff are provided, as well as streamflow using a river routing system. The NCA-LDAS data are archived at and distributed by the NASA Goddard Earth Sciences Data and Information Services Center (GES DISC). The data can be accessed via HTTP, OPeNDAP, Mirador search and download, and NASA Earth data Search. To further facilitate access and use, the NCA-LDAS data are integrated into the NASA Giovanni, for quick visualization and analysis, and into the Data Rods system, for retrieval of time series of long time periods. The temporal and spatial resolutions of the NCA-LDAS data are, respectively, daily-averages and 0.125x0.125 degree, covering North America (25N 53N; 125W 67W) and the period January 1979 to December 2015. The data files are in self-describing, machine-independent, CF-compliant netCDF-4 format.

data assimiliation↗
DOE OSTI2026

Stochastic Thermo-Hydro Modeling and Neural Network Surrogate Development for Thermal Resource Assessment of the Galleries-to-Calories Geobattery

The Galleries-to-Calories Geobattery concept explores the use of abandoned coal mine workings for large-scale thermal energy transport and storage. The system involves injecting waste heat from a supercomputing facility into flooded mine galleries, where groundwater flow can store and transport thermal energy for potential recovery in downgradient district heating and cooling applications. To evaluate the feasibility and performance of the Geobattery under geological and operational uncertainty, we developed a suite of stochastic thermo-hydrological (TH) simulations using Monte Carlo sampling of key uncertain parameters (e.g., permeability, porosity, thermal conductivity, specific heat capacity) and operating conditions (e.g., injection rate, injection temperature). Results identified injection rate and temperature as the most influential parameters governing thermal front propagation, while the geometry of the room-and-pillar structure played a critical role in directing the extent and orientation of thermal advancement. Optimal combinations of material properties for maximizing heat recovery were also determined. To address the high computational cost of coupled-process stochastic modeling, we trained a neural network surrogate model on 24,000 physics-based realizations, achieving an R² > 0.99 and MAE < 0.1 for temperature predictions at monitoring locations. This surrogate enabled an additional 100,000 realizations for global sensitivity analysis and probabilistic thermal resource assessment. The integrated stochastic physics–surrogate modeling framework offers a computationally efficient tool for quantifying uncertainty, identifying key drivers, and informing early-stage design decisions for Geobattery systems.

15 - GEOTHERMAL ENERGY↗
DOE OSTI2025

Modeling of multi-component precipitation and crystallization for zero-liquid-discharge desalination

This study proposes novel ZLD treatment trains that integrate multi-component chemical precipitation and multi-effect evaporative crystallization for efficient brine management. The methodology emphasizes sustainability by integrating CO₂ for chemical precipitation and reducing environmental impact, exemplified with two case studies: produced water that includes industrial waste heat utilization by adopting the emerging vacuum air-gapped membrane distillation (VAGMD) technology, and brackish groundwater that is abundant in sulphate which undergoes treatment by low-salt-rejection reverse osmosis (LSRRO) with interstage chemical precipitation. Here, this study is the first of its kind to simultaneously account for reducing risk of mineral scaling and effective recovery of valuable solids when incorporating VAGMD and LSRRO in ZLD treatment trains. Using Reaktoro and WaterTAP, both Python-based, open-source platforms, we model the recovery of high-purity magnesium, calcium, and sodium salts while optimizing energy consumption and operational efficiency in proposed ZLD pathways for case studies of produced water and brackish desalination brine management. Validation with experimental and reference data confirms the reliability of the models used. In both case studies, the optimized ZLD process achieves recovery rates of 97%, 99%, and 90% for Mg, Ca, and Na, with purities exceeding 99%, and brine volume reduced to less than 4% of the initial feed flow.

3D CT scans↗
DOE OSTI2022

Incorporating Heterogeneity into 3D Geologic Models for the Alluvial Basin Reference Case

The unsaturated zone alluvium reference case is one of several geologic systems under consideration by the U.S. Department of Energy Office of Nuclear Energy for hosting repositories for spent nuclear fuel and associated waste (Sevougian et al., 2019). As noted by Mariner et al. (2018), the generic alluvial basin offers positive attributes that merit its consideration as a reference case by the Spent Fuel and Waste Science and Technology (SFWST) campaign. Hundreds of alluvial basins and sub-basins are scattered across the arid western United States (Figure 1-1). Precipitation and infiltration rates are relatively low with high evapotranspiration, resulting in deep regional water tables. Deep regional water tables allow for significant vertical separation between proposed repository depths and groundwater and thus longer transport paths from a hypothetical repository to an aquifer. Accumulations of alluvial sediments within these basins are typically on the order of hundreds of meters, and locally may exceed 1,000 m in thickness, as is the case for the Deming sub-basin in southern New Mexico. In alluvial basins, alluvium creates a thick geologic host medium that serves as the natural barrier system (NBS) in the conceptual model framework of a geologic disposal system and isolates the waste packages from receptors in the biosphere. Further, alluvial basin fill is typically comprised of stacked playa and lacustrine deposits along the basin axis (Perry et al., 2018). Characterized by low permeability, these layers protect the biosphere above the repository and the groundwater resources below the repository by limiting flow and transport of radionuclides.

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

Characterization of a Uranium Groundwater Plume along the Columbia River

The 300 Area is adjacent to the Columbia River on the southeastern end of Hanford, approximately one mile north of the city of Richland. Uranium fuel fabrication took place in the 300 Area from 1943 to 1988. In this chapter, the historical waste releases that occurred within the 300 Area of the Hanford Site, and specifically the resulting uranium plume, will be used as a case study. There are 3 different operable units (OU) that were established in the 300 Area for cleanup purposes: The 300-FF-1 OU, 300-FF-2 OU and the 300-FF-5 groundwater OU. To determine the location of contaminant plumes within the Hanford Site’s 300 Area, several puzzle pieces must be put together using multiple resources. The first piece of the puzzle is determining the location and chemistry of waste disposed of in the area. Records beginning in the 1940’s can help piece this information together, and the known waste sites and waste disposed of are an example of that process knowledge. Characterization of the subsurface and groundwater is important to determine where the contaminants are now, and what form they are in. This information can be used to design a remediation strategy, as well as to inform modeling efforts to predict where the contaminants will move to next.

300-FF-5, uranium, 300 area, hanford book project↗
NASA NTRS2018

Emerging Trends in Global Freshwater Availability

Freshwater availability is changing worldwide. Here we quantify 34 trends in terrestrial water storage (TWS) observed by the Gravity Recovery and Climate Experiment (GRACE) satellites during 2002-2016 and categorize their drivers as natural interannual variability, unsustainable groundwater consumption, or climate change. Several of these trends had been lacking thorough investigation and attribution, including massive changes in northwestern China and the Okavango delta. Others are consistent with climate model predictions. This observation-based assessment of how the world's water landscape is responding to human impacts and climate variations provides a blueprint for evaluating and predicting emerging threats to water and food security.

Okavango delta↗
NASA NTRS2023

Global Assessment of Drought Characteristics in the Anthropocene

Contemporary understanding of the impacts of climate change on global drought characteristics (e.g., intensities, extents) is still limited and not well understood. This knowledge is critical because projected changes in climate are expected to impact on future water availability as well as influence decisions on how water resources are allocated. The main aim of this study is to improve understanding of drought characteristics (extents and duration) in the Anthropocene where rapid changes in the environment are caused by the composite influence of human activities and climate change. Multi-scale earth observation data (1980 - 2020) and the Coupled Model Intercomparison Project Phase 6 climate models, which incorporate the Shared Socioeconomic Pathways (2040 - 2070 and 2070 - 2100) are used to assess these characteristics as well as identify climatic hotspots where changes in drought characteristics could drive groundwater hydrology. Results show that towards the end of the 21st century, global land areas under drought will significantly decrease but their durations will not. Generally, there is evidence of significant decline in the proportion of areas that will experience various drought intensities (moderate, severe and extreme drought) in the future and for each category, drought affected areas will not reach 30% on average. Moreover, some regions are potential hotspots of climate–groundwater interactions where drought events could directly impact on groundwater. This is because of the varying degree of strong correlations (positive and negative) between climate and groundwater data in some areas (e.g., Australia, Europe, Southern Africa, Asia). The relatively strong negative correlations in some of these hotspots are indicative of the presence of considerable lags, that could be caused by aridity as well as human groundwater footprints.

Christopher E Ndehedehe↗
NASA NTRS2015

Development of a Complete Landsat Evapotranspiration and Energy Balance Archive to Support Agricultural Consumptive Water Use Reporting and Prediction in the Central Valley, CA

Mapping evapotranspiration (ET) from agricultural areas in Californias Central Valley is critical for understanding historical consumptive use of surface and groundwater. In addition, long histories of ET maps provide valuable training information for predictive studies of surface and groundwater demands. During times of drought, groundwater is commonly pumped to supplement reduced surface water supplies in the Central Valley. Due to the lack of extensive groundwater pumping records, mapping consumptive use using satellite imagery is an efficient and robust way for estimating agricultural consumptive use and assessing drought impacts. To this end, we have developed and implemented an algorithm for automated calibration of the METRIC remotely sensed surface energy balance model on NASAs Earth Exchange (NEX) to estimate ET at the field scale. Using automated calibration techniques on the NEX has allowed for the creation of spatially explicit historical ET estimates for the Landsat archive dating from 1984 to the near present. Further, our use of spatial NLDAS and CIMIS weather data, and spatial soil water balance simulations within the NEX METRIC workflow, has helped overcome challenges of time integration between satellite image dates. This historical and near present time archive of agricultural water consumption for the Central Valley will be an extremely useful dataset for water use and drought impact reporting, and predictive analyses of groundwater demands.

valleys↗
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↗
DOE OSTI2022

Long-term missing value imputation for time series data using deep neural networks

We present an approach that uses a deep learning model, in particular, a MultiLayer Perceptron, for estimating the missing values of a variable in multivariate time series data. We focus on filling a long continuous gap (e.g., multiple months of missing daily observations) rather than on individual randomly missing observations. Our proposed gap filling algorithm uses an automated method for determining the optimal MLP model architecture, thus allowing for optimal prediction performance for the given time series. We tested our approach by filling gaps of various lengths (three months to three years) in three environmental datasets with different time series characteristics, namely daily groundwater levels, daily soil moisture, and hourly Net Ecosystem Exchange. We compared the accuracy of the gap-filled values obtained with our approach to the widely used R-based time series gap filling methods ImputeTS and mtsdi. The results indicate that using an MLP for filling a large gap leads to better results, especially when the data behave nonlinearly. Thus, our approach enables the use of datasets that have a large gap in one variable, which is common in many long-term environmental monitoring observations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Modeling Green Roof Potential to Mitigate Urban Flooding in a Chinese City

The Middle and Lower Reaches of the Yangtze River (MLRYR) region, which has humid subtropical climate conditions and unique plum rain season, is characterized by a simultaneous high-frequency urban flooding and reduction in groundwater levels. Retrofitting the existing buildings into green roofs is a promising approach to combat urban flooding, especially for a densely developed city. Here, the application potential of the Green Roof System (GRS) and the Improved Green Roof System (IGRS) designed to divert overflowing water from green roofs to recharge groundwater were analyzed in a densely developed city, Nanchang, China. For the first time, the influence of GRS on the hydraulic condition of Combined Sewage System/Storm Water System (CSS/SWS) is analyzed, which is a direct reflection of the effect of GRS on alleviating urban flooding. The simulation results show that GRS can retain about 41–75% of precipitation in a 2-hour timescale and the flooding volumes in the GRS/IGRS region are 82% and 28% less than those of the Traditional Roof System (TRS) in 10- and 100-yr precipitation events, respectively. In the continuous simulations, GRS also enhances Evapotranspiration (ET), which accounts for 39% of annual precipitation, so that reduces the cumulative surface runoff. Considering the IGRS can provide more hydrological benefits than the GRS under the same climate conditions, we may conclude that the widespread implementation of both the GRS and the IGRS in Nanchang and other densely developed cities in the MLRYR region could significantly reduce surface and peak runoff rates.

54 ENVIRONMENTAL SCIENCES↗