Search NASA⌕ Search

SEARCH · Search NASA

Results for “liquid water”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Insights into the relation between vertical cloud structure and dynamics of three polar lows: Observations from COMBLE

The vertical structure of three polar lows is described using profiling radar, lidar, and passive remote sensors deployed at a coastal site in northern Norway. The data were collected as part of the Cold-air Outbreaks in the Marine Boundary Layer Experiment (COMBLE). These data are examined in the context of satellite imagery, operational weather radar reflectivity imagery, and output from the AROME-Arctic model. A comparison to satellite images shows that AROME-Arctic captured the polar lows well. All three polar lows form in marine cold-air outbreaks and are surrounded by open-cellular convection typical of such outbreaks. Two of the lows form in a forward shear environment, and one under reverse shear. Initially, the three polar lows are comma-shaped; two of them transition to be more spirali-form at their mature stage and have mostly cloud-free warm cores. The warm cores are the result of a warm-seclusion process. All three lows have stratiform precipitation bands, marked by little cloud liquid water, and rather high surface precipitation rate. The vertical drafts and turbulence in these stratiform clouds are generally weak. In conclusion, all three lows also have convective clouds, which have stronger vertical drafts, stronger turbulence, and pockets of high liquid water content.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Glaciation of mixed-phase clouds: insights from bulk model and bin-microphysics large-eddy simulation informed by laboratory experiment

Abstract. Mixed-phase clouds affect precipitation and radiation differently from liquid and ice clouds, posing greater challenges to their representation in numerical simulations. Recent laboratory experiments using the Pi Cloud Chamber explored cloud glaciation conditions based on increased injection of ice-nucleating particles. In this study, we use two approaches to reproduce the results of the laboratory experiments: a bulk scalar mixing model and large-eddy simulation (LES) with bin microphysics. The first approach assumes a well-mixed domain to provide an efficient assessment of the mean cloud properties for a wide range of conditions. The second approach resolves the energy-carrying turbulence, the particle size distribution, and their spatial distribution to provide more details. These modeling approaches enable a separate and detailed examination of liquid and ice properties, which is challenging in the laboratory. Both approaches demonstrate that, with an increased ice number concentration, the flow and microphysical properties exhibit the same changes in trends. Additionally, both approaches show that the ice integral radius reaches the theoretical glaciation threshold when the cloud is subsaturated with respect to liquid water. The main difference between the results of the two approaches is that the bulk model allows for the complete glaciation of the cloud. However, LES reveals that, in a dynamic system, the cloud is not completely glaciated as liquid water droplets are continuously produced near the warm lower boundary and subsequently mixed into the chamber interior. These results highlight the importance of the ice mass fraction in distinguishing the mixed-phase clouds and ice clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Multiscale ACI Satellite Database

The SATELLITE_EAGLES_PNNL NetCDF dataset contains a suite of satellite- and reanalysis-derived atmospheric and surface parameters on a regular latitude–longitude grid. The dataset includes core geophysical fields such as land fraction, aerosol optical depth at multiple wavelengths (465, 550, 667, and 865 nm), sea surface temperature, estimated inversion strength, and various thermodynamic and dynamic quantities (e.g., relative humidity, vertical velocity, boundary-layer height, and surface fluxes) from both MERRA and ERA reanalysis products, provided as daily-mean and instantaneous values. A major component of the dataset consists of MODIS-retrieved cloud microphysical properties, including cloud droplet number concentration, cloud effective radius, optical thickness, and liquid water path, provided for three compositing regimes (“All,” “Q06,” and “G18”). Corresponding cloud-top parameters—temperature, height, and pressure—along with total and domain-mean cloud fraction fields are also included. The file further integrates additional satellite data from AMSR-E (for cloud water, rain water, and surface precipitation retrievals) and CERES (for top-of-atmosphere radiative fluxes, cloud fractions, and albedo). This dataset is designed to evaluate aerosol–cloud interactions in warm clouds, emphasizing the use of MODIS for deriving cloud droplet number concentration and liquid water path statistics. The complementary satellite and reanalysis fields are co-located and time-matched to the same instantaneous MODIS observations, enabling consistent comparisons between cloud properties, aerosol loading, and large-scale meteorological conditions. The dataset is recently featured in Christensen et al. (2025), Machine Learning Reveals Strong Grid-Scale Dependence in the Satellite Nd–LWP Relationship, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2025-3850, 2025.

Christensen, Matthew [Pacific Northwest National L↗
DOE OSTI2025

Multiscale ACI Satellite Database

The SATELLITE_EAGLES_PNNL NetCDF dataset contains a suite of satellite- and reanalysis-derived atmospheric and surface parameters on a regular latitude–longitude grid. The dataset includes core geophysical fields such as land fraction, aerosol optical depth at multiple wavelengths (465, 550, 667, and 865 nm), sea surface temperature, estimated inversion strength, and various thermodynamic and dynamic quantities (e.g., relative humidity, vertical velocity, boundary-layer height, and surface fluxes) from both MERRA and ERA reanalysis products, provided as daily-mean and instantaneous values. A major component of the dataset consists of MODIS-retrieved cloud microphysical properties, including cloud droplet number concentration, cloud effective radius, optical thickness, and liquid water path, provided for three compositing regimes (“All,” “Q06,” and “G18”). Corresponding cloud-top parameters—temperature, height, and pressure—along with total and domain-mean cloud fraction fields are also included. The file further integrates additional satellite data from AMSR-E (for cloud water, rain water, and surface precipitation retrievals) and CERES (for top-of-atmosphere radiative fluxes, cloud fractions, and albedo). This dataset is designed to evaluate aerosol–cloud interactions in warm clouds, emphasizing the use of MODIS for deriving cloud droplet number concentration and liquid water path statistics. The complementary satellite and reanalysis fields are co-located and time-matched to the same instantaneous MODIS observations, enabling consistent comparisons between cloud properties, aerosol loading, and large-scale meteorological conditions. The dataset is recently featured in Christensen et al. (2025), Machine Learning Reveals Strong Grid-Scale Dependence in the Satellite Nd–LWP Relationship, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2025-3850, 2025.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Aerosol-induced closure of marine cloud cells: enhanced effects in the presence of precipitation

Abstract. The Weather Research Forecasting (WRF) version 4.3 model is configured within a Lagrangian framework to quantify the impact of aerosols on evolving cloud fields. Kilometer-scale simulations utilizing meteorological boundary conditions are based on 10 case study days offering diverse meteorology during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA). Measurements from aircraft, the ground-based Atmosphere Radiation Measurement (ARM) site at Graciosa Island in the Azores, and A-Train and geostationary satellites are utilized for validation, demonstrating good agreement with the WRF-simulated cloud and aerosol properties. Higher aerosol concentration leads to suppressed drizzle and increased cloud water content in all case study days. These changes lead to larger radiative cooling rates at cloud top, enhanced vertical velocity variance, and increased vertical and horizontal wind speed near the base of the lower-tropospheric inversion. As a result, marine cloud cell area expands, narrowing the gap between shallow clouds and increasing cloud optical thickness, liquid water content, and the top-of-atmosphere outgoing shortwave flux. While similar aerosol effects are observed in lightly to non-raining clouds, they tend to be smaller by comparison. These simulations show a relationship between cloud cell area expansion and the radiative adjustments caused by liquid water path and cloud fraction changes. The adjustments are positive and scale as 74 % and 51 %, respectively, relative to the Twomey effect. While higher-resolution large-eddy simulations may provide improved representation of cloud-top mixing processes, these results emphasize the importance of addressing mesoscale cloud-state transitions in the quantification of aerosol radiative forcing that cannot be attained from traditional climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Constraining the Twomey effect from satellite observations: issues and perspectives

Abstract. The Twomey effect describes the radiative forcing associated with a change in cloud albedo due to an increase in anthropogenic aerosol emissions. It is driven by the perturbation in cloud droplet number concentration (ΔNd, ant) in liquid-water clouds and is currently understood to exert a cooling effect on climate. The Twomey effect is the key driver in the effective radiative forcing due to aerosol–cloud interactions, but rapid adjustments also contribute. These adjustments are essentially the responses of cloud fraction and liquid water path to ΔNd, ant and thus scale approximately with it. While the fundamental physics of the influence of added aerosol particles on the droplet concentration (Nd) is well described by established theory at the particle scale (micrometres), how this relationship is expressed at the large-scale (hundreds of kilometres) perturbation, ΔNd, ant, remains uncertain. The discrepancy between process understanding at particle scale and insufficient quantification at the climate-relevant large scale is caused by co-variability of aerosol particles and updraught velocity and by droplet sink processes. These operate at scales on the order of tens of metres at which only localised observations are available and at which no approach yet exists to quantify the anthropogenic perturbation. Different atmospheric models suggest diverse magnitudes of the Twomey effect even when applying the same anthropogenic aerosol emission perturbation. Thus, observational data are needed to quantify and constrain the Twomey effect. At the global scale, this means satellite data. There are four key uncertainties in determining ΔNd, ant, namely the quantification of (i) the cloud-active aerosol – the cloud condensation nuclei (CCN) concentrations at or above cloud base, (ii) Nd, (iii) the statistical approach for inferring the sensitivity of Nd to aerosol particles from the satellite data and (iv) uncertainty in the anthropogenic perturbation to CCN concentrations, which is not easily accessible from observational data. This review discusses deficiencies of current approaches for the different aspects of the problem and proposes several ways forward: in terms of CCN, retrievals of optical quantities such as aerosol optical depth suffer from a lack of vertical resolution, size and hygroscopicity information, non-direct relation to the concentration of aerosols, difficulty to quantify it within or below clouds, and the problem of insufficient sensitivity at low concentrations, in addition to retrieval errors. A future path forward can include utilising co-located polarimeter and lidar instruments, ideally including high-spectral-resolution lidar capability at two wavelengths to maximise vertically resolved size distribution information content. In terms of Nd, a key problem is the lack of operational retrievals of this quantity and the inaccuracy of the retrieval especially in broken-cloud regimes. As for the Nd-to-CCN sensitivity, key issues are the updraught distributions and the role of Nd sink processes, for which empirical assessments for specific cloud regimes are currently the best solutions. These considerations point to the conclusion that past studies using existing approaches have likely underestimated the true sensitivity and, thus, the radiative forcing due to the Twomey effect.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Highly supercooled riming and unusual triple-frequency radar signatures over McMurdo Station, Antarctica

Abstract. Riming of ice crystals by supercooled water droplets is an efficient ice growth process, but its basic properties are still poorly known. While it has been shown to contribute significantly to surface precipitation at mid-latitudes, little is known about its occurrence at high latitudes. In Antarctica, two competing effects can influence the occurrence of riming: (i) the scarcity of supercooled liquid water clouds due to the extremely low tropospheric temperatures and (ii) the low aerosol concentration, which may lead to the formation of fewer and larger supercooled drops potentially resulting in an enhanced riming efficiency. In this work, by exploiting the deployment of an unprecedented number of multiwavelength remote sensing systems (including triple-frequency radar measurements) in West Antarctica, during the Atmospheric Radiation Measurements West Antarctic Radiation Experiment (AWARE) field campaign, we evaluate the riming incidence at McMurdo Station and find that riming occurs at much lower temperatures when compared to previous results in the mid-latitudes. This suggests the possible occurrence of a common atmospheric state over Antarctica that includes a rather stable atmosphere inhibiting turbulent mixing, and a high riming efficiency driven by large cloud droplets. We then focus on a peculiar case study featuring a persistent layer with a particularly pronounced riming signature in triple-frequency radar data but only a relatively modest amount of supercooled liquid water. In-depth analysis of the radar observations suggests that such signatures can only be explained by the combined effects of moderately rimed aggregates or similarly shaped florid polycrystals and a narrow particle size distribution (PSD). Simulations of this case study performed with a one-dimensional bin model indicate that similar triple-frequency radar observations can be reproduced when narrow PSDs are simulated. Such narrow PSDs can in turn be explained by two key factors: (i) the presence of a shallow homogeneous droplet or humidified aerosol freezing layer aloft seeding an underlying supercooled liquid layer, and (ii) the absence of turbulent mixing throughout a stable polar atmosphere that sustains narrow PSDs, as hydrometeors grow from the nucleation region aloft to ice particles of several millimeters in size, by vapor deposition and then riming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Visualization and analysis of coupling between plasmas self-organization and plasma-induced fluid circulation in 1 atm DC glows with liquid anode

The physical processes prevailing when plasma contacts liquid water are poorly understood; however, it is this very interaction that is the basis of new, exciting technologies for water purification, the treatment of cancer and disease and the production of new, high value chemical products. A chief impediment to widespread application of plasmas for the aforementioned applications is scale up and cost. Here we investigate plasma induced flows which facilitate transport of reactivity from the plasma to liquid water at rates much faster than diffusion. We find that this plasma induced motion varies depending on whether or not the plasma is self organized. The insight from this work paves the way to a better understanding of transport and ultimately informing plasma treatment technologies on how to best optimize the rate at which processes such as plasma based water purification occurs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Transferable Water Potentials Using Equivariant Neural Networks

Machine learning interatomic potentials (MLIPs) have emerged as a technique that promises quantum theory accuracy for reduced cost. It has been proposed [J. Chem. Phys. 2023, 158, 084111] that MLIPs trained on solely liquid water data cannot accurately transfer to the vapor–liquid equilibrium while recovering the many-body decomposition (MBD) analysis of gas-phase water clusters. This suggests that MLIPs do not directly learn the physically correct interactions of water molecules, limiting transferability. In this work, we show that MLIPs using equivariant architecture and trained on 3200 liquid water structures reproduces liquid-phase water properties (e.g., density within 0.003 g/cm 3 between 230 and 365 K), vapor–liquid equilibrium properties up to 550 K, the MBD analysis of gas-phase water cluster up to six-body interactions, and the relative energy and the vibrational density of states of ice phases. We show that potentials developed using equivariant MLIPs allow transferability for arbitrary phases of water that remain stable in nanosecond long simulations.

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

Large-Eddy Simulations of Marine Boundary Layer Clouds Associated with Cold-Air Outbreaks during the ACTIVATE Campaign. Part II: Aerosol–Meteorology–Cloud Interaction

Aerosol effects on micro/macrophysical properties of marine stratocumulus clouds over the western North Atlantic Ocean (WNAO) are investigated using in situ measurements and large-eddy simulations (LES) for two cold-air outbreak (CAO) cases (28 February and 1 March 2020) during the Aerosol Cloud Meteorology Interactions over the Western Atlantic Experiment (ACTIVATE). The LES is able to reproduce the vertical profiles of liquid water content (LWC), effective radius r eff and cloud droplet number concentration N c from fast cloud droplet probe (FCDP) in situ measurements for both cases. Furthermore, we show that aerosols affect cloud properties (N c , r eff , and LWC) via the prescribed bulk hygroscopicity of aerosols ($\bar{κ}$) and aerosol size distribution characteristics. N c , r eff , and liquid water path (LWP) are positively correlated to $\bar{κ}$ and aerosol number concentration (N a ) while cloud fractional cover (CFC) is insensitive to $\bar{κ}$ and aerosol size distributions for the two cases. The realistic changes to aerosol size distribution (number concentration, width, and the geometrical diameter) with the same meteorology state allow us to investigate aerosol effects on cloud properties without meteorological feedback. We also use the LES results to evaluate cloud properties from two reanalysis products, ERA5 and MERRA-2. Compared to LES, the ERA5 is able to capture the time evolution of LWP and total cloud coverage within the study domain during both CAO cases while MERRA-2 underestimates them.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Ultra-clean and smoky marine boundary layers frequently occur in the same season over the southeast Atlantic

We study 41 d with daily median surface accumulation mode aerosol particle concentrations below 50 cm -3 (ultra-clean conditions) observed at Ascension Island (ASI; 7.9° S, 14.4° W) between June 2016 and October 2017 as part of the Layered Atlantic Smoke Interactions with Clouds (LASIC) campaign. Interestingly, these days occur during a period of great relevance for aerosol–cloud–radiation interactions, the southeast Atlantic (SEATL)biomass-burning season (approximately June–October). That means that these critical months can feature both the highest surface aerosol numbers, from smoke intrusion into the marine boundary layer, as well as the lowest. While carbon monoxide and refractory black carbon concentrations on ultra-clean days do not approach those on days with heavy smoke, they also frequently exceed background concentrations calculated in the non-burning season from December 2016 to April 2017. This is evidence that even what become ultra-clean boundary layers can make contact with and entrain from an overlying SEATL smoke layer before undergoing a process of rapid aerosol removal. Because many ultra-clean and polluted boundary layers observed at Ascension Island during the biomass burning season follow similar isobaric back trajectories, the variability in this entrainment is likely more closely tied to the variability in the overlying smoke rather than large-scale horizontal circulation through the boundary layer. Since exceptionally low accumulation mode aerosol numbers at ASI do not necessarily indicate the relative lack of other trace pollutants, this suggests the importance of regional variations in what constitutes an “ultra-clean” marine boundary layer. Finally, surface drizzle rates, frequencies and accumulation – as well as retrievals of liquid water path– all consistently tend toward higher values on ultra-clean days. This implicates enhanced coalescence scavenging in low clouds as the key driver of ultra-clean events in the southeast Atlantic marine boundary layer. These enhancements occur against and are likely mediated by the backdrop of a seasonal increase in daily mean cloud fraction and daily median liquid water path over ASI, peaking in September and October in both LASIC years. Therefore the seasonality in ultra-clean day occurrence seems directly linked to the seasonality in SEATL cloud properties. These results highlight the importance of two-way aerosol–cloud interactions in the region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Identifying synoptic controls on boundary layer thermodynamic and cloud properties in a regional forecast model

Although most of our understanding of boundary layer cloudiness is based on idealized, subtropical, barotropic marine environments, boundary layer clouds exist across a range of conditions. In this study, we use the Naval Research Laboratory's Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS) and an automated cold-front-relative analysis framework to explore the boundary layer structure associated with low clouds across a transect through the cold front of a midlatitude synoptic cyclone. The model credibly captures boundary layer structure in line with conceptual models of midlatitude cyclones and ground-based observations at Graciosa Island in the Azores. The warm sector is conditionally unstable, with clouds that are too shallow and with too little liquid water, compared to cloud property retrievals from satellite and surface-based instruments. The cold-frontal region exhibits convection associated with weak stability and ascent. Northwest of the cold front, the boundary layer is well mixed, deeper, and capped by a strong inversion maintained by large-scale subsidence. Simulated clouds in frontal and post-frontal regions are mostly too thick, with too much liquid water and too little cloud base drizzle. The post-frontal clouds are associated with grid-scale updrafts, which appear to be the model's attempt to represent mesoscale organization of cellular convection typically observed in the cold sector of midlatitude cyclones. The deep, well-mixed post-frontal boundary layers and cloudiness are maintained by strong surface fluxes, as in cold air outbreaks. Our analysis framework serves as a unique approach to model verification, and our results offer insights into differences in boundary layer cloud behavior between subtropical and synoptic cold-sector regimes.

Eissner, Jordan M. [Univ. of Kansas, Lawrence, KS ↗
DOE OSTI2026

Substantial Diel Changes of Cloud Adjustments to Aerosols in Ship‐Tracks

Clouds adjust their albedo, amount and liquid water path to anthropogenic aerosols, and produce a climate forcing. Such cloud adjustments are complex functions of both environmental conditions and interaction time. However, direct observations of temporal changes in cloud adjustments are limited. Here, we develop a method to show observational evidence of substantial diel changes in cloud adjustments within ship-tracks. Polluted clouds experiencing extended nighttime exposure exhibit higher increases in cloud fraction while liquid-water path adjustments follow a non-monotonic evolution. The higher cloud fraction increase persists through daytime. Such diel variations likely stem from diel changes in solar heating and background cloud precipitation, and are important because nighttime changes persist and affect solar radiation during the day. Accounting for diel variations produces a stronger aerosol forcing, by −0.1 ∼−0.4 W⁢m −2 . Our approach reveals how the temporal evolution of cloud adjustments can have important implications for aerosol-cloud interactions.

Yuan, Tianle [NASA Goddard Space Flight Center (GS↗
DOE OSTI2024

Impacts of Synoptic-Scale Dynamics on Clouds and Radiation in High Southern Latitudes

High-latitudinal mixed-phase clouds significantly affect Earth's radiative balance. Observations of cloud and radiative properties from two field campaigns in the Southern Ocean and Antarctica were compared with two global climate model simulations. A cyclone compositing method was used to quantify “dynamics-cloud-radiation” relationships relative to the extratropical cyclone centers. Observations show larger asymmetry in cloud and radiative properties between western and eastern sectors at McMurdo compared with Macquarie Island. Most observed quantities at McMurdo are higher in the western (i.e., post-frontal) than the eastern (frontal) sector, including cloud fraction, liquid water path (LWP), net surface shortwave and longwave radiation (SW and LW), except for ice water path (IWP) being higher in the eastern sector. Here, the two models were found to overestimate cloud fraction and LWP at Macquarie Island but underestimate them at McMurdo Station. IWP is consistently underestimated at both locations, both sectors, and in all seasons. Biases of cloud fraction, LWP, and IWP are negatively correlated with SW biases and positively correlated with LW biases. The persistent negative IWP biases may have become one of the leading causes of radiative biases over the high southern latitudes, after correcting the underestimation of supercooled liquid water in the older model versions. By examining multi-scale factors from cloud microphysics to synoptic dynamics, this work will help increase the fidelity of climate simulations in this remote region.

58 GEOSCIENCES↗
DOE OSTI2025

The Snow Albedo Evolution (SALVO) Campaign at the North Slope of Alaska Field Campaign Report

The springtime surface-albedo transition in the Alaskan Arctic, along with the forces that determine its duration and nature, was the focus of the Snow ALbedo eVOlution (SALVO) campaign. The SALVO team investigated the reasons and durations of the stages of spring melt, during which albedo values decrease from 0.8 to 0.1, signifying the year’s largest and most significant radiative energy change. SALVO II (2022, 2024) built on findings from the successful SALVO I (2019) melt season campaign and previous research on the melt transition period conducted by Grenfell and Perovich (2004). SALVO used the ARM NSA observatory, thus aligning with a core principle of the ARM Decadal Vision to “provide comprehensive and impactful field measurements to support scientific advancement of atmospheric process understanding.” SALVO fieldwork was conducted in the spring near the ARM NSA Central Facility in Utqiagvik, Alaska (Figure 1). Comprehensive data sets were collected in 2019, 2022, and 2024, which included spectral and broadband albedos, snow and meltwater depths, snow stratigraphy, snow grain size data, and aerial imagery. In April of each project year, before the snow began to melt, we established survey lines at three or four locations: inland tundra at NSA E12, coastal tundra at the NSA C1 site, Elson Lagoon, and offshore on the Chukchi Sea sea ice (2022 only). We set up a 200-meter-long line, marked every five meters, to enable repeated measurements of surface albedo and snow depth at each location. During most site visits, the SALVO II team dug at least one snow pit to assess the characteristics of the vertical snow layers, including the types and sizes of snow grains, snow density, and the distribution of liquid water. We evaluated the characteristics of flowing or pooled water at the bottom of the snowpack. Watching snow melt in the Arctic has been compared by some to being as exciting as watching paint dry, but it is far more exciting and dynamic than that. Initially, changes happen slowly while the snow cover remains above 80%, but then reductions in the albedo of the landscape, from about 0.8 to 0.4, combined with rising spring temperatures, begin to accelerate the melt. It happens so quickly that no matter how hard the field team works, they cannot keep up with documenting the changes. At first, liquid water is scarce—only a few wet layers of snow in the snowpack. But then, water becomes ubiquitous, transporting melt energy and creating an extremely heterogeneous albedo and melt landscape. At the end of each field season, the SALVO team is exhausted, relieved to see the snow gone, and joyful to play in the melt ponds on the sea ice under the midnight sun. Overall, the 2024 SALVO field campaign was a highlight, reflecting numerous lessons learned from 2019 and 2022. Daily measurements were taken more than 10 times at each site, resulting in a comprehensive time series of snow conditions and albedo. The team deployed using snowmachines and then on foot when snowmachine travel was no longer permitted on the melting tundra. Instrument mounts and packaging were optimized for quick deployment, regardless of the mode of transportation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Improving the reliability of machine learned potentials for modeling inhomogeneous liquids

The atomic-scale response of inhomogeneous fluids at interfaces and surrounding solute particles plays a critical role in governing chemical, electrochemical, and biological processes. Classical molecular dynamics simulations have been applied extensively to simulate the response of fluids to inhomogeneities directly, but are limited by the accuracy of the underlying interatomic potentials. Here, we use neural network potentials (NNPs) trained to ab initio simulations to accurately predict the inhomogeneous responses of two distinct fluids: liquid water and molten NaCl. Although NNPs can be readily trained to model complex bulk systems across a range of state points, we show that to appropriately model a fluid's response at an interface, relevant inhomogeneous configurations must be included in the training data. In order to sufficiently sample appropriate configurations of such inhomogeneous fluids, we develop protocols based on molecular dynamics simulations in the presence of external potentials. We demonstrate that NNPs trained on inhomogeneous fluid configurations can more accurately predict several key properties of fluids—including the density response, surface tension and size-dependent cavitation free energies—for liquid water and molten NaCl, compared to both empirical interatomic potentials and NNPs that are not trained on such inhomogeneous configurations. This work therefore provides a first demonstration and framework to extract the response of inhomogeneous fluids from first principles for classical density-functional treatment of fluids free from empirical potentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗