Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geophysical”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Numerical Modeling of Anisotropic Elastic-Wave Sensitivity Propagation for Optimal Design of Time-Lapse Seismic Surveys

Reliable subsurface time-lapse seismic monitoring is crucial for many geophysical applications, such as enhanced geothermal system characterization, geologic carbon utilization and storage, and conventional and unconventional oil/gas reservoir characterization, etc. In this work, we develop an elastic-wave sensitivity propagation method for optimal design of cost-effective time-lapse seismic surveys considering the fact that most of subsurface geologic layers and fractured reservoirs are anisotropic instead of isotropic. For anisotropic media, we define monitoring criteria using qP- and qS-wave sensitivity energies after decomposing qP- and qS-wave components from the total elastic-wave sensitivity wavefield using a hybrid time- and frequency-domain approach. Geophones should therefore be placed at locations with significant qP- and qS-wave sensitivity energies for cost-effective time-lapse seismic monitoring in an anisotropic geology setting. Our numerical modeling results for a modified anisotropic Hess model demonstrate that, compared with the isotropic case, subsurface anisotropy changes the spatial distributions of elastic-wave sensitivity energies. Consequently, it is necessary to consider subsurface anisotropies when designing the spatial distribution of geophones for cost-effective time-lapse seismic monitoring. This finding suggests that it is essential to use our new anisotropic elastic-wave sensitivity modeling method for optimal design of time-lapse seismic surveys to reliably monitor the changes in subsurface reservoirs, fracture zones or target monitoring regions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Climate model projections from the Scenario Model Intercomparison Project (ScenarioMIP) of CMIP6

Abstract. The Scenario Model Intercomparison Project (ScenarioMIP) defines and coordinates the main set of future climate projections, based on concentration-driven simulations, within the Coupled Model Intercomparison Project phase 6 (CMIP6). This paper presents a range of its outcomes by synthesizing results from the participating global coupled Earth system models. We limit our scope to the analysis of strictly geophysical outcomes: mainly global averages and spatial patterns of change for surface air temperature and precipitation. We also compare CMIP6 projections to CMIP5 results, especially for those scenarios that were designed to provide continuity across the CMIP phases, at the same time highlighting important differences in forcing composition, as well as in results. The range of future temperature and precipitation changes by the end of the century (2081–2100) encompassing the Tier 1 experiments based on the Shared Socioeconomic Pathway (SSP) scenarios (SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5) and SSP1-1.9 spans a larger range of outcomes compared to CMIP5, due to higher warming (by close to 1.5 ∘C) reached at the upper end of the 5 %–95 % envelope of the highest scenario (SSP5-8.5). This is due to both the wider range of radiative forcing that the new scenarios cover and the higher climate sensitivities in some of the new models compared to their CMIP5 predecessors. Spatial patterns of change for temperature and precipitation averaged over models and scenarios have familiar features, and an analysis of their variations confirms model structural differences to be the dominant source of uncertainty. Models also differ with respect to the size and evolution of internal variability as measured by individual models' initial condition ensemble spreads, according to a set of initial condition ensemble simulations available under SSP3-7.0. These experiments suggest a tendency for internal variability to decrease along the course of the century in this scenario, a result that will benefit from further analysis over a larger set of models. Benefits of mitigation, all else being equal in terms of societal drivers, appear clearly when comparing scenarios developed under the same SSP but to which different degrees of mitigation have been applied. It is also found that a mild overshoot in temperature of a few decades around mid-century, as represented in SSP5-3.4OS, does not affect the end outcome of temperature and precipitation changes by 2100, which return to the same levels as those reached by the gradually increasing SSP4-3.4 (not erasing the possibility, however, that other aspects of the system may not be as easily reversible). Central estimates of the time at which the ensemble means of the different scenarios reach a given warming level might be biased by the inclusion of models that have shown faster warming in the historical period than the observed. Those estimates show all scenarios reaching 1.5 ∘C of warming compared to the 1850–1900 baseline in the second half of the current decade, with the time span between slow and fast warming covering between 20 and 27 years from present. The warming level of 2 ∘C of warming is reached as early as 2039 by the ensemble mean under SSP5-8.5 but as late as the mid-2060s under SSP1-2.6. The highest warming level considered (5 ∘C) is reached by the ensemble mean only under SSP5-8.5 and not until the mid-2090s.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

PAMTRA 1.0: the Passive and Active Microwave radiative TRAnsfer tool for simulating radiometer and radar measurements of the cloudy atmosphere

Forward models are a key tool to generate synthetic observations given knowledge of the atmospheric state. In this way, they are an integral part of inversion algorithms that aim to retrieve geophysical variables from observations or in data assimilation. Their application for the exploitation of the full information content of remote sensing observations becomes increasingly important when these are used to evaluate the performance of cloud-resolving models (CRMs). Herein, CRM profiles or fields provide the input to the forward model whose simulation results are subsequently compared to the observations. This paper introduces the freely available comprehensive microwave forward model PAMTRA (Passive and Active Microwave TRAnsfer), demonstrates its capabilities to simulate passive and active measurements across the microwave spectral region for upward- and downward-looking geometries, and illustrates how the forward simulations can be used to evaluate CRMs and to interpret measurements to improve our understanding of cloud processes. PAMTRA is unique as it treats passive and active radiative transfer (RT) in a consistent way with the passive forward model providing upwelling and downwelling polarized brightness temperatures and radiances for arbitrary observation angles. The active part is capable of simulating the full radar Doppler spectrum and its moments. PAMTRA is designed to be flexible with respect to instrument specifications and interfaces to many different formats of input and output, especially CRMs, spanning the range from bin-resolved microphysical output to one- and two-moment schemes, and to in situ measured hydrometeor properties. A specific highlight is the incorporation of the self-similar Rayleigh–Gans approximation (SSRGA) for both active and passive applications, which becomes especially important for the investigation of frozen hydrometeors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Testing the reliability of interpretable neural networks in geoscience using the Madden–Julian oscillation

Abstract. We test the reliability of two neural network interpretation techniques, backward optimization and layerwise relevance propagation, within geoscientific applications by applying them to a commonly studied geophysical phenomenon, the Madden–Julian oscillation. The Madden–Julian oscillation is a multi-scale pattern within the tropical atmosphere that has been extensively studied over the past decades, which makes it an ideal test case to ensure the interpretability methods can recover the current state of knowledge regarding its spatial structure. The neural networks can, indeed, reproduce the current state of knowledge and can also provide new insights into the seasonality of the Madden–Julian oscillation and its relationships with atmospheric state variables. The neural network identifies the phase of the Madden–Julian oscillation twice as accurately as a linear regression approach, which means that nonlinearities used by the neural network are important to the structure of the Madden–Julian oscillation. Interpretations of the neural network show that it accurately captures the spatial structures of the Madden–Julian oscillation, suggest that the nonlinearities of the Madden–Julian oscillation are manifested through the uniqueness of each event, and offer physically meaningful insights into its relationship with atmospheric state variables. We also use the interpretations to identify the seasonality of the Madden–Julian oscillation and find that the conventionally defined extended seasons should be shifted later by 1 month. More generally, this study suggests that neural networks can be reliably interpreted for geoscientific applications and may thereby serve as a dependable method for testing geoscientific hypotheses.

58 GEOSCIENCES↗
DOE OSTI2022

Efficient high-dimensional variational data assimilation with machine-learned reduced-order models

Abstract. Data assimilation (DA) in geophysical sciences remains the cornerstone of robust forecasts from numerical models. Indeed, DA plays a crucial role in the quality of numerical weather prediction and is a crucial building block that has allowed dramatic improvements in weather forecasting over the past few decades. DA is commonly framed in a variational setting, where one solves an optimization problem within a Bayesian formulation using raw model forecasts as a prior and observations as likelihood. This leads to a DA objective function that needs to be minimized, where the decision variables are the initial conditions specified to the model. In traditional DA, the forward model is numerically and computationally expensive. Here we replace the forward model with a low-dimensional, data-driven, and differentiable emulator. Consequently, gradients of our DA objective function with respect to the decision variables are obtained rapidly via automatic differentiation. We demonstrate our approach by performing an emulator-assisted DA forecast of geopotential height. Our results indicate that emulator-assisted DA is faster than traditional equation-based DA forecasts by 4 orders of magnitude, allowing computations to be performed on a workstation rather than a dedicated high-performance computer. In addition, we describe accuracy benefits of emulator-assisted DA when compared to simply using the emulator for forecasting (i.e., without DA). Our overall formulation is denoted AIEADA (Artificial Intelligence Emulator-Assisted Data Assimilation).

58 GEOSCIENCES↗
DOE OSTI2022

Lossy checkpoint compression in full waveform inversion: a case study with ZFPv0.5.5 and the overthrust model

This paper proposes a new method that combines checkpointing methods with error-controlled lossy compression for large-scale high-performance full-waveform inversion (FWI), an inverse problem commonly used in geophysical exploration. This combination can significantly reduce data movement, allowing a reduction in run time as well as peak memory. In the exascale computing era, frequent data transfer (e.g., memory bandwidth, PCIe bandwidth for GPUs, or network) is the performance bottleneck rather than the peak FLOPS of the processing unit. Like many other adjoint-based optimization problems, FWI is costly in terms of the number of floating-point operations, large memory footprint during backpropagation, and data transfer overheads. Past work for adjoint methods has developed checkpointing methods that reduce the peak memory requirements during backpropagation at the cost of additional floating-point computations. Combining this traditional checkpointing with error-controlled lossy compression, we explore the three-way tradeoff between memory, precision, and time to solution. We investigate how approximation errors introduced by lossy compression of the forward solution impact the objective function gradient and final inverted solution. Empirical results from these numerical experiments indicate that high lossy-compression rates (compression factors ranging up to 100) have a relatively minor impact on convergence rates and the quality of the final solution.

58 GEOSCIENCES↗
DOE OSTI2023

Massively parallel modeling and inversion of electrical resistivity tomography data using PFLOTRAN

Abstract. Electrical resistivity tomography (ERT) is a broadly accepted geophysical method for subsurface investigations. Interpretation of field ERT data usually requires the application of computationally intensive forward modeling and inversion algorithms. For large-scale ERT data, the efficiency of these algorithms depends on the robustness, accuracy, and scalability on high-performance computing resources. In this regard, we present a robust and highly scalable implementation of forward modeling and inversion algorithms for ERT data. The implementation is publicly available and developed within the framework of PFLOTRAN, an open-source, state-of-the-art massively parallel subsurface flow and transport simulation code. The forward modeling is based on a finite-volume discretization of the governing differential equations, and the inversion uses a Gauss–Newton optimization scheme. To evaluate the accuracy of the forward modeling, two examples are first presented by considering layered (1D) and 3D earth conductivity models. The computed numerical results show good agreement with the analytical solutions for the layered earth model and results from a well-established code for the 3D model. Inversion of ERT data, simulated for a 3D model, is then performed to demonstrate the inversion capability by recovering the conductivity of the model. To demonstrate the parallel performance of PFLOTRAN's ERT process model and inversion capabilities, large-scale scalability tests are performed by using up to 131 072 processes on a leadership class supercomputer. These tests are performed for the two most computationally intensive steps of the ERT inversion: forward modeling and Jacobian computation. For the forward modeling, we consider models with up to 122 ×106 degrees of freedom (DOFs) in the resulting system of linear equations and demonstrate that the code exhibits almost linear scalability on up to 10 000 DOFs per process. On the other hand, the code shows superlinear scalability for the Jacobian computation, mainly because all computations are fairly evenly distributed over each process with no parallel communication.

58 GEOSCIENCES↗
DOE OSTI2024

Identifying atmospheric rivers and their poleward latent heat transport with generalizable neural networks: ARCNNv1

Abstract. Atmospheric rivers (ARs) are extreme weather events that can alleviate drought or cause billions of US dollars in flood damage. By transporting significant amounts of latent energy towards the poles, they are crucial to maintaining the climate system's energy balance. Since there is no first-principle definition of an AR grounded in geophysical fluid mechanics, AR identification is currently performed by a multitude of expert-defined, threshold-based algorithms. The variety of AR detection algorithms has introduced uncertainty into the study of ARs, and the thresholds of the algorithms may not generalize to new climate datasets and resolutions. We train convolutional neural networks (CNNs) to detect ARs while representing this uncertainty; we name these models ARCNNs. To detect ARs without requiring new labeled data and labor-intensive AR detection campaigns, we present a semi-supervised learning framework based on image style transfer. This framework generalizes ARCNNs across climate datasets and input fields. Using idealized and realistic numerical models, together with observations, we assess the performance of the ARCNNs. We test the ARCNNs in an idealized simulation of a shallow-water fluid in which nearly all the tracer transport can be attributed to AR-like filamentary structures. In reanalysis and a high-resolution climate model, we use ARCNNs to calculate the contribution of ARs to meridional latent heat transport, and we demonstrate that this quantity varies considerably due to AR detection uncertainty.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Merged Observatory Data Files (MODFs): an integrated observational data product supporting process-oriented investigations and diagnostics

A large and ever-growing body of geophysical information is measured in campaigns and at specialized observatories as a part of scientific expeditions and experiments. These collections of observed data include many essential climate variables (as defined by the Global Climate Observing System) but are often distinguished by a wide range of additional non-routine measurements that are designed to not only document the state of the environment but also the drivers that contribute to that state. These field data are used not only to further understand environmental processes through observation-based studies but also to provide baseline data to test model performance and to codify understanding to improve predictive capabilities. To address the considerable barriers and difficulty in utilizing these diverse and complex data for observation–model research, the Merged Observatory Data File (MODF) concept has been developed. A MODF combines measurements from multiple instruments into a single file that complies with well-established data format and metadata practices and has been designed to parallel the development of corresponding Merged Model Data Files (MMDFs). Using the MODF and MMDF protocols will facilitate the evolution of model intercomparison projects into model intercomparison and improvement projects by putting observation and model data “on the same page” in a timely manner. The MODF concept was developed especially for weather forecast model studies in the Arctic. The surprisingly complex process of implementing MODFs in that context refined the concept itself. Thus, this article explains the concept of MODFs by providing details on the issues that were revealed and resolved during that first specific implementation. Detailed instructions are provided on how to make MODFs, and this article can be considered a MODF creation manual.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A strainmeter array as the fulcrum of novel observatory sites along the Alto Tiberina Near Fault Observatory

Fault slip is a complex natural phenomenon involving multiple spatiotemporal scales from seconds to days to weeks. To understand the physical and chemical processes responsible for the full fault slip spectrum, a multidisciplinary approach is highly recommended. The Near Fault Observatories (NFOs) aim at providing high-precision and spatiotemporally dense multidisciplinary near-fault data, enabling the generation of new original observations and innovative scientific products. The Alto Tiberina Near Fault Observatory is a permanent monitoring infrastructure established around the Alto Tiberina fault (ATF), a 60 km long low-angle normal fault (mean dip 20°), located along a sector of the Northern Apennines (central Italy) undergoing an extension at a rate of about 3 mm yr –1 . The presence of repeating earthquakes on the ATF and a steep gradient in crustal velocities measured across the ATF by GNSS stations suggest large and deep (5–12 km) portions of the ATF undergoing aseismic creep. Both laboratory and theoretical studies indicate that any given patch of a fault can creep, nucleate slow earthquakes, and host large earthquakes, as also documented in nature for certain ruptures (e.g., Iquique in 2014, Tōhoku in 2011, and Parkfield in 2004). Nonetheless, how a fault patch switches from one mode of slip to another, as well as the interaction between creep, slow slip, and regular earthquakes, is still poorly documented by near-field observation. With the strainmeter array along the Alto Tiberina fault system (STAR) project, we build a series of six geophysical observatory sites consisting of 80–160 m deep vertical boreholes instrumented with strainmeters and seismometers as well as meteorological and GNSS antennas and additional seismometers at the surface. By covering the portions of the ATF that exhibits repeated earthquakes at shallow depth (above 4 km) with these new observatory sites, we aim to collect unique open-access data to answer fundamental questions about the relationship between creep, slow slip, dynamic earthquake rupture, and tectonic faulting.

58 GEOSCIENCES↗
DOE OSTI2022

Resolving glacial isostatic adjustment (GIA) in response to modern and future ice loss at marine grounding lines in West Antarctica

Accurate glacial isostatic adjustment (GIA) modelling in the cryosphere is required for interpreting satellite, geophysical and geological records and for assessing the feedbacks of Earth deformation and sea-level change on marine ice-sheet grounding lines. GIA modelling in areas of active ice loss in West Antarctica is particularly challenging because the ice is underlain by laterally varying mantle viscosities that are up to several orders of magnitude lower than the global average, leading to a faster and more localised response of the solid Earth to ongoing and future ice-sheet retreat and necessitating GIA models that incorporate 3-D viscoelastic Earth structure. Improvements to GIA models allow for computation of the viscoelastic response of the Earth to surface ice loading at sub-kilometre resolution, and ice-sheet models and observational products now provide the inputs to GIA models at comparably unprecedented detail. However, the resolution required to accurately capture GIA in models remains poorly understood, and high-resolution calculations come at heavy computational expense. We adopt a 3-D GIA model with a range of Earth structure models based on recent seismic tomography and geodetic data to perform a comprehensive analysis of the influence of grid resolution on predictions of GIA in the Amundsen Sea Embayment (ASE) in West Antarctica. Through idealised sensitivity testing down to sub-kilometre resolution with spatially isolated ice loading changes, we find that a grid resolution of ~$\frac{1}{3}$ of the radius of the load or higher is required to accurately capture the elastic response of the Earth. However, when we consider more realistic, spatially coherent ice loss scenarios based on modern observational records and future ice-sheet model projections and adopt a viscoelastic Earth, we find that predicted deformation and sea-level change along the grounding line converge to within 5 % with grid resolutions of 7.5 km or higher, and to within 2% for grid resolutions of 3.75 km and higher, even when the input ice model is on a 1 km grid. Furthermore, we show that low mantle viscosities beneath the ASE lead to viscous deformation that contributes to the instrumental record on decadal timescales and equals or dominates over elastic effects by the end of the 21st century. Our findings suggest that for the range of resolutions of 1.9–15 km that we considered, the error due to adopting a coarser grid in this region is negligible compared to the effect of neglecting viscous effects and the uncertainty in the adopted mantle viscosity structure.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

US East Coast synthetic aperture radar wind atlas for offshore wind energy

We present the first synthetic aperture radar (SAR) offshore wind atlas of the US East Coast from Georgia to the Canadian border. Images from RADARSAT-1, Envisat, and Sentinel-1A/B are processed to wind maps using the geophysical model function (GMF) CMOD5.N. Extensive comparisons with 6008 collocated buoy observations of the wind speed reveal that biases of the individual systems range from -0.8 to 0.6 m s -1 . Unbiased wind retrievals are crucial for producing an accurate wind atlas, and intercalibration of the SAR observations is therefore applied. Wind retrievals from the intercalibrated SAR observations show biases in the range of to -0.2 to 0.0 m s -1 , while at the same time improving the root-mean-squared error from 1.67 to 1.46 m s -1 . The intercalibrated SAR observations are, for the first time, aggregated to create a wind atlas at the height 10 m a.s.l. (above sea level). The SAR wind atlas is used as a reference to study wind resources derived from the Wind Integration National Dataset Toolkit (WTK), which is based on 7 years of modelling output from the Weather Research and Forecasting (WRF) model. Comparisons focus on the spatial variation in wind resources and show that model outputs lead to lower coastal wind speed gradients than those derived from SAR. Areas designated for offshore wind development by the Bureau of Ocean Energy Management are investigated in more detail; the wind resources in terms of the mean wind speed show spatial variations within each designated area between 0.3 and 0.5 m s -1 for SAR and less than 0.2m s -1 for the WTK. Our findings indicate that wind speed gradients and variations might be underestimated in mesoscale model outputs along the US East Coast.

17 WIND ENERGY↗
DOE OSTI2025

First experimental determination of the ⁴⁰Ar(𝑛, 2𝑛)³⁹Ar reaction cross section and ³⁹Ar production in Earth’s atmosphere

The cosmogenic ³⁹Ar(t1∕2 = 268 years) isotope of argon is used for geophysical dating and tracing of underground and ocean water, as well as ice owing to its appropriate half-life and chemical inertness as a noble gas; ³⁹Ar serves also in nuclear weapon test monitoring. We measured for the first time the total cross section of the main ³⁹Ar cosmogenic production reaction in the atmosphere, namely ⁴⁰Ar(𝑛, 2𝑛)39 Ar, using 14.8 ± 0.3 MeV neutrons. The neutrons, produced by a deuterium-tritium generator, impinged on a stainless steel sphere filled with Ar gas highly enriched in the ⁴⁰Ar isotope and were monitored by a stack of fast-neutron activation foils. The reaction yield was measured by atom counting of long-lived ³⁹Ar with noble gas accelerator mass spectrometry and, independently, by decay counting relative to atmospheric argon (³⁹Ar/Ar= 8.12 × 10−16 ). A total ⁴⁰Ar(𝑛, 2𝑛)39 Ar cross section of 610 ± 100 mb was determined at 14.8 ± 0.3 MeV incident neutron energy. This result serves as a benchmark for recent theoretical calculations and evaluations, found to reproduce well the experimental total cross section. We use these energy-dependent theoretical cross sections together with experimental spectra of cosmogenic neutrons at different altitudes to calculate the global average rate of neutron-induced ³⁹Ar atmospheric production, resulting in 770 ± 240 39 Ar atoms/cm²/day. The secular equilibrium between the ³⁹Ar calculated production rate and radioactive decay rate leads to a partial isotopic abundance ³⁹Ar/Ar = (5.9 ± 1.8) × 10−16 , showing that ≈73% of atmospheric ³⁹Ar is produced by cosmogenic neutrons, the remaining part believed to be induced by muons and high-energy 𝛾 rays. The ⁴⁰Ar(𝑛, 2𝑛)³⁹Ar cross section at 14 MeV is also a key parameter for quantifying the anthropogenic contribution to atmospheric 39 Ar produced during the thermonuclear tests of the 1960s. We estimate that anthropogenic ³⁹Ar accounts for roughly 20% of the present atmospheric inventory.

39 Ar atmospheric production↗
DOE OSTI2021

CACTI ARM Mobile Facility (AMF) Measurements of Ice Nucleating Particles

The dataset comprises measures, using Colorado State University's Ice Spectrometer (IS, an immersion freezing device with a range from 0°C down to -26 to -29°C) of atmosperic ice nucleating particle (INP) concentrations taken on ARM's Mobile Facility (AMF-1) Aerosol Observing System, near Villa Yacanto in central Argentina, during the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) Experiment. A filter sampler was mounted on the AOS trailer (2 m below the AOS inlet). Single-use filter units open to the atmosphere fitted with pre-cleaned and pre-sterilized, 47-mm dia. Nuclepore polycarbonate filters were used. Filters were typically drawn for an 8-hour period, totaling 6000 liters sampled on average. After collection, filters were stored at -20°C until processed. Initial processing to obtain spectra of INP number concentration active via the immersion freezing mechanism versus temperature was conducted using CSU's IS instrument (McCluskey et al., 2018). For measurment of INPs, collected aerosol particles were re-suspended in 7-8 mL of 0.1 µm-filtered deionized (DI) water. Aliquots of each suspension, and serial dilutions, were dispensed into trays which were fit into aluminum blocks in the IS. Samples are cooled at 0.33°C min-1 and the freezing temperatures of wells recorded automatically. Cumulative INP concentrations were determined by first calculating the INPs per mL of suspension based on Vali (1971) and then converting to concentration per standard liter of air using the proportion of the total liquid sample dispensed and the air sample volumes. An aliquot of suspension fro selected samples were also heat treated (95°C for 20 min) to denature and deactivate biological INPs, and another aliquot digested in 10% H2O2 at 95°C under UV-B to remove all organic carbon INPs. McCluskey, C. S., J. Ovadnevaite, M. Rinaldi, J. Atkinson, F. Belosi, D. Ceburnis, S. Marullo, T. C. J. Hill, U. Lohmann, Z. A. Kanji, C. O’Dowd, S. M. Kreidenweis, P. J. DeMott, 2018: Marine and Terrestrial Organic Ice Nucleating Particles in Pristine Marine to Continentally-Influenced Northeast Atlantic Air Masses, Journal of Geophysical Research: Atmospheres, 123, 6196–6212, https://doi.org/10.1029/2017JD028033. Vali, G., 1971: Quantitative evaluation of experimental results on the heterogeneous freezing nucleation of supercooled liquids. J. Atmos. Sci., 28, 402–409.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

CACTI ARM Aerial Facility Measurements of Ice Nucleating Particles

The dataset comprises measures, using Colorado State University's Ice Spectrometer (IS, an immersion freezing device with a range from 0°C down to -26 to -29°C) of atmospheric ice nucleating particle (INP) concentrations taken on the Atmospheric Radiation Measurement (ARM) program Aerial Facilty (AAF) G-1 aircraft. INP measurements on the G-1 were collected from varied altitudes on different flights over the region of the Sierras de Córdoba mountain range of north-central Argentina, centred over ARM's Mobile Facility (AMF-1) near Villa Yacanto, where ground-based INP measures were being taken concurrently. Both studies took place as part of the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) Experiment. A filter sampling system was deployed on the G-1 to collect aerosol particles for post-processing of INPs collected on filters (measuring their immersion freezing ability) once returned to Colorado State University (CSU). Filter holders used were pre-cleaned aluminum in-line units loaded with pre-cleaned and pre-sterilized 47 mm dia. Nuclepore polycarbonate filters (0.2 mm pore size). Filters were drawn for varied times, resulting in varied volumes collected (51 to 1667 SL). Mass flow rate was recorded in real-time so that total sampled volume (at standard temperature and pressure) could be determined. A total of 34 sample filters were collected over the IOP, including 5 blanks. Filters were stored at -20°C freezer prior to frozen return to Colorado State University (CSU). . Processing to obtain spectra of INP number concentration active via the immersion freezing mechanism versus temperature was conducted using CSU's IS instrument (McCluskey et al., 2018). For measurment of INPs, collected aerosol particles were re-suspended in 7 mL of 0.02 µm-filtered deionized water. Aliquots of each suspension, and serial dilutions, were dispensed into trays which were fit into aluminum blocks in the IS. Samples are cooled at 0.33°C min-1 and the freezing temperatures of wells recorded automatically. Cumulative INP concentrations were determined by first calculating the INPs per mL of suspension based on Vali (1971) and then converting to concentration per standard liter of air using the proportion of the total liquid sample dispensed and the air sample volumes. Aliquots of suspensions from selected samples were also heat treated (95°C for 20 min) to denature and deactivate biological INPs, and digested in 10% H2O2 at 95°C under UV-B to remove all organic carbon INPs. McCluskey, C. S., J. Ovadnevaite, M. Rinaldi, J. Atkinson, F. Belosi, D. Ceburnis, S. Marullo, T. C. J. Hill, U. Lohmann, Z. A. Kanji, C. O’Dowd, S. M. Kreidenweis, P. J. DeMott, 2018: Marine and Terrestrial Organic Ice Nucleating Particles in Pristine Marine to Continentally-Influenced Northeast Atlantic Air Masses, Journal of Geophysical Research: Atmospheres, 123, 6196–6212, https://doi.org/10.1029/2017JD028033. Vali, G., 1971: Quantitative evaluation of experimental results on the heterogeneous freezing nucleation of supercooled liquids. J. Atmos. Sci., 28, 402–409.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

MOSAiC-Colorado State University Ice Spectrometer

This data set contains atmospheric ice nucleating particle (INP) measurements, using Colorado State University&rsquo;s (CSU) Ice Spectrometer (IS), of filter collections taken at the U.S. DOE ARM AMF2 site onboard the R/V Polarstern P-deck during the Multidisciplinary Drifting Observatory for the Study of Arctic Climate (MOSAiC) field campaign. Samples were collected from October 27, 2019 to September 24, 2020. A filter sampler was mounted approximately 15 m above ground level on a railing in proximity to (and approximately 3 m below) the Aerosol Observation System (AOS) inlet. Single-use filter units open to the atmosphere were pre-cleaned and pre-loaded with 47-mm diameter Nuclepore polycarbonate (0.2 &micro;m pore-diameter) filters. Filters were typically drawn for a three-day period, with an average volume of air filtered of 87,000 standard liters. Total volumes were calculated through recorded daily flow rates using a mass flow meter (TSI). After collection, filters were stored and transported frozen until analysis using CSU&rsquo;s IS instrument (McCluskey et al., 2018). Aerosol particles were first re-suspended in 8 mL of 0.1 &micro;m-filtered deionized (DI) water. Aliquots of each suspension, and corresponding 11-fold dilutions, were dispensed into polymerase chain reaction (PCR) trays and placed into the aluminum blocks of the IS. Samples were cooled at approximately 0.33 &deg;C min -1 and freezing detected optically with corresponding temperatures recorded. Cumulative INP concentrations were determined through calculating the number of INPs per mL of suspension (Vali, 1971) and converting to concentration per standard L of air by accounting for the proportion of liquid used and volume of air collected. All samples were corrected for the number of INPs on the average of four field blanks (cleaned, handled, transported, and analyzed in the same way without air flow). Two-tailed, 95% confidence intervals for binomial sampling are provided (Agresti and Coull, 1998). Select samples were also heat treated (95 &deg;C for 20 min) to denature and deactivate biological INPs present and digested in 10% H 2 O 2 at 95 &deg;C under UV-B for 20 min to remove any organic carbon INPs. Agresti, A, and BA Coull. 1998. "Approximate is better than &ldquo;exact&rdquo; for interval estimation of binomial proportions." American Statistics 52: 119&ndash;126. https://doi.org/10.2307/2685469 McCluskey, CS, J Ovadnevaite, M Rinaldi, J Atkinson, F Belosi, D Ceburnis, &hellip; and PJ DeMott. 2018. "Marine and Terrestrial Organic Ice-Nucleating Particles in Pristine Marine to Continentally Influenced Northeast Atlantic Air Masses." Journal of Geophysical Research: Atmospheres 123 (11): 6196&ndash;6212, https://doi.org/10.1029/2017JD028033 Vali, G. 1971. "Quantitative Evaluation of Experimental Results and the Heterogeneous Freezing Nucleation of Supercooled Liquids." Journal of the Atmospheric Sciences 28: 402-209. https://doi.org/10.1175/1520-0469(1971)028<0402:QEOERA>2.0.CO;2

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Ice nucleation measurements from DRUM impactors during TRACER campaign in the Houston TX region from July to September 2022

During TRACER, three Davis Rotating Uniform size-cut Monitors (DRUM; DRUMAir 4-DRUM) were used to collect aerosols for ice nucleation measurements in the Brooks laboratory at Texas A&M University. The three instruments were located at AMF1 in La Porte, Ancillary site in Guy (ANC), and onboard the Texas A&M Rapid Onsite Atmospheric Measurements Van (ROAM-V). ROAM-V was deployed to capture airmasses behind (maritime) and ahead (continental) of the passage of the sea-breeze front through Houston. On select sampling days, ROAM-V sampled in the morning/mid-day on the coast and then transited to a second inland site for the afternoon/evening. The suite of instruments deployed on ROAM-V included a Condensation Particle Counter (CPC; GRIMM Model 5.403 CPC), Scanning Mobility Particle Sizer (SMPS; TSI 3750 detector, TSI 3082 classifier, TSI 3088 neutralizer, TSI 3081A Differential Mobility Analyzer), Cloud Condensation Nuclei counter (Droplet Measurement Technologies CCN Counter), micro pulse lidar (Droplet Measurement Technologies Micro Pulse LiDAR (miniMPL)), and one of the Davis Rotating Uniform size-cut Monitors (DRUM; DRUMAir 4-DRUM). Before sampling at each location, the latitude and longitude were recorded using the GPS on the phone application “My Altitude”. Each DRUM sampler was operated at a flow rate of 23 LPM. Each DRUM has four stages where the aerodynamic diameter size cuts are as follows: stage 1 larger than 3 μm, stage 2 from 3 to 1.2 μm, stage 3 from 1.2 to 0.34 μm, and stage 4 from 0.34 to 0.15 μm. Pretreated aluminum foil was used as a substrate on all stages. At AMF1 and ANC sites, the DRUMs were operated on the shared aerosol inlet, with generous support of the DOE ARM site staff. These instruments rotated 24 mm per day and only contained aluminum foil substrates. Substrates were changed weekly and transported to Texas A&M for storage in -80C freezer until analysis. The DRUM onboard ROAM-V was operated at a faster rotation rate of 150 mm per day to clearly separate the multiple deployment locations for ROAM-V. For the ROAM-V DRUM, substrates were changed every deployment and transported to Texas A&M for storage in -80C freezer until analysis. At Texas A&M, ice nucleation experiments were conducted to measure the ice nucleation temperature of the ambient aerosol samples collected from the three DRUMs using our previously established procedures (Alsante et al., 2023; Fornea et al., 2009; Matthews et al., 2023). For ice nucleation, only samples collected on stage 3 were analyzed, given that these are the most relevant size (1.2 to 0.34 μm diameter) for potential ice nucleating particles. For the AMF1 and ANC sites, we cut and analyzed 2 mm (2-hour) samples. We analyzed the time periods of the AMF1 site instrument when the ROAM-V was deployed. A 72-hour period from July 11th at 23:49 through July 15th at 1:49 was analyzed from the ANC site instrument. For the ROAM-V instrument, we separated the daily samples by site location. Between 1- and 6-hour independent samples were analyzed at each location. We also cut the ROAM-V samples in half to allow for compositional analysis of the aerosol on the other half of the substrate. All viable samples from the ROAM-V were analyzed. Analysis was done using a custom-built ice nucleation apparatus, recently updated to include an array of 16 individual samples (Matthews et al., 2023). On our experimental setup, we used 100 μL of Ultra-High-Performance Liquid Chromatography (UHPLC) water (Sigma Aldrich, >99.9% purity) to wash off the aerosol from the DRUM substrate. Then, we micropipetted 2 μL droplet samples into each of the 16 wells of the array, using hydrophobically coated microscope slides. Experiments would cycle 28 times from 10 C to -40C over a 20-hour period. During analysis of the TRACER samples, nine experiments were conducted with UHPLC water process blanks that followed an identical preparation procedure. The average freezing temperature and standard deviation for these process blanks is -27.7±2.4 C. This data was collected for ARM Field Campaign AFC07023 and supported by DOE ASR grant DE-SC0021047. For any further questions, please feel free to contact the instrument PI, Sarah D. Brooks, sbrooks@tamu.edu. Alsante, A. N., Thornton, D. C., & Brooks, S. D. (2023). Ice nucleation catalyzed by the photosynthesis enzyme RuBisCO and other abundant biomolecules. Communications Earth & Environment, 4(1), 51. Fornea, A. P., Brooks, S. D., Dooley, J. B., & Saha, A. (2009). Heterogeneous freezing of ice on atmospheric aerosols containing ash, soot, and soil. Journal of Geophysical Research: Atmospheres, 114(D13). DOI:10.1029/2009JD011958 Matthews, B. H., Alsante, A. N., & Brooks, S. D. (2023). Pollen Emissions of Subpollen Particles and Ice Nucleating Particles. ACS Earth and Space Chemistry.

54 ENVIRONMENTAL SCIENCES↗