Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reservoir operation”

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.

810 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Simulation insights into wetting properties of hydrogen-brine-clay for hydrogen geo-storage

Hydrogen geo-storage is attracting substantial interdisciplinary interest as a cost-effective and sustainable option for medium- and long-term storage. Hydrogen can be stored underground in diverse formations, including aquifers, salt caverns, and depleted oil and gas reservoirs. The wetting dynamics of the hydrogen-brine-rock system are critical for assessing both structural and residual storage capacities, and ensuring containment safety. Through molecular dynamics simulations, we explore how varying concentrations of cushion gases (CO 2 or CH 4 ) influence the wetting properties of hydrogen-brine-clay systems under geological conditions (15 MPa and 333 K). We employed models of talc and the hydroxylated basal face of kaolinite (kaoOH) as clay substrates. Our findings reveal that the effect of cushion gases on hydrogen-brine-clay wettability is strongly dependent on the clay-brine interactions. Notably, CO 2 and CH 4 reduce the water wettability of talc in hydrogen-brine-talc systems, while exerting no influence on the wettability of hydrogen-brine-kaoOH systems. Detailed analysis of free energy of cavity formation near clay surfaces, clay-brine interfacial tensions, and the Willard-Chandler surface for gas-brine interfaces elucidate the molecular mechanisms underlying wettability changes. Our simulations identify empirical correlations between wetting properties and the average free energy required to perturb a flat interface when clay-brine interactions are less dominant. Here, our thorough thermodynamic analysis of rock-fluid and fluid-fluid interactions, aligning with key experimental observations, underscores the utility of simulated interfacial properties in refining contact angle measurements and predicting experimentally relevant properties. These insights significantly enhance the assessment of gas geo-storage potential. Prospectively, the approaches and findings obtained from this study could form a basis for more advanced multiscale simulations that consider a range of geological and operational variables, potentially guiding the development and improvement of geo-storage systems in general, with a particular focus on hydrogen storage.

25 ENERGY STORAGE↗
NASA NTRS2020

Encealdus Orbilander: A Flagship Mission Concept for Astrobiology

"Whether life exists beyond Earth remains a fundamental question driving our exploration of the Solar System. At Saturn’s moon Enceladus, plumes of oceanic material vented into space allow the investigation of the astrobiological potential of an ocean world, hinted at by Cassini, without the necessity of drilling through kilometers of ice crust. The Enceladus Orbilander is a flagship ($2.56B in fiscal year 2025 dollars) mission concept created for the 2023–2033 Planetary Science Decadal Survey. Orbilander takes full advantage of the opportunity provided by Enceladus’ plumes to search for signs of life. A single spacecraft both orbits and lands, capturing samples from four distinct reservoirs offered by the plumes. These samples, both particulate and vapor, are then analyzed by the Life Detection Suite (LDS), a set of five instruments conducting complementary and orthogonal biosignature-seeking measurements. To provide the context that specifically enhances interpretation of LDS measurements, geochemical and geophysical investigations are conducted both in orbit and on the surface. These reveal the physio-chemical state of the ocean and core as well as the processes involved in ejection of plume material and how these affect the ocean material analyzed by the LDS. The Orbilander can be delivered to the Saturn system via several launch vehicle and trajectory options, including a direct trajectory (7-year cruise), a ∆V-EGA trajectory (9-year cruise) and several options using an inner cruise with Venus and Earth flybys (10-year cruise). Upon Saturn Orbit Insertion, a 4-year moon tour pumps down the Orbilander’s orbit to intercept Enceladus. The most optimal arrival times balance the Jupiter flyby opportunities of the late 2030s and solar illumination at the Enceladus high southern latitudes where plume material is most abundant. This mission concept therefore targets project start in 2030. Upon Enceladus Orbit Insertion, the Orbilander begins a 1.5-year-long campaign of landing site reconnaissance, remote sensing science, and collecting sufficient plume sample to run all but one of the LDS measurements. After successful landing, the Orbilander spends 2 years on the surface conducting multiple LDS measurements with all five instruments on actively and passively collected plume material, as well as seismic investigations. The schedule laid out here is well-defined, but the mission also has operational and resource flexibility should additional reconnaissance be needed. As part of the design study, mission and development risks were identified and mitigation strategies proposed. Technologies key to achieving the life detection science objectives include instrumentation matured under programs like COLDTech and ICEE-2, such as aspects of the sampling system and microfluidic devices, as well as well-known techniques like high-resolution and separation-capable mass spectrometers. Autonomous onboard navigation is planned to maintain a halo orbit around Enceladus to enable passive sampling from orbit as well as reconnaissance measurements for use in site selection and landing. Terrain relative navigation is included to ensure safe landing, given that targeted areas may contain landing hazards. Continued development of radioisotope thermoelectric generator (RTG) technology and long-life batteries is essential for this long duration mission. The Enceladus Orbilander represents an optimal point in the trade space of science value versus cost, taking advantage of the extensive knowledge of Enceladus provided by Cassini, how well Enceladus lends itself to a search for life in material from its ocean, and the flexibility afforded by the innovative design developed by the APL team. By taking full advantage of Enceladus’ plumes both in orbit and on the surface, Orbilander represents a robust search for life with complementary and orthogonal biosignatures as well as contextual geophysical and geochemical measurements, determining not only whether Enceladus is inhabited (at levels up to 500,000× scarcer than in Earth’s oceans) but also why. "

Exobiology↗
NASA NTRS2012

The Discharging of Roving Objects in the Lunar Polar Regions

In 2007, the National Academy of Sciences identified the lunar polar regions as special environments: very cold locations where resources can be trapped and accumulated. These accumulated resources not only provide a natural reservoir for human explorers, but their very presence may provide a history of lunar impact events and possibly an indication of ongoing surface reactive chemistry. The recent LCROSS impacts confirm that polar crater floors are rich in material including approx 5%wt of water. An integral part of the special lunar polar environment is the solar wind plasma. Solar wind protons and electrons propagate outward from the Sun, and at the Moon's position have a nominal density of 5 el/cubic cm, flow speed of 400 km/sec, and temperature of 10 eV (approx. equal 116000K). At the sub-solar point, the flow of this plasma is effectively vertically incident at the surface. However, at the poles and along the lunar terminator region, the flow is effectively horizontal over the surface. As recently described, in these regions, local topography has a significant effect on the solar wind flow. Specifically, as the solar wind passes over topographic features like polar mountains and craters, the plasma flow is obstructed and creates a distinct plasma void in the downstream region behind the obstacle. An ion sonic wake structure forms behind the obstacle, not unlike that which forms behind a space shuttle. In the downstream region where flow is obstructed, the faster moving solar wind electrons move into the void region ahead of the more massive ions, thereby creating an ambipolar electric field pointing into the void region. This electric field then deflects ion trajectories into the void region by acting as a vertical inward force that draws ions to the surface. This solar wind 'orographic' effect is somewhat analogous to that occurring with terrestrial mountains. However, in the solar wind, the ambipolar E-field operating in the collision less plasma replaces the gradient in pressure that would act in a collisional neutral gas. Human systems (roving astronauts or robotic systems created by humans) may be required to gain access to the crater floor to collect resources such as water and other cold-trapped material. However, these human systems are also exposed to the above-described harsh thermal and electrical environments in the region. Thus, the objective of this work is to determine the nature of charging and discharging for a roving object in the cold, plasma-starved lunar polar regions. To accomplish this objective, we first define the electrical charging environment within polar craters. We then describe the subsequent charging of a moving object near and within such craters. We apply a model of an astronaut moving in periodic steps/cadence over a surface regolith. In fact the astronaut can be considered an analog for any kind of moving human system. An astronaut stepping over the surface accumulates charge via contact electrification (tribocharging) v.lith the lunar regolith. We present a model of this tribo-charge build-up. Given the environmental plasma in the region, we determine herein the dissipation time for the astronaut to bleed off its excess charge into the surrounding plasma.

Jackson, T. L.↗
NASA NTRS2007

Identification and Control of Aircrafts using Multiple Models and Adaptive Critics

We compared two possible implementations of local linear models for control: one approach is based on a self-organizing map (SOM) to cluster the dynamics followed by a set of linear models operating at each cluster. Therefore the gating function is hard (a single local model will represent the regional dynamics). This simplifies the controller design since there is a one to one mapping between controllers and local models. The second approach uses a soft gate using a probabilistic framework based on a Gaussian Mixture Model (also called a dynamic mixture of experts). In this approach several models may be active at a given time, we can expect a smaller number of models, but the controller design is more involved, with potentially better noise rejection characteristics. Our experiments showed that the SOM provides overall best performance in high SNRs, but the performance degrades faster than with the GMM for the same noise conditions. The SOM approach required about an order of magnitude more models than the GMM, so in terms of implementation cost, the GMM is preferable. The design of the SOM is straight forward, while the design of the GMM controllers, although still reasonable, is more involved and needs more care in the selection of the parameters. Either one of these locally linear approaches outperform global nonlinear controllers based on neural networks, such as the time delay neural network (TDNN). Therefore, in essence the local model approach warrants practical implementations. In order to call the attention of the control community for this design methodology we extended successfully the multiple model approach to PID controllers (still today the most widely used control scheme in the industry), and wrote a paper on this subject. The echo state network (ESN) is a recurrent neural network with the special characteristics that only the output parameters are trained. The recurrent connections are preset according to the problem domain and are fixed. In a nutshell, the states of the reservoir of recurrent processing elements implement a projection space, where the desired response is optimally projected. This architecture trades training efficiency by a large increase in the dimension of the recurrent layer. However, the power of the recurrent neural networks can be brought to bear on practical difficult problems. Our goal was to implement an adaptive critic architecture implementing Bellman s approach to optimal control. However, we could only characterize the ESN performance as a critic in value function evaluation, which is just one of the pieces of the overall adaptive critic controller. The results were very convincing, and the simplicity of the implementation was unparalleled.

Principe, Jose C.↗
DOE OSTI2021

Deep learning inversion of gravity data for detection of CO 2 plumes in overlying aquifers

In this work, we developed an effective U-Net based deep learning (DL) model for inversion of surface gravity data on a rectangular grid to predict 2-D high-resolution subsurface CO 2 distribution along a vertical cross-section due to CO 2 leakage through a wellbore within a deep CO 2 storage reservoir. We used synthetic data to model two types of CO 2 leakage scenarios: one CO 2 plume in a shallow aquifer (single plume case), and two plumes present at different depths (double plume case). The 3-D synthetic plume samples were created by sampling among predetermined CO 2 plume depths, saturations, and volumes. The corresponding surface gravity data on a rectangular grid were generated by a 3-D forward model. The U-Net model detected 72% of single-plume samples, and one or both plumes in 75% of double-plume samples. Most of the undetected single plumes have small gravity field strengths below the typical noise level of 5 μGal. This model generated reproducible, reliable predictions with acceptable errors and demonstrated improved spatial resolution over the conventional least-squares inversion. In contrast to the conventional least-squares inversion, which often overestimates the size of its target and underestimates its density, this U-Net model accurately delineated the boundary of a target. Furthermore, this DL inversion detected deep, small, or low saturation CO 2 plumes that are often more difficult to resolve with conventional gravity inversion methods. We note the limitations of this feasibility study, including the use of synthetic data with regular CO 2 plume shapes, and the prediction of a 2-D plume cross-section rather than the full 3-D plume, as well, we recognize the lower detection fraction for double-plume scenarios. Nevertheless, this study demonstrates that DL gravity inversion is a promising and potentially superior method to conventional least-squares inversion. Our U-Net based deep learning inversion approach may be adapted for inversion of other types of geophysical data. DL inversion can facilitate near real-time monitoring of geologic carbon sequestration to provide site operators with prompt information about subsurface CO 2 distribution for risk management and mitigation.

58 GEOSCIENCES↗
DOE OSTI2021

Mechanistic basis of post-treatment control of SIV after anti-α4β7 antibody therapy

Treating macaques with an anti-α4β7 antibody under the umbrella of combination antiretroviral therapy (cART) during early SIV infection can lead to viral remission, with viral loads maintained at < 50 SIV RNA copies/ml after removal of all treatment in a subset of animals. Depletion of CD8+ lymphocytes in controllers resulted in transient recrudescence of viremia, suggesting that the combination of cART and anti-α4β7 antibody treatment led to a state where ongoing immune responses kept the virus undetectable in the absence of treatment. A previous mathematical model of HIV infection and cART incorporates immune effector cell responses and exhibits the property of two different viral load set-points. While the lower set-point could correspond to the attainment of long-term viral remission, attaining the higher set-point may be the result of viral rebound. Here we expand that model to include possible mechanisms of action of an anti-α4β7 antibody operating in these treated animals. We show that the model can fit the longitudinal viral load data from both IgG control and anti-α4β7 antibody treated macaques, suggesting explanations for the viral control associated with cART and an anti-α4β7 antibody treatment. This effective perturbation to the virus-host interaction can also explain observations in other nonhuman primate experiments in which cART and immunotherapy have led to post-treatment control or resetting of the viral load set-point. Interestingly, because the viral kinetics in the various treated animals differed—some animals exhibited large fluctuations in viral load after cART cessation—the model suggests that anti-α4β7 treatment could act by different primary mechanisms in different animals and still lead to post-treatment viral control. This outcome is nonetheless in accordance with a model with two stable viral load set-points, in which therapy can perturb the system from one set-point to a lower one through different biological mechanisms.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2026

Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows

A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.

stagnation heat flux↗
NASA NTRS2010

Lidar Measurements of Methane and Applications for Aircraft and Spacecraft

Atmospheric methane levels have remained relatively constant over the last decade around 1.78 parts per million (ppm) but observations since 2007 show that levels may be increasing. This trend may be caused by increased fossil fuel production, rice farming, livestock and landfills, but the underlying causes are quite uncertain. One hypothesis is that reservoirs of carbon trapped in the permafrost regions of northern Canada, Europe, and Siberia thaw as global temperatures rise and are releasing increasing amounts of methane. Another hypothesis points to increased production of methane by microbes as the permafrost warms. Currently most observations of greenhouse gases are limited to in-situ (surface and tower sites) and limited airborne in-situ measurements. Space column density measurements are starting to become available from the GOSAT mission. Although methane survives for a shorter time in the atmosphere than CO2, its impact on climate change per molecule is about 23 times than that of CO2. Accurate global observations of several greenhouse gases, including methane, are urgently needed in order to better understand climate change processes and to reduce the uncertainty in the carbon budget. Differential absorption lidar is a well-established technique to measure atmospheric gases, and methane has optical absorption bands near 1.65,2.2,3.4 and 7.8 micron. The near infrared overtones lines of CH4 near 1650 nm are relatively free of interference from other species. There are absorption lines near 1651 nm which are both temperature insensitive and have line strengths well suited for lidar measurements. We have developed a laser and demonstrated lidar measurements of CH4 using lines in this band. Our laser uses a narrow linewidth 1064 nm laser pulse passing through a nonlinear crystal. We generate the tunable laser signals near 1651 nm by using the optical parametric amplification (OPA) process. Inside the crystal the 1064 nm beam overlaps with an injection seed laser near 1651 nm from a wavelength tunable diode laser. Incident photons from the pump laser pulse are converted into two photons, with one at the wavelength of the injection seeder. The wavelength of the OPA output is tuned via the wavelength of diode laser. Our laser is tunable, operates near 1651 nm and generates approximately 4 uJ/pulse at 6 KHz. We vary the emission wavelengths within this band by tuning the diode laser's wavelength. We have used this OPA transmitter to make measurements of CH4 at various pressures in a gas cell and over open outdoor horizontal paths. We have measured the lineshape of methane in a 6 cm long cell at various energy levels with this transmitter, with excellent agreement with the lineshape calculated by HITRAN. We have also measured the absorption lineshape of atmospheric methane in an open 3 km outdoor path. The agreement between the measurements and HITRAN, for 1746 ppb and 760 Torr was quite good. We have also made pulsed two wavelength lidar measurements of methane line absorption in the column to a tower at 1.5 km range. These used on- and off-line wavelengths of 1650.957 nm, and 1651.072 nm, and a 20 cm diameter receiver telescope with an infrared PMT detector. The absorption of the on-line photons was 30%. The methane column absorption was estimated via HITRAN, and was in good agreement with the expected methane absorption for a concentration of 1750 ppm. Finally we have calculated the measurement performance of an airborne methane lidar using this transmitter, as well as the energy and telescope scaling needed for a lidar for space. These results, and more details of our experiments will be described in the presentation.

Riris, Haris↗
NASA NTRS2022

An Empirical Algorithm to Map Perennial Firn Aquifers and Ice Slabs Within the Greenland Ice Sheet Using Satellite L-band Microwave Radiometry

Perennial firn aquifers are subsurface meltwater reservoirs consisting of a meters-thick water-saturated firn layer that can form on spatial scales as large as tens of kilometers. They have been observed within the percolation facies of glaciated regions experiencing intense seasonal surface melting and high snow accumulation. Widespread perennial firn aquifers have been identified within the Greenland Ice Sheet (GrIS) via field expeditions, airborne ice-penetrating radar surveys, and satellite microwave sensors. In contrast, ice slabs are nearly continuous ice layers that can also form on spatial scales as large as tens of kilometers as a result of surface and subsurface water-saturated snow and firn layers sequentially refreezing following multiple melting seasons. They have been observed within the percolation facies of glaciated regions experiencing intense seasonal surface melting but in areas where snow accumulation is at least 25 % lower as compared to perennial firn aquifer areas. Widespread ice slabs have recently been identified within the GrIS via field expeditions and airborne ice-penetrating radar surveys, specifically in areas where perennial firn aquifers typically do not form. However, ice slabs have yet to be identified from space. Together, these two ice sheet features represent distinct, but related, sub-facies within the broader percolation facies of the GrIS that can be defined primarily by differences in snow accumulation, which influences the englacial hydrology and thermal characteristics of firn layers at depth. Here, for the first time, we use enhanced-resolution vertically polarized L-band brightness temperature () imagery (2015–2019) generated using observations collected over the GrIS by NASA's Soil Moisture Active Passive (SMAP) satellite to map perennial firn aquifer and ice slab areas together as a continuous englacial hydrological system. We use an empirical algorithm previously developed to map the extent of Greenland's perennial firn aquifers via fitting exponentially decreasing temporal L-band signatures to a set of sigmoidal curves. This algorithm is recalibrated to also map the extent of ice slab areas using airborne ice-penetrating radar surveys collected by NASA's Operation IceBridge (OIB) campaigns (2010–2017). Our SMAP-derived maps show that between 2015 and 2019, perennial firn aquifer areas extended over 64 000 km2, and ice slab areas extended over 76 000 km2. Combined together, these sub-facies are the equivalent of 24 % of the percolation facies of the GrIS. As Greenland's climate continues to warm, seasonal surface melting will increase in extent, intensity, and duration. Quantifying the possible rapid expansion of these sub-facies using satellite L-band microwave radiometry has significant implications for understanding ice-sheet-wide variability in englacial hydrology that may drive meltwater-induced hydrofracturing and accelerated ice flow as well as high-elevation meltwater runoff that can impact the mass balance and stability of the GrIS.

Julie Z Miller↗
DOE OSTI2025

Environmental DNA as a tool for hydropower impact assessments: current status, special considerations, and future integration

Globally there is an urgent need to find sustainable solutions to balance energy production with the protection of vulnerable species and conservation of biodiversity. This is particularly critical for freshwater ecosystems, habitats, and species that may be impacted by hydropower development and operations needed to meet energy grid demands. Reliable and accurate environmental impact assessments (EIAs) that identify the biological, physical, or social impacts of hydropower are key to ensure biodiversity, ecosystem, and societal sustainability. The analysis of environmental DNA (eDNA) has the potential to transform hydropower EIAs, management and mitigation planning, and decision-making procedures. Further, the incorporation of eDNA surveys into EIAs during both hydropower planning and continued operations may streamline regulatory processes by improving our understanding of potentially impacted biota and habitats and evaluating environmental impacts mitigation. Here, we: (i) highlight current understanding and use of eDNA in freshwater environments; (ii) examine critical considerations for eDNA integration into hydropower EIAs and biological monitoring; (iii) identify knowledge gaps in eDNA analysis and applications unique to hydropower-regulated systems; and (iv) discuss future opportunities to bolster the incorporation of eDNA into hydropower research including regulatory acceptance and public engagement. While we acknowledge that there are several factors that may complicate the broad adoption of eDNA as a tool for assessing the impacts of hydropower, we anticipate that growing confidence in eDNA through hydropower-specific protocols, calibrations, and validations will overcome these inherent uncertainties.

aquatic biodiversity↗
DOE OSTI2021

DOE Award No.: DE-FE0023919 Phase 4 Scientific/Technical Report

This is the Phase 4 Report for the ‘Deepwater Methane Hydrate Characterization and Scientific Assessment or Genesis of Methane Hydrates in the Gulf of Mexico (GOM2)’ research project (DOE Award No. DE-FE0023919). The report summarizes activities from October 1, 2019 to September 30, 2020. The project is led by the University of Texas at Austin (UT). The project objective is to gain insight into the nature, formation, occurrence and physical properties of methane hydrate-bearing sediments for the purpose of methane hydrate resource appraisal through the planning and execution of drilling, coring, logging, testing and analytical activities that assess the geologic occurrence, regional context, and characteristics of marine methane hydrate deposits in the Gulf of Mexico outer continental shelf (OCS). We published a dedicated American Association of Petroleum Geologists Bulletin volume describing initial results from the UT-GOM2-1 expedition in Sept. 2020. This is part 1 of a multi-volume commitment by AAPG to this project. We further confirmed that the natural gas in hydrate at GC-955 was formed by primary microbial processes (>76.1 %). The in-situ effective permeability hydrate-bearing sandy silts at the GC-955 reservoir ranges from 0.1 md (1.0×10-16 m2) to 2.4 md (2.4×10-15m2) in cores with 83% to 93% hydrate saturation. The intrinsic permeability (the single phase permeability) is estimated from reconstituted samples to be ~12 md (1.2×10-14 m2) to ~41 md (4.1×10-14 m2). We used observation and models to interpret that the core degradation that is found in pressure cores is due to dissociation of the methane hydrate in the outer circumference of the core and dissolution of that methane into the fresh pore water that the core is stored with. We are designing approaches to minimize this core loss in the future. We spent an enormous amount of effort to further improve the ability of the pressure coring tool (the PCTB) to pressure seal correctly. We completed upgrading the upper section of the PCTB to address poor pressure. We successful tested the modifications at Geotek’s test facility in Salt Lake City (Bench Test II). We completed a Land Test of the PCTB at the Schlumberger Cameron Test and Training Facility (CTTF). The tool did not seal in 6 out of 7 tests and we clearly demonstrated that cuttings were wedging in the ball valve assembly, keeping the ball valve from sealing. We reproduced the failure mechanism observed during the land test at Salt Lake City and confirmed the sensitivity of the ball valve assembly to grit. Geotek designed and tested 9 modifications to address this issue and the PCTB is now 100% successfully sealing in the presence of grit. Our science expedition is scheduled for spring 2022 and we are fully focused on preparing for this. UT and Ohio State completed a Shallow Hazard Assessment report for each proposed UT-GOM2-2 drilling location, pursuant to 30 CFR 250.214(f) and 250.244 (f). The Shallow Hazard Reports will accompany the UT-GOM2-2 Exploration Plan that is submitted to BOEM, and completes the geological and geophysical analysis for UT-GOM2-2 permitting efforts. We updated the UT-GOM2-2 Operations Plan (Version 1). We completed the UT-GOM2-2 Science and Sample Distribution Plan (Version 1). We evaluated the scope, budget, and schedule that would result from using a commercial vessel. We developed detailed drilling schedule, mud volume, and resource estimates. We developed a vessel specification document and a well plan, and sent these documents to prospective vessel contractors.

03 NATURAL GAS↗
NASA NTRS2017

Constraints on Mechanisms of Chondrule Formation from Chondrule Precursors and Chronology of Transient Heating Events in the Protoplanetary Disk

The mineralogy, petrography, and oxygen-isotope compositions of porphyritic chondrules—dominant chondrule type in most chondrite groups—suggest formation by incomplete melting of isotopically diverse precursors during localized transient heating events in dust-rich regions of the protoplanetary disk characterized by 16O-poor compositions (D17Odust+gas~ –7‰ to +4‰) relative to the inferred Sun’s value (D17O ~ –28 ± 2‰). The chondrule precursors included Ca,Al-rich inclusions (CAIs), amoeboid olivine aggregates (AOAs), chondrules of earlier generations, fine-grained matrix-like material, and possibly fragments of pre-existing planetesimals. Like porphyritic chondrules, igneous CAIs formed by melting of isotopically diverse precursors during transient heating events, but in an isotopically distinct, solar-like reservoir of the protoplanetary disk (D17O dust + gas ~ –24‰), probably near the proto Sun. Based on a narrow range of the initial 26Al/27Alratios inferred from the internal Al-Mg isochrons in igneous CAIs, their melting started at the very beginning of Solar System formation (t0), defined by the CV CAIs with U-corrected Pb-Pb age of 4567.3 ± 0.16 Ma and the canonical 26Al/27Al ratio of (5.25 ± 0.02) ¥ 10–5, and lasted at least 0.3 Ma. The U-corrected Pb-Pb absolute and 26Al-26Mg relative ages of porphyritic chondrules from type 3 ordinary, CO, CV, and CR carbonaceous chondrites (assuming uniform distribution of 26Al in the disk at the canonical level) suggest chondrule formation started at t0 and lasted for about 4 Ma. These observations may preclude formation of the majority of porphyritic chondrules by splashing of differentiated planetesimals and by collisions between planetesimals; instead, they are consistent with melting of dust balls by bow shocks or magnetized turbulence in the disk. Some porphyritic chondrules in equilibrated (petrologic type 4–6) ordinary chondrites contain relict fragments of coarse-grained chromite, ilmenite, phosphates, and albitic plagioclase. The similar mineral assemblage is commonly observed in type 4–6 ordinary chondrites, but is absent in type 3 chondrites, suggesting these chondrules formed by incomplete melting of thermally metamorphosed ordinary chondrite material, possibly by impacts. The CB metal-rich carbonaceous chondrites contain exclusively magnesian non-porphyritic chondrules crystallized from complete melts. These chondrules formed in a gas-melt plume generated by a hypervelocity (≥20 km/s) collision between planetesimals ~4.8 Ma after t0 in a transition or a debris disk. One of the colliding bodies was probably differentiated. The CH metal-rich carbonaceous chondrites contain chondrules formed by different mechanisms. The magnesian non-porphyritic chondrules formed in the CB impact plume ~4.8 Ma after t0. The chemically diverse (magnesian, ferroan, and Al-rich) porphyritic chondrules formed by incomplete melting of isotopically diverse precursors in the protoplanetary disk, most likely prior the CB impact plume event. We conclude that there are multiple mechanisms of chondrule formation that operated over the entire life-time of the disk.

Chondrules↗
DOE OSTI2023

Hydroelectric Power and Hydrogen Production Integration

Hydropower-based hydrogen production could introduce opportunities for new revenue streams for hydropower plants, including from energy storage and regeneration as well as from sale of the hydrogen product to external markets. Hydrogen-based energy storage and regeneration could also help support Idaho Power’s decarbonization goals by decreasing dependence on fossil-based peaking power plants. Additionally, integration of hydrogen production with hydropower generation could help address the issue of low dissolved oxygen river water conditions that commonly accompany hydropower plant operations by utilizing the oxygen byproduct from an electrolytic hydrogen production process as a resource for mitigation of low dissolved oxygen water conditions. Comprehensive techno-economic analysis of the hybrid hydroelectric and hydrogen energy storage system has revealed critical insights into the pathways and considerations for optimizing the economic value and environmental benefits of such systems. Among the three identified pathways of natural gas blending, regeneration, and direct sale of hydrogen, the direct sale of hydrogen emerges as the most profitable, particularly given the current pricing dynamics of hydrogen and electricity. However, as we anticipate a future grid characterized by higher renewable energy penetration, the landscape may evolve, featuring lower average electricity prices, heightened fluctuations, and more significant seasonal variations. Consequently, the attractiveness of electricity regeneration through hydrogen and the benefits of long-duration hydrogen storage are expected to increase substantially in such a dynamic energy scenario. The careful selection of component sizes within the hybrid system proves to be paramount for ensuring cost-effectiveness. Notably, the size of the hydrogen market has emerged as a critical determinant for the optimal size of the electrolyzer. Hydrogen storage should be sized to meet energy shifting requirements. The appropriate size of the fuel cell/microturbine generator hinges on the shape of electricity prices and the available revenue streams derived from participating in grid services. Striking the right balance among these components is essential for maximizing the overall efficiency and profitability of the hydrogen facility. Furthermore, the by-product of electrolysis, namely oxygen, introduces an additional dimension to the system's functionality. The oxygen generated can be effectively utilized for dissolved oxygen (DO) mitigation, particularly with larger electrolyzer sizes capable of satisfying the complete oxygen demand for this purpose. While the economic benefits derived from saved oxygen purchase costs may be relatively modest compared to other revenue streams, the environmental advantages of repurposing oxygen for DO mitigation could help Idaho Power meet their environmental obligations. In addition to the identified factors shaping the viability of hybrid hydrogen production and hydroelectric generation, it is noteworthy that the integration of hydrogen energy storage offers a unique advantage during unusually wet years. In such periods of increased water inflow, the hydrogen storage capacity serves as a valuable supplement to the reservoir. By utilizing hydrogen energy storage as a complementary reservoir, the system gains flexibility in reservoir management when dealing with fluctuations in water availability.

08 HYDROGEN↗
NASA NTRS2004

ICESat Observations of Inland Surface Water Stage, Slope, and Extent: a New Method for Hydrologic Monitoring

River discharge and changes in lake, reservoir and wetland water storage are critical terms in the global surface water balance, yet they are poorly observed globally and the prospects for adequate observations from in-situ networks are poor (Alsdorf et al., 2003). The NASA-sponsored Surface Water Working Group has established a framework for advancing satellite observations of river discharge and water storage changes which focuses on obtaining measurements of water surface height (stage), slope, and extent. Satellite laser altimetry, which can achieve centimeter-level elevation precision for single, small laser footprints, provides a method to obtain these inland water parameters and contribute to global water balance monitoring. Since its launch in January, 2003 the Ice, Cloud, and land Elevation Satellite (ICESat), a NASA Earth Observing System mission, has achieved over 540 million laser pulse observations of ice sheet, ocean surface, land topography, and inland water elevations and cloud and aerosol height distributions. By recording the laser backscatter from 80 m diameter footprints spaced 175 m along track, ICESat acquires globally-distributed elevation profiles, using a 1064 nm laser altimeter channel, and cloud and aerosol profiles, using a 532 nm atmospheric lidar channel. The ICESat mission has demonstrated the following laser altimeter capabilities relevant to observations of inland water: (1) elevation measurements with a precision of 2 to 3 cm for flat surfaces, suitable for detecting river surface slopes along long river reaches or between multiple crossings of a meandering river channel, (2) from the laser backscatter waveform, detection of water surface elevations beneath vegetation canopies, suitable for measuring water stage in flooded forests, (3) single pulse absolute elevation accuracy of about 50 cm (1 sigma) for 1 degree sloped surfaces, with calibration work in progress indicating that a final accuracy of about 12 cm (1 sigma) will be achieved for clear atmosphere conditions, suitable for detection of stage changes through time, (4) ability to precisely point the spacecraft so as to position the laser profile on the Earth the surface with a cross-track accuracy of 50 m (1 sigma), enabling small water bodies and specific locations to be targeted and re-observed through time, (5) adequate signal levels from specular water surfaces up to 5 degrees off-nadir, enabling complete global access to any location on the Earth's surface from the ICESat repeat orbit by off-nadir pointing, and (6) day and night operation with successful laser ranging to the Earth's surface through thin to moderate cloud cover, enabling more frequent measurements than can be achieved by passive optical sensors. Here we illustrate these capabilities by showing ICESat observations through time for selected river and lake locations.

Harding, David J.↗
NASA NTRSDate not supplied

Development of Machine Learning Algorithms to Segment and Study Images of Astromaterial Samples

Introduction: Micrometer-scale chemical analyses of chondritic meteorites and mission-returned asteroid samples can reveal details of the physical and chemical processes operating in the early solar system, including processes that gave rise to planets, moons, and minor bodies. These primitive astromaterials are comprised of chondrules, calcium- and aluminum-rich inclusions (CAI), and many other silicates, oxides, metals, sulfides, and fine-grained materials. The chemical and mineralogical complexity of these samples, vast populations of different components, and heterogeneity across mm to km scales, all limit our understanding of the origin and evolution of these materials. Here, we describe recent efforts to use machine learning techniques to automate the segmentation of chemical maps of chondritic meteorites, designed to aid studies of asteroid samples returned by spacecraft. By automating the task of segmentation it will become possible to rapidly analyze and interpret the sizes, shapes, mineralogy, chemistry, and other properties of every chondrule, calcium- and aluminum-rich inclusion (CAI) and other clast within and between asteroid samples. Sample return missions significantly accelerate and heighten the need to develop such new data analysis techniques, and associated data repositories. Techniques: Neural networks require abundant training data, i.e. images which have been segmented by a human user. We have manually segmented data available from previous petrologic and chemical work at NASA Johnson Space Center and the American Museum of Natural History [1-4]. These data were derived from energy- and wavelength-dispersive X-ray spectroscopy (EDS, WDS) mapping of samples from many chondrite groups. The Deeplabv3+ [5] neural network architecture was trained on human-labeled masks and used to create machine-labeled masks. Several different algorithms were investigated, with inputs ranging from common RGB image formats through to hyperspectral datasets, with raw data comprising greyscale maps of Mg, Ca, and Al, with or without Si, Fe, Ti for both EDS and WDS data, and extending to other elements in EDS only. Each greyscale image was paired with a binary mask for each labelled particle type. Results: The trained algorithms can segment (Fig 1), classify, and measure the dimensions of thousands of particles in chemical maps of a standard 1-inch round petrographic section in seconds to minutes, rather than many hours needed by a human. Accuracy of the algorithms varied from chondrite to chondrite and across particle types. Further results and details of the algorithms will be presented at the workshop. Future directions: Machine learning has the potential to revolutionize our understanding of complex particle populations contained within primitive astromaterial, with segmentation being a critical first step. Example applications include better understanding of particle transport, nebular reservoirs, parent body accretion, and a deeper understanding of the relationships between particle populations and bulk rock elemental and isotopic compositions. In addition to benefits that machine learning can bring to individual researchers, building a community data repository of thousands to millions of particles across hundreds of samples will open up many other possibilities. For example, with a large enough dataset it will be possible to search for exceptionally closely matching particles across disparate samples. Such a capability would enable a single CAI from OSIRISREx or Hayabusa/II samples to be matched to chondritic CAIs that exhibit near-identical size, texture, and mineralogy, down to the level of similar core phenocrysts, zonation, and rim sequences. Such comparative analyses will help to disentangle precursor chemistry, chronology, gas/dust reservoirs during heating, and accretion. Such an endeavor would be impossible without machine learning and a large community data repository of astromaterial chemical/mineralogic maps.

Machine Learning↗
NASA NTRS1998

Deep Basalt Aquifers in Orcus Patera, Elysium Basin Mars: Perspectives for Exobiology Exploration

Direct indicators of shorelines, spillways, and terraces allowed to determine the extent of the Elysium Paleolake between the contour-lines 1000 and 500 m below the Martian datum. The Elysium Paleolake is bordered north by Orcus Patera (14N/181W), which lies west of the Tartarus Montes and Tartarus Colles. The Orcus Patera displays an ellipse-shaped collapsed caldera of 360-km long and 100-km wide. Viking topographic data show that the bottom of the caldera is located at 2500 below the Martian datum, and surrounded by a steep-walled ram art which crest is located at about 0 m elevation. Considering the localization of Orcus Patera in the Elysium paleolake, its altimetry, and the magmatic origin of this caldera, we propose the existence of a paleolake in Orcus Patera generated (a) by juvenile water from magma during the Noachian period, and (b) by intermittent influx of the Elysium Basin from Hesperian to Amazonian. Results are encouraging to consider this site as a potential high-energy source environment for microbial communities. are circumscribed by a 50-km wide lava field mapped as Noachian material. The structure of Orcus Patera represents the record of material erupted from a magmatic reservoir. The caldera is enclosed by steep inner walls (25% measured from topographic data), values which could be in agreement with the presence of a deep magmatic reservoir, as suggested by the typology of Crumpler et.al. The depth of the caldera might be due to the collapse of the magma reservoir, and the release of gases accompanying the magma thermal evolution. Origins of water for the paleolake(s): The water that generated a paleolake in Orcus Patera may have come from two origins: (1) Juvenile water: Plescia and Crips estimated a magma H20 content by weight between 0.5% and 1.5% using for the first value a comparison with terrestrial basalt, and for the second values from a Martian meteorite. The amount of H20 can be estimated by the volume of erupted lava, and the lava content of the caldera. In this study, we adopt a water content of 1%. The total volume of magma that has been contained in the caldera, and the volume of lava contained in the observed lava field is about 110 x 10(exp 6) cubic km, that gives a total volume of 1.10 x 10(exp 6) cubic km of water. The juvenile water expelled by the overpressure within the magma chamber charged with desolved water-vapor may have moved into the crust. The decrease in overburden pressure led to bubble formation. The ascent of these bubbles generated a pressurization of the magma, which was sufficient to fracture the overlaying magma layer, (2) Water from Elysium paleolake. During the Amazonian, the rise of the Elysium paleolake level generated an overspilling that supplied the caldera with water. The southern portion of the crest shows a deep gap 12-km wide at -1500 m elevation, locating the gap between 500 to 1000 in below the assumed water of Elysium paleolake, thus facilitating the influx of Elysium paleolake water into Orcus Patera. Bathymetric calculations give a floor area of 25,500 sq km at -2000 m elevation, and a water volume of 42,000 cubic km, with a lake-level at -1500 m. A substantial amount of water may have percolated through the fractured lava, and part of the volume may have overspilled the northern crest of Orcus Patera to debouch in the Tartarus Montes region. We envision the formation of a subsurface aqueous environment in basaltic rocks at the contact of the two water-source origins, possibly the percolating surface lake water, and more likely the juvenile water. Similarly to terrestrial calderas, Orcus Patera might be surrounded by ring-fractures caused by the collapse of the magma chamber that followed the release of gases. These ring-fractures may have been covered later by sedimentation in the caldera (lacustrine, aeolian, and volcanic), and by mass wasting. The detumescence of the magma in the caldera, and the vesiculation of the juvenile water may have operated simultaneously. Comparatively to terrestrial melts, Martian iron-rich melts are denser. This greater density implies greater effusion rates (eight-times terrestrial values), and larger fissuration widths (two-times terrestrial ones). With increasing vesiculation of magma, the bubbles interact with one-another because there are of similar pressure. They make a magma froth at the contact with the caldera surface, and on the walls of the fractures. In the saturated magma, froth, where the volume ratio of gases-to-liquid is about 4:1, the bubbles form a huge surface area of interconnected spaces. Bubbles near the caldera surface disrupt the magma, and fragmentation takes place, which moves downward through the magma column. On Earth, the bubbles are likely to grow between 1 and 50 mm in diameter due to the difference between the magma surface tension, and the bubble supersaturation pressure. The Martian low-pressure at surface level is likely to accelerate the expansion of the bubbles, and increase their final diameter and number, creating more voids in the magma. The strong magma froth with enclosed juvenile water bubbles interconnected with exsolved gas bubbles constitute a potential geothermal environment for geochemical energy production from basalt and water that does not require excessive temperatures. This process can start at +20C. Similar types of environments have been shown on Earth as potential energy sources for microbial metabolism, and could have provided deep aqueous basaltic niches for possible Martian microorganisms, even geologically recently. During the Amazonian, combination of volcanism and water activity still existed on Mars. Moreover, this type of potential niches open ways for investigation of possible oases of extinct or extant life, not only on paleolakes, and surface hydrothermalism spring areas, but also all large systems of fossae, which combine hydrologic and volcanic activities, and which provide an energy source, and an underground shelter to prevent surface UV bombardment. Additional information contained in the original.

Grin, E. A.↗
DOE OSTI2024

Idaho National Laboratory Integrated Multisite SSHAC Level 3: Probabilistic Volcanic Hazards Assessment

The Idaho National Laboratory (INL) resides on the eastern Snake River Plain (ESRP), part of the Snake River Plain (SRP) with a complex origin and geologic history of volcanism. Much of the Quaternary (last 2.58 million years) volcanism within 400 km of INL is genetically associated with a major thermal anomaly referred to as the Yellowstone hotspot, which is currently located more than 180 km northeast of INL. Near INL, local volcanic sources include silicic domes near its southern border, and numerous dike-fed basaltic vents of the ESRP, some of which are near or within the INL boundaries. Quantitative probabilistic assessments of screened-in volcanic hazardous phenomena for nine different facility complexes at the INL are presented for a Senior Seismic Hazard Analysis Committee (SSHAC) Level 3 (SL3) study. The comprehensive, integrated multisite SSHAC study consists of a single regional Probabilistic Volcanic Hazards Assessment (PVHA) that pertains to all nine INL facility complexes, with site-specific information developed at each respective area of interest (AOI), referred to as the "INL facility AOI" (Figure ES-1). Hazard products are generated for specified INL facility AOIs for use by multiple stakeholders from different agencies in risk-informed decision-making regarding site selection, operations, and design of nuclear facilities at INL, consistent with U.S. Department of Energy (DOE) and U.S. Nuclear Regulatory Commission (NRC) regulatory guidance. The study also serves as the basis for future periodic safety assessments required for existing DOE facilities, such as 10-year evaluations of natural-phenomena hazards. Elements of the INL PVHA including initial characterization, screening, quantitative assessments of eruption and hazard potential, and consideration of facility needs are conducted using the three phases of the SSHAC process: evaluation, integration, and documentation. As per regulatory guidance, the SSHAC framework provided the necessary processes and procedures for the PVHA Technical Integration (TI) team, at three workshops, five formal working meetings and many TI team remote meetings, to conduct initial characterization, screen the volcanic hazards, create PVHA model inputs, exercise those models in the PVHA, consider facility-specific volcanic hazard needs, and perform final hazard calculations. The Participatory Peer Review Panel (PPRP) provided independent oversight and performed process and technical reviews of the PVHA throughout its duration. The study included an extensive New Data Collection and Analyses (NDCA) program developed by the PVHA TI team to reduce uncertainties in hazard-significant elements in the PVHA model. NDCA activities generated 20 reports providing important contributory datasets and results to the project database for characterizing the ESRP. For example, a report compiling the dimensions of ESRP shield volcanoes and lava fields was used to construct volcanic footprints (areas of impact), was compared with data from INL subsurface cores, and was used to validate the results of lava-flow inundation modeling on the contemporary terrain. Another example is the acquisition of aeromagnetic data over INL and its surrounding area, with maps of buried magmatic features (e.g., subsurface volcanoes and swarms of feeder dikes) that informed the PVHA conceptual model of volcanism. The SSHAC evaluation process was used to conduct all elements of the PVHA including initial characterization and screening. Existing data and NDCA activities provided the foundation to develop the tectonomagmatic conceptual model of volcanism for the region of geographical interest in the SRP and Yellowstone hotspot volcanic system, and for considering volcanoes in the western US. Considering Quaternary volcanic sources active during this period, the screening approach identified and evaluated magma compositions, types of eruptive and intrusive phenomena, types of hazardous phenomena, and proximity of sources to INL. The PVHA TI team evaluated 18 types of magmatic sources in terms of 20 potentially hazardous phenomena, resulting in 360 screening decisions. The screening process resulted in 140 screened-in hazardous volcanic phenomena for Quaternary volcanic sources 1) proximal to INL facility complexes in the ESRP (64 basaltic and 54 silicic), 2) regional sources associated with the Yellowstone caldera system and Blackfoot Reservoir volcanic field (7), and 3) more distal sources from thirteen Cascade volcanoes and two volcanoes at Long Valley caldera (CA).

58 GEOSCIENCES↗
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↗