Search NASA⌕ Search

SEARCH · Search NASA

Results for “Atmospheric Model Parameterizations”

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

The contribution of Humboldt Glacier, northern Greenland, to sea-level rise through 2100 constrained by recent observations of speedup and retreat

Abstract. Humboldt Glacier, northern Greenland, has retreated and accelerated through the 21st century, raising concerns that it could be a significant contributor to future sea-level rise. We use a data-constrained ensemble of three-dimensional higher-order ice sheet model simulations to estimate the likely range of sea-level rise from the continued retreat of Humboldt Glacier. We first solve for basal traction using observed ice thickness, bed topography, and ice surface velocity from the year 2007 in a PDE-constrained (partial differential equation) optimization. Next, we impose calving rates to match mean observed retreat rates from winter 2007–2008 to winter 2017–2018 in a transient calibration of the exponent in the power-law basal friction relationship. We find that power-law exponents in the range of 1/7–1/5 – rather than the commonly used 1/3–1 – are necessary to reproduce the observed speedup over this period. We then tune an iceberg calving parameterization based on the von Mises stress yield criterion in another transient-calibration step to approximate both observed ice velocities and terminus position in 2017–2018. Finally, we use the range of basal friction relationship exponents and calving parameter values to generate the ensemble of model simulations from 2007–2100 under three climate forcing scenarios from CMIP5 (two RCP8.5 forcings, Representative Concentration Pathway) and CMIP6 (one SSP5-8.5 forcing, Shared Socioeconomic Pathway). Our simulations predict 5.2–8.7 mm of sea-level rise from Humboldt Glacier, significantly higher than a previous estimate (∼ 3.5 mm) and equivalent to a substantial fraction of the 40–140 mm predicted by ISMIP6 from the whole Greenland Ice Sheet. Our larger future sea-level rise prediction results from the transient calibration of our basal friction law to match the observed speedup, which requires a semi-plastic bed rheology. In many simulations, our model predicts the growth of a sizable ice shelf in the middle of the 21st century. Thus, atmospheric warming could lead to more retreat than predicted here if increased surface melt promotes hydrofracture of the ice shelf. Our data-constrained simulations of Humboldt Glacier underscore the sensitivity of model predictions of Greenland outlet glacier response to warming to choices of basal shear stress and iceberg calving parameterizations. Further, transient calibration of these parameterizations, which has not typically been performed, is necessary to reproduce observed behavior. Current estimates of future sea-level rise from the Greenland Ice Sheet could, therefore, contain significant biases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Assessment of Warm and Dry Bias over ARM SGP Site in E3SMv2 and E3SM-MMF

Abstract Many climate models exhibit a dry and warm bias over the central United States during the summer months, including the Energy Exascale Earth System Model (E3SM) and its Multiscale Modeling Framework (MMF) configuration. Understanding the causes of this bias is important to shine a light on this common model error and reduce the uncertainty in future projections. In this study, we use E3SMv2 and E3SM-MMF to assess how parameterized and resolved convection affect temperature and precipitation biases over the Southern Great Plains site of the Atmospheric Radiation Measurement program. Both configurations overestimate near-surface temperature and underestimate precipitation at the ARM SGP site. The bias is associated with a lack of low-level clouds during days without precipitation and too much incoming solar radiation causing the surface to warm. Low-level cloud fraction in E3SM-MMF during the nonprecipitating days is lower in comparison to E3SMv2 and observation, consistent with the larger warm bias. We also find that the underestimated precipitation can be characterized as “too frequent, too weak” in E3SMv2 and “too rare, too intense” in E3SM-MMF. These deficiencies conspire to sustain the warm and dry bias over the central United States.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Building an AI-enhanced modeling framework to address multiscale predictability challenges

Focal Area(s): Build an AI-enhanced modeling framework that integrates the three focus areas in the solicitation. Science Challenge: The 4M-2N complexities (Multiscale, Multiphysics, Multibody, Multidimension, Non-linearity, and Non-Gaussianality) of atmospheric aerosol-cloud-precipitation-turbulence-radiation system poses physical and computational challenges to further advance predictive models; We plan to address the challenges by developing an AI-enhanced modeling framework that facilitates automated calibration and improvement of subgrid parameterizations, enhance data assimilation of measurements to improve initial and boundary conditions used to drive the physical model, and optimally blends data-driven and physics-based forecasting models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Large-scale offshore wind farm effects on weather and climate in Puerto Rico (Final Technical Report)

Puerto Rico’s current electricity generation heavily relies on imported fossil fuels. This results in an average cost of electricity higher than the U.S. mainland average in all sectors (residential, commercial, and industrial), despite abundant local offshore wind resources, which have the potential to provide secure, low-cost energy generation and consequent economic prosperity. However, effects on atmospheric and oceanic circulation resulting from large-scale deployments of offshore wind farms have not been previously studied at tropical latitude. This project addressed this knowledge gap through a computational modeling effort designed to capture the coupled dynamics of the atmosphere and the ocean in presence of offshore wind farms. Results indicate that wind farm wakes can alter wind stress, generate Ekman-driven vertical transport, and potentially affect nutrient distribution. While full model coupling remains challenging, progress in parameterization and large-eddy simulations provides a foundation for future research. The project contributes to DOE’s Earth System modeling efforts and supports STEM workforce development.

17 WIND ENERGY↗
DOE OSTI2022

More Realistic Intermediate Depth Dry Firn Densification in the Energy Exascale Earth System Model (E3SM)

Abstract Earth system models account for seasonal snow cover, but many do not accommodate the deeper snowpack on ice sheets (aka firn) that slowly transforms to ice under accumulating snowfall. To accommodate and resolve firn depths of up to 60 m in the Energy Exascale Earth System Model's land surface model (ELM), we add 11 layers to its snowpack and evaluate three dry snow compaction equations in multi‐century simulations. After comparing results from ELM simulations (forced with atmospheric reanalysis) with empirical data, we find that implementing into ELM a two‐stage firn densification model produces more accurate dry firn densities at intermediate depths of 20–60 m. Compared to modeling firn using the equations in the (12 layer) Community Land Model (version 5), switching to the two‐stage firn densification model (with 16 layers) significantly decreases root‐mean‐square errors in upper 60 m dry firn densities by an average of 41 kg m −3 (31%). Simulations with three different firn density parameterizations show that the two‐stage firn densification model should be used for applications that prioritize accurate upper 60 m firn air content (FAC) in regions where the mean annual surface temperature is greater than roughly −31°C. Because snow metamorphism, firn density, and FAC are major components in modeling ice sheet surface albedo, melt water retention, and climatic mass balance, these developments advance broader efforts to simulate the response of land ice to atmospheric forcing in Earth system models.

58 GEOSCIENCES↗
DOE OSTI2022

Marine Shallow Cloud Adjustments to the Presence of Shortwave-Absorbing Aerosols: Advancing Understanding Through a Combined Analysis of Lasic Datasets and Process Modeling (Final Report)

This is a final report for this award, whose purpose was to support the analysis of data emanating from the LASIC (Layered Atlantic Smoke Interaction with Clouds) campaign. (Layered Atlantic Smoke Interactions with Clouds) is a strategy to improve our understanding of aged carbonaceous aerosol, its seasonal evolution, and the mechanisms by which clouds adjust to the presence of the aerosol. The observational strategy centers on deploying the AMF1 cloud, aerosol, and atmospheric profiling instrumentation to Ascension Island, located within the trade-wind shallow cumulus regime (14.50W, 80S) 3000 km offshore of continental Africa. The location is within the latitude zone of the maximum outflow of aerosol, with the deepening boundary layer known to entrain free-tropospheric smoke. The primary activities for LASIC are: 1) to improve current knowledge on aged biomass burning aerosol and its radiative properties as a function of the seasonal cycle; 2) to use surface-based remote sensing to sensitively interrogate the atmosphere for the relative vertical location of aerosol and clouds; 3) to improve our understanding of the cloud adjustments to the presence of shortwave-absorbing aerosol within the vertical column, both through aerosol-radiation and through aerosol-cloud interactions; 4) to aid low cloud parameterization efforts for climate models. This award serves to further these activities and includes a sub-award to Dr. Pablo Saide at UCLA for activities 3) and 4).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Climate and hydraulic traits interact to set thresholds for liana viability

Lianas, or woody vines, and trees dominate the canopy of tropical forests and comprise the majority of tropical aboveground carbon storage. These growth forms respond differently to contemporary variation in climate and resource availability, but their responses to future climate change are poorly understood because there are very few predictive ecosystem models representing lianas. We compile a database of liana functional traits (846 species) and use it to parameterize a mechanistic model of liana-tree competition. The substantial difference between liana and tree hydraulic conductivity represents a critical source of inter-growth form variation. Here, we show that lianas are many times more sensitive to drying atmospheric conditions than trees as a result of this trait difference. Further, we use our competition model and projections of tropical hydroclimate based on Representative Concentration Pathway 4.5 to show that lianas are more susceptible to reaching a hydraulic threshold for viability by 2100.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Light absorption enhancement of black carbon in a pyrocumulonimbus cloud

Abstract Pyrocumulonimbus (pyroCb) firestorm systems have been shown to inject significant amounts of black carbon (BC) to the stratosphere with a residence time of several months. Injected BC warms the local stratospheric air, consequently perturbing transport and hence spatial distributions of ozone and water vapor. A distinguishing feature of BC-containing particles residing within pyroCb smoke is their thick surface coatings made of condensed organic matter. When coated with non-refractory materials, BC’s absorption is enhanced, yet the absorption enhancement factor ( E abs ) for pyroCb BC is not well constrained. Here, we perform particle-scale measurements of BC mass, morphology, and coating thickness from inside a pyroCb cloud and quantify E abs using an established particle-resolved BC optics model. We find that the population-averaged E abs for BC asymptotes to 2.0 with increasing coating thickness. This value denotes the upper limit of E abs for thickly coated BC in the atmosphere. Our results provide observationally constrained parameterizations of BC absorption for improved radiative transfer calculations of pyroCb events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

ISMIP6-based projections of ocean-forced Antarctic Ice Sheet evolution using the Community Ice Sheet Model

Abstract. The future retreat rate for marine-based regions of the Antarctic Ice Sheet is one of the largest uncertainties in sea-level projections. The Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) aims to improve projections and quantify uncertainties by running an ensemble of ice sheet models with atmosphere and ocean forcing derived from global climate models. Here, the Community Ice Sheet Model (CISM) is used to run ISMIP6-based projections of ocean-forced Antarctic Ice Sheet evolution. Using multiple combinations of sub-ice-shelf melt parameterizations and calibrations, CISM is spun up to steady state over many millennia. During the spin-up, basal friction parameters and basin-scale thermal forcing corrections are adjusted to optimize agreement with the observed ice thickness. The model is then run forward for 550 years, from 1950–2500, applying ocean thermal forcing anomalies from six climate models. In all simulations, the ocean forcing triggers long-term retreat of the West Antarctic Ice Sheet, especially in the Filchner–Ronne and Ross sectors. Mass loss accelerates late in the 21st century and then rises steadily for several centuries without leveling off. The resulting ocean-forced sea-level rise at year 2500 varies from about 150 to 1300 mm, depending on the melt scheme and ocean forcing. Further experiments show relatively high sensitivity to the basal friction law, moderate sensitivity to grid resolution and the prescribed collapse of small ice shelves, and low sensitivity to the stress-balance approximation. The Amundsen sector exhibits threshold behavior, with modest retreat under many parameter settings but complete collapse under some combinations of low basal friction and high thermal forcing anomalies. Large uncertainties remain, as a result of parameterized sub-shelf melt rates, simplified treatments of calving and basal friction, and the lack of ice–ocean coupling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Scrutinizing entrainment and mass flux closures in shallow cumulus parameterizations using cloud-radar observations and large-eddy simulation (Final Report)

Representing warm, shallow clouds is a challenge for Global Climate Models, which leads to uncertainty in model projections of future climate. Our study was particularly concerned with how climate models represent the effects of shallow, non-precipitating cumulus clouds. Climate models must make certain assumptions when representing low clouds such as how much dry air from outside the cloud is mixed into the warm, moist cloud updraft. This process is called “entrainment” and tends to weaken the cloud updraft and dry out the cloud, compared to a theoretical cloud with no mixing. This project was a joint observational and modeling study to derive estimates of entrainment from field observations and use high-resolution modeling to evaluate the assumptions used in climate model representations of shallow cumulus. Our research took advantage of a wide variety of Department of Energy measurements from fixed climate-monitoring locations over the Southern Great Plains of the U.S. and the in the Eastern North Atlantic, along with short-term field projects (MAGIC, CAP–MBL, ACE–ENA). During the project, the PI maintained regular collaborations with Department of Energy Scientists at Brookhaven and Argonne National Laboratories. The project funded two students, who both obtained their M.S. in Atmospheric Science and are currently in the Atmospheric Science Ph.D. program. Both were lead authors on research papers funded by this grant.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Physics or Knob‐Tuning? Tropical Anvil Peak Is Captured by GCMs

Tropical anvil clouds peak near 200 hPa and significantly impact Earth's climate, yet its physical realism in coarse‐resolution General Circulation Models (GCMs) remains debated. We examine anvil cloud formation by performing simulations with a GCM with a hierarchy of cloud fraction schemes ranging from a complex prognostic Tiedtke scheme to a simple binary scheme. All schemes consistently reproduce the anvil peak. The robust anvil peak arises because extremely cold temperatures at the upper troposphere facilitate frequent saturation events, producing clouds that disproportionately influence mean cloud fraction. Sensitivity experiments with enhanced evaporation of cloud condensate unexpectedly show increased anvil coverage, highlighting how slight evaporative moistening reinforces local saturation in cold upper‐tropospheric conditions. These results demonstrate that the tropical anvil cloud peak emerges from fundamental thermodynamic constraints, rather than specific cloud fraction parameterization choices.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Using Radiative Convective Equilibrium to Explore Clouds and Climate in the Community Atmosphere Model

Abstract Characteristics of, and fundamental differences between, the radiative‐convective equilibrium (RCE) climate states following the Radiative‐Convective Equilibrium Model Intercomparison Project (RCEMIP) protocols in the Community Atmosphere Model version 5 (CAM5) and version 6 (CAM6) are presented. This paper explores the characteristics of clouds, moisture, precipitation and circulation in the RCE state, as well as the tropical response to surface warming, in CAM5 and CAM6 with different parameterizations. Overall, CAM5 simulates higher precipitation rates that result in larger global average precipitation, despite lower outgoing longwave radiation compared to CAM6. Differences in the structure of clouds, particularly the amount and vertical location of cloud liquid, exist between the CAM versions and can, in part, be related to distinct representations of shallow convection and boundary layer processes. Both CAM5 and CAM6 simulate similar peaks in cloud fraction, relative humidity, and cloud ice, linked to the usage of a similar deep convection parameterization. These anvil clouds rise and decrease in extent in response to surface warming. More generally, extreme precipitation, aggregation of convection, and climate sensitivity increase with warming in both CAM5 and CAM6. This analysis provides a benchmark for future studies that explore clouds, convection, and climate in CAM with the RCEMIP protocols now available in the Community Earth System Model. These results are discussed within the context of realistic climate simulations using CAM5 and CAM6, highlighting the usefulness of a hierarchical modeling approach to understanding model and parameterization sensitivities to inform model development efforts.

Reed, Kevin A.↗
DOE OSTI2023

A New Hybrid Mass-Flux/High-Order Turbulence Closure for Ocean Vertical Mixing

While various parameterizations of vertical turbulent fluxes at different levels of complexity have been proposed, each has its own limitations. For example, simple first-order closure schemes such as the K-Profile Parameterization (KPP) lack energetic constraints; two-equation models $k-ε$ like directly solve an equation for the turbulent kinetic energy but do not account for non-diffusive fluxes, and high-order closures that include the high-order transport terms are computationally expensive. To address these, we extend the Assumed-Distribution Higher-Order Closure (ADC) framework originally proposed for the atmospheric boundary layer and apply it to the ocean surface boundary layer. By assuming a probability distribution function relationship between the vertical velocity and tracers, all second-order and higher-order moments are exactly constructed and turbulence closure is achieved in the ADC scheme. In addition, this ADC parameterization has full energetic constraints and includes non-diffusive fluxes without the computational cost of a full higher-order closure scheme. We have tested the ADC scheme against a combination of large eddy simulation (LES), KPP, and $k-ε$ for surface buoyancy-driven convective mixing and found that the ADC scheme is robust with different vertical resolutions and compares well to the LES results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

ERF: Energy Research and Forecasting

The Energy Research and Forecasting (ERF) code is a new model that simulates the mesoscale and microscale dynamics of the atmosphere using the latest high-performance computing architectures. It employs hierarchical parallelism using an MPI+X model, where X may be OpenMP on multicore CPU-only systems, or CUDA, HIP, or SYCL on GPU-accelerated systems. ERF is built on AMReX (Zhang et al., 2019, 2021), a block-structured adaptive mesh refinement (AMR) software framework that provides the underlying performance-portable software infrastructure for block-structured mesh operations. The "energy" aspect of ERF indicates that the software has been developed with renewable energy applications in mind. In addition to being a numerical weather prediction model, ERF is designed to provide a flexible computational framework for the exploration and investigation of different physics parameterizations and numerical strategies, and to characterize the flow field that impacts the ability of wind turbines to extract wind energy. The ERF development is part of a broader effort led by the US Department of Energy's Wind Energy Technologies Office.

17 WIND ENERGY↗
DOE OSTI2021

Unprecedented cloud resolution in a GPU-enabled full-physics atmospheric climate simulation on OLCF’s summit supercomputer

Clouds represent a key uncertainty in future climate projection. While explicit cloud resolution remains beyond our computational grasp for global climate, we can incorporate important cloud effects through a computational middle ground called the Multi-scale Modeling Framework (MMF), also known as Super Parameterization. This algorithmic approach embeds high-resolution Cloud Resolving Models (CRMs) to represent moist convective processes within each grid column in a Global Climate Model (GCM). The MMF code requires no parallel data transfers and provides a self-contained target for acceleration. This study investigates the performance of the Energy Exascale Earth System Model-MMF (E3SM-MMF) code on the OLCF Summit supercomputer at an unprecedented scale of simulation. Hundreds of kernels in the roughly 10K lines of code in the E3SM-MMF CRM were ported to GPUs with OpenACC directives. A high-resolution benchmark using 4600 nodes on Summit demonstrates the computational capability of the GPU-enabled E3SM-MMF code in a full physics climate simulation.

58 GEOSCIENCES↗
DOE OSTI2023

Modeling the Impact of the Organic Aerosol Phase State on Multiphase OH Reactive Uptake Kinetics and the Resultant Heterogeneous Oxidation Timescale of Organic Aerosol in the Amazon Rainforest

Accurate predictions of chemical lifetime, i.e., oxidation timescale and change in the mass of organic aerosols (OAs) in atmospheric models, are critical to quantify the impacts of OA on aerosol–cloud interactions, radiative forcing, and air quality. The heterogeneous oxidation of OA by hydroxyl (OH) radicals is a key process governing OA mass changes during their oxidation. Recently, laboratory observations of the OH uptake coefficient (γ OH ) and heterogeneous reaction rate (k het,OH ) for different OA systems with varying phase states have been combined to develop a new parameterization of γ OH and k het,OH as a function of the OA phase state, i.e., viscosity (η OA ). Here, in this work, we use a recently published viscosity prediction framework to analyze the new γ OH and k het,OH parameterization with a box model. Subsequently, we implement this box model to predict γ OH and k het,OH over the Amazon rainforest in the dry-to-wet transition season with both significant biomass burning and biogenic influences. Relevant parameters within our box model are specified based on detailed regional model simulations over the Amazon using the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem). η OA is predicted as a function of species volatility, OA composition (including water uptake) along with ambient relative humidity (RH), and temperature. Based on ambient conditions simulated by WRF-Chem over the entire atmospheric column, we use the box model to predict the upper bounds of the heterogeneous oxidation timescale of OA. Based on previous laboratory measurements, we assume that this heterogeneous OH oxidation causes fragmentation (carbon loss) of OA, followed by the evaporation of the fragmented molecules that results in the exponential decay of OA. We predict that the oxidation timescale of OA is ~1 month near the Earth’s surface in the pristine Amazonian background because of low OH concentrations. But OA is oxidized more rapidly within urban and wildfire plumes near the Earth’s surface with an order of magnitude higher OH concentrations compared to the pristine background, causing the simulated oxidation timescale of OA to be much shorter, ~3–4 days. At 3–5 km altitudes where biomass burning OA is predicted to be semisolid, the heterogeneous oxidation timescale is estimated as ~6 days to 3 weeks and decreases with increasing OH concentrations within plumes. We show that the simulated mass loss of OA is strongly size-dependent, where smaller particles oxidize more rapidly compared to larger particles due to their greater per-particle surface area-to-volume ratio. Increasing urbanization and deforestation in the Amazon in the future might increase OH concentrations in the background Amazon, causing faster oxidation of OA. At higher altitudes above liquid clouds, especially in the upper troposphere where temperatures approach ~250 K, future measurements are needed to reduce the uncertainties related to the mass loss of OA at colder temperatures and low relative humidity conditions due to its heterogeneous OH oxidation.

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

Improving the particle dry deposition scheme in the CMAQ photochemical modeling system

Dry deposition of atmospheric aerosols in large-scale models is a critical, but highly uncertain, sink process with a strong dependence on particle size, meteorological conditions, and land surface properties. This study investigates the particle dry deposition scheme implemented in the standard Community Multiscale Air Quality (CMAQ) model v5.2.1, characterizes its underlying parameterized components with comparison to a similar scheme in a contemporary regional-scale model, and proposes two updated schemes that are then evaluated with available ambient particle deposition velocity (V d ) measurements. Both updated schemes reduce the surprisingly strong dependence of deposition velocity on the aerosol mode width, with one scheme further introducing a dependence on vegetation coverage that is broadly consistent with variability in observations between vegetated and non-vegetated surfaces. Compared to the base scheme, the updated scheme with vegetation dependence increases V d for submicron particles and decreases it for larger particles by an average of 37% and –66%, respectively. This scheme performs statistically better than the base scheme, reducing fractional biases by 56%–97% for vegetated land-use types and has roughly equivalent performance over water. The base and updated schemes are tested with three annual CMAQ (v5.2.1) simulations for the year 2011; predicted ambient aerosol concentrations are evaluated with routine monitoring network observations and predicted dry deposition fluxes are evaluated with data from the Clean Air Status and Trends Network (CASTNET). The updated scheme with vegetation dependence reduces negative fractional biases for PM 10 by 41% and positive fractional biases for PM 2.5 organic carbon by 15%. This scheme has been incorporated into the most recent publicly accessible versions of CMAQ (v5.3 and beyond) to replace the scheme used in previous versions of CMAQ (v4.5 through v5.2.1).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Challenging and Improving the Simulation of Mid-Level Mixed-Phase Clouds Over the High-Latitude Southern Ocean

Climate models exhibit major radiative biases over the Southern Ocean owing to a poor representation of mixed-phase clouds. In this study, we used the remote-sensing dataset from the Measurements of Aerosols, Radiation and Clouds over the Southern Ocean (MARCUS) campaign to assess the ability of the Weather Research and Forecasting (WRF) model to reproduce frontal clouds off Antarctica. It focuses on the modeling of thin mid-level supercooled liquid water layers which precipitate ice. The standard version of WRF produces almost fully glaciated clouds and cannot reproduce cloud top turbulence. Our work demonstrates the importance of adapting the ice nucleation parameterization to the pristine austral atmosphere to reproduce the supercooled liquid layers. Once simulated, droplets significantly impact the cloud radiative effect by increasing downwelling longwave fluxes and decreasing downwelling shortwave fluxes at the surface. The net radiative effect is a warming of snow and ice covered surfaces and a cooling of the ocean. Despite improvements in our simulations, the local turbulent circulation related to cloud-top radiative cooling is not properly reproduced, advocating for the need to develop a parameterization for top-down convection to capture the turbulence-microphysics interplay at cloud top.

54 ENVIRONMENTAL SCIENCES↗