Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “depth profiling techniques”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Development of a deep potential model for F and CF 2 etching of Si and SiO 2

An understanding of plasma-surface interactions at increasingly smaller scales is invaluable for the development of novel technologies and processing techniques. Molecular dynamics (MD) simulations can provide insights into atomic-scale interactions, though they are restricted by the availability of interatomic potentials. Machine learning methods, such as Deep Potential Molecular Dynamics (DeepMD), provide a systematic framework for the development of accurate and flexible ab initio-based models. In this work, we develop DeepMD models for the ion-enhanced etching of Si and SiO 2 by F and CF 2 radicals. We employ an active learning process to expand the data set on which the model is trained and demonstrate its effect on the model accuracy. The DeepMD results are compared to data from classical MD simulations and experiments. Physical sputtering yields of SiO 2 by Ar + ions show good agreement with previous simulation results using conventional interatomic potentials, though the predicted depth profiles are different. Etching yields are calculated as a function of ion energy and neutral to ion flux ratio for the Ar + ion-enhanced etching of SiO 2 and Si by F atoms, as well as for etching of SiO 2 by CF 2 radicals, showing reasonable agreement with experimental data. Finally, an ion-enhanced surface kinetic model is fitted to the DeepMD etch yields, and the fitted parameters are compared to quantities computed directly from DeepMD simulations. This study illustrates how molecular dynamics simulations using machine learning potentials can provide an accurate model of etching processes relevant to device manufacturing.

Kounis-Melas, Andreas [Princeton Univ., NJ (United↗
DOE OSTI2022

Mineralogical, magnetic and geochemical data constrain the pathways and extent of weathering of mineralized sedimentary rocks

The oxidative weathering of sulfidic rock can profoundly impact watersheds through the resulting export of acidity and metals. Weathering leaves a record of mineral transformation, particularly involving minor redox-sensitive phases, that can inform the development of conceptual and quantitative models. In sulfidic sedimentary rocks, however, variations in depositional history, diagenesis and mineralization can change or overprint the distributions of these trace minerals, complicating the interpretation of weathering signatures. Here we show that a combination of bulk mineralogical and geochemical techniques, micrometer-resolution X-ray fluorescence microprobe analysis and rock magnetic measurements, applied to drill core samples and single weathered fractures, can provide data that enable the development of a geochemically consistent weathering model. This work focused on one watershed in the Upper Colorado River Basin sitting within the Mesaverde Formation, a sedimentary sandstone bedrock with disseminated sulfide minerals, including pyrite and sphalerite, that were introduced during diagenesis and subsequent magmatic-hydrothermal mineralization. Combined analytical methods revealed the pathways of iron (Fe), carbonate and silicate mineral weathering and showed how pH controls element retention or release from the actively weathering fractured sandstone. Drill core logging, whole rock X-ray diffraction, and geochemical measurements document the progression from unweathered rock at depth to weathered rock at the surface. X-ray microprobe analyses of a 1-cm size weathering profile along a fracture surface are consistent with the mobilization of Fe(II) and Fe(III) into acidic pore water from the dissolution of primary pyrite, Fe-sphalerite, chlorite, and minor siderite and pyrrhotite. These reactions are followed by the precipitation of secondary minerals such as of goethite and jarosite, a Fe-(oxyhydr)oxide and hydrous Fe(III) sulfate, respectively. Microscale analyses also helped explain the weathering reactions responsible for the mineralogical transformations observed in the top and most weathered section of the drill core. For example, dissolution of feldspar and chlorite neutralizes the acidity generated by Fe and sulfide mineral oxidation, oversaturating the solution in both Fe-oxides. The combination of X-ray spectromicroscopy and magnetic measurements show that the Fe(III) product is goethite, mainly present either as a coating on fracture surfaces in the actively weathering region of the core or more homogeneously contained within the unconsolidated regolith at the top of the core. Low-temperature magnetic data reveal the presence of ferromagnetic Fe-sulfide pyrrhotite that, although it occurs at trace concentrations, could provide a qualitative proxy for unweathered sulfide minerals because the loss of pyrrhotite is associated with the onset of oxidative weathering. Pyrrhotite loss and goethite formation are detectable through room-temperature magnetic coercivity changes, suggesting that rock magnetic measurements can determine weathering intensity in rock samples at many scales. In conclusion, this work contributes evidence that the weathering of sulfidic sedimentary rocks follows a geochemical pattern in which the abundance of sulfide minerals controls the generation of acidity and dissolved elements, and the pH-dependent mobility of these elements controls their export to the ground- and surface-water.

58 GEOSCIENCES↗
DOE OSTI2024

Geophysical methods reveal the soil architecture and subsurface stratigraphic heterogeneities across land-lake interfaces along Lake Erie

Abstract Purpose The land-lake interface is a unique zone where terrestrial and aquatic ecosystems meet, forming part of the Earth’s most geochemically and biologically active zones. The unique characteristics of this interface are yet to be properly understood due to the inherently high spatiotemporal variability of subsurface properties, which are difficult to capture with the traditional soil sampling methods. Geophysical methods offer non-invasive techniques to capture variabilities in soil properties at a high resolution across various spatiotemporal scales. Methods We combined electromagnetic induction (EMI), electrical resistivity tomography (ERT), and ground penetrating radar (GPR) with data from soil cores and in situ sensors to investigate hydrostratigraphic heterogeneities across land-lake interfaces along the western basin of Lake Erie. Results EMI revealed high spatial heterogeneities in ECa distribution across the land-lake interfaces, with higher values in the wetland and transition zones compared to the upland zone. Soil ECa maps matched soil maps from a public database with the hydric soil units delineated as high conductivity zones (ECa > 40 mS/m). ERT and GPR showed vertical variation in soil properties with clear stratigraphic boundaries, and correlation of ERT profiles with lithologs from piezometers revealed the stratigraphic units of silt–clay and till sequence down to 3.5 m depth which are consistent with the surficial geology of the study area. Conclusions These results validate the use of multiple geophysical methods for extrapolating soil properties and mapping stratigraphic structures at land-lake interfaces, thereby providing the missing information required to improve the earth system model (ESM) of coastal interfaces.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Cybersecurity Threat Profile for a Connected Lighting System

In anticipation of improved energy performance and cost savings, cities and building owners are increasingly considering “smart lighting initiatives” that aim to convert their collection of simple luminaires (i.e., lighting fixtures) into an intelligent connected lighting system (CLS) capable of remotely monitoring energy consumption and fault conditions, and possibly implementing adaptive lighting schemes. The U.S. Department of Energy (DOE) has set an national goal of tripling the energy efficiency and demand flexibility of the buildings sector by 2030, relative to 2020 levels 1. It is forecast that connected lighting systems can contribute to that goal by delivering 125 TWh of annual energy savings by 2035 2, equivalent to the annual output of 50 typical (500 MW) power plants. However, these energy savings and the DOE goal are put at significant risk if connected technologies are not adopted due to real or perceived cybersecurity concerns. Connected IoT devices such as these have historically been rife with vulnerabilities which sometimes put security considerations secondary to functionality and operability. What are the cybersecurity threats that will impact these systems, as formerly banal luminaires transition into intelligent connected devices that collect information about themselves, their surrounding environment, and possibly us? In this paper we analyze a threat profile performed on a fault-detection use case for streetlights. A threat profile establishes security requirements, justifies security measures, yields actionable controls, and effectively communicates risk to stakeholders. This effort provides critical information for making threat-based decisions to increase security at a reasonable cost, and can effectively be used by development teams, software architects, and managers to make cybersecurity a part of their ongoing culture of awareness, training, and prevention. This leads to more secure systems and better-understood security. On-premise, cloud, and hybrid architectures with different authentication mechanisms were modeled and later categorized using the Microsoft STRIDE framework. An analysis of the recommended controls for each threat was performed to determine which controls could and should be put in place by manufacturers or third-party suppliers, and which controls need to be left up the end-user to implement. Fifty-seven threats were identified. Among our key findings: (1) 65% (37/57) of the threats did not involve the luminaires, but rather the other components needed to communicate with and manage them; (2) 63% (36/57) of the threats could have been mitigated through manufacturer-implemented defensive techniques or “controls”; and (3) 23% (13/57) of the threats were dependent on the network configuration. Recommendations based on the results of this work are made to key stakeholder groups. Notably, lighting technology developers are advised to address all threats that can be reasonably controlled with baked-in technology solutions (e.g., encryption or authentication controls), and employ some form of secure supply chain management and tracking where other parts (e.g., sensors, microprocessors) of a luminaire must also be built and manufactured with the proper security controls in place. Developers should also review threats involving assets not developed in-house to understand how connectivity with other devices will affect their product during system operation and determine if a compensating control for a defense-in-depth strategy will be needed. Finally, those interested in deploying CLS should compare the differences between cloud and on-premise models to determine which is more suitable for their needs and the abilities of their security team.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Play Fairway Analysis: Structurally Controlled Geothermal Systems in the Eastern Great Basin Extensional Regime, Utah

A research team with membership from the University of Utah/Energy & Geoscience Institute, the University of Utah/Dept. of Geology & Geophysics, and the Utah Geological Survey, undertook a play fairway analysis (PFA) for geothermal resources in the Eastern Great Basin (EGB) extensional tectonic regime of western Utah. This is a high-priority region for geothermal exploration because active Basin and Range (B&R) extension with volcanism having a N-S strike is superimposed upon pre-existing E-W belts of plutonic rocks and large-scale structural lineaments. Cumulative heat flow along the N-S strike of the state totals approximately 5 GWt above background stable interior. Three electricity producing power plants currently exist with substantial potential for increase. Succinctly, our PFA approach aims to resolve potential sources of heat and permeability in the region, which are the two principal criteria for establishing a geothermal resource. An initial Phase 1 was carried out using only existing geoscientific data in the area. Criteria selected for focusing heat potential include direct heat flow measurements in boreholes, magnetotelluric (MT) low resistivity anomalies, fluid/gas geochemistry, and proximity to recent volcanic eruptions. Permeability is established through geological structures (fault density, critically stressed areas, seismicity, gravity), and MT low resistivity anomalies. In Phase II of this PFA project, additional geological, geophysical and geochemical data were acquired and analysis carried out primarily over promising composite common risk (CCR) areas initially identified in Phase I in order to focus prospectivity and prepare for drilling recommendations. These prospects are near the Twin Peaks rhyolite field, high heat flow areas north of the producing Cove Fort system, and geophysical structure beneath the Crater Knoll area off the northeast flank of the Mineral Mountains. Additional data included MT site fill-in, structural mapping and analysis using high-resolution imagery, gravity and on-ground mapping, Nodal 3C passive seismic collection, and passive 3He surveying. Heat source and permeability potential are again expressed in terms of their individual common risk segment (CRS) maps, with a color scheme using green for most favorable (low risk) and red for least favorable (high risk). Diverse data types are united through the technique of probability kriging, which establishes prospectivity thresholds for each data type and then computes probability of exceeding that threshold over the PFA area. Modified CRS and CCRS maps are compared to those of Phase I to highlight tighter prospectivity focus. In doing so, the prospectivity threshold for heat was increased significantly to narrow the targeting. In the final Phase III of this project where a recommended deep thermal gradient hole was sited, additional geophysical, geochemical and geological field collection and analysis was carried out to refine drill hole targeting. This includes prospect-scale MT, gravity, structure, passive seismic deployment (Cove Fort area), and detailed 3He isotope profiling. It was the recommendation of the DOE Technical Monitoring Team (TMT) that one or more holes be sited in the north-ern Cove Fort area where legacy TG gradient holes showed high cumulative heat flow. These were to be of moderate depth, 2000-3000 feet, to reach the geothermal fluid table expected to start in excess of 1000 feet depth. The drilling organization stipulated by the DOE/GTO was that of the USGS Research Drilling Program (RDP) centered in Las Vegas, NV. A detailed well plan, appended to this report, was developed principally by Dr. Ben Barker consulting to University of Utah, Dr. Steve Pye on the DOE TMT, Mr. Steven Crawford of the USGS-RDP, and the project PI Phil Wannamaker. However, temperature and possible H 2 S at the systems lead to cancellation of the drilling last-minute as this appeared outside the experience base of the USGS-RDP. We hope to have the opportunity to revisit the test drilling and expand the Play Fairway Analysis of this region at some point in the future.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

High-Resolution Near-Surface Imaging at the Basin Scale Using Dark Fiber and Distributed Acoustic Sensing: Toward Site Effect Estimation in Urban Environments

Near-surface seismic structure, particularly the shear wave velocity (V s ), can strongly affect local site response, and should be accurately estimated for ground motion prediction during seismic hazard assessment. The Imperial Valley (California), occupying the southern end of the Salton Trough, is a seismically active basin with thick surficial lacustrine sedimentary deposits. In this study, we utilize ambient noise records and local earthquake events for high-resolution near-surface characterization and site effect estimation with an unlit fiber-optic telecommunication infrastructure (dark fiber) in Imperial Valley by using the distributed acoustic sensing (DAS) technique. We apply ambient noise interferometry to retrieve coherent surface waves from DAS records, and evaluate performances of three different surface wave methods on DAS ambient noise dispersion imaging. We develop a quality control workflow to improve the dispersion measurement of noisy portions of the DAS data set by using a data selection strategy. Using the joint inversion of both the fundamental mode and higher overtones of Rayleigh waves, a high resolution two-dimensional (2D) V s structure down to 70 m depth is obtained. We successfully achieve an improved V s 30 (the time-averaged shear-wave velocity in the top 30 m) model with higher spatial-resolution and reliability compared to the existing community model for the area. We also explore the potential for utilizing DAS earthquake events for site amplification estimation. The preliminary results reveal a clear anti-correlation between the approximated site response and the V s 30 profile. In conclusion, our results indicate the potential utility of DAS deployed on dark fiber for near-surface characterization in appropriate contexts.

58 GEOSCIENCES↗
DOE OSTI2021

Simultaneous x-ray and neutron 4D tomographic study of drying-driven hydro-mechanical behavior of cement-based materials at moderate temperatures

The drying of cement-based materials is intimately related to their durability, which has significant economic, social and environmental repercussions. The evolution of the saturation of the pore network and the associated drying shrinkage are in fact leading causes of cracking and of the ensuing penetration of aggressive chemicals. This process is highly heterogeneous, due to the thermo-hydric spatial gradients developing in the material from the exposed surfaces to its core and because of local effects, driven by the intrinsically heterogeneous micro-structure (e.g., by the distribution of pores and aggregates). It follows that macroscopic, sample-scale measurements cannot fully disclose the complexity of the underlying processes. In the last few decades, significant advances in full-field techniques have allowed an unprecedented insight into these local processes. For cement-based materials, x-ray and neutron tomography lend themselves as ideal, and highly complementary, tools for the study of their thermo-hydro-mechanical behavior. Notably, the high sensitivity to density variations of x-ray imaging gives access to the developments of fractures, in 4D (3D + time). On the other hand, neutron tomography allows the study of the evolution of the moisture field in 4D, thanks to its high hydrogen sensitivity. The combination of these two techniques provides a unique insight in thermo-hydro-mechanical couplings, e.g., the effect of cracks on the water content field. This contribution presents novel 5D datasets (3D tomographies along time, plus truly simultaneous x-ray and Neutron rapid acquisitions) in-operando of a cement paste and of a concrete sample heated at moderate temperatures (up to 140{sup ∘}C). The analysis of this 5D data-set (once aligned in time and across modalities) allows for example a quantification of the 4D moisture profiles which were found to predict an overall water loss at hydric equilibrium coherent with the corresponding analytical analysis. In the cement paste sample, the x-ray dataset captures the evolution of an extensive cracking network, opening and propagation toward the core of the sample. A novel analysis procedure is here proposed which allows the extraction of these fractures and the analysis of their interplay with local drying as captured through neutron imaging. This for example reveals the depth of penetration of drying in the vicinity of the fractures along time, which is essential for the assessment and calibration of hydro-mechanical coupled models.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modeling of streamflow in a 30 km long reach spanning 5 years using OpenFOAM 5.x

Abstract. Developing accurate and efficient modeling techniques for streamflow at the tens-of-kilometers spatial scale and multi-year temporal scale is critical for evaluating and predicting the impact of climate- and human-induced discharge variations on river hydrodynamics. However, achieving such a goal is challenging because of limited surveys of streambed hydraulic roughness, uncertain boundary condition specifications, and high computational costs. We demonstrate that accurate and efficient three-dimensional (3-D) hydrodynamic modeling of natural rivers at 30 km and 5-year scales is feasible using the following three techniques within OpenFOAM, an open-source computational fluid dynamics platform: (1) generating a distributed hydraulic roughness field for the streambed by integrating water-stage observation data, a rough wall theory, and a local roughness optimization and adjustment strategy; (2) prescribing the boundary condition for the inflow and outflow by integrating precomputed results of a one-dimensional (1-D) hydraulic model with the 3-D model; and (3) reducing computational time using multiple parallel runs constrained by 1-D inflow and outflow boundary conditions. Streamflow modeling for a 30 km long reach in the Columbia River (CR) over 58 months can be achieved in less than 6 d using 1.1 million CPU hours. The mean error between the modeled and the observed water stages for our simulated CR reach ranges from −16 to 9 cm (equivalent to approximately ±7 % relative to the average water depth) at seven locations during most of the years between 2011 and 2019. We can reproduce the velocity distribution measured by the acoustic Doppler current profiler (ADCP). The correlation coefficients of the depth-averaged velocity between the model and ADCP measurements are in the range between 0.71 and 0.83 at 75 % of the survey cross sections. With the validated model, we further show that the relative importance of dynamic pressure versus hydrostatic pressure varies with discharge variations and topography heterogeneity. Given the model's high accuracy and computational efficiency, the model framework provides a generic approach to evaluate and predict the impacts of climate- and human-induced discharge variations on river hydrodynamics at tens-of-kilometers and decadal scales.

58 GEOSCIENCES↗
DOE OSTI2021

Crystal plasticity model of residual stress in additive manufacturing using the element elimination and reactivation method

Selective laser melting is receiving increasing interest as an additive manufacturing technique. Residual stresses induced by the large temperature gradients and inhomogeneous cooling process can favour the generation of cracks. In this work, a crystal plasticity finite element model is developed to simulate the formation of residual stresses and to understand the correlation between plastic deformation, grain orientation and residual stresses in the additive manufacturing process. The temperature profile and grain structure from thermal-fluid flow and grain growth simulations are implemented into the crystal plasticity model. An element elimination and reactivation method is proposed to model the melting and solidification and to reinitialize state variables, such as the plastic deformation, in the reactivated elements. The accuracy of this method is judged against previous method based on the stiffness degradation of liquid regions by comparing the plastic deformation as a function of time induced by thermal stresses. The method is used to investigate residual stresses parallel and perpendicular to the laser scan direction, and the correlation with the maximum Schmid factor of the grains along those directions. The magnitude of the residual stress can be predicted as a function of the depth, grain orientation and position with respect to the molten pool. The simulation results are directly comparable to X-ray diffraction experiments and stress–strain curves.

316 stainless steel↗
DOE OSTI2026

Quasi-In-Situ Analysis of Electrode Top Atomic Layers via High-Sensitivity Low-Energy Ion Scattering and Potential-Controlled Sample Transfer

Electrocatalytic reactions involve interfacial interactions between the surfaces of electrodes and reactive species at an electrolyte interface. There are presently no universal or unambiguous methods to directly assay the active top atomic layer composition that influences the reactivity of these electrodes under relevant operating conditions. Low-energy ion scattering (LEIS) spectroscopy is a surface characterization technique that yields compositional analysis of the outermost atomic layer of a material, but it must be performed in ultrahigh vacuum (UHV). Application of LEIS measurements to electrochemical materials that are removed from ambient liquid-phase environments thus leaves an open question as to whether the surface that is transferred to UHV is truly the surface that manifested during the electrochemical reaction. Toward the goal of preserving the active surface state, we developed a sample transfer workflow for LEIS enabling air-free removal and drying of an electrode from an electrochemical cell while maintaining control of the potential using an auxiliary electrode. The potential-controlled emersion method was demonstrated to give distinct potential-dependent surface compositions for a Cu−Pd alloy relative to removal after uncontrolled return to open-circuit potential. A Cu-enriched surface was found at anodic potential and a Pd-enriched surface at cathodic potential, suggesting that the approach can be used to retain representative atomic configurations during transfer. Since adsorbates will often persist from the reaction environment, conventional sample pretreatment methods for removal, including atomic O and atomic H exposure, were also contrasted. Both methods were found to differ with results from incidental low-dose depth profiling by the LEIS primary ion source, which removes adventitious species and surface atoms during the course of repeated measurements. These depth profiles were found to be sensitive to sample history and thus qualitatively informative, despite the possible changes induced by ion damage. The results exhibit (i) the need for complete control over the polarization state of the sample at all times (no excursions to open circuit during transfer) and (ii) the utility of low-dose depth profiling to capture changes in the near-surface composition.

Alloys↗
DOE OSTI2020

A micromirror array with annular partitioning for high-speed random-access axial focusing

Dynamic axial focusing functionality has recently experienced widespread incorporation in microscopy, augmented/virtual reality (AR/VR), adaptive optics and material processing. However, the limitations of existing varifocal tools continue to beset the performance capabilities and operating overhead of the optical systems that mobilize such functionality. The varifocal tools that are the least burdensome to operate (e.g. liquid crystal, elastomeric or optofluidic lenses) suffer from low (≈100 Hz) refresh rates. Conversely, the fastest devices sacrifice either critical capabilities such as their dwelling capacity (e.g. acoustic gradient lenses or monolithic micromechanical mirrors) or low operating overhead (e.g. deformable mirrors). Here, we present a general-purpose random-access axial focusing device that bridges these previously conflicting features of high speed, dwelling capacity and lightweight drive by employing low-rigidity micromirrors that exploit the robustness of defocusing phase profiles. Geometrically, the device consists of an 8.2 mm diameter array of piston-motion and 48-μm-pitch micromirror pixels that provide 2π phase shifting for wavelengths shorter than 1100 nm with 10–90% settling in 64.8 μs (i.e., 15.44 kHz refresh rate). The pixels are electrically partitioned into 32 rings for a driving scheme that enables phase-wrapped operation with circular symmetry and requires <30 V per channel. Optical experiments demonstrated the array’s wide focusing range with a measured ability to target 29 distinct resolvable depth planes. Overall, the features of the proposed array offer the potential for compact, straightforward methods of tackling bottlenecked applications, including high-throughput single-cell targeting in neurobiology and the delivery of dense 3D visual information in AR/VR.

42 ENGINEERING↗
DOE OSTI2024

ToF-SIMS in material research: A view from nanoscale hydrogen detection

Hydrogen in materials has attracted tremendous interest as its incorporation leads to significant alterations in nanoscale structure, composition, and chemistry, impacting functional properties. It has also been integral to nuclear fusion reactors and is considered a future clean energy source. However, nanoscale characterization and manipulation of hydrogen in materials are challenging as only a selected few analytical techniques can readily detect hydrogen, among which time-of-flight secondary ion mass spectrometry (ToF-SIMS) is a unique and powerful one due to its excellent detection limit along with decent depth and lateral resolutions. Here, in this review, we discuss, using selected examples, how to detect and quantify hydrogen in materials by ToF-SIMS and its impact on revealing the hydrogenation/protonation-induced novel functional states in different classes of materials. In addition, we present our protocols on sample preparation and experimental conditions optimization, allowing us to achieve the best possible results. Finally, we highlight future research directions that can lead to the discovery of novel functional states and ultimately provide a deeper understanding of scientific questions in materials science.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Spatial evidence of cryptic methane cycling and methylotrophic metabolisms along a land–ocean transect in salt marsh sediment

Methylotrophic methanogenesis in the sulfate-rich zone of coastal and marine sediments couples with anaerobic oxidation of methane (AOM), forming the cryptic methane cycle. This study provides evidence of cryptic methane cycling in the sulfate-rich zone across a land–ocean transect of four stations–two brackish, one marine, and one hypersaline–within the Carpinteria Salt Marsh Reserve (CSMR), southern California, USA. Samples from the top 20 cm of sediment from the transect were analyzed through geochemical and molecular (16S rRNA) techniques, in-vitro methanogenesis incubations, and radiotracer incubations utilizing 35 S-SO 4 , 14 C-mono-methylamine, and 14 C-CH 4 . Sediment methane concentrations were consistently low (3 to 28 µM) at all stations, except for the marine station, where methane increased with depth reaching 665 µM. Methanogenesis from mono-methylamine was detected throughout the sediment at all stations with estimated CH 4 production rates in the sub-nanomolar to nanomolar range per cm 3 sediment and day. 16S rRNA analysis identified methanogenic archaea (Methanosarcinaceae, Methanomassiliicoccales, and Methanonatronarchaeacea) capable of producing methane from methylamines in sediment where methylotrophic methanogenesis was found to be active. Metabolomic analysis of porewater showed mono-methylamine was mostly undetectable (<3 µM) or present in trace amounts (<10 µM) suggesting rapid metabolic turnover. In-vitro methanogenesis incubations of natural sediment showed no linear methane buildup, suggesting a process limiting methane emissions. AOM activity, measured with 14 C-CH 4 , overlapped with methanogenesis from mono-methylamine activity at all stations, with rates ranging from 0.03 to 19.4 nmol cm −3 d −1 . Geochemical porewater analysis showed the CSMR sediments are rich in sulfate and iron. Porewater sulfate concentrations (9–91 mM) were non-limiting across the transect, supporting sulfate reduction activity (1.5–2,506 nmol cm −3 d −1 ). Porewater sulfide and iron (II) profiles indicated that the sediment transitioned from a predominantly iron-reducing environment at the two brackish stations to a predominantly sulfate-reducing environment at the marine and hypersaline stations, which coincided with the presence of phyla (Desulfobacterota) involved in these processes. AOM activity overlapped with sulfate reduction and porewater iron (II) concentrations suggesting that AOM is likely coupled to sulfate and possibly iron reduction at all stations. However, 16S rRNA analysis identified anaerobic methanotrophs (ANME-2) only at the marine and hypersaline stations while putative methanogens were found in sediment across all stations. In one sediment horizon at the marine station, methanogen families (Methanosarcinaceae, Methanosaetaceae, Methanomassiliicoccales, and Methanoregulaceae) and ANME 2a,2b, and 2c groups were found together. Collectively, our data suggest that at the brackish stations methanogens alone may be involved in cryptic methane cycling, while at the marine and hypersaline stations both groups may be involved in the process. Differences in rate constants from incubations with 14 C-labeled methane and mono-methylamine suggest a non-methanogenic process oxidizing mono-methylamine to inorganic carbon, likely mediated by sulfate-reducing bacteria. Understanding the potential competition of sulfate reducers with methanogens for mono-methylamine needs further investigation as it might be another important process responsible for low methane emissions in salt marshes.

Anaerobic oxidation of methane↗
DOE OSTI2022

Variational encoder geostatistical analysis (VEGAS) with an application to large scale riverine bathymetry

Estimation of riverbed profiles, also known as bathymetry, plays a vital role in many applications, such as safe and efficient inland navigation, prediction of bank erosion, land subsidence, and flood risk management. The high cost and complex logistics of direct bathymetry surveys, i.e, depth imaging, have encouraged the use of indirect measurements such as surface flow velocities. However, estimating high-resolution bathymetry from indirect measurements is an inverse problem that can be computationally challenging. Here, we propose a reduced-order model (ROM) based approach that utilizes a variational autoencoder (VAE), a type of deep neural network with a narrow layer in the middle, to compress bathymetry and flow velocity information and accelerate bathymetry inverse problems from flow velocity measurements. In our application, the shallow-water equations (SWE) with appropriate boundary conditions (BCs), e.g., the discharge and/or the free surface elevation, constitute the forward problem, to predict flow velocity. Then, ROMs of the SWEs are constructed on a nonlinear manifold of low dimensionality through a variational encoder and the bathymetry inversion problem is derived on the low-dimensional latent space in a Hierarchical Bayesian setting. Further, the reformulation allows variational inference with a small number (e.g., $\mathscr{O}$ (100) of ROM runs and efficient uncertainty quantification. We have tested our inversion approach on a one-mile reach of the Savannah River, GA, USA. Once the neural network is trained (offline stage), the proposed technique can perform the inversion operation orders of magnitude faster than traditional inversion methods that are commonly based on linear projections, such as principal component analysis (PCA), or the principal component geostatistical approach (PCGA). Furthermore, tests show that the algorithm can estimate the bathymetry with good accuracy even with sparse flow velocity measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Generation of Bright Soft X-ray Laser Beams (Final Report)

The project goal was to demonstrate new compact soft x-ray lasers emitting high energy ultrashort pulses at shorter wavelengths and increased repetition rates by amplification in atomic transitions in high density laser-created plasmas. To demonstrate these lasers we developed new high power solid state lasers. The project succeeded in demonstrating the shortest wavelength compact soft x-ray laser to date, 6.85 nm wavelength, and the highest repetition rate table-top soft x-ray laser, 400 Hz. We also applied these compact soft x-ray laser sources to determine the composition of materials with nanoscale resolution using a new technique that combines soft x-ray laser ablation with mass spectrometry, and we demonstrated that these lasers are capable of acquiring nanoscale resolution holograms with picosecond time resolution. The combination of an increased pumping efficiency with increased repetition rates can yield soft x-ray laser beams with an unprecedented average power on a table-top for applications demanding a high photon flux. We have demonstrated that the efficiency of plasma-based laser-pumped soft x-ray lasers can be significantly increased by detailed pulse shaping of the pump laser using a new pulse synthesizer. Results include the demonstration of 18.9 nm wavelength laser operation at a record 400 Hz repetition rate, and the extension of repetitive gain-saturated plasma-based table-top x-ray lasers to the shortest wavelength to date, 6.85 nm. Laser amplification was also observed at wavelengths as short at 5.8 nm by isolectronic scaling in Ni-like lanthanide ions. We have also conducted simulations to explore the possibility of demonstrating gain saturated atomic x-ray lasers near the carbon K-edge, in the “carbon window” for the first time. The work conducted under this program also includes the use of compact soft x-ray lasers in applications to problems of chemical interest, such as isotopic composition mapping of micron-size uranium particles with nanoscale resolution, and nanoscale resolution depth profiling analysis of CoNCN-coated electrodes for water oxidation catalysis. We also demonstrated single-shot high resolution soft x-ray Fourier transform holography over a broad 7 µm diameter field of view with ~ 5 ps temporal resolution.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Satellite-based evaluation of AeroCom model bias in biomass burning regions

Global models are widely used to simulate biomass burning aerosol (BBA). Exhaustive evaluations on model representation of aerosol distributions and properties are fundamental to assess health and climate impacts of BBA. Here we conducted a comprehensive comparison of Aerosol Comparisons between Observations and Models (AeroCom) project model simulations with satellite observations. A total of 59 runs by 18 models from three AeroCom Phase-III experiments (i.e., biomass burning emissions, CTRL16, and CTRL19) and 14 satellite products of aerosols were used in the study. Aerosol optical depth (AOD) at 550 nm was investigated during the fire season over three key fire regions reflecting different fire dynamics (i.e., deforestation-dominated Amazon, Southern Hemisphere Africa where savannas are the key source of emissions, and boreal forest burning in boreal North America). The 14 satellite products were first evaluated against AErosol RObotic NETwork (AERONET) observations, with large uncertainties found. But these uncertainties had small impacts on the model evaluation that was dominated by modeling bias. Through a comparison with Polarization and Directionality of the Earth’s Reflectances measurements with the Generalized Retrieval of Aerosol and Surface Properties algorithm (POLDER-GRASP), we found that the modeled AOD values were biased by –93 % to 152 %, with most models showing significant underestimations even for the state-of-the-art aerosol modeling techniques (i.e., CTRL19). By scaling up BBA emissions, the negative biases in modeled AOD were significantly mitigated, although it yielded only negligible improvements in the correlation between models and observations, and the spatial and temporal variations in AOD biases did not change much. For models in CTRL16 and CTRL19, the large diversity in modeled AOD was in almost equal measures caused by diversity in emissions, lifetime, and the mass extinction coefficient (MEC). We found that in the AeroCom ensemble, BBA lifetime correlated significantly with particle deposition (as expected) and in turn correlated strongly with precipitation. Additional analysis based on Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP) aerosol profiles suggested that the altitude of the aerosol layer in the current models was generally too low, which also contributed to the bias in modeled lifetime. Modeled MECs exhibited significant correlations with the Ångström exponent (AE, an indicator of particle size). Comparisons with the POLDER-GRASP-observed AE suggested that the models tended to overestimate the AE (underestimated particle size), indicating a possible underestimation of MECs in models. The hygroscopic growth in most models generally agreed with observations and might not explain the overall underestimation of modeled AOD. Our results imply that current global models contain biases in important aerosol processes for BBA (e.g., emissions, removal, and optical properties) that remain to be addressed in future research.

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

Predicting cutoff L-shells of solar protons using the GPPSn particle dataset

Solar energetic protons (SEPs) arriving at the Earth trigger severe radiation storms in the near-Earth space, directly impacting space missions operating at various altitudes. Therefore, monitoring SEP events and predicting the penetration depths of solar protons are critical for aerospace sectors. Building on previous efforts, here we demonstrate the feasibility of using proton measurements from the Global Prompt Proton Sensor network (GPPSn), enabled by Los Alamos National Laboratory developed combined X-ray dosimeters aboard GPS satellites, to characterize and predict the penetration of solar protons into the geomagnetic field. The inclined medium-Earth-orbits (MEOs) of the global GPS constellation offer a unique advantage of allowing simultaneous measurements of penetrating solar protons inside both open- and closed-field line regions. Therefore, the L-profiles of ∼10s–100 MeV solar protons and their associated cutoff L-shells can be determined from the GPPSn dataset, using predefined threshold proton flux values rather than traditional flux ratios. After examining a list of SEP event intervals across solar cycles 23, 24 and 25—including the 2024 Mother’s Day superstorm, we showcase how the latest GPPSn proton dataset (release v1.10), reprocessed and calibrated, can not only be used to monitor solar proton distributions inside the dynamic geomagnetic field for individual events, but also to derive a new empirical model linking cutoff L-shells with several key space weather parameters. This newly developed SEPCL-MEO model demonstrates high predictive performance; for example, predictions for > 30 MeV solar protons yield a correlation coefficient of 0.85 and performance efficiency of 0.67 when validated against GPPSn observations. Results from this pilot study underscores the scientific and operational value of the GPPSn dataset, and this dataset—when paired with machine-learning techniques—can play a critical role in observing and predicting the effects of future incoming SEP events, including extreme ones.

58 GEOSCIENCES↗
DOE OSTI2021

Realization of a high luminosity muon beam at the Jefferson Lab

The proton is the major component of the visible matter of our universe, so being able to determine its fundamental properties, such as its charge radius, is of vital importance. In recent years, several experiments have been performed for the accurate measurement of the proton radius with different experimental techniques. The main and most common are based on two different techniques: ep scattering and Lamb shift spectroscopy of regular hydrogen atoms. The latter can also be performed on muonic hydrogen atoms with a significant improvement in accuracy. Recent measurements that use this method have highlighted a discrepancy between the results obtained with the diffusion of electrons and those obtained with the spectroscopic technique leading to the so-called ?proton radius puzzle?. The explanation of why this difference between the measured values exists is still a matter of discussion. The last attempt to address the puzzle was the PRad Experiment. It was performed with a unique experimental setup that gave the scientists an excellent control on the systematic errors that are usually present in electron-proton scattering experiments. This has significantly increased the precision of the measurement, bringing it in the same area as the spectroscopic ones done on muonic hydrogen atoms, which are considered the most precise. However, to solve the puzzle once and for all, even more precise electron scattering measurements are needed. Alternatively, new methods can be used, such as eg. muon-proton scattering at low transferred momenta in order to lower the systematic errors that comes with the measurement. This thesis tries to address this problem by studying the set up of an intense muon beam at the Jefferson Laboratory (JLab) that can be used to carry out a high-precision m ? p scattering experiment. Using Montecarlo simulations I studied the muon production resulting from the interaction of the primary electron beam with the beam-dump at JLab Hall-A. The simulations were performed using the FLUKA package via the FLAIR graphical interface. In this study, I modified an existing geometry (built for the ?Beam Dump eXperiment?) and performed several studies in order to optimize the muon beam parameters. In order to estimate the beam intensity, I studied the muon generation and transportation through the concrete bunker surrounding the Hall-A beam-dump. Detailed studies on muon attenuation and dispersion at different depth allowed me to track the beam profile in space and momentum as a function of the distance from the beam-dump. Last, but not least, I studied a possible tagging system for the muon beam. The tagging system will allow us to know with good precision (O(1=1000)) the energy of each muon of the beam. Based on the results obtained, I created an event generator that incorporates the detailed description of the muon beam (energy, position, emittance) useful to define and optimize the experimental setup for a future measurement of the proton radius at the JLab.

Fulci, Antonino↗