Engineering topics

Shi, Xiangjun

Publications and source records attributed to Shi, Xiangjun.

11 records

DOE OSTI2022

Development of a triple-moment ice-phase cloud microphysics scheme and its application to the Single Column Atmosphere Model

Parameterization of cloud microphysics is critical for accurate simulation of weather and climate, in which the characteristics of cloud particle spectrum inevitably further affect climate simulation by changing the cloud evolution and cloud radiation effects. The popular currently used double-moment cloud microphysics schemes in numerical models can predict the intercept (N 0 ) and slope (λ) parameters of cloud particle spectrum but cannot predict the spectral shape parameter (μ), which hinders accurate description of cloud physical processes in climate models. Therefore, in the present study, we built upon the ideas of previously developed triple-moment cloud microphysics scheme, considered radar reflectivity factor as the third predictor in addition to number and mass concentration, and introduced its related prediction equation to the Single Column Atmosphere Model Version 5.3 (SCAM5.3). Moreover, the relevant microphysical process formulas in the model were revised, and a triple-moment ice-cloud microphysics scheme was constructed to predict the μ of ice particle spectrum. Based on this model, a 29-day case of the Atmospheric Radiation Measurement Program in the summer of 1997 (ARM97) was simulated, and the differences in cloud fraction and radiation simulations between the double- and triple-moment schemes were analysed. The μ of ice particle spectrum predicted by the triple-moment scheme mainly ranged from 0 to 4; the peak value was around 2, and at least 85% of the μ values were greater than the default 0. Therefore, the μ=0 setting in the double-moment scheme is unreasonable. Moreover, the developed triple-moment ice-cloud microphysics scheme yielded a narrower ice particle spectrum, which was closer to the observation results of previous studies, than the double-moment cloud microphysics scheme (e.g., height=233 hPa). Furthermore, compared with the double-moment scheme, the triple-moment scheme achieved closer simulations of the cloud fraction observations, particularly for ice clouds in the upper levels. Considering the close association between cloud fraction and radiation, the error with the triple-moment scheme was smaller than that with the double-moment scheme regardless of the shortwave (surface downward shortwave radiation, net downward shortwave radiation at the top of the atmosphere, and net shortwave radiation at the surface) or longwave (surface downward longwave radiation and net longwave radiation at the surface) flux density. Finally, the improvement mechanism of the triple-moment scheme on cloud fraction and radiation simulations was explored. The major reason is that the triple-moment scheme weakens the autoconversion of ice crystal to form snow and enhances the growth process of ice crystal deposition, thereby increasing the mass concentration of ice crystals and ice cloud fraction. Overall, the developed triple-moment ice-cloud microphysics scheme can help improve the simulation ability of models for cloud-climate feedback and other processes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

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

The effective radiative forcing of anthropogenic aerosols (ERF aer ) is an important measure of the anthropogenic aerosol effects simulated by a global climate model. Here we analyze ERF aer simulated by the E3SMv1 atmosphere model using both century-long free-running atmosphere-land simulations and short nudged simulations. We relate the simulated ERF aer to characteristics of the aerosol composition and optical properties, and 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. Consequently, although the global mean anthropogenic aerosol burden and optical depth increase from 1870 to 2014, the ERF aer magnitude does not increase after around year 1970. The relationships between key aerosol and cloud properties (relative changes between preindustrial and present-day conditions) also show evident changes after 1970, diverging from the linear relationships exhibited for the period from 1870 to 2014. The ERF aer 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 a stronger sensitivity of the cloud droplet effective radius to changes in the cloud droplet number concentration. 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. The ERF aer 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-atmosphere 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 ERF aer , it would be more effective to reduce the anthropogenic aerosol effect through liquid or mixed-phase clouds.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Analytical integration of the tractions induced by non-singular dislocations on an arbitrary shaped triangular quadratic element

Here, an analytical model is proposed to evaluate the nodal force induced by a segment of dislocation upon an arbitrary shaped triangular element. This calculation is required in hybrid methods that associate dislocation dynamics to boundary or finite element to solve simultaneously the evolution of large ensembles of dislocations with complex boundary conditions. Nodal forces are defined as the triple integration of the unbalanced traction field induced by a straight dislocation upon the surface of the element. Following our previous approach (Queyreau et al 2014 Modelling Simul. Mater. Sci. Eng. 22 035004) on a simpler geometry and in the case of linear isotropic elasticity, triple integrals are solved by sequences of integration by parts that exhibit recurrence relations. The traction field is defined and finite everywhere even at the core of dislocations, thanks to the use of the non-singular stress expression formulated by Cai et al (2006 J. Mech. Phys. Solids 54 561–587). The nodal force expressions can be used when considering both a single convolution or double convolution of the Green's function with the core distribution. A solution is also proposed for the case of a semi-infinite segment through the study of the asymptotic behavior of the analytical expressions. The proposed approach is exact and very computationally efficient. The choice of an arbitrary shaped triangular element and quadratic shape functions allow the consideration of complex geometries and comply with automatic meshing procedures. These analytical expressions could also be employed to estimate dislocation interactions with interfaces.

42 ENGINEERING↗

DOE OSTI2018

Impact of aerosols on ice crystal size

The interactions between aerosols and ice clouds represent one of the largest uncertainties in global radiative forcing from pre-industrial time to the present. In particular, the impact of aerosols on ice crystal effective radius (R ei ), which is a key parameter determining ice clouds' net radiative effect, is highly uncertain due to limited and conflicting observational evidence. Here we investigate the effects of aerosols on R ei under different meteorological conditions using 9-yearsatellite observations. We find that the responses of R ei to aerosol loadings are modulated by water vapor amount in conjunction with several other meteorological parameters. While there is a significant negative correlation between Rei and aerosol loading in moist conditions, consistent with the “Twomey effect” for liquid clouds, a strong positive correlation between the two occurs in dry conditions. Simulations based on a cloud parcel model suggest that water vapor modulates the relative importance of different ice nucleation modes, leading to the opposite aerosol impacts between moist and dry conditions. When ice clouds are decomposed into those generated from deep convection and formed in situ, the water vapor modulation remains in effect for both ice cloud types, although the sensitivities of R ei to aerosols differ noticeably between them due to distinct formation mechanisms. The water vapor modulation can largely explain the difference in the responses of R ei to aerosol loadings in various seasons. A proper representation of the water vapor modulation is essential for an accurate estimate of aerosol–cloud radiative forcing produced by ice clouds.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

Effects of pre-existing ice crystals on cirrus clouds and comparison between different ice nucleation parameterizations with the Community Atmosphere Model (CAM5)

In order to improve the treatment of ice nucleation in a more realistic manner in the Community Atmosphere Model version 5.3 (CAM5.3), the effects of pre-existing ice crystals on ice nucleation in cirrus clouds are considered. In addition, by considering the in-cloud variability in ice saturation ratio, homogeneous nucleation takes place spatially only in a portion of the cirrus cloud rather than in the whole area of the cirrus cloud. Compared to observations, the ice number concentrations and the probability distributions of ice number concentration are both improved with the updated treatment. The pre-existing ice crystals significantly reduce ice number concentrations in cirrus clouds, especially at mid- to high latitudes in the upper troposphere (by a factor of ~10). Furthermore, the contribution of heterogeneous ice nucleation to cirrus ice crystal number increases considerably. Besides the default ice nucleation parameterization of Liu and Penner (2005, hereafter LP) in CAM5.3, two other ice nucleation parameterizations of Barahona and Nenes (2009, hereafter BN) and Kärcher et al. (2006, hereafter KL) are implemented in CAM5.3 for the comparison. In-cloud ice crystal number concentration, percentage contribution from heterogeneous ice nucleation to total ice crystal number, and pre-existing ice effects simulated by the three ice nucleation parameterizations have similar patterns in the simulations with present-day aerosol emissions. However, the change (present-day minus pre-industrial times) in global annual mean column ice number concentration from the KL parameterization (3.24 × 10 6 m -2 ) is less than that from the LP (8.46 × 10 6 m -2 ) and BN (5.62 × 10 6 m -2 ) parameterizations. As a result, the experiment using the KL parameterization predicts a much smaller anthropogenic aerosol long-wave indirect forcing (0.24 W m -2 ) than that using the LP (0.46 W m −2 ) and BN (0.39 W m -2 ) parameterizations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2013

A Single Column Model Ensemble Approach Applied to the TWP-ICE Experiment

Single column models (SCM) are useful testbeds for investigating the parameterisation schemes of numerical weather prediction and climate models. The usefulness of SCM simulations are limited, however, by the accuracy of the best-estimate large-scale data prescribed. One method to address this uncertainty is to perform ensemble simulations of the SCM. This study first derives an ensemble of large-scale data for the Tropical Warm Pool International Cloud Experiment (TWP-ICE) based on an estimate of a possible source of error in the best-estimate product. This data is then used to carry out simulations with 11 SCM and 2 cloud-resolving models (CRM). Best-estimate simulations are also performed. All models show that moisture related variables are close to observations and there are limited differences between the best-estimate and ensemble mean values. The models, however, show different sensitivities to changes in the forcing particularly when weakly forced. The ensemble simulations highlight important differences in the moisture budget between the SCM and CRM. Systematic differences are also apparent in the ensemble mean vertical structure of cloud variables. The ensemble is further used to investigate relations between cloud variables and precipitation identifying large differences between CRM and SCM. This study highlights that additional information can be gained by performing ensemble simulations enhancing the information derived from models using the more traditional single best-estimate simulation.

applied↗

DOE OSTI2013

Two-moment Bulk Stratiform Cloud Microphysics in the Grid-point Atmospheric Model of IAP LASG (GAMIL)

A two-moment bulk stratiform microphysics scheme, including recently developed physically-based droplet activation/ice nucleation parameterizations has been implemented into the Grid-point Atmospheric Model of IAP LASG (GAMIL) as an effort to enhance the model capability for studying aerosol indirect effects. Unlike the previous one-moment cloud microphysics scheme, the new scheme produces reasonable representation of cloud particle size and number concentration. This scheme captures the observed spatial variations in cloud droplet number concentrations. Simulated ice crystal number concentrations in cirrus clouds qualitatively agree with in-situ observations. The longwave and shortwave cloud forcing are in better agreement with observations. Sensitivity tests show that the column cloud droplet number concentrations calculated from two different droplet activation parameterizations are similar. However, ice crystal number concentration in mixed-phased clouds is sensitive to different heterogeneous freezing formulations. The simulation with high ice crystal number concentration in mixed-phase clouds has less liquid water path and weaker cloud forcing. Furthermore, ice crystal number concentration in cirrus clouds is sensitive to different ice nucleation parameterizations. Sensitivity tests also suggest that impact of pre-existing ice crystals on homogeneous freezing in old clouds should be taken into account.

(GAMIL)↗

DOE OSTI2012

Sensitivity Studies of Dust Ice Nuclei Effect on Cirrus Clouds with the Community Atmosphere Model CAM5

In this study the effect of dust aerosol on upper tropospheric cirrus clouds through heterogeneous ice nucleation is investigated in the Community Atmospheric Model version 5 (CAM5) with two ice nucleation parameterizations. Both parameterizations consider homogeneous and heterogeneous nucleation and the competition between the two mechanisms in cirrus clouds, but differ significantly in the number concentration of heterogeneous ice nuclei (IN) from dust. Heterogeneous nucleation on dust aerosol reduces the occurrence frequency of homogeneous nucleation and thus the ice crystal number concentration in the Northern Hemisphere (NH) cirrus clouds compared to simulations with pure homogeneous nucleation. Global and annual mean shortwave and longwave cloud forcing are reduced by up to 2.0 ± 0.1 W m-2 (1σ uncertainty) and 2.4 ± 0.1 W m-2, respectively due to the presence of dust IN, with the net cloud forcing change of -0.40 ± 0.20 W m-2. Comparison of model simulations with in situ aircraft data obtained in NH mid-latitudes suggests that homogeneous ice nucleation may play an important role in the ice nucleation at these regions with temperatures of 205–230 K. However, simulations overestimate observed ice crystal number concentrations in the tropical tropopause regions with temperatures of 190–205 K, and overestimate the frequency of occurrence of high ice crystal number concentration (> 200 L-1) and underestimate the frequency of low ice crystal number concentration (< 30 L-1) at NH mid-latitudes. These results highlight the importance of quantifying the number concentrations and properties of heterogeneous IN (including dust aerosol) in the upper troposphere from the global perspective.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2012

Toward a Minimal Representation of Aerosols in Climate Models: Description and Evaluation in the Community Atmosphere Model CAM5

A modal aerosol module (MAM) has been developed for the Community Atmosphere Model version 5 (CAM5), the atmospheric component of the Community Earth System Model version 1 (CESM1). MAM is capable of simulating the aerosol size distribution and both internal and external mixing between aerosol components, treating numerous complicated aerosol processes and aerosol physical, chemical and optical properties in a physically based manner. Two MAM versions were developed: a more complete version with seven-lognormal modes (MAM7), and a three-lognormal mode version (MAM3) for the purpose of long-term (decades to centuries) simulations. Major approximations in MAM3 include assuming immediate mixing of primary organic matter (POM) and black carbon (BC) with other aerosol components, merging of the MAM7 fine dust and fine sea salt modes into the accumulation mode, merging of the MAM7 coarse dust and coarse sea salt modes into the single coarse mode, and neglecting the explicit treatment of ammonia and ammonium cycles. Simulated sulfate and secondary organic aerosol (SOA) mass concentrations are remarkably similar between MAM3 and MAM7 as most ({approx}90%) of these aerosol species are in the accumulation mode. Differences of POM and BC concentrations between MAM3 and MAM7 are also small (mostly within 10%) because of the assumed hygroscopic nature of POM, so that freshly emitted POM and BC are wet-removed before mixing internally with soluble aerosol species. Sensitivity tests with the POM assumed to be hydrophobic and with slower aging process increase the POM and BC concentrations, especially at high latitudes (by several times). The mineral dust global burden differs by 10% and sea salt burden by 30-40% between MAM3 and MAM7 mainly due to the different size ranges for dust and sea salt modes and different standard deviations of log-normal size distribution for sea salt modes between MAM3 and MAM7. The model is able to qualitatively capture the observed geographical and temporal variations of aerosol mass and number concentrations, size distributions, and aerosol optical properties. However, there are noticeable biases, e.g., simulated sulfate and mineral dust concentrations at surface over the oceans are too low. Simulated BC concentrations are significant low in the Arctic. There is a low bias in modeled aerosol optical depth on the global scale, especially in the developing counties. There biases in aerosol simulations clearly indicate the need for improvements of aerosol processes (e.g., emission fluxes of anthropogenic aerosols and precursor gases in developing countries, boundary layer nucleation) and properties (e.g., primary aerosol emission size, POM hygroscopicity). In addition the critical role of cloud properties (e.g., liquid water content, cloud fraction) responsible for the wet scavenging of aerosol is highlighted.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2011

Testing Cloud Microphysics Parameterizations in NCAR CAM5 with ISDAC and M-PACE Observations

Arctic clouds simulated by the NCAR Community Atmospheric Model version 5 (CAM5) are evaluated with observations from the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Indirect and Semi-Direct Aerosol Campaign (ISDAC) and Mixed-Phase Arctic Cloud Experiment (M-PACE), which were conducted at its North Slope of Alaska site in April 2008 and October 2004, respectively. Model forecasts for the Arctic Spring and Fall seasons performed under the Cloud- Associated Parameterizations Testbed (CAPT) framework generally reproduce the spatial distributions of cloud fraction for single-layer boundary layer mixed-phase stratocumulus, and multilayer or deep frontal clouds. However, for low-level clouds, the model significantly underestimates the observed cloud liquid water content in both seasons and cloud fraction in the Spring season. As a result, CAM5 significantly underestimates the surface downward longwave (LW) radiative fluxes by 20-40 W m-2. The model with a new ice nucleation parameterization moderately improves the model simulations by increasing cloud liquid water content in mixed-phase clouds through the reduction of the conversion rate from cloud liquid to ice by the Wegener-Bergeron- Findeisen (WBF) process. The CAM5 single column model testing shows that change in the homogeneous freezing temperature of rain to form snow from -5 C to -40 C has a substantial impact on the modeled liquid water content through the slowing-down of liquid and rain-related processes. In contrast, collections of cloud ice by snow and cloud liquid by rain are of minor importance for single-layer boundary layer mixed-phase clouds in the Arctic.

58 GEOSCIENCES↗

DOE OSTI2011

Droplet Nucleation: Physically-Based Parameterization and Comparative Evaluation

One of the greatest sources of uncertainty in simulations of climate and climate change is the influence of aerosols on the optical properties of clouds. The root of this influence is the droplet nucleation process. Droplet nucleation involves the simultaneous condensational growth of an aerosol population until maximum supersaturation is achieved and a subset of the particles are large enough to grow spontaneously into cloud droplets. Numerical models of droplet nucleation are capable of representing much of the complexity of the process, but at a computational cost that limits their application to simulations of hours or days. Parameterizations of droplet nucleation are designed to quickly estimate the number nucleated as a function of the primary controlling parameters: the cooling rate and the size distribution of aerosol number and hygroscopicity. Here we compare and contrast the key assumptions used in developing each of the most popular parameterizations and compare the performance of each under a variety of conditions. In general we find that the more complex parameterizations perform well under a wider variety of conditions, but all parameterizations perform well under the most common conditions. We then discuss the various applications of the parameterizations to cloud-resolving, regional and global models to study aerosol effects on clouds at a wide range of spatial and temporal scales. We conclude with a summary of the outstanding challenges remaining for further development.

54 ENVIRONMENTAL SCIENCES↗