Search

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.

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.

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.

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.

This study examines historical simulations of ENSO in the E3SM-1-0, CESM2, and GFDL-CM4 climate models, provided by three leading U.S. modeling centers as part of the Coupled Model Intercomparison Project phase 6 (CMIP6). These new models have made substantial progress in simulating ENSO’s key features, including: amplitude; timescale; spatial patterns; phase-locking; spring persistence barrier; and recharge oscillator dynamics. However, some important features of ENSO are still a challenge to simulate. In the central and eastern equatorial Pacific, the models’ weaker-than-observed subsurface zonal current anomalies and zonal temperature gradient anomalies serve to weaken the nonlinear zonal advection of subsurface temperatures, leading to insufficient warm/cold asymmetry of ENSO’s sea surface temperature anomalies (SSTA). In the western equatorial Pacific, the models’ excessive simulated zonal SST gradients amplify their zonal temperature advection, causing their SSTA to extend farther west than observed. The models underestimate both ENSO’s positive dynamic feedbacks (due to insufficient zonal wind stress responses to SSTA) and its thermodynamic damping (due to insufficient convective cloud shading of eastern Pacific SSTA during warm events); compensation between these biases leads to realistic linear growth rates for ENSO, but for somewhat unrealistic reasons. The models also exhibit stronger-than-observed feedbacks onto eastern equatorial Pacific SSTAs from thermocline depth anomalies, which accelerates the transitions between events and shortens the simulated ENSO period relative to observations. Implications for diagnosing and simulating ENSO in climate models are discussed.

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 and 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 timescales.

This study presents an experimental model for Seasonal Probabilistic Outlook for Tornadoes (SPOTter) in the contiguous United States for March, April, and May and evaluates its forecast skill. This forecast model uses the leading empirical orthogonal function modes of regional variability in tornadic environmental parameters (i.e., low-level vertical wind shear and convective available potential energy), derived from the NCEP Coupled Forecast System, version 2, as the primary predictors. A multiple linear regression is applied to the predicted modes of tornadic environmental parameters to estimate U.S. tornado activity, which is presented as the probability for above-, near-, and below-normal categories. The initial forecast is carried out in late February for March–April U.S. tornado activity and then is updated in late March for April–May activity. A series of reforecast skill tests, including the jackknife cross-validation test, shows that the probabilistic reforecast is overall skillful for predicting the above- and below-normal area-averaged activity in the contiguous United States for the target months of both March–April and April–May. The forecast model also successfully reforecasts the 2011 super-tornado-outbreak season and the other three most active U.S. tornado seasons in 1982, 1991, and 2008, and thus it may be suitable for an operational use for predicting future active and inactive U.S. tornado seasons. However, additional tests show that the regional reforecast is skillful for March–April activity only in the Ohio Valley and South and for April–May activity only in the Southeast and Upper Midwest. These and other limitations of the current model, along with the future advances needed to improve the U.S. regional-scale tornado forecast, are discussed.

Climate models often show errors in simulating and predicting tropical cyclone (TC) activity, but the sources of these errors are not well understood. This study proposes an evaluation framework and analyzes three sets of experiments conducted using a seasonal prediction model developed at the Geophysical Fluid Dynamics Laboratory (GFDL). These experiments apply the nudging technique to the model integration and/or initialization to estimate possible improvements from nearly perfect model conditions. The results suggest that reducing sea surface temperature (SST) errors remains important for better predicting TC activity at long forecast leads—even in a flux-adjusted model with reduced climatological biases. Other error sources also contribute to biases in simulated TC activity, with notable manifestations on regional scales. A novel finding is that the coupling and initialization of the land and atmosphere components can affect seasonal TC prediction skill. Simulated year-to-year variations in June land conditions over North America show a significant lead correlation with the North Atlantic large-scale environment and TC activity. Improved land–atmosphere initialization appears to improve the Atlantic TC predictions initialized in some summer months. For short-lead predictions initialized in June, the potential skill improvements attributable to land–atmosphere initialization might be comparable to those achievable with perfect SST predictions. Overall, this study delineates the SST and non-oceanic error sources in predicting TC activity and highlights avenues for improving predictions. The nudging-based evaluation framework can be applied to other models and help improve predictions of other weather extremes.

The majority of future projections in the Coupled Model Intercomparison Project (CMIP5) show more frequent exceedances of the 5 mm day−1 rainfall threshold in the eastern equatorial Pacific rainfall during El Niño, previously described in the literature as an increase in “extreme El Niño events”; however, these exceedance frequencies vary widely across models, and in some projections actually decrease. Here we combine single-model large ensemble simulations with phase 5 of the Coupled Model Intercomparison Project (CMIP5) to diagnose the mechanisms for these differences. The sensitivity of precipitation to local SST anomalies increases consistently across CMIP-class models, tending to amplify extreme El Niño occurrence; however, changes to the magnitude of ENSO-related SST variability can drastically influence the results, indicating that understanding changes to SST variability remains imperative. Future El Niño rainfall intensifies most in models with 1) larger historical cold SST biases in the central equatorial Pacific, which inhibit future increases in local convective cloud shading, enabling more local warming; and 2) smaller historical warm SST biases in the far eastern equatorial Pacific, which enhance future reductions in stratus cloud, enabling more local warming. These competing mechanisms complicate efforts to determine whether CMIP5 models under- or overestimate the future impacts of climate change on El Niño rainfall and its global impacts. However, the relation between future projections and historical biases suggests the possibility of using observable metrics as “emergent constraints” on future extreme El Niño, and a proof of concept using SSTA variance, precipitation sensitivity to SST, and regional SST trends is presented.

Positive precipitation biases over western North America have remained a pervasive problem in the current generation of coupled global climate models. These biases are substantially reduced, however, in a version of the Geophysical Fluid Dynamics Laboratory Forecast-Oriented Low Ocean Resolution (FLOR) coupled climate model with systematic sea surface temperature (SST) biases artificially corrected through flux adjustment. This study examines how the SST biases in the Atlantic and Pacific Oceans contribute to the North American precipitation biases. Experiments with the FLOR model in which SST biases are removed in the Atlantic and Pacific are carried out to determine the contribution of SST errors in each basin to precipitation statistics over North America. Tropical and North Pacific SST biases have a strong impact on northern North American precipitation, while tropical Atlantic SST biases have a dominant impact on precipitation biases in southern North America, including the western United States. Most notably, negative SST biases in the tropical Atlantic in boreal winter induce an anomalously strong Aleutian low and a southward bias in the North Pacific storm track. In boreal summer, the negative SST biases induce a strengthened North Atlantic subtropical high and Great Plains low-level jet. Each of these impacts contributes to positive annual mean precipitation biases over western North America. Both North Pacific and North Atlantic SST biases induce SST biases in remote basins through dynamical pathways, so a complete attribution of the effects of SST biases on precipitation must account for both the local and remote impacts.

Midlatitude baroclinic waves drive extratropical weather and climate variations, but their predictability beyond 2 weeks has been deemed low. Here we analyze a large ensemble of climate simulations forced by observed sea surface temperatures (SSTs) and demonstrate that seasonal variations of baroclinic wave activity (BWA) are potentially predictable. This potential seasonal predictability is denoted by robust BWA responses to SST forcings. To probe regional sources of the potential predictability, a regression analysis is applied to the SST-forced large ensemble simulations. By filtering out variability internal to the atmosphere and land, this analysis identifies both well-known and unfamiliar BWA responses to SST forcings across latitudes. Finally, we confirm the model-indicated predictability by showing that an operational seasonal prediction system can leverage some of the identified SST-BWA relationships to achieve skillful predictions of BWA. Our findings help to extend long-range predictions of the statistics of extratropical weather events and their impacts.

International audience; Large ensembles of model simulations require considerable resources, and thus defining an appropriate ensemble size for a particular application is an important experimental design criterion. We estimate the ensemble size (N) needed to assess a model’s ability to capture observed El Niño-Southern Oscillation (ENSO) behavior by utilizing the recently developed International CLIVAR ENSO Metrics Package. Using the larger ensembles available from CMIP6 and the US CLIVAR Large Ensemble Working Group, we find that larger ensembles are needed to robustly capture baseline ENSO characteristics (N > 50) and physical processes (N > 50) than the background climatology (N ≥ 12) and remote ENSO teleconnections (N ≥ 6). While these results vary somewhat across metrics and models, our study quantifies how larger ensembles are required to robustly evaluate simulated ENSO behavior, thereby providing some guidance for the design of model ensembles.

International audience; Climate variability in the tropical Pacific affects global climate on a wide range of time scales. On interannual time scales, the tropical Pacific is home to the El Niño–Southern Oscillation (ENSO). Decadal variations and changes in the tropical Pacific, referred to here collectively as tropical Pacific decadal variability (TPDV), also profoundly affect the climate system. Here, we use TPDV to refer to any form of decadal climate variability or change that occurs in the atmosphere, the ocean, and over land within the tropical Pacific. “Decadal,” which we use in a broad sense to encompass multiyear through multidecadal time scales, includes variability about the mean state on decadal time scales, externally forced mean-state changes that unfold on decadal time scales, and decadal variations in the behavior of higher-frequency modes like ENSO.

We describe the Geophysical Fluid Dynamics Laboratory's CM4.0 physical climate model, with emphasis on those aspects that may be of particular importance to users of this model and its simulations. The model is built with the AM4.0/LM4.0 atmosphere/land model and OM4.0 ocean model. Topics include the rationale for key choices made in the model formulation, the stability as well as drift of the preindustrial control simulation, and comparison of key aspects of the historical simulations with observations from recent decades. Notable achievements include the relatively small biases in seasonal spatial patterns of top‐of‐atmosphere fluxes, surface temperature, and precipitation; reduced double Intertropical Convergence Zone bias; dramatically improved representation of ocean boundary currents; a high‐quality simulation of climatological Arctic sea ice extent and its recent decline; and excellent simulation of the El Niño‐Southern Oscillation spectrum and structure. Areas of concern include inadequate deep convection in the Nordic Seas; an inaccurate Antarctic sea ice simulation; precipitation and wind composites still affected by the equatorial cold tongue bias; muted variability in the Atlantic Meridional Overturning Circulation; strong 100 year quasiperiodicity in Southern Ocean ventilation; and a lack of historical warming before 1990 and too rapid warming thereafter due to high climate sensitivity and strong aerosol forcing, in contrast to the observational record. Overall, CM4.0 scores very well in its fidelity against observations compared to the Coupled Model Intercomparison Project Phase 5 generation in terms of both mean state and modes of variability and should prove a valuable new addition for analysis across a broad array of applications.

We document the configuration and emergent simulation features from the Geophysical Fluid Dynamics Laboratory (GFDL) OM4.0 ocean/sea ice model. OM4 serves as the ocean/sea ice component for the GFDL climate and Earth system models. It is also used for climate science research and is contributing to the Coupled Model Intercomparison Project version 6 Ocean Model Intercomparison Project. The ocean component of OM4 uses version 6 of the Modular Ocean Model and the sea ice component uses version 2 of the Sea Ice Simulator, which have identical horizontal grid layouts (Arakawa C‐grid). We follow the Coordinated Ocean‐sea ice Reference Experiments protocol to assess simulation quality across a broad suite of climate‐relevant features. We present results from two versions differing by horizontal grid spacing and physical parameterizations: OM4p5 has nominal 0.5° spacing and includes mesoscale eddy parameterizations and OM4p25 has nominal 0.25° spacing with no mesoscale eddy parameterization. Modular Ocean Model version 6 makes use of a vertical Lagrangian‐remap algorithm that enables general vertical coordinates. We show that use of a hybrid depth‐isopycnal coordinate reduces the middepth ocean warming drift commonly found in pure z* vertical coordinate ocean models. To test the need for the mesoscale eddy parameterization used in OM4p5, we examine the results from a simulation that removes the eddy parameterization. The water mass structure and model drift are physically degraded relative to OM4p5, thus supporting the key role for a mesoscale closure at this resolution.

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.

The heat budget of the Pacific equatorial cold tongue (ECT) is explored using the GFDL-FLOR coupled GCM (the forecast-oriented low ocean resolution version of CM2.5) 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 that acts to cool the equatorial source waters; and 2) an excessive SST contrast between the ECT and off-equator areas, 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 either the simulated TIWs or 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 potential for explosive volcanism to affect the state of the El Niño-Southern Oscillation (ENSO) has been debated since the 1980s. Several observational studies, largely based on tree rings, have since found support for a positive ENSO phase in the year following large eruptions. Models of different complexities also simulate such a response, detectable above the backdrop of internal variability – though they disagree on the underlying mechanisms. In contrast, recent coral data from the heart of the tropical Pacific suggest no uniform ENSO response to all eruptions over the last millennium. Here we leverage paleoclimate data assimilation to integrate the latest paleoclimate evidence into a consistent dynamical framework and re-appraise this relationship. Our analysis finds only a weak statistical association between volcanism and ENSO, suggestive of either no causal association, or of an insufficient number of large volcanic events over the past millennium to obtain reliable statistics. While currently available observations do not support the model-based inference that tropical eruptions promote an ENSO response, there are hints of a response to hemispherically asymmetric forcing, consistent with the "ITCZ shift" mechanism. We discuss the difficulties of conclusively establishing a volcanic influence on ENSO given the many degrees of freedom affecting the response, including eruption season, spatial characteristics of the forcing, and ENSO phase preconditioning.

El Niño–Southern Oscillation (ENSO) is the dominant mode of interannual climate variability on the planet, with far-reaching global impacts. It is therefore key to evaluate ENSO simulations in state-of-the-art numerical models used to study past, present, and future climate. Recently, the Pacific Region Panel of the International Climate and Ocean: Variability, Predictability and Change (CLIVAR) Project, as a part of the World Climate Research Programme (WCRP), led a community-wide effort to evaluate the simulation of ENSO variability, teleconnections, and processes in climate models. The new CLIVAR 2020 ENSO metrics package enables model diagnosis, comparison, and evaluation to 1) highlight aspects that need improvement; 2) monitor progress across model generations; 3) help in selecting models that are well suited for particular analyses; 4) reveal links between various model biases, illuminating the impacts of those biases on ENSO and its sensitivity to climate change; and to 5) advance ENSO literacy. By interfacing with existing model evaluation tools, the ENSO metrics package enables rapid analysis of multipetabyte databases of simulations, such as those generated by the Coupled Model Intercomparison Project phases 5 (CMIP5) and 6 (CMIP6). The CMIP6 models are found to significantly outperform those from CMIP5 for 8 out of 24 ENSO-relevant metrics, with most CMIP6 models showing improved tropical Pacific seasonality and ENSO teleconnections. Only one ENSO metric is significantly degraded in CMIP6, namely, the coupling between the ocean surface and subsurface temperature anomalies, while the majority of metrics remain unchanged.

Future changes in the seasonal evolution of El Niño - Southern Oscillation (ENSO) during the onset and decay phases have received little attention by the research community. This work investigates the projected changes in the spatio-temporal evolution of El Niño events in the 21st Century (21C) using a large ensemble simulation of a couple general circulation model under anthropogenic forcing. Here we show that El Niño is projected to (1) initiate sooner in boreal spring, (2) to grow at a faster rate, (3) to persist longer over the eastern and far eastern Pacific, and (4) to have a broader impact on remote teleconnections. Significant changes in the tropical Pacific mean state, dominant feedback processes, and a projected increase in stochastic westerly wind burst forcing largely explain the fast growing and slow dissipating El Niño in the late 21C. Important implications of these findings are that the global climate impacts are projected to become more significant and persistent, owing to the extended persistence of El Niño.

The potential for explosive volcanism to affect the El Niño-Southern Oscillation (ENSO) has been debated since the 1980s. Several observational studies, based largely on tree-ring proxies, have since found support for a positive ENSO phase in the year following large eruptions. In contrast, recent coral data from the heart of the tropical Pacific suggest no uniform ENSO response to explosive volcanism over the last millennium. Here we leverage paleoclimate data assimilation to integrate both tree-ring and coral proxies into a reconstruction of ENSO state, and re-appraise this relationship. We find only a weak statistical association between volcanism and ENSO, and identify the selection of volcanic events as a key variable to the conclusion. We discuss the difficulties of conclusively establishing a volcanic influence on ENSO by empirical means, given the myriad factors affecting the response, including the spatiotemporal details of the forcing and ENSO phase preconditioning.

We describe the baseline coupled model configuration and simulation characteristics of GFDL's Earth System Model Version 4.1 (ESM4.1), which builds on component and coupled model developments at GFDL over 2013–2018 for coupled carbon‐chemistry‐climate simulation contributing to the sixth phase of the Coupled Model Intercomparison Project. In contrast with GFDL's CM4.0 development effort that focuses on ocean resolution for physical climate, ESM4.1 focuses on comprehensiveness of Earth system interactions. ESM4.1 features doubled horizontal resolution of both atmosphere (2° to 1°) and ocean (1° to 0.5°) relative to GFDL's previous‐generation coupled ESM2‐carbon and CM3‐chemistry models. ESM4.1 brings together key representational advances in CM4.0 dynamics and physics along with those in aerosols and their precursor emissions, land ecosystem vegetation and canopy competition, and multiday fire; ocean ecological and biogeochemical interactions, comprehensive land‐atmosphere‐ocean cycling of CO2, dust and iron, and interactive ocean‐atmosphere nitrogen cycling are described in detail across this volume of JAMES and presented here in terms of the overall coupling and resulting fidelity. ESM4.1 provides much improved fidelity in CO2 and chemistry over ESM2 and CM3, captures most of CM4.0's baseline simulations characteristics, and notably improves on CM4.0 in (1) Southern Ocean mode and intermediate water ventilation, (2) Southern Ocean aerosols, and (3) reduced spurious ocean heat uptake. ESM4.1 has reduced transient and equilibrium climate sensitivity compared to CM4.0. Fidelity concerns include (1) moderate degradation in sea surface temperature biases, (2) degradation in aerosols in some regions, and (3) strong centennial scale climate modulation by Southern Ocean convection.

To improve El Niño / Southern Oscillation (ENSO) predictions and projections in a changing climate, it is essential to better understand ENSO’s sensitivities to external radiative forcings. Strong volcanic eruptions can help to clarify ENSO’s sensitivities, mechanisms, and feedbacks. Strong explosive volcanic eruptions inject millions of tons of sulfur dioxide into the stratosphere, where they are converted into sulfate aerosols. For equatorial volcanoes, these aerosols can spread globally, scattering and absorbing incoming sunlight, and inducing a global-mean surface cooling. Despite this global-mean cooling effect, paleo data confirm remarkable warming of the eastern equatorial Pacific in the two years after a tropical eruption, with a shift towards an El Niño-like state. To illuminate this response and explain why it tends to occur during particular seasons and ENSO phases, we present a unified framework that includes the roles of the seasonal cycle, stochastic wind forcing, eruption magnitude, and various tropical Pacific climate feedbacks. Analyzing over 20,000 years of large-ensemble simulations from the GFDL-CM2.1 climate model forced by volcanic eruptions, we find that the ENSO response comprises both stochastic and deterministic components, which vary depending on the perturbation season and the ocean preconditioning. For boreal winter eruptions, stochastic dispersion largely obscures the deterministic response, being the largest for the strong El Niño preconditioning. Deterministic El Niño-like responses to summer eruptions are well seen on neutral ENSO and weak to moderate El Niño preconditioning and grow with the eruption magnitude. The relative balance of these components determines the predictability and strength of the ENSO response. The results clarify why previous studies obtained seemingly conflicting results.