Search NASA⌕ Search

SEARCH · Search NASA

Results for “map series”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Floating zone crystal growth, structure, and properties of a cubic Li 5.5 La 3 Nb 1.5 Zr 0.5 O 12 garnet-type lithium-ion conductor

As a promising candidate for solid-state electrolytes in Li-ion batteries, the garnet-type Li-ion conductor series Li 5+x La 3 Nb 2-x Zr x O 12 (LLNZO) (0 ≤ x ≤ 2) exhibits high ionic conductivity at room temperature. However, no previous single-crystal growth or characterization has been reported for LLNZO compositions 0 ≤ x ≤ 1. To obtain a complete understanding of the trend in the structure–property relationship in this class of materials, we used the floating zone (FZ) method to grow a single crystal of Li 5.5 La 3 Nb 1.5 Zr 0.5 O 12 that was 4 mm in diameter and 10 mm in length. Using Laue neutron single-crystal diffraction, two distinct Li sites were observed: tetrahedral 24d and octahedral 96h sites. The maximum entropy method (MEM) based on neutron single-crystal diffraction data was used to map Li nuclear density and estimate that the bottleneck of Li transport exists between neighboring tetrahedral and octahedral sites, and that Li is delocalized between split octahedral sites. Room-temperature Li-ion conductivity in Li 5.5 La 3 Nb 1.5 Zr 0.5 O 12 measured with electrochemical impedance spectroscopy (EIS) was 1.37 × 10 -4 S cm -1 . The Li migration activation energy was estimated to be 0.50 eV from EIS and 0.47 eV from dielectric relaxation measurements. The Li-ion jump attempt rate was estimated to be 1.47 × 10 12 Hz while the time scale of successful migration is 10 -7 to 10 -6 s.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Laboratory time series moisture manipulative experiment from sediment across the contiguous US: time series aerobic respiration and geochemistry (v2)

This dataset supports a broader study examining the effects of wetting and drying on hyporheic zone respiration across the contiguous United States (CONUS). The dataset provides data generated from a laboratory moisture manipulation experiment. The contents include time series aerobic respiration and moisture; dissolved oxygen; sediment geochemistry data; and field metadata (including qualitative information on instream and river corridor characteristics). Samples were collected as part of the WHONDRS CONUS-Scale Model-Sample Study (CM). This study was designed following ICON (integrated, coordinated, open, and networked) principles to facilitate a model-experiment (ModEx) iteration approach, leveraging crowdsourced sampling across the CONUS. The data package associated with the CM study is available at https://data.ess-dive.lbl.gov/view/doi:10.15485/1923689. CM sampling began in April 2022 and ended in October 2023. This study uses subsamples from a subset of CM samples collected between June 2022 and June 2023. The original field samples were labeled as CM_###. Subsequent subsamples for this study were labeled as EC_###. The labels from the field samples and the EC subsamples can be mapped directly based on the digits following the prefix and underscore (i.e., EC_001 is a subsample from CM_001). See the critical details section below for more details on sample naming. This data package was originally published in August 2024. It was updated in February 2026 (v2; new and modified files). See the change history section in the readme for more details. For details on how to navigate this data package, see this infographic from the River Corridor SFA https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. This dataset is comprised of one folder of raw Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) data and one main data folder containing (1) file-level metadata; (2) data dictionary; (3) field metadata; (4) readme; (5) field protocol; and a (6) a subfolder with sediment sample data from the incubation experiment. The sample data subfolder contains (1) dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC); (2) total nitrogen (TN); (3) adenosine triphosphate (ATP); (4) percent carbon and nitrogen; (5) effect size; (6) iron (II); (7) gravimetric moisture; (8) respiration rates and raw dissolved oxygen values; (9) specific conductance; (10) pH; (11) temperature; (12) a summary containing median values of each data type for each treatment (wet and dry); (13) methods codes; (14) FTICR-MS methods; and (15) a subfolder of 9.4 Tesla FTICR-MS data. This folder contains three subfolders, one containing the sediment .xml data files, one containing the sediment CoreMS output files, the other containing instructions and scripts for processing the files in CoreMS (https://github.com/EMSL-Computing/CoreMS). All files are .csv, .pdf, .R, .ref, or .xml.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Snow Distribution Patterns Revisited: A Physics-Based and Machine Learning Hybrid Approach to Snow Distribution Mapping in the Sub-Arctic

Snowpack distribution in Arctic and alpine landscapes often occurs in repeating, year-to-year patterns due to local topographic, weather, and vegetation characteristics. Previous studies have suggested that with years of observational data, these snow distribution patterns can be statistically integrated into a snow process modeling workflow. Recent advances in snow hydrology and machine learning (ML) have increased our ability to predict snowpack distribution using in-situ observations, remote sensing data sets, and simple landscape characteristics that can be easily obtained for most environments. Here, we propose a hybrid approach to couple a ML snow distribution pattern (MLSDP) map with a physics-based, snow process model. We trained a random forest ML algorithm on tens of thousands of snow survey observations from a subarctic study area on the Seward Peninsula, Alaska, collected during peak snow water equivalent (SWE). We validated hybrid model outputs using in-situ snow depth and SWE observations, as well as a light detection and ranging data set and a distributed temperature profiling sensor data set. When the hybrid results were compared with the physics-based method, the hybrid method more accurately depicted the spatial patterns of the snowpack, areas of drifting snow, and years when no in-situ observations were used in the random forest ML training data set. The hybrid method also showed improvements in root mean squared error at 61% of locations where time-series estimations of snow depth were observed. These results can be applied to any physics-based model to improve the snow distribution patterning to reflect observed conditions in high latitude and high elevation cold region environments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Evolution of magnetic field induced ordering in the layered quantum Heisenberg triangular-lattice antiferromagnet Ba 3 CoSb 2 O 9

Quantum fluctuations in the effective spin-12 layered triangular-lattice quantum Heisenberg antiferromagnet Ba 3 CoSb 2 O 9 lift the classical degeneracy of the antiferromagnetic ground state in magnetic field, producing a series of novel spin structures for magnetic fields applied within the crystallographic ab plane, including a celebrated collinear “up-up-down” spin ordering with magnetization equal to 1/3 of the saturation magnetization over an extended field range. Theoretically unresolved, however, are the effects of interlayer antiferromagnetic coupling and transverse magnetic fields on the ground states of this system. Additional magnetic field induced phase transitions are theoretically expected and in some cases have been experimentally observed, but details regarding their number, location, and physical character appear inconsistent with the predictions of existing models. Conversely, an absence of experimental measurements as a function of magnetic-field orientation has left other key predictions of these models untested. To address these issues, in this study we have used specific heat, neutron diffraction, thermal conductivity, and magnetic torque measurements to map out the phase diagram as a function of magnetic field intensity and orientation relative to the crystallographic ab plane. For H||ab, we have discovered an additional magnetic field induced phase transition at low temperature and an unexpected tetracritical point in the high-field phase diagram, which coupled with the apparent second-order nature of the phase transitions eliminates several theoretically proposed spin structures for the high-field phases. Our calorimetric measurements as a function of magnetic field orientation are in general agreement with theory for field-orientation angles close to plane parallel (H||a) but diverge at angles near plane perpendicular; a predicted convergence of two phase boundaries at finite angle and a corresponding change in the order of the field induced phase transition are not observed experimentally. Our results emphasize the role of interlayer coupling in selecting and stabilizing field induced phases, provide guidance on the nature of the magnetic order in each phase, and reveal the need for new physics to account for the nature of magnetic ordering in this archetypal two-dimensional spin-12 triangular-lattice quantum Heisenberg antiferromagnet.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Validating the galaxy and quasar catalog-level blinding scheme for the DESI 2024 analysis

In the era of precision cosmology, ensuring the integrity of data analysis through blinding techniques is paramount — a challenge particularly relevant for the Dark Energy Spectroscopic Instrument (DESI). DESI represents a monumental effort to map the cosmic web, with the goal to measure the redshifts of tens of millions of galaxies and quasars. Given the data volume and the impact of the findings, the potential for confirmation bias poses a significant challenge. To address this, we implement and validate a comprehensive blind analysis strategy for DESI Data Release 1 (DR1), tailored to the specific observables DESI is most sensitive to: Baryonic Acoustic Oscillations (BAO), Redshift-Space Distortion (RSD) and primordial non-Gaussianities (PNG). We carry out the blinding at the catalog level, implementing shifts in the redshifts of the observed galaxies to blind for BAO and RSD signals and weights to blind for PNG through a scale-dependent bias. We validate the blinding technique on mocks as well as on data by applying a second blinding layer to perform a series of sanity checks; the latter allows probing complexities in real data not captured in mocks. We find that the blinding strategy alters the data vector in a controlled way, and the BAO and RSD analysis choices are robust to blinding. The successful validation of the blinding strategy paves the way for the unblinded DESI DR1 analysis, alongside future blind analyses with DESI and other surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Hydrologic Impacts of a Strike-Slip Fault Zone: Insights from Joint 3D Body-Wave Tomography of Rock Valley

The Rock Valley fault zone (RVFZ), an intraplate strike-slip fault zone in the southern Nevada National Security Site (NNSS), hosted a series of very shallow (<3 km) earthquakes in 1993. The RVFZ may also have hydrological significance within the NNSS, potentially playing a role in regional groundwater flow, but there is a lack of local hydrological data. In the Spring of 2021, we collected active-source accelerated weight drop seismic data over part of the RVFZ to better characterize the shallow subsurface. Here, we manually picked ~17,000 P-wave travel times and over 14,000 S-wave travel times, which were inverted for P-wave velocity (V P ), S-wave velocity (V S ), and V P /V S ratio in a 3D joint tomographic inversion scheme. Seismic velocities are imaged as deep as ~700 m in areas and generally align with geologic and structural expectations. V P and V S are relatively reduced near mapped and inferred faults, with the most prominent lower V P and V S zone around the densest collection of faults. We image V P /V S ratios ranging from ~1.5 to ~2.4, the extremes of which occur at a depth of ~100 m and are juxtaposed across a fault. One possible interpretation of the imaged seismic velocities is enhanced fault damage near the densest collection of faults with relatively higher porosity and/or crack density at ~100 m depth, with patches of semiperched groundwater present in the sedimentary rock in higher V P /V S areas and drier rock in lower V P /V S areas. A relatively higher V P /V S area beneath the densest faults persists at depth, which suggests percolation of groundwater via the fault damage zone to the regionally connected lower carbonate aquifer. Potentially, the presence and movement of groundwater may have played a role in the 1993 earthquake aftershocks.

58 GEOSCIENCES↗
DOE OSTI2020

Centimeter-Scale Pillared-Layer Metal–Organic Framework Thin Films Mediated by Hydroxy Double Salt Intermediates for CO 2 Sensor Applications

Fabrication of metal–organic framework (MOF) thin films over macroscopic surface areas is a subject of great interest for gas sensor application platforms such as optics and microelectronics. However, a direct synthesis of MOF films at ambient conditions, in particular pillared-layer MOF films due to their anisotropic structures, remains a significant challenge. Herein, we demonstrate for the first time a facile construction of dense and continuous pillared-layer MOF thin films on a centimeter scale via an aluminum-doped zinc oxide template and hydroxy double salt (HDS) intermediates at room temperature. A series of Cu(II)-based pillared MOFs with different 1,4-benzenedicarboxylic acid (bdc) ligands were explored for optimizing MOF film formation for CO 2 sensor applications. Nonpolar ligands with lower water solubility preferentially formed crystalline pillared MOF structures from HDS intermediates. A Cu 2 (ndc) 2 (dabco) (ndc = 1,4-naphthalene-bdc; dabco = 1,4-diazabicyclo[2.2.2]octane) MOF demonstrated the most dense and uniform film growth with micrometer thickness over one square centimeter area. This synthetic approach for growing Cu 2 (ndc) 2 (dabco) MOF thin films was successfully translated toward two sensing platforms: a quartz crystal microbalance and an optical fiber sensor. These Cu 2 (ndc) 2 (dabco) MOF-coated sensors displayed sensitivity toward CO 2 and response/recovery time on the scale of seconds, even at moderate humidity levels. Finally, this work provides a road map for producing continuous and anisotropic crystalline MOF thin films over a centimeter scale area on various substrates, which will greatly facilitate their utilization in MOF-based sensor devices, among other applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Laboratory Testing and Performance Characterization of the High Efficiency Dehumidification System (HEDS) Technology

The high efficiency dehumidification system (HEDS) is a unique air handling unit (AHU) design with two water coils. Connected in series, these water coils can sensibly reheat the airstream after it is cooled and dehumidified, thereby mitigating (or sometimes eliminating) the need for conventional sensible reheat with additional energy from hot water, electric resistance or steam coils. Another unique design feature is the size of the HEDS primary cooling coil. While typical AHU cooling coils have a 6-row geometry, the HED's cooling coil has a deeper 10-row geometry enabling discharge air dew point temperatures at or below 55 degrees F despite warmer chilled water supply temperatures approaching (sometimes exceeding) 50 degrees F. The ability to meet acceptable discharge air dew points with warmer chilled water supply temperatures improves both operational efficiency and system resilience against chilled water plant operational or control issues. This report summarizes the 59 steady-state test results that form a comprehensive performance map of the HEDS technology including its reheat capability. The breadth of outdoor air conditions and flow rates over which the HEDS technology was tested and characterized encompasses the unit 1) operating as a dedicated outdoor air system or mixed air system and 2) controlled with constant air volume or variable air volume fan control.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Magnetotail Variability During Magnetospheric Substorms

Abstract In this work, we present a statistical study of substorms covering a five‐year period 2016–2020. Substorm phases were identified from time series of the SuperMAG AL (SML) index using a list of 5,077 previously identified substorm onsets, the SML peak value marking transition from expansion to recovery phase, and the recovery identified as return to activity less than −100 nT in the SML index. Magnetic field observations from THEMIS, RBSP, and MMS missions were used to study the magnetotail characteristics during the substorm evolution. A superposed epoch analysis indicates that the substorm onset occurs almost simultaneously with a few minutes of uncertainty throughout the magnetotail, ranging from geostationary orbit to 20 RE. The onset in the transition region precedes the ground onset by a few minutes. The peak SML time coincides with the peak of the outer transition region Δ B Z , which suggests that the field‐aligned currents driving the SML activity arise from the outer transition region. Analysis of 2D maps of the tail magnetic field shows that the magnetotail current changes are limited to the center of the tail within | Y | < 10 R E . The substorm recovery is fastest in the inner transition region and lasts longer when moving further out. We did not find major asymmetries in the substorm signatures associated with IMF B Y or B Z .

Kumar, S.↗
DOE OSTI2026

MnRhBi3: A Cleavable Antiferromagnetic Metal

This dataset contains DFT input and output files supporting the theoretical modeling in the associated publication (Chem. Mater. 2024, 36, 11306-11316). The calculations characterize MnRhBi3, an orthorhombic (Cmmm) van der Waals-layered intermetallic compound that cleaves easily between neighboring Bi layers. The dataset is organized into three calculation types: (i) Bulk: Structural relaxations of the periodic MnRhBi3 crystal in antiferromagnetic (AFM) and ferromagnetic (FM) configurations, using the vdW-DF-optB86b functional. These provide the equilibrium lattice constants, magnetic energy differences (AFM is 0.5 meV/f.u. lower than FM), and magnetic moments (4.4 µB/Mn, 0.17 µB/Rh, 0.18 µB/Bi) reported in Table 1 of the main text. (ii) Slab: Same magnetic configurations computed with an 18 Ang vacuum layer introduced between Bi layers, used to calculate the cleavage energy Ec = 0.56 J/m2 (AFM) and 0.57 J/m2 (FM), establishing MnRhBi3 as a van der Waals-layered material comparable to graphite, MoS2, and CrI3. (iii) ELF: Single-point calculation on the relaxed bulk AFM geometry with LELF=.TRUE., producing the ELFCAR file used to generate electron localization function isosurfaces and contour maps (Fig. 2, main text) showing Bi lone pairs directed into the van der Waals gaps. All folders contain CONTCAR, INCAR, KPOINTS, OUTCAR, and POSCAR. The ELF/ folder additionally contains ELFCAR. Calculations were performed using VASP 6.3.2 with PBE + vdW-DF-optB86b, PAW potentials, and an energy cutoff of 800 eV.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Spatial Study 2021: Sensor-Based Time Series of Surface Water Temperature, Specific Conductance, Total Dissolved Solids, pH, and Dissolved Oxygen from across Multiple Watersheds in the Yakima River Basin, Washington, USA (v3)

This dataset supports a broader study examining the drivers of spatial variability in sediment respiration rates in the Yakima River Basin. The dataset provides two-hour time series hydrological and water chemistry sensor data, manual chamber open channel respiration data, handheld sensor water chemistry data, river substrate grain size photos, general environmental context photos, and field metadata (including qualitative information on instream and river corridor characteristics) collected during the same two-week period at 47 sites within multiple rivers throughout the Yakima River Basin in Washington, USA. Grain size photos can be used to improve estimates of channel substrate D50 data. Related sample-based water chemistry data are published separately at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1898914.This dataset is comprised of four main folders, one containing three sensor-specific subfolders and the others containing photographs. The SFA_SpatialStudy_2021_SensorData main data folder includes file-level metadata (FLMD), data dictionary (dd), installation methods, field metadata, Ultrameter water chemistry data, field data collection protocols, international generic sample number (IGSN) mapping file, and a readme file. The “Sensor_Manual_Specifications” subfolder contains pdf files from the manufacturer of each sensor with details on the sensor specifications. Each sensor subfolder (BarotrollAtm, MantaRiver, and MinidotManualChamber) contains a sensor data subfolder for timeseries data and a subfolder for plots and summary statistics. The BarotrollAtm Data subfolder contains In Situ Rugged BaroTROLL pressure and temperature data. The MantaRiver Data subfolder contains Eureka Manta+ 35B multisonde temperature, specific conductance, and pH data. The MinidotManualChamber Data subfolder contains PME MiniDOT Logger dissolved oxygen (mg/L and percent saturation) and temperature data. The folder SFA_SpatialStudy_2021_EnvironmentalContextPhotos contains environmental context photographs and videos. The folders SFA_SpatialStudy_2021_SedimentQuadratPhotos_Part1 and SFA_SpatialStudy_2021_SedimentQuadratPhotos_Part2 contain sediment quadrat photographs. All files are .csv, .pdf, .R, .jpg, .jpeg, .mp4, or .mov. This data package was originally published September 2022. It was updated January 2023 (modified files) and June 2024 (new and modified files). See the change history in data package readme for more details.We acknowledge the Yakama Nation as owners and caretakers of the lands where we collected these data. We thank the Confederated Tribes and Bands of the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate sample collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Isolating the Contributions of Specific Network Sites to the Diffuse Vibrational Spectrum of Interfacial Water with Isotopomer-Selective Spectroscopy of Cold Clusters

Decoding the structural information contained in the interfacial vibrational spectrum of water requires understanding how the spectral signatures of individual water molecules respond to their local hydrogen bonding environments. In this study, we isolated the contributions for the five classes of sites that differ according to the number of donor (D) and acceptor (A) hydrogen bonds that characterize each site. These patterns were measured by exploiting the unique properties of the water cluster cage structures formed in the gas phase upon hydration of a series of cations M+·(H2O)n (M = Li, Na, Cs, NH4, CH3NH3, H3O, and n = 5, 20-22). This selection of ions was chosen to systematically express the A, AD, AAD, ADD, and AADD hydrogen bonding motifs. The spectral signatures of each site were measured using two-color, IR-IR isotopomer-selective photofragmentation vibrational spectroscopy of the cryogenically cooled, mass selected cluster ions in which a single intact H2O is introduced without isotopic scrambling, an important advantage afforded by the cluster regime. The resulting patterns provide an unprecedented picture of the intrinsic line shapes and spectral complexities associated with excitation of the individual OH groups, as well as the correlation between the frequencies of the two OH groups on the same water molecule, as a function of network site. The properties of the surrounding water network that govern this frequency map are evaluated by dissecting electronic structure calculations that explore how changes in the nearby network structures, both within and beyond the first hydration shell, affect the local frequency of an OH oscillator. The qualitative trends are recovered with a simple model that correlates the OH frequency with the network-modulated local electron density in the center of the OH bond.

Yang, Nan↗
DOE OSTI2020

Regional Assessment for the CO 2 Storage Potential in the Northern Niagaran Pinnacle Reef Trend

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. The report describes the regional analysis of the Northern Niagaran Pinnacle Reef Trend (NNPRT) in Northern Michigan as part of the MRCSP programs Regional Assessment goals. The goal of this task was to perform an initial assessment of the geologic storage capacity and injectivity of the Northern Fairway of Michigan’s Niagaran Reef Trend. This was achieved through three main objectives: 1) development of a reef atlas, 2) understanding the regional trends in geology, and 3) estimating CO 2 and CO 2 -EOR resources. Data collected across multiple subtasks were integrated into a mappable database that could be accessed from many of the leading software. The reef atlas was used to map the reefs and associated characteristics across the entire NNPRT. Several methodologies were explored to estimate the CO 2 and CO 2 -EOR resources across the trend which included three scenarios: 1) storage only using fluid substitution and volumetric estimations, 2) CO 2 -EOR which applied proximity analysis concepts and measured performance metrics from the Core Energy reefs to predict performance at all oil reefs, and 3) enhanced storage scenario which combined CO 2 -EOR with fluid substitution to represent maximized storage after completion of CO 2 -EOR.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

Photometry and astrometry with JWST – III. A NIRCam-Gaia DR3 analysis of the open cluster NGC 2506

ABSTRACT In the third paper of this series aimed at developing the tools for analysing resolved stellar populations using the cameras on board of the James Webb Space Telescope (JWST), we present a detailed multiband study of the 2 Gyr Galactic open cluster NGC 2506. We employ public calibration data sets collected in multiple filters to: (i) derive improved effective Point Spread Functions (ePSFs) for 10 NIRCam filters; (ii) extract high-precision photometry and astrometry for stars in the cluster, approaching the main sequence (MS) lower mass of ∼0.1 M⊙; and (iii) take advantage of the synergy between JWST and Gaia DR3 to perform a comprehensive analysis of the cluster’s global and local properties. We derived a MS binary fraction of ∼57.5 per cent, extending the Gaia limit (∼0.8 M⊙) to lower masses (∼0.4 M⊙) with JWST. We conducted a study on the mass functions (MFs) of NGC 2506, mapping the mass segregation with Gaia data, and extending MFs to lower masses with the JWST field. We also combined information on the derived MFs to infer an estimate of the cluster present-day total mass. Lastly, we investigated the presence of white dwarfs (WDs) and identified a strong candidate. However, to firmly establish its cluster membership, as well as that of four other WD candidates and of the majority of faint low-mass MS stars, further JWST equally deep observations will be required. We make publicly available catalogues, atlases, and the improved ePSFs.

Nardiello, D. (ORCID:0000000311493659)↗
DOE OSTI2021

Reproduction of melting behavior for vitrified hillforts based on amphibolite, granite, and basalt lithologies

Abstract European Bronze and Iron Age vitrified hillforts have been known since the 1700s, but archaeological interpretations regarding their function and use are still debated. We carried out a series of experiments to constrain conditions that led to the vitrification of the inner wall rocks in the hillfort at Broborg, Sweden. Potential source rocks were collected locally and heat treated in the laboratory, varying maximum temperature, cooling rate, and starting particle size. Crystalline and amorphous phases were quantified using X-ray diffraction both in situ, during heating and cooling, and ex situ, after heating and quenching. Textures, phases, and glass compositions obtained were compared with those for rock samples from the vitrified part of the wall, as well as with equilibrium crystallization calculations. ‘Dark glass’ and its associated minerals formed from amphibolite or dolerite rocks melted at 1000–1200 °C under reducing atmosphere then slow cooled. ‘Clear glass’ formed from non-equilibrium partial melting of feldspar in granitoid rocks. This study aids archaeological forensic investigation of vitrified hillforts and interpretation of source rock material by mapping mineralogical changes and glass production under various heating conditions.

42 ENGINEERING↗
DOE OSTI2022

Deep learning to estimate permeability using geophysical data

Time-lapse electrical resistivity tomography (ERT) is a popular geophysical method to estimate three-dimensional (3D) permeability fields from electrical potential difference measurements. Traditional inversion and data assimilation methods are used to ingest this ERT data into hydrogeophysical models to estimate permeability. Due to ill-posedness and the curse of dimensionality, existing inversion strategies provide poor estimates and low resolution of the 3D permeability field. Recent advances in deep learning provide us with powerful algorithms to overcome this challenge. This paper presents a deep learning (DL) framework to estimate the 3D subsurface permeability from time-lapse ERT data. To test the feasibility of the proposed framework, we train DL-enabled inverse models on simulation data. Each measurement in both synthetic and field data is standardized by removing the mean and scaling the time-series to unit variance. This pre-processing step is necessary to bring simulation data closer to field observations. Subsurface process models based on hydrogeophysics are used to generate this synthetic data. Training performed on limited simulation data resulted in the DL model over-fitting. An advanced data augmentation based on mixup is implemented to generate additional training samples to overcome this issue. This mixup technique creates weakly labeled (low-fidelity) samples from strongly labeled (high-fidelity) data. The weakly labeled training data is then used to develop DL-enabled inverse models and reduce over-fitting. As both time-lapse ERT (1133048 features/realization) and 3D permeability (585453 features/realization) data samples are from a high-dimensional space, principal component analysis (PCA) is employed to reduce dimensionality. Encoded ERT and encoded permeability are generated using the trained PCA estimators. A deep neural network is then trained to map the encoded ERT to encoded permeability. This mixup training and unsupervised learning allowed us to build a fast and reasonably accurate DL-based inverse model under limited simulation data. Results show that proposed weak supervised learning can capture salient spatial features in the 3D permeability field. Quantitatively, the average mean squared error (in terms of the natural log) on the strongly labeled training, validation, and test datasets is less than 0.5. The R 2 -score (global metric) is greater than 0.75, and the percent error in each cell (local metric) is less than 10%. Finally, an added benefit in terms of computational cost is that the proposed DL-based inverse model is at least O(10 4 ) times faster than running a forward model once it is trained. Data generation, DL model training, and hyperparameter tuning to identify optimal neural network architectures utilized high-performance computing resources while the DL inference is performed on a standard laptop. Approximately, O(10 5 ) processor hours are used for generating data and DL tuning and training. We acknowledge that the data generation and DL model development are expensive. But once a DL model is trained, it can be re-used for inversion rapidly for the given system, with set physics and domain. Note that traditional inversion may require multiple forward model simulations (e.g., in the order of 10 to 1000), which are very expensive. This computational savings ≈ O(10 5 ) – O(10 7 )) makes the proposed DL-based inverse model attractive for subsurface imaging and real-time ERT monitoring applications due to fast and yet reasonably accurate estimations of permeability field.

58 GEOSCIENCES↗
DOE OSTI2022

Krylov subspace recycling for evolving structures

Krylov subspace recycling is a powerful tool when solving a long series of large, sparse linear systems that change only slowly over time. In PDE constrained shape optimization, these series appear naturally, as typically hundreds or thousands of optimization steps are needed with only small changes in the geometry. In this setting, however, applying Krylov subspace recycling can be a difficult task. As the geometry evolves, in general, so does the finite element mesh defined on or representing this geometry, including the numbers of nodes and elements and element connectivity. This is especially the case if re-meshing techniques are used. As a result, the number of algebraic degrees of freedom in the system changes, and in general the linear system matrices resulting from the finite element discretization change size from one optimization step to the next. Changes in the mesh connectivity also lead to structural changes in the matrices. In the case of re-meshing, even if the geometry changes only a little, the corresponding mesh might differ substantially from the previous one. Obviously, this prevents any straightforward mapping of the approximate invariant subspace of the linear system matrix (the focus of recycling in this work) from one optimization step to the next; similar problems arise for other selected subspaces. In this paper, we present an algorithm to map an approximate invariant subspace of the linear system matrix for the previous optimization step to an approximate invariant subspace of the linear system matrix for the current optimization step, for general meshes. This is achieved by exploiting the map from coefficient vectors to finite element functions on the mesh, combined with interpolation or approximation of functions on the finite element mesh. We demonstrate the effectiveness of our approach numerically with several proof of concept studies for a specific meshing technique.

42 ENGINEERING↗