Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cloud parameterization”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Ice‐Nucleating Particles That Impact Clouds and Climate: Observational and Modeling Research Needs

Abstract Atmospheric ice‐nucleating particles (INPs) play a critical role in cloud freezing processes, with important implications for precipitation formation and cloud radiative properties, and thus for weather and climate. Additionally, INP emissions respond to changes in the Earth System and climate, for example, desertification, agricultural practices, and fires, and therefore may introduce climate feedbacks that are still poorly understood. As knowledge of the nature and origins of INPs has advanced, regional and global weather, climate, and Earth system models have increasingly begun to link cloud ice processes to model‐simulated aerosol abundance and types. While these recent advances are exciting, coupling cloud processes to simulated aerosol also makes cloud physics simulations increasingly susceptible to uncertainties in simulation of INPs, which are still poorly constrained by observations. Advancing the predictability of INP abundance with reasonable spatiotemporal resolution will require an increased focus on research that bridges the measurement and modeling communities. This review summarizes the current state of knowledge and identifies critical knowledge gaps from both observational and modeling perspectives. In particular, we emphasize needs in two key areas: (a) observational closure between aerosol and INP quantities and (b) skillful simulation of INPs within existing weather and climate models. We discuss the state of knowledge on various INP particle types and briefly discuss the challenges faced in understanding the cloud impacts of INPs with present‐day models. Finally, we identify priority research directions for both observations and models to improve understanding of INPs and their interactions with the Earth System.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Understanding the Microphysical Control and Spatial‐Temporal Variability of Warm Rain Probability Using CloudSat and MODIS Observations

Abstract By combining measurements from MODIS and the CloudSat radar, we develop a parameterization scheme to quantify the combined microphysical controls by liquid water path (LWP) and cloud droplet number concentration (CDNC) of the probability of precipitation (PoP) in marine low cloud over tropical oceans. We demonstrate that the spatial‐temporal variation of grid‐mean in‐cloud can be largely explained by the variation of the joint probability density function of LWP and CDNC in the phase space specified by the bivariate PoP (LWP and CDNC) function. Through a series of sensitivity tests guided by this understanding, we find that in the Southeastern Pacific and Atlantic the stratocumulus to cumulus transition of the is mainly due to the variation of CDNC while the annual cycle is mainly due to the variation of LWP. The results of this study provide a viable way to diagnose the root cause of warm rain problems in global climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The Influence of Shallow Cloud Populations on Transitions to Deep Convection in the Amazon

Abstract In this study, a pair of convection-permitting (2-km grid spacing), month-long, wet-season Weather Research and Forecasting (WRF) Model simulations with and without the eddy-diffusivity mass-flux (EDMF) scheme are performed for a portion of the Green Ocean Amazon (GoAmazon) 2014/15 field campaign period. EDMF produces an ensemble of subgrid-scale convective plumes that evolve in response to the boundary layer meteorological conditions and can develop into shallow clouds. The objective of this study is to determine how different treatments of shallow cumulus clouds (i.e., with and without EDMF) impact the total cloud population and precipitation across the Amazonian rain forest, with emphasis on impacts on the likelihood of shallow-to-deep convection transitions. Results indicate that the large-scale synoptic conditions in the EDMF and control simulations are nearly identical; however, on the local scale their rainfall patterns diverge drastically and the biases decrease in EDMF. The EDMF scheme significantly increases the frequency of shallow clouds, but the frequencies of deep clouds are similar between the simulations. Deep convective clouds are tracked using a cloud-tracking algorithm to examine the impact of shallow cumulus on the surrounding ambient environment where deep convective clouds initiate. Results suggest that a rapid increase of low-level cloudiness acts to cool and moisten the low to midtroposphere during the day, favoring the transition to deep convection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Interconnection of Aerosol‐Cloud Interactions and Cloud Feedback Through Warm Rain Process

Recent research has revealed a correlation within the Coupled Model Intercomparison Project Phase 6 (CMIP6) simulations: models exhibiting more pronounced warming due to positive cloud feedback simultaneously show enhanced cooling from aerosol‐cloud interactions (ACI). However, the interplay between ACI and cloud feedback is not well understood in these models. Our study addresses this gap by modifying the autoconversion rate in two Earth system models (ESMs), elucidating how ACI could connect with cloud feedback through the warm rain process. We find that higher autoconversion rates, which are associated with stronger ACI, suppress the increase in cloud liquid water path (LWP) and cloud optical depth with warming by enhancing the precipitation efficiency over the extratropical regions, leading to a larger positive cloud feedback. This study offers new insights into the compensatory mechanism between ACI and cloud feedback through the warm rain process and highlights the importance of constraining the autoconversion parameterization in models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A riming‐dependent parameterization of scattering by snowflakes using the self‐similar Rayleigh–Gans approximation

Abstract Riming is a key process of precipitation formation in ice‐containing clouds, but quantifying riming from observations is challenging, limiting our ability to evaluate the riming process in numerical weather models. One challenge for radar observations is that riming changes both the physical properties (mass, area cross‐section) and scattering properties of ice particles. These changes need to be implemented consistently as a function of riming in radar forward operators, which are required for retrievals and model evaluation in observation space. In this study, mass–size, cross‐section area–size, and backscattering cross‐section relations are developed as a function of the normalized rime mass for aggregates composed of various monomer types (columns, dendrites, needles, plates, and rosettes). The proposed framework allows us to simulate scattering properties of aggregated ice particles consistently as a function of riming in retrievals and radar forward operators. The parameterizations are developed from a large data set of simulated rimed aggregates of different sizes and monomer crystal types. The backscattering cross‐section parameterization (the “riming‐dependent parameterization”) is evaluated for radar frequencies of 35.6 and 94.0 GHz and is based on the Self‐Similar Rayleigh–Gans approximation (SSRGA), which is increasingly used to calculate microwave scattering of ice crystals and snowflakes. Compared with parameterizations from the literature that do not consider riming, the riming‐dependent parameterization leads to significantly smaller biases in terms of backscattering cross‐section. When using the particle masses and scattering properties of the individual particles simulated by the aggregation and riming model as a reference, the bias of our parameterization is below 1 dB when integrating over an exponential particle size distribution with sizes from 0.1–10 mm.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Horizontal Resolution Sensitivity of the Simple Convection-Permitting E3SM Atmosphere Model in a Doubly-Periodic Configuration

We develop a doubly periodic version of the Simple Convection-Permitting E3SM Atmosphere Model (SCREAM) to provide an efficient configuration for this global convection permitting model (GCPM), akin to a single column model often found in conventional general circulation models. The design details are explained, in addition to the extensive case library associated with the doubly periodic SCREAM (DP-SCREAM) configuration. We demonstrate that doubly periodic cloud resolving models are useful tools to explore the horizontal resolution sensitivity of GCPMs, in addition to replicating biases seen in the global models. Using DP-SCREAM, we show that SCREAM exhibits behaviors of a scale aware model as it is able to naturally partition between sub-grid scale (SGS) and resolved vertical transport across the gray zone of turbulence. We show that SCREAM is reasonably scale insensitive when run at resolutions from 1 to 5 km, but can exhibit sensitivity, particularly for the shallow convective regime, when run at resolutions approaching that of large eddy simulations. We conclude that SGS parameterization improvements are likely needed to reduce this scale sensitivity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Structurally flexible cloud microphysics, observationally constrained at all scales via ML-accelerated Bayesian inference

We discuss the challenge of developing observationally informed parameterizations of microphysics for use at a hierarchy of modeling scales. Our proposed approach is applicable to any domain that suffers from a two-fold parameterization problem, where physical processes are not resolved at the model scale (the first problem), and those processes are uncertain at any scale (the second problem). For such problems, a physical approach facilitates modeling across scales, as well as systematic observational inference accelerated by machine learning (ML) surrogate models, and quantification of physical uncertainties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Time Evolution and Diurnal Variability of the Parametric Sensitivity of Turbine-Height Winds in the MYNN-EDMF Parameterization

The Mellor-Yamada-Nakanishi-Niino (MYNN) parameterization applied in the Weather Research and Forecasting (WRF) model has been augmented to include the Eddy-Diffusion Mass-Flux (EDMF) approach to better represent transport by boundary-layer eddies. This change includes the addition of new parameters associated with convective updrafts and boundary-layer clouds. These additions lead to new parametric sensitivities in the hub-height wind speed compared to simulations using the standard MYNN parameterization that were explored previously. This work builds on efforts focused on WRF’s MYNN parameterization by examining the sensitivity of hub-height wind speed to parameters in the MYNN-EDMF parameterization as a function of simulation lead time. Summer and winter periods were selected from the second Wind Forecast Improvement Project (WFIP2). Five sets of simulations were completed for each season, with lead times ranging from 2 to 6 days. The results show that the sensitivity to the new parameters associated with the EDMF scheme is generally small compared to other parameters, with the largest impact in convective boundary layers. The spread in the perturbed parameter ensembles was found to grow quickly over the first 8 to 19 hours in the summer simulations and 17 to 24 hours in the winter simulations with little change after that, regardless of the simulation length. A strong diurnal cycle in the parameter sensitivity was also found associated with changes in the atmospheric stability, as well as an increase in the sensitivity to the entrainment parameter used in the EDMF parameterization that is associated with increasing fractional areas covered by plumes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Small field campaign: aerosol – ice formation closure pilot study. Final Report

Prediction of atmospheric ice formation from aerosol particles by heterogeneous nucleation represents one of the grand challenges in atmospheric science. Our insufficient predictive understanding of primary ice formation is the reason that climate models typically do not include heterogeneous ice nucleation with subsequent effects on climate uncertainty. Mixed-phase clouds, where supercooled water droplets and ice crystals coexist play globally an important role regulating climate. This is especially the case for the Arctic region that experiences the greatest warming due to climate change compared to other regions in the world. Immersion freezing initiated by ice-nucleating particles (INPs) in supercooled water droplets is recognized as the dominant primary ice formation pathway in mixed-phase cloud regimes. For this reason, it is crucial to evaluate our capability to predict immersion freezing for a given ambient aerosol population. The goal of this project is conducting a field-based pilot study at the U.S. DOE Atmospheric Radiation Measurement (ARM) user facility at Southern Great Plains (SGP) to evaluate our capability to predict the number concentration of aerosol particles that serve as INPs in the immersion freezing mode. Successful prediction of INP number concentrations is also termed “closure”. This field-observational approach represents a first-of-its kind attempt of an aerosol–ice formation closure study (AEROICESTUDY). Very few closure studies related to INPs have been conducted, and to our knowledge, none using robust size-resolved ambient aerosol composition measurements as a starting point. Achievement of aerosol–ice formation closure relies on our ability to characterize the ambient aerosol population with respect to particles size and composition and to determine INP number concentrations for specified freezing temperatures. This requires numerous online and offline instrumentation resulting in this pilot field campaign being a multi-institutional and community-collaborative effort. We chose the ARM SGP megasite for this first aerosol-ice formation closure pilot study due to its significant measurement capabilities available to obtain detailed physical characterization of the local aerosol population including size distribution, mass loading, and chemical composition of non-refractory aerosol particles. The overall objective of this project is to identify ice nucleation parameterizations that produce the most robust predictions of INP numbers and thus are best suited to be included in cloud and climate models. This objective includes the following goals for this field and laboratory-based project: i) What are the crucial aerosol physicochemical property measurements needed to accurately guide ice nucleation representations in models and long-term INP measurements? ii) What level of parameter details needs to be known to achieve aerosol–ice formation closure? iii) What are the leading causes for climate model bias in INP predictions? We found that the advances in our understanding of immersion freezing garnered over the last 20 years allowed us to yield partial and full closures of atmospheric immersion freezing from ambient aerosol particles. When the aerosol population is physicochemically complex and parameterizations for representative INP types are not yet available, we still struggle to accurately predict INP number concentrations. This project clearly demonstrates that with more laboratory and field measurements that are accompanied by particle composition analysis, the necessary datasets to achieve aerosol–ice formation closure for various locations will emerge, thus providing a robust foundation for guiding the representation of INPs in cloud and climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Tendencies of Soil Microbial NO Emissions During HI‐SCALE as Predicted by a Nitrification/Denitrification Scheme

Many atmospheric chemical processes, including the formation of secondary organic aerosol (SOA), are strongly modulated by the reactions of NO and NO 2 (NO x ). Though NO x is controlled by anthropogenic emissions near urban areas, in rural areas soil microbes can be a significant contribution to NO emission globally. The relative rates of emissions of different nitrogen-containing species (e.g., NO, N 2 O, HONO, and N 2 ) are strong functions of soil properties such as temperature, moisture content, pH, and soil carbon and nitrogen pools. However, typical large-scale biogeochemical models either express these emissions simplistically, or not at all. Here we investigate the potential impact of soil NO emissions on atmospheric chemistry and SOA formation over regional and monthly time scales, specifically the 2016 spring and summer Intensive Observational Periods (IOPs) of the Holistic Interactions of Shallow Clouds, Aerosols and Land Ecosystems (HI-SCALE) field campaign. We implement the soil NO nitrification/denitrification parameterization of Rasool et al. (2019) into the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem), supplemented by a 1-km soil moisture analysis. We then simulate both IOPs over the U.S. Great Plains, evaluating against ground stations and flight data. We show that soil NO emissions can account for a large fraction of total NO x and locally increase O 3 concentrations by up to 25%, while alleviating negative biases of gases and aerosols toward observations. Soil moisture and temperature changes between IOP1 and IOP2 lead to overall differences in emissions, but with large regional variability due to heterogeneous surface characteristics.

Atmoshpheric Chemistry↗
DOE OSTI2025

Dark magnetohydrodynamics: Black hole accretion in superradiant dark photon clouds

Black holes threaded by massive vector fields can be subject to a superradiant instability, growing a cloud of massive vector particles around it. In this work, we consider what happens if such a dark matter candidate field mimicking a dark photon interacts with an accretion flow onto the black hole. By including a kinetic mixing term with the standard model photon, we extend the commonly used equations of general-relativistic magnetohydrodynamics to a dark photon constituent. The coupling to the dark photon then appears as an effective dynamo term together with a dark Lorentz force acting on the accreting matter. We numerically study the interactions between the superradiant dark photon cloud and the inner accretion flow by solving the coupled system in full numerical relativity. By parameterically varying the mixing parameter between the dark and standard model sector, we provide a first investigation of how the accretion flow could be modified. In conclusion, depending on the coupling strength, our solutions exhibit increased wind launching, as well as oscillation modes in the disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Can general circulation models (GCMs) represent cloud liquid water path adjustments to aerosol–cloud interactions?

Abstract. General circulation models (GCMs), unlike other lines of evidence, indicate that anthropogenic aerosols cause a global-mean increase in cloud liquid water path (ℒ) and thus a negative adjustment to radiative forcing of the climate by aerosol–cloud interactions. In part 1 of this series of papers, we showed that this is true even in models that reproduce the negative correlation observed in present-day internal variability in ℒ and cloud droplet number concentration (Nd). We studied several possible confounding mechanisms that could explain the noncausal cloud–aerosol correlations in GCMs and that possibly contaminate observational estimates of radiative adjustments. Here, we perform single-column and full-atmosphere GCM experiments to investigate the causal model-physics mechanisms underlying the model radiative adjustment estimate. We find that both aerosol–cloud interaction mechanisms thought to be operating in real clouds – precipitation suppression and entrainment evaporation enhancement – are active in GCMs and behave qualitatively in agreement with physical process understanding. However, the modeled entrainment enhancement has a negligible global-mean effect. This raises the question of whether the GCM estimate is incorrect due to parametric or base-state representation errors or whether the process understanding gleaned from a limited set of canonical cloud cases is insufficiently representative of the diversity of clouds in the real climate. Regardless, even at limited resolution, the GCM physics appears able to parameterize the small-scale microphysics–turbulence interplay responsible for the entrainment enhancement mechanism. We suggest ways to resolve tension between current and future (storm-resolving) global modeling systems and other lines of evidence in synthesis climate projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Cloud microphysical response to entrainment and mixing is locally inhomogeneous and globally homogeneous: Evidence from the lab

Entrainment of dry air into clouds strongly influences cloud optical and precipitation properties and the response of clouds to aerosol perturbations. The response of cloud droplet size distributions to entrainment–mixing is examined in the Pi convection-cloud chamber that creates a turbulent, steady-state cloud. The experiments are conducted by injecting dry air with temperature (T e ) and flow rate (Q e ) through a flange in the top boundary, into the otherwise well-mixed cloud, to mimic the entrainment–mixing process. Due to the large-scale circulation, the downwind region is directly affected by entrained dry air, whereas the upwind region is representative of the background conditions. Droplet concentration (Cn) and liquid water content (L) decrease in the downwind region, but the difference in the mean diameter of droplets (D m ) is small. The shape of cloud droplet size distributions relative to the injection point is unchanged, to within statistical uncertainty, resulting in a signature of inhomogeneous mixing, as expected for droplet evaporation times small compared to mixing time scales. As T e and Q e of entrained air increase, however, Cn, L, and D m of the whole cloud system decrease, resulting in a signature of homogeneous mixing. The apparent contradiction is understood as the cloud microphysical responses to entrainment and mixing differing on local and global scales: locally inhomogeneous and globally homogeneous. This implies that global versus local sampling of clouds can lead to seemingly contradictory results for mixing, which informs the long-standing debate about the microphysical response to entrainment and the parameterization of this process for coarse-resolution models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Semi-Analytic Functions to Calculate the Deposition Coefficients for Ice Crystal Vapor Growth in Bin and Bulk Microphysical Models

Numerical cloud models require estimates of the vapor growth rate for ice crystals. Current bulk and bin microphysical parameterizations generally assume that vapor growth is diffusion limited, though some parameterizations include the influence of surface attachment kinetics through a constant deposition coefficient. A parameterization for variable deposition coefficients is provided herein. The parameterization is an explicit function of the ambient ice supersaturation and temperature, and an implicit function of crystal dimensions and pressure. The parameterization is valid for variable surface types including growth by dislocations and growth by step nucleation. Deposition coefficients are predicted for the two primary growth directions of crystals, allowing for the evolution of the primary habits. Comparisons with benchmark calculations of instantaneous mass growth indicate that the parameterization is accurate to within a relative error of 1%. Parcel model simulations using Lagrangian microphysics as a benchmark indicate that the bulk parameterization captures the evolution of mass mixing ratio and fall speed with typical relative errors of less than 10%, whereas the average axis lengths can have errors of up to 20%. The bin model produces greater accuracy with relative errors often less that 10%. The deposition coefficient parameterization can be used in any bulk and bin scheme, with low error, if an equivalent volume spherical radius is provided.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Warm-phase microphysical evolution in large-eddy simulations of tropical cumulus congestus: evaluating drop size distribution evolution using polarimetry retrievals, in situ measurements, and a thermal-based framework

Owing to uncertainties in convective microphysics processes, improving parameterizations in Earth system models (ESMs) can benefit from observationally constrained cases suitable for scaling between cloud-resolving models and ESMs. We propose a benchmark large-eddy simulation (LES) cumulus congestus case study from the NASA Cloud, Aerosol, and Monsoon Processes Philippines Experiment (CAMP 2 Ex) for evaluating and improving ESMs in single-column model (SCM) mode. We seek observational constraints using novel polarimetric retrievals and in situ cloud microphysics measurements. Simulations using bulk and bin microphysics initialized with observed aerosol profiles are compared to cloud-top retrievals of cloud droplet effective radius (R eff ), effective variance (ν eff ), and number concentration (N d ) from the airborne Research Scanning Polarimeter (RSP). Both schemes reproduce characteristics of cloud-top N d and R eff that increase and decrease with altitude, respectively. Cloud-top N d is low-biased relative to RSP retrievals in both schemes, potentially due to limitations in both simulations and retrieval assumptions. Cloud-top R eff is low-biased in the bulk scheme but reasonably reproduced by the bin scheme. Profiles of N d and R eff are sensitive to the collision–coalescence process and the vertical variation in aerosol size distribution. Comparison of simulated and in situ droplet size distributions (DSDs) shows that, to first order, integrated moments are always sensitive to sizes < ~ 30 µm and can also be sensitive to larger sizes if the DSDs are sufficiently broad, with implications for the assumed maximum observed size retrieved by the RSP. The bin scheme captures the observed extended tail of the DSD, while the bulk scheme is unable to due to parametric constraints. Differences in expected relationships between in situ measurements of cloud cores and cloud-top retrievals by RSP demonstrate difficulty in constraining well the case presented herein. Finally, a thermal-tracking framework demonstrates that the dilution of N d throughout a thermal's lifetime is heavily determined by collision–coalescence and the height-varying aerosol distribution and that, in the absence of these, the impact of entrainment on diluting N d is largely offset by secondary activation. Implications for evaluating warm-phase convective microphysics schemes in ESMs and translating results for use on global, space-based polarimetry platforms are discussed.

Stanford, McKenna Wallace [Columbia Univ., New Yor↗
DOE OSTI2022

OCEANFILMS (Organic Compounds from Ecosystems to Aerosols: Natural Films and Interfaces via Langmuir Molecular Surfactants) sea spray organic aerosol emissions – implementation in a global climate model and impacts on clouds

Sea spray aerosol is one of the major sources of atmospheric particulate matter globally. It has increasingly been recognized that organic matter derived from ocean biological precursors contributes significantly to the composition of submicron sea spray and may modify sea spray aerosol impacts on clouds and climate. This paper describes the implementation of the OCEANFILMS (Organic Compounds from Ecosystems to Aerosols: Natural Films and Interfaces via Langmuir Molecular Surfactants) parameterization for sea spray organic aerosol emissions in a global Earth system model, the Energy Exascale Earth System Model (E3SM). OCEANFILMS is a physically based model that links sea spray chemistry with ocean biogeochemistry using a Langmuir partitioning approach. We describe the implementation details of OCEANFILMS within E3SM, compare simulated aerosol fields with observations, and investigate impacts on simulated clouds and climate. Four sensitivity cases are tested, in which organic emissions either strictly add to or strictly replace sea salt emissions (in mass and number) and are either fully internally or fully externally mixed with sea salt. The simulation with internally mixed, added organics agrees reasonably well with observed seasonal cycles of organic matter in marine aerosol and has been selected as the default configuration of the E3SM. In this configuration, sea spray organic aerosol contributes an additional source of cloud condensation nuclei, adding up to 30 cm -3 to Southern Ocean boundary-layer cloud condensation nuclei concentrations (supersaturation = 0.1 %). The addition of this new aerosol source strengthens shortwave radiative cooling by clouds by -0.36 W m -2 in the global annual mean and contributes more than -3.5 W m -2 to summertime zonal mean cloud forcing in the Southern Ocean, with maximum zonal mean impacts of about -4 W m -2 around 50–60°S. This is consistent with a previous top-down, satellite-based empirical estimate of the radiative forcing by sea spray organic aerosol over the Southern Ocean. Through its mechanistic approach, OCEANFILMS offers a path towards improved understanding of the feedbacks between ocean biology, sea spray organic matter, and climate.

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

Ice Nucleating Particles, Aerosols and Clouds over the Higher Latitude Southern Ocean

Improved understanding of aerosol-cloud-precipitation interactions is a key goal of the DOE-ASR program. Ice Nucleating Particles (INPs) are rare aerosol particles that trigger ice formation in clouds. By doing so they may initiate precipitation in clouds, and this impacts cloud lifetime. Nowhere may the role of INPs be manifested so dramatically as in the Southern Ocean (SO). This vast region was the focus of the DOE-ARM MARCUS (Measurement of Aerosols, Radiation and CloUds over the Southern Oceans) and MICRE (Macquarie Island Cloud and Radiation Experiment) campaigns. In the SO, a deficit of INPs has been a leading hypothesis for why liquid clouds persist, which underlays the bias of global climate models in predicting excess shortwave radiation reaching the ocean surface between 55°S and Antarctica. In this study, we used the 6-month acquisition of INP data on four MARCUS cruises from Hobart, Tasmania, to Antarctica, and the single-location full annual INP cycle observed in MICRE, to gain a holistic, quantitative picture of the INP number concentrations as a function of temperature and their variability over SO latitudes from 43°S to Antarctica. Augmenting initial processing of filter collections of aerosol particles during MARCUS and MICRE, additional freezing studies following thermal treatment (to removes INPs associated with microbes) and peroxide treatment (to remove all organics), as well as ionic chemistry and total organic carbon analysis, were applied to samples to improve time and space resolution, and differentiate marine from terrestrial contributions. We also performed Next Generation Sequencing to determine that bacterial composition can tag periods of enhanced and degraded marine organic INPs and the relative lack of occurrence of land-sourced bio-particles. Through these and other analyses, INP data were categorized into representative source types (e.g., marine versus terrestrially-influenced). The products of our analyses are being used to parameterize ice formation for use in numerical modeling studies to determine if, and when, ocean-derived INPs control the microphysical composition and radiative balance of SO clouds in a manner not presently captured by global models. The specific objectives and approaches proposed were largely accomplished, including: 1) Additional and value-added chemical and biological analyses of archived samples of aerosols were completed to categorize and quantify INP types and concentrations over the SO. 2) Analyses were completed to place INP data in meteorological and aerosol context. 3) Sea spray aerosol INPs were demonstrated to dominate the SO region in all but episodic events. 4) New parameterizations for INP sources relevant to the SO were constructed. 5) Collaborative modeling studies demonstrating the crucial role of INPs in determining cloud radiative and precipitation properties were conducted with MARCUS and MICRE partners. Collaborative publication submissions in this regard ensued within the two year study. This work will advance the science of interactions of aerosols, clouds and precipitation, to improve their representation in regional and global climate models. Results have informed representation of primary ice crystal nucleation, give inference to the role of secondary processes, and improve investigations of aerosol influences on clouds via ice nucleation over SO high latitudes and similar vast ocean regions. Application will improve representation of such interactions for clouds in both regional and global climate models. The data base and analyses methods applied will continue to serve research studies in the future.

58 GEOSCIENCES↗
DOE OSTI2021

Nonparametric, data-based kernel interpolation for particle-tracking simulations and kernel density estimation

Traditional interpolation techniques for particle tracking include binning and convolutional formulas that use pre-determined (i.e., closed-form, parameteric) kernels. In many instances, the particles are introduced as point sources in time and space, so the cloud of particles (either in space or time) is a discrete representation of the Green’s function of an underlying PDE. As such, each particle is a sample from the Green’s function; therefore, each particle should be distributed according to the Green’s function. In short, the kernel of a convolutional interpolation of the particle sample “cloud” should be a replica of the cloud itself. This idea gives rise to an iterative method by which the form of the kernel may be discerned in the process of interpolating the Green’s function. When the Green’s function is a density, this method is broadly applicable to interpolating a kernel density estimate based on random data drawn from a single distribution. We formulate and construct the algorithm and demonstrate its ability to perform kernel density estimation of skewed and/or heavy-tailed data including breakthrough curves.

42 ENGINEERING↗