Matthew J. Harrison

GFDL's CM2 global coupled climate models. Part I: Formulation and simulation characteristics

Abstract The formulation and simulation characteristics of two new global coupled climate models developed at NOAAs Geophysical Fluid Dynamics Laboratory (GFDL) are described. The models were designed to simulate atmospheric and oceanic climate and variability from the diurnal time scale through multicentury climate change, given our computational constraints. In particular, an important goal was to use the same model for both experimental seasonal to interannual forecasting and the study of multicentury global climate change, and this goal has been achieved. Two versions of the coupled model are described, called CM2.0 and CM2.1. The versions differ primarily in the dynamical core used in the atmospheric component, along with the cloud tuning and some details of the land and ocean components. For both coupled models, the resolution of the land and atmospheric components is 2° latitude × 2.5° longitude; the atmospheric model has 24 vertical levels. The ocean resolution is 1° in latitude and longitude, wi...

Weakening of tropical Pacific atmospheric circulation due to anthropogenic forcing

Since the mid-nineteenth century the Earths surface has warmed, and models indicate that human activities have caused part of the warming by altering the radiative balance of the atmosphere. Simple theories suggest that global warming will reduce the strength of the mean tropical atmospheric circulation. An important aspect of this tropical circulation is a large-scale zonal (east–west) overturning of air across the equatorial Pacific Ocean—driven by convection to the west and subsidence to the east—known as the Walker circulation. Here we explore changes in tropical Pacific circulation since the mid-nineteenth century using observations and a suite of global climate model experiments. Observed Indo-Pacific sea level pressure reveals a weakening of the Walker circulation. The size of this trend is consistent with theoretical predictions, is accurately reproduced by climate model simulations and, within the climate models, is largely due to anthropogenic forcing. The climate model indicates that the weakened surface winds have altered the thermal structure and circulation of the tropical Pacific Ocean. These results support model projections of further weakening of tropical atmospheric circulation during the twenty-first century.

GFDL’s ESM2 Global Coupled Climate–Carbon Earth System Models. Part I: Physical Formulation and Baseline Simulation Characteristics

AbstractThe authors describe carbon system formulation and simulation characteristics of two new global coupled carbon–climate Earth System Models (ESM), ESM2M and ESM2G. These models demonstrate good climate fidelity as described in part I of this study while incorporating explicit and consistent carbon dynamics. The two models differ almost exclusively in the physical ocean component; ESM2M uses the Modular Ocean Model version 4.1 with vertical pressure layers, whereas ESM2G uses generalized ocean layer dynamics with a bulk mixed layer and interior isopycnal layers. On land, both ESMs include a revised land model to simulate competitive vegetation distributions and functioning, including carbon cycling among vegetation, soil, and atmosphere. In the ocean, both models include new biogeochemical algorithms including phytoplankton functional group dynamics with flexible stoichiometry. Preindustrial simulations are spun up to give stable, realistic carbon cycle means and variability. Significant differences...

GFDL's CM2 Global Coupled Climate Models. Part II: The Baseline Ocean Simulation

The current generation of coupled climate models run at the Geophysical Fluid Dynamics Laboratory (GFDL) as part of the Climate Change Science Program contains ocean components that differ in almost every respect from those contained in previous generations of GFDL climate models. This paper summarizes the new physical features of the models and examines the simulations that they produce. Of the two new coupled climate model versions 2.1 (CM2.1) and 2.0 (CM2.0), the CM2.1 model represents a major improvement over CM2.0 in most of the major oceanic features examined, with strikingly lower drifts in hydrographic fields such as temperature and salinity, more realistic ventilation of the deep ocean, and currents that are closer to their observed values. Regional analysis of the differences between the models highlights the importance of wind stress in determining the circulation, particularly in the Southern Ocean. At present, major errors in both models are associated with Northern Hemisphere Mode Waters and outflows from overflows, particularly the Mediterranean Sea and Red Sea.

System design and evaluation of coupled ensemble data assimilation for global oceanic climate studies

Abstract A fully coupled data assimilation (CDA) system, consisting of an ensemble filter applied to the Geophysical Fluid Dynamics Laboratory’s global fully coupled climate model (CM2), has been developed to facilitate the detection and prediction of seasonal-to-multidecadal climate variability and climate trends. The assimilation provides a self-consistent, temporally continuous estimate of the coupled model state and its uncertainty, in the form of discrete ensemble members, which can be used directly to initialize probabilistic climate forecasts. Here, the CDA is evaluated using a series of perfect model experiments, in which a particular twentieth-century simulation—with temporally varying greenhouse gas and natural aerosol radiative forcings—serves as a “truth” from which observations are drawn, according to the actual ocean observing network for the twentieth century. These observations are then assimilated into a coupled model ensemble that is subjected only to preindustrial forcings. By examining...

Impacts of Shortwave Penetration Depth on Large-Scale Ocean Circulation and Heat Transport

The impact of changes in shortwave radiation penetration depth on the global ocean circulation and heat transport is studied using the GFDL Modular Ocean Model (MOM4) with two independent parameterizations that use ocean color to estimate the penetration depth of shortwave radiation. Ten to eighteen percent increases in the depth of 1% downwelling surface irradiance levels results in an increase in mixed layer depths of 3–20 m in the subtropical and tropical regions with no change at higher latitudes. While 1D models have predicted that sea surface temperatures at the equator would decrease with deeper penetration of solar irradiance, this study shows a warming, resulting in a 10% decrease in the required restoring heat flux needed to maintain climatological sea surface temperatures in the eastern equatorial Atlantic and Pacific Oceans. The decrease in the restoring heat flux is attributed to a slowdown in heat transport (5%) from the Tropics and an increase in the temperature of submixed layer waters being transported into the equatorial regions. Calculations were made using a simple relationship between mixed layer depth and meridional mass transport. When compared with model diagnostics, these calculations suggest that the slowdown in heat transport is primarily due to off-equatorial increases in mixed layer depths. At higher latitudes (5°–40°), higher restoring heat fluxes are needed to maintain sea surface temperatures because of deeper mixed layers and an increase in storage of heat below the mixed layer. This study offers a way to evaluate the changes in irradiance penetration depths in coupled ocean–atmosphere GCMs and the potential effect that large-scale changes in chlorophyll a concentrations will have on ocean circulation.

Enhanced warming of the Northwest Atlantic Ocean under climate change

The Intergovernmental Panel on Climate Change (IPCC) fifth assessment of projected global and regional ocean temperature change is based on global climate models that have coarse (∼100 km) ocean and atmosphere resolutions. In the Northwest Atlantic, the ensemble of global climate models has a warm bias in sea surface temperature due to a misrepresentation of the Gulf Stream position; thus, existing climate change projections are based on unrealistic regional ocean circulation. Here we compare simulations and an atmospheric CO2 doubling response from four global climate models of varying ocean and atmosphere resolution. We find that the highest resolution climate model (∼10 km ocean, ∼50 km atmosphere) resolves Northwest Atlantic circulation and water mass distribution most accurately. The CO2 doubling response from this model shows that upper-ocean (0–300 m) temperature in the Northwest Atlantic Shelf warms at a rate nearly twice as fast as the coarser models and nearly three times faster than the global average. This enhanced warming is accompanied by an increase in salinity due to a change in water mass distribution that is related to a retreat of the Labrador Current and a northerly shift of the Gulf Stream. Both observations and the climate model demonstrate a robust relationship between a weakening Atlantic Meridional Overturning Circulation (AMOC) and an increase in the proportion of Warm-Temperate Slope Water entering the Northwest Atlantic Shelf. Therefore, prior climate change projections for the Northwest Atlantic may be far too conservative. These results point to the need to improve simulations of basin and regional-scale ocean circulation.

The Equatorial Thermocline Outcropping—A Seasonal Control on the Tropical Pacific Ocean–Atmosphere Instability Strength

Abstract One of the major factors determining the strength and extent of ENSO events is the instability state of the equatorial Pacific coupled ocean–atmosphere system and its seasonal variations. This study analyzes the coupled instability in a hybrid coupled model of the Indo–Pacific region, using the adjoint method for sensitivity studies. It is found that the seasonal changes in the ocean–atmosphere instability strength in the model used here are related to the outcropping of the thermocline in the east equatorial Pacific. From July to December, when the thermocline outcrops over a wide area in the east Pacific, there is a strong surface–thermocline connection and anomalies that arrive as Kelvin waves from the west along the thermocline can reach the surface and affect the SST and thus the coupled system. Conversely, from February to June, when the thermocline outcropping is minimal, the surface decouples from the thermocline and temperature anomalies in the thermocline depth range do not affect the s...

Initialization of an enso forecast system using a parallelized ensemble filter

Abstract As a first step toward coupled ocean–atmosphere data assimilation, a parallelized ensemble filter is implemented in a new stochastic hybrid coupled model. The model consists of a global version of the GFDL Modular Ocean Model Version 4 (MOM4), coupled to a statistical atmosphere based on a regression of National Centers for Environmental Prediction (NCEP) reanalysis surface wind stress, heat, and water flux anomalies onto analyzed tropical Pacific SST anomalies from 1979 to 2002. The residual part of the NCEP fluxes not captured by the regression is then treated as stochastic forcing, with different ensemble members feeling the residual fluxes from different years. The model provides a convenient test bed for coupled data assimilation, as well as a prototype for representing uncertainties in the surface forcing. A parallel ensemble adjustment Kalman filter (EAKF) has been designed and implemented in the hybrid model, using a local least squares framework. Comparison experiments demonstrate that t...

Sensitivity of Hybrid ENSO Models to Unresolved Atmospheric Variability

Abstract A common practice in the design of forecast models for ENSO is to couple ocean general circulation models to simple atmospheric models. Therefore, by construction these models (known as hybrid ENSO models) do not resolve various kinds of atmospheric variability [e.g., the Madden–Julian oscillation (MJO) and westerly wind bursts] that are often regarded as “unwanted noise.” In this work the sensitivity of three hybrid ENSO models to this unresolved atmospheric variability is studied. The hybrid coupled models were tuned to be asymptotically stable and the magnitude, and spatial and temporal structure of the unresolved variability was extracted from observations. The results suggest that this neglected variability can add an important piece of realism and forecast skill to the hybrid models. The models were found to respond linearly to the low-frequency part of the neglected atmospheric variability, in agreement with previous findings with intermediate models. While the wind anomalies associated wi...