A wide range of aerosol effects on precipitation have been proposed, from the scale of individual clouds to that of the globe.This presentation, based on the findings of an expert workshop under the umbrella of the GEWEX Aerosol Precipitation initiative, reviews the evidence and scientific consensus behind these effects and the underlying set of physical mechanisms, categorised into i) radiative effects via modification of radiative fluxes and the energy balance and ii) microphysical effects via modification of cloud droplets and ice crystals.There exists broad consensus and strong theoretical evidence that, because global mean precipitation is constrained by energetics and surface evaporation, aerosol radiative effects (aerosol-radiation interactions and aerosol-cloud interactions) act as drivers of precipitation changes. Likewise, aerosol radiative effects cause well-documented shifts of large-scale precipitation patterns, such as the Inter-Tropical Convergence Zone (ITCZ). The extent to which aerosol effects on precipitation are applicable at smaller scales and driven or buffered by compensating microphysical and dynamical mechanisms and budgetary constraints is less clear. Although there exists broad consensus and strong evidence that suitable aerosol perturbations increase cloud droplet numbers, reducing the efficiency of warm rain formation across cloud regimes, the overall aerosol effect on cloud microphysics and dynamics as well as the subsequent impact on local, regional and global precipitation is less constrained.This presentation provides a review of the physical mechanisms of aerosol effects on precipitation backed up by evidence from recent cloud-resolving and global modelling simulations as well as from satellite observations.

The objectives of 7-SEAS/BASELInE (Seven SouthEast Asian Studies/Biomass-burning Aerosols & Stratocumulus Environment: Lifecycles & Interactions Experiment) campaigns in spring 2013–2015 were to synergize measurements from uniquely distributed ground-based networks (e.g., AERONET, MPLNET) and sophisticated platforms (e.g., SMARTLabs, regional contributing instruments), along with satellite observations/retrievals and regional atmospheric transport/chemical models to establish a critically needed database, and to advance our understanding of biomass-burning aerosols and trace gases in Southeast Asia (SEA). We present a satellite-surface perspective of 7-SEAS/BASELInE and highlight scientific findings concerning: (1) regional meteorology of moisture fields conducive to the production and maintenance of low-level stratiform clouds over land, (2) atmospheric composition in a biomass-burning environment, particularly tracers/markers to serve as important indicators for assessing the state and evolution of atmospheric constituents, (3) applications of remote sensing to air quality and impact on radiative energetics, examining the effect of diurnal variability of boundary-layer height on aerosol loading, (4) aerosol hygroscopicity and ground-based cloud radar measurements in aerosol-cloud processes by advanced cloud ensemble models, and (5) implications of air quality, in terms of toxicity of nanoparticles and trace gases, to human health. This volume is the third 7-SEAS special issue (after Atmospheric Research, vol. 122, 2013; and Atmospheric Environment, vol. 78, 2013) and includes 27 papers published, with emphasis on air quality and aerosol-cloud effects on the environment. BASELInE observations of stratiform clouds over SEA are unique, such clouds are embedded in a heavy aerosol-laden environment and feature characteristically greater stability over land than over ocean, with minimal radar surface clutter at a high vertical spatial resolution. To facilitate an improved understanding of regional aerosol-cloud effects, we envision that future BASELInE-like measurement/modeling needs fall into two categories: (1) efficient yet critical in-situ profiling of the boundary layer for validating remote-sensing/retrievals and for initializing regional transport/chemical and cloud ensemble models, and (2) fully utilizing the high observing frequencies of geostationary satellites for resolving the diurnal cycle of the boundary-layer height as it affects the loading of biomass-burning aerosols, air quality and radiative energetics.

The 2013 7-SEAS/BASELInE campaign over northern Southeast Asia (SEA) provided, for the first time ever, comprehensive ground-based W-band radar measurements of the low-level stratocumulus (Sc) systems that often exist during the spring over northern Vietnam in the presence of biomass-burning aerosols. Although spatially limited, ground-based remote sensing observations are generally free of the surface contamination and signal attenuation effects that often hinder space-borne measurements of these low-level cloud systems. Such observations permit detailed measurements of structures and lifecycles of these clouds as part of a broader effort to study potential impacts of these coupled aerosol-cloud systems on local and regional weather and air quality. Introductory analyses of the W-band radar data show these Sc systems generally follow a diurnal cycle, with peak occurrences during the nighttime and early morning hours, often accompanied by light precipitation. Preliminary results from idealized simulations of Sc development over land based on the observations reveal the familiar response of increased numbers and smaller sizes of cloud droplets, along with suppressed drizzle formation, as aerosol concentrations increase. Slight reductions in simulated W-band reflectivity values also are seen with increasing aerosol concentrations and result primarily from decreased droplet sizes. As precipitation can play a large role in removing aerosol from the atmosphere, and thereby improving air quality locally, quantifying feedbacks between aerosols and cloud systems over this region are essential, particularly given the negative impacts of biomass burning on human health in SEA. Such an endeavor should involve improved modeling capabilities along with comprehensive measurements of time-dependent aerosol and cloud profiles.

The Goddard Convective-Stratiform Heating (CSH) algorithm has been used to retrieve latent heating (LH) associated with clouds and cloud systems in support of the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) mission. The CSH algorithm requires the use of a cloud-resolving model (CRM) to simulate LH profiles to build look-up tables. This paper describes the current V6 CSH and its differences/similarities versus the previous V5 CSH. Long-term CRM simulations were conducted to identify the impact of the CRM resolution and convective-stratiform separation method on LH structure/profiles. The TRMM and GPM Combined radar-radiometer algorithm-derived surface rain rates and their associated precipitation properties were the input to the CSH algorithm. CSH V6-retrieved regional LH profiles in the tropics and subtropics display the classic signatures of heating in the convective region and heating over cooling in the stratiform region. Because there is no direct measurement of LH structure, the performance of the CSH V6 algorithm is examined by comparing its vertically integrated heating (or equivalent surface rain rate) against the surface rain rate derived from the TRMM/GPM Combined algorithm. The CSH three-month and zonal mean equivalent surface rain rates are in good agreement with the Combined rain rates over the Inter Tropical Convergence Zone region; the agreement is best over the ocean. CSH three-month and zonal mean equivalent surface rain rates are larger than the Combined rain rates over land in both the tropics and subtropics. CSH three-month mean equivalent surface rain rates also have local differences with the Combined rain rates that can be smoothed by area averaging to larger horizontal resolutions (from the CSH standard grid of 0.25° × 0.25° to 0.5° × 0.5° or 1.0° × 1.0°). CSH equivalent surface rain rates have more light rain rates but less larger rates compared to the GPM Combined surface rain rates.

Fundamental relationships exist between cloud and precipitation development and their dynamic processes. Latent heat released by cloud/precipitation formation affects cloud vertical motions, which in turn affect convective cloud development. Here, a cloud‐resolving model is used to relate cloud properties and latent heating with cloud drafts using 15‐day simulations for an oceanic and continental environment. The results show condensation, deposition, and freezing occur mainly in moderate (3–5 m s −1 ) to strong (>10 m s −1 ) updrafts, evaporation and sublimation mainly in weak (1–2 m s −1 ) to moderate downdrafts, and melting in moderate updrafts and downdrafts. Active updrafts cover only a small percentage of the model domain but contribute significantly to the latent heat release and are associated with large proportions of the hydrometeors. Active updrafts with vertical velocities exceeding 1 and 2 m s −1 account for more than 75% and 50%, respectively, of the condensation, deposition, and freezing in both the oceanic and continental cases. However, active downdrafts with vertical velocity magnitudes exceeding |1 m s −1 | account for less than 40% and 25%, respectively, of the evaporation and sublimation. More evaporation and sublimation than condensation and deposition occur in the inactive cloud regions. Sensitivity tests are also conducted to assess the impact of model grid spacing (1,000 m vs. 250 m) and microphysical schemes (3 ice classes vs. 4 ice classes) on latent heat release and hydrometeor amount. The results show that model resolution had more impact than the microphysics on the simulated cloud properties in both cases.

This study contrasts midlatitude continental and tropical maritime deep convective cores using polarimetric radar observables and retrievals from selected deep convection episodes during the MC3E and TWPICE field campaigns. The continental convective cores produce stronger radar reflectivities throughout the profiles, while maritime convective cores produce more positive differential reflectivity Z dr and larger specific differential phase K dp above the melting level. Hydrometeor identification retrievals revealed the presence of large fractions of rimed ice particles (snow aggregates) in the continental (maritime) convective cores, consistent with the Z dr and K dp observations. The regional cloud-resolving model simulations with bulk and size-resolved bin microphysics are conducted for the selected cases, and the simulation outputs are converted into polarimetric radar signals and retrievals identical to the observational composites. Both the bulk and the bin microphysics reproduce realistic land and ocean (L-O) contrasts in reflectivity, polarimetric variables of rain drops, and hydrometeor profiles, but there are still large uncertainties in describing Z dr and K dp of ice crystals associated with the ice particle shapes/orientation assumptions. Sensitivity experiments are conducted by swapping background aerosols between the continental and maritime environments, revealing that background aerosols play a role in shaping the distinct L-O contrasts in radar reflectivity associated with raindrop sizes, in addition to the dominant role of background thermodynamics.

In this study, the impact of different configurations of the Goddard radiation scheme on convection-permitting simulations (CPSs) of the West African monsoon (WAM) is investigated using the NASA-Unified WRF (NU-WRF). These CPSs had 3 km grid spacing to explicitly simulate the evolution of mesoscale convective systems (MCSs) and their interaction with radiative processes across the WAM domain and were able to reproduce realistic precipitation and energy budget fields when compared with satellite data, although low clouds were overestimated. Sensitivity experiments reveal that (1) lowering the radiation update frequency (i.e., longer radiation update time) increases precipitation and cloudiness over the WAM region by enhancing the monsoon circulation, (2) deactivation of precipitation radiative forcing suppresses cloudiness over the WAM region, and (3) aggregating radiation columns reduces low clouds over ocean and tropical West Africa. The changes in radiation configuration immediately modulate the radiative heating and low clouds over ocean. On the 2nd day of the simulations, patterns of latitudinal air temperature profiles were already similar to the patterns of monthly composites for all radiation sensitivity experiments. Low cloud maintenance within the WAM system is tightly connected with radiation processes; thus, proper coupling between microphysics and radiation processes must be established for each modeling framework.

Cloud and precipitation systems are simulated with a multi-scale modeling framework (MMF) and compared over the Tropics and Subtropics against the Tropical Rainfall Measuring Mission (TRMM) Radar-defined Precipitation Features (RPFs) product. A methodology, in close analogy to the TRMM RPFs, is developed to produce simulated precipitation features (PFs) from the output of the embedded two-dimensional (2D) cloud-resolving models (CRMs) within an MMF. Despite the limitations of 2D CRMs, the simulated population distribution, horizontal and vertical structure of PFs, and the geographical location and local rainfall contribution of mesoscale convective systems (MCSs) are in good agreement with the TRMM observations. However, some model discrepancies are found and can be identified and quantified within the PF distributions. Using model biases in relative population and rainfall contributions, PFs can be characterized into four size categories: small, medium to large, very large, and extremely large. Four different major mechanisms might account for the model biases in each different category: (1) the two-dimensionality of the CRMs, (2) a positive convection-wind-evaporation feedback loop, (3) an artificial dynamic constraint in a bounded CRM domain with cyclic boundaries, and (4) the limited CRM domain size. The second and fourth mechanisms tend to contribute to the excessive tropical precipitation biases commonly found in most MMFs, whereas the other mechanisms reduce rainfall contributions from small and very large PFs. MMF sensitivity experiments with various CRM domain sizes and grid spacings showed that larger domains (higher resolutions) tend to shift PF populations toward larger (smaller) sizes.

Latent heat release itself is a consequence of phase changes between the vapor, liquid, and frozen states of water. This paper will highlight the retrieval of latent heat release from satellite measurements generated by the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Mission (GPM) satellite observatory. Both TRMM and GPM are providing four-dimensional account of rainfall and its associated precipitation properties over the global Tropics and mid-latitudes: information that can be used to estimate the space-time structure of latent heating. Goddard Convective-Stratiform or CSH retrieved LH is one of two standard LH products (the other one is from Japan Spectral Latent Heating or SLH). This paper will present (1) the new improvements of the CSH LH algorithm by better CRM simulated LH for its look-up table, and (2) the performance of CSH retrieved LH by comparison with surface rainfall rate. In addition, the similarities and differences of CSH retrieved LH obtained from TRMM and GPM measurements, respectively, will be presented. At the end of presentation, the further research on latent heating retrieval from satellites will be discussed.

With the launch of the Tropical Rainfall Measuring Mission satellite, equipped with high-functioning sensors to measure precipitation, four-dimensional latent heating retrievals from satellite observation have become available. With the launch of the Global Precipitation Measurement mission Core Observatory, the latent heating retrievals are now available at higher latitudes. Algorithms of two standard NASA/JAXA products, the Spectral Latent Heating, and the Convective Stratiform Heating and some results are introduced.

This study focuses on ocean surface fluxes, mainly the latent heat flux, and their impact on MJO propagation and associated precipitation structures over the Indian Ocean and Maritime Continent. The Coupled-Ocean-Atmosphere-Wave-Sediment Transport (COAWST) model is used to simulate two MJO events during the 2017-2018 season: the December 10 - January 20, 2017 case, which maintained its strong precipitation signal over the Maritime Continent, and the March 1 - 20, 2018 case, which was weaker and did not propagate through the Maritime Continent. Both simulated MJO events show positive biases in surface rainfall compared with GPM IMERG data. During the MJO suppressed phase, the simulations rain more often than the observations. During the active phase, the westward propagating precipitation structures are more organized and much stronger compared with the observations, sometimes forming westward propagating cyclones that weakened the eastward precipitation signals. Two aspects of the surface flux interactions are investigated: the impact of the domain mean surface fluxes, and the impact of storm scale circulations and their interactions with local surface fluxes. Both aspects affect water vapor budget, atmosphere instability and mean flow, through which convection initiation, organization, and propagation are influenced. Model sensitivity tests with different radiation, microphysics, PBL schemes and nudging schemes indicate that in the control simulations, higher SST and surface fluxes, especially during the suppressed period, are the main reason of rainfall overestimation compared with IMERG data. The strong westward propagating signals are caused by both increased atmosphere instability and reduced mean wind shear. Unfortunately, the small differences in mean SST and surface fluxes between different model sensitivity tests are all within the satellite observation error margin, and cannot be directly corroborated by observations. One of the advantages of CYGNSS satellites is that they observe ocean surface wind and heat fluxes underneath strong rainfall events such as the convective systems associated with MJO active phases. Currently we are comparing CYGNSS level 2 surface fluxes retrievals and the model simulations in order to better understand the second aspect of the MJO and surface fluxes interactions, and how this affects MJO strengths and propagations. The interactive atmosphere-ocean-wave model also provides cases that directly comparing satellite observables (the bistatic radar cross section) and the model simulations (through CYGNSS satellite simulator). These discrepancies are more prominent in coupled ocean simulations, mainly due to higher SST and enhanced surface fluxes.

In this paper, we have investigated the impact of radiation–cloud–convection–circulation interaction (RC3I) on structural changes and variability of the Inter-tropical Convergence Zone (ITCZ) using the Goddard Multi-scale Modeling Framework, where cloud processes are super-parameterized, i.e., explicitly resolved with 2-D cloud resolving models embedded in each coarse grid of the host Goddard Earth Observing System-Version 5 global climate model. Experiments have been conducted under prescribed sea surface temperature conditions for 10 years (2007–2016), with and without cloud radiation feedback in the atmosphere, respectively. Diagnostic analyses separately for January and July show that RC3I leads to an enhanced and expanded Hadley Circulation characterized by (1) a quasi-uniform warming and moistening of the tropical atmosphere and a sharpening of the ITCZ with enhanced deep convection, more intense precipitation and higher clouds, (2) extended drying of the tropical marginal convective zones, and extratropical mid- to lower troposphere, and (3) a cooling of the polar regions, with increased baroclinicity and midlatitude storm track activities. Computations based on the zonal mean thermodynamic energy balance equation show that the radiative warming and cooling are strongly balanced by local adiabatic processes associated with changes in large-scale vertical motions, as well as horizontal atmospheric heat transport. In the tropics, enhanced short-wave absorption and longwave water vapor greenhouse effects by high clouds play key roles in providing strong positive feedback to the tropospheric warming. In the extratropics, increased atmospheric heat transport associated with changes in the Hadley circulation is balanced by strong longwave cooling above, and warming below due to increased high clouds. We also find a strong positive correlation between daily and pentad heavy rain in the ITCZ core, and expansion of the drier zones coupled to a contraction of the highly convective zones in the ITCZ, indicating a strong tendency RC3I-induced convective aggregation in tropical clouds i.e., wet-regions-get-wetter and contracted, and dry-areas-get-drier and expanded.

25 A series of model simulations were conducted to investigate the effects of cloud 26 condensation nuclei (CCN) loading and convective available potential energy (CAPE) on 27 tropical maritime and mid-latitude continental deep convection. Dynamical downscaling from 28 global aerosol reanalysis was used to represent aerosol fields for the two convective regimes. 29 We describe a control run and multiple sensitivity experiments using a limited-area 30 model, employing spectral-bin cloud microphysics. The CCN loading is perturbed between the 31 target maritime and continental conditions, roughly 40 – 2000 cm at 850 hPa and 1% 32 supersaturation. Surface precipitation rates monotonically increase with increasing CCN loading 33 for both the maritime and continental situations, while these monotonic increases are disrupted in 34 the simulations with reduced CAPE. The increase in precipitation is in the form of convective 35 precipitation, at the expense of stratiform precipitation. 36 CCN increases promote increases in supercooled cloud water, in agreement with previous 37 modeling studies. However, in the simulations investigated herein, the changes in supercooled 38 water have different impacts on the cloud microphysics in the maritime and continental 39 simulations. Increased supercooled water contents lead to more hail and less graupel in the 40 continental simulation. For the maritime simulation, enhanced supercooled cloud water contents 41 promote an increase in graupel since little or no hail is produced. This distinction is due to the 42 difference in relative magnitudes and peak altitudes of supercooled water and snow amounts, 43 which is further attributable to moisture and dynamical differences in the two cases. 44 45 46 Confidential manuscript submitted to Journal of Geophysical Research: Atmospheres

A series of model simulations were conducted to investigate the effects of cloud condensation nuclei (CCN) loading and convective available potential energy (CAPE) on tropical maritime and midlatitude continental deep convection. Dynamical downscaling from global aerosol reanalysis was used to represent aerosol fields for the two convective regimes. We describe a control run and multiple sensitivity experiments using a limited‐area model, employing spectral‐bin cloud microphysics. The CCN loading is perturbed between the target maritime and continental conditions, roughly 40–2,000 cm−3 at 850 hPa and 1% supersaturation. Surface precipitation rates monotonically increase with increasing CCN loading for both the maritime and continental situations, while these monotonic increases are disrupted in the simulations with reduced CAPE. The increase in precipitation is in the form of convective precipitation, at the expense of stratiform precipitation. CCN increases promote increases in supercooled cloud water, in agreement with previous modeling studies. However, in the simulations investigated herein, the changes in supercooled water have different impacts on the cloud microphysics in the maritime and continental simulations. Increased supercooled water contents lead to more hail and less graupel in the continental simulation. For the maritime simulation, enhanced supercooled cloud water contents promote an increase in graupel since little or no hail is produced. This distinction is due to the difference in relative magnitudes and peak altitudes of supercooled water and snow amounts, which is further attributable to moisture and dynamical differences in the two cases.

Rainfall production is a fundamental process within the Earth's hydrological cycle because it represents both a principal forcing term in surface water budgets, and its energetics corollary, latent heating, is the principal source of atmospheric diabatic heating. Latent heat release itself is a consequence of phase changes between the vapor, liquid, and frozen states of water. The properties of the vertical distribution of latent heat release modulate large-scale meridional and zonal circulations within the Tropics - as well as modify the energetic efficiencies of mid-latitude weather systems. This paper highlights the retrieval of latent heat release from satellite measurements generated by the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) Mission satellite observatory, which were launched in November 1997 and February 2014, respectively. Both TRMM and GPM measurements have been providing an accurate four-dimensional account of rainfall over the global Tropics and mid-latitudes: information that can be used to estimate the space-time structure of latent heating.Two sets of latent heating (LH) retrieval algorithm methodologies (Goddard Convective-Stratiform Heating or CSH, and Japan's Spectral Latent Heating or SLH) have been developed to estimate latent heating based on rain rate profile retrievals obtained from TRMM and GPM measurements. The differences and similarities between the CSH and SLH algorithm designs will be presented. The presentation then will provide the results of LH structures derived from the TRMM and GPM. At the end of presentation, the further research on latent heating retrieval from satellites will be discussed.

The Goddard Convective-Stratiform Heating (CSH) algorithm has been used to retrieve latent heating (LH) associated with clouds and cloud systems in support of the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) mission. The CSH algorithm required the use of a cloud-resolving model (CRM) to simulate LH profiles to build look-up tables (LUTs). However, the current LUTs in the CSH algorithm are not suitable for retrieving LH profiles at high latitudes or winter conditions that are needed for GPM. The NASA Unified-Weather Research and Forecasting (NU-WRF) model is used to simulate three eastern continental US (CONUS) synoptic winter and three western coastal/offshore events. The relationship between LH structures (or profiles) and other precipitation properties (radar reflectivity, freezing level height, echo-top height, maximum radar reflectivity height and surface precipitation rate) is examined, and a new classification system is adopted with varying ranges for each of these precipitation properties to create LUTs representing high latitude/winter conditions. The performance of the new LUTs is examined using a self-consistency check for one CONUS and one West Coast offshore event by comparing LH profiles retrieved from the LUTs using model-simulated precipitation properties with those originally simulated by the model. The results of the self-consistency check validate the new classification and LUTs. High latitude retrievals from the new LUTs are merged with those from the CSH algorithm to retrieve LH profiles over the GPM domain using precipitation properties retrieved from the GPM combined algorithm.

This paper introduces a synthetic polarimetric radar simulator and retrieval package, POLArimetric Radar Retrieval and Instrument Simulator (POLARRIS), for evaluating cloud‐resolving models (CRMs). POLARRIS is composed of forward (POLARRIS‐f) and inverse (retrieval and diagnostic) components (iPOLARRIS) to generate not only polarimetric radar observables (Zh, Zdr, Kdp, ρhv) but also radar‐consistent geophysical parameters such as hydrometeor identification, vertical velocity, and rainfall rates retrieved from CRM data. To demonstrate its application and uncertainties, POLARRIS is applied to simulations of a mesoscale convective system over the Southern Great Plains on 23 May 2011, using the Weather Research and Forecasting model with both spectral bin microphysics (SBM) and the Goddard single‐moment bulk 4ICE microphysics. Statistical composites reveal a significant dependence of simulated polarimetric observables (Zdr, Kdp) on the assumptions of the particle axis ratio (oblateness) and orientation angle distributions. The simulated polarimetric variables differ considerably between the SBM and 4ICE microphysics in part due to the differences in their ice particle size distributions as revealed by comparisons with aircraft measurements. Regardless of these uncertainties, simulated hydrometeor identification distributions overestimate graupel and hail fractions, especially from the simulation with SBM. To minimize uncertainties in forward model, the particle shape and orientation angle distributions of frozen particles should be predicted in a microphysics scheme in addition to the size distributions and particle densities.

Advances in computing power allow atmospheric prediction and general circulation models to be run at progressively finer scales of resolution, using increasingly more sophisticated physical parameterizations. The representation of cloud microphysical processes is one of key components of these models. In addition, over the past decade both research and operational numerical weather prediction models have started using more complex microphysical schemes that were originally developed for high-resolution cloud-resolving models (CRMs). In the paper, we described different microphysics schemes that are used in Goddard Multi-scale Modeling System. There are three major models, Goddard Cumulus Ensemble (GCE), NASA Unified Weather Research Forecast (NU-WRF) and Multi-scale Modeling Framework (MMF) model, in this modeling system. The microphysics schemes are Goddard three class ice (3ICE) and four class (4ICE) scheme, Morrison two moments (2M) 3ICE, Colorado State University Regional Atmospheric Modeling System (RAMS) 2M five class ice (5ICE) and spectral bin microphysics schemes. The performance of these schemes are examined and compared with radar and satellite observation. In addition, the inter-comparison with different microphysics schemes are conducted. Current and future observations needed for microphysics schemes evaluation as well as major characteristics of current microphysics are discussed.