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 current GFDL seasonal prediction system, the Seamless System for Prediction and Earth System Research (SPEAR), has shown skillful prediction of Arctic sea ice extent with atmosphere and ocean constrained by observations. In this study we present improvements in subseasonal and seasonal predictions of Arctic sea ice by directly assimilating sea ice observations. The sea ice initial conditions from a data assimilation (DA) system that assimilates satellite sea ice concentration (SIC) observations are used to produce a set of reforecast experiments (IceDA) starting from the first day of each month from 1992 to 2017. Our evaluation of daily sea ice extent prediction skill concludes that the SPEAR system generally outperforms the anomaly persistence forecast at lead times beyond 1 month. We primarily focus our analysis on daily gridcell-level sea ice fields. SIC DA improves prediction skill of SIC forecasts prominently in the June-, July-, August-, and September-initialized reforecasts. We evaluate two additional user-oriented metrics: the ice-free probability (IFP) and ice-free date (IFD). IFP is the probability of a grid cell experiencing ice-free conditions in a given year, and IFD is the first date on which a grid cell is ice free. A combined analysis of IFP and IFD demonstrates that the SPEAR model can make skillful predictions of local ice melt as early as May, with modest improvements from SIC DA.

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.

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.

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.

Previous studies have shown the existence of internal multidecadal variability in the Southern Ocean using multiple climate models. This variability, associated with deep ocean convection, can have significant climate impacts. In this work, we use sensitivity studies based on Geophysical Fluid Dynamics Laboratory (GFDL) models to investigate the linkage of this internal variability with the background ocean mean state. We find that mean ocean stratification in the subpolar region that is dominated by mean salinity influences whether this variability occurs, as well as its time scale. The weakening of background stratification favors the occurrence of deep convection. For background stratification states in which the low-frequency variability occurs, weaker ocean stratification corresponds to shorter periods of variability and vice versa. The amplitude of convection variability is largely determined by the amount of heat that can accumulate in the subsurface ocean during periods of the oscillation without deep convection. A larger accumulation of heat in the subsurface reservoir corresponds to a larger amplitude of variability. The subsurface heat buildup is a balance between advection that supplies heat to the reservoir and vertical mixing/convection that depletes it. Subsurface heat accumulation can be intensified both by an enhanced horizontal temperature advection by the Weddell Gyre and by an enhanced ocean stratification leading to reduced vertical mixing and surface heat loss. The paleoclimate records over Antarctica indicate that this multidecadal variability has very likely happened in past climates and that the period of this variability may shift with different climate background mean state.

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 the future, and the related climate phenomena, for planning and mitigation purposes. Here, a suite of simulations from GFDL high-resolution coupled climate models are used 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 the future. The climate of twentieth-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 fails to show evidence for an 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 El Niño nor La Niña enhances the probability of extreme floods over Chennai. However, a warm Bay of Bengal tends to increase the odds of occurrence of extreme Chennai-like floods. In order to trigger a Chennai like-flood, a conducive weather event, such as a tropical depression over the Bay of Bengal with strong transport of moisture from a moist atmosphere over the warm Bay, is necessary for the intense precipitation.

This study proposes an integrated diagnostic framework based on atmospheric circulation regime spatial patterns and frequencies of occurrence to facilitate the identification of model systematic errors across multiple time scales. To illustrate the approach, three sets of 32-yr-long simulations are analyzed for northeastern North America and for the March–May season using the Geophysical Fluid Dynamics Laboratory’s Low Ocean–Atmosphere Resolution (LOAR) and Forecast-Oriented Low Ocean Resolution (FLOR) coupled models; the main difference between these two models is the horizontal resolution of the atmospheric model used. Regime-dependent biases are explored in the light of different atmospheric horizontal resolutions and under different nudging approaches. It is found that both models exhibit a fair representation of the observed circulation regime spatial patterns and frequencies of occurrence, although some biases are present independently of the horizontal resolution or the nudging approach and are associated with a misrepresentation of troughs centered north of the Great Lakes and deep coastal troughs. Moreover, the intraseasonal occurrence of certain model regimes is delayed with respect to observations. On the other hand, interexperiment differences in the mean frequencies of occurrence of the simulated weather types, and their variability across multiple time scales, tend to be negligible. This result suggests that low-resolution models could be of potential use to diagnose and predict physical variables via their simulated weather type characteristics.

This study explores the potential predictability of the Southern Ocean (SO) climate on decadal time scales as represented in the GFDL CM2.1 model using prognostic methods. Perfect model predictability experiments are conducted starting from 10 different initial states, showing potentially predictable variations of Antarctic bottom water (AABW) formation rates on time scales as long as 20 years. The associated Weddell Sea (WS) subsurface temperatures and Antarctic sea ice have potential predictability comparable to that of 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 this perfect predictability study, the authors also make hindcasts of SO decadal variability using the GFDL CM2.1 decadal prediction system. Significant predictive skill for SO SST on multiyear 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 some indication as to where and to what extent skillful decadal SO forecasts might be possible.

This study explores the role of the stratosphere as a source of seasonal predictability of surface climate over Northern Hemisphere extratropics both in the observations and climate model predictions. A suite of numerical experiments, including climate simulations and retrospective forecasts, are set up to isolate the role of the stratosphere in seasonal predictive skill of extratropical near-surface land temperature. It is shown that most of the lead-0-month spring predictive skill of land temperature over extratropics, particularly over northern Eurasia, stems from stratospheric initialization. It is further revealed that this predictive skill of extratropical land temperature arises from skillful prediction of the Arctic Oscillation (AO). The dynamical connection between the stratosphere and troposphere is also demonstrated by the significant correlation between the stratospheric polar vortex and sea level pressure anomalies, as well as the migration of the stratospheric zonal wind anomalies to the lower troposphere.

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 10 yr, 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 10 yr, 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). Additional simulations are used in which heat flux anomalies derived from observed NAO variations from 1901 to 2014 are applied to the ocean component of coupled models. It is shown that ocean dynamics allow models to reproduce important aspects of the observed AMO, mainly in the Subpolar Gyre.

Because of its persistence on seasonal time scales, Arctic sea ice thickness (SIT) is a potential source of predictability for summer sea ice extent (SIE). New satellite observations of SIT represent an opportunity to harness this potential predictability via improved thickness initialization in seasonal forecast systems. In this work, the evolution of Arctic sea ice volume anomalies is studied using a 700-yr control integration and a suite of initialized ensemble forecasts from a fully coupled global climate model. This analysis is focused on the September sea ice zone, as this is the region where thickness anomalies have the potential to impact the SIE minimum. The primary finding of this paper is that, in addition to a general decay with time, sea ice volume anomalies display a summer enhancement, in which anomalies tend to grow between the months of May and July. This summer enhancement is relatively symmetric for positive and negative volume anomalies and peaks in July regardless of the initial month. Analysis of the surface energy budget reveals that the summer volume anomaly enhancement is driven by a positive feedback between the SIT state and the surface albedo. The SIT state affects surface albedo through changes in the sea ice concentration field, melt-onset date, snow coverage, and ice thickness distribution, yielding an anomaly in the total absorbed shortwave radiation between May and August, which enhances the existing SIT anomaly. This phenomenon highlights the crucial importance of accurate SIT initialization and representation of ice–albedo feedback processes in seasonal forecast systems.

This study attempts to improve the prediction of western North Pacific (WNP) and East Asia (EA) landfalling tropical cyclones (TCs) using modes of large-scale climate variability [e.g., the Pacific meridional mode (PMM), the Atlantic meridional mode (AMM), and North Atlantic sea surface temperature anomalies (NASST)] as predictors in a hybrid statistical–dynamical scheme, based on dynamical model forecasts with the GFDL Forecast-Oriented Low Ocean Resolution version of CM2.5 with flux adjustments (FLOR-FA). Overall, the predictive skill of the hybrid model for the WNP TC frequency increases from lead month 5 (initialized in January) to lead month 0 (initialized in June) in terms of correlation coefficient and root-mean-square error (RMSE). The hybrid model outperforms FLOR-FA in predicting WNP TC frequency for all lead months. The predictive skill of the hybrid model improves as the forecast lead time decreases, with values of the correlation coefficient increasing from 0.56 for forecasts initialized in January to 0.69 in June. The hybrid models for landfalling TCs over the entire East Asian (EEA) coast and its three subregions [i.e., southern EA (SEA), middle EA (MEA), and northern EA (NEA)] dramatically outperform FLOR-FA. The correlation coefficient between predicted and observed TC landfall over SEA increases from 0.52 for forecasts initialized in January to 0.64 in June. The hybrid models substantially reduce the RMSE of landfalling TCs over SEA and EEA compared with FLOR-FA. This study suggests that the PMM and NASST/AMM can be used to improve statistical/hybrid forecast models for the frequencies of WNP or East Asia landfalling TCs.

The 2015 hurricane season in the eastern and central Pacific Ocean (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 tropical Pacific warming, but by warming in the subtropical Pacific Ocean. This warming is not typical of El Niño, but rather of 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 EPO, CPO, and Hawaii as a result of enhanced subtropical Pacific warming from anthropogenic greenhouse gas forcing.

The wintertime Arctic temperature (T; surface–400 hPa) decreased from 1979 to 1997 and increased rapidly from 1998 to 2012, in contrast to the global mean surface air temperature. Here aspects of circulation variability that are associated with these temperature changes are examined using the NCEP–NCAR reanalysis and ERA-Interim products. It is found that the Nordic–Siberia seesaw of meridional winds near 70°N is associated with two-thirds of the variance of the Arctic winter mean T, possibly contributing to the cooling and warming trends. It is suggested here that the seesaw accounts for much of the difference in Arctic amplification between observations and climate models. Growth of sea ice in winter is hindered by southerly winds over the Nordic region (0°–60°E). Through modulation of the wind seesaw, the eastern Atlantic (EA) pattern is found to be significantly associated with Arctic and East Asia winter climate variations. In one phase of the EA pattern, a midlatitude North Atlantic ridge anomaly is associated with a poleward shift of the mean storm track, a weakened eddy-driven jet over Eurasia, and above-normal sea level pressure (SLP) over Siberia, most significantly in the region to the northwest of Lake Baikal. The EA pattern is associated with two-thirds of the variance of winter-average SLP over Siberia.

Skillful seasonal forecasting of tropical cyclone (TC; wind speed ≥17.5 m s−1) activity is challenging, even more so when the focus is on major hurricanes (wind speed ≥49.4 m s−1), the most intense hurricanes (category 4 and 5; wind speed ≥58.1 m s–1), and landfalling TCs. This study shows that a 25-km-resolution global climate model [High-Resolution Forecast-Oriented Low Ocean Resolution (FLOR) model (HiFLOR)] developed at the Geophysical Fluid Dynamics Laboratory (GFDL) has improved skill in predicting the frequencies of major hurricanes and category 4 and 5 hurricanes in the North Atlantic as well as landfalling TCs over the United States and Caribbean islands a few months in advance, relative to its 50-km-resolution predecessor climate model (FLOR). HiFLOR also shows significant skill in predicting category 4 and 5 hurricanes in the western North Pacific and eastern North Pacific, while both models show comparable skills in predicting basin-total and landfalling TC frequency in the basins. The improved skillful forecasts of basin-total TCs, major hurricanes, and category 4 and 5 hurricane activity in the North Atlantic by HiFLOR are obtained mainly by improved representation of the TCs and their response to climate from the increased horizontal resolution rather than by improvements in large-scale parameters.

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 the key drivers. The atmospheric (land) initial conditions play a major (minor) role in the central and northwestern (eastern) United States. Because of changes in radiative forcing, the probability of areal-averaged summer temperature anomalies over the United States exceeding the observed 2012 anomaly increases with time over the early twenty-first century. La Niña (El Niño) events tend to increase (reduce) the occurrence rate of heat waves. The temperatures over the central United States 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.

Uncertainty in cumulus convection parameterization is one of the most important causes of model climate drift through interactions between large-scale background and local convection that use empirically set parameters. Without addressing the large-scale feedback, the calibrated parameter values within a convection scheme are usually not optimal for a climate model. This study first designs a multiple-column atmospheric model that includes large-scale feedbacks for cumulus convection and then explores the role of large-scale feedbacks in cumulus convection parameter estimation using an ensemble filter. The performance of convection parameter estimation with or without the presence of large-scale feedback is examined. It is found that including large-scale feedbacks in cumulus convection parameter estimation can significantly improve the estimation quality. This is because large-scale feedbacks help transform local convection uncertainties into global climate sensitivities, and including these feedbacks enhances the statistical representation of the relationship between parameters and state variables. The results of this study provide insights for further understanding of climate drift induced from imperfect cumulus convection parameterization, which may help improve climate modeling.

The seasonal predictability of extratropical storm tracks in the Geophysical Fluid Dynamics Laboratory’s (GFDL)’s high-resolution climate model has been investigated using an average predictability time analysis. The leading predictable components of extratropical storm tracks are the ENSO-related spatial patterns 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 2-m air temperature) in response to different ENSO phases. These results point toward the possibility of providing skillful seasonal predictions of the statistics of extratropical extremes over land using high-resolution coupled models.

This study demonstrates skillful seasonal prediction of 2-m air temperature and precipitation over land in a new high-resolution climate model developed by the Geophysical Fluid Dynamics Laboratory and explores the possible sources of the skill. The authors 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, the improved skill of the high-resolution model over the previous lower-resolution model in seasonal prediction of the Niño-3.4 index and other aspects of interest is shown. Then, the skill of temperature and precipitation in the high-resolution model for boreal winter and summer is measured, and the sources of the skill are diagnosed. Last, predictions are reconstructed using a few of the most predictable components to yield more skillful predictions than the raw model predictions. Over three decades of hindcasts, 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 2-m air temperature and precipitation over land.

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 basinwide, seasonally aggregated TC activity months, seasons, and even years in advance. It is shown that a newly developed high-resolution global climate model can produce skillful forecasts of seasonal TC activity on spatial scales finer than basinwide, 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 time scales, and comprises high-resolution (50 km × 50 km) atmosphere and land components as well as more moderate-resolution (~100 km) 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.” A suite of 12-month duration retrospective forecasts is performed 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 basinwide 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.

Decadal prediction experiments were conducted as part of phase 5 of the Coupled Model Intercomparison Project (CMIP5) using the GFDL Climate Model, version 2.1 (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 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-1990s shift, which is not captured by the uninitialized forecasts. All the hindcasts initialized in the early 1990s 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 predicting the mid-1990s 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 spindown 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 the central United States 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.

Given a biased coupled model and the atmospheric and oceanic observing system, maintaining a balanced and coherent climate estimation is of critical importance for producing accurate climate analysis and prediction initialization. However, because of limitations of the observing system (e.g., most of the oceanic measurements are only available for the upper ocean), directly evaluating climate estimation with real observations is difficult. With two coupled models that are biased with respect to each other, a biased twin experiment is designed to simulate the problem. To do that, the atmospheric and oceanic observations drawn from one model based on the modern climate observing system are assimilated into the other. The model that produces observations serves as the truth and the degree by which an assimilation recovers the truth steadily and coherently is an assessment of the impact of the data constraint scheme on climate estimation. Given the assimilation model bias of warmer atmosphere and colder ocean, where the atmospheric-only (oceanic only) data constraint produces an overcooling (overwarming) ocean through the atmosphere–ocean interaction, the constraints with both atmospheric and oceanic data create a balanced and coherent ocean estimate as the observational model. Moreover, the consistent atmosphere–ocean constraint produces the most accurate estimate for North Atlantic Deep Water (NADW), whereas NADW is too strong (weak) if the system is only constrained by atmospheric (oceanic) data. These twin experiment results provide insights that consistent data constraints of multiple components are very important when a coupled model is combined with the climate observing system for climate estimation and prediction initialization.

The decadal predictability of sea surface temperature (SST) and 2-m air temperature (T2m) in the Geophysical Fluid Dynamics Laboratory (GFDL) 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 Ensemble Coupled Data Assimilation system with uninitialized historical forcing simulations using the same model allows identification of the internal multidecadal pattern (IMP) for SST and T2m. The IMP of SST is characterized by an interhemisphere 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 bipolar seesaw, with warm anomalies centered in Greenland and cold anomalies centered in Antarctica. The retrospective prediction skill of the initialized system, verified against independent observational datasets, indicates that the IMP of SST may be predictable up to 4 (10) yr lead time at 95% (90%) significance level, and the IMP of T2m may be predictable up to 2 (10) yr at the 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 abruptly to 1961 values. These results point toward the possibility of meaningful decadal climate outlooks using dynamical coupled models if they are appropriately initialized from a sustained climate observing system.

This paper applies a new field significance test to establish the existence and consistency of ENSO teleconnection patterns across models and observations. An ENSO teleconnection pattern is defined as a field of regression coefficients between an index of the tropical Pacific sea surface temperature and a field of variables such as surface air temperature or precipitation. The test is applied to boreal winter and summer in six continents using observations and hindcasts from the Development of a European Multimodel Ensemble System for Seasonal-to-Interannual Prediction (DEMETER) and the ENSEMBLE-based predictions of climate changes and their impacts (ENSEMBLES) projects. This comparison represents one of the most comprehensive and up-to-date assessments of the extent to which ENSO teleconnection patterns exist and can be reproduced by coupled models. Statistically significant ENSO teleconnection patterns are detected in both observations and models and in all continents and in both winter and summer seasons, except in two cases: 1) Europe (both seasons and variables), and 2) North America (both variables in boreal summer). Despite many ENSO teleconnection patterns being significant, however, the patterns do not necessarily agree between observations and models. The degree of agreement between models and observations is characterized as “robust,” “moderate,” or “low.” Only Australia and South America are found to have robust agreement between ENSO teleconnection patterns, and then only for limited seasons and variables. Although many of our conclusions regarding teleconnection patterns conform to previous studies, there are exceptions, including the fact that the teleconnection for boreal winter precipitation is generally accepted to exist in Africa but in fact has only low agreement with model simulations, while that in Asia is not widely recognized to exist but is found to be significant and in moderate agreement with model teleconnections.

Regression patterns often are used to diagnose the relation between a field and a climate index, but a significance test for the pattern “as a whole” that accounts for the multiplicity and interdependence of the tests has not been widely available. This paper argues that field significance can be framed as a test of the hypothesis that all regression coefficients vanish in a suitable multivariate regression model. A test for this hypothesis can be derived from the generalized likelihood ratio test. The resulting statistic depends on relevant covariance matrices and accounts for the multiplicity and interdependence of the tests. It also depends only on the canonical correlations between the predictors and predictands, thereby revealing a fundamental connection to canonical correlation analysis. Remarkably, the test statistic is invariant to a reversal of the predictors and predictands, allowing the field significance test to be reduced to a standard univariate hypothesis test. In practice, the test cannot be applied when the number of coefficients exceeds the sample size, reflecting the fact that testing more hypotheses than data is ill conceived. To formulate a proper significance test, the data are represented by a small number of principal components, with the number chosen based on cross-validation experiments. However, instead of selecting the model that minimizes the cross-validated mean square error, a confidence interval for the cross-validated error is estimated and the most parsimonious model whose error is within the confidence interval of the minimum error is chosen. This procedure avoids selecting complex models whose error is close to much simpler models. The procedure is applied to diagnose long-term trends in annual average sea surface temperature and boreal winter 300-hPa zonal wind. In both cases a statistically significant 50-yr trend pattern is extracted. The resulting spatial filter can be used to monitor the evolution of the regression pattern without temporal filtering.