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.
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.
We describe GFDL's CM4.0 physical climate model, with emphasis on those aspects that may be of particular importance to users of this model and its simulations. The model is built with the AM4.0/LM4.0 atmosphere/land model and OM4.0 ocean model. Topics include the rationale for key choices made in the model formulation, the stability as well as drift of the pre‐industrial control simulation, and comparison of key aspects of the historical simulations with observations from recent decades. Notable achievements include the relatively small biases in seasonal spatial patterns of top‐of‐atmosphere fluxes, surface temperature, and precipitation; reduced double Intertropical Convergence Zone bias; dramatically improved representation of ocean boundary currents; a high quality simulation of climatological Arctic sea ice extent and its recent decline; and excellent simulation of the El Niño‐Southern Oscillation spectrum and structure. Areas of concern include inadequate deep convection in the Nordic Seas; an inaccurate Antarctic sea ice simulation; precipitation and wind composites still affected by the equatorial cold tongue bias; muted variability in the Atlantic Meridional Overturning Circulation; strong 100 year quasi‐periodicity in Southern Ocean ventilation; and a lack of historical warming before 1990 and too rapid warming thereafter due to high climate sensitivity and strong aerosol forcing, in contrast to the observational record. Overall, CM4.0 scores very well in its fidelity against observations compared to the Coupled Model Intercomparison Project Phase 5 generation in terms of both mean state and modes of variability and should prove a valuable new addition for analysis across a broad array of applications.
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.
Explosive volcanic eruptions have large climate impacts, and can serve as observable tests of the climatic response to radiative forcing. Using a high resolution climate model, we contrast the climate responses to Pinatubo, with symmetric forcing, and those to Santa Maria and Agung, which had meridionally asymmetric forcing. Although Pinatubo had larger global‐mean forcing, asymmetric forcing strongly shifts the latitude of tropical rainfall features, leading to larger local precipitation/TC changes. For example, North Atlantic TC activity over is enhanced/reduced by SH‐forcing (Agung)/NH‐forcing (Santa Maria), but changes little in response to the Pinatubo forcing. Moreover, the transient climate sensitivity estimated from the response to Santa Maria is 20% larger than that from Pinatubo or Agung. This spread in climatic impacts of volcanoes needs to be considered when evaluating the role of volcanoes in global and regional climate, and serves to contextualize the well‐observed response to Pinatubo.
As one of the first global coupled climate models to simulate and predict category 4 and 5 (Saffir–Simpson scale) tropical cyclones (TCs) and their interannual variations, the High-Resolution Forecast-Oriented Low Ocean Resolution (HiFLOR) model at the Geophysical Fluid Dynamics Laboratory (GFDL) represents a novel source of insight on how the entire TC intensification distribution could be transformed due to climate change. In this study, three 70-year HiFLOR experiments are performed to identify the effects of climate change on TC intensity and intensification. For each of the experiments, sea surface temperature (SST) is nudged to different climatological targets and atmospheric radiative forcing is specified, allowing us to explore the sensitivity of TCs to these conditions.
First, a control experiment, which uses prescribed climatological ocean and radiative forcing based on observations during the years 1986-2005, is compared to two observational records and evaluated for its ability to capture the mean TC behavior during these years. The simulated intensification distributions as well as the percentage of TCs that become major hurricanes show similarities with observations. The control experiment is then compared to two 21st century experiments, in which the climatological SSTs from the control experiment are perturbed by multimodel projected SST anomalies and atmospheric radiative forcing from either 2016-2035 or 2081-2100 (RCP4.5 scenario). The frequency, intensity, and intensification distribution of TCs all shift to higher values as the 21st century progresses. HiFLOR’s unique response to climate change and fidelity in simulating the present climate lays the groundwork for future studies involving models of this type.
Mountain snowpack in the western United States provides a natural reservoir for cold season precipitation; variations in snowpack influence warm season water supply, wildfire risk, ecology, and industries like agriculture dependent on snow and downstream water availability. Efforts to understand snowpack variability have predominantly been focused on either weekly (weather) or decadal to centennial (climate variability and change) timescales. We focus on a timescale between these ranges by demonstrating that a global climate model suite can provide snowpack predictions 8 months in advance. The predictions from climate models outperform statistical methods from observations alone. Our results show that seasonal hydroclimate predictions are possible and highlight areas for future prediction system improvements.
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 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.
Balaji, V, E Maisonnave, Niki Zadeh, Bryan N Lawrence, Joachim Biercamp, Uwe Fladrich, G Aloisio, Rusty Benson, Arnaud Caubel, Jeffrey W Durachta, M-A Foujols, Grenville Lister, S Mocavero, Seth D Underwood, and Garrett Wright, January 2017: CPMIP: Measurements of Real Computational Performance of Earth System Models. Geoscientific Model Development, 10(1), DOI:10.5194/gmd-10-19-2017. Abstract
A climate model represents a multitude of processes on a variety of time and space scales; a canonical example of multi-physics multi-scale modeling. The underlying climate system is physically characterized by sensitive dependence on initial conditions, and natural stochastic variability, so very long integrations are needed to extract signals of climate change. Algorithms generally possess weak scaling and can be I/O and/or memory bound. Such weak-scaling, I/O and memory-bound multi-physics codes present particular challenges to computational performance.
Traditional metrics of computational efficiency such as performance counters and scaling curves do not tell us enough about real sustained performance from climate models on different machines. They also do not provide a satisfactory basis for comparative information across models.
We introduce a set of metrics that can be used for the study of computational performance of climate (and Earth System) models. These measures do not require specialized software or specific hardware counters, and should be accessible to anyone. They are independent of platform, and underlying parallel programming models. We show how these metrics can be used to measure actually attained performance of Earth system models on different machines, and identify the most fruitful areas of research and development for performance engineering.
We present results for these measures for a diverse suite of models from several modeling centres, and propose to use these measures as a basis for a CPMIP, a computational performance MIP.
This study explores the role of the stratosphere as a source of seasonal predictability of surface climate over Northern Hemisphere extra-tropics both in the observations and climate model predictions. A suite of numerical experiments, including climate simulations and retrospective forecasts, are set up to isolate the role of the stratosphere in seasonal predictive skill of extra-tropical near surface land temperature. We show that most of the lead-0 month spring predictive skill of land temperature over extra-tropics, particularly over northern Eurasia, stems from stratospheric initialization. We further reveal that this predictive skill of extra-tropical land temperature arises from skillful prediction of the Arctic Oscillation (AO). The dynamical connection between the stratosphere and troposphere is also demonstrated by the significant correlation between the stratospheric polar vortex and sea level pressure anomalies, as well as the migration of the stratospheric zonal wind anomalies to the lower troposphere.
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.
In 2014 and 2015, post-monsoon extremely severe cyclonic storms (ESCS)—defined by the WMO as tropical storms with lifetime maximum winds greater than 46 m s−1—were first observed over the Arabian Sea (ARB), causing widespread damage. However, it is unknown to what extent this abrupt increase in post-monsoon ESCSs can be linked to anthropogenic warming, natural variability, or stochastic behaviour. Here, using a suite of high-resolution global coupled model experiments that accurately simulate the climatological distribution of ESCSs, we show that anthropogenic forcing has likely increased the probability of late-season ECSCs occurring in the ARB since the preindustrial era. However, the specific timing of observed late-season ESCSs in 2014 and 2015 was likely due to stochastic processes. It is further shown that natural variability played a minimal role in the observed increase of ESCSs. Thus, continued anthropogenic forcing will further amplify the risk of cyclones in the ARB, with corresponding socio-economic implications.
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.
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–5; wind speed ≥58.1 m s−1), and landfalling TCs. Here we show that a 25-km resolution global coupled model (HiFLOR) developed at the Geophysical Fluid Dynamics Laboratory (GFDL) has improved skill in predicting the frequencies of major hurricanes and Category 4–5 hurricanes in the North Atlantic, and 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–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–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 improvements in large-scale parameters.
The impact of atmosphere and ocean horizontal resolution on the climatology of North American Monsoon Gulf of California (GoC) moisture surges is examined in a suite of global circulation models (CM2.1, FLOR, CM2.5, CM2.6, HiFLOR) developed at the Geophysical Fluid Dynamics Laboratory (GFDL). These models feature essentially the same physical parameterizations, but differ in horizontal resolution in either the atmosphere (≃200, 50 and 25 km) or the ocean (≃1°, 0.25°, 0.1°). Increasing horizontal atmospheric resolution from 200 km to 50 km results in a drastic improvement in the model’s capability of accurately simulating surge events. The climatological near-surface flow and moisture and precipitation anomalies associated with GoC surges are overall satisfactorily simulated in all higher-resolution models. The number of surge events agrees well with reanalyses but models tend to underestimate July-August surge-related precipitation and overestimate September surge-related rainfall in the southwestern United States. Large-scale controls supporting the development of GoC surges, such as tropical easterly waves (TEWs), tropical cyclones (TCs) and trans-Pacific Rossby wave trains (RWTs), are also well captured, although models tend to underestimate the TEW/TC magnitude and number. Near-surface GoC surge features and their large-scale forcings (TEWs, TCs, RWTs) do not appear to be substantially affected by a finer representation of the GoC at higher ocean resolution. However, the substantial reduction of the eastern Pacific warm sea surface temperature bias through flux adjustment in the FLOR model leads to an overall improvement of tropical-extratropical controls on GoC moisture surges and the seasonal cycle of precipitation in the southwestern United States.
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 whether, and the extent to which, an increase in atmospheric resolution in versions of the Geophysical Fluid Dynamics Laboratory (GFDL) High-Resolution Forecast-oriented Low Ocean Resolution Version of CM2.5 (FLOR) with 50 km and HiFLOR with 25 km improves the simulation of the El Niño Southern Oscillation-tropical cyclone (ENSO-TC) connections in the western North Pacific (WNP). HiFLOR simulates better ENSO-TC connections in the WNP including TC track density, genesis and landfall than FLOR in both long-term control experiments and sea surface temperature (SST)- and sea surface salinity (SSS)-restoring historical runs (1971-2012). Restoring experiments are performed with SSS and SST restored to observational estimates of climatological SSS and interannually-varying monthly SST. In the control experiments of HiFLOR, an improved simulation of the Walker circulation arising from more realistic SST and precipitation is largely responsible for its better performance in simulating ENSO-TC connections in the WNP. In the SST-restoring experiments of HiFLOR, more realistic Walker circulation and steering flow during El Niño/La Niña are responsible for the improved simulation of ENSO-TC connections in the WNP. The improved simulation of ENSO-TC connections with HiFLOR arises from a better representation of SST and better responses of environmental large-scale circulation to SST anomalies associated with El Niño/La Niña. A better representation of ENSO-TC connections in HiFLOR can benefit the seasonal forecasting of TC genesis, track and landfall, improve our understanding of the interannual variation of TC activity, and provide better projection of TC activity under climate change.
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.
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.
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.
