The authors assess the uptake, transport, and storage of oceanic anthropogenic carbon and heat over the period 1861–2005 in a new set of coupled carbon–climate Earth system models conducted for the ...fifth phase of the Coupled Model Intercomparison Project (CMIP5), with a particular focus on the Southern Ocean. Simulations show that the Southern Ocean south of 30°S, occupying 30% of global surface ocean area, accounts for 43% ± 3% (42 ± 5 Pg C) of anthropogenic CO₂ and 75% ± 22% (23 ± 9 × 1022J) of heat uptake by the ocean over the historical period. Northward transport out of the Southern Ocean is vigorous, reducing the storage to 33 ± 6 Pg anthropogenic carbon and 12 ± 7 × 1022J heat in the region. The CMIP5 models, as a class, tend to underestimate the observation-based global anthropogenic carbon storage but simulate trends in global ocean heat storage over the last 50 years within uncertainties of observation-based estimates. CMIP5 models suggest global and Southern Ocean CO₂ uptake have been largely unaffected by recent climate variability and change. Anthropogenic carbon and heat storage show a common broad-scale pattern of change, but ocean heat storage is more structured than ocean carbon storage. The results highlight the significance of the Southern Ocean for the global climate and as the region where models differ the most in representation of anthropogenic CO₂ and, in particular, heat uptake.
Recent research highlights the role of land surface processes in heat waves, droughts, and other extreme events. Here we use an earth system model (ESM) from the Geophysical Fluid Dynamics Laboratory ...(GFDL) to investigate the regional impacts of historical anthropogenic land useland cover change (LULCC) on combined extremes of temperature and humidity. A bivariate assessment allows us to consider aridity and moist enthalpy extremes, quantities central to human experience of near-surface climate conditions. We show that according to this model, conversion of forests to cropland has contributed to much of the upper central US and central Europe experiencing extreme hot, dry summers every 2-3 years instead of every 10 years. In the tropics, historical patterns of wood harvesting, shifting cultivation and regrowth of secondary vegetation have enhanced near surface moist enthalpy, leading to extensive increases in the occurrence of humid conditions throughout the tropics year round. These critical land use processes and practices are not included in many current generation land models, yet these results identify them as critical factors in the energy and water cycles of the midlatitudes and tropics.
The 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 in simulation characteristics of these two models are described. Because of differences in oceanic ventilation rates, ESM2M has a stronger biological carbon pump but weaker northward implied atmospheric CO₂ transport than ESM2G. The major advantages of ESM2G over ESM2M are improved representation of surface chlorophyll in the Atlantic and Indian Oceans and thermocline nutrients and oxygen in the North Pacific. Improved tree mortality parameters in ESM2G produced more realistic carbon accumulation in vegetation pools. The major advantages of ESM2M over ESM2G are reduced nutrient and oxygen biases in the southern and tropical oceans.
The physical climate formulation and simulation characteristics of two new global coupled carbon–climate Earth System Models, ESM2M and ESM2G, are described. These models demonstrate similar climate ...fidelity as the Geophysical Fluid Dynamics Laboratory’s previous Climate Model version 2.1 (CM2.1) while incorporating explicit and consistent carbon dynamics. The two models differ exclusively in the physical ocean component; ESM2M uses Modular Ocean Model version 4p1 with vertical pressure layers while ESM2G uses Generalized Ocean Layer Dynamics with a bulk mixed layer and interior isopycnal layers. Differences in the ocean mean state include the thermocline depth being relatively deep in ESM2M and relatively shallow in ESM2G compared to observations. The crucial role of ocean dynamics on climate variability is highlighted in El Niño–Southern Oscillation being overly strong in ESM2M and overly weak in ESM2G relative to observations. Thus, while ESM2G might better represent climate changes relating to total heat content variability given its lack of long-term drift, gyre circulation, and ventilation in the North Pacific, tropical Atlantic, and Indian Oceans, and depth structure in the overturning and abyssal flows, ESM2M might better represent climate changes relating to surface circulation given its superior surface temperature, salinity, and height patterns, tropical Pacific circulation and variability, and Southern Ocean dynamics. The overall assessment is that neither model is fundamentally superior to the other, and that both models achieve sufficient fidelity to allow meaningful climate and earth system modeling applications. This affords the ability to assess the role of ocean configuration on earth system interactions in the context of two state-of-the-art coupled carbon–climate models.
We document the configuration and emergent simulation features from the Geophysical Fluid Dynamics Laboratory (GFDL) OM4.0 ocean/sea ice model. OM4 serves as the ocean/sea ice component for the GFDL ...climate and Earth system models. It is also used for climate science research and is contributing to the Coupled Model Intercomparison Project version 6 Ocean Model Intercomparison Project. The ocean component of OM4 uses version 6 of the Modular Ocean Model and the sea ice component uses version 2 of the Sea Ice Simulator, which have identical horizontal grid layouts (Arakawa C‐grid). We follow the Coordinated Ocean‐sea ice Reference Experiments protocol to assess simulation quality across a broad suite of climate‐relevant features. We present results from two versions differing by horizontal grid spacing and physical parameterizations: OM4p5 has nominal 0.5° spacing and includes mesoscale eddy parameterizations and OM4p25 has nominal 0.25° spacing with no mesoscale eddy parameterization. Modular Ocean Model version 6 makes use of a vertical Lagrangian‐remap algorithm that enables general vertical coordinates. We show that use of a hybrid depth‐isopycnal coordinate reduces the middepth ocean warming drift commonly found in pure z* vertical coordinate ocean models. To test the need for the mesoscale eddy parameterization used in OM4p5, we examine the results from a simulation that removes the eddy parameterization. The water mass structure and model drift are physically degraded relative to OM4p5, thus supporting the key role for a mesoscale closure at this resolution.
Key Points
Documentation is provided for a new generation of NOAA‐GFDL CMIP6/OMIP ocean ice climate models
Dynamical core and physical parameterizations are described and key features of interannual CORE simulations are assessed
Using hybrid vertical coordinates reduces spurious ocean heat drift
Using simulations performed with 18 coupled atmosphere–ocean global climate models from phase 5 of the Coupled Model Intercomparison Project (CMIP5), projections of the Northern Hemisphere snowfall ...under the representative concentration pathway (RCP4.5) scenario are analyzed for the period 2006–2100. These models perform well in simulating twentieth-century snowfall, although there is a positive bias in many regions. Annual snowfall is projected to decrease across much of the Northern Hemisphere during the twenty-first century, with increases projected at higher latitudes. On a seasonal basis, the transition zone between negative and positive snowfall trends corresponds approximately to the −10°C isotherm of the late twentieth-century mean surface air temperature, such that positive trends prevail in winter over large regions of Eurasia and North America. Redistributions of snowfall throughout the entire snow season are projected to occur—even in locations where there is little change in annual snowfall. Changes in the fraction of precipitation falling as snow contribute to decreases in snowfall across most Northern Hemisphere regions, while changes in total precipitation typically contribute to increases in snowfall. A signal-to-noise analysis reveals that the projected changes in snowfall, based on the RCP4.5 scenario, are likely to become apparent during the twenty-first century for most locations in the Northern Hemisphere. The snowfall signal emerges more slowly than the temperature signal, suggesting that changes in snowfall are not likely to be early indicators of regional climate change.
This contribution describes the ocean biogeochemical component of the Geophysical Fluid Dynamics Laboratory's Earth System Model 4.1 (GFDL‐ESM4.1), assesses GFDL‐ESM4.1's capacity to capture observed ...ocean biogeochemical patterns, and documents its response to increasing atmospheric CO2. Notable differences relative to the previous generation of GFDL ESM's include enhanced resolution of plankton food web dynamics, refined particle remineralization, and a larger number of exchanges of nutrients across Earth system components. During model spin‐up, the carbon drift rapidly fell below the 10 Pg C per century equilibration criterion established by the Coupled Climate‐Carbon Cycle Model Intercomparison Project (C4MIP). Simulations robustly captured large‐scale observed nutrient distributions, plankton dynamics, and characteristics of the biological pump. The model overexpressed phosphate limitation and open ocean hypoxia in some areas but still yielded realistic surface and deep carbon system properties, including cumulative carbon uptake since preindustrial times and over the last decades that is consistent with observation‐based estimates. The model's response to the direct and radiative effects of a 200% atmospheric CO2 increase from preindustrial conditions (i.e., years 101–120 of a 1% CO2 yr−1 simulation) included (a) a weakened, shoaling organic carbon pump leading to a 38% reduction in the sinking flux at 2,000 m; (b) a two‐thirds reduction in the calcium carbonate pump that nonetheless generated only weak calcite compensation on century time‐scales; and, in contrast to previous GFDL ESMs, (c) a moderate reduction in global net primary production that was amplified at higher trophic levels. We conclude with a discussion of model limitations and priority developments.
Plain Language Summary
This paper describes and evaluates the ocean biogeochemical component of the Geophysical Fluid Dynamics Laboratory's Earth System Model 4.1 (GFDL‐ESM4.1). GFDL‐ESM4.1 was developed to study the past, present, and future evolution of the Earth system under scenarios for natural and anthropogenic drivers of Earth system change, including greenhouse gases and aerosols. The response of the ocean's vast carbon and heat reservoirs to accumulating greenhouse gases greatly reduces their atmospheric and terrestrial impacts, but also puts ocean environments and the marine resources they support at risk. Relative to previous models, GFDL‐ESM4.1 improves the representation of (a) ocean food webs connecting plankton and fish; (b) biological processes influencing the sequestration of carbon in the deep ocean; and (c) land‐atmosphere‐ocean nutrient exchanges. While simulations have biases, they capture many critical aspects of the global ocean carbon cycle and ocean ecosystem, including the observed uptake of anthropogenic carbon over the last ~150 yr. Projections suggest that continued CO2 increases could significantly decrease ocean productivity and the ocean's capacity to sequester atmospheric carbon.
Key Points
Enhanced plankton food webs, remineralization, and Earth system linkages yield skillful global carbon cycle and ecosystem simulations
Ocean productivity estimates improved, but phosphate limitation and hypoxia overestimated in some areas
High CO2 and associated warming substantially reduce the biological pump and ocean productivity across trophic levels
We compare equilibrium climate sensitivity (ECS) estimates from pairs of long (≥ 800‐year) control and abruptly quadrupled CO2 simulations with shorter (150, 300 year) coupled atmosphere‐ocean ...simulations and Slab Ocean Models (SOM). Consistent with previous work, ECS estimates from shorter coupled simulations based on annual averages for years 1‐150 underestimate those from SOM (‐8% ± 13%) and long (‐14% ± 8%) simulations. Analysis of only years 21‐150 improved agreement with SOM (‐2% ± 14%) and long (‐8% ± 10%) estimates. Use of pentadal averages for years 51‐150 results in improved agreement with long simulations (‐4% ± 11%). While ECS estimates from current generation US models based on SOM and coupled annual averages of years 1‐150 range from 2.6°C to 5.3°C, estimates based longer simulations of the same models range from 3.2°C to 7.0°C. Such variations between methods argues for caution in comparison and interpretation of ECS estimates across models.
Previous studies have demonstrated the importance of enhanced vegetation growth under future elevated atmospheric CO ₂ for 21st century climate warming. Surprisingly no study has completed an ...analogous assessment for the historical period, during which emissions of greenhouse gases increased rapidly and land-use changes (LUC) dramatically altered terrestrial carbon sources and sinks. Using the Geophysical Fluid Dynamics Laboratory comprehensive Earth System Model ESM2G and a reconstruction of the LUC, we estimate that enhanced vegetation growth has lowered the historical atmospheric CO ₂ concentration by 85 ppm, avoiding an additional 0.31 ± 0.06 °C warming. We demonstrate that without enhanced vegetation growth the total residual terrestrial carbon flux (i.e., the net land flux minus LUC flux) would be a source of 65–82 Gt of carbon (GtC) to atmosphere instead of the historical residual carbon sink of 186–192 GtC, a carbon saving of 251–274 GtC.