Engineering topics

Ma, Po-Lun (ORCID:0000000331095316)

Publications and source records attributed to Ma, Po-Lun (ORCID:0000000331095316).

16 records

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 OSTI2024

Droplet collection efficiencies inferred from satellite retrievals constrain effective radiative forcing of aerosol–cloud interactions

Abstract. Process-oriented observational constraints for the anthropogenic effective radiative forcing due to aerosol–cloud interactions (ERFaci) are highly desirable because the uncertainty associated with ERFaci poses a significant challenge to climate prediction. The contoured frequency by optical depth diagram (CFODD) analysis supports the evaluation of model representation of cloud liquid-to-rain conversion processes because the slope of a CFODD, generated from joint MODerate Resolution Imaging Spectroradiometer (MODIS)-CloudSat cloud retrievals, provides an estimate of cloud droplet collection efficiency in single-layer warm liquid clouds. Here, we present an updated CFODD analysis as an observational constraint on the ERFaci due to warm rain processes and apply it to the U.S. Department of Energy's Energy Exascale Earth System Model version 2 (E3SMv2). A series of sensitivity experiments shows that E3SMv2 droplet collection efficiencies and ERFaci are highly sensitive to autoconversion, i.e., the rate of mass transfer from cloud liquid to rain, yielding a strong correlation between the CFODD slope and the shortwave component of ERFaci (ERFaciSW; Pearson's R=-0.91). E3SMv2's CFODD slope (0.20 ± 0.04) is in agreement with observations (0.20 ± 0.03). The strong sensitivity of ERFaciSW to the CFODD slope provides a useful constraint on highly uncertain warm rain processes, whereby ERFaciSW, constrained by MODIS-CloudSat, is estimated by calculating the intercept of the linear association between the ERFaciSW and the CFODD slopes, using the MODIS-CloudSat CFODD slope as a reference.

Beall, Charlotte M. (ORCID:0000000311370835)↗

DOE OSTI2024

Impacts of spatial heterogeneity of anthropogenic aerosol emissions in a regionally refined global aerosol–climate model

Abstract. Emissions of anthropogenic aerosol and their precursors are often prescribed in global aerosol models. Most of these emissions are spatially heterogeneous at model grid scales. When remapped from low-resolution data, the spatial heterogeneity in emissions can be lost, leading to large errors in the simulation. It can also cause the conservation problem if non-conservative remapping is used. The default anthropogenic emission treatment in the Energy Exascale Earth System Model (E3SM) is subject to both problems. In this study, we introduce a revised emission treatment for the E3SM Atmosphere Model (EAM) that ensures conservation of mass fluxes and preserves the original emission heterogeneity at the model-resolved grid scale. We assess the error estimates associated with the default emission treatment and the impact of improved heterogeneity and mass conservation in both globally uniform standard-resolution (∼ 165 km) and regionally refined high-resolution (∼ 42 km) simulations. The default treatment incurs significant errors near the surface, particularly over sharp emission gradient zones. Much larger errors are observed in high-resolution simulations. It substantially underestimates the aerosol burden, surface concentration, and aerosol sources over highly polluted regions, while it overestimates these quantities over less-polluted adjacent areas. Large errors can persist at higher elevation for daily mean estimates, which can affect aerosol extinction profiles and aerosol optical depth (AOD). We find that the revised treatment significantly improves the accuracy of the aerosol emissions from surface and elevated sources near sharp spatial gradient regions, with significant improvement in the spatial heterogeneity and variability of simulated surface concentration in high-resolution simulations. In the next-generation E3SM running at convection-permitting scales where the resolved spatial heterogeneity is significantly increased, the revised emission treatment is expected to better represent the aerosol emissions as well as their lifecycle and impacts on climate.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Assessing the sensitivity of aerosol mass budget and effective radiative forcing to horizontal grid spacing in E3SMv1 using a regional refinement approach

Abstract. Atmospheric aerosols have important impacts on air quality and the Earth–atmospheric energy balance. However, as computing power is limited, Earth system models generally use coarse spatial grids and parameterize finer-scale atmospheric processes. These parameterizations and the simulation of atmospheric aerosols are often sensitive to model horizontal resolutions. Understanding the sensitivities is necessary for the development of Earth system models at higher resolutions with the deployment of more powerful supercomputers. Using the Energy Exascale Earth System Model (E3SM) version 1, this study investigates the impact of horizontal grid spacing on the simulated aerosol mass budget, aerosol–cloud interactions, and the effective radiative forcing of anthropogenic aerosols (ERFaer) over the contiguous United States. We examine the resolution sensitivity by comparing the nudged simulation results for 2016 from the low-resolution model (LR) and the regional refinement model (RRM). As expected, the simulated emissions of natural dust, sea salt, and marine organic matter are substantially higher in the RRM than in the LR. In addition, RRM simulates stronger aqueous-phase production of sulfate through the enhanced oxidation of sulfur dioxide by hydrogen peroxide due to increased cloud liquid water content. In contrast, the gas-phase chemical production of sulfate is slightly suppressed. The RRM resolves more large-scale precipitation and produces less convective precipitation than the LR, leading to increased (decreased) aerosol wet scavenging by large-scale (convective) precipitation. Regarding aerosol effects on clouds, RRM produces larger temporal variabilities in the large-scale liquid cloud fractions than LR, resulting in increased microphysical cloud processing of aerosols (more interstitial aerosols are converted to cloud-borne aerosols via aerosol activation) in RRM. Water vapor condensation is also enhanced in RRM compared to LR. Consequently, the RRM simulation produces more cloud droplets, a larger cloud droplet radius, a higher liquid water path, and a larger cloud optical depth than the LR simulation. A comparison of the present-day and pre-industrial simulations indicates that, for this contiguous United States domain, the higher-resolution increases ERFaer at the top of the model by about 12 %, which is mainly attributed to the strengthened indirect effect associated with aerosol–cloud interactions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Understanding changes in cloud simulations from E3SM version 1 to version 2

Abstract. This study documents clouds simulated by the Energy Exascale Earth System Model (E3SM) version 2 (E3SMv2) and attempts to understand what causes the model behavior change in clouds relative to E3SMv1. This is done by analyzing the last 30-year (1985–2014) data from the 165-year historical simulations using E3SMv1 and v2 and four sensitivity tests to isolate the impact of changes in model parameter choices in its turbulence, shallow convection, and cloud macrophysics parameterization (Cloud Layers Unified By Binormals, CLUBB); microphysical parameterization (MG2); and deep-convection scheme (ZM), as well as model physics changes in convective triggering. It is shown that E3SMv2 significantly improves the simulation of subtropical coastal stratocumulus clouds and clouds with optical depth larger than 3.6 over the stratocumulus-to-cumulus transition regimes, where the shortwave cloud radiative effect (SWCRE) is also improved, and the Southern Ocean (SO) while seeing an overall slight degradation in low clouds over other tropical and subtropical oceans. The better performance in E3SMv1 over those regions is partially due to error compensation between its simulated optically thin and intermediate low clouds for which E3SMv2 actually improves simulation of optically intermediate low clouds. Sensitivity tests indicate that the changes in low clouds are primarily due to the tuning done in CLUBB. The impact of the ZM tuning is mainly on optically intermediate and thick high clouds, contributing to an improved SWCRE and longwave cloud radiative effect (LWCRE). The impact of the MG2 tuning and the new convective trigger is primarily on the high latitudes and the SO. They have a relatively smaller impact on clouds than CLUBB tuning and ZM tuning do. This study offers additional insights into clouds simulated in E3SMv2 by utilizing multiple data sets and the Cloud Feedback Model Intercomparison Project (CFMIP) Observation Simulator Package (COSP) diagnostic tool as well as sensitivity tests. The improved understanding will benefit future E3SM developments.

58 GEOSCIENCES↗

DOE OSTI2023

Earth System Model Aerosol–Cloud Diagnostics (ESMAC Diags) package, version 2: assessing aerosols, clouds, and aerosol–cloud interactions via field campaign and long-term observations

Abstract. Poor representations of aerosols, clouds, and aerosol–cloud interactions (ACIs) in Earth system models (ESMs) have long been the largest uncertainties in predicting global climate change. Huge efforts have been made to improve the representation of these processes in ESMs, and the key to these efforts is the evaluation of ESM simulations with observations. Most well-established ESM diagnostics packages focus on the climatological features; however, they lack process-level understanding and representations of aerosols, clouds, and ACIs. In this study, we developed the Earth System Model Aerosol–Cloud Diagnostics (ESMAC Diags) package to facilitate the routine evaluation of aerosols, clouds, and ACIs simulated the Energy Exascale Earth System Model (E3SM) from the US Department of Energy (DOE). This paper documents its version 2 functionality (ESMAC Diags v2), which has substantial updates compared with version 1 (Tang et al., 2022a). The simulated aerosol and cloud properties have been extensively compared with in situ and remote-sensing measurements from aircraft, ship, surface, and satellite platforms in ESMAC Diags v2. It currently includes six field campaigns and two permanent sites covering four geographical regions: the eastern North Atlantic, the central US, the northeastern Pacific, and the Southern Ocean. These regions produce frequent liquid- or mixed-phase clouds, with extensive measurements available from the DOE Atmospheric Radiation Measurement user facility and other agencies. ESMAC Diags v2 generates various types of single-variable and multivariable diagnostics, including percentiles, histograms, joint histograms, and heatmaps, to evaluate the model representation of aerosols, clouds, and ACIs. Select examples highlighting the capabilities of ESMAC Diags are shown using E3SM version 2 (E3SMv2). In general, E3SMv2 can reasonably reproduce many observed aerosol and cloud properties, with biases in some variables such as aerosol particle and cloud droplet sizes and number concentrations. The coupling of aerosol and cloud number concentrations may be too strong in E3SMv2, possibly indicating a bias in processes that control aerosol activation. Furthermore, the liquid water path response to a perturbed cloud droplet number concentration behaves differently in E3SMv2 and observations, which warrants further study to improve the cloud microphysics parameterizations in E3SMv2.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Long-term variability in immersion-mode marine ice-nucleating particles from climate model simulations and observations

Abstract. Ice-nucleating particles (INPs) in the Southern Ocean (SO) atmosphere have significant impacts on cloud radiative and microphysical properties. Yet, INP prediction skill in climate models remains poorly understood, in part because of the lack of long-term measurements. Here we show, for the first time, how model-simulated INP concentrations compare with year-round INP measurements during the Macquarie Island Cloud Radiation Experiment (MICRE) campaign from 2017–2018. We simulate immersion-mode INP concentrations using the Energy Exascale Earth System Model version 1 (E3SMv1) by combining simulated aerosols with recently developed deterministic INP parameterizations and the native classical nucleation theory (CNT) for mineral dust in E3SMv1. Because MICRE did not collect aerosol measurements of super-micron particles, which are more effective ice nucleators, we evaluate the model's aerosol fields at other high-latitude sites using long-term in situ observations of dust and sea spray aerosol. We find that the model underestimates dust and overestimates sea spray aerosol concentrations by 1 to 2 orders of magnitude for most of the high-latitude sites in the Southern Hemisphere. We next compare predicted INP concentrations with concentrations of INPs collected on filter samples (typically for 2 or 3 d) and processed offline using the Colorado State University ice spectrometer (IS) in immersion freezing mode. We find that when deterministic parameterizations for both dust and sea spray INPs are used, simulated INPs are within a factor of 10 of observed INPs more than 60 % of the time during summer. Our results also indicate that the E3SM's current treatment of mineral dust immersion freezing in the SO is impacted by compensating biases – an underprediction of dust amount was compensated by an overprediction of its effectiveness as INPs. We also perform idealized droplet freezing experiments to quantify the implications of the time-dependent behavior assumed by the E3SM's CNT-parameterization and compare with the ice spectrometer observations. We find that the E3SM CNT 10 s diagnostic used in this study is a reasonable approximation of the exact formulation of CNT, when applied to ice spectrometer measurements in low-INP conditions similar to Macquarie Island. However, the linearized 10 s diagnostic underestimates the exact formula by an order of magnitude or more in places with high-INP conditions like the Sahara. Overall, our findings suggest that it is important to correct the biases in E3SM's simulated dust life cycle and update E3SM's INP parameterizations. INP prediction errors of 2 to 3 orders of magnitude can have considerable impacts on the simulated cloud and radiative properties in global climate models. On comparing INP concentrations during MICRE against ship-based campaigns, Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS) and Antarctic Circumnavigation Expedition (ACE), we find that INPs from the latter are significantly higher only in regions closer to Macquarie Island. This alludes to the fact that physical, chemical and biological processes affecting INP concentrations as stimulated by the island could be partly responsible for the high INP concentrations observed at Macquarie Island during the MICRE campaign. Therefore, improvements to both aerosol simulation and INP parameterizations are required to adequately simulate INPs and their cloud impacts in E3SM. It will be helpful to include a parallel measurement of the size-resolved aerosol composition and explore opportunities for long-term measurement platforms in future field campaigns studying INP sources in remote marine regions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Emulating aerosol optics with randomly generated neural networks

Abstract. Atmospheric aerosols have a substantial impact on climate and remain one of the largest sources of uncertainty in climate prediction. Accurate representation of their direct radiative effects is a crucial component of modern climate models. However, direct computation of the radiative properties of aerosol populations is far too computationally expensive to perform in a climate model, so optical properties are typically approximated using a parameterization. This work develops artificial neural networks (ANNs) capable of replacing the current aerosol optics parameterization used in the Energy Exascale Earth System Model (E3SM). A large training dataset is generated by using Mie code to directly compute the optical properties of a range of atmospheric aerosol populations given a large variety of particle sizes, wavelengths, and refractive indices. Optimal neural architectures for shortwave and longwave bands are identified by evaluating ANNs with randomly generated wirings. Randomly generated deep ANNs are able to outperform conventional multilayer-perceptron-style architectures with comparable parameter counts. Finally, the ANN-based parameterization produces significantly more accurate bulk aerosol optical properties than the current parameterization when compared with direct Mie calculations using mean absolute error. The success of this approach makes possible the future inclusion of much more sophisticated representations of aerosol optics in climate models that cannot be captured by extension of the existing parameterization scheme and also demonstrates the potential of random-wiring-based neural architecture search in future applications in the Earth sciences.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Incorporation of aerosol into the COSPv2 satellite lidar simulator for climate model evaluation

Abstract. Atmospheric aerosol has substantial impacts on climate, air quality and biogeochemical cycles, and its concentrations are highly variable in space and time. A key variability to evaluate within models that simulate aerosol is the vertical distribution, which influences atmospheric heating profiles and aerosol–cloud interactions, to help constrain aerosol residence time and to better represent the magnitude of simulated impacts. To ensure a consistent comparison between modeled and observed vertical distribution of aerosol, we implemented an aerosol lidar simulator within the Cloud Feedback Model Intercomparison Project (CFMIP) Observation Simulator Package version 2 (COSPv2). We assessed the attenuated total backscattered (ATB) signal and the backscatter ratios (SRs) at 532 nm in the U.S. Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1). The simulator performs the computations at the same vertical resolution as the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP), making use of aerosol optics from the E3SMv1 model as inputs and assuming that aerosol is uniformly distributed horizontally within each model grid box. The simulator applies a cloud masking and an aerosol detection threshold to obtain the ATB and SR profiles that would be observed above clouds by CALIOP with its aerosol detection capability. Our analysis shows that the aerosol distribution simulated at a seasonal timescale is generally in good agreement with observations. Over the Southern Ocean, however, the model does not produce the SR maximum as observed in the real world. Comparison between clear-sky and all-sky SRs shows little differences, indicating that the cloud screening by potentially incorrect model clouds does not affect the mean aerosol signal averaged over a season. This indicates that the differences between observed and simulated SR values are due not to sampling errors, but to deficiencies in the representation of aerosol in models. Finally, we highlight the need for future applications of lidar observations at multiple wavelengths to provide insights into aerosol properties and distribution and their representation in Earth system models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Representing surface heterogeneity in land–atmosphere coupling in E3SMv1 single-column model over ARM SGP during summertime

Abstract. The Earth's land surface features spatial and temporal heterogeneity over a wide range of scales below those resolved by current Earth system models (ESMs). State-of-the-art land and atmosphere models employ parameterizations to represent their subgrid heterogeneity, but the land–atmosphere coupling in ESMs typically operates on the grid scale. Communicating the information on the land surface heterogeneity with the overlying atmospheric boundary layer (ABL) remains a challenge in modeling land–atmosphere interactions. In order to account for the subgrid-scale heterogeneity in land–atmosphere coupling, we implement a new coupling scheme in the Energy Exascale Earth system model version 1 (E3SMv1) that uses adjusted surface variances and covariance of potential temperature and specific water content as the lower boundary condition for the atmosphere model. The new lower boundary condition accounts for both the variability of individual subgrid land surface patches and the inter-patch variability. The E3SMv1 single-column model (SCM) simulations over the Atmospheric Radiation Measurement (ARM) Southern Great Plain (SGP) site were performed to assess the impacts. We find that the new coupling parameterization increases the magnitude and diurnal cycle of the temperature variance and humidity variance in the lower ABL on non-precipitating days. The impacts are primarily attributed to subgrid inter-patch variability rather than the variability of individual patches. These effects extend vertically from the surface to several levels in the lower ABL on clear days. We also find that accounting for surface heterogeneity increases low cloud cover and liquid water path (LWP). These cloud changes are associated with the change in cloud regime indicated by the skewness of the probability density function (PDF) of the subgrid vertical velocity. In precipitating days, the inter-patch variability reduces significantly so that the impact of accounting for surface heterogeneity vanishes. These results highlight the importance of accounting for subgrid heterogeneity in land–atmosphere coupling in next-generation ESMs.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Effective radiative forcing of anthropogenic aerosols in E3SM version 1: historical changes, causality, decomposition, and parameterization sensitivities

Abstract. The effective radiative forcing of anthropogenic aerosols (ERFaer) is an important measure of the anthropogenic aerosol effects simulated by a global climate model. Here we analyze ERFaer simulated by the E3SM version 1 (E3SMv1) atmospheric model using both century-long free-running atmosphere–land simulations and short nudged simulations. We relate the simulated ERFaer to characteristics of the aerosol composition and optical properties, and we evaluate the relationships between key aerosol and cloud properties. In terms of historical changes from the year 1870 to 2014, our results show that the global mean anthropogenic aerosol burden and optical depth increase during the simulation period as expected, but the regional averages show large differences in the temporal evolution. The largest regional differences are found in the emission-induced evolution of the burden and optical depth of the sulfate aerosol: a strong decreasing trend is seen in the Northern Hemisphere high-latitude region after around 1970, while a continued increase is simulated in the tropics. The relationships between key aerosol and cloud properties (relative changes between pre-industrial and present-day conditions) also show evident changes after 1970, diverging from the linear relationships exhibited for the period of 1870–1969. In addition to the regional differences in the simulated relationships, a reduced sensitivity in cloud droplet number and other cloud properties to aerosol perturbations is seen when the aerosol perturbation is large. Consequently, the global annual mean ERFaer magnitude does not increase after around 1970. The ERFaer in E3SMv1 is relatively large compared to the recently published multi-model estimates; the primary reason is the large indirect aerosol effect (i.e., through aerosol–cloud interactions). Compared to other models, E3SMv1 features large relative changes in the cloud droplet effective radius in response to aerosol perturbations. Large sensitivity is also seen in the liquid cloud optical depth, which is determined by changes in both the effective radius and liquid water path. Aerosol-induced changes in liquid and ice cloud properties in E3SMv1 are found to have a strong correlation, as the evolution of anthropogenic sulfate aerosols affects both the liquid cloud formation and the homogeneous ice nucleation in cirrus clouds (that causes a large effect on longwave ERFaer). As suggested by a previous study, the large ERFaer appears to be one of the reasons why the model cannot reproduce the observed global mean temperature evolution in the second half of the 20th century. Sensitivity simulations are performed to understand which parameterization and/or parameter changes have a large impact on the simulated ERFaer. The ERFaer estimates in E3SMv1 for the shortwave and longwave components are sensitive to the parameterization changes in both liquid and ice cloud processes. When the parameterization of ice cloud processes is modified, the top-of-model forcing changes in the shortwave and longwave components largely offset each other, so the net effect is negligible. This suggests that, to reduce the magnitude of the net ERFaer, it would be more effective to reduce the anthropogenic aerosol effect through liquid or mixed-phase clouds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Earth System Model Aerosol–Cloud Diagnostics (ESMAC Diags) package, version 1: assessing E3SM aerosol predictions using aircraft, ship, and surface measurements

Abstract. An Earth system model (ESM) aerosol–cloud diagnostics package is developed to facilitate the routine evaluation of aerosols, clouds, and aerosol–cloud interactions simulated by the Energy Exascale Earth System Model (E3SM) from the US Department of Energy (DOE). The first version focuses on comparing simulated aerosol properties with aircraft, ship, and surface measurements, which are mostly measured in situ. The diagnostics currently cover six field campaigns in four geographical regions: eastern North Atlantic (ENA), central US (CUS), northeastern Pacific (NEP), and Southern Ocean (SO). These regions produce frequent liquid- or mixed-phase clouds, with extensive measurements available from the Atmospheric Radiation Measurement (ARM) program and other agencies. Various types of diagnostics and metrics are performed for aerosol number, size distribution, chemical composition, cloud condensation nuclei (CCN) concentration, and various meteorological quantities to assess how well E3SM represents observed aerosol properties across spatial scales. Overall, E3SM qualitatively reproduces the observed aerosol number concentration, size distribution, and chemical composition reasonably well, but it overestimates Aitken-mode aerosols and underestimates accumulation-mode aerosols over the CUS and ENA regions, suggesting that processes related to particle growth or coagulation might be too weak in the model. The current version of E3SM struggles to reproduce the new particle formation events frequently observed over both the CUS and ENA regions, indicating missing processes in current parameterizations. The diagnostics package is coded and organized in a way that can be extended to other field campaign datasets and adapted to higher-resolution model simulations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Better calibration of cloud parameterizations and subgrid effects increases the fidelity of the E3SM Atmosphere Model version 1

Abstract. Realistic simulation of the Earth's mean-state climate remains a major challenge, and yet it is crucial for predicting the climate system in transition. Deficiencies in models' process representations, propagation of errors from one process to another, and associated compensating errors can often confound the interpretation and improvement of model simulations. These errors and biases can also lead to unrealistic climate projections and incorrect attribution of the physical mechanisms governing past and future climate change. Here we show that a significantly improved global atmospheric simulation can be achieved by focusing on the realism of process assumptions in cloud calibration and subgrid effects using the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1). The calibration of clouds and subgrid effects informed by our understanding of physical mechanisms leads to significant improvements in clouds and precipitation climatology, reducing common and long-standing biases across cloud regimes in the model. The improved cloud fidelity in turn reduces biases in other aspects of the system. Furthermore, even though the recalibration does not change the global mean aerosol and total anthropogenic effective radiative forcings (ERFs), the sensitivity of clouds, precipitation, and surface temperature to aerosol perturbations is significantly reduced. This suggests that it is possible to achieve improvements to the historical evolution of surface temperature over EAMv1 and that precise knowledge of global mean ERFs is not enough to constrain historical or future climate change. Cloud feedbacks are also significantly reduced in the recalibrated model, suggesting that there would be a lower climate sensitivity when it is run as part of the fully coupled E3SM. This study also compares results from incremental changes to cloud microphysics, turbulent mixing, deep convection, and subgrid effects to understand how assumptions in the representation of these processes affect different aspects of the simulated atmosphere as well as its response to forcings. We conclude that the spectral composition and geographical distribution of the ERFs and cloud feedback, as well as the fidelity of the simulated base climate state, are important for constraining the climate in the past and future.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Quantifying the structural uncertainty of the aerosol mixing state representation in a modal model

Abstract. Aerosol mixing state is an important emergent property that affects aerosol radiative forcing and aerosol–cloud interactions, but it has not been easy to constrain this property globally. This study aims to verify the global distribution of aerosol mixing state represented by modal models. To quantify the aerosol mixing state, we used the aerosol mixing state indices for submicron aerosol based on the mixing of optically absorbing and non-absorbing species (χo), the mixing of primary carbonaceous and non-primary carbonaceous species (χc), and the mixing of hygroscopic and non-hygroscopic species (χh). To achieve a spatiotemporal comparison, we calculated the mixing state indices using output from the Community Earth System Model with the four-mode version of the Modal Aerosol Module (MAM4) and compared the results with the mixing state indices from a benchmark machine-learned model trained on high-detail particle-resolved simulations from the particle-resolved stochastic aerosol model PartMC-MOSAIC. The two methods yielded very different spatial patterns of the mixing state indices. In some regions, the yearly averaged χ value computed by the MAM4 model differed by up to 70 percentage points from the benchmark values. These errors tended to be zonally structured, with the MAM4 model predicting a more internally mixed aerosol at low latitudes and a more externally mixed aerosol at high latitudes compared to the benchmark. Our study quantifies potential model bias in simulating mixing state in different regions and provides insights into potential improvements to model process representation for a more realistic simulation of aerosols towards better quantification of radiative forcing and aerosol–cloud interactions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

On the estimation of boundary layer heights: a machine learning approach

Abstract. The planetary boundary layer height (zi) is a key parameter used in atmospheric models for estimating the exchange of heat, momentum, and moisture between the surface and the free troposphere. Near-surface atmospheric and subsurface properties (such as soil temperature, relative humidity, etc.) are known to have an impact on zi. Nevertheless, precise relationships between these surface properties and zi are less well known and not easily discernible from the multi-year dataset. Machine learning approaches, such as random forest (RF), which use a multi-regression framework, help to decipher some of the physical processes linking surface-based characteristics to zi. In this study, a 4-year dataset from 2016 to 2019 at the Southern Great Plains site is used to develop and test a machine learning framework for estimating zi. Parameters derived from Doppler lidars are used in combination with over 20 different surface meteorological measurements as inputs to a RF model. The model is trained using radiosonde-derived zi values spanning the period from 2016 through 2018 and then evaluated using data from 2019. Results from 2019 showed significantly better agreement with the radiosonde compared to estimates derived from a thresholding technique using Doppler lidars only. Noteworthy improvements in daytime zi estimates were observed using the RF model, with a 50 % improvement in mean absolute error and an R2 of greater than 85 % compared to the Tucker method zi. We also explore the effect of zi uncertainty on convective velocity scaling and present preliminary comparisons between the RF model and zi estimates derived from atmospheric models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Physically regularized machine learning emulators of aerosol activation

Abstract. The activation of aerosol into cloud droplets is an important step in the formation of clouds and strongly influences the radiative budget of the Earth. Explicitly simulating aerosol activation in Earth system models is challenging due to the computational complexity required to resolve the necessary chemical and physical processes and their interactions. As such, various parameterizations have been developed to approximate these details at reduced computational cost and accuracy. Here, we explore how machine learning emulators can be used to bridge this gap in computational cost and parameterization accuracy. We evaluate a set of emulators of a detailed cloud parcel model using physically regularized machine learning regression techniques. We find that the emulators can reproduce the parcel model at higher accuracy than many existing parameterizations. Furthermore, physical regularization tends to improve emulator accuracy, most significantly when emulating very low activation fractions. This work demonstrates the value of physical constraints in machine learning model development and enables the implementation of improved hybrid physical and machine learning models of aerosol activation into next-generation Earth system models.

58 GEOSCIENCES↗