Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cloud 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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Quantifying and Modeling the Impact of Phase State on the Ice Nucleation Abilities of 2-Methyltetrols as a Key Component of Secondary Organic Aerosol Derived from Isoprene Epoxydiols

Organic aerosols (OAs) may serve as ice-nucleating particles (INPs), impacting the formation and properties of cirrus clouds when their phase state and viscosity are in the semisolid to glassy range. However, there is a lack of direct parameterization between aerosol viscosity and their ice nucleation capabilities. In this study, we experimentally measured the ice nucleation rate of 2-methyltetrols (2-MT) aerosols, a key component of isoprene-epoxydiol-derived secondary organic aerosols (IEPOX-SOA), at different viscosities. These results demonstrate that the phase state has a significant impact on the ice nucleation abilities of OA under typical cirrus cloud conditions, with the ice nucleation rate increasing by 2 to 3 orders of magnitude when the phase state changes from liquid to semisolid. An innovative parametric model based on classical nucleation theory was developed to directly quantify the impact of viscosity on the heterogeneous nucleation rate. This model accurately represents our laboratory measurement and can be implemented into climate models due to its simple, equation-based form. Based on data collected from the ACRIDICON-CHUVA field campaign, our model predicts that the INP concentration from IEPOX-SOA can reach the magnitude of 1 to tens per liter in the cirrus cloud region impacted by the Amazon rainforest, consistent with recent field observations and estimations. This novel parameterization framework can also be applied in regional and global climate models to further improve representations of cirrus cloud formation and associated climate impacts.

2-methyltetrol↗
DOE OSTI2020

The Strong Impact of Weak Horizontal Convergence on Continental Shallow Convection

Idealized large-eddy simulation (LES) is a basic tool for studying three-dimensional turbulence in the planetary boundary layer. LES is capable of providing benchmark solutions for parameterization development efforts. However, real small-scale atmospheric flows develop in heterogeneous and transient environments with locally varying vertical motions inherent to open multiscale interactive dynamical systems. These variations are often too subtle to detect them by state-of-the-art remote and in situ measurements, and are typically excluded from idealized simulations. The present study addresses the impact of weak [i.e., O(10-6) s-1] short-lived low-level large-scale convergence/divergence perturbations on continental shallow convection. The results show a strong response of shallow nonprecipitating convection to the applied weak large-scale dynamical forcing. Evolutions of CAPE, mean liquid water path, and cloud-top heights are significantly affected by the imposed convergence/divergence. In contrast, evolving cloud-base properties, such as the area coverage and mass flux, are only weakly affected. To contrast those impacts with microphysical sensitivity, the baseline simulations are perturbed assuming different observationally based cloud droplet number concentrations and thus different rainfall. For the tested range of microphysical perturbations, the imposed convergence/divergence provides significantly larger impact than changes in the cloud microphysics. Simulation results presented here provide a stringent test for convection parameterizations, especially important for large-scale models progressing toward resolving some nonhydrostatic effects.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

Arctic Cloud‐Base Ice Precipitation Properties Retrieved Using Bayesian Inference

Abstract Cloud‐climate feedbacks are still the greatest source of uncertainty in current climate projections. Arctic clouds, which are predominantly stratiform and supercooled, often long‐lived, and nearly continuously precipitate ice particles, contribute roughly 10% of the uncertainty attributed to the global cloud feedback. This Arctic cloud uncertainty is driven by incomplete observational and theoretical knowledge required to estimate and explain the state and active processes occurring in those clouds. A focus on ice precipitation properties at Arctic cloud base rather than the surface deconfounds the product of cloud condensate sink processes from the influence of the atmospheric thermodynamic state below cloud base, rendering cloud‐base properties a more appealing target for inference and evaluation of model simulations. Here I describe an inverse model for the estimation of cloud base ice precipitation properties over Utqiagvik, North Slope of Alaska, using the synthesis of ground‐based radar and lidar measurements. By leveraging a Markov Chain Monte Carlo algorithm as the core of the inverse model, a wide range of particle size distributions are sampled, and different combinations of ice habit models are examined, both of which are typically fixed in other retrieval methods. Results show intriguing links between different cloud base thermodynamic and ice precipitation properties. Apparent ice number concentration enhancements at temperatures of −5 and −15°C suggest possible secondary ice production (SIP). The analysis alludes to an overestimation of SIP occurrence and intensity, especially in studies relying only on radar or lidar measurements. Finally, reflectivity‐dependent ice precipitation rate and ice water content parameterizations are presented.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Vertical overlap of probability density functions of cloud and precipitation hydrometeors: CLOUD AND PRECIPITATION PDF OVERLAP

Coarse-resolution climate models increasingly rely on probability density functions (PDFs) to represent subgrid-scale variability of prognostic variables. While PDFs characterize the horizontal variability, a separate treatment is needed to account for the vertical structure of clouds and precipitation. When sub-columns are drawn from these PDFs for microphysics or radiation parameterizations, appropriate vertical correlations must be enforced via PDF overlap specifications. This study evaluates the representation of PDF overlap in the Subgrid Importance Latin Hypercube Sampler (SILHS) employed in the assumed PDF turbulence and cloud scheme called the Cloud Layers Unified By Binormals (CLUBB). PDF overlap in CLUBB-SILHS simulations of continental and tropical oceanic deep convection is compared with overlap of PDF of various microphysics variables in cloud-resolving model (CRM) simulations of the same cases that explicitly predict the 3D structure of cloud and precipitation fields. CRM results show that PDF overlap varies significantly between different hydrometeor types, as well as between PDFs of mass and number mixing ratios for each species, - a distinction that the current SILHS implementation does not make. In CRM simulations that explicitly resolve cloud and precipitation structures, faster falling species, such as rain and graupel, exhibit significantly higher coherence in their vertical distributions than slow falling cloud liquid and ice. These results suggest that to improve the overlap treatment in the sub-column generator, the PDF correlations need to depend on hydrometeor properties, such as fall speeds, in addition to the currently implemented dependency on the turbulent convective length scale.

58 GEOSCIENCES↗
DOE OSTI2021

“Experimental investigation of the governing parameters of atmospheric ice nucleation using field-collected and laboratory generated aerosol particles and its application in cloud resolving models” (Final Report)

The objective of this research project is to improve our understanding of the role of aerosol particles acting as ice-nucleating particles (INPs) which in turn define the mixed-phase and cirrus cloud radiative properties and thus climate. The focus is placed on how the physicochemical particle population properties determine the particles’ ability to initiate ice nucleation. This research project combined micro-spectroscopic particle analysis, experimental ice nucleation studies, and model sensitivity studies to advance our predictive capability of the formation of mixed-phase and cirrus clouds. These project activities have led to new ice nucleation data from laboratory generated and ambient (authentic) aerosol particles furthering process-level understanding, insights in the role of organic aerosol in ice formation, and advancements in the interpretation and parameterization of ice nucleation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Exposing Process‐Level Biases in a Global Cloud Permitting Model With ARM Observations

The emergence of global convective‐permitting models (GCPMs) represents a significant advancement in climate modeling, offering improved representation of deep convection and complex precipitation patterns. In this study, we evaluate the performance of the Simple Cloud‐Resolving E3SM Atmosphere Model (SCREAM) using its doubly periodic configuration (DP‐SCREAM) against large eddy simulations and modern observational data sets from the Atmospheric Radiation Measurement program. We introduce several new transitional cloud regime cases, such as the transition from shallow to deep convection and from stratocumulus to cumulus, as well as cold‐air outbreak scenarios. The results reveal both strengths and limitations of SCREAM, particularly in the accurate simulation of cloud transitions and midlevel convection, with varying degrees of sensitivity to horizontal and vertical resolution. Despite improvements at higher resolutions, key biases remain, including the abrupt transition from shallow to deep convection and the lack of congestus clouds. These findings underscore the need for further refinement in turbulence parameterizations and vertical grid resolution in GCPMs.

Bogenschutz, Peter A. [Lawrence Livermore National↗
DOE OSTI2026

Improving Process Level Understanding of Boundary Layer Winds over the Northeast U.S. Shelf: The Third Wind Forecast Improvement Project (WFIP3)

The third Wind Forecast Improvement Project (WFIP3), a U.S. Department of Energy and National Oceanic and Atmospheric Administration sponsored investigation, sought to improve understanding of the physical phenomena in the atmosphere and ocean that dictate the structure and variability of wind and thermodynamic fields within the marine atmospheric boundary layer. WFIP3 focused on mesoscale and submesoscale flows -- including sea breezes, low-level jets, low-level clouds, and coastal storms -- and the ability of advanced numerical model parameterizations to represent them within fully coupled oceanic and atmospheric modeling systems and foundational weather forecast models. WFIP3 conducted a comprehensive 18-month observational study over the Northeast U.S. outer continental shelf, a high-use coastal zone, using a 3D multiscale sensor array to highly resolve the temporal, vertical, and horizontal structure of the coupled atmospheric and oceanic boundary layers. Multiple land-based study sites adjacent to the coastal ocean observed surface meteorology and vertical profiles of atmospheric properties via passive infrared and microwave radiometers, active lidars and radars, and radiosondes. At sea, an array of surface flux buoys and two vertical profiling lidar buoys observed both atmospheric and oceanic properties, augmented by land-based oceanographic radar systems and routine ship-based surveys. Intensive observations of the marine atmospheric boundary layer over the ocean was done from an air-sea interaction flux tower and extended deployments of a large autonomous barge platform. Numerous critical forecasting phenomena were observed that are being evaluated within regional coupled and uncoupled modeling systems, including the National Oceanic and Atmospheric Administration's foundational High-Resolution Rapid Refresh forecast model.

Kirincich, Anthony↗
DOE OSTI2026

Improving High Resolution Offshore Wind Resource Assessments and Forecasts using Observations in the MA/RI Lease Areas

The third Wind Forecast Improvement Project (WFIP3) sought to improve understanding of the physical phenomena in the atmosphere and ocean that dictate the structure and variability of wind and thermodynamic fields within the Marine Atmospheric Boundary Layer (MABL). WFIP3 conducted a comprehensive 18-month observational study over the Northeast U.S. outer continental shelf, a high use coastal zone, using a 3-D multiscale sensor array to highly resolve the temporal, vertical, and horizontal structure of the coupled atmospheric and oceanic boundary layers. Multiple land-based study sites adjacent to the coastal ocean observed surface meteorology and vertical profiles of atmospheric properties via passive infrared and microwave radiometers, active lidars and radars, and radiosondes. At sea, an array of surface flux buoys and two vertical profiling lidar buoys observed both atmospheric and oceanic properties, augmented by land-based oceanographic radar systems and routine ship-based surveys. Intensive observations of the MABL over the ocean was done from an air-sea interaction flux tower and extended deployments of a large barge platform. WFIP3 focused on mesoscale and sub-mesoscale flows -- including sea breezes, low-level jets, low-level clouds, and coastal storms -- and the ability of advanced numerical model parameterizations to represent them within fully coupled oceanic and atmospheric modeling systems and foundational weather forecast models. Numerous critical forecasting phenomena were observed, however the project was terminated prior to the completion of the field observational period and the analysis period.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Modeling Condensation in Deep Convection

Cloud-scale models apply two drastically different methods to represent condensation of water vapor to form and grow cloud droplets. Maintenance of water saturation inside liquid clouds is assumed in the computationally efficient saturation adjustment approach used in most bulk microphysics schemes. When super- or subsaturations are allowed, condensation/evaporation can be calculated using the predicted saturation ratio and (either predicted or prescribed) mean droplet radius and concentration. Here, we investigate differences between simulations of deep unorganized convection applying a saturation adjustment condensation scheme (SADJ) and a scheme with supersaturation prediction (SPRE). A double-moment microphysics scheme with CCN activation parameterized as a function of the local vertical velocity is applied to compare cloud fields simulated applying SPRE and SADJ. Clean CCN conditions are assumed to demonstrate upper limits of the SPRE and SADJ difference. Microphysical piggybacking is used to extract the impacts with confidence. Results show a significant impact on deep convection dynamics, with SADJ featuring more cloud buoyancy and thus stronger updrafts. This leads to around a 3% increase of the surface rain accumulation in SADJ. Upper-tropospheric anvil cloud fractions are much larger in SPRE than in SADJ because of the higher ice concentrations and thus longer residence times of anvil particles in SPRE, as demonstrated by sensitivity tests. Higher ice concentrations in SPRE come from significantly larger ice supersaturations in strong convective updrafts that feature water supersaturations of several percent.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Shallow Cumulus in WRF Parameterizations Evaluated against LASSO Large-Eddy Simulations

Abstract Representation of shallow cumulus is a challenge for mesoscale numerical weather prediction models. These cloud fields have important effects on temperature, solar irradiance, convective initiation, and pollutant transport, among other processes. Recent improvements to physics schemes available in the Weather Research and Forecasting (WRF) Model aim to improve representation of shallow cumulus, in particular over land. The DOE LES ARM Symbiotic Simulation and Observation Workflow (LASSO) project provides several cases that we use here to test the new physics improvements. The LASSO cases use multiple large-scale forcings to drive large-eddy simulations (LES), and the LES output is easily compared to output from WRF single-column simulations driven with the same initial conditions and forcings. The new Mellor–Yamada–Nakanishi–Niino (MYNN) eddy diffusivity mass-flux (EDMF) boundary layer and shallow cloud scheme produces clouds with timing, liquid water path (LWP), and cloud fraction that agree well with LES over a wide range of those variables. Here we examine those variables and test the scheme’s sensitivity to perturbations of a few key parameters. We also discuss the challenges and uncertainties of single-column tests. The older, simpler total energy mass-flux (TEMF) scheme is included for comparison, and its tuning is improved. This is the first published use of the LASSO cases for parameterization development, and the first published study to use such a large number of cases with varying cloud amount. This is also the first study to use a more precise combined infrared and microwave retrieval of LWP to evaluate modeled clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

ICON-A, the Atmosphere Component of the ICON Earth System Model: I. Model Description

ICON–A is the new icosahedral nonhydrostatic (ICON) atmospheric general circulation model in a configuration using the Max Planck Institute physics package, which originates from the ECHAM6 general circulation model, and has been adapted to account for the changed dynamical core framework. The coupling scheme between dynamics and physics employs a sequential updating by dynamics and physics, and a fixed sequence of the physical processes similar to ECHAM6. To allow a meaningful initial comparison between ICON–A and the established ECHAM6–LR model, a setup with similar, low resolution in terms of number of grid points and levels is chosen. The ICON–A model is tuned on the base of the Atmospheric Model Intercomparison Project (AMIP) experiment aiming primarily at a well balanced top–of atmosphere energy budget to make the model suitable for coupled climate and Earth system modeling. The tuning addresses first the moisture and cloud distribution to achieve the top–of–atmosphere energy balance, followed by the tuning of the parameterized dynamic drag aiming at reduced wind errors in the troposphere. The resulting version of ICON–A has overall biases, which are comparable to those of ECHAM6. Problematic specific biases remain in the vertical distribution of clouds and in the stratospheric circulation, where the winter vortices are too weak. Biases in precipitable water and tropospheric temperature are, however, reduced compared to the ECHAM6. As a result, ICON–A will serve as the basis of further development and as the atmosphere component to the coupled model, ICON–Earth system model (ESM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

A flexible importance sampling method for integrating subgrid processes

Numerical models of weather and climate need to compute grid-box-averaged rates of physical processes such as microphysics. These averages are computed by integrating subgrid variability over a grid box. For this reason, an important aspect of atmospheric modeling is spatial integration over subgrid scales. The needed integrals can be estimated by Monte Carlo integration. Monte Carlo integration is simple and general but requires many evaluations of the physical process rate. To reduce the number of function evaluations, this paper describes a new, flexible method of importance sampling. It divides the domain of integration into eight categories, such as the portion that contains both precipitation and cloud, or the portion that contains precipitation but no cloud. It then allows the modeler to prescribe the density of sample points within each of the eight categories. The new method is incorporated into the Subgrid Importance Latin Hypercube Sampler (SILHS). The resulting method is tested on drizzling cumulus and stratocumulus cases. In the cumulus case, the sampling error can be considerably reduced by drawing more sample points from the region of rain evaporation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A 1D Model for Nucleation of Ice From Aerosol Particles: An Application to a Mixed‐Phase Arctic Stratus Cloud Layer

Abstract Mixed‐phase clouds (MPCs) have been identified as significant contributors to uncertainties in climate projections, attributable to model representation of processes controlling the formation and loss of supercooled water droplets and ice particles from the atmosphere. Arctic MPCs are commonly widespread and long‐lived, with sustained ice crystal formation processes that challenge current understanding. This study examines the ice‐nucleating particle (INP) reservoir dynamics governing immersion‐mode heterogeneous freezing in an observed case of Arctic MPCs using a simplified 1D aerosol‐cloud model. The model setup includes prescribed dynamical forcings and thermodynamic profiles, and represents INPs as multicomponent and polydisperse particle size distributions. Diagnostic and prognostic approaches to immersion freezing parameterization are compared, including time‐independent (singular) number‐ and surface area‐based descriptions and a time‐dependent description following classical nucleation theory (CNT). The choice of freezing parameterization defines the size of the INP reservoir. The CNT‐based description yields an orders of magnitude larger INP reservoir than the singular parameterizations, which is the dominant factor for sustained ice crystal formation. The efficiency of the freezing process and cloud cooling are of secondary importance. A diagnostic treatment neglecting INP loss is only accurate when the INP reservoir size is large and INP depletion weak. Since a larger INP reservoir sustains ice crystal formation substantially longer, and ice water path scales with ice crystal concentrations for the conditions considered, resolving the source of differences in INP reservoir dynamics due to model implementation is a high priority for advancing climate model physics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Long-term observations of cloud condensation nuclei over the Amazon rain forest - Part 2: Variability and characteristics of biomass burning, long-range transport, and pristine rain forest aerosols

For this study, size-resolved measurements of atmospheric aerosol and cloud condensation nuclei (CCN) concentrations and hygroscopicity were conducted over a full seasonal cycle at the remote Amazon Tall Tower Observatory (ATTO, March 2014–February 2015). In a preceding companion paper, we presented annually and seasonally averaged data and parametrizations (Part 1; Pöhlker et al., 2016a). In the present study (Part 2), we analyze key features and implications of aerosol and CCN properties for the following characteristic atmospheric conditions: Empirically pristine rain forest (PR) conditions, where no influence of pollution was detectable, as observed during parts of the wet season from March to May. The PR episodes are characterized by a bimodal aerosol size distribution (strong Aitken mode with D Ait ≈ 70nm and N Ait ≈ 160cm -3 , weak accumulation mode with D acc ≈ 160nm and N acc ≈ 90cm -3 ), a chemical composition dominated by organic compounds, and relatively low particle hygroscopicity (κ Ait ≈ 0.12, κ acc ≈ 0.18). Long-range-transport (LRT) events, which frequently bring Saharan dust, African biomass smoke, and sea spray aerosols into the Amazon Basin, mostly during February to April. The LRT episodes are characterized by a dominant accumulation mode ( DAit ≈ 80nm, N Ait ≈ 120cm -3 vs. D acc ≈ 180nm, N acc ≈ 310cm -3 ), an increased abundance of dust and salt, and relatively high hygroscopicity (κ Ait ≈ 0.18, κ acc ≈ 0.35). The coarse mode is also significantly enhanced during these events. Biomass burning (BB) conditions characteristic for the Amazonian dry season from August to November. The BB episodes show a very strong accumulation mode (D Ait ≈ 70nm, N Ait ≈ 140cm -3 vs. D acc ≈ 170nm, N acc ≈ 3400cm -3 ), very high organic mass fractions (~90%), and correspondingly low hygroscopicity (κ Ait ≈ 0.14, κ acc ≈ 0.17). Mixed-pollution (MPOL) conditions with a superposition of African and Amazonian aerosol emissions during the dry season. During the MPOL episode presented here as a case study, we observed African aerosols with a broad monomodal distribution (D ≈ 130nm, N CN, 10 ≈ 1300cm -3 ), with high sulfate mass fractions (~20%) from volcanic sources and correspondingly high hygroscopicity (κ < 100 nm ≈ 0.14, κ > 10 nm ≈ 0.22), which were periodically mixed with fresh smoke from nearby fires (D ≈ 110nm, N CN, 10 ≈ 2800cm -3 ) with an organic-dominated composition and sharply decreased hygroscopicity (κ < 150 nm ≈ 0.10, κ > 150 nm ≈ 0.20). Insights into the aerosol mixing state are provided by particle hygroscopicity (κ) distribution plots, which indicate largely internal mixing for the PR aerosols (narrow κ distribution) and more external mixing for the BB, LRT, and MPOL aerosols (broad κ distributions). The CCN spectra (CCN concentration plotted against water vapor supersaturation) obtained for the different case studies indicate distinctly different regimes of cloud formation and microphysics depending on aerosol properties and meteorological conditions. The measurement results suggest that CCN activation and droplet formation in convective clouds are mostly aerosol-limited under PR and LRT conditions and updraft-limited under BB and MPOL conditions. Normalized CCN efficiency spectra (CCN divided by aerosol number concentration plotted against water vapor supersaturation) and corresponding parameterizations (Gaussian error function fits) provide a basis for further analysis and model studies of aerosol–cloud interactions in the Amazon.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

A variance-based decomposition and global sensitivity index method for uncertainty quantification: Application to retrieved ice cloud properties

This study develops a novel uncertainty quantification (UQ) method for cloud microphysical property retrievals using variance-based decomposition and global sensitivity index. In this UQ framework, empirical orthogonal function (EOF) analysis is applied to the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) ground-based observations, which are the inputs for the cloud retrieval studied here. The principal components (PCs) in the EOF expansion are parameterized as random input variables, and hence, the input dimension is greatly reduced (up to a factor of 50), allowing large ensemble of random samplings. The EOF expansion improves the accuracy of the uncertainty estimation by taking into account the cross correlations in the input data profiles. Here, this method enables a probabilistic representation of a retrieval process by adding normally distributed perturbations into PCs of sample means of input data profiles within a time window. Therefore, it effectively facilitates objective validation of climate models against cloud retrievals under a probabilistic framework for rigorous statistical inferences. Moreover, the variance-based global sensitivity index analysis, part of this method, attributes the output uncertainties to each individual source, thus providing directions for improving retrieval algorithms and observation strategies. For demonstration, we apply this method to quantify the uncertainties of the ARM program's baseline cloud retrieval algorithm for an ice cloud case observed at the Southern Great Plains site on 9 March 2000.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Assessing impacts of PBL and surface layer schemes in simulating the surface–atmosphere interactions and precipitation over the tropical ocean using observations from AMIE/DYNAMO

Accuracy of turbulence parameterization in representing Planetary Boundary Layer (PBL) processes in climate models is critical for predicting the initiation and development of clouds, air quality issues, and underlying surface-atmosphere-cloud interactions. In this study, we 1) evaluate WRF model-simulated spatial patterns of precipitation and surface fluxes, as well as vertical profiles of potential temperature, humidity, moist static energy and moisture tendency terms as simulated by WRF at various spatial resolutions and with PBL, surface layer and shallow convection schemes against measurements, 2) identify model biases by examining the moisture tendency terms contributed by PBL and convection processes through nudging experiments, and 3) evaluate the dependence of modeled surface latent heat (LH) fluxes onPBL and surface layer schemes over the tropical ocean. The results show that PBL and surface parameterizations have surprisingly large impacts on precipitation, convection initiation and surface moisture fluxes over tropical oceans. All of the parameterizations tested tend to overpredict moisture in PBL and free atmosphere, and consequently result in larger moist static energy and precipitation. Moisture nudging tends to suppress the initiation of convection and reduces the excess precipitation. The reduction in precipitation bias in turn reduces the surface wind and LH flux biases, which suggests that the model drifts at least partly because of a positive feedback between precipitation and surface fluxes. The updated shallow convection scheme KF-CuP tends to suppress the initiation and development of deep convection, consequently decreasing precipitation. The Eta surface layer scheme predicts more reasonable LH fluxes and the LH-Wind Speed relationship than the MM5 scheme, especially when coupled with the MYJ scheme. By examining various parameterization schemes in WRF, we identify sources of biases and weaknesses of current PBL, surface layer and shallow convection schemes in reproducing PBL processes, the initiation of convection and intra-seasonal variability of precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Mesoscale Convective Systems in a Superparameterized E3SM Simulation at High Resolution

Accurately representing mesoscale convective systems (MCSs) is crucial to simulating the energy and water cycles in global climate models. Using a novel MCS identification and tracking algorithm applied to observations and model simulations, we evaluate how well the Energy Exascale Earth System Model simulates MCSs over the central US. Simulations performed by E3SM at 25 km grid spacing with (SP-E3SM) and without (E3SM) superparameterization are compared and evaluated against observations using multiple metrics highlighting important MCS characteristics. Compared to E3SM, SP-E3SM better simulates MCS number and MCS precipitation amount, diurnal cycle, propagation, and the probability distribution of precipitation rate in both spring and summer. The improvement from SP is partly contributed by improvement in simulating the large-scale environments, featuring enhanced atmospheric low-level moisture and larger moisture transport to the central US relative to E3SM. However, SP-E3SM still underestimates MCS precipitation amount, particularly in summer. This underestimation is closely related to the negative bias in MCS precipitation intensity, although the drier environments simulated during summer also contributes to underestimation of MCS frequency. Without SP, the larger bias in MCS precipitation amount is closely related to the negative bias in MCS frequency, while the substantial dry bias in the large-scale environments also contributes to underestimation of MCS intensity. Our results suggest that SP improves MCS simulation by improving modeling of the large-scale environments and convection initiation, which are both major limiting factors in E3SM even at 25 km grid spacing where deep convection is represented by a cumulus parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Linking Large-Scale Double-ITCZ Bias to Local-Scale Drizzling Bias in Climate Models

Tropical precipitation in climate models presents significant biases in both the large-scale pattern (i.e., double intertropical convergence zone bias) and local-scale characteristics (i.e., drizzling bias with too frequent drizzle/convection and reduced occurrences of no and heavy precipitation). By untangling the coupled system and analyzing the biases in precipitation, cloud, and radiation, this study shows that local-scale drizzling bias in atmospheric models can lead to large-scale double-ITCZ bias in coupled models by inducing convective-regime-dependent biases in precipitation and cloud radiative effects (CRE). The double-ITCZ bias consists of a hemispherically asymmetric component that arises from the asymmetric SST bias and a nearly symmetric component that exists in atmospheric models without the SST bias. By increasing light rain but reducing heavy rain, local-scale drizzling bias induces positive (negative) precipitation bias in the moderate (strong) convective regime, leading to the nearly symmetric wet bias in atmospheric models. By affecting the cloud profile, local-scale drizzling bias induces positive (negative) CRE bias in the stratocumulus (convective) regime in atmospheric models. Because the stratocumulus (convective) region is climatologically more pronounced in the southern (northern) tropics, the CRE bias is deemed to be hemispherically asymmetric and drives warm and wet (cold and dry) biases in the southern (northern) tropics when coupled to ocean. Our results suggest that correcting local-scale drizzling bias is critical for fixing large-scale double-ITCZ bias. Finally, the drizzling and double-ITCZ biases are not alleviated in models with mesoscale (0.25°–0.5°) or even storm-resolving (~3 km) resolution, implying that either large-eddy simulation or fundamental improvement in small-scale subgrid parameterizations is needed.

54 ENVIRONMENTAL SCIENCES↗