Beadling, Rebecca L., Pu Lin, John P Krasting, William Ellinger, Anna Coomans, James Milward, Katherine Turner, Xiaoqi Xu, Torge Martin, and Maria J Molina, November 2024: From the surface to the stratosphere: Large-scale atmospheric response to Antarctic meltwater. Geophysical Research Letters, 51(21), DOI:10.1029/2024GL110157. Abstract
The ocean response to Antarctic Ice Sheet (AIS) mass loss has been extensively studied using numerical models, but less attention has been given to the atmosphere. We examine the global atmospheric response to AIS meltwater in an ensemble of experiments performed using two fully coupled climate models under a pre-industrial climate. In response to AIS meltwater, the experiments yield cooling from the surface to the tropopause over the subpolar Southern Ocean, warming in the Southern Hemisphere polar stratosphere, and cooling in the upper tropical troposphere. Positive feedbacks, initiated by disrupted ocean-atmosphere heat exchange, result in a change in the top-of-atmosphere radiative balance caused primarily through surface and near-surface albedo changes. Changes in the atmospheric thermal structure alter the jet streams aloft. The results highlight the global influence of AIS melting on the climate system and the potential for impacts on mid-latitude climate patterns and delayed regional warming signals.
Clark, Joseph P., Pu Lin, and Spencer A Hill, June 2024: ITCZ response to disabling parameterized convection in global fixed-SST GFDL-AM4 aquaplanet simulations at 50 and 6 km resolutions. Journal of Advances in Modeling Earth Systems, 16(6), DOI:10.1029/2023MS003968. Abstract
As the community increases climate model horizontal resolutions and experiments with removing moist convective parameterizations entirely, it is imperative to understand how these advances affect the InterTropical Convergence Zone (ITCZ). We investigate how the ITCZ responds to deactivating parameterized convection at two resolutions, 50 and 6 km, in fixed sea surface temperature, aquaplanet simulations with the NOAA GFDL AM4 atmospheric model. Disabling parameterized convection at 50 km resolution narrows the ITCZ and increases its precipitation minus evaporation (P–E) maximum by ∼78%, whereas at 6 km resolution doing so widens the ITCZ and decreases its P–E maximum by ∼50%. Using the column-integrated moist static energy budget, we decompose these tropical P–E responses into contributions from changes in atmospheric energy input (AEI), gross moist stability, and gross moisture stratification. At 6 km, the ITCZ weakens due to increased gross moist stability. Disabling the convective parameterization at this finer resolution deepens the circulation, favoring more efficient poleward energy transport out of the deep tropics and reduced precipitation in the core of the ITCZ. Conversely, at 50 km the ITCZ strengthening is primarily driven by AEI, which in turn stems primarily from increased low cloud amount and thus longwave cloud radiative cooling in the Hadley cell subsiding branch. The Hadley circulation overturning intensifies to produce poleward energy fluxes that compensate the longwave cooling, yielding a stronger ITCZ. We further show that the low level diabatic heating profiles over the descending region are instrumental in understanding such diverse responses.
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.
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.
Martin, Zane K., Isla Simpson, and Pu Lin, et al., September 2023: The lack of a QBO-MJO connection in climate models with a nudged stratosphere. JGR Atmospheres, 128(17), DOI:10.1029/2023JD038722. Abstract
The observed stratospheric quasi-biennial oscillation (QBO) and the tropospheric Madden-Julian oscillation (MJO) are strongly connected in boreal winter, with stronger MJO activity when lower-stratospheric winds are easterly. However, the current generation of climate models with internally generated representations of the QBO and MJO do not simulate the observed QBO-MJO connection, for reasons that remain unclear. This study builds on prior work exploring the QBO-MJO link in climate models whose stratospheric winds are relaxed toward reanalysis, reducing stratospheric biases in the model and imposing a realistic QBO. A series of ensemble experiments are performed using four state-of-the-art climate models capable of representing the MJO over the period 1980–2015, each with similar nudging in the stratosphere. In these four models, nudging leads to a good representation of QBO wind and temperature signals, however no model simulates the observed QBO-MJO relationship. Biases in MJO vertical structure and cloud-radiative feedbacks are investigated, but no conclusive model bias or mechanism is identified that explains the lack of a QBO-MJO connection.
Andrews, Timothy, Alejandro Bodas-Salcedo, Jonathan M Gregory, Yue Dong, Kyle Armour, David J Paynter, and Pu Lin, et al., September 2022: On the effect of historical SST patterns on radiative feedback. JGR Atmospheres, 127(18), DOI:10.1029/2022JD036675. Abstract
We investigate the dependence of radiative feedback on the pattern of sea-surface temperature (SST) change in 14 Atmospheric General Circulation Models (AGCMs) forced with observed variations in SST and sea-ice over the historical record from 1871 to near-present. We find that over 1871–1980, the Earth warmed with feedbacks largely consistent and strongly correlated with long-term climate sensitivity feedbacks (diagnosed from corresponding atmosphere-ocean GCM abrupt-4xCO2 simulations). Post 1980, however, the Earth warmed with unusual trends in tropical Pacific SSTs (enhanced warming in the west, cooling in the east) and cooling in the Southern Ocean that drove climate feedback to be uncorrelated with—and indicating much lower climate sensitivity than—that expected for long-term CO2 increase. We show that these conclusions are not strongly dependent on the Atmospheric Model Intercomparison Project (AMIP) II SST data set used to force the AGCMs, though the magnitude of feedback post 1980 is generally smaller in nine AGCMs forced with alternative HadISST1 SST boundary conditions. We quantify a “pattern effect” (defined as the difference between historical and long-term CO2 feedback) equal to 0.48 ± 0.47 [5%–95%] W m−2 K−1 for the time-period 1871–2010 when the AGCMs are forced with HadISST1 SSTs, or 0.70 ± 0.47 [5%–95%] W m−2 K−1 when forced with AMIP II SSTs. Assessed changes in the Earth's historical energy budget agree with the AGCM feedback estimates. Furthermore satellite observations of changes in top-of-atmosphere radiative fluxes since 1985 suggest that the pattern effect was particularly strong over recent decades but may be waning post 2014.
Lawrence, Zachary D., Marta Abalos, Blanca Ayarzagüena, D Barriopedro, Amy Butler, Natalia Calvo, Alvaro de la Cámara, Andrew J Charlton-Perez, Daniela I V Domeisen, Etienne Dunn-Sigouin, J García-Serrano, Chaim I Garfinkel, Neil P Hindley, Liwei Jia, Martin Jucker, Alexey Y Karpechko, Hera Kim, Andrea L Lang, Simon H Lee, and Pu Lin, et al., August 2022: Quantifying stratospheric biases and identifying their potential sources in subseasonal forecast systems. Weather and Climate Dynamics, 3(3), DOI:10.5194/wcd-3-977-2022977-1001. Abstract
The stratosphere can be a source of predictability for surface weather on timescales of several weeks to months. However, the potential predictive skill gained from stratospheric variability can be limited by biases in the representation of stratospheric processes and the coupling of the stratosphere with surface climate in forecast systems. This study provides a first systematic identification of model biases in the stratosphere across a wide range of subseasonal forecast systems.
It is found that many of the forecast systems considered exhibit warm global-mean temperature biases from the lower to middle stratosphere, too strong/cold wintertime polar vortices, and too cold extratropical upper-troposphere/lower-stratosphere regions. Furthermore, tropical stratospheric anomalies associated with the Quasi-Biennial Oscillation tend to decay toward each system's climatology with lead time. In the Northern Hemisphere (NH), most systems do not capture the seasonal cycle of extreme-vortex-event probabilities, with an underestimation of sudden stratospheric warming events and an overestimation of strong vortex events in January. In the Southern Hemisphere (SH), springtime interannual variability in the polar vortex is generally underestimated, but the timing of the final breakdown of the polar vortex often happens too early in many of the prediction systems.
These stratospheric biases tend to be considerably worse in systems with lower model lid heights. In both hemispheres, most systems with low-top atmospheric models also consistently underestimate the upward wave driving that affects the strength of the stratospheric polar vortex. We expect that the biases identified here will help guide model development for subseasonal-to-seasonal forecast systems and further our understanding of the role of the stratosphere in predictive skill in the troposphere.
Morgenstern, Olaf, Douglas Kinnison, Michael J Mills, Martine Michou, Larry W Horowitz, and Pu Lin, et al., October 2022: Comparison of Arctic and Antarctic stratospheric climates in chemistry versus no-chemistry climate models. JGR Atmospheres, 127(20), DOI:10.1029/2022JD037123. Abstract
Using nine chemistry-climate and eight associated no-chemistry models, we investigate the persistence and timing of cold episodes occurring in the Arctic and Antarctic stratosphere during the period 1980–2014. We find systematic differences in behavior between members of these model pairs. In a first group of chemistry models whose dynamical configurations mirror their no-chemistry counterparts, we find an increased persistence of such cold polar vortices, such that these cold episodes often start earlier and last longer, relative to the times of occurrence of the lowest temperatures. Also the date of occurrence of the lowest temperatures, both in the Arctic and the Antarctic, is often delayed by 1–3 weeks in chemistry models, versus their no-chemistry counterparts. This behavior exacerbates a widespread problem occurring in most or all models, a delayed occurrence, in the median, of the most anomalously cold day during such cold winters. In a second group of model pairs there are differences beyond just ozone chemistry. In particular, here the chemistry models feature more levels in the stratosphere, a raised model top, and differences in non-orographic gravity wave drag versus their no-chemistry counterparts. Such additional dynamical differences can completely mask the above influence of ozone chemistry. The results point toward a need to retune chemistry-climate models versus their no-chemistry counterparts.
Abalos, Marta, Natalia Calvo, Samuel Benito-Barca, Hella Garny, Steven C Hardiman, and Pu Lin, et al., September 2021: The Brewer-Dobson circulation in CMIP6. Atmospheric Chemistry and Physics, 21(17), DOI:10.5194/acp-21-13571-202113571-13591. Abstract
The Brewer–Dobson circulation (BDC) is a key feature of the stratosphere that models need to accurately represent in order to simulate surface climate variability and change adequately. For the first time, the Climate Model Intercomparison Project includes in its phase 6 (CMIP6) a set of diagnostics that allow for careful evaluation of the BDC. Here, the BDC is evaluated against observations and reanalyses using historical simulations. CMIP6 results confirm the well-known inconsistency in the sign of BDC trends between observations and models in the middle and upper stratosphere. Nevertheless, the large uncertainty in the observational trend estimates opens the door to compatibility. In particular, when accounting for the limited sampling of the observations, model and observational trend error bars overlap in 40 % of the simulations with available output. The increasing CO2 simulations feature an acceleration of the BDC but reveal a large spread in the middle-to-upper stratospheric trends, possibly related to the parameterized gravity wave forcing. The very close connection between the shallow branch of the residual circulation and surface temperature is highlighted, which is absent in the deep branch. The trends in mean age of air are shown to be more robust throughout the stratosphere than those in the residual circulation.
Freidenreich, Stuart, David J Paynter, Pu Lin, V Ramaswamy, Alexandra L Jones, Daniel Feldman, and William D Collins, June 2021: An investigation into biases in instantaneous aerosol radiative effects calculated by shortwave parameterizations in two Earth system models. JGR Atmospheres, 126(11), DOI:10.1029/2019JD032323. Abstract
Because the forcings to which Coupled Model Intercomparison Project - Phase 5 (CMIP5) models were subjected were poorly quantified, recent efforts from the Radiative Forcing Model Intercomparison Project (RFMIP) have focused on developing and testing models with exacting benchmarks. Here, we focus on aerosol forcing to understand if for a given distribution of aerosols, participating models are producing a radiometrically-accurate forcing. We apply the RFMIP experimental protocol for assessing flux biases in aerosol instantaneous radiative effect (IRE) on two participating models, GFDL AM4 and CESM 1.2.2. The latter model contains the RRTMG radiation code which is widely used among CMIP6 GCM's. We conduct a series of calculations that test different potential sources of error in these models relative to line-by-line benchmarks. We find two primary sources of error: two-stream solution methods and the techniques to resolve spectral dependencies of absorption and scattering across the solar spectrum. The former is the dominant source of error for both models but the latter is more significant as a contributing factor for CESM 1.2.2. Either source of error can be addressed in future model development, and these results both demonstrate how the RFMIP protocol can help determine the origins of parameterized errors relative to their equivalent benchmark calculations for participating models, as well as highlight a viable path towards a more rigorous quantification and control of forcings for future CMIP exercises.
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.
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.
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.
Ayarzagüena, Blanca, Andrew J Charlton-Perez, Amy Butler, Peter Hitchcock, Isla Simpson, Lorenzo M Polvani, Neal Butchart, Edwin P Gerber, Lesley Gray, Birgit Hassler, and Pu Lin, et al., March 2020: Uncertainty in the Response of Sudden Stratospheric Warmings and Stratosphere‐Troposphere Coupling to Quadrupled CO2 Concentrations in CMIP6 Models. JGR Atmospheres, 125(6), DOI:10.1029/2019JD032345. Abstract
Major sudden stratospheric warmings (SSWs), vortex formation, and final breakdown dates are key highlight points of the stratospheric polar vortex. These phenomena are relevant for stratosphere‐troposphere coupling, which explains the interest in understanding their future changes. However, up to now, there is not a clear consensus on which projected changes to the polar vortex are robust, particularly in the Northern Hemisphere, possibly due to short data record or relatively moderate CO2 forcing. The new simulations performed under the Coupled Model Intercomparison Project, Phase 6, together with the long daily data requirements of the DynVarMIP project in preindustrial and quadrupled CO2 (4xCO2) forcing simulations provide a new opportunity to revisit this topic by overcoming the limitations mentioned above. In this study, we analyze this new model output to document the change, if any, in the frequency of SSWs under 4xCO2 forcing. Our analysis reveals a large disagreement across the models as to the sign of this change, even though most models show a statistically significant change. As for the near‐surface response to SSWs, the models, however, are in good agreement as to this signal over the North Atlantic: There is no indication of a change under 4xCO2 forcing. Over the Pacific, however, the change is more uncertain, with some indication that there will be a larger mean response. Finally, the models show robust changes to the seasonal cycle in the stratosphere. Specifically, we find a longer duration of the stratospheric polar vortex and thus a longer season of stratosphere‐troposphere coupling.
Fu, Qiang, Rachel H White, Mingcheng Wang, B Alexander, S Solomon, Andrew Gettelman, David S Battisti, and Pu Lin, March 2020: The Brewer‐Dobson Circulation during the Last Glacial Maximum. Geophysical Research Letters, 47(5), DOI:10.1029/2019GL086271. Abstract
The Brewer‐Dobson circulation (BDC) during the last glacial maximum (LGM) is investigated in simulations using the Whole Atmosphere Community Climate Model version 6 (WACCM6). We examine vertical mass fluxes, age of stratospheric air, and the transformed Eulerian mean stream function and find that the modeled annual‐mean BDC during the LGM is almost everywhere slower than that in the modern climate (with or without anthropogenic ozone depleting substances). Compared to the modern climate, the annual‐mean tropical upwelling in the LGM is 11.3‐16.9%, 11.2‐15.8%, and 4.4‐10.2% weaker, respectively, at 100, 70, and 30 hPa. Simulated decreases in annual‐mean mass fluxes at 70 and 100 hPa are caused by a weaker parameterized orographic gravity wave drag and resolved wave drag, respectively.
Guan, Weina, Xianan Jiang, Xuejuan Ren, Gang Chen, Pu Lin, and Hai Lin, November 2020: The leading intraseasonal variability mode of wintertime surface air temperature over the North American sector. Journal of Climate, 33(21), DOI:10.1175/JCLI-D-20-0096.1. Abstract
In this study, detailed characteristics of the leading intraseasonal variability mode of boreal winter surface air temperature (SAT) over the North American (NA) sector are investigated. This intraseasonal SAT mode, characterized by two anomalous centers with an opposite sign—one over central NA and another over east Siberia (ES)/Alaska—bears a great resemblance to the “warm Arctic–cold continent” pattern of the interannual SAT variability over NA. This intraseasonal SAT mode and associated circulation exert a pronounced influence on regional weather extremes, including precipitation over the northwest coast of NA, sea ice concentration over the Chukchi and Bering Seas, and extreme warm and cold events over the NA continent and Arctic region. Surface warming and cooling signals of the intraseasonal SAT mode are connected to temperature anomalies in a deep-tropospheric layer up to 300 hPa with a decreasing amplitude with altitude. Particularly, a coupling between the troposphere and stratosphere is found during evolution of the intraseasonal SAT variability, although whether the stratospheric processes are essential in sustaining the leading intraseasonal SAT mode is difficult to determine based on observations alone. Two origins of wave sources are identified in contributing to vertically propagating planetary waves near Alaska: one over ES/Alaska associated with local intraseasonal variability and another from the subtropical North Pacific via Rossby wave trains induced by tropical convective activity over the western Pacific, possibly associated with the Madden–Julian oscillation.
We describe the baseline model configuration and simulation characteristics of the Geophysical Fluid Dynamics Laboratory (GFDL)'s Atmosphere Model version 4.1 (AM4.1), which builds on developments at GFDL over 2013–2018 for coupled carbon‐chemistry‐climate simulation as part of the sixth phase of the Coupled Model Intercomparison Project. In contrast with GFDL's AM4.0 development effort, which focused on physical and aerosol interactions and which is used as the atmospheric component of CM4.0, AM4.1 focuses on comprehensiveness of Earth system interactions. Key features of this model include doubled horizontal resolution of the atmosphere (~200 to ~100 km) with revised dynamics and physics from GFDL's previous‐generation AM3 atmospheric chemistry‐climate model. AM4.1 features improved representation of atmospheric chemical composition, including aerosol and aerosol precursor emissions, key land‐atmosphere interactions, comprehensive land‐atmosphere‐ocean cycling of dust and iron, and interactive ocean‐atmosphere cycling of reactive nitrogen. AM4.1 provides vast improvements in fidelity over AM3, captures most of AM4.0's baseline simulations characteristics, and notably improves on AM4.0 in the representation of aerosols over the Southern Ocean, India, and China—even with its interactive chemistry representation—and in its manifestation of sudden stratospheric warmings in the coldest months. Distributions of reactive nitrogen and sulfur species, carbon monoxide, and ozone are all substantially improved over AM3. Fidelity concerns include degradation of upper atmosphere equatorial winds and of aerosols in some regions.
Fu, Qiang, S Solomon, H A Pahlavan, and Pu Lin, November 2019: Observed changes in Brewer–Dobson circulation for 1980–2018. Environmental Research Letters, 14(11), DOI:10.1088/1748-9326/ab4de7. Abstract
Previous work has examined the Brewer–Dobson circulation (BDC) changes for 1980–2009 based on satellite Microwave Sounding Unit (MSU/AMSU) lower-stratospheric temperature (T LS ) observations and ERA-Interim reanalysis data. Here we examine the BDC changes for the longer period now available (1980–2018), which also allows analysis of both the ozone depletion (1980–1999) and ozone healing (2000–2018) periods. We provide observational evidence that the annual mean BDC accelerated for 1980–1999 but decelerated for 2000–2018, with the changes largely driven by the Southern Hemisphere (SH), which might be partly contributed by the effects of ozone depletion and healing. We also show that the annual mean BDC has accelerated in the last 40 years (at the 90% confidence level) with a relative strengthening of ~1.7% per decade. This overall acceleration was driven by both Northern Hemisphere (40%) and SH (60%) cells. Significant SH radiative warming is also identified in September for 2000–2018 after excluding the year 2002 when a very rare SH stratospheric sudden warming occurred, supporting the view that healing of the Antarctic ozone layer has now begun to occur during the month of September.
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.
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.
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.
Geng, L, Lee T Murray, L J Mickley, and Pu Lin, et al., June 2017: Isotopic evidence of multiple controls on atmospheric oxidants over climate transitions. Nature, 546(7656), DOI:10.1038/nature22340. Abstract
The abundance of tropospheric oxidants, such as ozone (O3) and hydroxyl (OH) and peroxy radicals (HO2 + RO2), determines the lifetimes of reduced trace gases such as methane and the production of particulate matter important for climate and human health. The response of tropospheric oxidants to climate change is poorly constrained owing to large uncertainties in the degree to which processes that influence oxidants may change with climate1 and owing to a lack of palaeo-records with which to constrain levels of atmospheric oxidants during past climate transitions2. At present, it is thought that temperature-dependent emissions of tropospheric O3 precursors and water vapour abundance determine the climate response of oxidants, resulting in lower tropospheric O3 in cold climates while HOx (= OH + HO2 + RO2) remains relatively buffered3. Here we report observations of oxygen-17 excess of nitrate (a proxy for the relative abundance of atmospheric O3 and HOx) from a Greenland ice core over the most recent glacial–interglacial cycle and for two Dansgaard–Oeschger events. We find that tropospheric oxidants are sensitive to climate change with an increase in the O3/HOx ratio in cold climates, the opposite of current expectations. We hypothesize that the observed increase in O3/HOx in cold climates is driven by enhanced stratosphere-to-troposphere transport of O3, and that reactive halogen chemistry is also enhanced in cold climates. Reactive halogens influence the oxidative capacity of the troposphere directly as oxidants themselves and indirectly4 via their influence on O3 and HOx. The strength of stratosphere-to-troposphere transport is largely controlled by the Brewer–Dobson circulation5, which may be enhanced in colder climates owing to a stronger meridional gradient of sea surface temperatures6, with implications for the response of tropospheric oxidants7 and stratospheric thermal and mass balance8. These two processes may represent important, yet relatively unexplored, climate feedback mechanisms during major climate transitions.
Hardiman, Steven C., and Pu Lin, et al., March 2017: The influence of dynamical variability on the observed Brewer-Dobson circulation trend. Geophysical Research Letters, 44(6), DOI:10.1002/2017GL072706. Abstract
The strength of the Brewer-Dobson circulation (BDC) is predicted to increase due to climate change. However, this increase has yet to be robustly detected in observational analyses. In this study a long control simulation is used to calculate the Time of Emergence of the BDC trend and how much of that trend may be masked by dynamical variability in current observations. A Time of Emergence of around 30 years is found (assuming a 2%/decade trend in the BDC), similar to the length of current reanalysis data sets. However, the discrepancies in vertical velocities between different reanalysis products remain far larger than any predicted trend. Furthermore, dynamical variability can completely mask the BDC trend on time scales shorter than around 12 years. Thus, more robust observational analyses of vertical velocity are likely to be needed for at least the next decade before detection of a statistically significant trend can be expected.
This study explores the role of the stratosphere as a source of seasonal predictability of surface climate over Northern Hemisphere extra-tropics both in the observations and climate model predictions. A suite of numerical experiments, including climate simulations and retrospective forecasts, are set up to isolate the role of the stratosphere in seasonal predictive skill of extra-tropical near surface land temperature. We show that most of the lead-0 month spring predictive skill of land temperature over extra-tropics, particularly over northern Eurasia, stems from stratospheric initialization. We further reveal that this predictive skill of extra-tropical land temperature arises from skillful prediction of the Arctic Oscillation (AO). The dynamical connection between the stratosphere and troposphere is also demonstrated by the significant correlation between the stratospheric polar vortex and sea level pressure anomalies, as well as the migration of the stratospheric zonal wind anomalies to the lower troposphere.
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.
Solomon, S, D Ivy, M Gupta, J Bandoro, B D Santer, Qiang Fu, and Pu Lin, et al., August 2017: Mirrored Changes in Antarctic Ozone and Stratospheric Temperature in the late 20th versus early 21st Centuries. Journal of Geophysical Research: Atmospheres, 122(16), DOI:10.1002/2017JD026719. Abstract
Observed and modeled patterns of lower stratospheric seasonal trends in Antarctic ozone and temperature in the late 20th (1979-2000) and the early 21st (2000-2014) centuries are compared. Patterns of pre-2000 observed Antarctic ozone decreases and stratospheric cooling as a function of month and pressure are followed by opposite-signed (i.e., “mirrored”) patterns of ozone increases and warming post-2000. An interactive chemistry-climate model forced by changes in anthropogenic ozone depleting substances produces broadly similar mirrored features. Statistical analysis of unforced model simulations (from long term model control simulations of a few centuries up to 1000 years) suggests that internal and solar natural variability alone is unable to account for the pattern of observed ozone trend mirroring, implying that forcing is the dominant driver of this behavior. Radiative calculations indicate that ozone increases have contributed to Antarctic warming of the lower stratosphere over 2000-2014, but dynamical changes that are likely due to internal variability over this relatively short period also appear to be important. Overall, the results support the recent finding that the healing of the Antarctic ozone hole is underway, and that coupling between dynamics, chemistry, and radiation is important for a full understanding of the causes of observed stratospheric temperature and ozone changes.
Guan, Xiaodan, Jianping Huang, Ruixia Guo, and Pu Lin, July 2015: The role of dynamically induced variability in the recent warming trend slowdown over the Northern Hemisphere. Scientific Reports, 5, 12669, DOI:10.1038/srep12669. Abstract
Since the slowing of the trend of increasing surface air temperature (SAT) in the late 1990 s, intense interest and debate have arisen concerning the contribution of human activities to the warming observed in previous decades. Although several explanations have been proposed for the warming-trend slowdown (WTS), none has been generally accepted. We investigate the WTS using a recently developed methodology that can successfully identify and separate the dynamically induced and radiatively forced SAT changes from raw SAT data. The dynamically induced SAT changes exhibited an obvious cooling effect relative to the warming effect of the adjusted SAT in the hiatus process. A correlation analysis suggests that the changes are dominated primarily by the North Atlantic Oscillation (NAO), Pacific Decadal Oscillation (PDO) and Atlantic Multidecadal Oscillation (AMO). Our results confirm that dynamically induced variability caused the WTS. The radiatively forced SAT changes are determined mainly by anthropogenic forcing, indicating the warming influence of greenhouse gases (GHGs), which reached levels of 400 ppm during the hiatus period. Therefore, the global SAT will not remain permanently neutral. The increased radiatively forced SAT will be amplified by increased dynamically induced SAT when the natural mode returns to a warming phase in the next period.
Guan, Xiaodan, Jianping Huang, Ruixia Guo, Haipeng Yu, Pu Lin, and Yanting Zhang, December 2015: Role of radiatively forced temperature changes in enhanced semi-arid warming in the cold season over east Asia. Atmospheric Chemistry and Physics, 15, DOI:10.5194/acp-15-13777-201513777-13786. Abstract
As climate change has occurred over east Asia since the 1950s, intense interest and debate have arisen concerning the contribution of human activities to the observed warming in past decades. In this study, we investigate regional surface temperature change during the boreal cold season using a recently developed methodology that can successfully identify and separate the dynamically induced temperature (DIT) and radiatively forced temperature (RFT) changes in raw surface air temperature (SAT) data. For regional averages, DIT and RFT contribute 44 and 56 % to the SAT over east Asia, respectively. The DIT changes dominate the SAT decadal variability and are mainly determined by internal climate variability, represented by the North Atlantic Oscillation (NAO), Pacific Decadal Oscillation (PDO), and Atlantic Multi-decadal Oscillation (AMO). Radiatively forced SAT changes have made a major contribution to the global-scale warming trend and the regional-scale enhanced semi-arid warming (ESAW). Such enhanced warming is also found in radiatively forced daily maximum and minimum SAT. The long-term global-mean SAT warming trend is mainly related to radiative forcing produced by global well-mixed greenhouse gases. The regional anthropogenic radiative forcing, however, caused the enhanced warming in the semi-arid region, which may be closely associated with local human activities. Finally, the relationship between the so-called "global warming hiatus" and regional enhanced warming is discussed.
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.
Fueglistaler, Stephan, Marta Abalos, Thomas J Flannaghan, Pu Lin, and W J Randel, December 2014: Variability and trends in dynamical forcing of tropical lower stratospheric temperatures. Atmospheric Chemistry and Physics, 14(24), DOI:10.5194/acp-14-13439-2014. Abstract
We analyse the relation between tropical lower stratospheric temperatures and dynamical forcing over the period 1980–2011 using NCEP, MERRA and ERA-Interim reanalyses. The tropical mean thermodynamic energy equation with Newtonian cooling for radiation is forced with two dynamical predictors: (i) the average eddy heat flux of both hemispheres; and (ii) tropical upwelling estimated from momentum balance following Randel et al. (2002). The correlation (1995–2011) for deseasonalised tropical average temperatures at 70 hPa with the eddy heat flux based predictor is 0.84 for ERA-Interim (0.77 for the momentum balance calculation), and 0.87 for MERRA. The eddy heat flux based predictor indicates a dynamically forced cooling of the tropics of ∼−0.1 K decade−1 (∼−0.2 K decade−1 excluding volcanic periods) for the period 1980–2011 in MERRA and ERA-Interim. ERA-Interim eddy heat fluxes drift slightly relative to MERRA in the 2000's, possibly due to onset of GPS temperature data assimilation. While NCEP gives a small warming trend, all 3 reanalyses show a similar seasonality, with strongest cooling in January/February (∼−0.4 K decade−1, from northern hemispheric forcing) and October (∼−0.3 K decade−1, from southern hemispheric forcing). Months preceding and following the peaks in cooling trends show pronounced smaller, or even warming, trends. Consequently, the seasonality in the trends arises in part due to a temporal shift in eddy activity. Over all months, the Southern Hemisphere contributes more to the tropical cooling in both MERRA and ERA-Interim. The residual time series (observed minus estimate of dynamically forced temperature) are well correlated between ERA-Interim and MERRA, with differences largely due to temperature differences. The residual time series is dominated by the modification of the radiative balance by volcanic aerosol following the eruption of El Chichon (maximum warming of ∼3 K at 70 hPa) and Pinatubo (maximum warming of ∼4 K at 70 hPa), with a strong dynamical response during Pinatubo partially masking the aerosol heating.