SEARCH · Engineering Papers

Results for “reference concentration”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,027 records · Page 57

All sources NASA NTRS DOE OSTI

NASA NTRS2004

Sensor for Monitoring Nanodevice-Fabrication Plasmas

The term plasma process diagnostics (PPD) refers to a spectroscopic technique and sensing hardware that have been proposed for monitoring plasma processes used to fabricate electronic devices that feature sizes as small as several nanometers. Nanometer dimensions are characteristic of the quantum level of miniaturization, where single impurity atoms or molecules can drastically change the local properties of the nanostructures. Such changes may be purposely used in nanoscale design but may also be extremely damaging or cause improper operation of the fabricated devices. Determination of temperature and densities of reactants near the developing features is important, since the structural synthesis is affected by characteristics of the local microenvironment. Consequently, sensors capable of nonintrusive monitoring with high sensitivity and high resolution are essential for real-time atomistic control of reaction kinetics and minimizing trace contamination in plasma processes used to fabricate electronic nanodevices. Such process-monitoring sensors are required to be compact, multiparametric, and immune to the harsh environments of processing plasmas. PPD is intended to satisfy these requirements. The specific technique used to implement plasma diagnostics with a PPD sensor would be an advanced version of continuous-wave cavity-ringdown spectroscopy (CW-CRDS) capable of profiling spectral line broadenings in order to derive both Doppler and Stark components. CRDS is based on measurements of the rate of absorption of laser light in an optical resonator. The ultimate sensitivity results from a very long absorption path length within the cavity and immunity to variations in incident laser intensity. The proposed version of this technique would involve the use of multiplexing tunable laser diodes and an actively modulated high-reflectivity optical resonator, thus offering a synergistic combination of simplicity, compactness, high sensitivity, and high resolution. The multiplexing capabilities of diode lasers could be utilized to make the PPD sensor a single, simple, compact, and inexpensive tool for the acquisition of multiparametric data. A PPD sensor would be capable of continuous measurement of such physical parameters as gas temperature, gas velocity, electron number density, and absolute densities of reacting chemical species. A laser beam can be easily adjusted to analyze the immediate vicinity of the growing nanostructures (or features etched down) in real time. The absorption enhancement in an optical cavity would afford the sensitivity needed for measurement of the temperature and densities of species at concentrations significantly lower than measurable by other nonintrusive techniques. It is anticipated that fully developed PPD sensors would enable simultaneous measurement of local temperature and determination of plasma species responsible for the synthesis and functionalization of nanodevices. These sensors would also enable tracking the pathways and origins of damaging contaminants, thereby providing feedback for adjustment of processes to optimize them and reduce contamination. The PPD sensors should also be useful for optimization of conventional microelectronics manufacturing plasma processes. Going beyond plasma processes for fabrication of electronic devices, PPD sensors could be used for monitoring of atoms, molecules, ions, radicals, clusters, and particles in a variety of other settings, including outer space. Because of their high sensitivity, such sensors could also prove useful for detecting traces of illegal drugs and explosives.

Bolshakov, Alexander↗

DOE OSTI2024

Voltage Optimization

Voltage optimization refers to a volt-var optimization technique which was originally designed to minimize energy consumption and improve end-use efficiency on the distribution system. This reduces source voltage at a substation, which lowers generation demand, reduces system losses, and improves system stability. This has been the conventional implementation of voltage optimization: focusing on regulation of voltages throughout the distribution system via coordinated adjustments of load tap changers, line-voltage regulators, switched shunt compensation, and the targeted placement of new shunt compensation. These methods have received increased focus following a growth in distributed-energy resources (DERs) interconnecting at the distribution level, given their developing capability to participate in voltage regulation. This is paired with the intermittency challenges they pose, load imbalances, and voltage fluctuations. Voltage optimization has long been an area of extensive research within the distribution system, with various strategies and philosophies demonstrating a range of effects. The conventional focus on the distribution system owes to the typical operation of transmission systems, which assumes firm local resources and flows strictly from generation resources to loads. The strategies employed in distribution have also been optimized for this typically radial configuration and operation of the distribution system.

24 - POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRSDate not supplied

Additional Trials & Tribulations in Synthesizing a Sulfide Standard

Well-characterized sulfide reference materials that can serve as matrix-matched calibrants for in-situ trace element analyses via laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) remain elusive. Here, we describe the creation an in-house sulfide standard at NASA JSC that will be used to measure siderophile and chalcophile trace elements in low pressure experimental products, specifically pentlandite ([FeNi]9S8) and pyrrhotite ([FeNi]1-xS), as well as natural sulfides in terrestrial and meteoritic mineral assemblages. Recent methods in creating a homogenous sulfide standard include pressed pellets, synthesized chips, or fused glass [e.g. 1-3]. Other studies, however, have had success in forming homogenous FeS standards via synthesization methods [e.g. 4-6]. We present a variation of the method described in [4], in which we create a pyrrhotite standard doped with a variety of trace elements (Zn, As, Se, Mo, Ru, Rh, Pd, Sn, Sb, Te, W, Os, Ir, Pt, Au) at ~10-40 ppm and Cu at ~200 ppm. Sulfide compositions were constructed using Fe and Ni metal and elemental sulfur powders. Trace elements were added to elemental sulfur from Atomic Absorption Spectroscopy (AAS) elemental standards as either nitrate or chloride solutions to prevent oxidation of the metal powders. The dried sulfur mixture was mixed with Fe and Ni powders and mechanically mixed before two aliquots were placed in separate SiO2 tubes. Each tube was held under vacuum for ~30 minutes, sealed under vacuum, and then heated at 800C for 48 hours. Like [4], the synthesis products were composed of pourous sulfide crystals. Major element analyses of both experimental aliquots, collected using an Electron Probe Microanalyser (EPMA), yield indistinguishable major element compositions (uncertainties in 2SE & 2RSE[%]), with an average of 57.90  0.09 (0.16 %), 4.95  0.03 (0.55 %), and 39.70  0.14 (0.34 %), for Fe, Ni, and S, respectively. Trace element data were measured using a Photon Machines 193nm laser ablation system coupled to a Thermo-Scientific Element-XR ICP-MS. Spot sizes were limited to 50 μm due to the porous nature of the sulfide target material. Trace element abundances, normalized to Fe as an internal standard, were also found to be homogenous between the two aliquots, with weighted mean 2RSE (%) values of <3.0 for all trace elements. Synthesized products were re-powdered and absolute concentrations measured via solution ICP-MS. Although the sulfide appears to be homogenous, sintering experiments will be performed to more closely match the standard density to natural sulfides and minimize differences in ablation behavior. Improved density also allows for higher sensitivity (i.e. more compact target material) and larger spot sizes or traverses, as void space is eliminated. Additional major and trace element analyses on the products of the sintering experiments will be undertaken.

Jacob B Setera↗

DOE OSTI2022

IM3/HyperFACETS Thermodynamic Global Warming (TGW) Simulation Datasets

Publication For a thorough description of the methods, see the peer-reviewed paper: Jones, A.D., Rastogi, D., Vahmani, P. et al. Continental United States climate projections based on thermodynamic modification of historical weather. Sci Data 10, 664 (2023). https://doi.org/10.1038/s41597-023-02485-5 Overview The IM3 / HyperFACETS climate simulations provide 40-year historical (1980-2019) as well as four 80-year future simulations (2020-2099) over the U.S. The future simulations are split into near (2020-2059) and far future (2060-2099) segments. The future scenarios span a range of plausible changes in future climate (both Global Circulation Model (GCM) and Representative Concentration Pathways/Shared Socioeconomic Pathway (RCP/SSP) dimensions). The simulations provide climate variables with high spatiotemporal resolution (25 hourly variables and 207 3-hourly variables at 12 km2). The datasets are generated using dynamical downscaling with the WRF (Weather Research and Forecasting) model (version 4.2.1) and therefore preserve physical consistency across variables. WRF is a state-of-the-art, fully compressible, non-hydrostatic, mesoscale numerical weather prediction model. WRF is coupled with an urban canopy model (UCM), which resolves urban surfaces. The future scenarios were developed using a thermodynamic global warming approach where past events are replayed under a range of future warming conditions. These scenarios therefore provide a perspective on potential increases in extreme event intensity, geographic scope, and duration, with previously non-extreme conditions potentially crossing new thresholds to be considered extreme by today's standards. This approach is not intended to estimate future changes in extreme event frequency that might result from changes in large-scale atmospheric dynamics. This dataset has NOT been bias corrected. A bias corrected version of selected variables is under development and will be released here when available. Scenarios Files Data for each scenario is provided in weekly NetCDF files. 25 variables are available at hourly resolution, and 207 variables are available at three-hourly resolution. Spatial resolution is 12km and spans the conterminous United States (CONUS), including some areas of Canada and Mexico, resulting in a grid of 424 by 299 cells. The spatial projection is a Lambert Conformal Conic with the following proj-string: "+proj=lcc +lat_0=40.0000076293945 +lon_0=-97 +lat_1=30 +lat_2=45 +x_0=0 +y_0=0 +R=6370000 +units=m +no_defs". The available scenarios and simulation periods are listed below: historical | 1980 - 2019 rcp45cooler | 2020 - 2059 rcp45cooler | 2060 - 2099 rcp45hotter | 2020 - 2059 rcp45hotter | 2060 - 2099 rcp85cooler | 2020 - 2059 rcp85cooler | 2060 - 2099 rcp85hotter | 2020 - 2059 rcp85hotter | 2060 - 2099 * The first year (1979, 2019, and 2059) of data within each scenario represents a model warmup period and should not be used. These are located in the `spinup_files` directory. Historical year 2020 is considered an extra year of data beyond the simulation period and can be found in the `additional_files` directory. For information on specific variables and a more in-depth discussion of methodology, please refer to the data landing page at https://tgw-data.msdlive.org. Delta Warming Files The global and CONUS warming deltas for each scenario are provided in degrees Celsius annually and monthly. Restart Files Yearly restart files are provided for each scenario which can be used to restart the WRF model at a particular point in time. Spinup Files The first year of data within each simulation period represents a model warmup period and should not be used. The files are provided here for the sake of reproducibility. Additional Files Additional years of data are provided as an extension of the historic simulation.

Jones, Andrew D.↗

NASA NTRS1993

Composition and maturity of the 60013/14 core

The 60013/14 double drive tube (62 cm deep) is one of three regolith cores taken 35-40 m apart in a triangular array on the Cayley plains at station 10' (LM/ALSEP), Apollo 16. This trio, which includes double drive tube 60009/10 (59 cm deep) and deep drill core 60001-7 (220 cm), is the only such array of cores returned from the Moon. The top 45 cm of 60013/14 is mature, as is surface reference soil 60601 taken nearby. Maturity generally decreases with depth, with soil below 45 cm being submature. The zone of lowest maturity (34 is less than or equal to I(sub s)/FeO is less than 50) extends from 46 to 58 cm depth, and corresponds to the distinct region of light-colored soil observed during core processing. In the other two cores, most of the compositional variation results from mixing between fine-grained, mature soil with 10-11 micro-g/g Sc and coarse-grained ferroan anorthosite consisting of greater than 99% plagioclase with less than 0.5 micro-g/g Sc. This is most evident in 60009/10 which contains a high abundance of plagioclase at about 54 cm depth (minimum Sc: 3-4 micro-g/g); a similar zone occurs in 60001-7 at 17-22 cm (MPU-C), although it is not as rich in plagioclase (minimum Sc: 6-7 micro-g/g). Compositional variations are less in 60013/14 than in the other two cores (range: 7.9-10.0 micro-g/g Sc), but are generally consistent with the 'plagioclase dilution' effect seen in 60009/10, i.e., most 60013/14 samples plot along the mixing line of 60009/10. However, a plagioclase component is not the cause of the lower maturity and lighter color of the unit at 46-58 cm depth in 60013/14. Many of the samples in this zone have distinctly lower Sm/Sc ratios than typical LM-area soils and plot off the mixing trend defined by 60009/10. This requires a component with moderately high Sc, but low-Sm/Sc, such as feldspathic fragmental breccia (FFB) or granulitic breccia. A component of Descartes regolith, such as occurs at North Ray Crater (NRC) and which is rich in FFB, could account for the composition of these soils (i.e., a 3:1 mixture of 60601 and NRC soil). It seems unlikely that NRC ejecta would occur half a meter deep at the LM station, thus this low-Sm/Sc component may result from an older, local crater that penetrated the Cayley surface layer and excavated underlying Descartes material, as did North Ray Crater. There is no evidence for such a unit or component in the other two cores. Soil below the light-colored unit (58-62) cm has 'typical' Sm/Sc ratios, but the lowest absolute Sc concentrations, i.e., it is compositionally equivalent to a mixture of surface soil and plagioclase such as that in ferroan anorthosite. This is the only soil that might be related to the plagioclase-rich units in the other two cores. Except for the mature soil at the top of each core and, perhaps, the plagioclase-rich layers, there is little compositional evidence for any common unit among the three cores. Soil corresponding to the mare-glass-bearing unit (MPU-B) and regolith-breccia-bearing unit (MPU-A) of 60001-7 do not occur in 60013/14 or 60009/10.

Korotev, Randy L.↗

NASA NTRSDate not supplied

Enumeration and Fluorescence In Situ Hybridization of Microbial Bioburden on Cleanroom Surfaces

Introduction: Microorganisms are everywhere on Earth, even in the cleanest of places. Spacecraft assembly cleanrooms can harbor low levels of living and dead microbial cells (e.g., [1,2]), and cleanroom bioburden can also include organic molecules from industrial sources and in situ biomass. Life detection missions require careful attention to avoid contaminants that can be easily convoluted with analytical targets. We are evaluating epifluorescent microscopy and fluorescence in situ hybridization (FISH) as methods to complement organic contamination detection techniques. Epifluorescent cell counting offers an accurate and cost-effective way to quantify low levels of surface biomass. FISH could allow for the identification of residual organisms, and can be targeted to detect active populations of specific organisms such as bacteria known to resist cleaning procedures. This effort is part of a larger study that is concentrated on characterizing the surface and airborne molecular organic contamination background in Johnson Space Center (JSC) Astromaterials curation laboratories and Goddard Space Flight Center (GSFC) spacecraft assembly rooms, and understanding contaminants in the context of cleaning procedures and residual bioburden. Methods: Samples were collected by swabbing surfaces in ISO 5 and ISO 7 equivalent cleanrooms at JSC. Swabs for FISH were fixed in 4% paraformaldehyde (PFA) for 3 hours and then stored in 1:1 ethanol:PBS, while swabs for cell counting were stored in 4% PFA until analysis to avoid any cell loss during centrifugation that could impact quantification of very low biomass samples. Cell counting was performed with SYBR Gold as in [3], but adapted for very low biomass. FISH was performed as in [4], using DAPI as a counterstain for all DNA-containing cells. Negative controls included wells with no probe applied, to test for natural fluorescence, as well as the nonsense probe NONEUB (reverse complement of EUB338) to evaluate non-specific probe binding. Results and Discussion: Cleanroom surfaces had 102-103 cells cm-2. The extremely low biomass of these samples was challenging for enumeration, and required careful and routine use of “field” and laboratory blanks. FISH was performed with the general archaeal and bacterial probes ARCH915 and EUB338 (EUBMIX, [4]), probe GAMBET ([4]), and PSE227, which targets the genus Pseudomonas [5]). The latter two probes were selected because Pseudomonas spp. and other Gammaproteobacteria have not been isolated from cleanroom surfaces but do appear frequently in rRNA gene libraries from these surfaces. While some active bacteria were identified (Fig. 1c), most cells detectable by DAPI did not have a strong or any fluorescent signal (e.g., Fig. 1d), indicating that the vast majority of cells are dead or inactive. This suggests that cleaning protocols are effective at inactivating microbial contaminants, but that dead or inactive cells can remain on surfaces. Cells were often clumped in a weakly autofluorescent matrix, possibly biofilm material (Fig. 1c,d). We also observed other particulate material that was collected by the swabs, including apparent textile fibers (Fig. 1b). Our results are consistent with other studies that show that the bioburden present in clean rooms includes active, dormant, and dead cells. We will discuss how FISH and epifluorescent cell counting could be applied in planetary protection protocols, including the advantages and disadvantages of FISH and cell counting for routine use, as well as different possible applications for more specialized FISH procedures. References: [1] Moissl-Eichinger et al. (2015) Sci Rep, 5, 9156 [2] Hendrickson et al. (2021) Microbiome, 9, 238 [3] Jones et al. (2017) Appl Environ Microbiol, 83, e00909-17 [4] Jones et al. (2015) Appl Environ Microbiol, 81, 1242-1250. [5] Watt et al. (2006) Environ Microbiol, 8, 871-884

C J Huff↗

NASA NTRS2025

Lunar Terrain Vehicle (LTV) Remote Teleoperation Studies Under Four Lunar Communication Latencies

Remotely operating a lunar rover from Earth while subject to an Earth-Moon time delay of multiple seconds could result in a dangerous state where the roving vehicle is either damaged or lost, thereby potentially compromising an entire mission or series of missions. Providing the right capabilities to the remote operator to manage inherent communication latencies will be important for remote driving to be successful. The National Aeronautics and Space Administration (NASA) conducted two studies to investigate the average speed and number of kilometers per day that an operator on Earth could teleoperate a notional Artemis unpressurized rover with minimal remote operator capabilities under 0- and 4-second communication delays (April 2023 study) and 6- and 8-second delays (August 2023 study). A primary goal of these studies was to understand if an Artemis Lunar Terrain Vehicle (LTV) could cover 6 kilometers (km) in 24 hours when operated remotely. During the April 2023 evaluation, eight test operators used an in-house simulation of the lunar surface South Pole to teleoperate a NASA government reference LTV. Each operator received approximately 30 minutes of remote driving familiarization/training prior to their test run. Operators viewed the surrounding terrain via a single, rover mast-mounted, high-resolution camera with pan/tilt/zoom capabilities; continuous communication was provided throughout all testing. In the August 2023 evaluation, remote operators received approximately 3 hours of familiarization training in each latency, and the simulation environment provided remote operators with an operator-selected rate limiter to enable finer sensitivity in the hand controller and a predictive circle function to better assist operators with predicting the path the vehicle could take. All test operators were able to successfully navigate and drive through six different types of terrain and five planned traverse scenarios using natural lighting under all communication delays. Results for average speeds for each communication delay, computed by averaging the data from all test conditions for that latency and all operators, are shown in the table below. The average speed data was then used to derive the total time needed to cover 6 km, 8 km, and 20 km. Remote operators drove slower and used the brake more frequently when subject to a communication latency as opposed to no communication latency. Subjective workload assessments revealed that while operating in a latency the overall workload significantly increased when compared to a 0-s delay with mental demand, frustration, and performance being the primary contributing factors. Driving strategies in the 0-s delay did not vary significantly among subjects; however, in the 4-s delay condition, three different driving strategies were identified. In the 6-s and 8-s latency conditions the operator’s use of the cruise control to maintain speed was more apparent. Additionally, over the course of the August study, the operator took advantage of the predictive circle indicator on the navigation display and over 95% of the operator’s navigation used the mast camera 180-degree panning function for ground truthing in terms of boulders and craters. Operators started to define more specific parameters in driving strategies for general operations. This consisted of setting the vehicle into a low-speed cruise mode of approximately 1–1.5 kph and noticing driving performance of the vehicle seemed to be much harder at slower speeds 0.4–0.8 kph; however, the vehicle was more responsive at speeds of 2.9–3.6 kph. Regardless of communication delay, operators used both the horizontal translation rails and the vehicle fenders as guides to predict a path for the vehicle through heavily concentrated terrain features. Test operators acknowledged that the teleoperations training for this study was substantially less than what an actual LTV remote operator will ultimately receive. They estimated a minimum of 20 to 100 hours spread across multiple days and weeks (e.g., strategies included immersion training over a 3-day period, to a short 8-week starter program) would be needed to get an operator ~ 60% proficient (i.e., able to complete a subset of remote driving tasks), to a yearlong program for full proficiency in remote driving tasks under all terrain types and natural lighting conditions. Remotely operating a vehicle on another planetary body while subject to communication latency is a complex task. Speed, distance covered, time spent driving, time spent navigating, brake usage and rock contacts are all affected by operator workload, driving strategies, workstation ergonomics and training. These studies provided a “firstlook” answer to a potential system requirement (namely if a remote operator could cover a given distance in a given amount of time); however, considerable general knowledge was gained to begin to understand what it will take to make a successful lunar rover teleoperator.

LTV↗

NASA NTRS1999

Experimental Investigation of Pool Boiling Heat Transfer Enhancement in Microgravity in the Presence of Electric Fields

In boiling high heat fluxes are possible driven by relatively small temperature differences, which make its use increasingly attractive in aerospace applications. The objective of the research is to develop ways to overcome specific problems associated with boiling in the low gravity environment by substituting the buoyancy force with the electric force to enhance bubble removal from the heated surface. Previous studies indicate that in terrestrial applications nucleate boiling heat transfer can be increased by a factor of 50, as compared to values obtained for the same system without electric fields. The goal of our research is to experimentally explore the mechanisms responsible for EHD heat transfer enhancement in boiling in low gravity conditions, by visualizing the temperature distributions in the vicinity of the heated surface and around the bubble during boiling using real-time holographic interferometry (HI) combined with high-speed cinematography. In the first phase of the project the influence of the electric field on a single bubble is investigated. Pool boiling is simulated by injecting a single bubble through a nozzle into the subcooled liquid or into the thermal boundary layer developed along the flat heater surface. Since the exact location of bubble formation is known, the optical equipment can be aligned and focused accurately, which is an essential requirement for precision measurements of bubble shape, size and deformation, as well as the visualization of temperature fields by HI. The size of the bubble and the frequency of bubble departure can be controlled by suitable selection of nozzle diameter and mass flow rate of vapor. In this approach effects due to the presence of the electric field can be separated from effects caused by the temperature gradients in the thermal boundary layer. The influence of the thermal boundary layer can be investigated after activating the heater at a later stage of the research. For the visualization experiments a test cell was developed. All four vertical walls of the test cell are transparent, and they allow transillumination with laser light for visualization experiments by HI. The bottom electrode is a copper cylinder, which is electrically grounded. The copper block is heated with a resistive heater and it is equipped with 6 thermocouples that provide reference temperatures for the measurements with HI. The top electrode is a mesh electrode. Bubbles are injected with a syringe into the test cell through the bottom electrode. The working fluids presently used in the interferometric visualization experiments, water and PF 5052, satisfy requirements regarding thermophysical, optical and electrical properties. A 30kV power supply equipped with a voltmeter allows to apply the electric field to the electrodes during the experiments. The magnitude of the applied voltage can be adjusted either manually or through the LabVIEW data acquisition and control system connected to a PC. Temperatures of the heated block are recorded using type-T thermocouples, whose output is read by a data acquisition system. Images of the bubbles are recorded with 35mm photographic and 16mm high-speed cameras, scanned and analyzed using various software packages. Visualized temperature fields HI allows the visualization of temperature fields in the vicinity of bubbles during boiling in the form of fringes. Typical visualized temperature distributions around the air bubbles injected into the thermal boundary layer in PF5052 are shown. The temperature of the heated surface is 35 C. The temperature difference for a pair of fringes is approximately 0.05 C. The heat flux applied to the bottom surface is moderate, and the fringe patterns are regular. In the image a bubble penetrating the thermal boundary layer is visible. Because of the axial symmetry of the problem, simplified reconstruction techniques can be applied to recover the temperature field. The thermal plume developing above the heated surface for more intensive heating is shown. The temperature distribution in the liquid is clearly 3D, and tomographic techniques have to be applied to recover the temperature distribution in such a physical situation. A sequence of interferometric images showing the temperature distribution around the rising bubble, recorded with a high-speed camera is shown. Again, the temperature distribution is 3D, and a more complex approach to the evaluation, the tomographic reconstruction has to be taken. Measurement of the temperature distribution from the fringe pattern temperature distributions that yield important information regarding heat transfer are determined. Two algorithms that allow the quantitative evaluation of interferometric fringe patterns and the reconstruction of temperature fields during boiling have been developed at the Heat Transfer Laboratory of the Johns Hopkins University. In the first algorithm the bubble is assumed to be axially symmetrical, which significantly reduces the computational effort for quantifying the temperature distribution around the bubble. For this purpose the thermal boundary layer around the bubble is divided into equidistant concentric shells, and the refractive index is assumed to be constant in each of the shells. Since large temperature gradients are expected in the vicinity of the bubble during boiling, the deflection of the light beam cannot be neglected in boiling experiments. Since the exit angle of the light beam is known, this allows to account for the deflections and phase shifts outside the boundary layer (in the bulk fluid and in the windows of the test cell). Three dimensional temperature distributions in the vicinity of the bubble are reconstructed using tomographic techniques. In tomography, the measurement volume is sliced into 2D planes. In the present study these planes are parallel to the heated surface. The objective is to determine the values of the field parameter of interest in form of the field function in these 2D planes. The field parameter is the change of the refractive index of the liquid in the measurement volume caused by temperature changes. By superimposing data for many 2D planes recorded at the same time instant, the 3D temperature distribution in the measurement volume is recovered.

Herman, Cila↗

DOE OSTI2020

Dynamic projection of anthropogenic emissions in China: methodology and 2015–2050 emission pathways under a range of socio-economic, climate policy, and pollution control scenarios

Abstract. Future trends in air pollution and greenhouse gas (GHG)emissions for China are of great concern to the community. A set of globalscenarios regarding future socio-economic and climate developments, combiningshared socio-economic pathways (SSPs) with climate forcing outcomes asdescribed by the Representative Concentration Pathways (RCPs), was createdby the Intergovernmental Panel on Climate Change (IPCC). Chinese researchers have also developed various emission scenarios by considering detailed local environmental and climate policies. However, a comprehensive scenario set connecting SSP–RCP scenarios with local policies and representing dynamic emission changes under local policies is still missing. In this work, to fill this gap, we developed a dynamic projection model, the Dynamic Projection model for Emissions in China (DPEC), to explore China'sfuture anthropogenic emission pathways. The DPEC is designed tointegrate the energy system model, emission inventory model, dynamicprojection model, and parameterized scheme of Chinese policies. The modelcontains two main modules, an energy-model-driven activity rate projectionmodule and a sector-based emission projection module. The activity rateprojection module provides the standardized and unified future energyscenarios after reorganizing and refining the outputs from the energy systemmodel. Here we use a new China-focused version of the Global ChangeAssessment Model (GCAM-China) to project future energy demand and supply inChina under different SSP–RCP scenarios at the provincial level. Theemission projection module links a bottom-up emission inventory model, theMulti-resolution Emission Inventory for China (MEIC), to GCAM-China andaccurately tracks the evolution of future combustion and production technologiesand control measures under different environmental policies. We developedtechnology-based turnover models for several key emitting sectors (e.g.coal-fired power plants, key industries, and on-road transportationsectors), which can simulate the dynamic changes in the unit/vehicle fleetturnover process by tracking the lifespan of each unit/vehicle on an annualbasis. With the integrated modelling framework, we connected five SSP scenarios(SSP1–5), five RCP scenarios (RCP8.5, 7.0, 6.0, 4.5, and 2.6), and threepollution control scenarios (business as usual, BAU; enhanced controlpolicy, ECP; and best health effect, BHE) to produce six combined emissionscenarios. With those scenarios, we presented a wide range of China's futureemissions to 2050 under different development and policy pathways. We foundthat, with a combination of strong low-carbon policy and air pollutioncontrol policy (i.e. SSP1-26-BHE scenario), emissions of major airpollutants (i.e. SO 2 , NO x , PM 2.5 , and non-methane volatile organic compounds – NMVOCs) in China willbe reduced by 34%–66% in 2030 and 58%–87% in 2050 compared to 2015. End-of-pipe control measures are more effective for reducing air pollutant emissions before 2030, while low-carbon policy will play a more important rolein continuous emission reduction until 2050. In contrast, China's emissionswill remain at a high level until 2050 under a reference scenario without activeactions (i.e. SSP3-70-BAU). Compared to similar scenarios set from theCMIP6 (Coupled Model Intercomparison Project Phase 6), our estimates ofemission ranges are much lower than the estimates from the harmonized CMIP6 emissions dataset in2020–2030, but their emission ranges become similar in the year 2050.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2024

Lunar Terrain Vehicle (LTV) Remote Teleoperation Studies Under Four Lunar Communication Latencies

Remotely operating a lunar rover from Earth while subject to an Earth-Moon time delay of multiple seconds could result in a dangerous state where the roving vehicle is either damaged or lost, thereby potentially compromising an entire mission or series of missions. Providing the right capabilities to the remote operator to manage inherent communication latencies will be important for remote driving to be successful. NASA conducted two studies to investigate the average speed and number of kilometers per day that an operator on Earth could teleoperate a notional Artemis unpressurized rover with minimal remote operator capabilities under 0- and 4-second communication delays (April 2023 study) and 6- and 8-second delays (August 2023 study). A primary goal of these studies was to understand if an Artemis Lunar Terrain Vehicle (LTV) could cover 6 kilometers (km) in 24 hours when operated remotely. During the April 2023 evaluation, eight test operators used an in-house simulation of the lunar surface South Pole to teleoperate a NASA government reference LTV. Each operator received approximately 30 minutes of remote driving familiarization/training prior to their test run. Operators viewed the surrounding terrain via a single, rover mast-mounted, high-resolution camera with pan/tilt/zoom capabilities; continuous communication was provided throughout all testing. In the August 2023 evaluation, remote operators received approximately 3 hours of familiarization training in each latency, and the simulation environment provided remote operators with an operator-selected rate limiter to enable finer sensitivity in the hand controller and a predictive circle function to better assist operators with predicting the path the vehicle could take. All test operators were able to successfully navigate and drive through six different types of terrain and five planned traverse scenarios using natural lighting under all communication delays. Results for average speeds for each communication delay, computed by averaging the data from all test conditions for that latency and all operators, are shown in the table below. The average speed data was then used to derive the total time needed to cover 6 km, 8 km, and 20 km (distances relevant to LTV-SYS071 and -029 requirements). Remote operators drove slower and used the brake more frequently when subject to a communication latency as opposed to no communication latency. Subjective workload assessments revealed that while operating in a latency the overall workload significantly increased when compared to a 0-s delay with mental demand, frustration, and performance being the primary contributing factors. Driving strategies in the 0-s delay did not vary significantly among subjects; however, in the 4-s delay condition, three different driving strategies were identified. In the 6-s and 8-s latency conditions the operator’s use of the cruise control to maintain speed was more apparent. Additionally, over the course of the August study, the operator took advantage of the predictive circle indicator on the navigation display and over 95% of the operator’s navigation used the mast camera 180-degree panning function for ground truthing in terms of boulders and craters. Operators started to define more specific parameters in driving strategies for general operations. This consisted of setting the vehicle into a low-speed cruise mode of approximately 11.5 kph and noticing driving performance of the vehicle seemed to be much harder at slower speeds 0.4–0.8 kph; however, the vehicle was more responsive at speeds of 2.9–3.6 kph. Regardless of communication delay, operators used both the horizontal translation rails and the vehicle fenders as guides to predict a path for the vehicle through heavily concentrated terrain features. Test operators acknowledged that the teleoperations training for this study was substantially less than what an actual LTV remote operator will ultimately receive. They estimated a minimum of 20 to 100 hours spread across multiple days and weeks (e.g., strategies included immersion training over a 3-day period, to a short 8-week starter program) would be needed to get an operator ~ 60% proficient (i.e., able to complete a subset of remote driving tasks), to a yearlong program for full proficiency in remote driving tasks under all terrain types and natural lighting conditions. Remotely operating a vehicle on another planetary body while subject to communication latency is a complex task. Speed, distance covered, time spent driving, time spent navigating, brake usage and rock contacts are all affected by operator workload, driving strategies, workstation ergonomics and training. These studies provided a “first-look” answer to a potential system requirement (namely if a remote operator could cover a given distance in a given amount of time); however, considerable general knowledge was gained to begin to understand what it will take to make a successful lunar rover teleoperator.

LTV↗

NASA NTRSDate not supplied

Does Collection Time Bias the Ecology of Cleanroom Air Samples?

Microbial monitoring of astromaterials collections has taken on increased importance with the return of biologically sensitive samples from the asteroids Ryugu and Bennu and the initiation of the Mars Sample Return Program. Terrestrial bacteria and fungi can alter the mineralogy and organic composition of our collections causing irreversible contamination of pristine samples and increasing the risk of false positives for life detection measurements. NASA has conducted routine microbial monitoring of its existing collections since 20181. Initial monitoring focused on surface samples collected with foam swabs. Although, airborne microbiology is often decoupled from surface microbiology in the built environment2 culture-based air sampling techniques like impactors were not compliant with existing contamination control requirements. Bringing organic rich media, gelatin or liquids into curation cleanrooms presents an unacceptable risk to pristine samples. In 2022 NASA purchased a materials complaint air sampler and began collecting air samples from the cleanrooms in addition to surface samples3. The new instrument uses an electret filter to collect samples that are suitable for cultivating organisms or for direct DNA sequencing. Preliminary DNA sequencing results appeared to indicate that longer sampling times biased the microbial community in favor of hearty, spore-forming bacteria3. We present the results of a study comparing overnight sampling (17 hours) to short (1 hour) sampling of unoccupied curation cleanrooms. The results will help us optimize our monitoring protocols and develop a more detailed inventory of the ecology of astromaterials curation cleanrooms. Methods: We analyzed 72 paired air samples from six different cleanrooms including the meteorite processing lab (ISO 7 equivalent, 16 samples), the lunar lab (ISO 6 equivalent, 10 samples), the stardust lab (ISO 5 equivalent 14 samples), the OSIRIS-REx lab (ISO 5 equivalent, 12 samples), the Hayabusa2 lab (ISO 5 equivalent, 14 samples), and the Genesis lab (ISO 4 equivalent, 6 samples). All the samples were collected with an InnovaPrep Bobcat air sampler operating at a sampling rate of 200 L/min. The sampler operates for 5 minutes out of every 20 minute period. Half of the samples were collected by filtering 3,000L (15 min. of active sampling) of air across an electret filter for one hour. The rest of the samples were collected by filtering approximately 51,000 L air across the filter overnight (~17 hours, 255 min. of active sampling). Cells were eluted from the filter using 6-7 ml of pressurized 0.15% tween 20 in PBS (phosphate buffered saline). This liquid was used to cultivate bacteria according to previously published methods1,4,5 and for DNA extraction and next generation sequencing. DNA was extracted with a Qiagen MagAttract PowerMicrobiome kit6. To identify bacteria and archaea, the 16S rRNA gene was amplified using Earth Microbiome primers for the V4 region 7. The amplified DNA was sequenced on an Illumina MiSeq using a V3 reagent kit. The resulting sequences were processed using DADA2 and QIIME2 as implemented on the EDGE bioinformatics platform8–10. Results: Only two of the 72 samples had no amplifiable DNA. Amplified DNA concentrations ranged from 2.67 – 0.272 ng/µl. The median concentration of amplified DNA for the 1 hour samples was 0.770 ± 0.368 ng/µl. The median concentration of amplified DNA for the overnight samples was 0.877 ± 0.434 ng/µl. On average the overnight samples had slightly more sequences (58,960 vs. 59,456) and ASV’s (amplicon sequence variants) (60 vs 64.5) than the one hour samples, but these differences are not statistically significant. The most abundant ASV in every sample mapped to the genus Cupravidus. ASV’s mapping to the genuses Bacillus, Schlegelella, Thermus, and Staphylococcus were also common. Discussion and Future Work: Alpha diversity statistics like Shannon Entropy and Faith Phylogenetic Diversity are used to describe the diversity of organisms in a single sample. If a longer sampling time was biasing the data, we would expect to see a change in these diversity statistics vs. sample time. However, we did not observe this in our data. The median Shannon entropy was slightly higher for the overnight samples (3.773 vs 3.611) as was the Faith Phylogenetic Diversity (4.042 vs 3.596), but both values were within a standard deviation of each other for the two sampling times (Fig. 1). It is unlikely, that the longer sampling time is introducing bias into our data. We do observe a significant decrease in diversity when comparing the air samples by lab. The Genesis lab (ISO 4 equivalent) has a lower median number of ASV’s (45.5) than the other labs (62). Median values for Shannon Entropy (3.717 vs. 3.430) and Faith Phylogenetic Diversity (3.796 vs. 3.548) are also lower for Genesis, but those values are with one standard deviation of each other for the different sampling times. This is consistent with previous culture-based results suggesting that the environment in cleanrooms tends to select for a core group of organisms capable of surviving under dry, low nutrient, conditions. The presence of the ASV’s mapping to Cupravidus and Thermus in our sequencing blanks and controls suggests that several of the most common organisms in our samples represent contaminants from the reagents used to perform the DNA extractions and sequencing. Further work is needed to identify these contaminants, remove them from our data and recalculate the diversity statistics. This is a systematic error. Therefore, we do not expect removing the sequencing contaminants to change our conclusions. Longer air sample collection times appear to result in slightly higher diversity and do not bias the results towards “hardy” bacteria like spore-formers. Based on these preliminary results we conclude that sampling at least 3,000 liters of air is sufficient to capture the microbial diversity of cleanrooms, and that air samples can also be collected overnight without negatively impacting diversity. These results allow us to be flexible when designing microbial monitoring plans so that they do not interfere with routine lab activity. References: 1. Regberg, A. B. et al. 49th Lunar and Planetary Science Conference (2018). 2. The United States Pharmacopeial Convention. USP General Chapter <1116> (2013). 3. Regberg, A. B., et al. 54th Lunar and Planetary Science Conference (2023). 4. Regberg, A. B. et al. 53rd Lunar and Planetary Science Conference ( 2022). 5. Davis, R. E.,et al. 50th Lunar and Planetary Science Conference (2019). 6. Qiagen. MagAttract® PowerMicrobiome® DNA/RNA EP Kit Handbook. (2018). 7. Walters, W. et al. mSystems 1, (2015). 8. Callahan, B. J. et al. Nat. Methods 13, 581–583 (2016). 9. Hall, M. & Beiko, R. G. Microbiome Analysis: Methods and Protocols113–129 (Springer, 2018). 10. Philipson, C. et al. Bio-Protoc. 7, e2622 (2017).

A. B. Regberg↗

DOE OSTI2023

Direct observation of C 3 S particle dissolution using fast nano X-ray computed tomography

Tricalcium silicate (C 3 S) occupies 50 % to 70 % of ordinary portland cement (OPC) by mass and it is an important component affecting the hydration of OPC [1], [2], [3], [4], [5], [6], [7]. Generally, the hydration of C 3 S is described by two processes: the dissolution of C 3 S particles and the precipitation of hydration products. While it is understood that the dissolution rates of C 3 S vary with time, more precise measurements are needed to understand this process. Many mechanisms have been proposed to explain the time-evolving dissolution rates of C 3 S [2]. The metastable barrier hypothesis suggests that a thin metastable layer of hydrates forms around the C 3 S particle surface and prohibits the access of grains to the aqueous solution [8], [9], [10], [11], [12]. The slow dissolution step hypothesis suggests that the increased ion concentration from the initial reaction delays the C 3 S dissolution [2], [13], [14], [15], [16], [17]. More recent publications suggest that C 3 S may react differently depending on the existence of crystallographic defects [18], [19], [20]. Etch pits are thought to open on the particle surface during the initial reaction and this contributes to the C 3 S dissolution [21], [22]. As hydrates precipitate and cover these highly reactive surfaces, hydration slows down and the induction period starts [18], [23], [24], [25], [26]. Many experiments have been conducted to test the aforementioned mechanisms. Some hydration studies utilize bulk measurements, such as isothermal calorimetry [27], [28], [29], pore solution analysis [30], quasi-elastic neutron scattering [31], and nuclear magnetic resonance spectroscopy [32], [33]. One limitation of these measurements is that they do not provide direct and detailed information on the individual C 3 S particles. Some other studies utilize imaging techniques, such as scanning electron microscopy (SEM) [34], [35], [36] and transmission electron microscopy (TEM) [37]. However, SEM/TEM cannot track the evolution of individual particles throughout hydration [34], [35], [38], [39], [40] and they do not give insights into the microstructure of materials before hydration [34], [35], [39]. This makes it challenging to draw strong conclusions from only SEM or TEM observations. Synchrotron X-ray tomography techniques have been used more broadly in recent years to study cement hydration. They are not only non-destructive but also able to image a sample in full 3D with resolutions that can reach from micron to nanoscale. Nano computed tomography (nCT) is one technique that has been applied to study cement hydration at the nanoscale [26], [41]. A typical nCT can reach a pixel size from 15 to 65 nm, providing enough detail for observing features <1 μm. However, nCT often takes >0.5 h to finish one scan. This makes the application of this technique on continuous scans for in-situ observations challenging. Fast X-ray computerized tomography (fCT) is another technique that has shown success in studying the time-evolving cement microstructures [20], [42], [43], [44], [45], [46], [47], [48], [49]. Due to the high flux of the X-ray beam from the synchrotron ring, fCT allows a scan to be captured within 1 min at a pixel size of 1 μm. This allows a paste sample to be continuously scanned during the hydration process. However, the micron-sized resolutions limit does not provide detailed insights for particles <5 μm [20], [49]. Fortunately, the combination of nCT and fCT has allowed the development of fast nano X-ray computed tomography (fnCT). fnCT can capture a 3D data set in <2 min at a pixel size of 50 nm. This makes this procedure an exciting method to evaluate hydrating pastes. fnCT collects multiple X-ray radiographs at various rotation angles and generates a 3D model of the scanned sample, which is also referred to as a 3D tomography [50], [51]. In one tomography, the X-ray absorptions of different components (e.g., C 3 S and hydrates) differ as functions of density and chemistry [52], [53]. These X-ray absorption contrasts can be used to extract detailed information about the 3D microstructure [26], [54], [55]. In this paper, fnCT is used to collect time-lapse tomographs of hydrating C 3 S paste from 18 min after mixing to 7 h of hydration. The bulk measurements of anhydrous C 3 S, as well as the microstructural changes of individual C 3 S particles, are directly observed, quantified, and discussed. The dissolution behavior of C 3 S particles at various size scales is systematically analyzed and compared. This work aims to find the relationship between the size of C 3 S particle sizes and their dissolution rates. This provides significant insights into the early-age hydration of C 3 S on length and time scales not previously possible. Because of the magnitude of the data and the substantial amount of observations, this work will solely focus on the change in the anhydrous particles. Changes in the hydration products will be reported in future work.

42 ENGINEERING↗

DOE OSTI2024

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗

DOE OSTI2022

Vitrification of Hanford Tank 241-AP-101 Waste and Simulant

Hanford tank 241-AP-101 (referred to herein as AP-101) is the second Hanford radioactive tank waste planned to be processed and vitrified. A simulant version of AP 101 waste was formulated from the best-basis inventory (BBI) for the Hanford Tank 241-AP-101 liquid with an assumed target dilution of the waste from the BBI sodium molarity of 8.61 M to the desired 5.5 M Na. After the addition of glass-forming chemicals (GFCs), the simulant melter feed was processed in a non-radioactive, continuous laboratory-scale melter (CLSM) system. The AP-101 simulant melter feed was charged into the CLSM for 6.11 h of processing, which produced 6.55 kg of glass, for an average glass production rate of 2275 kg m 2 d -1 . Since there were no processing issues with the AP-101 simulant melter feed, AP-101 melter feed made with actual waste was then processed in a CLSM system built into a contamination area in a radioactive environment. The melting behavior characteristics appeared similar for both the simulant and waste melter feeds. The AP-101 waste melter feed was charged into the CLSM for 12.14 h of processing, which produced 8.75 kg of glass, for an average glass production rate of 1530 kg m 2 d -1 . During the AP-101 waste melter feed charging, the pump used to move the feed reached a maximum and it is believed that if the pump had a greater capacity, a greater average glass production rate could have been achieved. A constituent of interest present in low quantities in the AP-101 waste is 99 Tc or its non-radioactive surrogate, Re, added to the AP-101 simulant. Analysis for the quantities of 99 Tc and Re in the AP-101 glass product resulted in an average single-pass retention from the melter feed during relative chemical steady state of 55 ± 2 % for 99 Tc and 45 ± 2 % for Re. Compared to the processing of other melter feeds, the retention of 99 Tc in the AP-101 glass was greater than in both AP-107 and AP-105 glass, while the retention of Re in the AP-101 was less than in the AP-107 glass, but greater than in the AP-105 glass. A spike of I was added into the AP-101 melter feed that could be detected above the analysis detection limits. However, the iodine was only detectable above the ~6 ppm limit in one glass pour: the pour immediately following the burn off of the cold cap, where the I level reached ~30 ppm. This event was significant because the glass was poured immediately after burn off and thus it is presumed that the iodine had yet to volatilize from the glass melt while idling. It is recommended to perform future tests with I spikes at greater levels so that it can be detected in additional glass pours to determine if the expected 50 % retention of I used in the Kim et al. glass models can be confirmed. Offgas liquid samples were analyzed for acetonitrile, which was present at greater concentrations in CLSM liquids than in other scaled melter systems. This result was expected based on unique conditions with the CLSM system including a small plenum space leading to low residence time for offgas and the rapidity of offgas cooling upon exiting the CLSM vessel due to the location and environment. About 90 % of the total acetonitrile captured during both the AP-101 simulant and waste CLSM runs was found in the offgas condensate and demister liquids, thus it is recommended that only those liquids be sent for analysis if future testing to study the presence of acetonitrile in offgas products is desired.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRSDate not supplied

Trace Elemental Abundances in Calcium-Aluminum-Rich Inclusions in CV Chondrites

Introduction: Calcium-aluminum-rich inclusions (CAIs), are the first formed solids that define the age of the Solar System [1,2]. CAIs are thought to have condensed from nebular gas [3,4] within the first <1 Ma of Solar System formation [5,6]. CAIs have experienced numerous early Solar System processes including condensation, evaporation, melting, recrystallization, and aqueous alteration [e.g., 7]. The chemical, mineralogical, and textural diversity among CAIs results from a range of chemical and physical processes recorded during nebular and parent body epoch. This study aims to explore the mineralogical, textural, and chemical compositions of CAIs including the trace elemental abundances in CAI phases to determine the early Solar System processes recorded in them. Samples and Analytical Methods: We analyzed one CAI each from CV3 chondrites Northwest Africa (NWA) 5508 designated as ‘Saguaro’, and Northwest Africa (NWA) 12772 designated as ‘Hoopoe’. Back-scatter electron (BSE) images were collected using a Phenom XL scanning electron microscope (SEM) at the Lunar and Planetary Institute (LPI) and the JEOL JXA-8530F electron probe microanalyzer (EPMA) at Johnson Space Center (JSC)-NASA. Additionally, energy dispersive X-ray spectrometry (EDS) elemental maps of select areas for these samples were collected using a 15.0kV beam energy and a 40µA emission current. Using the EPMA, wavelength-dispersive X-ray spectroscopy (WDS) quantitative data were collected. In-situ trace element measurements for both CAIs were determined at JSC-NASA using a Photon Machines 193nm laser ablation system and a Thermo-Scientific Element-XR inductively coupled plasma mass spectrometer (ICP-MS). Analyses consisted of 30s ablations at 10Hz, spot sizes of 20-25µm, and a fluence of 6.0 J/cm2 for anorthite and melilite, and a 3.5 J/cm2 fluence for all other phases. NIST612 was used to correct for instrument drift, while BHVO-2g was used as a primary calibration standard. BCR-2g and in-house mineral standards were regularly measured as unknowns to ensure accuracy. Results: Saguaro is a coarse-grained CAI, ~11 x 6 mm in dimensions. Saguaro contains spinel, Al-rich pyroxene, anorthite, Mg-rich melilite, and minor perovskite in its interior and is therefore classified as a Type B CAI. Individual melilite grains shows normal compositional zoning with an Ak content ranging from ~24 to 54 with no apparent trend from the core to the edge of the CAI. The spinel appears euhedral and occurs both as clusters and as spinel palisades [8]. Two rim sequences surround most of the sample: the inner rim being a Wark-Lovering (WL) rim (~10-35 µm) containing pyroxene, spinel, and melilite (or anorthite), and the outer rim is a finer-grained, thicker (~100 µm), accretionary rim (Fig. 1). The mineral phases in Saguaro record an overall flat REE pattern with an average negative Eu anomaly in pyroxene, and an average positive Eu anomaly in anorthite and melilite respectively. Anorthite, melilite, and pyroxene have a minor depletion in Tm (Fig. 2). The Hoopoe CAI is a compact, coarse-grained ~6 × 4 mm in size. The major mineralogy includes hibonite, spinel, melilite, anorthite, and perovskite. Therefore, it is classified as a compact transitional type A and B. (?)zoning was observed in some hibonites. Individual melilite grains show both reverse and normal zoning, where the Ak content ranges from ~6- to 28. Melilite shows two distinct textures. One texture consisted of smooth melilite that appeared homogenous, while the second appeared to consist of many fine fractures. The spinel also often appears clustered. The WL-rim sequence surrounding Hoopoe is ~25 µm thick and composed of spinel, perovskite, hibonite, and melilite/anorthite. It is then partially surrounded by an outer accretionary rim (~75µm). Like before, refractory metal nuggets appeared concentrated near the WL rims. Other metal assemblages rich in Fe and Ni were also observed. All major mineral phases in Hoopoe display relatively flat REE patterns, except for varying Eu and Tm between phases (Fig. 2). There is a prominent negative Eu anomaly in perovskite and an average positive Eu anomaly in spinel, anorthite, and hibonite respectively (Fig. 2, 3). The mixed phases along the rim of the CAI also display a negative Eu anomaly, and all phases the CAI were depleted in Pb. Discussion: The CV3 CAIs analyzed in this study were classified based on their mineralogy and textures into Type A versus Type B CAIs [10]. Hibonite appears to be pseudomorphically replacing the spinel, (i.e., is hibonite in composition, but appears in the shape of spinel). Spinel palisades. The presence of spinel palisades present in Saguaro are consistent with the melting and recrystallization experienced by this CAI. Trace elemental analyses. Saguaro and Hoopoe display similar trace element patterns to each other, with both appearing generally flat, with anomalies in Eu, and Tm. Melilite and anorthite display positive Eu anomalies in both CAIs, in addition to the hibonite in Hoopoe (Fig. 2). The phases that are depleted in Eu are pyroxene and perovskite in both Saguaro and Hoopoe, respectively (Fig. 3). Given that Eu is volatile in reducing environments [11], this could possibly indicate reducing conditions at the time anorthite and melilite crystallized, with the gas they formed from containing Eu. As these CAIs continued to form, this gas as a result would become depleted in Eu. This also could be supported by the propensity of anorthite and melilite to take up Eu from its surroundings and incorporate it into their structure [12]. In addition, analyzing the assemblage of the phases in the Saguaro, melilite and anorthite (Eu enriched) often surround the pyroxene (Eu depleted) as they are crystallized. This intergrowth of phases and the proximity of the phases would support that the Eu is being incorporated into some phases, preventing it from incorporating into other. Trace elemental analyses of the CAI rims will be evaluated in more detail, as they are complicated by the transient signal being composed of a mixture of mineral phases. Broadly, however, the patterns in the rims of both CAIs are comparable to each other, and for Hoopoe, to the mixed phase patterns in the core (Fig. 3). Other studies have found that CAI rims can be depleted in Ce and Yb [13], however we did not observe these anomalies in the two CAIs discussed here. Given their similarity, the trace elemental analyses of the mixed interior (i.e. core) and rim phases could be interpreted as forming from similar, if not the same, reservoirs. The REE abundance between the rim and core of Hoopoe are also similar, indicating they may have formed from a gas of the same or similar composition. Acknowledgments: We thank the ASU Center for Meteorite Studies for loaning the samples used in this work and Tabb Prissel for his assistance with the analysis. Mouti Al-Hashimi thanks Sam Crossley and Cyrena Goodrich for their help with the LPI SEM training. This work was supported by the LPI Summer Intern Program in Planetary Science and the LPI Cooperative Agreement. References: [1] Connelly J.N. (2012) Science, 338, 651-655. [2] MacPherson G. J. (2014) Treatise on Geochem., 2, 139-179. [3] Grossman L. (1972) GCA, 36, 597-619. [4] Ebel, D.S. (2006) Meteorites and the Early Solar System II (D. S. Lauretta & H. Y. McSween, Eds.) 253-277. [5] MacPherson G. J. (2012) Earth Planet. Sci. Lett., 331-332, 43-54. [6] MacPherson G.J. (2017) GCA, 201, 65-82. [7] Krot A.N. (1995) Meteoritics & Planet. Sci., 30, 748-775. [8] Wark and Lovering (1982) GCA, 46, 2595-2607. [9] Palme H. and Jones A. (2003) Treatise on Geochemistry (H. D. Holland and K. K. Turekian Eds.), 1, 41-61. [10] Grossman L. (1980) Ann. Rev. Earth Planet. Sci., 8, 559-608. [11] Floss C. et al. (1996) GCA, 60, 1975-1997. [12] Mason B. and Martin P. M. (1974) Earth Planet. Sci. Lett., 22, 141-144. [13] Wark B. and Boynton W. V. (2001) Meteoritics & Planet. Sci., 36, 1135-1166.

X Mouti↗

DOE OSTI2022

Emissions mitigation technology for advanced water-lean solvent-based CO 2 capture processes

This technical final report submitted to DOE/NETL presents all the research activities performed during the entirety of DE-FE0031660 project-Emissions Mitigation Technology for Advanced Water-Lean Solvent-Based CO 2 Capture Processes which spans from October 2018 through March 2022. RTI International has been conducting studies from fundamental and operational aspects to reduce the overall amine emissions from the advanced Water-Lean Solvent (WLS) systems, specifically RTI’s Non-Aqueous Solvent (NAS). This technical final report will highlight the key findings from project which align closely to the project objectives which are: Identify the contribution of vapor loss, entrainment, and aerosols to the overall emissions of water-lean systems; Determine the significance of CO 2 capture system operating parameters to the amine emissions; Develop an emissions model based on critical operating parameters; Evaluate the effectiveness of emissions mitigation devices to reduce the amine emissions to <1 ppm under flue coal-fired flue gas; and, Determine the contribution of the ECTs to the overall CO 2 capture cost. The following are the key findings based on numerous tests using both lab-scale setups and parametric testing performed at RTI’s Bench-scale Gas Absorption System (BsGAS). During the BP1, the aerosol generation system and monitoring equipment were installed at BsGAS to produce and determine the aerosol characteristics during the NAS CO 2 capture process. The aerosol produced by this setup produced aerosols with the peak diameter and concentration of 50 micron and 1.2E10 7 cm -3 , respectively. These particle sizes and concentrations are matched to those observed in the actual coal-fired power plant flue gases and expected to be found at the absorber inlet of the CO 2 capture system. Over 1,300 hours of parametric testing have been conducted to evaluate the impact of the aerosols and operating conditions during the CO 2 capture with NAS on the overall amine emissions in the treated flue gas. At the worse condition tested, the presence of the aerosols in the flue gas could increase the overall emissions by 10X compared to the baseline emissions from NAS’s vapor pressure. CO 2 capture rate was found to be a main factor impacting the overall emissions as well as aerosol size and concentrations in the absorber off-gas. The higher CO 2 capture rate, the higher amine emissions in the treated gas. The temperature difference between the temperature bulge seen in the absorber and the water wash temperature also impacts the particle growth where the larger the temperature difference, the more amine emissions from aerosols in the treated gas. The majority of the aerosols did not grow substantially in the system, and the particle concentrations remained nearly constant between the absorber inlet and wash outlet. Only a small portion of the particles were found to grow significantly. The high efficiency demister with mesh size of 5-10 micron can be installed to remove a portion of the aerosols from the gas stream leaving the water wash. Overall, these results from parametric testing have established the emission baseline and validate our assumption on the need of emission control technologies (ECT) in order to minimize the emissions from the baseline NAS CO 2 capture process. Over 2,000 of BsGAS operating hours was used to investigate a handful of process improvements which led to a selection of the vital few changes that effectively control the amine emissions. These process improvements are lime-coated-filters for absorber gas inlet, advanced demister at the top of the absorber, a second water wash with amine recovery unit were designed, installed, and tested at BsGAS at the end of BP1. The result showed that the NAS CO 2 capture process with these additional emission control devices could lower the amine emission in the treated gas to about 1 ppm using a simulated coal-fire flue gas stream. The main contributor in lowering the amine emission came from the second water wash with amine recovery unit where the amine concentration in the scrubbing water was kept below 2 wt% through a continuous amine removal via an adsorbent bed, resulting in a low amine vapor pressure. The adsorbent bed was regenerated via a direct steam regeneration and the recover amine was returned to the absorber to minimize wastewater and makeup amine. A flue gas generation system was designed and installed during the first half of BP2 to support the emission testing using a real coal-derived flue gas. The system is capable of generating both coal- and natural gas- derived flue gases with the composition of the gaseous species highly resemble to that of the power plant flue gases. The particulates detected in the coal-derived flue gas showed the mean diameter of 1 micron. The CO 2 capture operating was then proceed using the real coal-derived flue gas where the amine emission was controlled to be about 0-3 ppm for the total run time of about 200 hours. Similar testing was conducted with natural gas-derived flue gas and the result showed a highly amine emission of 30 ppm under the total run time of 200 hours. The Principal Component Analysis (PCA) and the Partial Least Squares Projection to Latent Structures (PLS) techniques were applied to the parametric testing data to derive a multivariate statistical model. The model was validated and trained with half of the data collected, and the predictive ability of the model was evaluated using the remaining half of the data. The resulting empirical model was capable of predicting the overall emissions from the NAS process without the ECTs with ±15% accuracy (average absolute deviation, AAD) in BP1. As more emission data were obtained under the real coal-flue gas in the BP2, the model incorporated these new set of data to reflect the final process configuration, operating parameters, and amine emission. This results in the updated empirical model predicting the amine emission from the NAS CO 2 capture process with 84% goodness-of-fit (R 2 ), 85% predictability (Q 2 ), and 15% AAD. The study evaluates the use of RTI’s Non-Aqueous Solvent technology for 90% CO 2 capture from a net 650 MWe pulverized coal power plant, downstream of the flue-gas desulfurization unit. The captured CO 2 has a purity of > 95% CO 2 , and is dried, compressed to 15.3 MPa (2,215 psia), ready for sequestration. The analysis uses Case B12B from the DOE Baseline study on Bituminous Coal, Revision 4 where the Cansolv CO 2 capture plant is replaced by the RTI CO 2 Capture plant. The CO 2 capture plant has been sized to capture >90% CO 2 from flue gas derived from a net 650 MWe supercritical pulverized coal power plant. The CO 2 capture plant is equipped with emission control technologies that limits the amine emissions to < 1 ppm. Two different cases were evaluated for the technoeconomic study. The key difference between the two cases is the regenerator pressure. In Case 1, the regenerator operates at 0.195 MPa (28.3 psia), whereas in Case 2, the regenerator pressure is 0.44 MPa (64 psia) thus removing the need for the first stage of compression of the eight-stage compression train. Results from the TEA are compared against the DOE reference cases for SCPC plant with and without CO 2 Capture (Case B12A and Case B12B of the DOE Baseline study, respectively). Case 2 with CO 2 regeneration at higher pressure results in the lower cost of CO 2 capture. The total capital cost of the capture process has been estimated using 2018 dollars in Aspen Process Economic Analyzer and was estimated to be $579 MM. The capture plant operation leads to a total parasitic power loss rate of 96 MWe, resulting in a decrease in pulverized coal power plant efficiency of 7.8% points. The resulting cost of electric power increases from 64.4 mills/kWh, for no capture, to 97.5 mills/kWh, with 90% capture, an increase of 51% in the COE. The cost of capturing 90% CO 2 was estimated to be $38.2/tonne-CO 2 , and meets the DOE target of $40/t-CO 2 . Emission control technologies (ECT) investigated in this project includes a second water wash with use of activated carbon beds for removal of amine from the wash water prior to recirculation in the water wash. These ECT allow operation of the CO 2 capture plant with < 1 ppm amine emissions with the treated flue gas and contributes to $2.4/t-CO 2 captured. Amine emissions derived from thermal and oxidative degradations were investigated under this project along with the emissions derived from aerosols for the NAS system. The thermally degraded of the lean NAS showed less than 4% decreased of the original total amine content in the NAS at 150 °C while the result obtained at 120 °C showed no drop in total amine content, suggesting that thermal degradation of the NAS is minimal. These results also suggested that the thermally degraded species are not likely formed and contributed to the emissions due to the low regeneration temperature of the NAS at 90-105 °C. The oxidative degradation, on the other hand, could become problematic as some of these oxidative degraded species were observed during the NAS-5 testing at National Carbon Capture Center (NCCC) and SINTEF in our previous project. The rapid screening of selected inhibitors suggested that oxidative degradation of NAS can be suppressed using thiol containing compounds in amounts of at least 1 mol%. The detailed mechanistic degradation pathway was conceived for a specific amine used in NAS formulation during BP2. he reduction of the nitrosamines caused by the NO x present in the flue gas was also examined. The study suggested that the thermo-chemical treatment of the NAS solvent would be a more effective and economically viable compared to removing NO x at the DCC.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2022

Quasars with large proper motions: A selection from the LQAC-5 catalogue combined with Gaia EDR3: Focusing on astrometric and photometric properties

Thanks to the cross-identification between the Gaia EDR3 (E-Data Release 3) catalogue and the quasars’ compiled catalogue LQAC-5 (Large Quasar Astrometric Catalogue-5), accurate astrometric parameters as well as photometric measurements were available for a set of 416 113 quasars. After analysing the astrometric and photometric properties of these quasars, we concentrate our study on objects characterized by significant proper motion, which contradicts the postulate that they are fixed in the celestial sphere. We derived statistical properties of the positions, parallaxes, and proper motions of the quasars cross-matched with Gaia EDR3. We carried out a specific investigation of quasars with a proper motion larger than 10 mas yr -1 , and used colour-colour diagrams in mid-infrared to discriminate these quasars from stellar contamination. We discuss, in detail, a few cases where the abnormal proper motion could be caused by misidentifications. We note that 43% of LQAC-5 quasars were not detected by Gaia EDR3, mainly because of the cut-off magnitude at G = 21. The comparison between r(G) (Gaia) and r (SDSS) shows a normal distribution. The normalized distributions of parallaxes and proper motions are identical to that of the Gaia CRF2, as well as the distribution of the G BP -G RP colour index. We have identified 937 LQAC-5 objects (0.27% of the sample) with Gaia EDR3 proper motion which is significant at the 5σ level. We have clearly shown that below z = 0.35, the integrated magnitude G(r) obtained from SDSS spectra is smaller than G (Gaia), leading to the statement that the quasars are generally resolved under this threshold, which should lead to a bad determination of proper motion. Moreover, to discriminate quasars from stellar contamination, we have plotted them in a mid-infrared colour–colour diagram. Finally, we have made a final selection of 41 objects with μ > 10 mas yr -1 , satisfying drastic constraints, which makes them potential candidates for considering their proper motion as real. In this paper we carried out a complete analysis of the quasars of the LQAC-5 which were cross-identified with Gaia EDR3 ; this serves as an alternative and complementary study with respect to Gaia CRF2, involving a different population of quasars. We selected a set of 41 quasars with a proper motion exceeding 10 mas yr -1 , which can be considered as very high for objects which are a priori fixed in the celestial sphere.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRSDate not supplied

4.32 Billion Year Old Impact Melts at Apollo 14: Dating the Procellarum Basin?

Introduction: Recent work has improved our understanding of lunar crustal structure and basin evolution on the Moon, but its early impact history and implications for solar system dynamics remain poorly established [1,2]. In particular, there is a derth of absolute ages, especially for the period >4.0 Ga, that could be related to basin-forming events. Here we present U-Pb isotopic data for Zr-rich minerals found in impact-melt fragments from Apollo 14 soil sample 14163 that, when combined with previously published data, dates an impact event at 4324±15 Ma. The compositions and ages of these impact-melt fragments provide unique information about the timing of early impact events and the composition of the lunar crust. Results: Fourteen rocklets ranging in size from 1 to 3 mm were extracted from soil 14163. Most are impactmelt rocks with 10-30% clasts of mostly pyroxene and plagioclase (100 to 500 μm) in a crystalline matrix. Textures of the matrix vary from subophitic, formed by intergrowths of 10-20 μm plagioclase and pyroxene crystals, to poikilitic, with plagioclase and pyroxene reaching 20-50 μm size. Some fragments that contain no visible clasts have similar textures are also interpreted as impact-melt rocks. All fragments contain notably large proportions of euhedral to subhedral ilmenite, zircon, apatite/merrillite and less abundant zirconolite and baddeleyite in the melt matrices. These grains often form intricate intergrowths with each other and rock-forming minerals, indicating their crystallization from the melt (Fig. 1). However, some slightly larger (~50 μm) zircon and phosphate grains can be interpreted as relict clasts based on their granular textures and relationships with the surrounding phases (Fig. 1). This implies that zircon and phosphate minerals were present in the target rocks. Some of these grains were profoundly remelted during the impact, which resulted in an oversaturation of the melt in Zr and P and crystallization of new grains of Zrrich minerals and phosphates during melt solidification. The presence of zircon and phosphates indicates that the melts were significantly enriched in KREEP components and ties their provenance to the Procellarum- KREEP Terrane (PKT) [3]. Combined U-Pb zircon data obtained for all fragments form two main clusters on a concordia diagram (Fig. 2), between about 4.3 and 3.9 Ga. Phosphate data mostly concentrate on the younger end of this age range but some analyses are almost as old as the older zircon grains (Fig. 2). Figure 1: Zircon (Zr) and phosphate (Apt; Mer) grains in impact melt fragments from Apollo 14 soil 14161. A-zircon grains crystallized from impact melt; B-phosphate grain inherited from the target; C- granular zircon grains Our interpretation of these data is that all U-rich minerals experienced variable resetting of the U-Pb system, with phosphate, where closure temperature is significantly lower than that in zircon, affected more profoundly by Pb loss than zircon. Our best estimate of the time of formation of zircon and phosphate is based on statistically valid analysis of grains interpreted texturally as grown from the impact melt (Fig. 2). Ten of these analyses define an age of 4324±15 Ma (MSWD=3.0, probability of fit P=0.002). Our best estimate for the time of resetting is obtained by combining data from phosphates that are statistically indistinguishable from 3.9 Ga within the analytical uncertainties. This group is represented by 18 analyses of 15 phosphate grains from different fragments (Fig. 2) and defines an age of 3922±6 Ma (MSWD=1.2, P=0.23). Conversely, a minimum age of the target lithologies, remelted in the impact that produced the rocklets studied here (Fig. 2), can be determined from the five oldest analyses of zircon clasts at 4338±13 Ma (MSWD=1.5, P=0.2), which is indistinguishable from the age of the impact melt within the uncertainties. The obtained ages, combined with textural evidence, imply that the impact melt was formed at 4324±15 Ma, and that it occurred in a zircon-rich target with a minimum age of 4338±13 Ma. Further reworking occurred during a second impact event at 3922±6 Ma. Discussion: The younger age of 3922±6 Ma can be interpreted as the time of the Imbrium impact. The older age of 4324±15 Ma would then be the time of formation of the impact melt, which was then caught in the Imbrium ejecta either at the Apollo 14 landing site or within the target rocks of the Imbrium impact. The presence of abundant accessory phases such as zircons and phosphates is consistent with a substantial KREEP component in the analyzed particles. The current distribution of KREEP on the lunar surface appears to be strongly influenced by Imbrum ejecta [4] so interpretation of the 4.32 Ga age depends in part on assumed structure of the pre-impact crust and distribution of KREEP within the crust. If KREEP was present only in the deep crust at 4.32 Ga, then a basin-scale impact possibly analogous to Imbrium or larger may be necessary to excavate a KREEPy impact melt at this time. However, if KREEPy materials were present closer to the surface perhaps due to redistribution related to Mg-suite magmatism, then smaller impacts might be able to rework KREEPy compositions at shallower depths. The coherence of the data on the particles analysed here suggests a large volume of melt that has been preserved since 4.32 Ga, consistent with a large impact event. Relict zircons and mineral clasts suggests that the igneous crust in the vicinity of this impact was well developed by at least 4.34 Ga, similar to the model age of KREEP and older than the isochron ages of many Mg-suite cumulates [5]. A problem that confronts all lunar sample studies using the current collection is that the pre-Imbrium geology of the PKT (the source of these 4.32 Ga impact melts) is not well constrained. The South Pole-Aitken basin contains regions that are moderately enriched in Th, but its ejecta is Th-poor [6]; therefore these A14 fragments are probably not SPA ejecta. Alternatively, these fragments may represent formation of a hypothesized Procellarum basin [7] although the lack of a clearly defined basin ejecta signature is a potential problem with that interpretation. In any case, they provide a younger limit on the age of lunar differentiation and formation of KREEP within the lunar crust. Figure 2: U-Pb data for zircon and phosphate grains from impact melt fragments. A-all data; B-data used for age calculations. References: [1] Orgel C., Michael G., Fassett C. I., van der Bogert C. H., Riedel C., Kneissl T., and Hiesinger H. (2018) J. Geophys. Res. Planets 123, 748- 762. [2] Evans A. J., Andrews-Hanna J. C., Head J. W., Soderblom J. M., Solomon S. C., and Zuber M. T. (2018) J. Geophys. Planets. 123, 1596-1617. [3] Jolliff, B.L., Gillis, J.J., Haskin, L.A., Korotev, R.L. and Wieczorek, M.A. (2000) J. Geophys. Res: Planets 105, 4197-4216. [4] Haskin L. A. (1998) J. Geophys. Res. Planets 103, 1679-1689. [5] Borg L.E., Gaffney A.M., and Shearer C.K. (2015) MAPS 50, 715-732. [6] Moriarty, D.P., Watkins, R.N., Valencia, S.N., Kendall, J.D., Evans, A.J., Dygert, N. and Petro, N.E. (2021) J. Geophys. Res. Planets 126. [7] Zhu, M.H., Wünnemann, K., Potter, R.W., Kleine, T. and Morbidelli, A. (2019) J. Geophys. Res: Planets 124, 2117-2140

M D Norman↗