Search NASA⌕ Search

SEARCH · Search NASA

Results for “reservoir lining systems”

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.

68 records · Page 4

Detection of young (≤20 Myr) stellar populations in apparently quenched low-mass galaxies using red spectral line indices

ABSTRACT We report on the detection of a small contribution (around and below 1 per cent in mass) from young stellar components with ages ≤20 Myr in low-mass galaxies purposely selected from the MaNGA survey to be already-quenched systems. Among the sample of 28 galaxies, 8 of them show signatures of having suffered a very recent burst of star formation. The detection has been done through the analysis of line-strength indices in the red spectral range [5700,8800] Å. The increasing contribution of red supergiants to this red regime is responsible for a deviation of the index measurements with respect to their position within the model grids in the standard spectral range [3600,5700] Å. We demonstrate that a combination of red indices, as well as a qualitative assessment of the mean luminosity-weighted underlying stellar population, is required in order to distinguish between a true superyoung population and other possible causes of this deviation, such as abundance ratio variations. Our result implies that many presumably quenched low-mass galaxies actually contain gas that is triggering some level of star formation. They have, therefore, either accreted external gas, internally recycled enough gas from stellar evolution to trigger new star formation, or they kept a gas reservoir after the harassment or stripping process that quenched them in the first place. Internal processes are favoured since we find no particular trends between our non-quenched galaxies and their environment, although more work is needed to fully discard an external influence.

de Lorenzo-Cáceres, A.↗

On the origin and early evolution of biological catalysis and other studies on chemical evolution

One of the lines of research in molecular evolution which we have developed for the past three years is related to the experimental and theoretical study of the origin and early evolution of biological catalysis. In an attempt to understand the nature of the first peptidic catalysts and coenzymes, we have achieved the non-enzymatic synthesis of the coenzymes ADPG, GDPG, and CDP-ethanolamine, under conditions considered to have been prevalent on the primitive Earth. We have also accomplished the prebiotic synthesis of histidine, as well as histidyl-histidine, and we have measured the enhancing effects of this catalytic dipeptide on the dephosphorylation of deoxyribonucleotide monophosphates, the hydrolysis of oligo A, and the oligomerization 2', 3' cAMP. We reviewed and further developed the hypothesis that RNA preceded double stranded DNA molecules as a reservoir of cellular genetic information. This led us to undertake the study of extant RNA polymerases in an attempt to discover vestigial sequences preserved from early Archean times. In addition, we continued our studies of on the chemical evolution of organic compounds in the solar system and beyond.

Oro, J.↗

Coordinated analyses of hydrated Kuiper Belt dust

We report on coordinated mineralogical, chemical, and isotopic analyses of hydrated carbon-rich interplanetary dust particles (IDPs). We focus on hydrated IDPs that contain high solar flare track densities consistent with an origin in outer Solar System bodies such as Edgeworth-Kuiper belt objects (EKBOs) [1]. We are searching for evidence of the heavy water reservoir predicted by self-shielding models [2]. The IDPs were embedded in sulfur and microtome thin sections were obtained for Fourier-transform infrared (FTIR), transmission electron microscopy (TEM), and nanoSIMS analyses. The IDPs were only partly microtomed and the remainder was extracted from the sulfur and pressed into In or Au for high precision O isotopic analyses. The hydrated IDPs in this study are compact low-porosity objects with a mineralogy dominated by fine-grained saponite and lesser serpentine, Fe-Ni sulfides, Mg-Fe carbonates, rare crystaline silicates, and abundant organic matter. The C in the IDPs is abundant (typically >4X CI) and occurs as isolated nanoglobules and finely disseminated within the phyllosilicate matrix of the IDPs. TEM analyses of the nanoglobules show O-rich, S-depleted rims with higher carbonyl functionality than their cores. FTIR spectra show strong aliphatic C-H stretches consistent with the presence of short-chain aliphatic hydrocarbons in these IDPs. The oxygen isotopic compositions measured in 6 IDPs all plot within error of the CCAM line [3]. For several of the IDPs we were able to obtain multiple analyses. The ∆^(17)O values range from -23‰ up to +3.4‰ with most analyses clustering above the terrestrial fractionation line. The O isotopic compositions and the high C contents of these IDPs are distinct from known hydrated carbonaceous chondrites. The track-rich hydrated IDPs record parent body aqueous alteration effects through interactions with 16O-poor H2O on a EKBO parent body. The aqueous alteration activity likely resulted from the heat generated by collisional processes among EKBOs given that these objects are believed to have accreted late after the decay of 26Al [4].

Lindsay P. Keller↗

The Iodine Satellite (iSat) Propellant Feed System - Design and Demonstration

CUBESATS are relatively new spacecraft platforms that are typically deployed from a launch vehicle as a secondary payload, providing low-cost access to space for a wide range of end-users. These satellites are comprised of building blocks having dimensions of 10x10x10 cm3 and a mass of 1.33 kg (a 1-U size). While providing low-cost access to space, a major operational limitation is the lack of a propulsion system that can fit within a CubeSat and is capable of executing high Delta V maneuvers. This makes it difficult to use CubeSats on missions requiring certain types of maneuvers (i.e. formation flying, spacecraft rendezvous). Work has been performed investigating the use of iodine as a propellant for Hall-effect thrusters (HETs) that could subsequently be used to provide a high specific impulse path to CubeSat propulsion. One of the systems under development to support such a technology is the propellant feed system, which must be capable of storing solid iodine propellant, applying heat to sublime the stored solid into the vapor phase, and then control the flow of low-pressure gaseous iodine to both the thruster and cathode. In a test conducted in 2016, a first-generation iodine propellant feed system was integrated with a cathode and Hall thruster. While this test had to be terminated, the feed system in this first test was able to support both cathode and integrated cathode and thruster operation prior to the termination of the test. In the present paper, we describe work performed since that initial integrated test. The effort uses lessons learned from the previous integrated test, retiring risk associated with the iodine propellant feed system, answering open design-space questions, and demonstrating iodine flow control in an integrated system. The work is undertaken at both the component level and then at the integrated subsystem level to systematically improve the feed system design, improving the hardware fidelity so the appearance and operation of the system are as flight-like as possible. At the component level, the work focuses on the propellant tank, the feed system tubing, the valves used to control the flow to the cathode and thruster, and the heaters that maintain the temperature of the flowpaths and keep iodine from redepositing and clogging the system. Work on the propellant reservoir focuses on fabricating a tank that matches the geometry of the flight design, which allows for the identification of flight tank fabrication issues that may arise and permits thermal testing of a tank possessing the same size and thermal mass as the flight design, which can be used to anchor thermal modeling of the component. This is critical for finalizing the tank heater power requirements that feed into the heater design. All metallic materials in the feed system are hastelloy or Inconel, as these materials are resistant to chemical attack by the highly-reactive iodine vapor. The tubing in the iodine feed system must possess ports to permit a neutral gas purge of the system that clear impurities after iodine is loaded into the propellant tank. A procedure is discussed whereby these ports are crimped and sealed after the purge process is completed so as to not re-expose the iodine system to air. The valves are a critical component for control of the flow to the thruster and the cathode. Significant effort has gone into upgrading the materials of the valves to make them more resistant to chemical attack and into developing an understanding of the use of these valves during the startup and operation of the cathode and thruster. The heaters that line the entire feed system are designed to draw minimal power from the power processing unit (PPU) while still having the capacity to maintain all the feed system components at the temperatures required to discourage iodine deposition inside components downstream of the propellant tank exit. The heaters possess two separate resistive traces, giving the design redundancy should a failure occur in the primary heater circuit of one of the heater zones. The task of operating a feed system in conjunction with a thruster and cathode is undertaken in a series of sub-steps. The system is first assembled and operated on xenon gas, using the valves for cathode startup and thruster control based on measurement of the discharge current. After startup and control on xenon are demonstrated, the thruster will be transitioned to iodine operation, demonstrating thruster startup and feed system control while using a xenon-fed cathode. Finally, the last step is to integrate an iodine-compatible cathode with the system, demonstrate autonomous cathode start-up with open-loop control and thruster start-up with closed-loop control for multiple cycles.

Polzin, Kurt A.↗

Simplex GPS and InSAR Inversion Software

Changes in the shape of the Earth's surface can be routinely measured with precisions better than centimeters. Processes below the surface often drive these changes and as a result, investigators require models with inversion methods to characterize the sources. Simplex inverts any combination of GPS (global positioning system), UAVSAR (uninhabited aerial vehicle synthetic aperture radar), and InSAR (interferometric synthetic aperture radar) data simultaneously for elastic response from fault and fluid motions. It can be used to solve for multiple faults and parameters, all of which can be specified or allowed to vary. The software can be used to study long-term tectonic motions and the faults responsible for those motions, or can be used to invert for co-seismic slip from earthquakes. Solutions involving estimation of fault motion and changes in fluid reservoirs such as magma or water are possible. Any arbitrary number of faults or parameters can be considered. Simplex specifically solves for any of location, geometry, fault slip, and expansion/contraction of a single or multiple faults. It inverts GPS and InSAR data for elastic dislocations in a half-space. Slip parameters include strike slip, dip slip, and tensile dislocations. It includes a map interface for both setting up the models and viewing the results. Results, including faults, and observed, computed, and residual displacements, are output in text format, a map interface, and can be exported to KML. The software interfaces with the QuakeTables database allowing a user to select existing fault parameters or data. Simplex can be accessed through the QuakeSim portal graphical user interface or run from a UNIX command line.

Donnellan, Andrea↗

Comparison of the earliest NC and CC planetesimals: Evidence from ungrouped iron meteorites

Isotope anomalies in meteorites reveal a fundamental dichotomy between Non-Carbonaceous- (NC) and Carbonaceous-type (CC) planetary bodies. Until now, this dichotomy is established for the major meteorite groups, representing about 36 distinct parent bodies. Ungrouped meteorites represent an even larger number of additional parent bodies, but whether they conform to the overall NC-CC dichotomy is unknown. Here, the genetics and chronology of 26 ungrouped iron meteorites are considered through nucleosynthetic Mo and radiogenic W isotopic compositions. Secondary cosmic ray-induced modifications of these isotope compositions are corrected using Pt isotope measurements on the same samples. We find that all of the ungrouped irons have Mo isotope anomalies within the range of the major meteorite groups and confirm the NC-CC dichotomy for Mo, where NC and CC meteorites define two distinct, subparallel s-process mixing lines. All ungrouped NC irons fall on the NC-line, which is now precisely defined for 41 distinct parent bodies. The ungrouped CC irons show scatter around the CC-line indicative of small r-process Mo heterogeneities among these samples. These r-process Mo isotope variations correlate with O isotope anomalies, most likely reflecting mixing of CI chondrite-like matrix, chondrule precursors, and Ca-Al-rich inclusions. This implies that CC iron meteorite parent bodies accreted the same nebular components as the later-formed carbonaceous chondrites. The Hf-W model ages of core formation for the ungrouped irons overlap with those of the iron meteorite groups from each reservoir and reveal a narrow age peak at ∼3.3 Ma after Ca-Al-rich inclusions for the CC irons. By contrast, the NC irons display more variable ages, including younger ages indicative of impact-induced melting events, which seem absent among the CC irons. This is attributed to the more fragile and porous nature of the CC bodies, making impact-induced melting on their surfaces difficult. The chemical characteristics of all iron meteorites together reveal slightly more oxidizing conditions during core formation for CC compared to NC irons. More strikingly, strong depletions in moderately volatile elements, typical of many iron meteorite parent bodies, predominantly occur among CC irons, for reasons that remain unclear at present.

Chondrite components↗

Online monitoring of lanthanide species with combined spectroscopy in flowing aqueous aerosol systems

Combined spectroscopic analysis through absorption and laser-induced breakdown spectroscopy (LIBS) was used to monitor Nd and Pr concentrations in a flowing aqueous system. A unique, online sampling approach was employed that allowed for an approximately closed analysis loop of liquids from a reservoir. Absorption spectroscopy was performed with an optical flow cell, and LIBS was performed on an aerosol stream. Multivariate calibrations based on combined absorption and LIBS signals were built for Nd and Pr and then used to monitor concentrations in mixed solutions in a series of spiking tests. In these tests, the concentrations of Nd and Pr in solution were intermittently changed while spectroscopic signals were monitored in real-time. The combined spectroscopic signals and multivariate models were successful in monitoring changing concentrations of lanthanide species with high accuracy and minimal latency. Root-mean squared error of predictions were 0.015 mol L −1 and 0.019 mol L −1 for Nd and Pr respectively, and these lanthanides were able to be monitored to an accuracy of < 0.2 wt%. This work demonstrates both the capabilities of a sampling system for online analysis of liquids and the capabilities of multimodal spectroscopic characterization for real-time, continuous tracking of species in potentially hazardous liquid systems.

absorbance spectroscopy↗

Kimberlina 1.2 CCUS Geophysical Models and Synthetic Data Sets

This synthetic multi-scale and multi-physics data set was produced in collaboration with teams at the Lawrence Berkeley National Laboratory, National Energy Technology Laboratory, Los Alamos National Laboratory, and Colorado School of Mines through the Science-informed Machine Learning for Accelerating Real-Time Decisions in Subsurface Applications (SMART) Initiative. Data are associated with the following publication: Alumbaugh, D., Gasperikova, E., Crandall, D., Commer, M., Feng, S., Harbert, W., Li, Y., Lin, Y., and Samarasinghe, S., “The Kimberlina Synthetic Geophysical Model and Data Set for CO2 Monitoring Investigations”, The Geoscience Data Journal, 2023, DOI: 10.1002/gdj3.191. The dataset uses the Kimberlina 1.2 CO2 reservoir flow model simulations based on a hypothetical CO2 storage site in California (Birkholzer et al., 2011; Wainwright et al., 2013). Geophysical properties models (P- and S-wave seismic velocities, saturated density, and electrical resistivity) were produced with an approach similar to that of Yang et al. (2019) and Gasperikova et al. (2022) for 100 Kimberlina 1.2 reservoir models. Links to individual resources are provided below: [CO2 Saturation Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-co2-saturation-models); Resistivity Models – [part 1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-1), [part 2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-2), and [part 3](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-3); [Vp Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vp-velocity-models); [Vs Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vs-velocity-models); [Density Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-density-models). The 3D distributions of geophysical properties for the 33 time stamps of the SIM001 model were used to generate synthetic seismic, gravity, and electromagnetic (EM) responses for 33 times between zero and 200 years. Synthetic surface seismic data were generated using 2D and 3D finite-difference codes that simulate the acoustic wave equation (Moczo et al., 2007). 2D data were simulated for six point-pressure sources along a 2D line with 10 m receiver spacing and a time spacing of 0.0005 s. 3D simulations were completed for 25 surface pressure sources using a source separation of 1 km in both the x and y directions and a time spacing of 0.001 s. Links to individual resources are provided below: [2D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-velocity-models) and [2D surface seismic data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-surface-seismic-data). [3D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-velocity-models), and 3D seismic data [year0](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year0), [year1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year1), [year2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year2), [year5](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year5), [year10](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year10), [year15](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year15), [year20](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year20), [year25](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year25), [year30](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year30), [year35](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year35), [year40](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year40), [year45](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year45), [year49](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year49), [year50](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year50), [year51](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year51), [year52](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year52), [year55](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year55), [year60](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year60), [year65](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year65), [year70](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year70), [year75](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year75), [year80](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year80), [year85](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year85), [year90](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year90), [year95](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year95), [year100](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year100), [year110](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year110), [year120](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year120), [year130](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year130), [year140](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year140), [year150](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year150), [year175](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year175), [year200](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year200). The Python scripts to read these models and data are provided [here](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-python-scripts). EM simulations used a borehole-to-surface survey configuration, with the source located near the reservoir level and receivers on the surface using the code developed by Commer and Newman (2008). Pseudo-2D data for the source at [2500 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz2500m) and [3025 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz3025m), used a 2D inline receiver configuration to simulate a response over 3D resistivity models. The [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-csem-data) contain electric fields generated by borehole sources at monitoring well locations and measured over a surface receiver grid. Vector gravity data, both on the surface and in boreholes, were simulated using a modeling code developed by Rim and Li (2015). The simulation scenarios were parallel to those used for the EM: [pseudo-2D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were calculated along the same lines and within the same boreholes, and [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were simulated over 3D models on the surface and in three monitoring wells. A series of [synthetic well logs](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-well-logs) of CO2 saturation, acoustic velocity, density, and induction resistivity in the injection well and three monitoring wells are also provided at 0, 1, 2, 5, 10, 15, and 20 years after the initiation of injection. These were constructed by combining the low-frequency trend of the geophysical models with the high-frequency variations of actual well logs collected in the Kimberlina 1 well that was drilled at the proposed site. Measurements of permeability and pore connectivity were made on cores of Vedder Sandstone, which forms the primary reservoir unit: [CT micro scans](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-ct-micro-scans-of-vedder-formation) and [Industrial CT Images](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-industrial-ct-images-vedder-formation). These measurements provide the range of scales in the otherwise synthetic data set to be as close to a real-world situation as possible. References: Birkholzer, J.T., Zhou, Q., Cortis, A. and Finsterle, S., 2011. A sensitivity study on regional pressure buildup from large-scale CO2 storage projects. Energy Procedia, 4, 4371-4378. Commer, M., and Newman, G.A., 2008. New advances in three-dimensional controlled-source electromagnetic inversion, Geophysical Journal International, 172, 513-535. Gasperikova, E., Appriou, D., Bonneville, A., Feng, Z., Huang, L., Gao, K., Yang, X., Daley, T., 2022, Sensitivity of geophysical techniques for monitoring secondary CO2 storage plumes, Int. J. Greenh. Gas Control, Volume 114, 103585, ISSN 1750-5836, https://doi.org/10.1016/j.ijggc.2022.103585. Moczo, P., J.O. Robertsson and L. Eisner, 2007, The finite-difference time-domain method for modeling of seismic wave propagation: Advances in geophysics, 48, 421-516. Rim, H., and Y. Li, 2015, Advantages of borehole vector gravity in density imaging, Geophysics, 80, G1-G13. Wainwright, H. M.; Finsterle, S.; Zhou, Q.; Birkholzer, J. T., 2013. Modeling the Performance of Large-Scale CO2 Storage Systems: A Comparison of Different Sensitivity Analysis Methods. International Journal of Greenhouse Gas Control, 17, 189205. https://doi.org/10.1016/j.ijggc.2013.05.007, DOI: 10.18141/1603331. Yang, X., Buscheck, T.A., Mansoor, K., Wang, Z., Gao, K., Huang, L., Appriou, D., and Carroll, S.A., 2019. Assessment of geophysical monitoring methods for detection of brine and CO2 leakage in drinking water aquifers, International Journal of Greenhouse Gas Control, 90, 102803, https://doi.org/10.1016/j.ijggc.2019.102803.

CCUS↗

Precision Oxygen Isotope Measurements of Two C-Rich Hydrated Interplanetary Dust Particles

Introduction: Chondritic-smooth IDPs (Interplanetary Dust Particles) are low porosity objects whose mineralogy is dominated by aqueous alteration products such as Mg-rich phyllosilicates (smectite and serpentine group) and Mg-Fe carbonate minerals. Their hydrated mineralogy combined with low atmospheric entry velocities have been used to infer an origin largely from asteroidal sources. Spectroscopic studies show that the types and abundance of organic matter in CS IDPs is similar to that in CP IDPs. Although CS IDPs show broad similarities to primitive carbonaceous chondrites, only a few particles have been directly linked to specific meteorite groups such as CM and CI chondrites based on the presence of diagnostic minerals. Many CS IDPs however, have carbon contents that greatly exceed that of known meteorite groups suggesting that they either may derive from comets or represent samples of more primitive parent bodies than do meteorites. It is now recognized that many large, dark primitive asteroids in the outer main belt, as well as some trans-Neptunian objects, show spectroscopic evidence for aqueous alteration products on their surfaces. Some CS IDPs exhibit large bulk D enrichments similar to those observed in the cometary CP IDPs. While hydrated minerals in comets have not been unambiguously identified to date, the presence of the smectite group mineral nontronite has been inferred from infrared spectra obtained from the ejecta from comet 9P/Tempel 1 during the Deep Impact mission. Recent observations of low temperature sulfide minerals in Stardust mission samples suggest that limited aqueous activity occurred on comet Wild-2. All of these observations, taken together, suggest that the high-carbon hydrated IDPs are abundant and important samples of primitive solar system objects not represented in meteorite collections. Oxygen isotopic compositions of chondrites reflect mixing between a 16O-rich reservoir and a 17O,18O-rich reservoir produced via mass-independent fractionation. The composition of the 16O-rich reservoir is well constrained but material representing the 17O,18O-rich end-member is rare. Self-shielding models predict that cometary water, presumed to represent this reservoir, should be enriched in 17O and 18O by greater than 200 per mille. The high-carbon hydrated IDPs may be among the best materials available to search for preserved "cometary" H2O signatures. In order to better understand the origin and evolution of these particles, we have obtained 10 hydrated interplanetary dust particles for coordinated mineralogical, isotopic and organic analyses. We have previously reported the results of mineralogical and O isotopic measurements of two hydrated IDPs; here we present results of O isotopic measurements of three additional IDPs. Samples and Methods: Three interplanetary dust particles (L2079C35, L2083D46 and L2083E46) were embedded in S and partially ultramicrotomed into approximately 70 nanometer sections for analysis via transmission electron microscopy (TEM). The remainders of the unsliced particles were removed from S and pressed into high purity Au foil that was cleaned with HF acid and annealed at 800 degrees Centigrade. The pressed IDPs were analyzed via electron microprobe analysis (EPMA) for quantitative bulk chemical analysis. After EPMA analysis, the IDPs were subjected to precision O isotope analysis with the UCLA Cameca IMS-1270 ion probe. A 20 kiloelectronvolt, 0.5 nanoangstrom Cs+ primary beam of approximately 15 micrometers diameter was used for each measurement. Small particles of San Carlos olivine and Burma spinel were pressed into the Au foil for use as standards to correct for instrumental mass fractionation. The detection system was configured for multicollection, with 16O measured on a Faraday cup, and 17O and 18O measured on electron multipliers (EMs). Individual analyses consisted of 15 cycles of 10 seconds per cycle. Additionally, two microtome thin sections were measured for H isotopic compositions with the JSC NanoSIMS 50L ion probe. An 8 picoangstrom, 16 kiloelectronvolt Cs plus primary beam was used. Measurements consisted of H, C, 12C, 16O, and 18O collected with EMs in multicollection. Terrestrial biotite and kerogen were used for isotopic standards. A significant challenge in O isotope measurement of hydrated minerals is the interference from 16OH at mass 17O. We ensured that the 17O and 16OH peaks were fully resolved by using a mass resolution of greater than 7000 and by careful analyses of San Carlos olivine, Burma spinel and chlorite hydrated mineral standards. The hydride was further suppressed with a cold finger attached to an LN2 dewar to trap volatiles in the sample chamber. All sputtered ions were counted (i.e. presputtering was not used); after applying background, yield and deadtime corrections, we performed a change-point analysis on our data via R in order to determine when the sample reached sputtering equilibrium; data points collected prior to the change point were excluded. Change-point analysis was also used to determine whether the IDP had completely sputtered. Results: Mineralogy. IDPs C35 and E46 exhibited hydrated mineralogies, Fe-Ni sulfide grains, nanoglobules and occasional enstatite grains distributed throughout a fine-grained Mg-Fe saponite matrix. C35 also contained breunnerite (Mg,Fe)CO3; solar flare tracks were observed in enstatite, indicating minimal atmospheric entry heating. The mineralogy of D46 is dominated by a large FeS grain with a minor component of adhering silicate material. D46 was strongly heated during atmorpheric entry as evidenced by a well-developed magnetite rim. EPMA analyses show that both C35 and E46 have high carbon contents of 20 weight percentage (approximately 6X CI). D46 contains approximately 6 wt.weight percentage C. A significant amount of the carbon is present as carbon nanoglobules. Results: Hydrogen isotopes: Although the bulk delta D values of both sections of L2079C35 were within error of SMOW (minus 33 plus or minus 19 per mille, 1 plus or minus 14 per mille 1 sigma), several delta D-rich hotspots were also identified, reaching 2000 per mille. As shown in Fig. 1, these hotspots are clearly associated with discrete carbonaceous inclusions that are akin to nanoglobules found in many meteorites and other IDPs. Oxygen Isotopes. Results of the oxygen isotope measurements are shown in Figure 1. The oxygen isotope composition for L2079C35 was delta 18O equals plus11.6 plus or minus 1.9per mille, delta 17O equals plus 7.9 plus or minus 1.9per mille (2 standard errors). The oxygen isotope composition for L2083D46 was delta 18O equals minus 8.1 plus or minus 1.9 per mille, delta 17O equals minus 6.4 plus or minus 3.1 per mille (2 standard errors). The oxygen isotope composition for L2083E46 was delta 18O equals plus 12.0 plus or minus 1.9 per mille, delta 17O equals plus 9.2 plus or minus 2.0 per mille (2 standard errors). Discussion: Despite mineralogical similarities to highly aqueously altered carbonaceous chondrites, the hydrated IDPs we analyzed have oxygen isotopic compositions that are distinct from matrix materials in the CI, CM, and CR chondrites. The IDPs plot along the Young-Russell line, with delta 17O values for C35 and E46 suggestive of interaction with a 16O-poor reservoir. However, we have thus far not observed evidence of extreme 16O-poor reservoirs expected from self-shielding models and observed in Acfer 094 simplectite. The high carbon contents of the IDPs also set them apart from known meteoritic samples. The lack of atmospheric entry heating effects are consistent with low encounter velocities and suggest either an asteroidal source, or a low inclination, low eccentricity cometary origin. Conclusions: The unusual oxygen isotopic compositions, high carbon contents, and the abundance of Drich nanoglobules, together, suggest that the high-carbon, hydrated IDPs are derived from a primitive source that is not yet represented in meteorite collections.

Snead, C. J.↗

SUBTASK 1.6 – BASIN ELECTRIC CARBON STORAGE RESEARCH PROJECT: NOVEL MONITORING TECHNIQUES

The Energy & Environmental Research Center (EERC) conducted baseline activities associated with an applied research project at Basin Electric Power Cooperative’s (Basin’s) carbon capture and storage (CCS) site in Beulah, North Dakota, to establish novel carbon storage-monitoring techniques as commercial methods under Cooperative Agreement No. DE-FE0024233, Subtask 1.6. The following report summarizes the baseline activities performed and briefly describes the subsequent (operational monitoring) activities that have been proposed to the U.S. Department of Energy (DOE) as part of the overall project to develop and demonstrate novel monitoring techniques at North America’s largest permitted CCS operation. Dakota Gasification Company (DGC), a wholly owned subsidiary of Basin, owns and operates the Great Plains Synfuels Plant (GPSP) approximately 5 miles northwest of the town of Beulah, North Dakota (Figure 1). In 2023, DGC received approval from the North Dakota Industrial Commission (NDIC) to develop a storage facility on-site for injecting a stream of carbon dioxide (CO2) captured from GPSP. DGC will transport the captured CO2 stream with approximately 6.8 miles of transmission lines that extend north of GPSP and inject >1 million tonnes (MMt) of CO2 annually (>1 MMt/yr) over a 12-year period with up to six underground injection control (UIC) Class VI-compliant injection wells completed in the Broom Creek Formation, a predominantly sandstone reservoir and saline aquifer underlying GPSP. The Broom Creek Formation lies approximately 5900 feet (ft) below ground surface (bgs) at GPSP. The commercial scale (i.e., >1 MMt/yr) of DGC’s permitted carbon storage project is ideal for developing and testing the novel monitoring techniques included within Subtask 1.6. The goals of this project are to demonstrate 1) the cost-effectiveness of novel monitoring technologies included as part of this research, 2) technology capability for tracking the CO2 plume and/or associated pressure response in the subsurface and monitoring out-of-zone migration, and 3) compliance with UIC Class VI program requirements. The research activities proposed for the overall project include 1) design of an automated, integrated, modular (AIM) monitoring station; 2) time-lapse electromagnetic (EM) field surveys; 3) drone-based surveillance studies; 4) time-lapse monitoring with seismic methods; 5) advanced wellbore-monitoring methods; 6) deployment of an AIM monitoring network; 7) EM monitoring of CO2 with real-time data processing; 8) continued seasonal drone-based surveillance studies; 9) seismic monitoring with passive and active surveys; and 10) wellbore monitoring with nuclear magnetic resonance (NMR) for near-surface characterization. Completion of Activities 1.0–5.0 (baseline activities) are described in this report. Upon authorization of funding by DOE, the EERC will initiate Activities 6.0– 10.0 (operational monitoring activities). Current state-of-the-art (SOA) carbon storage-monitoring techniques require countless labor hours dedicated to the acquisition of data. Once data are gathered, these SOA techniques often rely on commercial facilities to process raw data from the field. However, it is anticipated that next-generation monitoring techniques, such as those being demonstrated, will lower acquisition footprints, be less operationally intensive, and improve data acquisition efficiencies. These new techniques are more conducive to the application of machine learning, artificial intelligence, and automation, thus providing a pathway for integration into active control systems, informing site operability, and improving the integration of data for future CCS projects across the United States. Additionally, reclaimed and active mining lands are present within the project site, creating a unique opportunity to demonstrate the effectiveness of remote sensing and surface-based geophysics monitoring techniques at similar project sites that may include disturbed, unconsolidated, or actively excavated near-surface environments. The efforts included in the overall project will produce necessary designs, learnings, and data acquired during the baseline and operational monitoring periods that are necessary for time-lapse demonstration and validation of the described monitoring techniques. In addition, it is anticipated that the monitoring technologies included in this study will be compliant with UIC Class VI requirements to enable the potential for implementation at other CCS sites across the United States.

42 ENGINEERING↗

Comet C2012 S1 (ISON): Observations of the Dust Grains From SOFIA and of the Atomic Gas From NSO Dunn and Mcmath-Pierce Solar Telescopes

Comet C/2012 S1 (ISON) is unique in that it is a dynamically new comet derived from the Oort cloud reservoir of comets with a sun-grazing orbit. Infrared (IR) and visible wavelength observing campaigns were planned on NASA's Stratospheric Observatory For Infrared Astronomy (SOFIA) and on National Solar Observatory Dunn (DST) and McMath-Pierce Solar Telescopes, respectively. We highlight our SOFIA (+FORCAST) mid- to far-IR images and spectroscopy (approx. 5-35 microns) of the dust in the coma of ISON are to be obtained by the ISON-SOFIA Team during a flight window 2013 Oct 21-23 UT (r_h approx. = 1.18 AU). Dust characteristics, identified through the 10 micron silicate emission feature and its strength, as well as spectral features from cometary crystalline silicates (Forsterite) at 11.05-11.2 microns, and near 16, 19, 23.5, 27.5, and 33 microns are compared with other Oort cloud comets that span the range of small and/or highly porous grains (e.g., C/1995 O1 (Hale-Bopp) and C/2001 Q4 (NEAT) to large and/or compact grains (e.g., C/2007 N4 (Lulin) and C/2006 P1 (McNaught)). Measurement of the crystalline peaks in contrast to the broad 10 and 20 micron amorphous silicate features yields the cometary silicate crystalline mass fraction, which is a benchmark for radial transport in our protoplanetary disk. The central wavelength positions, relative intensities, and feature asymmetries for the crystalline peaks may constrain the shapes of the crystals. Only SOFIA can look for cometary organics in the 5-8 micron region. Spatially resolved measurements of atoms and simple molecules from when comet ISON is near the Sun (r_h< 0.4 AU, near Nov-20-Dec-03 UT) were proposed for by the ISON-DST Team. Comet ISON is the first comet since comet Ikeya-Seki (1965f) suitable for studying the alkalai metals Na and K and the atoms specifically attributed to dust grains including Mg, Si, Fe, as well as Ca. DST's Horizontal Grating Spectrometer (HGS) measures 4 settings: Na I, K, C2 to sample cometary organics (along with Mg I), and [OI] as a proxy for activity from water (along with Si I and Fe I). State-of-the-art instruments that will also be employed include IBIS, which is a Fabry-Perot spectral imaging system that concurrently measures lines of Na, K, Ca II, or Fe, and ROSA (CSUN/QUB), which is a rapid imager that simultaneously monitors Ca II or CN. From McMath-Pierce, the Solar-Stellar Spectrograph also will target ISON (320-900 nm, R approx. 21,000, r_h<0.3 AU). Assuming survival, the intent is to target ISON over r_h<0.4 AU, characteristic of prior Na detections.

Oort cloud comets↗

Thermocapillary Flow and Aggregation of Bubbles on a Solid Wall

During the electrolytic evolution of oxygen bubbles forming on a vertically oriented transparent tin oxide electrode, bubbles were found to be mutually attractive. The mechanism of the aggregation had never been explained satisfactorily until Guelcher et al. attributed it to thermocapillary flow. The gradient of surface tension of the liquid at the bubble's surface, which was established because of reaction heat and ohmic heat loss at the electrode wall, drives flow of the liquid adjacent to each bubble; the bubble "pumps" fluid along its surface away from the wall. Fluid flows toward the bubble to conserve mass and entrains nearby bubbles in the flow pattern. The same logic would apply when two bubbles of equal size are adjacent to each other on a warm wall. Each bubble drives thermocapillary flow and hence entrains the other in its flow pattern, which drives the aggregation. Our objective here is to perform experiments where the temperature gradient at the wall is well known and controlled. The theory can be quantitatively tested by studying aggregation of bubble pairs of equal size, and by varying system parameters such as temperature gradient, bubble size and fluid viscosity. The results are then compared with the theory in a quantitatively rigorous manner. We demonstrate that the theory without adjustable parameters is capable of quantitatively modeling the rate of aggregation of two bubbles. The equations governing the thermocapillary flow around a single stationary bubble on a heated or cooled wall in a semi-infinite domain were solved. Both Reynolds number and Marangoni number were much less than unity. The critical result is that liquid in the vicinity of a warm wall flows toward a stationary collector bubble. Consequently the thermocapillary flow around the stationary bubble entrains another bubble toward itself. The bubbles undergo hindered translation parallel to the wall with velocity U while the fluid flow field is described with u. Two velocities were equated by using a wall hindrance parameter q: U = qu [1] which shows the velocity of bubble is proportional to the entraining velocity. The hindrance parameter q can experimentally be measured independently. q can also be calculated by solving the equations of motion for a bubble translating parallel to a solid wall. The experimental cell is cylindrical with an ID of 10 cm and consists of a 1 cm deep main cell filled with silicone oil and flanked by two thermal reservoirs. The upper thermal reservoir was heated and the lower thermal reservoir was cooled so that the bubbles aggregate. Two types of silicone oil (eta = 0.02 and 0.50 Pa s) were used. Two equal sized air bubbles were injected into the cell with a syringe. The center-to-center distance of bubbles was observed through a microscope. Bubble radius ranged from 0.40 mm to 0.65 mm and the temperature gradients along with the cell ranged from 1400 to 5000 K/m. The bubbles aggregated when heat flows from the wall to the fluid. The velocities of bubbles were in the range of 1 - 10 microns/s. The separation r decreased more quickly when the temperature gradient was higher, bubble size was larger, and the oil viscosity was lower. r decreased more rapidly as the bubbles approached each other. Dimensionless time was arbitrarily set to be zero when the dimensionless center-to-center distance between the bubbles was 4. All the bubble trajectories fall onto one line, especially in the range of dimensionless distance from 4 to 3. This means the relative movement of the bubble pair is proportional to the temperature gradient and bubble size and it is inversely proportional to the viscosity of the oil. This result strongly suggests that the thermocapillary flow-based aggregation mechanism is correct. A value of q can be estimated by fitting the scaled data to Eq. [1]. A best fit value of q was obtained as q = 0.26 with a standard deviation of 0.03. Independent experimental results for q for a 0.5 mm radius bubble, give values of q in the range 0.11 to 0.23. The value of q obtained from solving the equations of motion reveals q has values in the range 0.23 - 0.30. Since the full scale of possible values of q is zero to one, the maximum deviation of independently determined values of q from the best fit value was 15% of this full scale. Thus reasonable quantitative agreement between theory and experiment has been obtained.

Kasumi, Hiroki↗

Heterogeneity of Bulk Oxygen Isotopic Compositions in Anhydrous Interplanetary Dust Particles

Introduction: Anhydrous interplanetary dust particles are one of the least altered ancient solar systemmaterials. Their bulk chemical compositions match those of CI chondrites within a factor of 2-3, except for carbonwhich is enriched in both anhydrous and hydrated IDPs by ~4x CI [1]. Hydrous IDPs often show 16O-poor isotopiccompositions, likely due to the interaction with isotopically heavy H2O [2]. Anhydrous IDPs are interpreted tooriginate from comets from the outer solar system [3]. Oxygen isotopic data for these materials is scarce, but resultspoint to a wider range of compositions [4]. This heterogeneity might be a result of mixing between a 16O-rich and a16O-poor reservoir formed via self-shielding and photodissociation of CO [5]. Here, we report bulk oxygen isotopiccompositions of three anhydrous IDPs in context of their mineralogical composition and discuss possible origins ofisotopic heterogeneity. Experimental: Bulk chemical composition and mineralogy of 70 nm thin ultramicrotomed sections of threeanhydrous IDPs (L2099A7, L2099A8, and L2071AB1) were analyzed using a Thermo Scientific Titan Themis G360-300 TEM equipped with a four-quadrant energy-dispersive X-ray detector (Super-X G2) at University ofMünster. Bulk oxygen isotopic compositions were measured using a NanoSIMS 50 at the Max-Planck-Institut fürChemie in Mainz. A focused Cs+ ion beam (~1 pA, ~100 nm) was rastered over a field of view varying from 5 x 5 to6 x 6 μm2, depending on particle size for 15 layers. Negative secondary ions of 16O, 17O, 18O, 12C14N, and 28Si werecollected on electron multipliers. The 16OH contribution to the 17O peak was 1-2 ‰. As a standard, a matrix regionof meteorite CR2 Queen Alexandra Range (QUE) 99177 was used, whose oxygen isotopic composition wasmeasured by [6]. Results and Discussion: Particles A7 and A8 are both fine-grained, consisting mostly of small, 100 nm-sizedequilibrated aggregates (EA). Both particles exhibit discontinuous magnetite rims with thicknesses up to 50 nm inA7 and 100 nm in A8 on the outsides. In Particle AB1, the upper part of the particle consists of a few EAs, while thelower part contains lots of small GEMS grains with 100-200 nm diameter clustered together. Magnetite rims arenearly absent, the rare rims reaching maximum 20 nm thickness, indicating that AB1 is of more primitive naturethan the other two IDPs [1]. Furthermore, it contains several diffuse GEMS-like areas, identified by an amorphoussilicate groundmass with small nano-inclusions of Fe,Ni-metal and Fe-sulfides. All three particles are subsolar forall major element/Si ratios. Mean S/Si of particles A7 and A8 is 0.046 and 0.048, respectively, while it is an order ofmagnitude higher in AB1 (0.184). This reflects the higher degree of thermal alteration in A7 and A8, resulting in theloss of volatile S and oxidation of Fe-sulfides and FeNi-metal to magnetite as present in particle rims. A8 is rather16O-poor with δ17OSMOW = 7.8 ± 3.4 ‰ and δ18OSMOW = 9.5 ± 3.1 ‰ (1σ), followed by A7 with δ17OSMOW = − 0.8 ±4.6 ‰ and δ18OSMOW = −2.2 ± 3.4 ‰. In contrast to that, AB1 has the most 16O-rich bulk composition with δ17OSMOW= −24.2 ± 5.4 ‰ and δ18OSMOW = −25.3 ± 3.6 ‰. All investigated IDPs plot slightly above the CCAM line aspreviously reported for other anhydrous IDPs [4], maybe due to a contribution from circumstellar dust from AGBstars enriched in 17O [7]. The 16O-rich composition of AB1 is similar to the most 16O-rich anhydrous IDP U2015D21measured by [8] which is dominated by GEMS [9] and the GEMS-rich IDP GM4-2 reported by [10]. Nevertheless,it is unlikely that GEMS are the carrier for the 16O-rich composition of AB1 because most GEMS have O isotopiccompositions indistinguishable from terrestrial values, although errors on these analyses are large [11]. The isotopicheterogeneity is best explained by contribution of grains from different regions in the protoplanetary disk, samplingdifferent O isotope reservoirs, created by self-shielding and photodissociation of CO. This process resulted in 16O-rich CO gas and 16O-poor H2O that froze as ice-mantles onto dust grains [12]. Alternatively, the higher degree ofheating in A7 and A8 could have influenced isotopic fractionation, because higher δ17,18O values are reported fromthe regions of the atmosphere where small particles are heated and oxidized, as observed for the thermal alteration ofcosmic spherules [13]. Acknowledgements: We would like to thank NASA Astromaterials Acquisition & Curation Office for providing IDPsamples and DFG for funding this project (VO1816/5-1)

B Schulz↗

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.↗