Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ice physics”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Using distributions to understand neutron and x-ray production in ICF ignition capsules and other high energy density plasmas

Here this paper describes how x-ray and neutron distribution functions can be useful tools to visualize the conditions measured in many types of plasma physics experiments. In particular, we model a standard inertial confinement fusion ignition capsule that consists of a Si doped plastic ablator surrounding a layer of deuterium–tritium (DT) ice as the yield varies from 18 kJ to 16.7 MJ and use the distribution functions to show that neutrons and high energy x rays (15 keV) are produced under similar conditions when the yield is low. However, as the capsule starts to support a propagating burn due to alpha heating, the x rays and neutrons are produced under somewhat different conditions in different parts of the plasma. In particular, the x-ray production takes place mainly in the hot plastic ablator for the full yield ignition capsule under quite different plasma conditions from the DT region producing the 14 MeV neutrons, which results in x-ray images with larger radii than the corresponding neutron images. These same distribution functions can be applied to many other plasma physics experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Dark energy survey year 3 results: Cosmological constraints from cluster abundances, weak lensing, and galaxy clustering

Galaxy clusters provide a unique probe of the late-time cosmic structure and serve as a powerful independent test of the Λ⁢CDM model. This work presents the first set of cosmological constraints derived with ∼16,000 optically selected redMaPPer clusters across nearly 5000 deg 2 using DES year 3 datasets. Our analysis leverages a consistent modeling framework for galaxy cluster cosmology and DES-Y3 joint analyses of galaxy clustering and weak lensing (3 × 2⁢pt), ensuring direct comparability with the DES-Y3 3 × 2⁢pt analysis. Here, we obtain constraints of 𝑆 8 = 0.864 ± 0.035 and Ω m = $0.26⁢5^{+0.019}_{−0.031}$ from the cluster-based data vector. We find that cluster constraints and 3 × 2⁢pt constraints are consistent under the Λ⁢CDM model with a posterior predictive distribution (PPD) value of 0.53. The consistency between clusters and 3 × 2⁢pt provides a stringent test of Λ⁢CDM across different mass and spatial scales. Jointly analyzing clusters with 3 × 2⁢pt further improves cosmological constraints, yielding 𝑆 8 = $0.81⁢1^{+0.022}_{−0.020}$ and Ω m = $0.29⁢4^{+0.022}_{−0.033}$, a 24% improvement in the Ω m − 𝑆 8 figure of merit over 3 × 2⁢pt alone. Moreover, we find no significant deviation from the Planck CMB constraints with a probability to exceed (PTE) value of 0.6, significantly reducing previous 𝑆 8 tension claims. Finally, combining DES 3 × 2⁢pt, DES clusters, and Planck CMB places an upper limit on the sum of neutrino masses of ∑𝑚 𝜈 < 0.26 eV at 95% confidence under the Λ⁢CDM model. These results establish optically selected clusters as a key cosmological probe and pave the way for cluster-based analyses in upcoming stage-IV surveys such as LSST, Euclid, and Roman.

Abbott, T. M. C. [NSF’s National Optical-Infrared ↗
DOE OSTI2024

A multi-model assessment of the early last deglaciation (PMIP4 LDv1): a meltwater perspective

Abstract. The last deglaciation (∼20–11 ka BP) is a period of a major, long-term climate transition from a glacial to interglacial state that features multiple centennial- to decadal-scale abrupt climate variations whose root cause is still not fully understood. To better understand this time period, the Paleoclimate Modelling Intercomparison Project (PMIP) has provided a framework for an internationally coordinated endeavour in simulating the last deglaciation whilst encompassing a broad range of models. Here, we present a multi-model intercomparison of 17 transient simulations of the early part of the last deglaciation (∼20–15 ka BP) from nine different climate models spanning a range of model complexities and uncertain boundary conditions and forcings. The numerous simulations available provide the opportunity to better understand the chain of events and mechanisms of climate changes between 20 and 15 ka BP and our collective ability to simulate them. We conclude that the amount of freshwater forcing and whether it follows the ice sheet reconstruction or induces an inferred Atlantic meridional overturning circulation (AMOC) history, heavily impacts the deglacial climate evolution for each simulation rather than differences in the model physics. The course of the deglaciation is consistent between simulations except when the freshwater forcing is above 0.1 Sv – at least 70 % of the simulations agree that there is warming by 15 ka BP in most places excluding the location of meltwater input. For simulations with freshwater forcings that exceed 0.1 Sv from 18 ka BP, warming is delayed in the North Atlantic and surface air temperature correlations with AMOC strength are much higher. However, we find that the state of the AMOC coming out of the Last Glacial Maximum (LGM) also plays a key role in the AMOC sensitivity to model forcings. In addition, we show that the response of each model to the chosen meltwater scenario depends largely on the sensitivity of the model to the freshwater forcing and other aspects of the experimental design (e.g. CO2 forcing or ice sheet reconstruction). The results provide insight into the ability of our models to simulate the first part of the deglaciation and how choices between uncertain boundary conditions and forcings, with a focus on freshwater fluxes, can impact model outputs. We can use these findings as helpful insight in the design of future simulations of this time period.

Snoll, Brooke (ORCID:0000000230914394)↗
DOE OSTI2024

Refactoring the elastic–viscous–plastic solver from the sea ice model CICE v6.5.1 for improved performance

This study focuses on the performance of the elastic–viscous–plastic (EVP) dynamical solver within the sea ice model, CICE v6.5.1. The study has been conducted in two steps. First, the standard EVP solver was extracted from CICE for experiments with refactored versions, which are used for performance testing. Second, one refactored version was integrated and tested in the full CICE model to demonstrate that the new algorithms do not significantly impact the physical results. The study reveals two dominant bottlenecks, namely (1) the number of Message Parsing Interface (MPI) and Open Multi-Processing (OpenMP) synchronization points required for halo exchanges during each time step combined with the irregular domain of active sea ice points and (2) the lack of single-instruction, multiple-data (SIMD) code generation. The standard EVP solver has been refactored based on two generic patterns. The first pattern exposes how general finite differences on masked multi-dimensional arrays can be expressed in order to produce significantly better code generation by changing the memory access pattern from random access to direct access. The second pattern takes an alternative approach to handle static grid properties. The measured single-core performance improvement is more than a factor of 5 compared to the standard implementation. The refactored implementation of strong scales on the Intel® Xeon® Scalable Processors series node until the available bandwidth of the node is used. For the Intel® Xeon® CPU Max series, there is sufficient bandwidth to allow the strong scaling to continue for all the cores on the node, resulting in a single-node improvement factor of 35 over the standard implementation. This study also demonstrates improved performance on GPU processors.

58 GEOSCIENCES↗
DOE OSTI2020

Analysis of combustion acoustic phenomena in compression–ignition engines using large eddy simulation

As computational capabilities continue to grow, exploring the limits of computational fluid dynamics to capture complex and elusive phenomena, which are otherwise difficult to study by experimental techniques, is one of the main targets for the research community. This paper presents a detailed analysis of the physical processes that lead to combustion noise emissions in internal combustion engines. In particular, diesel combustion in a compression-ignition (CI) engine is studied in order to understand the singular behavior of the in-cylinder flow field responsible for the acoustic emissions. The main objective is, therefore, to improve the understanding of the phenomena involved in CI engine noise using large eddy simulations. Several visualization methods are employed to investigate the connection between combustion behavior and its effects on the pressure field. In addition, proper orthogonal decomposition is used to analyze the modal energy distribution among all the acoustic modes. The results show that the acoustic signature is fundamentally conditioned by the intensity of the premixed combustion rather than by the pressure oscillations generated by turbulent fluctuations in the flame surface established during the diffusion stage.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

High-speed impact of micron-sized diesel drop trains—Splashing dynamics, secondary droplet formation, and effects of pre-existing film thickness

This work investigates surface impingement of monodispersed diesel drop trains using computational simulations. The three-dimensional simulations include high-speed impact of micrometer-sized fuel drops onto initially dry and wetted stainless steel substrates. The drop size and impact velocity are representative of fuel injection conditions in internal combustion engines (ICEs). The drop trains serve as a simplified representation of fuel spray. To accurately capture the film that naturally develops on initially dry substrates, a dynamic diesel-stainless steel contact angle model was developed using experiments of single drop impact under ICE representative temperature. Simulations are presented for a highly splashing case and the effects of a pre-existing film on splashing dynamics were investigated, including the temporal evolution of splashed mass and film thickness. It was concluded that for thin films, the effect of pre-existing film thickness is less significant than in single drop impingement. Secondary droplet characterization was performed on simulation results, leading to instantaneous and time-averaged distributions of secondary droplet size, velocity magnitude, and trajectory angle. It was found that for each drop impingement, approximately 58% of the splashed mass is from the impinging drop itself, while the remainder of the splashed mass is composed of film liquid. From a high resolution case, details of secondary droplet formation are observed and three distinct phases of secondary droplet formation are identified. The detailed analysis of drop train impingement under engine-relevant conditions serves as a first step toward a robust understanding of fuel impingement and the development of cleaner and more efficient ICEs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Characterization of Orography-Influenced Riming and Secondary Ice Production and Their Effects on Precipitation Rates Using Radar Polarimetry and Doppler Spectra (CORSIPP-SAIL)

The Characterization of Orography-Influenced Riming and Secondary Ice Production and Their Effects on Precipitation Rates Using Radar Polarimetry and Doppler Spectra (CORSIPP) project was conducted to help improve the understanding of precipitation formation in orographically influenced terrain. Special focus is put on the two processes of riming and secondary ice production and their external drivers. Two instruments, a polarimetric W-band simultaneous transmission simultaneous reception (STSR) Doppler cloud radar manufactured by Radiometer Physics GmbH (RPG, instrument type RPG-FMCW-94-DP), from now on named LIMRAD94, and the video in situ snowfall sensor (VISSS), were deployed at the U. S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s Surface Atmosphere Integrated Field Laboratory (SAIL) site in Gothic, Colorado between November 2022 and June 2023 during the second SAIL winter. Note that the exact dates of data availability differ between the instruments. Both instruments arrived at Gothic on November 2, 2022, 09:20 local time. VISSS, described by Maahn et al., is equipped with two camera systems with telecentric lenses. The two cameras are at a 90° angle to each other. This configuration allows for size-independent measurements by capturing images of hydrometeors from two sides at a high frame rate of 250 Hz. With a minimum detection size of 200 μm, VISSS provides valuable insights into particle size, number, shape, complexity, and fall velocity. The VISSS was deployed on the grassland next to the ARM facility with the amazing help of the ARM employees on site. The setup started on November 2, 2022, and was finished on November 5, 2022, without major problems. VISSS measurements were started on November 6, 2022. LIMRAD94 was installed on a scaffolding platform near Orehouse (Gothic) on November 9, 2022, with the great help of RMBL staff. LIMRAD94 was mounted on a cold temperature scanner prototype. After a short test of the setup on November 9, 2022, the digital control of the scanner elevation stopped working for (at that time) unknown reasons. All attempts to resolve the problem failed. This malfunction made it impossible to operate LIMRAD94 in scanning mode. The scanner was then manually moved to zenith pointing mode and between November 10 and November 15, 2022, vertical observations for polarimetric calibration were performed. On November 15, 2022, after the polarimetric calibration was applied, the scanner was manually moved to 40° elevation with azimuthal view towards the Ka-band ARM Zenith Radar (KAZR) and measurements were continued at constant elevation. Investigation of the scanner malfunction on February 6, 2023, by Benn Schmatz revealed a disengagement between the cogwheel of the elevation motor and the cogwheel moving the scanner in elevation. This mechanical issue was temporarily solved by re-engaging the cogwheels. This made the scanner operational again for about four weeks, until mechanical force disengaged the cogwheels again on March 15, 2023. This repeated scanner failure remained undetected for about three weeks until April 8, 2023; during this time the scanner was stuck at 72° elevation. However, the radar software continued to produce data files falsely indicating that the scanner was still operational. After the scanner failure was noticed, scanning was stopped again, and we returned to constant elevation measurements. On May 15, 2023, the cogwheels of the elevation motor were secured with additional screws sent by the manufacturer. At some point in May, the cogwheels of the azimuth motor were also disengaged by mechanical force, which still allowed for range height indicator (RHI) but no plan position indicator (PPI) scans in the last weeks of the campaign. The azimuth motor was repaired in Germany after the end of the campaign. Throughout the campaign, RMBL and ARM staff kept the radar and the VISSS free of snow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Density functional theory of water with the machine-learned DM21 functional

The delicate interplay between functional-driven and density-driven errors in density functional theory (DFT) has hindered traditional density functional approximations (DFAs) from providing an accurate description of water for over 30 years. Recently, the deep-learned DeepMind 21 (DM21) functional has been shown to overcome the limitations of traditional DFAs as it is free of delocalization error. To determine if DM21 can enable a molecular-level description of the physical properties of aqueous systems within Kohn–Sham DFT, we assess the accuracy of the DM21 functional for neutral, protonated, and deprotonated water clusters. Furthermore, we find that the ability of DM21 to accurately predict the energetics of aqueous clusters varies significantly with cluster size. Additionally, we introduce the many-body MB-DM21 potential derived from DM21 data within the many-body expansion of the energy and use it in simulations of liquid water as a function of temperature at ambient pressure. We find that size-dependent functional-driven errors identified in the analysis of the energetics of small clusters calculated with the DM21 functional result in the MB-DM21 potential systematically overestimating the hydrogen-bond strength and, consequently, predicting a more ice-like local structure of water at room temperature.

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

The Climatic Impact‐Driver Framework for Assessment of Risk‐Relevant Climate Information

Abstract The climate science and applications communities need a broad and demand‐driven concept to assess physical climate conditions that are relevant for impacts on human and natural systems. Here, we augment the description of the “climatic impact‐driver” (CID) approach adopted in the Working Group I (WGI) contribution to the Intergovernmental Panel on Climate Change (IPCC) Sixth Assessment Report. CIDs are broadly defined as “physical climate system conditions (e.g., means, events, and extremes) that affect an element of society or ecosystems. Depending on system tolerance, CIDs and their changes can be detrimental, beneficial, neutral, or a mixture of each across interacting system elements and regions.” We give background information on the IPCC Report process that led to the development of the 7 CID types (heat and cold, wet and dry, wind, snow and ice, coastal, open ocean, and other) and 33 distinct CID categories, each of which may be evaluated using a variety of CID indices. This inventory of CIDs was co‐developed with WGII to provide a useful collaboration point between physical climate scientists and impacts/risk experts to assess the specific climatic phenomena driving sectoral responses and identify relevant CID indices within each sector. The CID Framework ensures that a comprehensive set of climatic conditions informs adaptation planning and risk management and may also help prioritize improvements in modeling sectoral dynamics that depend on climatic conditions. CIDs contribute to climate services by increasing coherence and neutrality when identifying and communicating relevant findings from physical climate research to risk assessment and planning activities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Biasing from galaxy trough and peak profiles with the DES Y3 redMaGiC galaxies and the weak lensing mass map

We measure the correspondence between the distribution of galaxies and matter around troughs and peaks in the projected galaxy density, by comparing redMaGiC galaxies (0.15 < z < 0.65) to weak lensing mass maps from the Dark Energy Survey (DES) Y3 data release. We obtain stacked profiles, as a function of angle θ, of the galaxy density contrast δ g and the weak lensing convergence κ, in the vicinity of these identified troughs and peaks, referred to as ‘void’ and ‘cluster’ superstructures. The ratio of the profiles depend mildly on θ, indicating good consistency between the profile shapes. We model the amplitude of this ratio using a function $F(\boldsymbol{\eta }, \theta )$ that depends on cosmological parameters $\boldsymbol{\eta }$, scaled by the galaxy bias. We construct templates of $F(\boldsymbol{\eta }, \theta )$ using a suite of N-body (‘Gower Street’) simulations forward-modelled with DES Y3-like noise and systematics. We discuss and quantify the caveats of using a linear bias model to create galaxy maps from the simulation dark matter shells. We measure the galaxy bias in three lens tomographic bins (near to far): $2.32^{+0.86}_{-0.27}, 2.18^{+0.86}_{-0.23}, 1.86^{+0.82}_{-0.23}$ for voids, and $2.46^{+0.73}_{-0.27}, 3.55^{+0.96}_{-0.55}, 4.27^{+0.36}_{-1.14}$ for clusters, assuming the best-fit Planck cosmology. Similar values with ∼0.1σ shifts are obtained assuming the mean DES Y3 cosmology. The biases from troughs and peaks are broadly consistent, although a larger bias is derived for peaks, which is also larger than those measured from the DES Y3 3 × 2-point analysis. This method shows an interesting avenue for measuring field-level bias that can be applied to future lensing surveys.

cosmology: observations↗
DOE OSTI2024

The importance of sub-meter-scale snow roughness on conductive heat flux of Arctic sea ice

Abstract The conductive heat flux through the snow and ice is a critical component of the mass and energy budgets in the Arctic sea ice system. We use high horizontal resolution (3–15 cm) measurements of snow topography to explore the impacts of sub-meter-scale snow surface roughness on heat flux as simulated by the Finite Element method. Simulating horizontal heat flux in a variable snow cover modestly increases the total simulated heat flux. With horizontal heat flux, as opposed to simple 1D-vertical heat flux modeling, the simulated heat flux is 10% greater than that for uniform snow with the same mean snow thickness for a 31.5 × 21 m region of sea ice (the largest region we studied). Vertical-only (1D) heat flux simulates just a 6% increase for the same region. However, this is highly dependent on observation resolution. Had we measured the snow cover at 1 m horizontal spacing or greater, simulating horizontal heat flux would not have changed the net heat flux from that simulated with vertical-only heat flux. These findings suggest that measuring and modeling snow roughness at sub-meter horizontal scales may be necessary to accurately represent horizontal heat flux on level Arctic sea ice.

Geology↗
DOE OSTI2023

A first-principles machine-learning force field for heterogeneous ice nucleation on microcline feldspar

The formation of ice in the atmosphere affects precipitation and cloud properties, and plays a key role in the climate of our planet. Although ice can form directly from liquid water at deeply supercooled conditions, the presence of foreign particles can aid ice formation at much warmer temperatures. Over the past decade, experiments have highlighted the remarkable efficiency of feldspar minerals as ice nuclei compared to other particles present in the atmosphere. However, the exact mechanism of ice formation on feldspar surfaces has yet to be fully understood. Here, we develop a first-principles machine-learning model for the potential energy surface aimed at studying ice nucleation at microcline feldspar surfaces. The model is able to reproduce with high-fidelity the energies and forces derived from density-functional theory (DFT) based on the SCAN exchange and correlation functional. Our training set includes configurations of bulk supercooled water, hexagonal and cubic ice, microcline, and fully-hydroxylated feldspar surfaces exposed to vacuum, liquid water, and ice. In this work, we apply the machine-learning force field to study different fully-hydroxylated terminations of the (100), (010), and (001) surfaces of microcline exposed to vacuum. Our calculations suggest that terminations that do not minimize the number of broken bonds are preferred in vacuum. We also study the structure of supercooled liquid water in contact with microcline surfaces, and find that water density correlations extend up to around 10 Å from the surfaces. Finally, we show that the force field maintains a high accuracy during the simulation of ice formation at microcline surfaces, even for large systems of around 30,000 atoms. Future work will be directed towards the calculation of nucleation free energy barriers and rates using the force field developed herein, and understanding the role of different microcline surfaces on ice nucleation.

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

Analysis of Dynamic Behavior of the Target and the Deuterium-Tritium Ice in Magnetic-Field Assisted Implosions

The use of strong magnetic fields to augment the output energy of inertial confinement fusion experiments at the National Ignition Facility is of high interest. It offers the potential of reducing electron thermal conduction and increasing hot-spot alpha heating with little to no change in hohlraum behavior. In these magnetically assisted ignition experiments, the ultimate goal is to add a B-field in the form of a pulse ranging from 25 to 60 T to a high-performing hohlraum implosion several microseconds before impingement of the laser beams. This requires eliminating metallic components in the target and replacing them with electrically nonconducting materials. However, the strong eddy currents generated by the rapidly increasing high B-field, which were calculated to be as high as 2000 K, can heat the hohlraum. Here in this paper, we examine the transient effects of this rapid temperature change on the behavior of the target as well as the fuel layer composed typically of deuterium and tritium. Using simulations and calculations for limiting case scenarios, we find that the effect of the heating is not restrictive toward the performance of the target or the quality of the deuterium-tritium ice.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Quantum Theory of Semiconductor Photo-Catalysis and Solar Water Splitting / Theory and Computation for Semiconductor Catalysis of Solar Water Splitting / Static and Dynamic Thermal Effects in Solids. Final Report

This is the final report. The finished research was published in 27 journal papers listed below. The topics fall roughly into seven categories: Photochemistry of water at the GaN ($10\bar{1}0$) surface; Quasiparticles: Band electrons and Anharmonic phonons; Quasiharmonic theory; Anharmonic lattice thermal conductivity; Water and Ice; Electron-phonon effects; and, Alloy theory.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Reduction of the type Ia supernova host galaxy step in the outer regions of galaxies

Using 1533 type Ia supernovae (SNe Ia) from the 5-yr sample of the Dark Energy Survey (DES), we investigate the relationship between the projected galactocentric separation of the SNe and their host galaxies and their light curves and standardization. We show, for the first time, that the difference in SN Ia post-standardization brightnesses between high- and low-mass hosts reduces from $0.078\pm 0.011$ mag in the full sample to $0.036 \pm 0.018$ mag for SNe Ia located in the outer regions of their host galaxies, while increasing to $0.100 \pm 0.014$ mag for SNe in the inner regions. The difference in the size of the mass step between inner and outer regions is $0.064\pm 0.023$ mag. In these inner regions, the step can be reduced (but not removed) using a model where the $R_V$ of dust along the line of sight to the SN changes as a function of galaxy properties. We investigate the remaining difference using the distributions of the SN Ia stretch parameter to test the inferred age of SN progenitors. Comparing red (older) environments only, outer regions have a higher proportion of high-stretch SNe and a more homogeneous stretch distribution. However, this effect cannot explain the reduction in significance of any Hubble residual step in outer regions. We conclude that the standardized distances of SNe Ia located in the outer regions of galaxies are less affected by their global host galaxy properties than those in the inner regions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The DOE E3SM Model Version 2: Overview of the Physical Model and Initial Model Evaluation

Abstract This work documents version two of the Department of Energy's Energy Exascale Earth System Model (E3SM). E3SMv2 is a significant evolution from its predecessor E3SMv1, resulting in a model that is nearly twice as fast and with a simulated climate that is improved in many metrics. We describe the physical climate model in its lower horizontal resolution configuration consisting of 110 km atmosphere, 165 km land, 0.5° river routing model, and an ocean and sea ice with mesh spacing varying between 60 km in the mid‐latitudes and 30 km at the equator and poles. The model performance is evaluated with Coupled Model Intercomparison Project Phase 6 Diagnosis, Evaluation, and Characterization of Klima simulations augmented with historical simulations as well as simulations to evaluate impacts of different forcing agents. The simulated climate has many realistic features of the climate system, with notable improvements in clouds and precipitation compared to E3SMv1. E3SMv1 suffered from an excessively high equilibrium climate sensitivity (ECS) of 5.3 K. In E3SMv2, ECS is reduced to 4.0 K which is now within the plausible range based on a recent World Climate Research Program assessment. However, a number of important biases remain including a weak Atlantic Meridional Overturning Circulation, deficiencies in the characteristics and spectral distribution of tropical atmospheric variability, and a significant underestimation of the observed warming in the second half of the historical period. An analysis of single‐forcing simulations indicates that correcting the historical temperature bias would require a substantial reduction in the magnitude of the aerosol‐related forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Mapping Glacier Basal Sliding Applying Machine Learning

During the RESOLVE project (“High-resolution imaging in subsurface geophysics: development of a multi-instrument platform for interdisciplinary research”), continuous surface displacement and seismic array observations were obtained on Glacier d’Argentière in the French Alps for 35 days in May 2018. The data set is used to perform a detailed study of targeted processes within the highly dynamic cryospheric environment. In particular, the physical processes controlling glacial basal motion are poorly understood and remain challenging to observe directly. Especially in the Alpine region for temperate based glaciers where the ice rapidly responds to changing climatic conditions and thus, processes are strongly intermittent in time and heterogeneous in space. Spatially dense seismic and Global Positioning System (GPS) measurements are analyzed applying machine learning to gain insight into the processes controlling glacial motions of Glacier d’Argentière. Using multiple bandpass-filtered copies of the continuous seismic waveforms, we compute energy-based features, develop a matched field beamforming catalog and include meteorological observations. Features describing the data are analyzed with a gradient boosting decision tree model to directly estimate the GPS displacements from the seismic noise. We posit that features of the seismic noise provide direct access to the dominant parameters that drive displacement on the highly variable and unsteady surface of the glacier. The machine learning model infers daily fluctuations and longer term trends. The results show on-ice displacement rates are strongly modulated by activity at the base of the glacier. The techniques presented provide a new approach to study glacial basal sliding and discover its full complexity.

58 GEOSCIENCES↗