Cooper, Vincent T., Kyle Armour, Gregory J Hakim, Jessica E Tierney, Matthew B Osman, Cristian Proistosescu, Yue Dong, Natalie J Burls, Timothy Andrews, Daniel E Amrhein, Jiang Zhu, Wenhao Dong, Yi Ming, and Philip Chmielowiec, April 2024: Last Glacial Maximum pattern effects reduce climate sensitivity estimates. Science Advances, 10(16), DOI:10.1126/sciadv.adk9461. Abstract
Here, we show that the Last Glacial Maximum (LGM) provides a stronger constraint on equilibrium climate sensitivity (ECS), the global warming from increasing greenhouse gases, after accounting for temperature patterns. Feedbacks governing ECS depend on spatial patterns of surface temperature (“pattern effects”); hence, using the LGM to constrain future warming requires quantifying how temperature patterns produce different feedbacks during LGM cooling versus modern-day warming. Combining data assimilation reconstructions with atmospheric models, we show that the climate is more sensitive to LGM forcing because ice sheets amplify extratropical cooling where feedbacks are destabilizing. Accounting for LGM pattern effects yields a median modern-day ECS of 2.4°C, 66% range 1.7° to 3.5°C (1.4° to 5.0°C, 5 to 95%), from LGM evidence alone. Combining the LGM with other lines of evidence, the best estimate becomes 2.9°C, 66% range 2.4° to 3.5°C (2.1° to 4.1°C, 5 to 95%), substantially narrowing uncertainty compared to recent assessments.
The Taklamakan and Gobi Desert (TGD) region has experienced a pronounced increase in summer precipitation, including high-impact extreme events, over recent decades. Despite identifying large-scale circulation changes as a key driver of the wetting trend, understanding the relative contributions of internal variability and external forcings remains limited. Here, we approach this problem by using a hierarchy of numerical simulations, complemented by diverse statistical analysis tools. Our results offer strong evidence that the atmospheric internal variations primarily drive this observed trend. Specifically, recent changes in the North Atlantic Oscillation have redirected the storm track, leading to increased extratropical storms entering TGD and subsequently more precipitation. A clustering analysis further demonstrates that these linkages predominantly operate at the synoptic scale, with larger contributions from large precipitation events. Our analysis highlights the crucial role of internal variability, in addition to anthropogenic forcing, when seeking a comprehensive understanding of future precipitation trends in TGD.
Coupled ocean and prescribed sea surface temperature (SST) experiments are performed to investigate the drivers of Northern Hemisphere (NH) midlatitude winter circulation and blocking changes in warmer climates. In coupled experiments, a historical simulation is compared to a simulation following an end of the twenty-first-century shared socioeconomic pathway (SSP5-8.5) emission scenario. The SSP5-8.5 simulation yields poleward-shifted jets and an enhanced stationary wave pattern compared to the historical simulation. In terms of blocking, a reduction is found across North America and over the Pacific Ocean with the suggestion of more blocking over parts of Eurasia. Separately, prescribed SST experiments are performed decomposing the SSP5-8.5 SST response into a uniform warming component plus a spatially dependent change in SST pattern. SSP5-8.5 changes in circulation are primarily driven by a uniform warming of SST. Uniform warming is also found to account for most of the SSP5-8.5 blocking reduction over North America and the Pacific Ocean, but not over Eurasia. El Niño–like changes to the SST pattern also yield less blocking over the Pacific and North America. However, adding the responses of uniform and pattern experiments yields a nonlinear overreduction of blocking compared to the SSP5-8.5 experiment. Regional analyses of block energetics suggest that much of the reductions in blocking in warming simulations are driven by decreased baroclinic conversion in some regions and enhanced dissipation from diabatic sources in others.
Nikumbh, Akshaya C., Pu Lin, David J Paynter, and Yi Ming, June 2024: Does increasing horizontal resolution improve the simulation of intense tropical rainfall in GFDL's AM4 model?Geophysical Research Letters, 51(12), DOI:10.1029/2023GL106708. Abstract
We examine tropical rainfall from the Geophysical Fluid Dynamics Laboratory's Atmosphere Model version 4 (GFDL AM4) at three horizontal resolutions of 100 km, 50 km, and 25 km. The model produces more intense rainfall at finer resolutions, but a large discrepancy still exists between the simulated and the observed frequency distribution. We use a theoretical precipitation scaling diagnostic to examine the frequency distribution of the simulated rainfall. The scaling accurately produces the frequency distribution at moderate-to-high intensity (≥10 mm day−1). Intense tropical rainfall at finer resolutions is produced primarily from the increased contribution of resolved precipitation and enhanced updrafts. The model becomes more sensitive to the grid-scale updrafts than local thermodynamics at high rain rates as the contribution from the resolved precipitation increases.
Stier, Philip, Susan C van den Heever, Matthew W Christensen, Edward Gryspeerdt, Guy Dagan, Stephen M Saleeby, Massimo Bollasina, Leo J Donner, Kerry A Emanuel, Annica M L Ekman, Graham Feingold, Paul Field, Piers M Forster, Jim M Haywood, Ralph A Kahn, Ilan Koren, Christian Kummerow, Tristan L'Ecuyer, Ülrike Lohmann, Yi Ming, Gunnar Myhre, Johannes Quaas, Daniel Rosenfeld, Bjørn H Samset, Axel Seifert, Graeme L Stephens, and Wei-Kuo Tao, August 2024: Multifaceted aerosol effects on precipitation. Nature Geoscience, 17, DOI:10.1038/s41561-024-01482-6. Abstract
Aerosols have been proposed to influence precipitation rates and spatial patterns from scales of individual clouds to the globe. However, large uncertainty remains regarding the underlying mechanisms and importance of multiple effects across spatial and temporal scales. Here we review the evidence and scientific consensus behind these effects, categorized into radiative effects via modification of radiative fluxes and the energy balance, and microphysical effects via modification of cloud droplets and ice crystals. Broad consensus and strong theoretical evidence exist that aerosol radiative effects (aerosol–radiation interactions and aerosol–cloud interactions) act as drivers of precipitation changes because global mean precipitation is constrained by energetics and surface evaporation. Likewise, aerosol radiative effects cause well-documented shifts of large-scale precipitation patterns, such as the intertropical convergence zone. The extent of aerosol effects on precipitation at smaller scales is less clear. Although there is broad consensus and strong evidence that aerosol perturbations microphysically increase cloud droplet numbers and decrease droplet sizes, thereby slowing precipitation droplet formation, the overall aerosol effect on precipitation across scales remains highly uncertain. Global cloud-resolving models provide opportunities to investigate mechanisms that are currently not well represented in global climate models and to robustly connect local effects with larger scales. This will increase our confidence in predicted impacts of climate change.
You, Zhenyu, Yi Deng, Yi Ming, and Wenhao Dong, February 2024: A multiscale assessment of the springtime U.S. mesoscale convective systems in the NOAA GFDL AM4. Climate Dynamics, DOI:10.1007/s00382-024-07114-4. Abstract
This study presents a multiscale assessment of the springtime U.S. Mesoscale Convective Systems (MCSs) in the NOAA Geophysical Fluid Dynamics Laboratory (GFDL)’s Atmosphere Model version 4 (AM4). In AM4, MCSs exhibit lower intensity but longer duration, producing more precipitation compared to observation. The overall MCS activity demonstrates a “location bias” with its peak shifting from the Southern Great Plains to the Midwest in AM4, causing an eastward shift in associated precipitation. However, the dry bias of MCS precipitation over the Great Plains due to this shift is compensated by additional precipitation from amplified extratropical cyclone activities. Further analysis reveals that AM4 effectively reproduces the spatiotemporal distribution and relative frequency contribution of large-scale forcing patterns driving MCS genesis. The MCS location bias emerges under all forms of large-scale forcing patterns and is further attributed to local dynamic and thermodynamic factors including weaker surface lows, eastward-shifted fronts, and suppressed low-level jets (LLJs). Here we argue that the MCS location bias results from AM4 biases in both synoptic-mesoscale anomalies (i.e., fronts and LLJs) and seasonal mean circulations. The lack of two-way air-sea interaction in AM4 creates a hemispheric-scale sea level pressure bias, which is ultimately responsible for a seasonal mean northerly bias in lower-tropospheric winds and the subsequent weakening of LLJs. The existence of such biases in prescribed sea surface temperature (SST) experiments implies the need for extra caution when utilizing extended-range forecasts for MCSs over the continental U.S.
Accurate representation of mesoscale scale convective systems (MCSs) in climate models is of vital importance to understanding global energy, water cycles, and extreme weather. In this study, we evaluate the simulated MCS features over the United States from the newly developed GFDL global high-resolution (∼50 km) AM4 model by comparing them with the observations during spring to early summer (April–June) and late summer (July–August). The results show that the spatial distribution and seasonality of occurrence and genesis frequency of MCSs are reasonably simulated over the central United States in both seasons. The model reliably reproduces the observed features of MCS duration, translation speed, and size over the central United States, as well as the favorable large-scale circulation pattern associated with MCS development over the central United States during spring and early summer. However, the model misrepresents the amplitude and the phase of the diurnal cycle of MCSs during both seasons. In addition, the spatial distribution of occurrence and genesis frequency of MCSs over the eastern United States is substantially overestimated, with larger biases in early spring and summer. Furthermore, while large-scale circulation patterns are reasonably simulated in spring and early summer, they are misrepresented in the model during summer. Finally, we examine MCS-related precipitation, finding that the model overestimates MCS-related precipitation during spring and early summer, but this bias is insufficient to explain the significant dry bias observed in total precipitation over the central United States. Nonetheless, the dry biases in MCS-associated precipitation during late summer likely contribute to the overall precipitation deficit in the model.
We performed a series of aquaplanet simulations at the horizontal resolution from 50 to 6 km with identical parameterization settings using the Geophysical Fluid Dynamics Laboratory's Atmosphere Model version 4 implemented with the two-moment Morrison-Gettelman cloud microphysics with prognostic precipitation (GFDL AM4-MG2). At the finer resolution, the global mean resolved-scale precipitation increases and that from cumulus parameterization decreases. The model also simulates less/thinner clouds over the low latitudes and more/thicker clouds over the high latitudes as the resolution increases. The precipitation over the deep tropics is investigated in detail. We find little resolution sensitivity in the daily mean precipitation extremes. Changes of the equatorial resolved precipitation with resolution cannot be fully explained by the resolution dependence in the vertical velocity amplitude. We report a robust sensitivity in the convective organization over the deep tropics to the model resolution. In simulations of finer resolution, the localized convection is suppressed, and the organized convective system associated with large-scale circulations becomes more prominent.
Marine cloud brightening is a proposal to counteract global warming by increasing sea salt aerosol emissions. In theory, this increases the cloud droplet number concentration of subtropical marine stratocumulus decks, increasing cloud brightness and longevity. However, this theoretical progression remains uncertain in coupled climate models, especially the response of liquid water path and cloud fraction to aerosol seeding. We use the GFDL CM4 climate model to simulate marine cloud brightening following the published G4sea-salt protocol, in which sea salt aerosol emissions are uniformly increased over 30 S–30 N in addition to standard forcings from a SSP2-4.5 future warming scenario. The perturbed radiative and cloud responses are temporally stable though spatially heterogeneous, and direct scattering by the added sea salt predominates over changes to cloud reflectance. In fact, feedbacks in the coupled simulation lead to a net warming, rather than cooling, response by clouds.
Narinesingh, Veeshan, James F Booth, and Yi Ming, January 2023: Northern hemisphere heat extremes in a warmer climate: More probable but less colocated with blocking. Geophysical Research Letters, 50(2), DOI:10.1029/2022GL101211. Abstract
This work uses reanalysis and NOAA Geophysical Fluid Dynamics Laboratory's Coupled Model Intercomparison Project 6 model, CM4, to investigate the colocation of heat extremes and atmospheric blocking in the current climate and an end of 21st century, extreme-emissions projection. In the present day, the colocation of heat events and blocking is greatest for the strongest heat events. Block-heat extreme colocation is found to be less prevalent over ocean than land, exhibiting regional variation throughout the Northern Hemisphere. Over North America, colocation is greatest near the northwestern and northeastern coasts, minimizing near the center; over Eurasia, colocation is most prevalent in northern regions. In an RCP 8.5 projection, the historical 90th percentile temperature decreases to 0–70th percentile, depending on the region. This is primarily driven by mean state warming. Blocking is found to decrease along with the colocation of blocking and heat extremes, suggesting that in some regions, the mechanisms driving heat extremes will change in future climates.
Neelin, J D., John P Krasting, Aparna Radhakrishnan, Jessica Liptak, Thomas Jackson, Yi Ming, Wenhao Dong, Andrew Gettelman, Danielle R Coleman, Eric Maloney, Allison A Wing, and Yi-Hung Kuo, et al., August 2023: Process-oriented diagnostics: Principles, practice, community development and common standards. Bulletin of the American Meteorological Society, 104(8), DOI:10.1175/BAMS-D-21-0268.1E1452-E1468. Abstract
Process-oriented diagnostics (PODs) aim to provide feedback for model developers through model analysis based on physical hypotheses. However, the step from a diagnostic based on relationships among variables, even when hypothesis driven, to specific guidance for revising model formulation or parameterizations can be substantial. The POD may provide more information than a purely performance-based metric, but a gap between POD principles and providing actionable information for specific model revisions can remain. Furthermore, in coordinating diagnostics development, there is a trade-off between freedom for the developer, aiming to capture innovation, and near-term utility to the modeling center. Best practices that allow for the former, while conforming to specifications that aid the latter, are important for community diagnostics development that leads to tangible model improvements. Promising directions to close the gap between principles and practice include the interaction of PODs with perturbed physics experiments and with more quantitative process models as well as the inclusion of personnel from modeling centers in diagnostics development groups for immediate feedback during climate model revisions. Examples are provided, along with best-practice recommendations, based on practical experience from the NOAA Model Diagnostics Task Force (MDTF). Common standards for metrics and diagnostics that have arisen from a collaboration between the MDTF and the Department of Energy’s Coordinated Model Evaluation Capability are advocated as a means of uniting community diagnostics efforts.
Interpreting behaviors of low-level clouds (LLCs) in a climate model is often not straightforward. This is particularly so over polar oceans where frozen and unfrozen surfaces coexist, with horizontal winds streaming across them, shaping LLCs. To add clarity to this interpretation issue, we conduct budget analyses of LLCs using a global atmosphere model with a fully prognostic cloud scheme. After substantiating the model’s skill in reproducing observed LLCs, we use the modeled budgets of cloud fraction and water content to elucidate physics governing changes of LLCs across sea ice edges. Contrasting LLC regimes between open water and sea ice are found. LLCs over sea ice are primarily maintained by large-scale condensation: intermittent intrusions of maritime humid air and surface radiative cooling jointly sustain high relative humidity near the surface, forming extensive but tenuous stratus. This contrasts with the LLCs over open water where the convection and boundary layer condensation sustain the LLCs on top of deeper boundary layers. Such contrasting LLC regimes are influenced by the direction of horizontal advection. During on-ice flow, large-scale condensation dominates the regions, both open water and sea ice regions, forming clouds throughout the lowest several kilometers of the troposphere. During off-ice flow, as cold air masses travel over the open water, the cloud layer lifts and becomes denser, driven by increased surface fluxes that generate LLCs through boundary layer condensation and convective detrainment. These results hold in all seasons except summer when the atmosphere–surface decoupling substantially reduces the footprints of surface type changes.
Dong, Wenhao, and Yi Ming, September 2022: Seasonality and variability of snowfall to total precipitation ratio over high mountain Asia simulated by the GFDL high-resolution AM4. Journal of Climate, 35(17), DOI:10.1175/JCLI-D-22-0026.15573-5589. Abstract
The ratio of snowfall to total precipitation (S/P ratio) is an important metric that is widely used to detect and monitor hydrologic responses to climate change over mountainous areas. Changes in the S/P ratio over time have proved to be reliable indicators of climatic warming. In this study, the seasonality and interannual variability of monthly S/P ratios over High Mountain Asia (HMA) have been examined during the period 1950–2014 based on a three-member ensemble of simulations using the latest GFDL AM4 model. The results show a significant decreasing trend in S/P ratios during the analysis period, which has mainly resulted from reductions in snowfall, with increases in total precipitation playing a secondary role. Significant regime shifts in S/P ratios are detected around the mid-1990s, with rainfall becoming the dominant form of precipitation over HMA after the changepoints. Attribution analysis demonstrates that increases in rainfall during recent decades were primarily caused by a transformation of snowfall to rainfall as temperature warmed. A logistic equation is used to explore the relationship between the S/P ratio and surface temperature, allowing calculation of a threshold temperature at which the S/P ratio equals 50% (i.e., precipitation is equally likely to take the form of rainfall or snowfall). These temperature thresholds are higher over higher elevations. This study provides an extensive evaluation of simulated S/P ratios over the HMA that helps clarify the seasonality and interannual variability of this metric over the past several decades. The results have important socioeconomic and environmental implications, particularly with respect to water management in Asia under climate change.
An event-based assessment of the sea surface temperature (SST) threshold at the genesis of tropical mesoscale convective systems (MCSs) is performed in this study. We show that this threshold (SSTG) has undergone a significant warming trend at a rate of ∼0.2°C per decade. The SSTG shows a remarkable correspondence with the tropical mean SST and upper-tropospheric temperature on interannual and longer timescales. Using a high-resolution global climate model that permits realistic simulations of tropical MCSs, we find that the observed features of SSTG are well simulated. Both observation and model simulations demonstrate that the upward tendency in SSTG primarily results from the environmental SST warming over MCS genesis regions rather than the changes in MCS genesis location. A continuous increase in SSTG is projected in a warming simulation, but the relationship between SSTG and upper-tropospheric temperature remains unchanged, suggesting that the tropical tropospheric temperature generally follows a moist-adiabatic adjustment.
Frazer, Michelle E., and Yi Ming, April 2022: Understanding the extratropical liquid water path feedback in mixed-phase clouds with an idealized global climate model. Journal of Climate, 35(8), DOI:10.1175/JCLI-D-21-0334.12391-2406. Abstract
A negative shortwave cloud feedback associated with higher extratropical liquid water content in mixed-phase clouds is a common feature of global warming simulations, and multiple mechanisms have been hypothesized. A set of process-level experiments performed with an idealized global climate model (a dynamical core with passive water and cloud tracers and full Rotstayn–Klein single-moment microphysics) show that the common picture of the liquid water path (LWP) feedback in mixed-phase clouds being controlled by the amount of ice susceptible to phase change is not robust. Dynamic condensate processes—rather than static phase partitioning—directly change with warming, with varied impacts on liquid and ice amounts. Here, three principal mechanisms are responsible for the LWP response, namely higher adiabatic cloud water content, weaker liquid-to-ice conversion through the Bergeron–Findeisen process, and faster melting of ice and snow to rain. Only melting is accompanied by a substantial loss of ice, while the adiabatic cloud water content increase gives rise to a net increase in ice water path (IWP) such that total cloud water also increases without an accompanying decrease in precipitation efficiency. Perturbed parameter experiments with a wide range of climatological LWP and IWP demonstrate a strong dependence of the LWP feedback on the climatological LWP and independence from the climatological IWP and supercooled liquid fraction. This idealized setup allows for a clean isolation of mechanisms and paints a more nuanced picture of the extratropical mixed-phase cloud water feedback than simple phase change.
Frazer, Michelle E., and Yi Ming, August 2022: Understanding controlling factors of extratropical humidity and clouds with an idealized general circulation model. Journal of Climate, 35(16), DOI:10.1175/JCLI-D-21-0580.15321-5337. Abstract
This paper examines the physical controls of extratropical humidity and clouds by isolating the effects of cloud physics factors in an idealized model. The Held–Suarez dynamical core is used with the addition of passive water vapor and cloud tracers, allowing cloud processes to be explored cleanly. Separate saturation adjustment and full cloud scheme controls are used to consider the strength of advection–condensation theory. Three sets of perturbations to the cloud scheme are designed to test the model’s sensitivity to the physics of condensation, sedimentation, and precipitation formation. The condensation and sedimentation perturbations isolate two key differences between the control cases. First, the sub-grid-scale relative humidity distribution assumed for the cloud macrophysics influences the location and magnitude of the extratropical cloud maxima, which interrupt the isentropic transport of moisture to the polar troposphere. Second, within the model’s explicit treatment of cloud microphysics, re-evaporation of hydrometeors moistens and increases clouds in the lower troposphere. In contrast, microphysical processes of precipitation formation (specifically, the ratio of accretion to autoconversion) have negligible effects on humidity, cloudiness, and precipitation apart from the strength of the large-scale condensation and formation cycle. In addition, counterintuitive relationships—such as cloud condensate and cloud fraction responding in opposing directions—emphasize the need for careful dissection of physical mechanisms. In keeping with advection–condensation theory, circulation sets the patterns of humidity, clouds, and precipitation to first order, with factors explored herein providing secondary controls. The results substantiate the utility of such idealized modeling and highlight key cloud processes to constrain.
We describe the model performance of a new global coupled climate model configuration, CM4-MG2. Beginning with the Geophysical Fluid Dynamics Laboratory's fourth-generation physical climate model (CM4.0), we incorporate a two-moment Morrison-Gettelman bulk stratiform microphysics scheme with prognostic precipitation (MG2), and a mineral dust and temperature-dependent cloud ice nucleation scheme. We then conduct and analyze a set of fully coupled atmosphere-ocean-land-sea ice simulations, following Coupled Model Intercomparison Project Phase 6 protocols. CM4-MG2 generally captures CM4.0's baseline simulation characteristics, but with several improvements, including better marine stratocumulus clouds off the west coasts of Africa and North and South America, a reduced bias toward “double” Intertropical Convergence Zones south of the equator, and a stronger Madden-Julian Oscillation (MJO). Some degraded features are also identified, including excessive Arctic sea ice extent and a stronger-than-observed El Nino-Southern Oscillation. Compared to CM4.0, the climate sensitivity is reduced by about 10% in CM4-MG2.
MacDonald, Cameron G., and Yi Ming, December 2022: Tropical intraseasonal variability response to zonally asymmetric forcing in an idealized moist GCM. Journal of Climate, 35(24), DOI:10.1175/JCLI-D-22-0344.14479-4501. Abstract
The tropical intraseasonal variability in an idealized moist general circulation model (GCM) coupled to a slab ocean is investigated. The model has a simple moist convection scheme and realistic radiative transfer, but no parameterization of cloud processes. In a zonally symmetric aquaplanet state, variability is dominated by westward-propagating Rossby waves. Enforcing zonal asymmetry through the application of a prescribed ocean heat flux in the bottom boundary leads to the development of a slow, eastward propagating mode that bears some of the characteristics of the observed Madden–Julian oscillation (MJO). When the ocean heat flux is made stronger, high-frequency Kelvin waves exist alongside the MJO mode. The strength of the disturbances and the spatial distribution of their precipitation anomalies are sensitive to the strength of intraseasonal sea surface temperature (SST) anomalies. The greatest resemblance to the MJO is observed when shallow slab ocean depths (1 m) are used, but the mode still exists at deeper slabs. Sensitivity experiments to the parameters of the convection scheme suggest that the simulated MJO mode couples to convection in a way that is distinct from both Kelvin and Rossby waves generated by the model. Analysis of the column moist static energy (CMSE) budget of the MJO mode suggests that column radiative heating plays only a weak role in destabilizing the mode relative to the stabilizing contribution of vertical advection. The CMSE budget analysis highlights the importance of the life cycle of horizontal advection for the destabilization and propagation of the MJO. Synergies between the generated MJO mode and linear theories of the MJO are discussed as well.
This study examines the climatology and dynamics of atmospheric blocking, and the general circulation features that influence blocks in GFDL’s atmosphere-only (AM4) and coupled atmosphere–ocean (CM4) comprehensive models. We compare AM4 and CM4 with reanalysis, focusing on winter in the Northern Hemisphere. Both models generate the correct blocking climatology and planetary-scale signatures of the stationary wave. However, at regional scales some biases exist. In the eastern Pacific and over western North America, both models generate excessive blocking frequency and too strong of a stationary wave. In the Atlantic, the models generate too little blocking and a weakened stationary wave. A block-centered compositing analysis of block-onset dynamics reveals that the models 1) produce realistic patterns of high-frequency (1–6-day) eddy forcing and 2) capture the notable differences in the 500-hPa geopotential height field between Pacific and Atlantic blocking. However, the models fail to reproduce stronger wave activity flux convergence in the Atlantic compared to the Pacific. Overall, biases in the blocking climatology in terms of location, frequency, duration, and area are quite similar between AM4 and CM4 despite the models having large differences in sea surface temperatures and climatological zonal circulation. This could suggest that other factors could be more dominant in generating blocking biases for these GCMs.
The characteristics of tropical mesoscale convective systems (MCSs) simulated with a finer-resolution (~50 km) version of the Geophysical Fluid Dynamics Laboratory (GFDL) AM4 model are evaluated by comparing with a comprehensive long-term observational dataset. It is shown that the model can capture the various aspects of MCSs reasonably well. The simulated spatial distribution of MCSs is broadly in agreement with the observations. This is also true for seasonality and interannual variability over different land and oceanic regions. The simulated MCSs are generally longer-lived, weaker, and larger than observed. Despite these biases, an event-scale analysis suggests that their duration, intensity, and size are strongly correlated. Specifically, longer-lived and stronger events tend to be bigger, which is consistent with the observations. The same model is used to investigate the response of tropical MCSs to global warming using time-slice simulations forced by prescribed sea surface temperatures and sea ice. There is an overall decrease in occurrence frequency, and the reduction over land is more prominent than over ocean.
A two-moment Morrison-Gettelman bulk cloud microphysics with prognostic precipitation (MG2), together with a mineral dust and temperature-dependent ice nucleation scheme, have been implemented into the Geophysical Fluid Dynamics Laboratory's Atmosphere Model version 4.0 (AM4.0). We refer to this configuration as AM4-MG2. This paper describes the configuration of AM4-MG2, evaluates its performance, and compares it with AM4.0. It is shown that the global simulations with AM4-MG2 compare favorably with observations and reanalyses. The model skill scores are close to AM4.0. Compared to AM4.0, improvements in AM4-MG2 include (a) better coastal marine stratocumulus and seasonal cycles, (b) more realistic ice fraction, and (c) dominant accretion over autoconversion. Sensitivity tests indicate that nucleation and sedimentation schemes have significant impacts on cloud liquid and ice water fields, but higher horizontal resolution (about 50 km instead of 100 km) does not.
The effect of stratospheric ozone depletion is simulated in GFDL AM4 model with three ozone schemes: Prescribing monthly zonal mean ozone concentration, full interactive stratospheric chemistry, and a simplified linear ozone chemistry scheme but with full dynamical interactions. While similar amounts of ozone loss are simulated by the three schemes, the two interactive ozone schemes produce significantly stronger stratospheric cooling than the prescribed one. We find that this temperature difference is driven by the dynamical responses to ozone depletion. In particular, the existence of ozone hole leads to strong ozone eddies that are in‐phase with the temperature eddies. The coherence between ozone and temperature anomalies leads to a weaker radiative damping as ozone absorbs shortwave radiation that compensates for the longwave cooling. As a result, less wave dissipates at the lower stratosphere, leading to a weaker descending and dynamical heating over the polar lower stratosphere, and hence a stronger net cooling there. The covariance between ozone and temperature is largely suppressed when ozone is prescribed as monthly zonal mean time series, as is the reduction in the radiative damping following ozone depletion. With much lower computational cost, the simplified ozone scheme is capable of producing similar magnitude of ozone loss and the consequent dynamical responses to those simulated by the full chemistry.
Loeb, Norman G., Wenying Su, Nicolas Bellouin, and Yi Ming, March 2021: Changes in clear-sky shortwave aerosol direct radiative effects since 2002. JGR Atmospheres, 126(5), DOI:10.1029/2020JD034090. Abstract
A new method for determining clear-sky shortwave aerosol direct radiative effects (ADRE) from the Clouds and the Earth's Radiant Energy System is used to examine changes in ADRE since 2002 alongside changes in aerosol optical depth (AOD) from the Moderate Resolution Spectroradiometer. At global scales, neither ADRE nor AOD show a significant trend. Over the northern hemisphere (NH), ADRE increases by 0.18 ± 0.17 Wm−2 per decade (less reflection to space) but shows no significant change over the southern hemisphere. The increase in the NH is primarily due to emission reductions in China, the United States, and Europe. The COVID-19 shutdown shows no noticeable impact on either global ADRE or AOD, but there is a substantial influence over northeastern China in March 2020. In contrast, February 2020 anomalies in ADRE and AOD are within natural variability even though the impact of the shutdown on industry was more pronounced in February than March. The reason is because February 2020 was exceptionally hot and humid over China, which compensated for reduced emissions. After accounting for meteorology and normalizing by incident solar flux, February ADRE anomalies increase substantially, exceeding the climatological mean ADRE by 23%. February and March 2020 correspond to the only period in which adjusted anomalies exceed the 95% confidence interval for 2 consecutive months. Distinct water-land differences over northeastern China are observed in ADRE but not in AOD. This is likely due to the influence of surface albedo on ADRE in the presence of absorbing aerosols.
The Coronavirus Disease 2019 (COVID‐19) pandemic led to a widespread reduction in aerosol emissions. Using satellite observations and climate model simulations, we study the underlying mechanisms of the large decreases in solar clear‐sky reflection (3.8 W m−2 or 7%) and aerosol optical depth (0.16 W m−2 or 32%) observed over the East Asian Marginal Seas in March 2020. By separating the impacts from meteorology and emissions in the model simulations, we find that about one‐third of the clear‐sky anomalies can be attributed to pandemic‐related emission reductions, and the rest to weather variability and long‐term emission trends. The model is skillful at reproducing the observed interannual variations in solar all‐sky reflection, but no COVID‐19 signal is discerned. The current observational and modeling capabilities will be critical for monitoring, understanding, and predicting the radiative forcing and climate impacts of the ongoing crisis.
Ramaswamy, V, Yi Ming, and M Daniel Schwarzkopf, April 2021: Forcing of global hydrological changes in the twentieth and twenty-first centuries In Hydrological Aspects of Climate Change [Pandey, A., S. Kumar, and A. Kumar (eds.)], Springer, Singapore, Springer Transactions in Civil and Environmental Engineering, DOI:10.1007/978-981-16-0394-561-76. Abstract
The Earth’s climate system in the twentieth century has experienced significant effects due to human-influenced factors. In this paper, we focus on the manner in which anthropogenic aerosols have radiatively forced changes in temperature and precipitation and contrast the effects with that due to the influence of well-mixed greenhouse gases. We employ the NOAA/Geophysical Fluid Dynamics Laboratory 3rd generation global climate model to simulate and derive a mechanistic understanding of the response to the forcings. We find that, over the twentieth century, anthropogenic aerosols have counteracted greenhouse gas effects to a substantial extent with regards to climate forcing, temperature and precipitation. The manner in which this comes about is traced through the effects on the atmosphere and surface heat balance, with resultant effects on the hydrologic cycle. Understanding of the twentieth century precipitation change is a prerequisite for confidence in model-based projections of the effects in the twenty-first century in response to emission scenarios of greenhouse gases and aerosols.
Smyth, Jane E., and Yi Ming, October 2021: Investigating the impact of land surface characteristics on monsoon dynamics with idealized model simulations and theories. Journal of Climate, 34(19), DOI:10.1175/JCLI-D-20-0954.17943-7958. Abstract
Monsoons emerge over a range of land surface conditions and exhibit varying physical characteristics over the seasonal cycle, from onset to withdrawal. Systematically varying the moisture and albedo parameters over land in an idealized modeling framework allows one to analyze the physics underlying the successive stages of monsoon development. To this end, we implement an isolated South American continent with reduced heat capacity but no topography in an idealized moist general circulation model. Irrespective of the local moisture availability, the seasonal cycles of precipitation and circulation over the South American monsoon sector are distinctly monsoonal with the default surface albedo. The dry land case (zero evaporation) is characterized by a shallow overturning circulation with vigorous lower-tropospheric ascent, transporting water vapor from the ocean. By contrast, with bucket hydrology or unlimited land moisture, the monsoon features deep moist convection that penetrates the upper troposphere. A series of land albedo perturbation experiments indicates that the monsoon strengthens with the net column energy flux and the near-surface moist static energy with all land moisture conditions. When the land–ocean thermal contrast is strong enough, inertial instability alone is sufficient for producing a shallow but vigorous circulation and converging a large amount of moisture from the ocean even in the absence of land moisture. Once the land is sufficiently moist, convective instability takes hold and the shallow circulation deepens. These results have implications for monsoon onset and intensification, and may elucidate the seasonal variations in how surface warming impacts tropical precipitation over land.
Zhang, Pengfei, Gang Chen, and Yi Ming, September 2021: Quantifying the mechanisms of atmospheric circulation response to greenhouse gas increases in a forcing–feedback framework. Journal of Climate, 34(17), DOI:10.1175/JCLI-D-20-0778.17005-7022. Abstract
While there is substantial evidence for tropospheric jet shift and Hadley cell expansion in response to greenhouse gas increases, quantitative assessments of individual mechanisms and feedback for atmospheric circulation changes remain lacking. We present a new forcing–feedback analysis on circulation response to increasing CO2 concentration in an aquaplanet atmospheric model. This forcing–feedback framework explicitly identifies a direct zonal wind response by holding the zonal mean zonal wind exerting on the zonal advection of eddies unchanged, in comparison with the additional feedback induced by the direct response in zonal mean zonal wind. It is shown that the zonal advection feedback accounts for nearly half of the changes to the eddy-driven jet shift and Hadley cell expansion, largely contributing to the subtropical precipitation decline, when the CO2 concentration varies over a range of climates. The direct response in temperature displays the well-known tropospheric warming pattern to CO2 increases, but the feedback exhibits negative signals. The direct response in eddies is characterized by a reduction in upward wave propagation and a poleward shift of midlatitude eddy momentum flux (EMF) convergence, likely due to an increase in static stability from moist thermodynamic adjustment. In contrast, the feedback features a dipole pattern in EMF that further shifts and strengthens midlatitude EMF convergence, resulting from the upper-level zonal wind increase seen in the direct response. Interestingly, the direct response produces an increase in eddy kinetic energy (EKE), but the feedback weakens EKE. Thus, the forcing–feedback framework highlights the distinct effect of zonal mean advecting wind from direct thermodynamic effects in atmospheric response to greenhouse gas increases.
Zhang, Pengfei, Gang Chen, Weiming Ma, Yi Ming, and Zheng Wu, September 2021: Robust atmospheric river response to global warming in idealized and comprehensive climate models. Journal of Climate, 34(18), DOI:10.1175/JCLI-D-20-1005.17717-7734. Abstract
Atmospheric rivers (ARs), narrow intense moisture transport, account for much of the poleward moisture transport in midlatitudes. While studies have characterized AR features and the associated hydrological impacts in a warming climate in observations and comprehensive climate models, the fundamental dynamics for changes in AR statistics (e.g., frequency, length, width) are not well understood. Here we investigate AR response to global warming with a combination of idealized and comprehensive climate models. To that end, we developed an idealized atmospheric GCM with Earth-like global circulation and hydrological cycle, in which water vapor and clouds are modeled as passive tracers with simple cloud microphysics and precipitation processes. Despite the simplicity of the model physics, it reasonably reproduces observed dynamical structures for individual ARs, statistical characteristics of ARs, and spatial distributions of AR climatology. Under climate warming, the idealized model produces robust AR changes similar to CESM large ensemble simulations under RCP8.5, including AR size expansion, intensified landfall moisture transport, and an increased AR frequency, corroborating previously reported AR changes under global warming by climate models. In addition, the latitude of AR frequency maximum shifts poleward with climate warming. Further analysis suggests that the thermodynamic effect (i.e., an increase in water vapor) dominates the AR statistics and frequency changes while both the dynamic and thermodynamic effects contribute to the AR poleward shift. These results demonstrate that AR changes in a warming climate can be understood as passive water vapor and cloud tracers regulated by large-scale atmospheric circulation, whereas convection and latent heat feedback are of secondary importance.
Zhou, Xiaoli, Rachel Atlas, Isabel L McCoy, Christopher S Bretherton, Charles Bardeen, Andrew Gettelman, Pu Lin, and Yi Ming, February 2021: Evaluation of cloud and precipitation simulations in CAM6 and AM4 using observations over the Southern Ocean. Earth and Space Science, 8(2), DOI:10.1029/2020EA001241. Abstract
This study uses cloud and radiative properties collected from in situ and remote sensing instruments during two coordinated campaigns over the Southern Ocean between Tasmania and Antarctica in January–February 2018 to evaluate the simulations of clouds and precipitation in nudged-meteorology simulations with the CAM6 and AM4 global climate models sampled at the times and locations of the observations. Fifteen SOCRATES research flights sampled cloud water content, cloud droplet number concentration, and particle size distributions in mixed-phase boundary layer clouds at temperatures down to −25°C. The 6-week CAPRICORN2 research cruise encountered all cloud regimes across the region. Data from vertically pointing 94 GHz radars deployed was compared with radar simulator output from both models. Satellite data were compared with simulated top-of-atmosphere (TOA) radiative fluxes. Both models simulate observed cloud properties fairly well within the variability of observations. Cloud base and top in both models are generally biased low. CAM6 overestimates cloud occurrence and optical thickness while cloud droplet number concentrations are biased low, leading to excessive TOA reflected shortwave radiation. In general, low clouds in CAM6 precipitate at the same frequency but are more homogeneous compared to observations. Deep clouds are better simulated but produce snow too frequently. AM4 underestimates cloud occurrence but overestimates cloud optical thickness even more than CAM6, causing excessive outgoing longwave radiation fluxes but comparable reflected shortwave radiation. AM4 cloud droplet number concentrations match observations better than CAM6. Precipitating low and deep clouds in AM4 have too little snow. Further investigation of these microphysical biases is needed for both models.
Atlas, Rachel, Christopher S Bretherton, Peter N Blossey, Andrew Gettelman, Charles Bardeen, Pu Lin, and Yi Ming, November 2020: How well do large-eddy simulations and global climate models represent observed boundary layer structures and low clouds over the summertime Southern Ocean?Journal of Advances in Modeling Earth Systems, 12(11), DOI:10.1029/2020MS002205. Abstract
Climate models struggle to accurately represent the highly reflective boundary layer clouds overlying the remote and stormy Southern Ocean. We use in situ aircraft observations from the Southern Ocean Clouds, Radiation and Aerosol Transport Experimental Study (SOCRATES) to evaluate Southern Ocean clouds in a cloud-resolving large-eddy simulation (LES) and two coarse resolution global atmospheric models, the CESM Community Atmosphere Model (CAM6) and the GFDL Atmosphere Model (AM4), run in a nudged hindcast framework. We develop six case studies from SOCRATES data which span the range of observed cloud and boundary layer properties. For each case, the LES is run once forced purely using reanalysis data (fifth generation European Centre for Medium-Range Weather Forecasts atmospheric reanalysis, “ERA5 based”) and once strongly nudged to an aircraft profile(“Obs based”). The ERA5-based LES can be compared with the global models, which are also nudged to reanalysis data and are better for simulating cumulus. The Obs-based LES closely matches an observed cloud profile and is useful for microphysical comparisons and sensitivity tests and simulating multilayer stratiform clouds. We use two-moment Morrison microphysics in the LES and find that it simulates too few frozen particles in clouds occurring within the Hallett-Mossop temperature range. We tweak the Hallett-Mossop parameterization so that it activates within boundary layer clouds, and we achieve better agreement between observed and simulated microphysics. The nudged global climate models (GCMs) simulate liquid-dominated mixed-phase clouds in the stratiform cases but excessively glaciate cumulus clouds. Both GCMs struggle to represent two-layer clouds, and CAM6 has low droplet concentrations in all cases and underpredicts stratiform cloud-driven turbulence.
Chua, Xin Rong, and Yi Ming, December 2020: Convective invigoration traced to warm-rain microphysics. Geophysical Research Letters, 47(23), DOI:10.1029/2020GL089134. Abstract
Aerosols are postulated to alter moist convection by increasing cloud droplet number concentration (Nd). Cloud-resolving model simulations of radiative-convective equilibrium show that higher Nd leads to stronger convective mass flux, seemingly in line with a hypothesis that links the convective invigoration to delayed rain formation allowing more cloud liquid condensate to be frozen. Yet, the invigoration is also present in an alternative model configuration with warm-rain microphysics only, suggesting that ice microphysics is not central to the phenomenon. The key dynamical mechanism lies in the different vertical distributions of the increases in water vapor condensation and in cloud liquid reevaporation, causing a dipole pattern favoring convection. This is further supported by a pair of mechanism-denial experiments in which an imposed weakening of cloud liquid reevaporation tends to mute invigoration.
In this paper, it is shown that westward-propagating monsoon-low-pressure-system-like disturbances in the South Asian monsoon region can be simulated in an idealized moist general circulation model through the addition of a simplified parameterization of land. Land is parameterized as having one-tenth the heat capacity of the surrounding slab ocean, with evaporation limited by a bucket hydrology model. In this model, the prominent topography of the Tibetan Plateau does not appear to be necessary for these storm systems to form or propagate; therefore focus is placed on the simulation with land but no topography.
The properties of the simulated storms are elucidated using regression analysis and compared to results from composites of storms from comprehensive GCMs in prior literature and reanalysis. The storms share a similar vertical profile in anomalous Ertel potential vorticity to those in reanalysis. Propagation, however, does not seem to be strongly dictated by beta-drift. Rather, it seems to be more closely consistent with linear moisture vortex instability theory, with the exception of the importance of the vertical advection term in the Ertel potential vorticity budget toward the growth and maintenance of disturbances. The results presented here suggest that a simplified GCM configuration might be able to be used to gain a clearer understanding of the sensitivity of monsoon low pressure systems to changes in the mean state climate.
Monsoon low-pressure systems (MLPSs) are among the most important synoptic-scale disturbances of the South Asian summer monsoon. Potential changes in their characteristics in a warmer climate would have broad societal impacts. Yet, the findings from a few existing studies are inconclusive. We use the Geophysical Fluid Dynamics Laboratory (GFDL) coupled climate model CM4.0 to examine the projected changes in the simulated MLPS activity under a future emission scenario. It is shown that CM4.0 can skillfully simulate the number, genesis location, intensity and lifetime of MLPSs. Global warming gives rise to a significant decrease in MLPS activity. An analysis of several large-scale environmental variables, both dynamic and thermodynamic, suggests that the decrease in MLPS activity can be attributed mainly to a reduction in low-level relative vorticity over the core genesis region. The decreased vorticity is consistent with weaker large-scale ascent, which leads to less vorticity production through the stretching term in the vorticity equation. Assuming a fixed radius of influence, the projected reduction in MLPSs would significantly lower the associated precipitation over the north central India, despite an overall increase in mean precipitation.
https://doi.org/10.1175/JCLI-D-20-0168.1
Kuo, Yi-Hung, J David Neelin, C-C Chen, W-T Chen, Leo J Donner, Andrew Gettelman, Xianan Jiang, K-T Kuo, Eric Maloney, C R Mechoso, Yi Ming, K A Schiro, Charles J Seman, Chien-Ming Wu, and Ming Zhao, January 2020: Convective transition statistics over tropical oceans for climate model diagnostics: GCM evaluation. Journal of the Atmospheric Sciences, 77(1), DOI:10.1175/JAS-D-19-0132.1. Abstract
To assess deep-convective parameterizations in a variety of GCMs and examine the fast-timescale convective transition, a set of statistics characterizing the pickup of precipitation as a function of column water vapor (CWV), PDFs and joint-PDFs of CWV and precipitation, and the dependence of the moisture-precipitation relation on tropospheric temperature is evaluated using the hourly output of two versions of GFDL AM4, NCAR CAM5 and superparameterized CAM (SPCAM). The 6-hourly output from the MJOTF/GASS project is also analyzed. Contrasting statistics produced from individual models that primarily differ in representations of moist convection suggest that convective transition statistics can substantially distinguish differences in convective representation and its interaction with the large-scale flow, while models that differ only in spatial-temporal resolution, microphysics, or ocean-atmosphere coupling result in similar statistics. Most of the models simulate some version of the observed sharp increase in precipitation as CWV exceeds a critical value, as well as that convective onset occurs at higher CWV but at lower column RH as temperature increases. While some models quantitatively capture these observed features and associated probability distributions, considerable intermodel spread and departures from observations in various aspects of the precipitation-CWV relationship are noted. For instance, in many of the models, the transition from the low-CWV, non-precipitating regime to the moist regime for CWV around and above critical is less abrupt than in observations. Additionally, some models overproduce drizzle at low CWV, and some require CWV higher than observed for strong precipitation. For many of the models, it is particularly challenging to simulate the probability distributions of CWV at high temperature.
Li, Xiaoqiong, Mingfang Ting, Yujia You, Dong-Eun Lee, Daniel M Westervelt, and Yi Ming, January 2020: South Asian summer monsoon response to aerosol‐forced sea surface temperatures. Geophysical Research Letters, 47(1), DOI:10.1029/2019GL085329. Abstract
Climate models suggest that anthropogenic aerosol‐induced drying dominates the historical rainfall changes over the heavily populated south Asian monsoon region. The regional response depends on both the aerosol fast radiative effect and the slow process through sea surface temperature (SST) cooling. Two atmospheric general circulation models, NCAR‐CAM5 and GFDL‐AM3, are used to investigate the monsoon response to prescribed aerosol‐forced SSTs. The total SST is separated into uniform cooling and a spatially‐varying component characterized by interhemispheric asymmetry. The monsoon rainfall is predominantly controlled by the non‐uniform SSTs, in the local Indian Ocean, South and East China Seas (IO‐CSs). The reduced meridional SST gradient in the IO‐CSs leads to weakened monsoon circulation, which drives a north‐south dipole rainfall change. The latitudinal location of the dipole shows model dependence due to differences in local SSTs and their meridional gradient, which determines the latitudinal location of the meridional overturning circulation responses.
Liu, Zhen, Yi Ming, Chun Zhao, and Ngar-Cheung Lau, et al., January 2020: Contribution of local and remote anthropogenic aerosols to a record-breaking torrential rainfall event in Guangdong Province, China. Atmospheric Chemistry and Physics, 20(1), DOI:10.5194/acp-20-223-2020. Abstract
A torrential rainfall case, which happened in Guangdong Province during 14–16 December 2013, broke the historical rainfall record in the province in terms of duration, affected area, and accumulative precipitation. The influence of anthropogenic aerosols on this extreme rainfall event is examined using a coupled meteorology–chemistry–aerosol model. Up to 33.7 mm precipitation enhancement in the estuary and near the coast is mainly attributed to aerosol–cloud interactions (ACI), whereas aerosol–radiation interaction partially offsets 14 % of the precipitation increase. Our further analysis of changes in hydrometeors and latent heat sources suggests that the ACI effects on the intensification of precipitation can be divided into two stages: cold rain enhancement in the former stage followed by warm rain enhancement in the latter. Responses of precipitation to the changes in anthropogenic aerosol concentration from local (i.e., Guangdong Province) and remote (i.e., outside Guangdong Province) sources are also investigated through simulations with reduced aerosol emissions from either local or remote sources. Accumulated aerosol concentration from local sources aggregates mainly near the ground surface and dilutes quickly after the precipitation initiated. By contrast, the aerosols from remote emissions extend up to 8 km above ground and last much longer before decreasing until peak rainfall begins, because aerosols are continuously transported by the strong northerly winds. The patterns of precipitation response to remote and local aerosol concentrations resemble each other. However, compared with local aerosols through warm rain enhancement, remote aerosols contribute more than twice the precipitation increase by intensifying both cold and warm rain, occupying a predominant role. A 10-time emission sensitivity test shows about 10 times the PM2.5 concentration compared with the control run. Cold (warm) rain is drastically enhanced (suppressed) in the 10× run. In response to 10× aerosol emissions, the pattern of precipitation and cloud property changes resembles the differences between CTL and CLEAN, but with a much greater magnitude. The precipitation average over Guangdong decreases by 1.0 mm in the 10× run but increases by 1.4 mm in the control run compared with the CLEAN run. We note that the precipitation increase is concentrated within a more narrowed downstream region of the aerosol source, whereas the precipitation decrease is more dispersed across the upstream region. This indicates that the excessive aerosols not only suppress rainfall, but also change the spatial distribution of precipitation, increasing the rainfall range, thereby potentially exacerbating flood and drought elsewhere. This study highlights the importance of considering aerosols in meteorology to improve extreme weather forecasting. Furthermore, aerosols from remote emissions may outweigh those from local emissions in the convective invigoration effect.
Loeb, Norman G., Hailan Wang, Richard P Allan, Timothy Andrews, Kyle Armour, Jason N S Cole, J-L Dufresne, Piers M Forster, Andrew Gettelman, Huan Guo, T Mauritsen, Yi Ming, and David J Paynter, et al., March 2020: New Generation of Climate Models Track Recent Unprecedented Changes in Earth's Radiation Budget Observed by CERES. Geophysical Research Letters, 47(5), DOI:10.1029/2019GL086705. Abstract
We compare top‐of‐atmosphere (TOA) radiative fluxes observed by the Clouds and the Earth's Radiant Energy System (CERES) and simulated by seven general circulation models forced with observed sea‐surface temperature (SST) and sea‐ice boundary conditions. In response to increased SSTs along the equator and over the eastern Pacific (EP) following the so‐called global warming “hiatus” of the early 21st century, simulated TOA flux changes are remarkably similar to CERES. Both show outgoing shortwave and longwave TOA flux changes that largely cancel over the west and central tropical Pacific, and large reductions in shortwave flux for EP low‐cloud regions. A model's ability to represent changes in the relationship between global mean net TOA flux and surface temperature depends upon how well it represents shortwave flux changes in low‐cloud regions, with most showing too little sensitivity to EP SST changes, suggesting a “pattern effect” that may be too weak compared to observations.
Many fundamental questions remain about the roles and effects of stationary forcing on atmospheric blocking. As such, this work utilizes an idealized moist general circulation model (GCM) to investigate atmospheric blocking in terms of dynamics, geographical location, and duration. The model is first configured as an aquaplanet, then orography is added in separate integrations. Block-centered composites of wave activity fluxes and height show that blocks in the aquaplanet undergo a realistic dynamical evolution when compared to reanalysis. Blocks in the aquaplanet are also found to have similar life cycles to blocks in model integrations with orography. These results affirm the usefulness of both zonally symmetric and asymmetric idealized model configurations for studying blocking. Adding orography to the model leads to an increase in blocking. This mirrors what is observed when comparing the Northern Hemisphere (NH) and Southern Hemisphere (SH), where the NH contains more orography and thus more blocking. As the prescribed mountain height increases, so do the magnitude and size of climatological stationary waves, resulting in more blocking overall. Increases in blocking, however, are not spatially uniform. Orography is found to induce regions of enhanced block frequency just upstream of mountains, near high pressure anomalies in the stationary waves, which is poleward of climatological minima in upper-level zonal wind, while block frequency minima and jet maxima occur eastward of the wave trough. This result matches what is observed near the Rocky Mountains. Finally, an analysis of block duration suggests blocks generated near stationary wave maxima last slightly longer than blocks that form far from or without orography. Overall, the results of this work help to explain some of the observed similarities and differences in blocking between the NH and SH and emphasize the importance of general circulation features in setting where blocks most frequently occur.
Anthropogenic aerosols have been postulated to have a cooling effect on climate, but its magnitude remains uncertain. Using atmospheric general circulation model simulations, we separate the land temperature response into a fast response to radiative forcings and a slow response to changing oceanic conditions and find that the former accounts for about one fifth of the observed warming of the Northern Hemisphere land during summer and autumn since the 1960s. While small, this fast response can be constrained by observations. Spatially varying aerosol effects can be detected on the regional scale, specifically warming over Europe and cooling over Asia. These results provide empirical evidence for the important role of aerosols in setting regional land temperature trends and point to an emergent constraint that suggests strong global aerosol forcing and high transient climate response.
Smyth, Jane E., and Yi Ming, November 2020: Characterizing drying in the South American monsoon onset season with the moist static energy budget. Journal of Climate, 33(22), DOI:10.1175/JCLI-D-20-0217.19735-9748. Abstract
The tropical atmospheric circulation and attendant rainfall exhibit seasonally dependent responses to increasing temperatures. Understanding changes in the South American monsoon system is of particular interest given the sensitivity of the southern Amazon rainforest to changes in dry season length. We utilize the latest Geophysical Fluid Dynamics Laboratory Atmospheric Model (GFDL AM4) to analyze the response of the South American monsoon to uniform sea surface temperature (SST) warming. SST warming is a poorly understood yet impactful component of greenhouse gas–induced climate change. Region-mean rainfall declines by 11%, and net precipitation (precipitation minus evaporation) declines by 40%, during the monsoon onset season (September–November), producing a more severe dry season. The column-integrated moist static energy (MSE) budget helps elucidate the physical mechanisms of the simulated drying. Based on the seasonal analysis, precipitation reductions tend to occur when 1) a convecting region’s climatological MSE export is dominated by horizontal rather than vertical advection, and 2) the horizontal MSE advection increases in the perturbed climate, impeding ascent. On a synoptic scale, the South American low-level jet strengthens and exports more moisture from the monsoon sector, exacerbating spring drying.
Westervelt, Daniel M., Yujia You, Xiaoqiong Li, Mingfang Ting, Dong-Eun Lee, and Yi Ming, July 2020: Relative importance of greenhouse gases, sulfate, organic carbon, and black carbon aerosol for South Asian monsoon rainfall changes. Geophysical Research Letters, 47(13), DOI:10.1029/2020GL088363. Abstract
The contribution of individual aerosol species and greenhouse gases to precipitation changes during the South Asian summer monsoon is uncertain. Mechanisms driving responses to anthropogenic forcings need further characterization. We use an atmosphere‐only climate model to simulate the fast response of the summer monsoon to different anthropogenic aerosol types and to anthropogenic greenhouse gases. Without normalization, sulfate is the largest driver of precipitation change between 1850 and 2000, followed by black carbon and greenhouse gases. Normalized by radiative forcing, the most effective driver is black carbon. The precipitation and moisture budget responses to combinations of aerosol species perturbed together scale as a linear superposition of their individual responses. We use both a circulation‐based and moisture budget‐based argument to identify mechanisms of aerosol and greenhouse gas induced changes to precipitation, and find that in all cases the dynamic contribution is the dominant driver to precipitation change in the monsoon region.
Adames, A F., D Kim, Spencer K Clark, Yi Ming, and Kuniaki Inoue, December 2019: Scale analysis of moist thermodynamics in a simple model and the relationship between moisture modes and gravity waves. Journal of the Atmospheric Sciences, 76(12), DOI:10.1175/JAS-D-19-0121.1. Abstract
Observations and theory of convectively-coupled equatorial waves suggest that they can be categorized into two distinct groups. Moisture modes are waves whose thermodynamics are governed by moisture fluctuations. The thermodynamics of the gravity wave group, on the other hand, are rooted in buoyancy (temperature) fluctuations. On the basis of scale analysis it is found that a simple nondimensional parameter –akin to the Rossby number– can explain the processes that lead to the existence of these two groups. This parameter, defined as Nmode, indicates that moisture modes arise when anomalous convection lasts sufficiently long so that dry gravity waves eliminate the temperature anomalies in the convective region, satisfying weak temperature gradient (WTG) balance. This process causes moisture anomalies to dominate the distribution of moist enthalpy (or moist static energy), and hence the evolution of the wave. Conversely, convectively-coupled gravity waves arise when anomalous convection eliminates the moisture anomalies more rapidly than dry gravity waves can adjust the troposphere towards WTG balance, causing temperature to govern the moist enthalpy distribution and evolution. Spectral analysis of reanalysis data indicates that slowly-propagating waves (cp ~ 3 m s-1) are likely to be moisture modes while fast waves (cp ~ 30 m s-1) exhibit gravity wave behavior, with “mixed moisture-gravity” waves existing in between. While these findings are obtained from a highly idealized framework, it is hypothesized that they can be extended to understand simulations of convectively-coupled waves in GCMs and the thermodynamics of more complex phenomena.
Chua, Xin Rong, Yi Ming, and Nadir Jeevanjee, August 2019: Investigating the Fast Response of Precipitation Intensity and Boundary Layer Temperature to Atmospheric Heating Using a Cloud‐Resolving Model. Geophysical Research Letters, 46(15), DOI:10.1029/2019GL082408. Abstract
Coarse‐resolution global climate models cannot explicitly resolve the intensity distribution of tropical precipitation and how it responds to a forcing. We use a cloud‐resolving model to study how imposed atmospheric radiative heating (such as that caused by greenhouse gases or absorbing aerosols) may alter precipitation intensity in the setting of radiative‐convective equilibrium. It is found that the decrease in total precipitation is realized through reducing weak events. The intensity of strong precipitation events is maintained by a cancellation between the moistening of air parcels and weakening of updrafts. A boundary layer energy budget analysis suggests that free‐tropospheric heating raises boundary layer temperatures mainly through a reduction in rain re‐evaporation. This insight leads to a predictive scaling for the surface sensible and latent flux changes. The results imply that cloud microphysical processes play a key role in shaping the temperature and precipitation responses to atmospheric heating.
We describe GFDL's CM4.0 physical climate model, with emphasis on those aspects that may be of particular importance to users of this model and its simulations. The model is built with the AM4.0/LM4.0 atmosphere/land model and OM4.0 ocean model. Topics include the rationale for key choices made in the model formulation, the stability as well as drift of the pre‐industrial control simulation, and comparison of key aspects of the historical simulations with observations from recent decades. Notable achievements include the relatively small biases in seasonal spatial patterns of top‐of‐atmosphere fluxes, surface temperature, and precipitation; reduced double Intertropical Convergence Zone bias; dramatically improved representation of ocean boundary currents; a high quality simulation of climatological Arctic sea ice extent and its recent decline; and excellent simulation of the El Niño‐Southern Oscillation spectrum and structure. Areas of concern include inadequate deep convection in the Nordic Seas; an inaccurate Antarctic sea ice simulation; precipitation and wind composites still affected by the equatorial cold tongue bias; muted variability in the Atlantic Meridional Overturning Circulation; strong 100 year quasi‐periodicity in Southern Ocean ventilation; and a lack of historical warming before 1990 and too rapid warming thereafter due to high climate sensitivity and strong aerosol forcing, in contrast to the observational record. Overall, CM4.0 scores very well in its fidelity against observations compared to the Coupled Model Intercomparison Project Phase 5 generation in terms of both mean state and modes of variability and should prove a valuable new addition for analysis across a broad array of applications.
Lee, J E., Baylor Fox-Kemper, C Horvat, and Yi Ming, October 2019: The response of East Asian monsoon to the precessional cycle: A new study using the Geophysical Fluid Dynamics Laboratory model. Geophysical Research Letters, 46(20), DOI:10.1029/2019GL082661. Abstract
Speleothem oxygen isotopes have been shown to exhibit a close relationship with summer insolation in the northern hemisphere, leading to the hypothesis that East Asian monsoon intensity is proportional to the summer insolation. This hypothesis, however, has been questioned because previous climate model simulations have been unable to simulate the observed large variation in precipitation or the precipitation isotope values, about a half of the variation in the entire modern tropical regions, in response to the insolation change due to the precession cycle. Here we show new results, using the fully‐coupled Geophysical Fluid Dynamics Laboratory (GFDL) model, that it is dynamically possible to have much higher precipitation during the high summer insolation period compared with the low summer insolation period in the East Asian monsoon region. We conclude that past East Asian monsoon intensity probably increased with increasing northern hemispheric insolation, given a large change in speleothem oxygen isotopes.
Li, Y, Yi Deng, S Yang, Henian Zhang, Yi Ming, and Zhaoyi Shen, June 2019: Multi-scale temporal-spatial variability of the East Asian summer monsoon frontal system: observation versus its representation in the GFDL HiRAM. Climate Dynamics, 52(11), DOI:10.1007/s00382-018-4546-z. Abstract
This study examines the representation of the multi-scale temporospatial variability of the East Asian summer monsoon stationary front (MSF) in the High-Resolution Atmospheric Model (HiRAM) of the National Oceanic and Atmospheric Administration (NOAA) Geophysical Fluid Dynamics Laboratory. Compared with the observed variability of the MSF in the European Center for Medium-Range Weather Forecasts Reanalysis Interim (ERA-Interim), HiRAM reproduces reasonably well the seasonal mean precipitation pattern and the seasonal migration of MSF. However, wet biases are found over the northern and eastern China and northern Japan, and dry biases extend from the southern China to the western North Pacific. These rainfall biases are directly tied to a northwestward bias in the model simulated seasonal mean location of MSF and this location bias is most pronounced in the month of May. In general, the MSF in HiRAM is more intense, located more northwestward, and more stationary with weaker interannual variations compared to the observed. A pronounced positive bias in the ocean-land sea level pressure contrast over East Asia, largely manifested as the westward expansion of the western North Pacific subtropical high, is hypothesized to be the main cause of the northwestward location bias of MSF in HiRAM. This bias in sea level pressure contrast likely results from the missing of realistic air-sea interactions in the HiRAM simulations.
An unprecedented disruption of the Quasi-Biennial Oscillation (QBO) started to develop from late 2015. The early development of this event is analyzed using the space-time spectra of eddies from reanalysis data. While the extratropical waves propagating horizontally into the tropics were assumed to be the main driver for the disruption, it was not clear why these waves dissipated near the jet core instead of jet edge as linear theory predicts. This study shows that the drastic deceleration of the equatorial jet was largely brought about by a single strong wave packet, and the local winds experienced by the wave packet served as a better indicator of the wave breaking latitude than the zonal mean winds.
Surprisingly, tropical mixed Rossby gravity waves also made an appreciable contribution to the deceleration of the equatorial westerly jet by the horizontal eddy momentum fluxes, especially before January 2016. The horizontal eddy momentum fluxes associated with the tropical waves arise from the deformation of the wave structure when background westerlies increase with height. These horizontal eddy momentum anomalies from the tropical waves are commonly observed in the reanalysis data, but are typically much weaker than those in the 2015/2016 winter. The possibility exists that exceptionally strong equatorially trapped waves precondition the flow to disruption by an extratropical disturbance.
Liu, Zhen, Yi Ming, Lin Wang, Massimo Bollasina, M Luo, and Ngar-Cheung Lau, et al., August 2019: A Model Investigation of Aerosol‐Induced Changes in the East Asian Winter Monsoon. Geophysical Research Letters, 46(16), DOI:10.1029/2019GL084228. Abstract
The response of the East Asian winter monsoon (EAWM) circulation to aerosols is studied using a coupled atmosphere‐slab ocean general circulation model. In the extratropics, the aerosol‐induced cooling in the mid‐latitudes leads to an intensified subtropical jet stream, a deepened East Asian trough, and thus an enhanced EAWM. In the tropics, the local Hadley circulation shifts southward to compensate for the interhemispheric asymmetry in aerosol radiative cooling. Anomalous subsidence at around 10°N leads to a salient anticyclone to the southwest of the Philippines. The associated southwesterlies advect abundant moisture to South China, resulting in local precipitation increase and suggesting a weaker EAWM. The EAWM response to aerosol forcing is thus driven by a competition between tropical and extratropical mechanisms, which has important implications for the future monsoon evolution as aerosol changes may follow different regional‐dependent trajectories.
Maloney, Eric, Andrew Gettelman, Yi Ming, J David Neelin, D Barrie, Annarita Mariotti, C-C Chen, D B Coleman, Yi-Hung Kuo, B Singh, H Annamalai, Alexis Berg, James F Booth, Suzana J Camargo, A Dai, A Gonzalez, J Hafner, Xianan Jiang, X Jing, D Kim, Arun Kumar, Yumin Moon, C M Naud, Adam H Sobel, K Suzuki, F Wang, J Wang, Allison A Wing, X Xu, and Ming Zhao, September 2019: Process-Oriented Evaluation of Climate and Weather Forecasting Models. Bulletin of the American Meteorological Society, 100(9), DOI:10.1175/BAMS-D-18-0042.1. Abstract
Outcomes of NOAA MAPP Model Diagnostics Task Force activities to promote process-oriented diagnosis of models to accelerate development are described.
Realistic climate and weather prediction models are necessary to produce confidence in projections of future climate over many decades and predictions for days to seasons. These models must be physically justified and validated for multiple weather and climate processes. A key opportunity to accelerate model improvement is greater incorporation of process-oriented diagnostics (PODs) into standard packages that can be applied during the model development process, allowing the application of diagnostics to be repeatable across multiple model versions and used as a benchmark for model improvement. A POD characterizes a specific physical process or emergent behavior that is related to the ability to simulate an observed phenomenon. This paper describes the outcomes of activities by the Model Diagnostics Task Force (MDTF) under the NOAA Climate Program Office (CPO) Modeling, Analysis, Predictions and Projections (MAPP) program to promote development of PODs and their application to climate and weather prediction models. MDTF and modeling center perspectives on the need for expanded process-oriented diagnosis of models are presented. Multiple PODs developed by the MDTF are summarized, and an open-source software framework developed by the MDTF to aid application of PODs to centers’ model development is presented in the context of other relevant community activities. The paper closes by discussing paths forward for the MDTF effort and for community process-oriented diagnosis.
The clouds in southern hemisphere extratropical cyclones generated by the GFDL climate model are analyzed against MODIS, CloudSat and CALIPSO cloud and precipitation observations. Two model versions are used: one is a developmental version of AM4, a model GFDL will utilize for CMIP6, the other is the same model with a different parameterization of moist convection. Both model versions predict a realistic top-of-atmosphere cloud cover in the southern oceans, within 5% of the observations. However, an examination of cloud cover transects in extratropical cyclones reveals a tendency in the models to overestimate high-level clouds (by differing amounts) and underestimate cloud cover at low-levels (again by differing amounts), especially in the post-cold frontal (PCF) region, when compared to observations. Focusing on only the models, their differences in high and mid-level clouds are consistent with their differences in convective activity and relative humidity (RH), but the same is not true for the PCF region. In this region, RH is higher in the model with less cloud fraction. These seemingly contradictory cloud and RH differences can be explained by differences in the cloud parameterization tuning parameters that ensure radiative balance. In the PCF region, the model cloud differences are smaller than either of the model biases with respect to observations, suggesting other physics changes are needed to address the bias. The process-oriented analysis used to assess these model differences will soon be automated and shared.
Surface layer (SL) variables [e.g., 2‐m temperature (T2) and 10‐m wind (U10)] are diagnosed by applying the flux‐profile relationships based on Monin‐Obukhov similarity theory to the lowest model height (LMH). This assumes that the LMH is in the SL, which is approximately the bottom 10% of the boundary layer, but atmospheric general circulation models rarely satisfy this in stable boundary layers (SBLs). To assess errors in the diagnostic variables due to the LMH solely linked to the diagnostic algorithm, offline tests of the flux‐profile relationships are performed with LMH from a few meters to 60 m for three SBL regimes: weakly stable, very stable, and transition stability regimes. The results show that T2 and U10 are underestimated by O(0.1–1 °C) and O(0.1–1 m s−1), respectively, if the LMH is higher than the SL height. The stronger the SL stability is, the larger the temperature biases are. The negative wind biases increase with the surface stress. Based on these findings, we analyze the impacts of the LMH on the climatologies of the diagnostic parameters in the GFDL AM4.0/LM4.0. The results show reduced negative biases in T2 and U10 by lowering the LMH. The decrease of the overall bias over land is mainly due to the sensitivity of the diagnostic method to the LMH in SBLs, as shown in the offline tests. The overall increase in T2 and U10 over the oceans results from the increase in the actual near‐surface temperature and wind rather than from the diagnostic method.
Xie, Yongkun, Jianping Huang, and Yi Ming, April 2019: Robust Regional Warming Amplifications Directly Following the Anthropogenic Emission. Earth's Future, 7(4), DOI:10.1029/2018EF001068. Abstract
Amplified regional warming suggests that some regions suffer more from global warming. However, present knowledge is still not enough to understand this paramount aspect of global warming. By using climate model simulations, we first highlight that regional warming amplifications directly following anthropogenic emission would torment many undeveloped countries. Then, we find that cloud has the largest contribution to polar amplification, rather than the frequently blamed sea ice decline. The summertime sea ice decline plays the second role in generating strong wintertime polar amplification via lagged energy release. Our findings also propose that heat uptake by the deeper ocean greatly favors the terrestrial amplification. Additionally, the cloud also shows its importance in helping dryland amplification. Our findings suggest that these robust regional warming amplifications are very likely inevitable under anthropogenic emission, therefore strictest global mitigation is required to prevent unacceptable warming over the high risk regions.
Adames, A F., and Yi Ming, June 2018: Interactions between water vapor and potential vorticity in synoptic-scale monsoonal disturbances: Moisture vortex instability. Journal of the Atmospheric Sciences, 75(6), DOI:10.1175/JAS-D-17-0310.1. Abstract
South Asian monsoon low pressure systems, referred to as synoptic-scale monsoonal disturbances (SMDs), are convectively-coupled cyclonic disturbances that are responsible for up to half of the total monsoon rainfall. In spite of their importance, the mechanisms that lead to the growth of these systems have remained elusive. It has long been thought that SMDs grow due to a variant of baroclinic instability that includes the effects of convection. Recent work, however, has shown that this framework is inconsistent with the observed structure and dynamics of SMDs.
Here we present an alternative framework that may explain the growth of SMDs and may also be applicable to other modes of tropical variability. Moisture is prognostic and is coupled to precipitation through a simplified Betts-Miller scheme. Interactions between moisture and potential vorticity (PV ) in the presence of a moist static energy gradient can be understood in terms of a “gross” PV (qG) equation. qG summarizes the dynamics of SMDs and reveals the relative role that moist and dry dynamics play in these disturbances, which is largely determined by the gross moist stability.
Linear solutions to the coupled PV and moisture equations reveal Rossby-like modes that grow due to a moisture-vortex instability. Meridional temperature and moisture advection to the west of the PV maximum moisten and destabilize the column, which results in enhanced convection and SMD intensification through vortex stretching. This instability occurs only if the moistening is in the direction of propagation of the SMD and is strongest at the synoptic scale.
Adames, A F., and Yi Ming, June 2018: Moisture and moist static energy budgets of South Asian monsoon low pressure systems in GFDL AM4.0. Journal of the Atmospheric Sciences, 75(6), DOI:10.1175/JAS-D-17-0309.1. Abstract
The mechanisms that lead to the propagation of anomalous moisture and moist static energy (MSE) in monsoon low and high pressure systems, collectively referred to as synoptic-scale monsoonal disturbances (SMDs), are investigated using daily output fields from GFDL's atmospheric general circulation model (AM4.0). On the basis of linear regression analysis of westward-propagating rainfall anomalies of timescales shorter than 15 days, it is found that SMDs are organized into wavetrains of 3-4 individual cyclones and anticyclones. These events amplify over the Bay of Bengal, reach a maximum amplitude over the eastern coast of India and dissipate as they approach the Arabian Sea. The structure and propagation of the simulated SMDs resemble those documented in observations.
It is found that moisture and MSE anomalies exhibit similar horizontal structures in the simulated SMDs, indicating that moisture is the leading contributor to MSE. Propagation of the moisture anomalies is governed by vertical moisture advection while the MSE anomalies propagate due to horizontal advection of dry static energy by the anomalous winds. By combining the budgets, we interpret the propagation of the moisture anomalies in terms of lifting that is forced by horizontal dry static energy advection, that is, ascent along sloping isentropes. This process moistens the lower free-troposphere, producing an environment that is more favorable to deep convection. Ascent driven by radiative heating is of primary importance to the maintenance of the moisture anomalies.
In comprehensive and idealized general circulation models, hemispherically asymmetric forcings lead to shifts in the latitude of the Intertropical Convergence Zone (ITCZ). Prior studies using comprehensive GCMs (with complicated parameterizations of radiation, clouds, and convection) suggest that the water vapor feedback tends to amplify the movement of the ITCZ in response to a given hemispherically asymmetric forcing, but this effect has yet to be elucidated in isolation. This study uses an idealized moist model, coupled to a full radiative transfer code, but without clouds, to examine the role of the water vapor feedback in a targeted manner.
In experiments with interactive water vapor and radiation, the ITCZ latitude shifts roughly twice as much off the equator as in cases with the water vapor field seen by the radiation code prescribed to a static hemisperically-symmetric control distribution. Using energy flux equator theory for the latitude of the ITCZ, the amplification of the ITCZ shift is attributed primarily to the longwave water vapor absorption associated with the movement of the ITCZ into the warmer hemisphere, further increasing the net column heating asymmetry. Local amplification of the imposed forcing by the shortwave water vapor feedback plays a secondary role. Experiments varying the convective relaxation time, an important parameter in the convection scheme used in the idealized moist model, yield qualitatively similar results, suggesting some degree of robustness to the model physics; however, the sensitivity experiments do not preclude that more extreme modifications to the convection scheme could lead to qualitatively different behavior.
Dong, Wenhao, Yanluan Lin, Jonathon S Wright, Yuanyu Xie, and Yi Ming, et al., August 2018: Regional disparities in warm season rainfall changes over arid eastern-central Asia. Scientific Reports, 8, 13051, DOI:10.1038/s41598-018-31246-3. Abstract
Multiple studies have reported a shift in the trend of warm season rainfall over arid eastern–central Asia (AECA) around the turn of the new century, from increasing over the second half of the twentieth century to decreasing during the early years of the twenty-first. Here, a closer look based on multiple precipitation datasets reveals important regional disparities in these changes. Warm-season rainfall increased over both basin areas and mountain ranges during 1961–1998 due to enhanced moisture flux convergence associated with changes in the large-scale circulation and increases in atmospheric moisture content. Despite a significant decrease in warm-season precipitation over the high mountain ranges after the year 1998, warm season rainfall has remained large over low-lying basin areas. This discrepancy, which is also reflected in changes in river flow, soil moisture, and vegetation, primarily results from disparate responses to enhanced warming in the mountain and basin areas of AECA. In addition to changes in the prevailing circulation and moisture transport patterns, the decrease in precipitation over the mountains has occurred mainly because increases in local water vapor saturation capacity (which scales with temperature) have outpaced the available moisture supply, reducing relative humidity and suppressing precipitation. By contrast, rainfall over basin areas has been maintained by accelerated moisture recycling driven by rapid glacier retreat, snow melt, and irrigation expansion. This trend is unsustainable and is likely to reverse as these cryospheric buffers disappear, with potentially catastrophic implications for local agriculture and ecology.
How the globally uniform component of sea surface temperature (SST) warming influences rainfall in the African Sahel remains under-studied, despite mean SST warming being among the most robustly simulated and theoretically grounded features of anthropogenic climate change. A prior study using the NOAA Geophysical Fluid Dynamics Laboratory (GFDL) AM2.1 atmospheric general circulation model (AGCM) demonstrated that uniform SST warming strengthens the prevailing northerly advection of dry Saharan air into the Sahel. The present study uses uniform SST warming simulations performed with seven GFDL and ten CMIP5 AGCMs to assess the robustness of this drying mechanism across models and uses observations to assess the physical credibility of the severe drying response in AM2.1.
In all seventeen AGCMs, mean SST warming enhances the free-tropospheric meridional moisture gradient spanning the Sahel and with it the Saharan dry air advection. Energetically, this is partially balanced by anomalous subsidence, yielding decreased precipitation in fourteen of the seventeen models. Anomalous subsidence and precipitation are tightly linked across the GFDL models but not the CMIP5 models, precluding the use of this relationship as the start of a causal chain ending in an emergent observational constraint. For AM2.1, cloud-rainfall covariances generate radiative feedbacks on drying through the subsidence mechanism and through surface hydrology that are excessive compared to observations at the interannual timescale. These feedbacks also act in the equilibrium response to uniform warming, calling into question the Sahel’s severe drying response to warming in all coupled models using AM2.1.
This paper introduces an idealized general circulation model (GCM) in which water vapor and clouds are tracked as tracers, but are not allowed to affect circulation either through latent heat release or cloud radiative effects. The cloud scheme includes an explicit treatment of cloud microphysics and diagnoses cloud fraction from a prescribed sub-grid distribution of total water. The model is capable of qualitatively capturing many large-scale features of water vapor and cloud distributions outside of the boundary layer and deep tropics. The subtropical dry zones, mid-latitude storm tracks and upper-tropospheric cirrus are simulated reasonably well. The inclusion of cloud microphysics (namely rain re-evaporation) has a modest but significant effect of moistening the lower troposphere in this model. When being subjected to a uniform fractional increase of saturated water vapor pressure, the model produces little change in cloud fraction. A more realistic perturbation, which considers the non-linearity of the Clausius-Clapeyron relation and spatial structure of CO2-induced warming, results in a substantial reduction in the free-tropospheric cloud fraction. This is reconciled with an increase of relative humidity by analyzing the probability distributions of both quantities, and may help explain partly similar decreases in cloud fraction in full GCMs. The model provides a means to isolate individual processes or model components for studying their influences on cloud simulation in the extratropical free troposphere.
Despite distinct geographic distributions of top-of-the-atmosphere radiative forcing, anthropogenic greenhouse gases and aerosols have been found to produce similar patterns of climate response in atmosphere-and-ocean coupled climate model simulations. Understanding surface energy flux changes, a crucial pathway by which atmospheric forcing is communicated to the ocean, is a vital bridge to explaining the similar full atmosphere-and-ocean responses to these disparate forcings. Here we analyze the fast, atmosphere-driven change in surface energy flux caused by present-day greenhouse gases vs aerosols to elucidate its role in shaping the subsequent slow, coupled response. We find that the surface energy flux response patterns achieve roughly two-thirds of the anti-correlation seen in the fully coupled response, driven by Rossby waves excited by symmetric changes to the land–sea contrast. Our results suggest that atmosphere and land surface processes are capable of achieving substantial within-hemisphere homogenization in the climate response to disparate forcers on fast, societally-relevant timescales.
This study examines how aerosol absorption affects the extratropical circulation by analyzing the response to a globally uniform increase in black carbon (BC) simulated with an atmospheric general circulation model forced by prescribed sea surface temperatures. The model includes aerosol direct and semidirect effects, but not indirect or cloud-absorption effects. BC-induced heating in the free troposphere stabilizes the mid-latitude atmospheric column, which results in less energetic baroclinic eddies and thus reduced meridional energy transport at mid-latitudes. Upper tropospheric BC also decreases the meridional temperature gradient on the equatorward flank of the tropospheric jet and yields a weakening and poleward shift of the jet, while boundary layer BC has no significant influence on the large-scale circulation since most of the heating is diffused by turbulence in the boundary layer. The effectiveness of BC in altering circulation generally increases with height.
Dry baroclinic eddy theories can explain most of the extratropical response to free troposphere BC. Specifically, the decrease in vertical eddy heat flux related to a more stable atmosphere is the main mechanism for re-establishing atmospheric energy balance in the presence of BC-induced heating. Similar temperature responses are found in a dry idealized model, which further confirms the dominant role of baroclinic eddies in driving the extratropical circulation changes. The strong atmospheric-only response to BC suggests that absorbing aerosols are capable of altering synopic-scale weather patterns. Its height dependence highlights the importance of better constraining model-simulated aerosol vertical distributions with satellite and field measurements.
This study describes the performance of two Geophysical Fluid Dynamics Laboratory (GFDL) atmospheric general circulation models (AGCMs) in simulating the climatologies of planetary boundary layer (PBL) parameters, with a particular focus on the diurnal cycles. The two models differ solely in the PBL parameterization: one uses a prescribed K-profile PBL (KPP) scheme with an entrainment parameterization, and the other employs a turbulence kinetic energy (TKE) scheme. The models are evaluated through the comparison to the reanalysis ensemble, which is generated from ERA-20C, ERA-Interim, NCEP-CFSR and NASA-MERRA, and the following systematic biases are identified. The models exhibit wide-spread cold biases in the high latitudes, and the biases are smaller when the KPP scheme is used. The diurnal cycle amplitudes are underestimated in most dry regions, and the model with the TKE scheme simulates larger amplitudes. For the near-surface winds, the models underestimate both the daily means and the diurnal amplitudes. The differences between the models are relatively small compared to the biases.
The role of the PBL schemes in simulating the PBL parameters is investigated through the analysis of vertical profiles. The Sahara, which is suitable for focusing on the role of vertical mixing in dry PBLs, is selected for a detailed analysis. It reveals that compared to the KPP scheme, the heat transport is weaker with the TKE scheme in both convective and stable PBLs due to weaker vertical mixing, resulting in larger diurnal amplitudes. Lack of non-local momentum transport from the nocturnal low-level jets to the surfaces appears to explain the underestimation of the near-surface winds in the models.
Smyth, Jane E., Spencer A Hill, and Yi Ming, November 2018: Simulated Responses of the West African Monsoon and Zonal‐Mean Tropical Precipitation to Early Holocene Orbital Forcing. Geophysical Research Letters, 45(21), DOI:10.1029/2018GL080494. Abstract
This study seeks to improve our mechanistic understanding of how the insolation changes associated with orbital forcing impact the West African monsoon and zonal‐mean tropical precipitation. We impose early Holocene orbital parameters in simulations with the Geophysical Fluid Dynamics Laboratory AM2.1 atmospheric general circulation model, either with fixed sea surface temperatures, a 50‐m thermodynamic slab ocean, or coupled to a dynamic ocean (CM2.1). In all cases, West African Monsoon rainfall expands northward, but the summer zonal‐mean Intertropical Convergence Zone does not—there is drying near 10°N, and in the slab ocean experiment a southward shift of rainfall. This contradicts expectations from the conventional energetic framework for the Intertropical Convergence Zone location, given anomalous southward energy fluxes in the deep tropics. These anomalous energy fluxes are not accomplished by a stronger Hadley circulation; instead, they arise from an increase in total gross moist stability in the northern tropics.
Undorf, S, D Polson, Massimo Bollasina, and Yi Ming, et al., May 2018: Detectable impact of local and remote anthropogenic aerosols on the 20th century changes of West African and South Asian monsoon precipitation. Journal of Geophysical Research: Atmospheres, 123(10), DOI:10.1029/2017JD027711. Abstract
Anthropogenic aerosols are a key driver of changes in summer monsoon precipitation in the Northern Hemisphere during the 20th century. Here we apply detection and attribution methods to investigate causes of change in the West African and South Asian monsoons separately and identify the aerosol source regions that are most important for explaining the observed changes during 1920‐2005. Historical simulations with the GFDL‐CM3 model are used to derive fingerprints of aerosol forcing from different regions. For West Africa, remote aerosol emissions from North America and Europe (NAEU) are essential in order to detect the anthropogenic signal in observed monsoon precipitation changes. The changes are significantly underestimated in the model, however. While natural (volcanic) forcing seems to also play a role, the dominant contribution is found to come from aerosol‐induced changes in the inter‐hemispheric temperature gradient and associated meridional shifts of the Inter‐tropical Convergence Zone. For South Asia, in contrast, changes in observed monsoon precipitation can not be explained without local emissions. Here the findings show a weakening of the monsoon circulation, driven by the increase of remote NAEU aerosol emissions until 1975, and since then by the increase in local emissions offsetting the decrease of NAEU emissions. The results show that the aerosol forcing from individual emission regions is strong enough to be detected over internal variability. They also underscore the importance of the spatial pattern of global aerosol emissions, which is likely to continue to change throughout the expected near‐future decline in global emissions.
Xiang, Baoqiang, Ming Zhao, and Yi Ming, et al., July 2018: Contrasting Impacts of radiative forcing in the Southern Ocean versus Southern Tropics on ITCZ position and energy transport in one GFDL climate model. Journal of Climate, 31(14), DOI:10.1175/JCLI-D-17-0566.1. Abstract
Most current climate models suffer from pronounced cloud and radiation biases in the Southern Ocean (SO) and in the tropics. Using one GFDL climate model, this study investigates the migration of the Inter-tropical Convergence Zone (ITCZ) with prescribed Top of Atmosphere (TOA) shortwave radiative heating in the SO (50°S-80°S) versus the Southern Tropics (ST, 0-20°S). Results demonstrate that the ITCZ position response to the ST forcing is twice as strong as the SO forcing, which is primarily driven by the contrasting sea surface temperature (SST) gradient over the tropics; however, the mechanism for the formation of the SST pattern remains elusive.
Energy budget analysis reveals that the conventional energetic constraint framework is inadequate in explaining the ITCZ shift in these two perturbed experiments. For both cases, the anomalous Hadley circulation does not contribute to transport the imposed energy from the Southern Hemisphere to the Northern Hemisphere, given a positive mean gross moist stability in the equatorial region. Changes in the cross-equatorial atmospheric energy are primarily transported by atmospheric transient eddies when the anomalous ITCZ shift is most pronounced during December-May.
The partitioning of energy transport between the atmosphere and ocean shows latitudinal dependence: the atmosphere and ocean play an overall equivalent role in transporting the imposed energy for the extratropical SO forcing, while for the ST forcing, the imposed energy is nearly completely transported by the atmosphere. This contrast originates from the different ocean heat uptake and also the different meridional scale of the anomalous ocean circulation.
In this two-part paper, a description is provided of a version of the AM4.0/LM4.0 atmosphere/land model that will serve as a base for a new set of climate and Earth system models (CM4 and ESM4) under development at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). This version, with roughly 100km horizontal resolution and 33 levels in the vertical, contains an aerosol model that generates aerosol fields from emissions and a “light” chemistry mechanism designed to support the aerosol model but with prescribed ozone. In Part I, the quality of the simulation in AMIP (Atmospheric Model Intercomparison Project) mode – with prescribed sea surface temperatures (SSTs) and sea ice distribution – is described and compared with previous GFDL models and with the CMIP5 archive of AMIP simulations. The model's Cess sensitivity (response in the top-of-atmosphere radiative flux to uniform warming of SSTs) and effective radiative forcing are also presented. In Part II, the model formulation is described more fully and key sensitivities to aspects of the model formulation are discussed, along with the approach to model tuning.
In Part II of this two-part paper, documentation is provided of key aspects of a version of the AM4.0/LM4.0 atmosphere/land model that will serve as a base for a new set of climate and Earth system models (CM4 and ESM4) under development at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). The quality of the simulation in AMIP (Atmospheric Model Intercomparison Project) mode has been provided in Part I. Part II provides documentation of key components and some sensitivities to choices of model formulation and values of parameters, highlighting the convection parameterization and orographic gravity wave drag. The approach taken to tune the model's clouds to observations is a particular focal point. Care is taken to describe the extent to which aerosol effective forcing and Cess sensitivity have been tuned through the model development process, both of which are relevant to the ability of the model to simulate the evolution of temperatures over the last century when coupled to an ocean model.
Brown, Patrick T., Yi Ming, Wenhong Li, and Spencer A Hill, October 2017: Change in the magnitude and mechanisms of global temperature variability with warming. Nature Climate Change, 7(10), DOI:10.1038/nclimate3381. Abstract
Natural unforced variability in global mean surface air temperature (GMST) can mask or exaggerate human-caused global warming, and thus a complete understanding of this variability is highly desirable. Significant progress has been made in elucidating the magnitude and physical origins of present-day unforced GMST variability, but it has remained unclear how such variability may change as the climate warms. Here we present modelling evidence that indicates that the magnitude of low-frequency GMST variability is likely to decline in a warmer climate and that its generating mechanisms may be fundamentally altered. In particular, a warmer climate results in lower albedo at high latitudes, which yields a weaker albedo feedback on unforced GMST variability. These results imply that unforced GMST variability is dependent on the background climatological conditions, and thus climate model control simulations run under perpetual pre-industrial conditions may have only limited relevance for understanding the unforced GMST variability of the future.
Climate models generate a wide range of precipitation responses to global warming in the African Sahel, but all that use the NOAA Geophysical Fluid Dynamics Laboratory AM2.1 model as their atmospheric component dry the region sharply. This study compares the Sahel’s wet season response to uniform 2 K SST warming in AM2.1 using either its default convective parameterization, Relaxed Arakawa-Schubert (RAS), or an alternate, the University of Washington (UW) parameterization, using the moist static energy (MSE) budget to diagnose the relevant mechanisms.
UW generates a drier, cooler control Sahel climate than does RAS and a modest rainfall increase with SST warming rather than a sharp decrease. Horizontal advection of dry, low-MSE air from the Sahara Desert – a leading-order term in the control MSE budget with either parameterization – is enhanced with oceanic warming, driven by enhanced meridional MSE and moisture gradients spanning the Sahel. With RAS, this occurs throughout the free troposphere and is balanced by anomalous MSE import through anomalous subsidence, which must be especially large in the mid-troposphere where the moist static stability is small. With UW, the strengthening of the meridional MSE gradient is mostly confined to the lower troposphere, due in part to comparatively shallow prevailing convection. This necessitates less subsidence, enabling convective and total precipitation to increase with UW, although both large-scale precipitation and precipitation minus evaporation decrease. This broad set of hydrological and energetic responses persists in simulations with SSTs varied over a wide range.
Lee, J E., A Shen, Baylor Fox-Kemper, and Yi Ming, January 2017: Hemispheric sea ice distribution sets the glacial tempo. Geophysical Research Letters, 44(2), DOI:10.1002/2016GL071307. Abstract
The proxy record of global temperature shows that the dominant periodicity of the glacial cycle shifts from 40 kyr (obliquity) to 100 kyr (eccentricity) about a million years ago. Using climate model simulations, here we show that the pace of the glacial cycle depends on the pattern of hemispheric sea ice growth. In a cold climate the sea ice grows asymmetrically between two hemispheres under changes to Earth's orbital precession, because sea ice growth potential outside of the Arctic Circle is limited. This difference in hemispheric sea ice growth leads to an asymmetry in absorbed solar energy for the two hemispheres, particularly when eccentricity is high, even if the annual average insolation is similar. In a warmer climate, the hemispheric asymmetry of the sea ice decreases as mean Arctic and Antarctic sea ice decreases, diminishing the precession and eccentricity signals and explaining the dominant obliquity signal (40 kyr) before the mid-Pleistocene transition.
This paper investigates changes in the tropical tropopause layer (TTL) in response to carbon dioxide increase and surface warming separately in an atmospheric general circulation model, finding that both effects lead to a warmer tropical tropopause. Surface warming also results in an upward shift of the tropopause. A detailed heat budget analysis is performed to quantify the contributions from different radiative and dynamic processes to changes in the TTL temperature. When carbon dioxide increases with fixed surface temperature, a warmer TTL mainly results from the direct radiative effect of carbon dioxide increase. With surface warming, the largest contribution to the TTL warming comes from the radiative effect of the warmer troposphere, which is partly canceled by the radiative effect of the moistening at the TTL. Strengthening of the stratospheric circulation following surface warming cools the lower stratosphere dynamically and radiatively via changes in ozone. These two effects are of comparable magnitudes. This circulation change is the main cause of temperature changes near 63 hPa but is weak near 100 hPa. Contributions from changes in convection and clouds are also quantified. These results illustrate the heat budget analysis as a useful tool to disentangle the radiative–dynamical–chemical–convective coupling at the TTL and to facilitate an understanding of intermodel difference.
We contrast the responses to ozone depletion in two climate models: CAM3 and GFDL AM3. Although both models are forced with identical ozone concentration changes, the stratospheric cooling simulated in CAM3 is 30% stronger than in AM3 in annual mean, and twice as strong in December. We find that this difference originates from the dynamical response to ozone depletion, and its strength can be linked to the timing of the climatological springtime polar vortex breakdown. This mechanism is further supported by a variant of the AM3 simulation in which the Southern stratospheric zonal wind climatology is nudged to be CAM3-like. Given that the delayed breakdown of the Southern polar vortex is a common bias among many climate models, previous model-based assessments of the forced responses to ozone depletion may have been somewhat overestimated.
Pan, Fang, X Huang, S S Leroy, Pu Lin, L Larrabee Strow, Yi Ming, and V Ramaswamy, August 2017: The stratospheric changes inferred from 10 years of AIRS and AMSU-A radiances. Journal of Climate, 30(15), DOI:10.1175/JCLI-D-17-0037.1. Abstract
We analyze global-mean radiances observed by AIRS (Atmospheric Infrared Sounder) and AMSU-A (Advanced Microwave Sounding Unit) from 2003 to 2012. We focus on channels sensitive to emission and absorption in the stratosphere. Optimal fingerprinting is used to obtain estimates of changes of stratospheric temperature in five vertical layers due to external forcing in the presence of natural variability. Natural variability is estimated using synthetic radiances based on the 500-year GFDL CM3 and 240-year HADGEM2-CC control runs. The results show a cooling rate of 0.65±0.11(2σ) K decade-1 in the upper stratosphere above 6hPa, ~0.46±0.24 K decade-1 in two middle stratospheric layers between 6hPa and 30hPa, and 0.39±0.32 K decade-1 in the lower stratosphere (30-60hPa). The cooling rate in the lowest part of the stratosphere (60-100hPa) is -0.014±0.22 K decade-1, which is smallest among all five layers and statistically insignificant. The synergistic use of well-calibrated passive infrared and microwave radiances permits disambiguation of trends of carbon dioxide and stratospheric temperature, increases vertical resolution of detected stratospheric temperature trends, and effectively reduces uncertainties of estimated temperature trends.
East Asia has some of the largest concentrations of absorbing aerosols globally, and these, along with the region’s scattering aerosols, have both reduced the amount of solar radiation reaching the Earth’s surface regionally (“solar dimming”) and increased shortwave absorption within the atmosphere, particularly during the peak months of the East Asian Summer Monsoon (EASM). This study analyzes how atmospheric absorption and surface solar dimming compete in driving the response of regional summertime climate to anthropogenic aerosols, which dominates, and why—issues of particular importance for predicting how East Asian climate will respond to projected changes in absorbing and scattering aerosol emissions in the future. These questions are probed in a state-of-the-art general circulation model using a combination of realistic and novel idealized aerosol perturbations that allow analysis of the relative influence of absorbing aerosols’ atmospheric and surface-driven impacts on regional circulation and climate. Results show that even purely absorption-driven dimming decreases EASM precipitation by cooling the land surface, counteracting climatological land-sea contrast and reducing ascending atmospheric motion and on-shore winds, despite the associated positive top-of-atmosphere regional radiative forcing. Absorption-driven atmospheric heating does partially offset the precipitation and surface evaporation reduction from surface dimming, but the overall response to aerosol absorption more closely resembles the response to its surface dimming than to its atmospheric heating. These findings provide a novel decomposition of absorbing aerosol’s impacts on regional climate and demonstrate that the response cannot be expected to follow the sign of absorption’s top-of-atmosphere or even atmospheric radiative perturbation.
Arctic haze has a distinct seasonal cycle with peak concentrations in winter but pristine conditions in summer. It is demonstrated that the Geophysical Fluid Dynamics Laboratory (GFDL) atmospheric general circulation model AM3 can reproduce the observed seasonality of Arctic black carbon (BC), an important component of Arctic haze. We use the model to study how large-scale circulation and removal drive the seasonal cycle of Arctic BC. It is found that despite large seasonal shifts in the general circulation pattern, the transport of BC into the Arctic varies little throughout the year. The seasonal cycle of Arctic BC is attributed mostly to variations in the controlling factors of wet removal, namely the hydrophilic fraction of BC and wet deposition efficiency of hydrophilic BC. Specifically, a confluence of low hydrophilic fraction and weak wet deposition, owing to slower aging process and less efficient mixed-phase cloud scavenging, respectively, is responsible for the wintertime peak of BC. The transition to low BC in summer is the consequence of a gradual increase in the wet deposition efficiency, while the increase of BC in late fall can be explained by a sharp decrease in the hydrophilic fraction. The results presented here suggest that future changes in the aging and wet deposition processes can potentially alter the concentrations of Arctic aerosols and their climate effects.
Wang, Yan, Yuanyu Xie, Wenhao Dong, and Yi Ming, et al., October 2017: Adverse Effects of Increasing Drought on Air Quality via Natural Processes. Atmospheric Chemistry and Physics, 17(20), DOI:10.5194/acp-17-12827-2017. Abstract
Drought is a recurring extreme of the climate system with well-documented impacts on agriculture and water resources. The strong perturbation of drought to the land biosphere and atmospheric water cycle will affect atmospheric composition, the nature and extent of which are not well understood. Here we present observational evidence that surface ozone and PM2.5 in the US are significantly correlated with drought severity, with 3.5 ppbv (8 %) and 1.6 μg m−3 (17 %) increases respectively under severe drought. These enhancements show little sensitivity to the decreasing trend of US anthropogenic emissions, indicating natural processes as the primary cause. Elevated ozone and PM2.5 are attributed to the combined effects of drought on natural emissions (wildfires, biogenic VOCs and dust), deposition, and chemistry. Most climate-chemistry models are not able to reproduce the observed responses of ozone and PM2.5 to drought severity, suggesting a lack of mechanistic understanding of drought effects on atmospheric composition. The model deficiencies are partly attributed to the lack of drought-induced changes in land-atmosphere exchanges of reactive gases and particles and misrepresentation of cloud changes under drought conditions. By applying the observed relationships between drought and air pollutants to climate model projected drought occurrences, we estimate an increase of 1–6 % for ground-level O3 and 1–16 % for PM2.5 in the US by 2100 compared to the 2000s due to increasing drought alone. Drought thus poses another aspect of climate change penalty on air quality not recognized before. Improvements in the models are imperative to facilitate better prediction of air quality challenges due to changing hydroclimate and atmospheric composition feedback to climate.
Dong, Wenhao, Yanluan Lin, Jonathon S Wright, and Yi Ming, et al., March 2016: Summer rainfall over the southwestern Tibetan Plateau controlled by deep convection over the Indian subcontinent. Nature Communications, 7, 10925, DOI:10.1038/ncomms10925. Abstract
Despite the importance of precipitation and moisture transport over the Tibetan Plateau for
glacier mass balance, river runoff and local ecology, changes in these quantities remain highly
uncertain and poorly understood. Here we use observational data and model simulations to
explore the close relationship between summer rainfall variability over the southwestern
Tibetan Plateau (SWTP) and that over central-eastern India (CEI), which exists despite the
separation of these two regions by the Himalayas. We show that this relationship is maintained
primarily by ‘up-and-over’ moisture transport, in which hydrometeors and moisture are
lifted by convective storms over CEI and the Himalayan foothills and then swept over the
SWTP by the mid-tropospheric circulation, rather than by upslope flow over the Himalayas.
Sensitivity simulations confirm the importance of up-and-over transport at event scales, and
an objective storm classification indicates that this pathway accounts for approximately half
of total summer rainfall over the SWTP.
Li, Z, K M Lau, V Ramanathan, Guoxiong Wu, Y Ding, M G Manoj, J Liu, Yitian Qian, J Li, Tianjun Zhou, J Fan, Daniel Rosenfeld, and Yi Ming, et al., December 2016: Aerosol and Monsoon Climate Interactions over Asia. Reviews of Geophysics, 54(4), DOI:10.1002/2015RG000500. Abstract
The increasing severity of droughts/floods and worsening air quality from increasing aerosols in Asia monsoon regions are the two gravest threats facing over 60% of the world population living in Asian monsoon regions. These dual threats have fueled a large body of research in the last decade on the roles of aerosols in impacting Asian monsoon weather and climate. This paper provides a comprehensive review of studies on Asian aerosols, monsoons, and their interactions. The Asian monsoon region is a primary source of emissions of diverse species of aerosols from both anthropogenic and natural origins. The distributions of aerosol loading are strongly influenced by distinct weather and climatic regimes, which are, in turn, modulated by aerosol effects. On a continental scale, aerosols reduce surface insolation and weaken the land-ocean thermal contrast, thus inhibiting the development of monsoons. Locally, aerosol radiative effects alter the thermodynamic stability and convective potential of the lower atmosphere leading to reduced temperatures, increased atmospheric stability, and weakened wind and atmospheric circulations. The atmospheric thermodynamic state, which determines the formation of clouds, convection, and precipitation, may also be altered by aerosols serving as cloud condensation nuclei or ice nuclei. Absorbing aerosols such as black carbon and desert dust in Asian monsoon regions may also induce dynamical feedback processes, leading to a strengthening of the early monsoon and affecting the subsequent evolution of the monsoon. Many mechanisms have been put forth regarding how aerosols modulate the amplitude, frequency, intensity, and phase of different monsoon climate variables. A wide range of theoretical, observational, and modeling findings on the Asian monsoon, aerosols, and their interactions are synthesized. A new paradigm is proposed on investigating aerosol-monsoon interactions, in which natural aerosols such as desert dust, black carbon from biomass burning, and biogenic aerosols from vegetation are considered integral components of an intrinsic aerosol-monsoon climate system, subject to external forcing of global warming, anthropogenic aerosols, and land use and change. Future research on aerosol-monsoon interactions calls for an integrated approach and international collaborations based on long-term sustained observations, process measurements, and improved models, as well as using observations to constrain model simulations and projections.
Uncertainty in equilibrium climate sensitivity impedes accurate climate projections. While the inter-model spread is known to arise primarily from differences in cloud feedback, the exact processes responsible for the spread remain unclear. To help identify some key sources of uncertainty, we use a developmental version of the next generation Geophysical Fluid Dynamics Laboratory global climate model (GCM) to construct a tightly controlled set of GCMs where only the formulation of convective precipitation is changed. The different models provide simulation of present-day climatology of comparable quality compared to the CMIP5 model ensemble. We demonstrate that model estimates of climate sensitivity can be strongly affected by the manner through which cumulus cloud condensate is converted into precipitation in a model’s convection parameterization, processes that are only crudely accounted for in GCMs. In particular, two commonly used methods for converting cumulus condensate into precipitation can lead to drastically different climate sensitivity, as estimated here with an atmosphere/land model by increasing sea surface temperatures uniformly and examining the response in the top-of-atmosphere energy balance. The effect can be quantified through a bulk convective detrainment efficiency, which measures the ability of cumulus convection to generate condensate per unit precipitation. The model differences, dominated by shortwave feedbacks, come from broad regimes ranging from large-scale ascent to subsidence regions. Given current uncertainties in representing convective precipitation microphysics and our current inability to find a clear observational constraint that favors one version of our model over the others, the implications of this ability to engineer climate sensitivity needs to be considered when estimating the uncertainty in climate projections.
Anthropogenically forced changes to the mean and spatial pattern of sea surface temperatures (SSTs) alter tropical atmospheric meridional energy transport throughout the seasonal cycle – in total, its partitioning between the Hadley cells and eddies, and, for the Hadley cells, the relative roles of the mass flux and the gross moist stability (GMS). We investigate this behavior using an atmospheric general circulation model forced with SST anomalies caused by either historical greenhouse gas or aerosol forcing, dividing the SST anomalies into two components: the tropical mean SST anomaly applied uniformly, and the full SST anomalies minus the tropical mean.
For greenhouse gases, the polar-amplified SST spatial pattern partially negates enhanced eddy poleward energy transport driven by mean warming. Both SST components weaken winter Hadley cell circulation and alter GMS. The Northern Hemisphere-focused aerosol cooling induces northward energy flux anomalies in the deep tropics, which manifest partially via strengthened northern and weakened southern Hadley cell overturning. Aerosol-induced GMS changes also contribute to the northward energy fluxes. A simple thermodynamic scaling qualitatively captures these changes, though it performs less well for the greenhouse gas simulations. The scaling provides an explanation for the tight correlation demonstrated in previous studies between shifts in the Intertropical Convergence Zone position and cross-equatorial energy fluxes.
The behavior of the Brewer-Dobson circulation is investigated using a suite of global climate model simulations with different forcing agents, in conjunction with observation-based analysis. We find that the variations in the Brewer-Dobson circulationare strongly correlated with those in the tropical-mean surface temperature through changes in the upper tropospheric temperature and zonal winds. This correlation is seen on both interannual and multi-decadal timescales, and holds for natural and forced variations alike. The circulation change is relatively insensitive to the spatial pattern of the forcings. Consistent changes in the Brewer-Dobson circulation with respect to those in the tropical-mean surface temperature prevail across timescales and forcings, and constitute an important attribution element of the atmospheric adjustment to global climate change.
Ban-Weiss, G A., Ling Jin, Susanne E Bauer, R Bennartz, Xiaoping Liu, Kai Zhang, Yi Ming, Huan Guo, and J H Jiang, September 2014: Evaluating clouds, aerosols, and their interactions in three global climate models using satellite simulators and observations. Journal of Geophysical Research: Atmospheres, 119(18), DOI:10.1002/2014JD021722. Abstract
Accurately representing aerosol-cloud interactions in global climate models is challenging. As parameterizations evolve, it is important to evaluate their performance with appropriate use of observations. In this investigation we compare aerosols, clouds, and their interactions in three global climate models (GFDL-AM3, NCAR-CAM5, GISS-ModelE2) to MODIS satellite observations. Modeled cloud properties are diagnosed using a MODIS simulator. Cloud droplet number concentrations (N) are computed identically from satellite-simulated and MODIS-observed values of liquid cloud optical depth and droplet effective radius. We find that aerosol optical depth (τa) simulated by models is similar to observations in many regions around the globe. For N, AM3 and CAM5 capture the observed spatial pattern of higher values in coastal marine stratocumulus versus remote ocean regions, though modeled values in general are higher than observed. Aerosol-cloud interactions were computed as the sensitivity of ln(N) to ln(τa) for coastal marine liquid clouds near South Africa (SAF) and Southeast Asia (SEA) where τa varies in time. AM3 and CAM5 are more sensitive than observations, while the sensitivity for ModelE2 is statistically insignificant. This widely used sensitivity could be subject to misinterpretation due to the confounding influence of meteorology on both aerosols and clouds. A simple framework for assessing the sensitivity of ln(N) to ln(τa) at constant meteorology illustrates that observed sensitivity can change from positive to statistically insignificant when including the confounding influence of relative humidity. Satellite-simulated versus standard model values of N from CAM5 are compared in SAF; standard model values are significantly lower with a bias of 83 cm−3.
The late 20th century response of the South Asian monsoon to changes in anthropogenic aerosols from local (i.e., South Asia) and remote (i.e., outside South Asia) sources was investigated using historical simulations with a state-of-the-art climate model. The observed summertime drying over India is replaced by widespread wettening once local aerosol emissions are kept at pre-industrial levels while all the other forcings evolve. Constant remote aerosol emissions partially suppress the precipitation decrease. While predominant precipitation changes over India are thus associated with local aerosols, remote aerosols contribute as well, especially in favoring an earlier monsoon onset in June and enhancing summertime rainfall over the northwestern regions. Conversely, temperature and near-surface circulation changes over South Asia are more effectively driven by remote aerosols. These changes are reflected into northward cross-equatorial anomalies in the atmospheric energy transport induced by both local and, to a greater extent, remote aerosols.
Brown, Patrick T., Wenhong Li, L Li, and Yi Ming, July 2014: Top-of-atmosphere radiative contribution to unforced decadal global temperature variability in climate models. Geophysical Research Letters, 41(14), DOI:10.1002/2014GL060625. Abstract
Much recent work has focused on unforced global mean surface air temperature (T) variability associated with the efficiency of heat transport into the deep ocean. Here the relationship between unforced variability in T and the Earth's top-of-atmosphere (TOA) energy balance is explored in preindustrial control runs of the Coupled Model Intercomparison Project Phase 5 multimodel ensemble. It is found that large decadal scale variations in T tend to be significantly enhanced by the net energy flux at the TOA. This indicates that unforced decadal variability in T is not only caused by a redistribution of heat within the climate system but can also be associated with unforced changes in the total amount of heat in the climate system. It is found that the net TOA radiation imbalances result mostly from changes in albedo associated with the Interdecadal Pacific Oscillation that temporarily counteracts the climate system's outgoing longwave (i.e., Stefan-Boltzmann) response to T change.
Ocko, I B., V Ramaswamy, and Yi Ming, July 2014: Contrasting Climate Responses to the Scattering and Absorbing Features of Anthropogenic Aerosol Forcings. Journal of Climate, 27(14), DOI:10.1175/JCLI-D-13-00401.1. Abstract
Anthropogenic aerosols are comprised of optically scattering and absorbing particles, with the principal concentrations being in the Northern Hemisphere, and yielding negative and positive global mean radiative forcings, respectively. Aerosols also influence cloud albedo, yielding additional negative radiative forcings. Climate responses to a comprehensive set of isolated aerosol forcing simulations are investigated in a coupled atmosphere-ocean framework, forced by preindustrial to present-day aerosol-induced radiative perturbations. Atmospheric and oceanic climate responses (including precipitation, atmospheric circulation, atmospheric and oceanic heat transport, sea surface temperature, and salinity) to positive and negative particulate forcings are consistently anti-correlated. The striking effects include distinct patterns of changes north and south of the equator that are governed by the sign of the aerosol forcing and its initiation of an interhemispheric forcing asymmetry. The presence of opposing signs of the forcings between the aerosol scatterers and absorbers, and the resulting contrast in climate responses, thus dilutes the effects of individual aerosol types on influencing global and regional climate conditions. The aerosol-induced changes in the variables also have a distinct fingerprint when compared to the responses of the more globally uniform and interhemispherically symmetric well-mixed greenhouse gas forcing. The significance of employing a full ocean model is demonstrated in this study by the ability to partition how individual aerosols influence atmospheric and oceanic conditions separately.
Persad, Geeta, Yi Ming, and V Ramaswamy, September 2014: The Role of Aerosol Absorption in Driving Clear-Sky Solar Dimming over East Asia. Journal of Geophysical Research: Atmospheres, 119(17), DOI:10.1002/2014JD021577. Abstract
Surface-based observations indicate a significant decreasing trend in clear-sky downward surface solar radiation (SSR) over East Asia since the 1960s. This “dimming" is thought to be driven by the region's long-term increase in aerosol emissions, butlittle work has been done to quantify the underlying physical mechanisms or the contribution from aerosol absorption within the atmospheric column. Given the distinct climate impacts that absorption-driven dimming may produce, this constitutes an important, but thus far rather neglected, line of inquiry.
We examine experiments conducted in the Geophysical Fluid Dynamics Laboratory's Atmospheric General Circulation Models, AM2.1 and AM3, in order to analyze the model-simulated East Asian clear-sky SSR trends. We also use the models’ standalone radiation module to examine the contribution from various aerosol characteristics in the two models (such as burden, mixing state, hygroscopicity, and seasonal distribution) to the trends. Both models produce trends in clear-sky SSR that are comparable tothat observed, but via disparate mechanisms. Despite their different aerosol characteristics, the models produce nearly identical increases in aerosol absorption since the 1960s, constituting as much as half of the modeled clear-sky dimming. This is due to a compensation between the differences in aerosol column burden and mixing state assumed in the two models, i.e. plausible clear-sky SSR simulations can be achieved via drastically different aerosol parameterizations. Our novel results indicatethat trends in aerosol absorption drive a large portion of East Asian clear-sky solar dimming in the models presented here and for the time periods analyzed, and that mechanistic analysis of the factors involved in aerosol absorption is an important diagnostic in evaluating modeled clear-sky solar dimming trends.
Rotstayn, L D., E L Plymin, Mark A Collier, Olivier Boucher, J-L Dufresne, J-J Luo, K von Salzen, S J Jeffrey, M-A Foujols, Yi Ming, and Larry W Horowitz, September 2014: Declining aerosols in CMIP5 projections: effects on atmospheric temperature structure and midlatitude jets. Journal of Climate, 27(18), DOI:10.1175/JCLI-D-14-00258.1. Abstract
We assess the effects of declining anthropogenic aerosols in Representative Concentration Pathway 4.5 (RCP4.5) in four models from the Coupled Model Intercomparison Project Phase 5 (CMIP5), with a focus on annual, zonal-mean atmospheric temperature structure and zonal winds. For each model, we diagnose the effect of declining aerosols from the difference between a projection forced by RCP4.5 for 2006–2100 and another that has identical forcing, except that anthropogenic aerosols are fixed at early 21st century levels. The response to declining aerosols is interpreted in terms of the meridional structure of aerosol radiative forcing, which peaks near 40°N and vanishes at the South Pole.
Increasing greenhouse gases cause amplified warming in the tropical upper troposphere and strengthening midlatitude jets in both hemispheres. However, for declining aerosols the vertically averaged tropospheric temperature response peaks near 40°N, rather than in the tropics. This implies that for declining aerosols the tropospheric meridional temperature gradient generally increases in the Southern Hemisphere (SH), but in the Northern Hemisphere (NH) it decreases in the tropics and subtropics. Consistent with thermal wind balance, the NH jet then strengthens on its poleward side and weakens on its equatorward side, whereas the SH jet strengthens more than the NH jet. The asymmetric response of the jets is thus consistent with the meridional structure of aerosol radiative forcing and the associated tropospheric warming: in the NH the latitude of maximum warming is roughly collocated with the jet, whereas in the SH warming is strongest in the tropics and weakest at high latitudes.
Bollasina, Massimo, and Yi Ming, February 2013: The general circulation model precipitation bias over the southwestern equatorial Indian Ocean and its implications for simulating the South Asian monsoon. Climate Dynamics, 40(3-4), DOI:10.1007/s00382-012-1347-7. Abstract
Most of current general circulation models
(GCMs) show a remarkable positive precipitation bias over
the southwestern equatorial Indian Ocean (SWEIO), which
can be thought of as a westward expansion of the simulated
IO convergence zone toward the coast of Africa. The bias
is common to both coupled and uncoupled models, suggesting
that its origin does not stem from the way boundary
conditions are specified. The spatio-temporal evolution of
the precipitation and associated three-dimensional atmospheric
circulation biases is comprehensively characterized
by comparing the GFDL AM3 atmospheric model to
observations. It is shown that the oceanic bias, which
develops in spring and reduces during the monsoon season,
is associated to a consistent precipitation and circulation
anomalous pattern over the whole Indian region. In the
vertical, the areas are linked by an anomalous Hadley-type
meridional circulation, whose northern branch subsides
over northeastern India significantly affecting the monsoon
evolution (e.g., delaying its onset). This study makes the
case that the precipitation bias over the SWEIO is forced
by the model excess response to the local meridional sea
surface temperature (SST) gradient through enhanced nearsurface
meridional wind convergence. This is suggested by
observational evidence and supported by AM3 sensitivity
experiments. The latter show that relaxing the magnitude
of the meridional SST gradient in the SWEIO can lead to a
significant reduction of both local and large-scale
precipitation and circulation biases. The ability of local
anomalies over the SWEIO to force a large-scale remote
response to the north is further supported by numerical
experiments with the GFDL spectral dry dynamical core
model. By imposing a realistic anomalous heating source
over the SWEIO the model is able to reproduce the main
dynamical features of the AM3 bias. These results indicate
that improved GCM simulations of the South Asian summer
monsoon could be achieved by reducing the springtime
model bias over the SWEIO. Deficiencies in the
atmospheric model, and in particular in the convective
parameterization, are suggested to play a key role. Finally,
the important mechanism controlling the simulated precipitation
distribution over South Asia found here should
be considered in the interpretation and attribution
Bollasina, Massimo, and Yi Ming, November 2013: The role of land-surface processes in modulating the Indian monsoon annual cycle. Climate Dynamics, 41(9-10), DOI:10.1007/s00382-012-1634-3. Abstract
The annual cycle of solar radiation, together with the resulting land–ocean differential heating, is traditionally considered the dominant forcing controlling the northward progression of the Indian monsoon. This study makes use of a state-of-the-art atmospheric general circulation model in a realistic configuration to conduct “perpetual†experiments aimed at providing new insights into the role of land–atmosphere processes in modulating the annual cycle of precipitation over India. The simulations are carried out at three important stages of the monsoon cycle: March, May, and July. Insolation and SSTs are held fixed at their respective monthly mean values, thus eliminating any external seasonal forcing. In the perpetual May experiment both precipitation and circulation are able to considerably evolve only by regional internal land–atmosphere processes and the mediation of soil hydrology. A large-scale equilibrium state is reached after approximately 270 days, closely resembling mid-summer climatological conditions. As a result, despite the absence of external forcing, intense and widespread rains over India are able to develop in the May-like state. The interaction between soil moisture and circulation, modulated by surface heating over the northwestern semi-arid areas, determines a slow northwestward migration of the monsoon, a crucial feature for the existence of desert regions to the west. This also implies that the land–atmosphere system in May is far from being in equilibrium with the external forcing. The inland migration of the precipitation front comprises a succession of large-scale 35–50 day coupled oscillations between soil moisture, precipitation, and circulation. The oscillatory regime is self-sustained and entirely due to the internal dynamics of the system. In contrast to the May case, minor changes in the land–atmosphere system are found when the model is initialized in March and, more surprisingly, in July, the latter case further emphasizing the role of northwestern surface heating.
Bollasina, Massimo, Yi Ming, and V Ramaswamy, July 2013: Earlier onset of the Indian Monsoon in the late 20th century: The role of anthropogenic aerosols. Geophysical Research Letters, 40(14), DOI:10.1002/grl.50719. Abstract
The impact of the late 20th century increase of anthropogenic aerosols on the Indian monsoon onset was investigated with a state-of-the-art climate model with fully-interactive aerosols and chemistry. We find that aerosols are likely responsible for the observed earlier onset, resulting in enhanced June precipitation over most of India. This shift is preceded by strong aerosol forcing over the Bay of Bengal and Indochina, mostly attributable to the direct effect, resulting in increased atmospheric stability that inhibits the monsoon migration in May. The adjusted atmospheric circulation leads to thermodynamical changes over the northwestern continental region, including increased surface temperature and near-surface moist static energy, which support a stronger June flow and, facilitated by a relative warming of the Indian Ocean, a vigorous northwestward precipitation shift. These findings underscore the importance of dynamical feedbacks and of regional land-surface processes for the aerosol-monsoon link.
Huang, X, H-W Chuang, A Dessler, X Chen, K Minschwaner, Yi Ming, and V Ramaswamy, March 2013: A radiative-convective equilibrium perspective of the weakening of tropical Walker circulation in response to global warming. Journal of Climate, 26(5), DOI:10.1175/JCLI-D-12-00288.1. Abstract
Both observational analysis and GCM model simulations indicate that the tropical
Walker circulation is becoming weaker and may continue to weaken as a consequence of
climate change. Here we use a conceptual radiative-convective-equilibrium (RCE)
framework to interpret the weakening of the Walker circulation as simulated by the
GFDL coupled-GCM. Based on the modeled lapse rate and clear-sky cooling rate profiles,
the RCE framework can directly compute the change of vertical velocity in the
descending branch of the Walker circulation, which agrees with the counterpart simulated
by the GFDL model. Our results show that the vertical structure of clear-sky radiative
cooling rate (QR) will change in response to the increased water vapor as the globe warms.
We explain why the change of QR is positive in the upper most part of the troposphere
(<300 hPa) and is negative for the rest part of the troposphere. As a result, both the
change of clear-sky cooling rate and the change of tropospheric lapse rate contribute to
the weakening of circulation. The vertical velocity changes due to the two factors are
comparable to each other from the top of planetary boundary layer to 600hPa. From
600hPa to 300hPa, lapse rate changes are the dominant cause of the weakening
circulation. Above 300hPa, the change due to QR is opposite to the change due to lapse
rate, which forces a slight increase in vertical velocity that is seen in the model simulation.
Employing the Geophysical Fluid Dynamics Laboratory (GFDL)'s fully-coupled chemistry-climate (ocean/atmosphere/land/sea ice) model (CM3) with an explicit physical representation of aerosol indirect effects (cloud-water droplet activation), we find that the dramatic emission reductions (35–80%) in anthropogenic aerosols and their precursors projected by Representative Concentration Pathway (RCP) 4.5 result in ~1°C of additional warming and ~0.1 mm day−1 of additional precipitation, both globally averaged, by the end of the 21st century. The impact of these reductions in aerosol emissions on simulated global mean surface temperature and precipitation becomes apparent by mid-21st century. Furthermore, we find that the aerosol emission reductions cause precipitation to increase in East and South Asia by ~1.0 mm day−1 through the 2nd half of the 21st century. Both the simulated temperature and precipitation responses in CM3 are significantly stronger than the previously simulated responses in our earlier climate model (CM2.1) that only considered direct radiative forcing by aerosols. We conclude that sulfate aerosol indirect effects greatly enhance the impacts of aerosols on surface temperature in CM3, while both direct and indirect effects from sulfate aerosols dominate the strong precipitation response, possibly with a small contribution from carbonaceous aerosols. Just as we found with the previous GFDL model, CM3 produces surface warming patterns that are uncorrelated with the spatial distribution of 21stcentury changes in aerosol loading. However, the largest precipitation increases in CM3 are co-located with the region of greatest aerosol decrease, in and downwind of Asia.
A set of GFDL AM2 sensitivity simulations by varying an entrainment threshold rate to control deep convection occurrence are used to investigate how cumulus parameterization impacts tropical cloud and precipitation characteristics. In the Tropics, model convective precipitation (CP) is frequent and light, while large-scale precipitation (LSP) is intermittent and strong. With deep convection inhibited, CP decreases significantly over land and LSP increases prominently over ocean. This results in an overall redistribution of precipitation from land to ocean. A composite analysis reveals that cloud fraction (low and middle) and cloud condensate associated with LSP is substantially larger than those associated with CP. With about the same total precipitation and precipitation frequency distribution over the Tropics, simulations having greater LSP fraction tend to have larger cloud condensate and low and middle cloud fraction.
Simulations having greater LSP fraction tend to be drier and colder in the upper-troposphere. The induced unstable stratification supports strong transient wind perturbations and LSP. Greater LSP also contributes to greater intraseasonal (20-100 day) precipitation variability. Model LSP has a close connection to the low level convergence via the resolved grid-scale dynamics and thus a close coupling with the surface heat flux. Such wind-evaporation feedback is essential to the development and maintenance of LSP and enhances model precipitation variability. LSP has stronger dependence and sensitivity on column moisture than CP. The moisture-convection feedback, critical to tropical intraseasonal variability, is enhanced in simulations with large LSP. Strong precipitation variability accompanied by the worse mean state implies that an optimal precipitation partitioning is critical to model tropical climate simulation.
Identifying the prime drivers of the twentieth-century multidecadal variability in the Atlantic Ocean is crucial for predicting how the Atlantic will evolve in the coming decades and the resulting broad impacts on weather and precipitation patterns around the globe. Recently Booth et al (2012) showed that the HadGEM2-ES climate model closely reproduces the observed multidecadal variations of area-averaged North Atlantic sea surface temperature in the 20th century. The multidecadal variations simulated in HadGEM2-ES are primarily driven by aerosol indirect effects that modify net surface shortwave radiation. On the basis of these results, Booth et al (2012) concluded that aerosols are a prime driver of twentieth-century North Atlantic climate variability. However, here it is shown that there are major discrepancies between the HadGEM2-ES simulations and observations in the North Atlantic upper ocean heat content, in the spatial pattern of multidecadal SST changes within and outside the North Atlantic, and in the subpolar North Atlantic sea surface salinity. These discrepancies may be strongly influenced by, and indeed in large part caused by, aerosol effects. It is also shown that the aerosol effects simulated in HadGEM2-ES cannot account for the observed anti-correlation between detrended multidecadal surface and subsurface temperature variations in the tropical North Atlantic. These discrepancies cast considerable doubt on the claim that aerosol forcing drives the bulk of this multidecadal variability.
Fan, Songmiao, J P Schwarz, Junfeng Liu, D W Fahey, Paul Ginoux, Larry W Horowitz, Hiram Levy II, Yi Ming, and J R Spackman, December 2012: Inferring ice formation processes from global-scale black carbon profiles observed in the remote atmosphere and model simulations. Journal of Geophysical Research: Atmospheres, 117, D23205, DOI:10.1029/2012JD018126. Abstract
Black carbon (BC) aerosol absorbs solar radiation and can act as cloud condensation nucleus and ice formation nucleus. The current generation of climate models have difficulty in accurately predicting global-scale BC concentrations. Previously, an ensemble of such models was compared to measurements, revealing model biases in the tropical troposphere and in the polar troposphere. Here, global aerosol distributions are simulated using different parameterizations of wet removal and model results are compared to BC profiles observed in the remote atmosphere to explore the possible sources of these biases. The model-data comparison suggests a slow removal of BC aerosol during transport to the Arctic in winter and spring, because ice crystal growth causes evaporation of liquid cloud via the Bergron process and, hence, release of BC aerosol back to ambient air. By contrast, more efficient model wet removal is needed in the cold upper troposphere over the tropical Pacific. Parcel model simulations with detailed droplet and ice nucleation and growth processes suggest that ice formation in this region may be suppressed due to a lack of ice nuclei (mainly insoluble dust particles) in the remote atmosphere, allowing liquid and mixed-phase clouds to persist under freezing temperatures, and forming liquid precipitation capable of removing aerosol incorporated in cloud water. Falling ice crystals can scavenge droplets in lower clouds, which also results in efficient removal of cloud condensation nuclei. The combination of models with global-scale BC measurements in this study has provided new, latitude-dependent information on ice formation processes in the atmosphere, and highlights the importance of a consistent treatment of aerosol and moist physics in climate models.
To counteract global warming, there have been suggestions to increase the albedo of low-level marine clouds through the aerosol indirect effects by injecting them with sea salt. However, the full climate response to this geoengineering scheme is currently poorly understood. We simulate cloud seeding in a coupled mixed-layer ocean-atmosphere general circulation model in order to identify the specific physical mechanisms through which seeding could perturb the climate system's radiative balance, and cause temperature and precipitation changes. Seeding stratocumulus decks over three tropical maritime regions in the North Pacific, South Pacific and South Atlantic produces strong local reductions in solar absorption. Over half of the radiative cooling is due to direct scattering of solar radiation by the added sea salt aerosols, while the rest comes from enhancement of the local cloud albedo. The oceanic cooling due to the seeding over the southeastern equatorial Pacific induces a La Ni\~na-like response, with tropical precipitation changes resembling La Ni\~na anomalies and teleconnections occurring in the mid-latitude North Pacific and North America. Additionally, model runs in which only one of the three regions is seeded indicate nonlinearity in the climate response. We identify dynamical and thermodynamical constraints respectively on the temperature and hydrological cycle responses to cloud seeding, but the full response to such geoengineering remains poorly constrained.
Radiative flux perturbation (RFP) is defined as the top-of-the-atmosphere (TOA) radiative imbalance after the atmosphere-land system adjusts fully to an external perturbation, and serves as a useful metric for quantifying climate forcing. This paper presents an effort to address the issue whether a forcing imposed initially over a specific region may alter the radiative balance elsewhere through atmospheric circulation, thus giving rise to a non-local component of RFP. We perturb an atmospheric general circulation model (AGCM) by increasing the cloud droplet number concentration exclusively over land, and observe widespread positive (warming) RFP over ocean, along with the expected negative (cooling) RFP over land. A detailed analysis suggests that the oceanic (or non-local) RFP is closely associated with a reduction in low cloud amount, which can be attributed primarily to the horizontal advection of drier land boundary layer air and to the oceanic boundary layer top entrainment of drier free troposphere air. By examining how the land surface and atmospheric energy balances are re-established, we are able to identify the physical mechanisms behind the strong hydrological impact of a tropical land shortwave (SW) forcing (the multiplier effect). In contrast, the hydrological cycle is relatively insensitive to an extratropical forcing.
The direct radiative forcing of the climate system includes effects due to scattering and absorbing aerosols. This study explores how important physical climate characteristics contribute to the magnitudes of the direct radiative forcings (DRF) from anthropogenic sulfate, black carbon, and organic carbon. For this purpose, we employ the GFDL CM2.1 global climate model, which has reasonable aerosol concentrations and reconstruction of 20th Century climate change. Sulfate and carbonaceous aerosols constitute the most important anthropogenic aerosol perturbations to the climate system, and provide striking contrasts between primarily scattering (sulfate and organic carbon) and primarily absorbing (black carbon) species. The quantitative roles of cloud coverage, surface albedo, and relative humidity in governing the sign and magnitude of all-sky top-of-atmosphere forcings are examined. Clouds reduce the global-mean sulfate TOA DRF by almost 50%, reduce the global-mean organic carbon TOA DRF by more than 30%, and increase the global-mean black carbon TOA DRF by almost 80%. Sulfate forcing is increased by over 50% as a result of hygroscopic growth, while high-albedo surfaces are found to have only a minor (less than 10%) impact on all global-mean forcings. Although the radiative forcing magnitudes are subject to uncertainties in the state of mixing of the aerosol species, it is clear that fundamental physical climate characteristics play a large role in governing aerosol direct radiative forcing magnitudes.
Absorbing aerosols affect the Earth's climate through direct radiative heating of the troposphere. We analyze the tropical tropospheric-only response to a globally uniform increase in black carbon, simulated with an atmospheric general circulation model, in order to gain insight into the interactions that determine the radiative flux perturbation. Over the convective regions, heating in the free troposphere hinders the vertical development of deep cumulus clouds, resulting in the detrainment of more cloudy air into the large-scale environment and stronger cloud reflection. A different mechanism operates over the subsidence regions, where heating near the boundary layer top causes a substantial reduction in low cloud amount thermodynamically by decreasing relative humidity and dynamically by lowering cloud top. These findings, which align well with previous general circulation model and large eddy simulation calculations for black carbon, provide physically based explanations for the main characteristics of the tropical tropospheric adjustment. The implications for quantifying the climate perturbation posed by absorbing aerosols are discussed.
Zhou, C, Joyce Penner, Yi Ming, and X Huang, October 2012: Aerosol forcing based on CAM5 and AM3 meteorological fields. Atmospheric Chemistry and Physics, 12(20), DOI:10.5194/acp-12-9629-2012. Abstract
We use a single aerosol model to explore the effects of the differing meteorological fields from the NCAR CAM5 and GFDL AM3 models. We simulate the global distributions of sulfate, black carbon, organic matter, dust and sea salt using the University of Michigan IMPACT model and use these fields to calculate aerosol direct and indirect forcing, thereby isolating the impacts of the differing meteorological fields.
Over all, the IMPACT-AM3 model predicts larger burdens and longer aerosol lifetimes than the IMPACT-CAM5 model. However, the IMPACT-CAM5 simulations transport more black carbon to the polar regions and more dust from Asia towards North America. These differences can mainly be attributed to differences in: (1) the vertical cloud mass flux and large-scale precipitation fields which determine the wet deposition of aerosols; (2) the in-cloud liquid water content and the cloud coverage which determine the wet aqueous phase production of sulfate. The burden, lifetime and global distribution, especially black carbon in polar regions, are strongly affected by choice of the parameters used for wet deposition.
The total annual mean aerosol optical depth (AOD) at 550 nm ranges from 0.087 to 0.122 for the IMPACT-AM3 model and from 0.138 to 0.186 for the IMPACT-CAM5 model (range is due to different parameters used for wet deposition). Even though IMPACT-CAM5 has smaller aerosol burdens, its AOD is larger due to the much higher relative humidity in CAM5 which leads to more hygroscopic growth. The corresponding global annual average anthropogenic and all-sky aerosol direct forcing at the top of the atmosphere ranges from −0.25 W m−2 to −0.30 W m−2 for IMPACT-AM3 and from −0.48 W m−2 to −0.64 W m−2 for IMPACT-CAM5. The global annual average anthropogenic 1st aerosol indirect effect at the top of the atmosphere ranges from −1.26 W m−2 to −1.44 W m−2 for IMPACT-AM3 and from −1.74 W m−2 to −1.77 W m−2 for IMPACT-CAM5.
Observations show that South Asia underwent a widespread summertime drying during the
second half of the 20th century, but it is unclear whether this trend was due to natural variations or
human activities. We used a series of climate model experiments to investigate the South Asian
monsoon response to natural and anthropogenic forcings. We find that the observed precipitation
decrease can be attributed mainly to human-influenced aerosol emissions. The drying is a
robust outcome of a slowdown of the tropical meridional overturning circulation, which
compensates for the aerosol-induced energy imbalance between the Northern and Southern
Hemispheres. These results provide compelling evidence of the prominent role of aerosols in
shaping regional climate change over South Asia.
Chen, Gang, Yi Ming, N D Singer, and Jian Lu, February 2011: Testing the Clausius-Clapeyron constraint on the aerosol-induced changes in mean and extreme precipitation. Geophysical Research Letters, 38, L04807, DOI:10.1029/2010GL046435. Abstract
The impacts of aerosol and greenhouse gas forcing of the 20th century on the climatological mean and extremes of precipitation are compared in an atmospheric GCM with improved parameterizations of aerosol direct and indirect effects. In spite of different forcing patterns, the thermodynamic effects of aerosol cooling and greenhouse gas warming on the zonally averaged precipitation have similar latitudinal patterns but opposite signs, plausibly due to the effects of temperature on atmospheric water vapor content and moisture flux. The fractional thermodynamic change, for both the moisture convergence in mid- and high latitudes and precipitation extremes at all latitudes, scales linearly with the abundance of atmospheric moisture at a rate of ∼5%/K, somewhat less than the expectation from the Clausius-Clapeyron relation.
The Geophysical Fluid Dynamics Laboratory (GFDL) has developed a coupled general circulation model (CM3) for atmosphere, oceans, land, and sea ice. The goal of CM3 is to address emerging issues in climate change, including aerosol-cloud interactions, chemistry-climate interactions, and coupling between the troposphere and stratosphere. The model is also designed to serve as the physical-system component of earth-system models and models for decadal prediction in the near-term future, for example, through improved simulations in tropical land precipitation relative to earlier-generation GFDL models. This paper describes the dynamical core, physical parameterizations, and basic simulation characteristics of the atmospheric component (AM3) of this model.
Relative to GFDL AM2, AM3 includes new treatments of deep and shallow cumulus convection, cloud-droplet activation by aerosols, sub-grid variability of stratiform vertical velocities for droplet activation, and atmospheric chemistry driven by emissions with advective, convective, and turbulent transport. AM3 employs a cubed-sphere implementation of a finite-volume dynamical core and is coupled to LM3, a new land model with eco-system dynamics and hydrology.
Most basic circulation features in AM3 are simulated as realistically, or more so, than in AM2. In particular, dry biases have been reduced over South America. In coupled mode, the simulation of Arctic sea ice concentration has improved. AM3 aerosol optical depths, scattering properties, and surface clear-sky downward shortwave radiation are more realistic than in AM2. The simulation of marine stratocumulus decks and the intensity distributions of precipitation remain problematic, as in AM2.
The last two decades of the 20th century warm in CM3 by .32°C relative to 1881-1920. The Climate Research Unit (CRU) and Goddard Institute for Space Studies analyses of observations show warming of .56°C and .52°C, respectively, over this period. CM3 includes anthropogenic cooling by aerosol cloud interactions, and its warming by late 20th century is somewhat less realistic than in CM2.1, which warmed .66°C but did not include aerosol cloud interactions. The improved simulation of the direct aerosol effect (apparent in surface clear-sky downward radiation) in CM3 evidently acts in concert with its simulation of cloud-aerosol interactions to limit greenhouse gas warming in a way that is consistent with observed global temperature changes.
Ghan, S, and Yi Ming, et al., October 2011: Droplet nucleation: Physically-based parameterizations and comparative evaluation. Journal of Advances in Modeling Earth Systems, 3, M10001, DOI:10.1029/2011MS000074. Abstract
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, which involves the spontaneous growth of aerosol into cloud droplets at cloud edges, during the early stages of cloud formation, and in some cases within the interior of mature clouds. Numerical models of droplet nucleation represent much of the complexity of the process, but at a computational cost that limits their application to simulations of hours or days. Physically-based parameterizations of droplet nucleation are designed to quickly estimate the number nucleated as a function of the primary controlling parameters: the aerosol number size distribution, hygroscopicity and cooling rate. Here we compare and contrast the key assumptions used in developing each of the most popular parameterizations and compare their performances under a variety of conditions. We find that the more complex parameterizations perform well under a wider variety of nucleation 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 compare estimates of anthropogenic aerosol indirect effects using two different parameterizations applied to the same global climate model, and find that the estimates of indirect effects differ by only 10%. We conclude with a summary of the outstanding challenges remaining for further development and application.
The recently developed GFDL Atmospheric Model version 3 (AM3), an atmospheric general circulation model (GCM), incorporates a prognostic treatment of cloud drop number to simulate the aerosol indirect effect. Since cloud drop activation depends on cloud-scale vertical velocities, which are not reproduced in present-day GCMs, additional assumptions on the subgrid variability are required to implement a local activation parameterization into a GCM.
This paper describes the subgrid activation assumptions in AM3 and explores sensitivities by constructing alternate configurations. These alternate model configurations exhibit only small differences in their present-day climatology. However, the total anthropogenic radiative flux perturbation (RFP) between present-day and preindustrial conditions varies by ±50% from the reference, because of a large difference in the magnitude of the aerosol indirect effect. The spread in RFP does not originate directly from the subgrid assumptions but indirectly through the cloud retuning necessary to maintain a realistic radiation balance. In particular, the paper shows a linear correlation between the choice of autoconversion threshold radius and the RFP.
Climate sensitivity changes only minimally between the reference and alternate configurations. If implemented in a fully coupled model, these alternate configurations would therefore likely produce substantially different warming from preindustrial to present day.
We study how anthropogenic aerosols, alone or in conjunction with radiatively active gases, affect the tropical circulation with an atmosphere/mixed layer ocean general circulation model. Aerosol-induced cooling gives rise to a substantial increase in the overall strength of the tropical circulation, a robust outcome consistent with a thermodynamical scaling argument. Owing to the interhemispheric asymmetry in aerosol forcing, the zonal-mean and zonally asymmetrical components of the tropical circulation respond differently. The Hadley circulation weakens in the Northern Hemisphere, but strengthens in the Southern Hemisphere. The resulting northward cross-equatorial moist static energy flux compensates partly for the aerosol radiative cooling in the Northern Hemisphere. In contrast, the less restricted zonally asymmetrical circulation does not show sensitivity to the spatial structure of aerosols, and strengthens in both hemispheres. Our results also point to the possible role of aerosols in driving the observed reduction in the equatorial sea level pressure gradient.
These circulation changes have profound implications for the hydrological cycle. We find that aerosols alone make the subtropical dry zones in both hemispheres wetter, as the local hydrological response is controlled thermodynamically by atmospheric moisture content. The deep tropical rainfall undergoes a dynamically induced southward shift, a robust pattern consistent with the adjustments in the zonal-mean circulation and in the meridional moist static energy transport. Less certain is the magnitude of the shift. The nonlinearity exhibited by the combined hydrological response to aerosols and radiatively active gases is dynamical in nature.
Ming, Yi, V Ramaswamy, and Gang Chen, December 2011: A model investigation of aerosol-induced changes in boreal winter extratropical circulation. Journal of Climate, 24(23), DOI:10.1175/2011JCLI4111.1. Abstract
We examine the key characteristics of the boreal winter extratropical circulation changes in response to anthropogenic aerosols, simulated with a coupled atmosphere-slab ocean general circulation model. The zonal-mean response features a pronounced equatorward shift of the Northern Hemisphere subtropical jet owing to the mid-latitude aerosol cooling. The circulation changes also show strong zonal asymmetry. In particular, the cooling is more concentrated over the North Pacific than over the North Atlantic despite similar regional forcings. With the help of an idealized model, we demonstrate that the zonally asymmetrical response is linked tightly to the stationary Rossby waves excited by the anomalous diabatic heating over the tropical East Pacific. The altered wave pattern leads to a southeastward shift of the Aleutian low (and associated changes in winds and precipitation), while leaving the North Atlantic circulation relatively unchanged.
Despite the rich circulation changes, the variations in the extratropical meridional latent heat transport are controlled strongly by the dependence of atmospheric moisture content on temperature. This suggests that one can project reliably the changes in extratropical zonal-mean precipitation solely from the global-mean temperature change, even without a good knowledge of the detailed circulation changes caused by aerosols. On the other hand, such knowledge is indispensable for understanding zonally asymmetrical (regional) precipitation changes.
Ming, Yi, March 2011: Aerosols In Encyclopedia of Climate and Weather, second edition, Oxford University Press, 34-38.
Ming, Yi, V Ramaswamy, and Geeta Persad, July 2010: Two opposing effects of absorbing aerosols on global-mean precipitation. Geophysical Research Letters, 37, L13701, DOI:10.1029/2010GL042895. Abstract
Absorbing aerosols affect global-mean precipitation primarily in two ways. They give rise to stronger shortwave atmospheric heating, which acts to suppress precipitation. Depending on the top-of-the-atmosphere radiative flux change, they can also warm up the surface with a tendency to increase precipitation. Here, we present a theoretical framework that takes into account both effects, and apply it to analyze the hydrological responses to increased black carbon burden simulated with a general circulation model. It is found that the damping effect of atmospheric heating can outweigh the enhancing effect of surface warming, resulting in a net decrease in precipitation. The implications for moist convection and general circulation are discussed.
A new stratiform cloud scheme including a two-moment bulk microphysics module, a cloud cover parameterization allowing ice supersaturation, and an ice nucleation parameterization has been implemented into the recently developed GFDL AM3 general circulation model (GCM) as part of an effort to treat aerosol-cloud-radiation interactions more realistically. Unlike the original scheme, the new scheme facilitates the study of cloud-ice-aerosol interactions via influences of dust and sulfate on ice nucleation. While liquid and cloud ice water path associated with stratiform clouds are similar for the new and the original scheme, column integrated droplet numbers and global frequency distributions (PDFs) of droplet effective radii differ significantly. This difference is in part due to a difference in the implementation of the Wegener-Bergeron-Findeisen (WBF) mechanism, which leads to a larger contribution from super-cooled droplets in the original scheme. Clouds are more likely to be either completely glaciated or liquid due to the WBF mechanism in the new scheme. Super-saturations over ice simulated with the new scheme are in qualitative agreement with observations, and PDFs of ice numbers and effective radii appear reasonable in the light of observations. Especially, the temperature dependence of ice numbers qualitatively agrees with in-situ observations. The global average long-wave cloud forcing decreases in comparison to the original scheme as expected when super-saturation over ice is allowed. Anthropogenic aerosols lead to a larger decrease in short-wave absorption (SWABS) in the new model setup, but outgoing long-wave radiation (OLR) decreases as well, so that the net effect of including anthropogenic aerosols on the net radiation at the top of the atmosphere (netradTOA = SWABS-OLR) is of similar magnitude for the new and the original scheme.
Shindell, Drew, M Schulz, Yi Ming, Toshihiko Takemura, G Faluvegi, and V Ramaswamy, October 2010: Spatial scales of climate response to inhomogeneous radiative forcing. Journal of Geophysical Research: Atmospheres, 115, D19110, DOI:10.1029/2010JD014108. Abstract
The distances over which localized radiative forcing influences surface temperature have
not been well characterized. We present a general methodology to analyze the spatial
scales of the forcing/response relationship, and apply it to simulations of historical
aerosol forcing and response in four climate models. We find that the surface temperature
response is not strongly sensitive to the longitude of forcing, but is fairly sensitive to
latitude. Surface temperature responses in the Arctic and the Southern Hemisphere
extratropics, where forcing was small, show little relationship to local forcing. Restricting
the analysis to 30S-60N, where nearly all the forcing was applied, shows that forcing
strongly influences response out to ~4500 km away examining all directions. The
meridional length of influence is somewhat shorter (~3500 km or 30 degrees), while it
extends out to at least 12,000 km in the zonal direction. Substantial divergences between
the models are seen over the oceans, whose physical representations differ greatly among
the models. Length scales are quite consistent over 30S-60N land areas, however, despite
differences in both the forcing applied and the physics of the models themselves. The
results suggest that better understanding of regionally inhomogeneous radiative forcing
would lead to improved projections of regional climate change over land areas. They also
provide quantitative estimates of the spatial extent of the climate impacts of pollutants,
which can extend thousands of kilometers beyond polluted areas, especially in the zonal
direction.
Magi, B I., Paul Ginoux, Yi Ming, and V Ramaswamy, July 2009: Evaluation of tropical and extratropical Southern Hemisphere African aerosol properties simulated by a climate model. Journal of Geophysical Research, D14204, DOI:10.1029/2008JD011128. Abstract
We compare aerosol optical depth (AOD) and single scattering albedo (SSA) simulated by updated configurations of a version of the atmospheric model (AM2) component of the NOAA Geophysical Fluid Dynamics Laboratory general circulation model over Southern Hemisphere Africa with AOD and SSA derived from research aircraft measurements and NASA Aerosol Robotic Network (AERONET) stations and with regional AOD from the NASA Moderate Resolution Imaging Spectroradiometer satellite. The results of the comparisons suggest that AM2 AOD is biased low by 30–40% in the tropics and 0–20% in the extratropics, while AM2 SSA is biased high by 4–8%. The AM2 SSA bias is higher during the biomass burning season, and the monthly variations in AM2 SSA are poorly correlated with AERONET. On the basis of a comparison of aerosol mass in the models with measurements from southern Africa, and a detailed analysis of aerosol treatment in AM2, we suggest that the low bias in AOD and high bias in SSA are related to an underestimate of carbonaceous aerosol emissions in the biomass burning inventories used by AM2. Increases in organic matter and black carbon emissions by factors of 1.6 and 3.8 over southern Africa improve the biases in AOD and especially SSA. We estimate that the AM2 biases in AOD and SSA imply that the magnitude of annual top of the atmosphere radiative forcing in clear-sky conditions over southern Africa is overestimated (too negative) by ∼8% while surface radiative forcing is underestimated (not negative enough) by ∼20%.
The equilibrium temperature and hydrological responses to the total aerosol effects (i.e., direct, semidirect, and indirect effects) are studied using a modified version of the Geophysical Fluid Dynamics Laboratory atmosphere general circulation model (AM2.1) coupled to a mixed layer ocean model. The treatment of aerosol–liquid cloud interactions and associated indirect effects is based upon a prognostic scheme of cloud droplet number concentration, with an explicit representation of cloud condensation nuclei activation involving sulfate, organic carbon, and sea salt aerosols. Increasing aerosols from preindustrial (1860) to presentday (1990) levels leads to a decrease of 1.9 K in the global annual mean surface temperature. The cooling is relatively strong over the Northern Hemisphere midlatitude land owing to the high aerosol burden there, while being amplified at high latitudes. When being subject to aerosols and radiatively active gases (i.e., wellmixed greenhouse gases and ozone) simultaneously, the model climate behaves nonlinearly; the simulated increase in surface temperature (0.55 K) is considerably less than the arithmetic sum of separate aerosol and gas effects (0.86 K). The thermal responses are accompanied by the nonlinear changes in cloud fields, which are amplified owing to the surface albedo feedback at high latitudes. The two effects completely offset each other in the Northern Hemisphere, while gas effect is dominant in the Southern Hemisphere. Both factors are crucial in shaping the regional responses. Interhemispheric asymmetry in aerosol-induced cooling yields a southward shift of the intertropical convergence zone, thus giving rise to a significant reduction in precipitation north of the equator, and an increase to the south. The simulations show that the change of precipitation in response to the simultaneous increases in aerosols and gases not only largely follows the same pattern as that for aerosols alone, but that it is also substantially strengthened in terms of magnitude south of 10°N. This is quite different from the damping expected from adding up individual responses, and further indicates the nonlinearity in the model's hydrological response.
Quaas, Johannes, Yi Ming, Leo J Donner, and Paul Ginoux, et al., November 2009: Aerosol indirect effects – general circulation model intercomparison and evaluation with satellite data. Atmospheric Chemistry and Physics, 9(22), DOI:10.5194/acp-9-8697-2009. Abstract
Aerosol indirect effects continue to constitute one of the most important uncertainties for anthropogenic climate perturbations. Within the international AEROCOM initiative, the representation of aerosol-cloud-radiation interactions in ten different general circulation models (GCMs) is evaluated using three satellite datasets. The focus is on stratiform liquid water clouds since most GCMs do not include ice nucleation effects, and none of the model explicitly parameterises aerosol effects on convective clouds. We compute statistical relationships between aerosol optical depth (τa) and various cloud and radiation quantities in a manner that is consistent between the models and the satellite data. It is found that the model-simulated influence of aerosols on cloud droplet number concentration (Nd) compares relatively well to the satellite data at least over the ocean. The relationship between τa and liquid water path is simulated much too strongly by the models. This suggests that the implementation of the second aerosol indirect effect mainly in terms of an autoconversion parameterisation has to be revisited in the GCMs. A positive relationship between total cloud fraction (fcld) and τa as found in the satellite data is simulated by the majority of the models, albeit less strongly than that in the satellite data in most of them. In a discussion of the hypotheses proposed in the literature to explain the satellite-derived strong fcld–τa relationship, our results indicate that none can be identified as a unique explanation. Relationships similar to the ones found in satellite data between τa and cloud top temperature or outgoing long-wave radiation (OLR) are simulated by only a few GCMs. The GCMs that simulate a negative OLR–τa relationship show a strong positive correlation between τa and fcld. The short-wave total aerosol radiative forcing as simulated by the GCMs is strongly influenced by the simulated anthropogenic fraction of τa, and parameterisation assumptions such as a lower bound on Nd. Nevertheless, the strengths of the statistical relationships are good predictors for the aerosol forcings in the models. An estimate of the total short-wave aerosol forcing inferred from the combination of these predictors for the modelled forcings with the satellite-derived statistical relationships yields a global annual mean value of −1.5±0.5 Wm−2. In an alternative approach, the radiative flux perturbation due to anthropogenic aerosols can be broken down into a component over the cloud-free portion of the globe (approximately the aerosol direct effect) and a component over the cloudy portion of the globe (approximately the aerosol indirect effect). An estimate obtained by scaling these simulated clear- and cloudy-sky forcings with estimates of anthropogenic τa and satellite-retrieved Nd–τa regression slopes, respectively, yields a global, annual-mean aerosol direct effect estimate of −0.4±0.2 Wm−2 and a cloudy-sky (aerosol indirect effect) estimate of −0.7±0.5 Wm−2, with a total estimate of −1.2±0.4 Wm−2.
Lee, S S., Leo J Donner, Vaughan T J Phillips, and Yi Ming, 2008: The dependence of aerosol effects on clouds and precipitation on cloud-system organization, shear and stability. Journal of Geophysical Research, 113, D16202, DOI:10.1029/2007JD009224. Abstract
Precipitation suppression due to an increase of aerosol number concentration in stratiform cloud is well-known. It is not certain whether the suppression applies for deep convection. Recent studies have suggested increasing precipitation from deep convection with increasing aerosols under some, but not all, conditions. Increasing precipitation with increasing aerosols can result from strong interactions in deep convection between dynamics and microphysics. High cloud liquid, due to delayed autoconversion, provides more evaporation, leading to more active downdrafts, convergence fields, condensation, collection of cloud liquid by precipitable hydrometeors, and precipitation. Evaporation of cloud liquid is a primary determinant of the intensity of the interactions. It is partly controlled by wind shear modulating the entrainment of dry air into clouds and transport of cloud liquid into unsaturated areas. Downdraft-induced convergence, crucial to the interaction, is weak for shallow clouds, generally associated with low convective available potential energy (CAPE). Aerosol effects on cloud and precipitation can vary with CAPE and wind shear. Pairs of idealized numerical experiments for high and low aerosol cases were run for five different environmental conditions to investigate the dependence of aerosol effect on stability and wind shear. In the environment of high CAPE and strong wind shear, cumulonimbus- and cumulus-type clouds were dominant. Transport of cloud liquid to unsaturated areas was larger at high aerosol, leading to stronger downdrafts. Because of the large vertical extent of those clouds, strong downdrafts and convergence developed for strong interactions between dynamics and microphysics. These led to larger precipitation at high aerosol. Detrainment of cloud liquid and associated evaporation were less with lower CAPE and wind shear, where dynamically weaker clouds dominated. Transport of cloud liquid to unsaturated areas was not as active as in the environment of high CAPE and strong shear. Also, evaporatively driven differences in downdrafts at their level of initial descent were not magnified in clouds with shallow depth as much as in deep convective clouds as they accelerated to the surface over shorter distances. Hence the interaction between dynamics and microphysics was reduced, leading to precipitation suppression at high aerosol. These results demonstrate that increasing aerosol can either decrease or increase precipitation for an imposed large-scale environment supporting cloud development. The implications for larger-scale aspects of the hydrological cycle will require further study with larger-domain models and cumulus parameterizations with advanced microphysics.
Lee, S S., Leo J Donner, Vaughan T J Phillips, and Yi Ming, 2008: Examination of aerosol effects on precipitation in deep convective clouds during the 1997 ARM summer experiment. Quarterly Journal of the Royal Meteorological Society, 134(634), DOI:10.1002/qj.287. Abstract
It has been generally accepted that increasing aerosols suppress precipitation. The aerosol-induced precipitation suppression was suggested by the study of shallow stratiform clouds. Recent studies of convective clouds showed increasing aerosols could increase precipitation. Those studies showed that intense feedbacks between aerosols and cloud dynamics led to increased precipitation in some cases of convective clouds. This study expanded those studies by analyzing detailed microphysical and dynamical modifications by aerosols leading to increased precipitation. This study focused on three observed cases of mesoscale cloud ensemble (MCE) driven by deep convective clouds, since MCE accounts for a large proportion of the Earth's precipitation and the study of aerosol effects on MCE is at its incipient stage. Those MCEs were observed during the 1997 Atmospheric Radiation Measurement (ARM) summer experiment. Two numerical experiments were performed for each of the MCEs to simulate aerosol effects on deep convection. The first was with high aerosol number concentration, and the second was with low concentration. The results showed an increased precipitation at high aerosol, due to stronger, more numerous updraughts, initiated by stronger convergence lines at the surface in convective regions of the MCE. The stronger convergence lines were triggered by increased evaporation of cloud liquid in the high-aerosol case, made possible by higher values of cloud liquid necessary for autoconversion.
The generality of these results requires further investigation. However, they demonstrate that the response of precipitation to increased aerosols in deep convection can be different from that in shallow cloud systems, at least for the cases studied here.
Ming, Yi, V Ramaswamy, Leo J Donner, Vaughan T J Phillips, Stephen A Klein, Paul Ginoux, and Larry W Horowitz, February 2007: Modeling the interactions between aerosols and liquid water clouds with a self-consistent cloud scheme in a general circulation model. Journal of the Atmospheric Sciences, 64(4), DOI:10.1175/JAS3874.1. Abstract
To model aerosol-cloud interactions in general circulation
models (GCMs), a prognostic cloud scheme of cloud liquid water and amount is expanded to include droplet number concentration (Nd) in a way that allows them to be calculated using the same large-scale and convective updraft velocity field. In the scheme, the evolution of droplets fully interacts with the model meteorology. An explicit treatment of cloud condensation nuclei (CCN) activation enables the scheme to take into account the contributions to Nd of multiple aerosol species (i.e., sulfate, organic, and sea-salt aerosols) and to consider kinetic limitations of the activation process. An implementation of the prognostic scheme in the Geophysical Fluid Dynamics Laboratory (GFDL) AM2 GCM yields a vertical distribution of Nd with a characteristic maximum in the lower troposphere; this feature differs from the profile that would be obtained if Ndis diagnosed from the sulfate mass concentration based on an often-used empirical relationship. Prognosticated Nd exhibits large variations with respect to the sulfate mass concentration. The mean values are generally consistent with the empirical relationship over ocean, but show negative biases over the Northern Hemisphere midlatitude land, perhaps owing to the neglect of subgrid variations of large-scale ascents and inadequate convective sources. The prognostic scheme leads to a substantial improvement in the agreement of model-predicted present-day liquid water path (LWP) and cloud forcing with satellite measurements compared to using the empirical relationship.
The simulations with preindustrial and present-day aerosols show that the
combined first and second indirect effects of anthropogenic sulfate and organic aerosols give rise to a steady-state global annual mean flux change of -1.8 W m-2, consisting of -2.0 W m-2 in shortwave and 0.2 W m-2 in longwave. The ratios of the flux changes in the Northern Hemisphere (NH) to that in Southern Hemisphere (SH) and of the flux changes over ocean to that over land are 2.9 and 0.73, respectively. These estimates are consistent with the averages of values from previous studies stated in a recent review. The model response to higher Nd alters the cloud field; LWP and total cloud amount increase by 19% and 0.6%, respectively. Largely owing to high sulfate concentrations from fossil fuel burning, the NH midlatitude land and oceans experience strong radiative cooling. So does the tropical land, which is dominated by biomass burning-derived organic aerosol. The computed annual, zonal-mean flux changes are determined to be statistically significant, exceeding the model's natural variations in the NH low and midlatitudes and in the SH low latitudes. This study reaffirms the major role of sulfate in providing CCN for cloud formation.
Ming, Yi, V Ramaswamy, Leo J Donner, and Vaughan T J Phillips, 2006: A new parameterization of cloud droplet activation applicable to general circulation models. Journal of the Atmospheric Sciences, 63(4), DOI:10.1175/JAS3686.1. Abstract
A new parameterization is proposed to link the droplet number concentration to the size distribution and chemical composition of aerosol and updraft velocity. Except for an empirical assumption of droplet growth, the parameterization is formulated almost entirely on first principles to allow for satisfactory performance under a variety of conditions. For a series of updraft velocity ranging from 0.03 to 10.0 m s−1, the droplet number concentrations predicted with the parameterization are in good agreement with the detailed parcel model simulations with an average error of −4 ± 26% (one standard deviation). The accuracy is comparable to or better than some existing parameterizations. The parameterization is able to account for the effects of droplet surface tension and mass accommodation coefficient on activation without adjusting the empirical parameter. These desirable attributes make the parameterization suitable for being used in the prognostic determination of the cloud droplet number concentration in general circulation models (GCMs).
This study simulates the direct radiative forcing of organic aerosol using the GFDL AM2 GCM. The aerosol climatology is provided by the MOZART chemical transport model (CTM). The approach to calculating aerosol optical properties explicitly considers relative humidity–dependent hygroscopic growth by employing a functional group–based thermodynamic model, and makes use of the size distribution derived from AERONET measurements. The preindustrial (PI) and present-day (PD) global burdens of organic carbon are 0.17 and 1.36 Tg OC, respectively. The annual global mean total-sky and clear-sky top-of-the atmosphere (TOA) forcings (PI to PD) are estimated as −0.34 and −0.71 W m−2, respectively. Geographically the radiative cooling largely lies over the source regions, namely part of South America, Central Africa, Europe and South and East Asia. The annual global mean total-sky and clear-sky surface forcings are −0.63 and −0.98 W m−2, respectively. A series of sensitivity analyses shows that the treatments of hygroscopic growth and optical properties of organic aerosol are intertwined in the determination of the global organic aerosol forcing. For example, complete deprivation of water uptake by hydrophilic organic particles reduces the standard (total-sky) and clear-sky TOA forcing estimates by 18% and 20%, respectively, while the uptake by a highly soluble organic compound (malonic acid) enhances them by 18% and 32%, respectively. Treating particles as non-absorbing enhances aerosol reflection and increases the total-sky and clear-sky TOA forcing by 47% and 18%, respectively, while neglecting the scattering brought about by the water associated with particles reduces them by 24% and 7%, respectively.
Ming, Yi, V Ramaswamy, Paul Ginoux, Larry W Horowitz, and L M Russell, 2005: Geophysical Fluid Dynamics Laboratory general circulation model investigation of the indirect radiative effects of anthropogenic sulfate aerosol. Journal of Geophysical Research, 110, D22206, DOI:10.1029/2005JD006161. Abstract
The Geophysical Fluid Dynamics Laboratory (GFDL) atmosphere general circulation model, with its new cloud scheme, is employed to study the indirect radiative effect of anthropogenic sulfate aerosol during the industrial period. The preindustrial and present-day monthly mean aerosol climatologies are generated from running the Model for Ozone And Related chemical Tracers (MOZART) chemistry-transport model. The respective global annual mean sulfate burdens are 0.22 and 0.81 Tg S. Cloud droplet number concentrations are related to sulfate mass concentrations using an empirical relationship (Boucher and Lohmann, 1995). A distinction is made between "forcing" and flux change at the top of the atmosphere in this study. The simulations, performed with prescribed sea surface temperature, show that the first indirect "forcing" ("Twomey" effect) amounts to an annual mean of -1.5 W m-2, concentrated largely over the oceans in the Northern Hemisphere (NH). The annual mean flux change owing to the response of the model to the first indirect effect is -1.4 W m-2, similar to the annual mean forcing. However, the model's response causes a rearrangement of cloud distribution as well as changes in longwave flux (smaller than solar flux changes). There is thus a differing geographical nature of the radiation field than for the forcing even though the global means are similar. The second indirect effect, which is necessarily an estimate made in terms of the model's response, amounts to -0.9 W m-2, but the statistical significance of the simulated geographical distribution of this effect is relatively low owing to the model's natural variability. Both the first and second effects are approximately linearly additive, giving rise to a combined annual mean flux change of -2.3 W m-2, with the NH responsible for 77% of the total flux change. Statistically significant model responses are obtained for the zonal mean total indirect effect in the entire NH and in the Southern Hemisphere low latitudes and midlatitudes (north of 45°S). The area of significance extends more than for the first and second effects considered separately. A comparison with a number of previous studies based on the same sulfate-droplet relationship shows that, after distinguishing between forcing and flux change, the global mean change in watts per square meter for the total effect computed in this study is comparable to existing studies in spite of the differences in cloud schemes.
Ming, Yi, L M Russell, and D F Bradford, 2005: Health and climate policy impacts on sulfur emission control. Reviews of Geophysics, 43, RG4001, DOI:10.1029/2004RG000167. Abstract
Sulfate aerosol from burning fossil fuels not only has strong cooling effects on the Earth's climate but also imposes substantial costs on human health. To assess the impact of addressing air pollution on climate policy, we incorporate both the climate and health effects of sulfate aerosol into an integrated-assessment model of fossil fuel emission control. Our simulations show that a policy that adjusts fossil fuel and sulfur emissions to address both warming and health simultaneously will support more stringent fossil fuel and sulfur controls. The combination of both climate and health objectives leads to an acceleration of global warming in the 21st century as a result of the short-term climate response to the decreased cooling from the immediate removal of short-lived sulfate aerosol. In the long term (more than 100 years), reducing sulfate aerosol emissions requires that we decrease fossil fuel combustion in general, thereby removing some of the coemitted carbon emissions and leading to a reduction in global warming.
Ming, Yi, and L M Russell, 2004: Organic aerosol effects on fog droplet spectra. Journal of Geophysical Research, 109, D10206, DOI:10.1029/2003JD004427. Abstract
Organic aerosol alters cloud and fog properties through surface tension and solubility effects. This study characterizes the role of organic compounds in affecting fog droplet number concentration by initializing and comparing detailed particle microphysical simulations with two field campaigns in the Po Valley. The size distribution and chemical composition of aerosol were based on the measurements made in the Po Valley Fog Experiments in 1989 and 1998–1999. Two types of aerosol with different hygroscopicity were considered: the less hygroscopic particles, composed mainly of organic compounds, and the more hygroscopic particles, composed mainly of inorganic salts. The organic fraction of aerosol mass was explicitly modeled as a mixture of seven soluble compounds [Fuzzi et al., 2001] by employing a functional group-based thermodynamic model [Ming and Russell, 2002]. Condensable gases in the vapor phase included nitric acid, sulfuric acid, and ammonia. The maximum supersaturation in the simulation is 0.030% and is comparable to the calculation by Noone et al. [1992] inferred from measured residual particle fractions. The minimum activation diameters of the less and more hygroscopic particles are 0.49 μm and 0.40 μm, respectively. The predicted residual particle fractions are in agreement with measurements. The organic components of aerosol account for 34% of the droplet residual particle mass and change the average droplet number concentration by −10–6%, depending on the lowering of droplet surface tension and the interactions among dissolving ions. The hygroscopic growth of particles due to the presence of water-soluble organic compounds enhances the condensation of nitric acid and ammonia due to the increased surface area, resulting in a 9% increase in the average droplet number concentration. Assuming ideal behavior of aqueous solutions of water-soluble organic compounds overestimates the hygroscopic growth of particles and increases droplet numbers by 6%. The results are sensitive to microphysical processes such as condensation of soluble gases, which increases the average droplet number concentration by 26%. Wet deposition plays an important role in controlling liquid water content in this shallow fog.
Ming, Yi, et al., 2004: Free energy perturbation study of water dimer dissociation kinetics. Journal of Chemical Physics, 121, 773-777.
Ming, Yi, and L M Russell, 2002: Thermodynamic equilibrium of aqueous solutions of organic-electrolyte mixtures in aerosol particles. AIChE Journal, 48, 1331.
Russell, L M., and Yi Ming, 2002: Deliquescence of small particles. Journal of Chemical Physics, 116, 311-321.
Ming, Yi, and L M Russell, 2001: Predicted hygroscopic growth of sea salt aerosol. Journal of Geophysical Research, 106(D22), 28,259-28,274. Abstract
Organic species in sea salt particles can significantly reduce hygroscopic growth in subsaturated conditions, an important uncertainty in the radiative effect of aerosol particles on the atmosphere. This hygroscopic behavior is predicted with a numerical model of the the organic‐water, electrolyte‐water, and organicelectrolyte interactions in complex mixtures of organic species and inorganic ions. The results show a 15% decrease in hygroscopic growth above 75% relative humidity for particles that include as little as 30% organic mass. Organic compositions of 50% organic mass reduce hygroscopic growth by 25%. This prediction relies on particle chemical composition estimated from measurements of insoluble organic species in marine‐derived particles and of soluble organic species measured in seawater. Twenty insoluble and four soluble organic species are used to represent the behavior of sea salt organic composition. The hygroscopic growth is strongly sensitive to the organic fraction that is soluble or slightly soluble, although variations among different soluble or insoluble species are small above the sodium chloride deliquescence point. Interactions between organic and electrolyte species depend primarily on the “salting out” behavior of NaCl with alkanes, carboxylic acids, and alcohols, although interactions with other inorganic ions in sea salt were estimated to cause small changes in the hygroscopic growth. The predicted growth factors for sea salt with < 30% organic species are consistent with growth factors measured for ambient marine‐derived particles by another group [Berg et al., 1998; Swietlicki et al, 2000; Zhou et al, 2001]. This coincidence suggests that the less‐hygroscopic particles could indicate the presence of marine organic compounds, although multiple combinations of inorganic and anthropogenic organic species would also satisfy the measured behavior.
Prenni, A J., P J DeMott, Marta Abalos, D E Sherman, L M Russell, and Yi Ming, 2001: The effects of low molecular weight dicarboxylic acids on cloud formation. Journal of Physical Chemistry A, 105, 11240-11248.
