We describe the baseline model configuration and simulation characteristics of the Geophysical Fluid Dynamics Laboratory (GFDL)'s Land Model version 4.1 (LM4.1), which builds on component and coupled model developments over 2013–2019 for the coupled carbon-chemistry-climate Earth System Model Version 4.1 (ESM4.1) simulation as part of the sixth phase of the Coupled Model Intercomparison Project. Analysis of ESM4.1/LM4.1 is focused on biophysical and biogeochemical processes and interactions with climate. Key features include advanced vegetation dynamics and multi-layer canopy energy and moisture exchanges, daily fire, land use representation, and dynamic atmospheric dust coupling. We compare LM4.1 performance in the GFDL Earth System Model (ESM) configuration ESM4.1 to the previous generation component LM3.0 in the ESM2G configuration. ESM4.1/LM4.1 provides significant improvement in the treatment of ecological processes from GFDL's previous generation models. However, ESM4.1/LM4.1 likely overestimates the influence of land use and land cover change on vegetation characteristics, particularly on pasturelands, as it overestimates the competitiveness of grasses versus trees in the tropics, and as a result, underestimates present-day biomass and carbon uptake in comparison to observations.
The rate of sea level rise (SLR) along the Southeast Coast of the U.S. increased significantly after 2010. While anthropogenic radiative forcing causes an acceleration of global mean SLR, regional changes in the rate of SLR are strongly influenced by internal variability. Here we use observations and climate models to show that the rapid increase in the rate of SLR along the U.S. Southeast Coast after 2010 is due in part to multidecadal buoyancy-driven Atlantic meridional overturning circulation (AMOC) variations, along with heat transport convergence from wind-driven ocean circulation changes. We show that an initialized decadal prediction system can provide skillful regional SLR predictions induced by AMOC variations 5 years in advance, while wind-driven sea level variations are predictable 2 years in advance. Our results suggest that the rate of coastal SLR and its associated flooding risk along the U.S. southeastern seaboard are potentially predictable on multiyear timescales.
Using a state-of-the-art coupled general circulation model, physical processes underlying Antarctic sea ice multidecadal variability and predictability are investigated. Our model simulations constrained by atmospheric reanalysis and observed sea surface temperature broadly capture a multidecadal variability in the observed sea ice extent (SIE) with a low sea ice state (late 1970s–1990s) and a high sea ice state (2000s–early 2010s), although the model overestimates the SIE decrease in the Weddell Sea around the 1980s. The low sea ice state is largely due to the deepening of the mixed layer and the associated deep convection that brings subsurface warm water to the surface. During the high sea ice period (post-2000s), the deep convection substantially weakens, so surface wind variability plays a greater role in the SIE variability. Decadal retrospective forecasts started from the above model simulations demonstrate that the Antarctic sea ice multidecadal variability can be skillfully predicted 6–10 years in advance, showing a moderate correlation with the observation. Ensemble members with a deeper mixed layer and stronger deep convection tend to predict a larger sea ice decrease in the 1980s, whereas members with a larger surface wind variability tend to predict a larger sea ice increase after the 2000s. Therefore, skillful simulation and prediction of the Antarctic sea ice multidecadal variability require accurate simulation and prediction of the mixed layer, deep convection, and surface wind variability in the model.
The Model-Analogs technique is used in the present study to assess the decadal sea surface temperature (SST) prediction skill over the Southern Ocean (SO). The Model-Analogs here is based on reanalysis products and model control simulations that have ∼1° ocean/ice (refined to 0.5° at high latitudes) components and 100 km atmosphere/land components. It is found that the model analog hindcasts show comparable skills with the initialized retrospective decadal hindcasts south of 50°S, with even higher skills over the Weddell Sea at longer lead years. The high SST skills primarily arise from the successful capture of SO deep convection states. This deep ocean memory and the associated decadal predictability are also clearly seen when we assess the Model-Analogs technique in a perfect model context. Within 30°S–50°S latitudinal band, the model analog hindcasts show low skills. When we include the externally forced signals estimated from the large ensemble simulations, the model analog hindcasts and initialized decadal hindcasts show identical skills. The Model-Analogs method therefore provides a great baseline for developing future decadal forecast systems. It is unclear whether such analog techniques would also be successful with models that explicitly resolve ocean mesoscale eddies or other small-scale processes. This area of research needs to be explored further.
Long-term sea-level rise and multiyear to decadal sea level variations pose substantial risks for flooding and erosion in coastal communities. Here we use observations and climate model predictions to show that sea level variations along the U.S. East Coast are skillfully predictable 3 to 10 years in advance. The most predictable component of sea level is a basin scale upward trend, predictable a decade in advance and primarily a response to increasing greenhouse gases. Significant additional predictability comes from multidecadal variations of the Atlantic Meridional Overturning Circulation (AMOC). While perfect model simulations show AMOC-related sea level predictability of 5-7 years, model biases and initialization uncertainties reduce the realized predictive skill to 3-5 years, depending on location. Overall, greenhouse gas warming and predictable AMOC variations lead to multiyear to decadal prediction skill for sea level along the U.S. East Coast. Such skill could have significant societal benefit for planning and adaptation.
Research over the past decade has demonstrated that dynamical forecast systems can skillfully predict pan-Arctic sea ice extent (SIE) on the seasonal time scale; however, there have been fewer assessments of prediction skill on user-relevant spatial scales. In this work, we evaluate regional Arctic SIE predictions made with the Forecast-Oriented Low Ocean Resolution (FLOR) and Seamless System for Prediction and Earth System Research (SPEAR_MED) dynamical seasonal forecast systems developed at the NOAA/Geophysical Fluid Dynamics Laboratory. Compared to FLOR, we find that the recently developed SPEAR_MED system displays improved skill in predicting regional detrended SIE anomalies, partially owing to improvements in sea ice concentration (SIC) and thickness (SIT) initial conditions. In both systems, winter SIE is skillfully predicted up to 11 months in advance, whereas summer minimum SIE predictions are limited by the Arctic spring predictability barrier, with typical skill horizons of roughly 4 months. We construct a parsimonious set of simple statistical prediction models to investigate the mechanisms of sea ice predictability in these systems. Three distinct predictability regimes are identified: a summer regime dominated by SIE and SIT anomaly persistence; a winter regime dominated by SIE and upper-ocean heat content (uOHC) anomaly persistence; and a combined regime in the Chukchi Sea, characterized by a trade-off between uOHC-based and SIT-based predictability that occurs as the sea ice edge position evolves seasonally. The combination of regional SIE, SIT, and uOHC predictors is able to reproduce the SIE skill of the dynamical models in nearly all regions, suggesting that these statistical predictors provide a stringent skill benchmark for assessing seasonal sea ice prediction systems.
Hermanson, Leon, Doug Smith, Melissa Seabrook, Roberto Bilbao, Francisco J Doblas-Reyes, Etienne Tourigny, Vladimir Lapin, Viatcheslav Kharin, William J Merryfield, Reinel Sospedra-Alfonso, Panos Athanasiadis, Dario Nicolí, Silvio Gualdi, Nick Dunstone, Rosie Eade, Adam A Scaife, Mark A Collier, Terence O'Kane, Vassili Kitsios, Paul Sandery, Klaus Pankatz, Barbara Früh, Holger Pohlmann, Wolfgang A Müller, Takahito Kataoka, Hiroaki Tatebe, Masayoshi Ishii, Yukiko Imada, Tim Kruschke, Torben Koenigk, Mehdi Pasha Karami, Shuting Yang, Tian Tian, Liping Zhang, Thomas L Delworth, Xiaosong Yang, and Fanrong Zeng, et al., April 2022: WMO global annual to decadal climate update: A prediction for 2021–25. Bulletin of the American Meteorological Society, 103(4), DOI:10.1175/BAMS-D-20-0311.1E1117-E1129. Abstract
As climate change accelerates, societies and climate-sensitive socioeconomic sectors cannot continue to rely on the past as a guide to possible future climate hazards. Operational decadal predictions offer the potential to inform current adaptation and increase resilience by filling the important gap between seasonal forecasts and climate projections. The World Meteorological Organization (WMO) has recognized this and in 2017 established the WMO Lead Centre for Annual to Decadal Climate Predictions (shortened to “Lead Centre” below), which annually provides a large multimodel ensemble of predictions covering the next 5 years. This international collaboration produces a prediction that is more skillful and useful than any single center can achieve. One of the main outputs of the Lead Centre is the Global Annual to Decadal Climate Update (GADCU), a consensus forecast based on these predictions. This update includes maps showing key variables, discussion on forecast skill, and predictions of climate indices such as the global mean near-surface temperature and Atlantic multidecadal variability. it also estimates the probability of the global mean temperature exceeding 1.5°C above preindustrial levels for at least 1 year in the next 5 years, which helps policy-makers understand how closely the world is approaching this goal of the Paris Agreement. This paper, written by the authors of the GADCU, introduces the GADCU, presents its key outputs, and briefly discusses its role in providing vital climate information for society now and in the future.
This study shows that the frequency of North American summertime (June–August) heat extremes is skillfully predicted several months in advance in the newly developed Geophysical Fluid Dynamics Laboratory (GFDL) Seamless System for Prediction and Earth System Research (SPEAR) seasonal forecast system. Using a statistical optimization method, the average predictability time, we identify three large-scale components of the frequency of North American summer heat extremes that are predictable with significant correlation skill. One component, which is related to a secular warming trend, shows a continent-wide increase in the frequency of summer heat extremes and is highly predictable at least 9 months in advance. This trend component is likely a response to external radiative forcing. The second component is largely driven by the sea surface temperatures in the North Pacific and North Atlantic and is significantly correlated with the central U.S. soil moisture. The second component shows largest loadings over the central United States and is significantly predictable 9 months in advance. The third component, which is related to the central Pacific El Niño, displays a dipole structure over North America and is predictable up to 4 months in advance. Potential implications for advancing seasonal predictions of North American summertime heat extremes are discussed.
The Kuroshio Extension (KE), an eastward-flowing jet located in the Pacific western boundary current system, exhibits prominent seasonal-to-decadal variability, which is crucial for understanding climate variations in the northern midlatitudes. We explore the representation and prediction skill for the KE in the GFDL SPEAR (Seamless System for Prediction and Earth System Research) coupled model. Two different approaches are used to generate coupled reanalyses and forecasts: 1) restoring the coupled model’s SST and atmospheric variables toward existing reanalyses, or 2) assimilating SST and subsurface observations into the coupled model without atmospheric assimilation. Both systems use an ocean model with 1° resolution and capture the largest sea surface height (SSH) variability over the KE region. Assimilating subsurface observations appears to be essential to reproduce the narrow front and related oceanic variability of the KE jet in the coupled reanalysis. We demonstrate skillful retrospective predictions of KE SSH variability in monthly (up to 1 year) and annual-mean (up to 5 years) KE forecasts in the seasonal and decadal prediction systems, respectively. The prediction skill varies seasonally, peaking for forecasts initialized in January and verifying in September due to the winter intensification of North Pacific atmospheric forcing. We show that strong large-scale atmospheric anomalies generate deterministic oceanic forcing (i.e., Rossby waves), leading to skillful long-lead KE forecasts. These atmospheric anomalies also drive Ekman convergence and divergence, which forms ocean memory, by sequestering thermal anomalies deep into the winter mixed layer that re-emerge in the subsequent autumn. The SPEAR forecasts capture the recent negative-to-positive transition of the KE phase in 2017, projecting a continued positive phase through 2022.
Quantifying the response of atmospheric rivers (ARs) to radiative forcing is challenging due to uncertainties caused by internal climate variability, differences in shared socioeconomic pathways (SSPs), and methods used in AR detection algorithms. In addition, the requirement of medium-to-high model resolution and ensemble sizes to explicitly simulate ARs and their statistics can be computationally expensive. In this study, we leverage the unique 50-km large ensembles generated by a Geophysical Fluid Dynamics Laboratory next-generation global climate model, Seamless system for Prediction and EArth system Research, to explore the warming response in ARs. Under both moderate and high emissions scenarios, increases in AR-day frequency emerge from the noise of internal variability by 2060. This signal is robust across different SSPs and time-independent detection criteria. We further examine an alternative approach proposed by Thompson et al. (2015), showing that unforced AR variability can be approximated by a first-order autoregressive process. The confidence intervals of the projected response can be analytically derived with a single ensemble member.
A subseasonal-to-seasonal (S2S) prediction system was recently developed using the GFDL Seamless System for Prediction and Earth System Research (SPEAR) global coupled model. Based on 20-yr hindcast results (2000–19), the boreal wintertime (November–April) Madden–Julian oscillation (MJO) prediction skill is revealed to reach 30 days measured before the anomaly correlation coefficient of the real-time multivariate (RMM) index drops to 0.5. However, when the MJO is partitioned into four distinct propagation patterns, the prediction range extends to 38, 31, and 31 days for the fast-propagating, slow-propagating, and jumping MJO patterns, respectively, but falls to 23 days for the standing MJO. A further improvement of MJO prediction requires attention to the standing MJO given its large gap with its potential predictability (38 days). The slow-propagating MJO detours southward when traversing the Maritime Continent (MC), and confronts the MC prediction barrier in the model, while the fast-propagating MJO moves across the central MC without this prediction barrier. The MJO diversity is modulated by stratospheric quasi-biennial oscillation (QBO): the standing (slow-propagating) MJO coincides with significant westerly (easterly) phases of QBO, partially explaining the contrasting MJO prediction skill between these two QBO phases. The SPEAR model shows its capability, beyond the propagation, in predicting their initiation for different types of MJO along with discrete precursory convection anomalies. The SPEAR model skillfully predicts the observed distinct teleconnections over the North Pacific and North America related to the standing, jumping, and fast-propagating MJO, but not the slow-propagating MJO. These findings highlight the complexities and challenges of incorporating MJO prediction into the operational prediction of meteorological variables.
One of the most puzzling observed features of recent climate has been a multidecadal surface cooling trend over the subpolar Southern Ocean (SO). In this study we use large ensembles of simulations with multiple climate models to study the role of the SO meridional overturning circulation (MOC) in these sea surface temperature (SST) trends. We find that multiple competing processes play prominent roles, consistent with multiple mechanisms proposed in the literature for the observed cooling. Early in the simulations (twentieth century and early twenty-first century) internal variability of the MOC can have a large impact, in part due to substantial simulated multidecadal variability of the MOC. Ensemble members with initially strong convection (and related surface warming due to convective mixing of subsurface warmth to the surface) tend to subsequently cool at the surface as convection associated with internal variability weakens. A second process occurs in the late-twentieth and twenty-first centuries, as weakening of oceanic convection associated with global warming and high-latitude freshening can contribute to the surface cooling trend by suppressing convection and associated vertical mixing of subsurface heat. As the simulations progress, the multidecadal SO variability is suppressed due to forced changes in the mean state and increased oceanic stratification. As a third process, the shallower mixed layers can then rapidly warm due to increasing forcing from greenhouse gas warming. Also, during this period the ensemble spread of SO SST trend partly arises from the spread of the wind-driven Deacon cell strength. Thus, different processes could conceivably have led to the observed cooling trend, consistent with the range of possibilities presented in the literature. To better understand the causes of the observed trend, it is important to better understand the characteristics of internal low-frequency variability in the SO and the response of that variability to global warming.
The low Antarctic sea ice extent following its dramatic decline in late 2016 has persisted over a multiyear period. However, it remains unclear to what extent this low sea ice extent can be attributed to changing ocean conditions. Here, we investigate the causes of this period of low Antarctic sea ice extent using a coupled climate model partially constrained by observations. We find that the subsurface Southern Ocean played a smaller role than the atmosphere in the extreme sea ice extent low in 2016, but was critical for the persistence of negative anomalies over 2016–2021. Prior to 2016, the subsurface Southern Ocean warmed in response to enhanced westerly winds. Decadal hindcasts show that subsurface warming has persisted and gradually destabilized the ocean from below, reducing sea ice extent over several years. The simultaneous variations in the atmosphere and ocean after 2016 have further amplified the decline in Antarctic sea ice extent.
Compared to the Arctic, seasonal predictions of Antarctic sea ice have received relatively little attention. In this work, we utilize three coupled dynamical prediction systems developed at the Geophysical Fluid Dynamics Laboratory to assess the seasonal prediction skill and predictability of Antarctic sea ice. These systems, based on the FLOR, SPEAR_LO, and SPEAR_MED dynamical models, differ in their coupled model components, initialization techniques, atmospheric resolution, and model biases. Using suites of retrospective initialized seasonal predictions spanning 1992–2018, we investigate the role of these factors in determining Antarctic sea ice prediction skill and examine the mechanisms of regional sea ice predictability. We find that each system is capable of skillfully predicting regional Antarctic sea ice extent (SIE) with skill that exceeds a persistence forecast. Winter SIE is skillfully predicted 11 months in advance in the Weddell, Amundsen/Bellingshausen, Indian, and west Pacific sectors, whereas winter skill is notably lower in the Ross sector. Zonally advected upper-ocean heat content anomalies are found to provide the crucial source of prediction skill for the winter sea ice edge position. The recently developed SPEAR systems are more skillful than FLOR for summer sea ice predictions, owing to improvements in sea ice concentration and sea ice thickness initialization. Summer Weddell SIE is skillfully predicted up to 9 months in advance in SPEAR_MED, due to the persistence and drift of initialized sea ice thickness anomalies from the previous winter. Overall, these results suggest a promising potential for providing operational Antarctic sea ice predictions on seasonal time scales.
Atmospheric rivers (ARs) exert significant socioeconomic impacts in western North America, where 30% of the annual precipitation is determined by ARs that occur in less than 15% of wintertime. ARs are thus beneficial to water supply but can produce extreme precipitation hazards when making landfall. While most prevailing research has focused on the subseasonal (<5 weeks) prediction of ARs, only limited efforts have been made for AR forecasts on multiseasonal timescales (>3 months) that are crucial for water resource management and disaster preparedness. Through the analysis of reanalysis data and retrospective predictions from a new seasonal-to-decadal forecast system, this research shows the existing potential of multiseasonal AR frequency forecasts with predictive skills 9 months in advance. Additional analysis explores the dominant predictability sources and challenges for multiseasonal AR prediction.
Using GFDL's new coupled model SPEAR, we have developed a decadal coupled reanalysis/initialization system (DCIS) that does not use subsurface ocean observations. In DCIS, the winds and temperature in the atmosphere, along with sea surface temperature (SST), are restored to observations. Under this approach the ocean component of the coupled model experiences a sequence of surface heat and momentum fluxes that are similar to observations. DCIS offers two initialization approaches, called A1 and A2, which differ only in the atmospheric forcing from observations. In A1, the atmospheric winds/temperature are restored toward the JRA reanalysis; in A2, surface pressure observations are assimilated in the model. Two sets of coupled reanalyses have been completed during 1961–2019 using A1 and A2, and they show very similar multi-decadal variations of the Atlantic Meridional Overturning Circulation (AMOC). Two sets of retrospective decadal forecasts were then conducted using initial conditions from the A1 and A2 reanalyses. In comparison with previous prediction system CM2.1, SPEAR-A1/A2 shows comparable skill of predicting the North Atlantic subpolar gyre SST, which is highly correlated with initial values of AMOC at all lead years. SPEAR-A1 significantly outperforms CM2.1 in predicting multi-decadal SST trends in the Southern Ocean (SO). Both A1 and A2 have skillful prediction of Sahel precipitation and the associated ITCZ shift. The prediction skill of SST is generally lower in A2 than A1 especially over SO presumably due to the sparse surface pressure observations.
We document the development and simulation characteristics of the next generation modeling system for seasonal to decadal prediction and projection at the Geophysical Fluid Dynamics Laboratory (GFDL). SPEAR (Seamless System for Prediction and EArth System Research) is built from component models recently developed at GFDL ‐ the AM4 atmosphere model, MOM6 ocean code, LM4 land model and SIS2 sea ice model. The SPEAR models are specifically designed with attributes needed for a prediction model for seasonal to decadal time scales, including the ability to run large ensembles of simulations with available computational resources. For computational speed SPEAR uses a coarse ocean resolution of approximately 1.0o (with tropical refinement). SPEAR can use differing atmospheric horizontal resolutions ranging from 1o to 0.25o. The higher atmospheric resolution facilitates improved simulation of regional climate and extremes. SPEAR is built from the same components as the GFDL CM4 and ESM 4 models, but with design choices geared toward seasonal to multidecadal physical climate prediction and projection. We document simulation characteristics for the time‐mean climate, aspects of internal variability, and the response to both idealized and realistic radiative forcing change. We describe in greater detail one focus of the model development process that was motivated by the importance of the Southern Ocean to the global climate system. We present sensitivity tests that document the influence of the Antarctic surface heat budget on Southern Ocean ventilation and deep global ocean circulation. These findings were also useful in the development processes for the GFDL CM4 and ESM 4 models.
The Caribbean low-level jet (CLLJ) is an important component of the atmospheric circulation over the Intra-Americas Sea (IAS) which impacts the weather and climate both locally and remotely. It influences the rainfall variability in the Caribbean, Central America, northern South America, the tropical Pacific and the continental Unites States through the transport of moisture. We make use of high-resolution coupled and uncoupled models from the Geophysical Fluid Dynamics Laboratory (GFDL) to investigate the simulation of the CLLJ and its teleconnections and further compare with low-resolution models. The high-resolution coupled model FLOR shows improvements in the simulation of the CLLJ and its teleconnections with rainfall and SST over the IAS compared to the low-resolution coupled model CM2.1. The CLLJ is better represented in uncoupled models (AM2.1 and AM2.5) forced with observed sea-surface temperatures (SSTs), emphasizing the role of SSTs in the simulation of the CLLJ. Further, we determine the forecast skill for observed rainfall using both high- and low-resolution predictions of rainfall and SSTs for the July–August–September season. We determine the role of statistical correction of model biases, coupling and horizontal resolution on the forecast skill. Statistical correction dramatically improves area-averaged forecast skill. But the analysis of spatial distribution in skill indicates that the improvement in skill after statistical correction is region dependent. Forecast skill is sensitive to coupling in parts of the Caribbean, Central and northern South America, and it is mostly insensitive over North America. Comparison of forecast skill between high and low-resolution coupled models does not show any dramatic difference. However, uncoupled models show improvement in the area-averaged skill in the high-resolution atmospheric model compared to lower resolution model. Understanding and improving the forecast skill over the IAS has important implications for highly vulnerable nations in the region.
Responses of tropical cyclones (TCs) to CO2 doubling are explored using coupled global climate models (GCMs) with increasingly refined atmospheric/land horizontal grids (~ 200 km, ~ 50 km and ~ 25 km). The three models exhibit similar changes in background climate fields thought to regulate TC activity, such as relative sea surface temperature (SST), potential intensity, and wind shear. However, global TC frequency decreases substantially in the 50 km model, while the 25 km model shows no significant change. The ~ 25 km model also has a substantial and spatially-ubiquitous increase of Category 3–4–5 hurricanes. Idealized perturbation experiments are performed to understand the TC response. Each model’s transient fully-coupled 2 × CO2 TC activity response is largely recovered by “time-slice” experiments using time-invariant SST perturbations added to each model’s own SST climatology. The TC response to SST forcing depends on each model’s background climatological SST biases: removing these biases leads to a global TC intensity increase in the ~ 50 km model, and a global TC frequency increase in the ~ 25 km model, in response to CO2-induced warming patterns and CO2 doubling. Isolated CO2 doubling leads to a significant TC frequency decrease, while isolated uniform SST warming leads to a significant global TC frequency increase; the ~ 25 km model has a greater tendency for frequency increase. Global TC frequency responds to both (1) changes in TC “seeds”, which increase due to warming (more so in the ~ 25 km model) and decrease due to higher CO2 concentrations, and (2) less efficient development of these“seeds” into TCs, largely due to the nonlinear relation between temperature and saturation specific humidity.
Kam, Jonghun, Thomas R Knutson, Fanrong Zeng, and Andrew T Wittenberg, January 2018: CMIP5 Model-Based Assessment of Anthropogenic Influence on Highly Anomalous Arctic Warmth During November-December 2016, [in “Explaining Extreme Events of 2016 from a Climate Perspective”]. Bulletin of the American Meteorological Society, 99(1), DOI:10.1175/BAMS-D-17-0116.1S34-S38.
Knutson, Thomas R., Jonghun Kam, Fanrong Zeng, and Andrew T Wittenberg, January 2018: CMIP5 Model-Based Assessment of Anthropogenic Influence on Record Global Warmth During 2016, [in “Explaining Extreme Events of 2016 from a Climate Perspective”]. Bulletin of the American Meteorological Society, 99(1), DOI:10.1175/BAMS-D-17-0104.1S11-S15.
Precipitation trends for 1901-2010, 1951-2010 and 1981- 2010 over relatively well-observed global land regions are assessed for detectable anthropogenic influences and for consistency with Coupled Model Intercomparison Project 5 (CMIP5) historical simulations. The CMIP5 historical All-Forcing runs are broadly consistent with the observed trend pattern (1901-2010), but with an apparent low trend bias tendency in the simulations. Despite this bias, observed and modeled trends are statistically consistent over 59% of the analyzed area. Over 20% (9%) of the analyzed area, increased (decreased) precipitation is partly attributable to anthropogenic forcing. These inferred human-induced changes include: increases over regions of the north-central U.S., southern Canada, Europe, and southern South America; and decreases over parts of the Mediterranean region and northern tropical Africa. Trends for the shorter periods (1951-2010 and 1981-2010) do not indicate a prominent low trend bias in the models, as found for the 1901-2010 trends. An atmosphere-only model, forced with observed sea surface temperatures and other climate forcing agents, also under-predicts the observed precipitation increase in the northern hemisphere extratropics since 1901. The CMIP5 All-Forcing ensemble’s low bias in simulated trends since 1901 is a tentative finding which, if borne out in further studies, suggests that precipitation projections using these regions/models could overestimate future drought risk, and underestimate future flooding risk.
Unprecedented high intensity flooding induced by extreme precipitation was reported over Chennai in India during November-December of 2015, which led to extensive damage to human life and property. It is of utmost importance to determine the odds of occurrence of such extreme floods in future and the related climate phenomena, for planning and mitigation purposes. Here, we make use of a suite of simulations from GFDL high-resolution coupled climate models to investigate the odds of occurrence of extreme floods induced by extreme precipitation over Chennai and the role of radiative forcing and/or large-scale SST forcing in enhancing the probability of such events in future. Climate of 20th century experiments with large ensembles suggest that the radiative forcing may not enhance the probability of extreme floods over Chennai. Doubling of CO2 experiments also fail to show evidence for increase of such events in a global warming scenario. Further, this study explores the role of SST forcing from the Indian and Pacific Oceans on the odds of occurrence of Chennai-like floods. Neither an El Niño nor La Niña enhances the probability of extreme floods over Chennai. However, warm Bay of Bengal tends to increase the odds of occurrence of extreme Chennai-like floods. The atmospheric condition such as a tropical depression over Bay of Bengal favoring the transport of moisture from warm Bay of Bengal is conducive for intense precipitation.
The Pacific equatorial cold tongue plays a leading role in Earth’s strongest and most predictable climate signals. To illuminate the processes governing cold tongue temperatures, the upper-ocean heat budget is explored using the GFDL FLOR coupled GCM. Starting from the exact temperature budget for layers of time-varying thickness, the layer temperature tendency terms are studied using hourly-, daily-, and monthly-mean output from a 30-year simulation driven by present-day radiative forcings. The budget is then applied to (1) a surface mixed layer whose temperature is highly correlated with SST, in which the air-sea heat flux is balanced mainly by downward diffusion of heat across the layer base; and (2) a thicker advective layer that subsumes most of the vertical mixing, in which the air-sea heat flux is balanced mainly by monthly-scale advection. The surface warming from shortwave fluxes and submonthly meridional advection, and the subsurface cooling from monthly vertical advection, are both shown to be essential to maintain the cold tongue thermal stratification against the destratifying effects of vertical mixing. Although layer undulations strongly mediate the tendency terms on diurnal-to-interannual scales, the 30-year-mean tendencies are found to be well summarized by analogous budgets developed for stationary but spatially-varying layers. The results are used to derive practical simplifications of the exact budget, to support the analyses in Part II of this paper and to facilitate broader application of heat budget analyses when evaluating and comparing climate simulations.
The heat budget of the Pacific equatorial cold tongue (ECT) is explored using the GFDL FLOR coupled GCM and ocean reanalyses, leveraging the two-layer framework developed in Part I. Despite FLOR’s relatively weak meridional stirring by tropical instability waves (TIWs), the model maintains a reasonable SST and thermocline depth in the ECT via two compensating biases: (1) enhanced monthly-scale vertical advective cooling below the surface mixed layer (SML), due to overly cyclonic off-equatorial wind stress which acts to cool the equatorial source waters; and (2) an excessive SST contrast between the ECT and off-equator, which boosts the equatorward heat transport by TIWs. FLOR’s strong advective cooling at the SML base is compensated by strong downward diffusion of heat out of the SML, which then allows FLOR’s ECT to take up a realistic heat flux from the atmosphere. Correcting FLOR’s climatological SST and wind stress biases via flux adjustment (FA) leads to weaker deep advective cooling of the ECT, which then erodes the upper-ocean thermal stratification, enhances vertical mixing, and excessively deepens the thermocline. FA does strengthen FLOR’s meridional shear of the zonal currents in the east Pacific, but this does not amplify the simulated TIWs nor their equatorward heat transport, likely due to FLOR’s coarse zonal ocean resolution. The analysis suggests that to advance coupled simulations of the ECT, improved winds and surface heat fluxes must go hand in hand with improved subseasonal and parameterized ocean processes. Implications for model development and the tropical Pacific observing system are discussed.
The Geophysical Fluid Dynamics Laboratory (GFDL) has recently developed two global coupled GCMs, FLOR and HiFLOR, which are now being utilized for climate research and seasonal predictions. Compared to their predecessor CM2.1, the new versions have improved ocean/atmosphere physics and numerics, and refinement of the atmospheric horizontal grid from 220 km (CM2.1) to 55 km (FLOR) and 26 km (HiFLOR). Both FLOR and HiFLOR demonstrate greatly improved simulations of the tropical Pacific annual‐mean climatology, with FLOR practically eliminating any equatorial cold bias in sea surface temperature. An additional model experiment (LOAR1) using FLOR's ocean/atmosphere physics, but with the atmospheric grid coarsened toward that of CM2.1, is used to further isolate the impacts of the refined atmospheric grid versus the improved physics and numerics. The improved ocean/atmosphere formulations are found to produce more realistic tropical Pacific patterns of sea surface temperature and rainfall, surface heat fluxes, ocean mixed layer depths, surface currents, and tropical instability wave (TIW) activity; enhance the near‐surface equatorial upwelling; and reduce the inter‐centennial warm drift of the tropical Pacific upper ocean. The atmospheric grid refinement further improves these features, and also improves the tropical Pacific surface wind stress, implied Ekman and Sverdrup transports, subsurface temperature and salinity structure, and heat advection in the equatorial upper ocean. The results highlight the importance of nonlocal air‐sea interactions in the tropical Pacific climate system, including the influence of off‐equatorial surface fluxes on the equatorial annual‐mean state. Implications are discussed for improving future simulations, observations, and predictions of tropical Pacific climate.
A “typical” El Niño leads to wet (dry) wintertime anomalies over the southern (northern) half of the Western United States (WUS). However, during the strong El Niño of 2015/16, the WUS winter precipitation pattern was roughly opposite to this canonical (average of the record) anomaly pattern. To understand why this happened, and whether it was predictable, we use a suite of high-resolution seasonal prediction experiments with coupled climate models. We find that the unusual 2015/16 precipitation pattern was predictable at zero-lead time horizon when the ocean/atmosphere/land components were initialized with observations. However, when the ocean alone is initialized the coupled model fails to predict the 2015/16 pattern, although ocean initial conditions alone can reproduce the observed WUS precipitation during the 1997/98 strong El Niño. Further observational analysis shows that the amplitudes of the El Niño induced tropical circulation anomalies during 2015/16 were weakened by about 50% relative to those of 1997/98. This was caused by relative cold (warm) anomalies in the eastern (western) tropical Pacific suppressing (enhancing) deep convection anomalies in the eastern (western) tropical Pacific during 2015/16. The reduced El Niño teleconnection led to a weakening of the subtropical westerly jet over the southeast North Pacific and southern WUS, resulting in the unusual 2015/16 winter precipitation pattern over the WUS. This study highlights the importance of initial conditions not only in the ocean, but in the land and atmosphere as well, for predicting the unusual El Niño teleconnection and its influence on the winter WUS precipitation anomalies during 2015/16.
The relationship between the North Atlantic Oscillation (NAO) and Atlantic sea surface temperature (SST) variability is investigated using models and observations. Coupled climate models are used in which the ocean component is either a fully dynamic ocean, or a slab ocean with no resolved ocean heat transport. On time scales less than ten years NAO variations drive a tripole pattern of SST anomalies in both observations and models. This SST pattern is a direct response of the ocean mixed layer to turbulent surface heat flux anomalies associated with the NAO.
On time scales longer than ten years a similar relationship exists between the NAO and the tripole pattern of SST anomalies in models with a slab ocean. A different relationship exists both for the observations and for models with a dynamic ocean. In these models a positive (negative) NAO anomaly leads, after a decadal-scale lag, to a monopole pattern of warming (cooling) that resembles the Atlantic Multidecadal Oscillation (AMO), although with smaller than observed amplitudes of tropical SST anomalies. Ocean dynamics are critical to this decadal scale response in the models. The simulated Atlantic Meridional Overturning Circulation (AMOC) strengthens (weakens) in response to a prolonged positive (negative) phase of the NAO, thereby enhancing (decreasing) poleward heat transport, leading to broad scale warming (cooling).
We use additional simulations in which heat flux anomalies derived from observed NAO variations from 1901 to 2014 are applied to the ocean component of coupled models. We show that ocean dynamics allow models to reproduce important aspects of the observed AMO, mainly in the subpolar gyre.
The 2015 hurricane season in the Eastern and Central Pacific Oceans (EPO and CPO), particularly around Hawaii, was extremely active – including a record number of tropical cyclones (TCs) and the first instance of three simultaneous Category 4 hurricanes in the EPO and CPO. A strong El Niño developed during the 2015 boreal summer season, and was attributed by some to be the cause of the extreme number of TCs. However, according to a suite of targeted high-resolution model experiments, the extreme 2015 EPO and CPO hurricane season was not primarily induced by the 2015 El Niño’s tropical Pacific warming, but by warming in the subtropical Pacific Ocean. This warming is not typical of El Niño, but rather the “Pacific Meridional Mode (PMM)” superimposed on long-term anthropogenic warming. Although the likelihood of such an extreme year depends on the phase of natural variability, the coupled GCM projects an increase in the frequency of such extremely active TC years over the next few decades for the EPO, CPO, and Hawaii due enhanced subtropical Pacific warming from anthropogenic greenhouse forcing.
Observations and model simulations of the climate responses to strong explosive low-latitude volcanic eruptions suggest a significant increase in the likelihood of El Niño during the eruption and posteruption years, though model results have been inconclusive and have varied in magnitude and even sign. In this study, we test how this spread of responses depends on the initial phase of El Niño-Southern Oscillation (ENSO) in the eruption year and on the eruption's seasonal timing. We employ the Geophysical Fluid Dynamics Laboratory CM2.1 global coupled general circulation model to investigate the impact of the Pinatubo 1991 eruption, assuming that in 1991 ENSO would otherwise be in central or eastern Pacific El Niño, La Niña, or neutral phases. We obtain statistically significant El Niño responses in a year after the eruption for all cases except La Niña, which shows no response in the eastern equatorial Pacific. The eruption has a weaker impact on eastern Pacific El Niños than on central Pacific El Niños. We find that the ocean dynamical thermostat and (to a lesser extent) wind changes due to land-ocean temperature gradients are the main feedbacks affecting El Niño development after the eruption. The El Niño responses to eruptions occurring in summer are more pronounced than for winter and spring eruptions. That the climate response depends on eruption season and initial ENSO phase may help to reconcile apparent inconsistencies among previous studies.
This study examines the year-to-year modulation of the western North Pacific (WNP) tropical cyclones (TC) activity by the Atlantic Meridional Mode (AMM) using both observations and the Geophysical Fluid Dynamics Laboratory Forecast-oriented Low Ocean Resolution Version of CM2.5 (FLOR) global coupled model. 1. The positive (negative) AMM phase suppresses (enhances) WNP TC activity in observations. The anomalous occurrence of WNP TCs results mainly from changes in TC genesis in the southeastern part of the WNP. 2. The observed responses of WNP TC activity to the AMM are connected to the anomalous zonal vertical wind shear (ZVWS) caused by AMM-induced changes to the Walker circulation. During the positive AMM phase, the warming in the North Atlantic induces strong descending flow in the tropical eastern and central Pacific, which intensifies the Walker cell in the WNP. The intensified Walker cell is responsible for the suppressed (enhanced) TC genesis in the eastern (western) part of the WNP by strengthening (weakening) ZVWS. 3. The observed WNPTC–AMM linkage is examined by the long-term control and idealized perturbations experiment with FLOR-FA. A suite of sensitivity experiments strongly corroborate the observed WNPTC–AMM linkage and underlying physical mechanisms.
The impact of multidecadal variations of the Atlantic meridional overturning circulation (AMOC) on the Southern Ocean (SO) is investigated in the current paper using a coupled ocean–atmosphere model. We find that the AMOC can influence the SO via fast atmosphere teleconnections and subsequent ocean adjustments. A stronger than normal AMOC induces an anomalous warm SST over the North Atlantic, which leads to a warming of the Northern Hemisphere troposphere extending into the tropics. This induces an increased equator-to-pole temperature gradient in the Southern Hemisphere (SH) upper troposphere and lower stratosphere due to an amplified tropical upper tropospheric warming as a result of increased latent heat release. This altered gradients leads to a poleward displacement of the SH westerly jet. The wind change over the SO then cools the SST at high latitudes by anomalous northward Ekman transports. The wind change also weakens the Antarctic bottom water (AABW) cell through changes in surface heat flux forcing. The poleward shifted westerly wind decreases the long term mean easterly winds over the Weddell Sea, thereby reducing the turbulent heat flux loss, decreasing surface density and therefore leading to a weakening of the AABW cell. The weakened AABW cell produces a temperature dipole in the SO, with a warm anomaly in the subsurface and a cold anomaly in the surface that corresponds to an increase of Antarctic sea ice. Opposite conditions occur for a weaker than normal AMOC. Our study here suggests that efforts to attribute the recent observed SO variability to various factors should take into consideration not only local process but also remote forcing from the North Atlantic.
This study explores the potential predictability of the Southern Ocean (SO) climate on decadal timescales as represented in the GFDL CM2.1 model using prognostic methods. We conduct perfect model predictability experiments starting from ten different initial states, and show potentially predictable variations of Antarctic bottom water formation (AABW) rates on time scales as long as twenty years. The associated Weddell Sea (WS) subsurface temperatures and Antarctic sea ice have comparable potential predictability as the AABW cell. The predictability of sea surface temperature (SST) variations over the WS and the SO is somewhat smaller, with predictable scales out to a decade. This reduced predictability is likely associated with stronger damping from air-sea interaction. As a complement to our perfect predictability study, we also make hindcasts of SO decadal variability using the GFDL CM2.1 decadal prediction system. Significant predictive skill for SO SST on multi-year time scales is found in the hindcast system. The success of the hindcasts, especially in reproducing observed surface cooling trends, is largely due to initializing the state of the AABW cell. A weak state of the AABW cell leads to cooler surface conditions and more extensive sea ice. Although there are considerable uncertainties regarding the observational data used to initialize the hindcasts, the consistency between the perfect model experiments and the decadal hindcasts at least gives us some indication as to where and to what extent skillful decadal SO forecasts might be possible.
The impact of the North Atlantic Oscillation (NAO) on the Atlantic Meridional Overturning Circulation (AMOC) and large-scale climate is assessed using simulations with three different climate models. Perturbation experiments are conducted in which a pattern of anomalous heat flux corresponding to the NAO is added to the model ocean. Differences between the perturbation experiments and a control illustrate how the model ocean and climate system respond to the NAO. A positive phase of the NAO strengthens the AMOC by extracting heat from the subpolar gyre, thereby increasing deepwater formation, horizontal density gradients, and the AMOC. The flux forcings have the spatial structure of the observed NAO, but the amplitude of the forcing varies in time with distinct periods varying from 2 to 100 years. The response of the AMOC to NAO variations is small at short time scales, but increases up to the dominant time scale of internal AMOC variability (20-30 years for the models used). The amplitude of the AMOC response, and associated oceanic heat transport, is approximately constant as the timescale of the forcing is increased further. In contrast, the response of other properties, such as hemispheric temperature or Arctic sea ice, continues to increase as the time scale of the forcing becomes progressively longer. The larger response is associated with the time integral of the anomalous oceanic heat transport at longer time scales, combined with an increased impact of radiative feedback processes. We show that NAO fluctuations, similar in amplitude to those observed over the last century, can modulate hemispheric temperature by several tenths of a degree.
Pronounced climate changes have occurred since the 1970s, including rapid loss of Arctic sea ice1, large-scale warming2 and increased tropical storm activity3 in the Atlantic. Anthropogenic radiative forcing is likely to have played a major role in these changes4, but the relative influence of anthropogenic forcing and natural variability is not well established. The above changes have also occurred during a period in which the North Atlantic Oscillation has shown marked multidecadal variations5. Here we investigate the role of the North Atlantic Oscillation in these rapid changes through its influence on the Atlantic meridional overturning circulation and ocean heat transport. We use climate models to show that observed multidecadal variations of the North Atlantic Oscillation can induce multidecadal variations in the Atlantic meridional overturning circulation and poleward ocean heat transport in the Atlantic, extending to the Arctic. Our results suggest that these variations have contributed to the rapid loss of Arctic sea ice, Northern Hemisphere warming, and changing Atlantic tropical storm activity, especially in the late 1990s and early 2000s. These multidecadal variations are superimposed on long-term anthropogenic forcing trends that are the dominant factor in long-term Arctic sea ice loss and hemispheric warming.
This study investigates the roles of radiative forcing, sea surface temperatures (SSTs), and atmospheric and land initial conditions in the summer warming episodes of the United States. The summer warming episodes are defined as the significantly above normal (1983-2012) June-August 2-m temperature anomalies, and are referred to as heat waves in this study. Two contrasting cases, the summers of 2006 and 2012, are explored in detail to illustrate the distinct roles of SSTs, direct radiative forcing, and atmospheric and land initial conditions in driving U.S. summer heat waves. For 2012, simulations with the GFDL atmospheric general circulation model reveal that SSTs play a critical role. Further sensitivity experiments reveal the contributions of uniform global SST warming, SSTs in individual ocean basins and direct radiative forcing to the geographic distribution and magnitudes of warm temperature anomalies. In contrast, for 2006, the atmospheric and land initial conditions are key drivers. The atmospheric (land) initial conditions play a major (minor) role in the central and northwestern (eastern) U.S.. Due to changes in radiative forcing, the probability of areal-averaged summer temperature anomalies over U.S. exceeding the observed 2012 anomaly increases with time over the early 21st century. La Niña (El Niño) events tend to increase (reduce) the occurrence rate of heat waves. The temperatures over the central U.S. are mostly influenced by El Niño/La Niña, with the central tropical Pacific playing a more important role than the eastern tropical Pacific. Thus, atmospheric and land initial conditions, SSTs and radiative forcing are all important drivers of, and sources of predictability for U.S. summer heat waves.
Kam, Jonghun, Thomas R Knutson, Fanrong Zeng, and Andrew T Wittenberg, December 2016: Multimodel Assessment of Anthropogenic Influence on Record Global and Regional Warmth During 2015, Section 2 of “[Explaining extreme events of 2015 from a climate perspective]”. Bulletin of the American Meteorological Society, 97(12), DOI:10.1175/BAMS-D-16-0138.1S4-S8.
Tropical cyclone (TC) activity in the North Pacific and North Atlantic Oceans is known to be affected by the El Niño Southern Oscillation (ENSO). This study uses GFDL FLOR model, which has relatively high-resolution in the atmosphere, as a tool to investigate the sensitivity of TC activity to the strength of ENSO events. We show that TCs exhibit a non-linear response to the strength of ENSO in the tropical eastern North Pacific (ENP) but a quasi-linear response in the tropical western North Pacific (WNP) and tropical North Atlantic. Specifically, stronger El Niño results in disproportionate inhibition of TCs in the ENP and North Atlantic, and leads to an eastward shift in the location of TCs in the southeast of the WNP. However, the character of the response of TCs in the Pacific is insensitive to the amplitude of La Niña events. The eastward shift of TCs in the southeast of the WNP in response to a strong El Niño is due to an eastward shift of the convection and of the associated environmental conditions favorable for TCs. The inhibition of TC activity in the ENP and Atlantic during El Niño is attributed to the increase in the number of days with strong vertical wind shear during stronger El Niño events. These results are further substantiated with coupled model experiments. Understanding of the impact of strong ENSO on TC activity is important for present and future climate as the frequency of occurrence of extreme ENSO events is projected to increase in future.
Precipitation extremes have a widespread impact on societies and ecosystems; it is therefore important to understand current and future patterns of extreme precipitation. Here, a set of new global coupled climate models with varying atmospheric resolution has been used to investigate the ability of these models to reproduce observed patterns of precipitation extremes and to investigate changes in these extremes in response to increased atmospheric CO2 concentrations. The atmospheric resolution was increased from 2°×2° grid cells (typical resolution in the CMIP5 archive) to 0.25°×.25° (tropical cyclone-permitting). Analysis has been confined to the contiguous United States (CONUS). It is shown that, for these models, integrating at higher atmospheric resolution improves all aspects of simulated extreme precipitation: spatial patterns, intensities and seasonal timing. In response to 2×CO2 concentrations, all models show a mean intensification of precipitation rates during extreme events of approximately 3-4% K−1. However, projected regional patterns of changes in extremes are dependent on model resolution. For example, the highest-resolution models show increased precipitation rates during extreme events in the hurricane season in the CONUS southeast, this increase is not found in the low-resolution model. These results emphasize that, for the study of extreme precipitation there is a minimum model resolution that is needed to capture the weather phenomena generating the extremes. Finally, the observed record and historical model experiments were used to investigate changes in the recent past. In part because of large intrinsic variability, no evidence was found for changes in extreme precipitation attributable to climate change in the available observed record.
This study aims to assess the connections between the El Niño Southern Oscillation (ENSO) and tropical cyclones near Guam (GuamTC) using the state-of-the-art Geophysical Fluid Dynamics Laboratory (GFDL) Forecast-oriented Low Ocean Resolution Version of CM2.5 (FLOR). In observations, more (less) GuamTCs occur in El Niño (La Niña) years and the ENSO-GuamTC connections arise from TC genesis locations in ENSO phases. The observed ENSO-GuamTC connections are realistically simulated in the two control experiments that use two versions of FLOR, the standard version and another with flux adjustments (FLOR-FA). The ENSO-GuamTC connections in FLOR-FA are closer to observations than those in FLOR because of a better representation of TC genesis during ENSO phases. The physical mechanisms underlying the observed ENSO-GuamTC connections are further supported in the long-term control experiments with FLOR/FLOR-FA. The ENSO-GuamTC connections in sea surface temperature (SST)- and sea surface salinity (SSS)-restoring experiments with FLOR 1990 strongly resemble the observations, suggesting the ENSO-GuamTC connections arise substantially from the forcing of SST. The prediction skill of FLOR-FA for GuamTC frequency is quite promising in terms of correlation and root mean square error and is higher than that of FLOR for the period 1980-2014. This study shows the capability of global climate models (FLOR/FLOR-FA) in simulating the linkage between ENSO and TC activity near a highly localized region (i.e., Guam) and in predicting the frequency of TCs at the sub-basin scale.
Observed austral summertime (November through April) rainfall in southeastern South America (SESA)—including northern Argentina, Uruguay, southern Brazil and Paraguay—has exhibited substantial low-frequency variations with a multi-decadal moistening trend during the 20th century and a subsequent decadal drying trend during the current century. Understanding the mechanisms responsible for these variations is essential for predicting long-term rainfall changes. Here with a suite of attribution experiments using a pair of high-resolution global climate models—GFDL CM2.5 and FLOR_FA, we investigate the causes of these regional rainfall variations. Both models reproduce the 20th-century moistening trend, albeit with a weaker magnitude than observed, in response to the radiative forcing associated with increasing greenhouse gases. The increasing greenhouse gases drive tropical expansion; consequently, the subtropical dry branch of Hadley cell moves away from SESA, leading to the rainfall increase. The amplitude discrepancy between the observed and simulated rainfall changes suggests a possible underestimation by the models of the atmospheric response to the radiative forcing, as well as an important role for low-frequency internal variability in the observed moistening trend. Over the current century, increasing greenhouse gases drive a continuous SESA rainfall increase in the models. However, the observed decadal rainfall decline is largely (~60%) reproduced in response to the observed Pacific trade wind strengthening, which is likely associated with natural Pacific decadal variability. These results suggest that the recent summertime rainfall decline in SESA is temporary and the positive trend will resume in response to both increasing greenhouse gases and a return of Pacific trade winds to normal conditions.
Portions of western North America have experienced prolonged drought over the last decade. This drought has occurred at the same time as the global warming hiatus – a decadal period with little increase in global mean surface temperature. We use climate models and observational analyses to clarify the dual role of recent tropical Pacific changes in driving both the global warming hiatus and North American drought. When we insert observed tropical Pacific wind stress anomalies into coupled models, the simulations produce persistent negative sea surface temperature anomalies in the eastern tropical Pacific, a hiatus in global warming, and drought over North America driven by SST-induced atmospheric circulation anomalies. In our simulations the tropical wind anomalies account for 92% of the simulated North American drought during the recent decade, with 8% from anthropogenic radiative forcing changes. This suggests that anthropogenic radiative forcing is not the dominant driver of the current drought, unless the wind changes themselves are driven by anthropogenic radiative forcing. The anomalous tropical winds could also originate from coupled interactions in the tropical Pacific or from forcing outside the tropical Pacific. The model experiments suggest that if the tropical winds were to return to climatological conditions, then the recent tendency toward North American drought would diminish. Alternatively, if the tropical winds were to persist, then the impact on North American drought would continue; however, the impact of the enhanced Pacific easterlies on global temperature diminishes after a decade or two due to a surface reemergence of warmer water that was initially subducted into the ocean interior.
This study demonstrates skillful seasonal prediction of 2m air temperature and precipitation over land in a new high-resolution climate model developed by Geophysical Fluid Dynamics Laboratory, and explores the possible sources of the skill. We employ a statistical optimization approach to identify the most predictable components of seasonal mean temperature and precipitation over land, and demonstrate the predictive skill of these components. First, we show improved skill of the high-resolution model over the previous lower-resolution model in seasonal prediction of NINO3.4 index and other aspects of interest. Then we measure the skill of temperature and precipitation in the high-resolution model for boreal winter and summer, and diagnose the sources of the skill. Lastly, we reconstruct predictions using a few most predictable components to yield more skillful predictions than the raw model predictions. Over three decades of hindcasts, we find that the two most predictable components of temperature are characterized by a component that is likely due to changes in external radiative forcing in boreal winter and summer, and an ENSO-related pattern in boreal winter. The most predictable components of precipitation in both seasons are very likely ENSO-related. These components of temperature and precipitation can be predicted with significant correlation skill at least 9 months in advance. The reconstructed predictions using only the first few predictable components from the model show considerably better skill relative to observations than raw model predictions. This study shows that the use of refined statistical analysis and a high-resolution dynamical model leads to significant skill in seasonal predictions of 2m air temperature and precipitation over land.
Kam, Jonghun, Thomas R Knutson, Fanrong Zeng, and Andrew T Wittenberg, December 2015: Record annual mean warmth over Europe, the northeast Pacific, and the northwest Atlantic during 2014: Assessment of anthropogenic influence. Bulletin of the American Meteorological Society, 96(12), DOI:10.1175/BAMS-EEE_2014_ch13.1.
This study investigates the seasonality of the relationship between the Great Plains low-level jet (GPLLJ) and the Pacific Ocean from spring to summer, using observational analysis and coupled model experiments. The observed GPLLJ and El Niño-Southern Oscillation (ENSO) relation undergoes seasonal changes with a stronger GPLLJ associated with La Niña in boreal spring and El Niño in boreal summer. The ability of the GFDL FLOR global coupled climate model, which has the high-resolution atmospheric and land components, to simulate the observed seasonality in the GPLLJ-ENSO relationship is assessed. The importance of simulating the magnitude and phase-locking of ENSO accurately in order to better simulate its seasonal teleconnections with the Intra-Americas Seas (IAS) is demonstrated. This study explores the mechanisms for seasonal changes in the GPLLJ-ENSO relation in model and observations. It is hypothesized that ENSO affects the GPLLJ variability through the Caribbean low-level jet (CLLJ) during the summer and spring seasons. These results suggest that climate models with improved ENSO variability would advance our ability to simulate and predict seasonal variations of the GPLLJ and their associated impacts on the United States.
A new high-resolution Geophysical Fluid Dynamics Laboratory (GFDL) coupled model (HiFLOR) has been developed and used to investigate potential skill in simulation and prediction of tropical cyclone (TC) activity. HiFLOR comprises of high-resolution (~25-km mesh) atmosphere and land components and a more moderate-resolution (~100-km mesh) sea ice and ocean components. HiFLOR was developed from the Forecast Oriented Low Resolution Ocean model (FLOR) by decreasing the horizontal grid spacing of the atmospheric component from 50-km to 25-km, while leaving most of the sub-gridscale physical parameterizations unchanged. Compared with FLOR, HiFLOR yields a more realistic simulation of the structure, global distribution, and seasonal and interannual variations of TCs, and a comparable simulation of storm-induced cold wakes and TC-genesis modulation induced by the Madden Julian Oscillation (MJO). Moreover, HiFLOR is able to simulate and predict extremely intense TCs (categories 4 and 5) and their interannual variations, which represents the first time a global coupled model has been able to simulate such extremely intense TCs in a multi-century simulation, sea surface temperature restoring simulations, and retrospective seasonal predictions.
The seasonal predictability of extratropical storm tracks in Geophysical Fluid Dynamics Laboratory (GFDL)’s high-resolution climate model has been investigated using an average predictability time analysis. The leading predictable components of extratropical storm tracks are ENSO-related spatial pattern for both boreal winter and summer, and the second predictable components are mostly due to changes in external radiative forcing and multidecadal oceanic variability. These two predictable components for both seasons show significant correlation skill for all leads from 0 to 9 months, while the skill of predicting the boreal winter storm track is consistently higher than that of the austral winter. The predictable components of extratropical storm tracks are dynamically consistent with the predictable components of the upper troposphere jet flow for both seasons. Over the region with strong storm track signals in North America, the model is able to predict the changes in statistics of extremes connected to storm track changes (e.g., extreme low and high sea level pressure and extreme 2m air temperature) in response to different ENSO phases. These results point towards the possibility of providing skillful seasonal predictions of the statistics of extratropical extremes over land using high-resolution coupled models.
This study examines two sets of high-resolution coupled model forecasts starting from no-tropical cyclone (TC) and correct-TC-statistics initial conditions to understand the role of TC events on climate prediction. While the model with no-TC initial conditions can quickly spin up TCs within a week, the initial conditions with a corrected TC distribution can produce more accurate forecast of sea surface temperature up to one and half months and maintain larger ocean heat content up to 6 months due to enhanced mixing from continuous interactions between initialized and forecasted TCs and the evolving ocean states. The TC-enhanced tropical ocean mixing strengthens the meridional heat transport in the Southern Hemisphere driven primarily by Southern Ocean surface Ekman fluxes but weakens the Northern Hemisphere poleward transport in this model. This study suggests a future plausible initialization procedure for seamless weather-climate prediction when individual convection-permitting cyclone initialization is incorporated into this TC-statistics-permitting framework.
Precipitation in austral autumn and winter has declined over parts of southern and especially southwestern Australia in the past few decades. According to observations and climate models, at least part of this decline is associated with changes in large-scale atmospheric circulation, including a poleward movement of the westerly winds and increasing atmospheric surface pressure over parts of southern Australia. Here we use a high-resolution global climate model to analyse the causes of this rainfall decline. In our simulations, many aspects of the observed regional rainfall decline over southern and southwest Australia are reproduced in response to anthropogenic changes in levels of greenhouse gases and ozone in the atmosphere, whereas anthropogenic aerosols do not contribute to the simulated precipitation decline. Simulations of future climate with this model suggest amplified winter drying over most parts of southern Australia in the coming decades in response to a high-end scenario of changes in radiative forcing. The drying is most pronounced over southwest Australia, with total reductions in austral autumn and winter precipitation of approximately 40% by the late twenty-first century.
Global tropical cyclone (TC) activity is simulated by the Geophysical Fluid Dynamics Laboratory (GFDL) CM2.5, which is a fully coupled global climate model with horizontal resolution of about 50km for atmosphere and 25 km for ocean. The present climate simulation shows fairly realistic global TC frequency, seasonal cycle, and geographical distribution. The model has some notable biases in regional TC activity, including simulating too few TCs in the North Atlantic. The regional biases in TC activity are associated with simulation biases in the large-scale environment such as sea surface temperature, vertical wind shear, and vertical velocity. Despite these biases, the model simulates the large-scale variations of TC activity induced by El Nino/Southern Oscillation fairly realistically.
The response of TC activity in the model to global warming is investigated by comparing the present climate with a CO2 doubling experiment. Globally, TC frequency decreases (-19%) while the intensity increases (+2.7%) in response to CO2 doubling, consistent with previous studies. The average TC lifetime decreases by -4.6%, while the TC size and rainfall increase by about 3% and 12%, respectively. These changes are generally reproduced across the different basins in terms of the sign of the change, although the percent changes vary from basin to basin and within individual basins. For the Atlantic basin, although there is an overall reduction in frequency from CO2 doubling, the warmed climate exhibits increased interannual hurricane frequency variability so that the simulated Atlantic TC activity is enhanced more during unusually warm years in the CO2-warmed climate relative to that in unusually warm years in the control climate.
In this extended abstract, we report on progress in two areas of research at GFDL relating to Indian Ocean regional climate and climate change. The first topic is an assessment of regional surface temperature trends in the Indian Ocean and surrounding region. Here we illustrate the use of a multi-model approach (CMIP3 or CMIP5 model ensembles) to assess whether an anthropogenic warming signal has emerged in the historical data, including identification of where the observed trends are consistent or not with current climate models. Trends that are consistent with All Forcing runs but inconsistent with Natural Forcing Only runs are ones which we can attribute, at least in part, to anthropogenic forcing.
Knutson, Thomas R., Fanrong Zeng, and Andrew T Wittenberg, September 2014: Seasonal and Annual Mean Precipitation Extremes Occurring During 2013: A U.S. Focused Analysis [in "Explaining Extremes of 2013 from a Climate Perspective"]. Bulletin of the American Meteorological Society, 95(9), S19-S23. Abstract
Explaining Extreme Events of 2013 from a Climate Perspective. Bull. Amer. Meteor. Soc., 95, S1–S104.
doi: http://dx.doi.org/10.1175/1520-0477-95.9.S1.1
Knutson, Thomas R., Fanrong Zeng, and Andrew T Wittenberg, September 2014: Multimodel Assessment of Extreme Annual-Mean Warm Anomalies During 2013 over Regions of Australia and the Western Tropical Pacific [in "Explaining Extremes of 2013 from a Climate Perspective"]. Bulletin of the American Meteorological Society, 95(9), S26-S30. Abstract
Explaining Extreme Events of 2013 from a Climate Perspective. Bull. Amer. Meteor. Soc., 95, S1–S104.
doi: http://dx.doi.org/10.1175/1520-0477-95.9.S1.1
Coupled climate model simulations of volcanic eruptions and abrupt changes in CO2 concentration are compared in multiple realizations of Geophysical Fluid Dynamics Laboratory’s (GFDL) CM2.1. The change in global-mean surface temperature (GMST) is analyzed to determine whether a fast component of the climate sensitivity of relevance to the transient climate response (TCR, defined with the 1% yr−1 CO2-increase scenario) can be estimated from shorter-timescale climate changes. The fast component of the climate sensitivity estimated from the response of the climate model to volcanic forcing is similar to that of the simulations forced by abrupt CO2 changes, but is 5–15% smaller than the TCR. In addition, the partition between the top-of-atmosphere radiative restoring and ocean heat uptake is similar across radiative forcing agents. The possible asymmetry between warming and cooling climate perturbations, which may affect the utility of volcanic eruptions for estimating the TCR, is assessed by comparing simulations of abrupt CO2 doubling to abrupt CO2 halving. There is slightly less (~5%) GMST change in 0.5 × CO2 simulations than in 2 × CO2 simulations on the short (~10-yr) timescales relevant to the fast component of the volcanic signal. However, inferring the TCR from volcanic eruptions is more sensitive to uncertainties from internal climate variability and the estimation procedure.
The response of the GMST to volcanic eruptions is similar in GFDL’s CM2.1 and CM3, even though the latter has a higher TCR associated with a multidecadal timescale in its response. This is consistent with the expectation that the fast component of the climate sensitivity inferred from volcanic eruptions is a lower bound for the TCR.
Decadal prediction experiments were conducted as part of CMIP5 using the GFDL-CM2.1 forecast system. The abrupt warming of the North Atlantic subpolar gyre (SPG) that was observed in the mid 1990s is considered as a case study to evaluate our forecast capabilities and better understand the reasons for the observed changes. Initializing the CM2.1 coupled system produces high skill in retrospectively predicting the mid-90s shift, which is not captured by the uninitialized forecasts. All the hindcasts initialized in the early 90s show a warming of the SPG, however, only the ensemble mean hindcasts initialized in 1995 and 1996 are able to reproduce the observed abrupt warming and the associated decrease and contraction of the SPG. Examination of the physical mechanisms responsible for the successful retrospective predictions indicates that initializing the ocean is key to predict the mid 90s warming. The successful initialized forecasts show an increased Atlantic Meridional Overturning Circulation and North Atlantic current transport, which drive an increased advection of warm saline subtropical waters northward, leading to a westward shift of the subpolar front and subsequently a warming and spin down of the SPG. Significant seasonal climate impacts are predicted as the SPG warms, including a reduced sea-ice concentration over the Arctic, an enhanced warming over central US during summer and fall, and a northward shift of the mean ITCZ. These climate anomalies are similar to those observed during a warm phase of the Atlantic Multidecadal Oscillation, which is encouraging for future predictions of North Atlantic climate.
In our original paper (Vecchi et al., 2013, hereafter V13) we stated “the skill in the initialized forecasts comes in large part from the persistence of the mid-1990s shift by the initialized forecasts, rather than from predicting its evolution”. Smith et al (2013, hereafter S13) challenge that assertion, contending that DePreSys was able to make a successful retrospective forecast of that shift. We stand by our original assertion, and present additional analyses using output from DePreSys retrospective forecasts to support our assessment.
Tropical cyclones (TCs) are a hazard to life and property and a prominent element of the global climate system, therefore understanding and predicting TC location, intensity and frequency is of both societal and scientific significance. Methodologies exist to predict basin-wide, seasonally-aggregated TC activity months, seasons and even years in advance. We show that a newly developed high-resolution global climate model can produce skillful forecasts of seasonal TC activity on spatial scales finer than basin-wide, from months and seasons in advance of the TC season. The climate model used here is targeted at predicting regional climate and the statistics of weather extremes on seasonal to decadal timescales, and is comprised of high-resolution (50km×50km) atmosphere and land components, and more moderate resolution (~100km) sea ice and ocean components. The simulation of TC climatology and interannual variations in this climate model is substantially improved by correcting systematic ocean biases through “flux-adjustment.” We perform a suite of 12-month duration retrospective forecasts over the 1981-2012 period, after initializing the climate model to observationally-constrained conditions at the start of each forecast period – using both the standard and flux-adjusted versions of the model. The standard and flux-adjusted forecasts exhibit equivalent skill at predicting Northern Hemisphere TC season sea surface temperature, but the flux-adjusted model exhibits substantially improved basin-wide and regional TC activity forecasts, highlighting the role of systematic biases in limiting the quality of TC forecasts. These results suggest that dynamical forecasts of seasonally-aggregated regional TC activity months in advance are feasible.
Observations and climate simulations exhibit epochs of extreme El Niño/Southern Oscillation (ENSO) behavior that can persist for decades. Previous studies have revealed a wide range of ENSO responses to forcings from greenhouse gases, aerosols, and orbital variations – but they have also shown that interdecadal modulation of ENSO can arise even without such forcings. The present study examines the predictability of this intrinsically-generated component of ENSO modulation, using a 4000-year unforced control run from a global coupled GCM (GFDL-CM2.1) with a fairly realistic representation of ENSO. Extreme ENSO epochs from the unforced simulation are reforecast using the same (“perfect”) model, but slightly-perturbed initial conditions. These 40-member reforecast ensembles display potential predictability of the ENSO trajectory, extending up to several years ahead. However, no decadal-scale predictability of ENSO behavior is found. This indicates that multidecadal epochs of extreme ENSO behavior can arise not only intrinsically, but delicately, and entirely at random. Previous work had shown that CM2.1 generates strong, reasonably-realistic, decadally-predictable high-latitude climate signals, as well as tropical and extratropical decadal signals that interact with ENSO. However, those slow variations appear not to lend significant decadal predictability to this model’s ENSO behavior, at least in the absence of external forcings. While the potential implications of these results are sobering for decadal predictability, they also suggest an expedited approach to model evaluation and development – in which large ensembles of short runs are executed in parallel, to quickly and robustly evaluate simulations of ENSO. Further implications are discussed for decadal prediction, attribution of past and future ENSO variations, and societal vulnerability.
Regional surface temperature trends from the CMIP3 and CMIP5 20th century runs are compared with observations -- at spatial scales ranging from global averages to individual grid points -- using simulated intrinsic climate variability from pre-industrial control runs to assess whether observed trends are detectable and/or consistent with the models’ historical run trends. The CMIP5 models are also used to detect anthropogenic components of the observed trends, by assessing alternative hypotheses based on scenarios driven with either anthropogenic plus natural forcings combined, or with natural forcings only. Modeled variability is assessed via inspection of control run time series, standard deviation maps, spectral analyses, and low-frequency variance consistency tests. The models are found to provide plausible representations of internal climate variability, though there is room for improvement. The influence of observational uncertainty on the trends is assessed, and found to be generally small compared to intrinsic climate variability.
Observed temperature trends over 1901-2010 are found to contain detectable anthropogenic warming components over a large fraction (about 80%) of the analyzed global area. In several regions, the observed warming is significantly underestimated by the models, including parts of the southern Ocean, south Atlantic, far eastern Atlantic, and far west Pacific. Regions without detectable warming signals include the high latitude North Atlantic, the eastern U.S., and parts of the eastern Pacific. For 1981-2010, the observed warming trends over about 45% of the globe are found to contain a detectable anthropogenic warming; this includes much of the globe within about 40-45 degrees of the equator, except for the eastern Pacific.
Knutson, Thomas R., Fanrong Zeng, and Andrew T Wittenberg, September 2013: The Extreme March-2012 Warm Anomaly over the Eastern United States: Global Context and Multimodel Trend Analysis [in “Explaining Extreme Events of 2012 from a Climate Perspective”]. Bulletin of the American Meteorological Society, 94(9), S13-S17.
The impact of climate warming on the upper layer of the Bering Sea is investigated by using a high-resolution coupled global climate model. The model is forced by increasing atmospheric CO2 at a rate of 1% per year until CO2 reaches double its initial value (after 70 years), after which it is held constant. In response to this forcing, the upper layer of the Bering Sea warms by about 2�C in the southeastern shelf and by a little more than 1�C in the western basin. The wintertime ventilation to the permanent thermocline weakens in the western Bering Sea. After CO2 doubling, the southeastern shelf of the Bering Sea becomes almost ice-free in March, and the stratification of the upper layer strengthens in May and June. Changes of physical condition due to the climate warming would impact the pre-condition of spring bio-productivity in the southeastern shelf.
Retrospective predictions of multi-year North Atlantic hurricane frequency are explored, by applying a hybrid statistical-dynamical forecast system to initialized and non-initialized multi-year forecasts of tropical Atlantic and tropical mean sea surface temperatures (SSTs) from two global climate model forecast systems. By accounting for impacts of initialization and radiative forcing, retrospective predictions of five-year mean and nine-year mean tropical Atlantic hurricane frequency show significant correlation relative to a null hypothesis of zero correlation. The retrospective correlations are increased in a two-model average forecast and by using a lagged-ensemble approach, with the two-model ensemble decadal forecasts hurricane frequency over 1961-2011 yielding correlation coefficients that approach 0.9.
These encouraging retrospective multi-year hurricane predictions, however, should be interpreted with care: although initialized forecasts have higher nominal skill than uninitialized ones, the relatively short record and large autocorrelation of the time series limits our confidence in distinguishing between the skill due to external forcing and that added by initialization. The nominal increase in correlation in the initialized forecasts relative to the uninitialized experiments is due to improved representation of the multi-year tropical Atlantic SST anomalies. The skill in the initialized forecasts comes in large part from the persistence of a mid-1990s shift by the initialized forecasts, rather than from predicting its evolution. Predicting shifts like that observed in 1994-1995 remains a critical issue for the success of multi-year forecasts of Atlantic hurricane frequency. The retrospective forecasts highlight the possibility that changes in observing system impact forecast performance.
The decadal predictability of sea surface temperature (SST) and 2m air temperature (T2m) in Geophysical Fluid Dynamics Laboratory (GFDL)'s decadal hindcasts, which are part of the Fifth Coupled Model Intercomparison Project experiments, has been investigated using an average predictability time (APT) analysis. Comparison of retrospective forecasts initialized using the GFDL's Ensemble Coupled Data Assimilation system with uninitialized historical forcing simulations using the same model, allows identification of internal multidecadal pattern (IMP) for SST and T2m. The IMP of SST is characterized by an inter-hemisphere dipole, with warm anomalies centered in the North Atlantic subpolar gyre region and North Pacific subpolar gyre region, and cold anomalies centered in the Antarctic Circumpolar Current region. The IMP of T2m is characterized by a general bi-polar seesaw, with warm anomalies centered in Greenland, and cold anomalies centered in Antarctica. The retrospective prediction skill of the initialized system, verified against independent observations, indicates that the IMP of SST may be predictable up to 4 (10) year lead time at 95% (90%) significance level, and the IMP of T2m may be predictable up to 2 (10) years at 95% (90%) significance level. The initialization of multidecadal variations of northward oceanic heat transport in the North Atlantic significantly improves the predictive skill of the IMP. The dominant roles of oceanic internal dynamics in decadal prediction are further elucidated by fixed-forcing experiments, in which radiative forcing is returned to 1961 values. These results point towards the possibility of meaningful decadal climate outlooks using dynamical coupled models, if they are appropriately initialized from a sustained climate observing system.
We present results for simulated climate and climate change from a newly developed high-resolution global climate model (GFDL CM2.5). The GFDL CM2.5 model has an atmospheric resolution of approximately 50 Km in the horizontal, with 32 vertical levels. The horizontal resolution in the ocean ranges from 28 Km in the tropics to 8 Km at high latitudes, with 50 vertical levels. This resolution allows the explicit simulation of some mesoscale eddies in the ocean, particularly at lower latitudes.
We present analyses based on the output of a 280 year control simulation; we also present results based on a 140 year simulation in which atmospheric CO2 increases at 1% per year until doubling after 70 years.
Results are compared to the GFDL CM2.1 climate model, which has somewhat similar physics but coarser resolution. The simulated climate in CM2.5 shows marked improvement over many regions, especially the tropics, including a reduction in the double ITCZ and an improved simulation of ENSO. Regional precipitation features are much improved. The Indian monsoon and Amazonian rainfall are also substantially more realistic in CM2.5.
The response of CM2.5 to a doubling of atmospheric CO2 has many features in common with CM2.1, with some notable differences. For example, rainfall changes over the Mediterranean appear to be tightly linked to topography in CM2.5, in contrast to CM2.1 where the response is more spatially homogeneous. In addition, in CM2.5 the near-surface ocean warms substantially in the high latitudes of the Southern Ocean, in contrast to simulations using CM2.1.
Delworth, Thomas L., and Fanrong Zeng, July 2012: Multicentennial variability of the Atlantic Meridional Overturning Circulation and its climatic influence in a 4000 year simulation of the GFDL CM2.1 climate model. Geophysical Research Letters, 39, L13702, DOI:10.1029/2012GL052107. Abstract
We investigate decadal to multicentennial variability of Northern Hemisphere surface air temperature in a 4000-year control simulation of the GFDL CM2.1 climate model. Spectral analysis shows the presence of a distinct multicentennial timescle of temperature variability. The associated spatial pattern is broad, covering the entire Northern Hemisphere extratropics, but with enhanced amplitude in the Atlantic and Arctic sectors. This variability appears to be driven by interhemispheric fluctuations in oceanic heat transport associated with the Atlantic Meridional Overturning Circulation (AMOC). The AMOC variability is associated with century-scale propagation of salinity anomalies from the Southern Ocean to the subpolar North Atlantic, with out of phase transport variations between the upper ocean and deeper layers of the Atlantic. When positive (negative) upper ocean salinity anomalies reach the subpolar North Atlantic they strengthen (weaken) the AMOC by modulating upper ocean density and vertical stratification. The large-scale warming also appears to be enhanced by reductions in surface albedo associated with reduced sea-ice and low-level cloudiness, thereby increasing the absorption of shortwave radiation and amplifying the warming from AMOC changes. We speculate that such multicentennial variations in the AMOC could contribute to long-time scale climate fluctuations in the observed paleo record. This could arise purely as internal variability of the climate system, or through radiatively-induced changes to atmospheric circulation patterns, such as the NAO, that would in turn influence the AMOC.
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.
The sensitivity of the North Atlantic Ocean Circulation to an abrupt change in the Nordic Sea overflow is investigated for the first time using a high resolution eddy-permitting global coupled ocean-atmosphere model (GFDL CM2.5). The Nordic Sea overflow is perturbed through the change of the bathymetry in GFDL CM2.5. We analyze the Atlantic Meridional Overturning Circulation (AMOC) adjustment process and the downstream oceanic response to the perturbation. The results suggest that north of 34N, AMOC changes induced by changes in the Nordic Sea overflow propagate on the slow tracer advection time scale, instead of the fast Kelvin wave time scale, resulting in a time lead of several years between subpolar and subtropical AMOC changes. The results also show that a stronger and deeper-penetrating Nordic Sea overflow leads to stronger and deeper AMOC, stronger northward ocean heat transport, reduced Labrador Sea deep convection, stronger cyclonic Northern Recirculation Gyre (NRG), westward shift of the North Atlantic Current (NAC) and southward shift of the Gulf Stream, warmer sea surface temperature (SST) east of Newfoundland and colder SST south of the Grand Banks, stronger and deeper NAC and Gulf Stream, and stronger oceanic eddy activities along the NAC and the Gulf Stream paths. A stronger/weaker Nordic Sea overflow also leads to a contracted/expanded subpolar gyre (SPG). This sensitivity study points to the important role of the Nordic Sea overflow in the large scale North Atlantic ocean circulation, and it is crucial for climate models to have a correct representation of the Nordic Sea overflow.
Simulations from a fine-resolution global coupled model, the Geophysical Fluid Dynamics Laboratory
Climate Model, version 2.4 (CM2.4), are presented, and the results are compared with a coarse version of the
same coupled model, CM2.1, under idealized climate change scenarios. A particular focus is given to the
dynamical response of the Southern Ocean and the role played by the eddies—parameterized or permitted—
in setting the residual circulation and meridional density structure. Compared to the case in which eddies are
parameterized and consistent with recent observational and idealized modeling studies, the eddy-permitting
integrations of CM2.4 show that eddy activity is greatly energized with increasing mechanical and buoyancy
forcings, buffering the ocean to atmospheric changes, and the magnitude of the residual oceanic circulation
response is thus greatly reduced. Although compensation is far from being perfect, changes in poleward eddy
fluxes partially compensate for the enhanced equatorward Ekman transport, leading to weak modifications in
local isopycnal slopes, transport by the Antarctic Circumpolar Current, and overturning circulation. Since the
presence of active ocean eddy dynamics buffers the oceanic response to atmospheric changes, the associated
atmospheric response to those reduced ocean changes is also weakened. Further, it is hypothesized that
present numerical approaches for the parameterization of eddy-induced transports could be too restrictive
and prevent coarse-resolution models from faithfully representing the eddy response to variability and change
in the forcing fields.
The fast and slow components of global warming in a comprehensive climate model are isolated by examining the response to an instantaneous return to pre-industrial forcing. The response is characterized by an initial fast exponential decay with an e-folding time smaller than 5 years, leaving behind a remnant that evolves more slowly. The slow component is estimated to be small at present, as measured by the global mean near-surface air temperature, and, in the model examined, grows to 0.4C by 2100 in the A1B SRES scenario and then to 1.4C by 2300 if one holds radiative forcing fixed after 2100. The dominance of the fast component at present is supported by examining the response to an instantaneous doubling of CO2 and by the excellent fit to the model's ensemble mean 20th century evolution with a simple one-box model with no long times scales.
Understanding the plausible causes for the observed high dust concentrations in Antarctic ice cores during
the Last Glacial Maximum (LGM) is crucial for interpreting the Antarctic dust records in the past climates
and could provide insights into dust variability in future climates. Using the Geophysical Fluid Dynamics
Laboratory (GFDL) General Circulation Models, we conduct an investigation into the various factors
modulating dust emission, transport and deposition, with a view towards an improved quantification of the
LGM dust enhancements in the Antarctic ice cores. The model simulations show that the expansion of
source areas and changes in the Antarctic ice accumulation rates together can account for most of the
observed increase of dust concentrations in the Vostok, Dome C and Taylor Dome cores, but there is an
overestimate of the LGM-to-present ratio in the case of the Byrd core. The source expansion due to the
lowering of sea level yields a factor of 2–3 higher contribution than that due to the reduction of continental
vegetation. The changes in other climate parameters (e.g., SH precipitation change) are estimated to be
relatively less important within the context of this sensitivity study, while the model-simulated LGM
surface winds yield a 20–30 % reduction rather than an increase in dust deposition in Antarctica. This
research yields insights towards a fundamental understanding of the causes for the significant enhancement
of the dust deposition in the Antarctic ice cores during the LGM.
Sulfate aerosols resulting from strong volcanic explosions last for 2–3 years in the lower stratosphere. Therefore it was traditionally believed that volcanic impacts produce mainly short-term, transient climate perturbations. However, the ocean integrates volcanic radiative cooling and responds over a wide range of time scales. The associated processes, especially ocean heat uptake, play a key role in ongoing climate change. However, they are not well constrained by observations, and attempts to simulate them in current climate models used for climate predictions yield a range of uncertainty. Volcanic impacts on the ocean provide an independent means of assessing these processes. This study focuses on quantification of the seasonal to multidecadal time scale response of the ocean to explosive volcanism. It employs the coupled climate model CM2.1, developed recently at the National Oceanic and Atmospheric Administration's Geophysical Fluid Dynamics Laboratory, to simulate the response to the 1991 Pinatubo and the 1815 Tambora eruptions, which were the largest in the 20th and 19th centuries, respectively. The simulated climate perturbations compare well with available observations for the Pinatubo period. The stronger Tambora forcing produces responses with higher signal-to-noise ratio. Volcanic cooling tends to strengthen the Atlantic meridional overturning circulation. Sea ice extent appears to be sensitive to volcanic forcing, especially during the warm season. Because of the extremely long relaxation time of ocean subsurface temperature and sea level, the perturbations caused by the Tambora eruption could have lasted well into the 20th century.nd sea level, the perturbations caused by the Tambora eruption could last well into the 20th century.
Previous work has suggested that the strength and latitudinal position of the Southern Hemisphere (SH) mid-latitude westerly winds has an important impact on climate and the Atlantic Meridional Overturning Circulation (AMOC). We probe this hypothesis by conducting ensembles of experiments using the GFDL CM2.1 coupled ocean-atmosphere model with altered SH wind stress. We find, consistent with previous work, that enhanced (reduced) and poleward (equatorward) displaced SH westerly winds lead to an AMOC intensification (weakening). While the AMOC takes more than a century to respond fully to the altered SH winds, initial effects in the North Atlantic can occur within a few decades. The AMOC changes generate SST and surface air temperature responses in the North Atlantic and adjacent continental regions. In the Southern Hemisphere, the atmosphere responds to the altered ocean circulation with a further strengthening and poleward movement of the SH winds, thereby constituting a modest positive feedback.
Gnanadesikan, Anand, Joellen L Russell, and Fanrong Zeng, 2007: How does ocean ventilation change under global warming?Ocean Science, 3(1), 43-53. Abstract PDF
Since the upper ocean takes up much of the heat added to the earth system by anthropogenic global warming, one would expect that global warming would lead to an increase in stratification and a decrease in the ventilation of the ocean interior. However, multiple simulations in global coupled climate models using an ideal age tracer which is set to zero in the mixed layer and ages at 1 yr/yr outside this layer show that the intermediate depths in the low latitudes, Northwest Atlantic, and parts of the Arctic Ocean become younger under global warming. This paper reconciles these apparently contradictory trends, showing that the decreases result from changes in the relative contributions of old deep waters and younger surface waters. Implications for the tropical oxygen minimum zones, which play a critical role in global biogeochemical cycling are considered in detail.
The formulation and simulation characteristics of two new global coupled climate models developed at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL) are described. The models were designed to simulate atmospheric and oceanic climate and variability from the diurnal time scale through multicentury climate change, given our computational constraints. In particular, an important goal was to use the same model for both experimental seasonal to interannual forecasting and the study of multicentury global climate change, and this goal has been achieved.
Two versions of the coupled model are described, called CM2.0 and CM2.1. The versions differ primarily in the dynamical core used in the atmospheric component, along with the cloud tuning and some details of the land and ocean components. For both coupled models, the resolution of the land and atmospheric components is 2° latitude × 2.5° longitude; the atmospheric model has 24 vertical levels. The ocean resolution is 1° in latitude and longitude, with meridional resolution equatorward of 30° becoming progressively finer, such that the meridional resolution is 1/3° at the equator. There are 50 vertical levels in the ocean, with 22 evenly spaced levels within the top 220 m. The ocean component has poles over North America and Eurasia to avoid polar filtering. Neither coupled model employs flux adjustments.
The control simulations have stable, realistic climates when integrated over multiple centuries. Both models have simulations of ENSO that are substantially improved relative to previous GFDL coupled models. The CM2.0 model has been further evaluated as an ENSO forecast model and has good skill (CM2.1 has not been evaluated as an ENSO forecast model). Generally reduced temperature and salinity biases exist in CM2.1 relative to CM2.0. These reductions are associated with 1) improved simulations of surface wind stress in CM2.1 and associated changes in oceanic gyre circulations; 2) changes in cloud tuning and the land model, both of which act to increase the net surface shortwave radiation in CM2.1, thereby reducing an overall cold bias present in CM2.0; and 3) a reduction of ocean lateral viscosity in the extratropics in CM2.1, which reduces sea ice biases in the North Atlantic.
Both models have been used to conduct a suite of climate change simulations for the 2007 Intergovernmental Panel on Climate Change (IPCC) assessment report and are able to simulate the main features of the observed warming of the twentieth century. The climate sensitivities of the CM2.0 and CM2.1 models are 2.9 and 3.4 K, respectively. These sensitivities are defined by coupling the atmospheric components of CM2.0 and CM2.1 to a slab ocean model and allowing the model to come into equilibrium with a doubling of atmospheric CO2. The output from a suite of integrations conducted with these models is freely available online (see http://nomads.gfdl.noaa.gov/).
Manuscript received 8 December 2004, in final form 18 March 2005
The current generation of coupled climate models run at the Geophysical Fluid Dynamics Laboratory (GFDL) as part of the Climate Change Science Program contains ocean components that differ in almost every respect from those contained in previous generations of GFDL climate models. This paper summarizes the new physical features of the models and examines the simulations that they produce. Of the two new coupled climate model versions 2.1 (CM2.1) and 2.0 (CM2.0), the CM2.1 model represents a major improvement over CM2.0 in most of the major oceanic features examined, with strikingly lower drifts in hydrographic fields such as temperature and salinity, more realistic ventilation of the deep ocean, and currents that are closer to their observed values. Regional analysis of the differences between the models highlights the importance of wind stress in determining the circulation, particularly in the Southern Ocean. At present, major errors in both models are associated with Northern Hemisphere Mode Waters and outflows from overflows, particularly the Mediterranean Sea and Red Sea.
The climate response to idealized changes in the atmospheric CO2 concentration by the new GFDL climate model (CM2) is documented. This new model is very different from earlier GFDL models in its parameterizations of subgrid-scale physical processes, numerical algorithms, and resolution. The model was constructed to be useful for both seasonal-to-interannual predictions and climate change research. Unlike previous versions of the global coupled GFDL climate models, CM2 does not use flux adjustments to maintain a stable control climate. Results from two model versions, Climate Model versions 2.0 (CM2.0) and 2.1 (CM2.1), are presented.
Two atmosphere–mixed layer ocean or slab models, Slab Model versions 2.0 (SM2.0) and 2.1 (SM2.1), are constructed corresponding to CM2.0 and CM2.1. Using the SM2 models to estimate the climate sensitivity, it is found that the equilibrium globally averaged surface air temperature increases 2.9 (SM2.0) and 3.4 K (SM2.1) for a doubling of the atmospheric CO2 concentration. When forced by a 1% per year CO2 increase, the surface air temperature difference around the time of CO2 doubling [transient climate response (TCR)] is about 1.6 K for both coupled model versions (CM2.0 and CM2.1). The simulated warming is near the median of the responses documented for the climate models used in the 2001 Intergovernmental Panel on Climate Change (IPCC) Working Group I Third Assessment Report (TAR).
The thermohaline circulation (THC) weakened in response to increasing atmospheric CO2. By the time of CO2 doubling, the weakening in CM2.1 is larger than that found in CM2.0: 7 and 4 Sv (1 Sv 106 m3 s−1), respectively. However, the THC in the control integration of CM2.1 is stronger than in CM2.0, so that the percentage change in the THC between the two versions is more similar. The average THC change for the models presented in the TAR is about 3 or 4 Sv; however, the range across the model results is very large, varying from a slight increase (+2 Sv) to a large decrease (−10 Sv).
We present results from a series of ensemble integrations of a global coupled atmosphere-ocean model for the period 1865-1997. Each ensemble consists of three integrations initialized from different points in a long-running GFDL R30 coupled model control simulation. The first ensemble includes time-varying forcing from greenhouse gases only. In the remaining three ensembles, forcings from anthropogenic sulfate aerosols, solar variability, and volcanic aerosols in the stratosphere are added progressively, such that the fourth ensemble uses all four of these forcings. The effects of anthropogenic sulfate aerosols are represented by changes in surface albedo, and the effects of volcanic aerosols are represented by latitude-dependent perturbations in incident solar radiation. Comparisons with observations reveal that the addition of the natural forcings (solar and volcanic) improves the simulation of global multidecadal trends in temperature, precipitation, and ocean heat content. Solar and volcanic forcings are important contributors to early twentieth century warming. Volcanic forcing reduces the warming simulated for the late twentieth century. Interdecadal variations in global mean surface air temperature from the ensemble of experiments with all four forcings are very similar to observed variations during most of the twentieth century. The improved agreement of simulated and observed temperature trends when natural climate forcings are included supports the climatic importance of variations in radiative forcing during the twentieth century.
