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
The GFDL CM3 Coupled Climate Model Griffies, Stephen M.; Winton, Michael; Donner, Leo J. ...
Journal of climate,
07/2011, Volume:
24, Issue:
13
Journal Article
Peer reviewed
Open access
This paper documents time mean simulation characteristics from the ocean and sea ice components in a new coupled climate model developed at the NOAA Geophysical Fluid Dynamics Laboratory (GFDL). The ...GFDL Climate Model version 3 (CM3) is formulated with effectively the same ocean and sea ice components as the earlier CM2.1 yet with extensive developments made to the atmosphere and land model components. Both CM2.1 and CM3 show stable mean climate indices, such as large-scale circulation and sea surface temperatures (SSTs). There are notable improvements in the CM3 climate simulation relative to CM2.1, including a modified SST bias pattern and reduced biases in the Arctic sea ice cover. The authors anticipate SST differences between CM2.1 and CM3 in lower latitudes through analysis of the atmospheric fluxes at the ocean surface in corresponding Atmospheric Model Intercomparison Project (AMIP) simulations. In contrast, SST changes in the high latitudes are dominated by ocean and sea ice effects absent in AMIP simulations. The ocean interior simulation in CM3 is generally warmer than in CM2.1, which adversely impacts the interior biases.
The influence of changing ocean currents on climate change is evaluated by comparing an earth system model’s response to increased CO₂ with and without an ocean circulation response. Inhibiting the ...ocean circulation response, by specifying a seasonally varying preindustrial climatology of currents, has a much larger influence on the heat storage pattern than on the carbon storage pattern. The heat storage pattern without circulation changes resembles carbon storage (either with or without circulation changes) more than it resembles the heat storage when currents are allowed to respond. This is shown to be due to the larger magnitude of the redistribution transport—the change in transport due to circulation anomalies acting on control climate gradients—for heat than for carbon. The net ocean heat and carbon uptake are slightly reduced when currents are allowed to respond. Hence, ocean circulation changes potentially act to warm the surface climate. However, the impact of the reduced carbon uptake on radiative forcing is estimated to be small while the redistribution heat transport shifts ocean heat uptake from low to high latitudes, increasing its cooling power. Consequently, global surface warming is significantly reduced by circulation changes. Circulation changes also shift the pattern of warming from broad Northern Hemisphere amplification to a more structured pattern with reduced warming at subpolar latitudes in both hemispheres and enhanced warming near the equator.
•Phase II of the Coordinated Ocean-ice Reference Experiments (CORE-II) is introduced.•Solutions from CORE-II simulations from eighteen participating models are presented.•Mean states in the North ...Atlantic with a focus on AMOC are examined.•The North Atlantic solutions differ substantially among the models.•Many factors, including parameterization choices, contribute to these differences.
Simulation characteristics from eighteen global ocean–sea-ice coupled models are presented with a focus on the mean Atlantic meridional overturning circulation (AMOC) and other related fields in the North Atlantic. These experiments use inter-annually varying atmospheric forcing data sets for the 60-year period from 1948 to 2007 and are performed as contributions to the second phase of the Coordinated Ocean-ice Reference Experiments (CORE-II). The protocol for conducting such CORE-II experiments is summarized. Despite using the same atmospheric forcing, the solutions show significant differences. As most models also differ from available observations, biases in the Labrador Sea region in upper-ocean potential temperature and salinity distributions, mixed layer depths, and sea-ice cover are identified as contributors to differences in AMOC. These differences in the solutions do not suggest an obvious grouping of the models based on their ocean model lineage, their vertical coordinate representations, or surface salinity restoring strengths. Thus, the solution differences among the models are attributed primarily to use of different subgrid scale parameterizations and parameter choices as well as to differences in vertical and horizontal grid resolutions in the ocean models. Use of a wide variety of sea-ice models with diverse snow and sea-ice albedo treatments also contributes to these differences. Based on the diagnostics considered, the majority of the models appear suitable for use in studies involving the North Atlantic, but some models require dedicated development effort.
•Inter-annual to decadal variability in AMOC from CORE-II simulations is presented.•AMOC variability shows three stages, with maximum transports in mid- to late-1990s.•North Atlantic temporal ...variability features are in good agreement among simulations.•Such agreements suggest variability is dictated by the atmospheric data sets.•Simulations differ in spatial structures of variability due to ocean dynamics.
Simulated inter-annual to decadal variability and trends in the North Atlantic for the 1958–2007 period from twenty global ocean – sea-ice coupled models are presented. These simulations are performed as contributions to the second phase of the Coordinated Ocean-ice Reference Experiments (CORE-II). The study is Part II of our companion paper (Danabasoglu et al., 2014) which documented the mean states in the North Atlantic from the same models. A major focus of the present study is the representation of Atlantic meridional overturning circulation (AMOC) variability in the participating models. Relationships between AMOC variability and those of some other related variables, such as subpolar mixed layer depths, the North Atlantic Oscillation (NAO), and the Labrador Sea upper-ocean hydrographic properties, are also investigated. In general, AMOC variability shows three distinct stages. During the first stage that lasts until the mid- to late-1970s, AMOC is relatively steady, remaining lower than its long-term (1958–2007) mean. Thereafter, AMOC intensifies with maximum transports achieved in the mid- to late-1990s. This enhancement is then followed by a weakening trend until the end of our integration period. This sequence of low frequency AMOC variability is consistent with previous studies. Regarding strengthening of AMOC between about the mid-1970s and the mid-1990s, our results support a previously identified variability mechanism where AMOC intensification is connected to increased deep water formation in the subpolar North Atlantic, driven by NAO-related surface fluxes. The simulations tend to show general agreement in their temporal representations of, for example, AMOC, sea surface temperature (SST), and subpolar mixed layer depth variabilities. In particular, the observed variability of the North Atlantic SSTs is captured well by all models. These findings indicate that simulated variability and trends are primarily dictated by the atmospheric datasets which include the influence of ocean dynamics from nature superimposed onto anthropogenic effects. Despite these general agreements, there are many differences among the model solutions, particularly in the spatial structures of variability patterns. For example, the location of the maximum AMOC variability differs among the models between Northern and Southern Hemispheres.
Coordinated Ocean-ice Reference Experiments (COREs) are presented as a tool to explore the behaviour of global ocean-ice models under forcing from a common atmospheric dataset. We highlight issues ...arising when designing coupled global ocean and sea ice experiments, such as difficulties formulating a consistent forcing methodology and experimental protocol. Particular focus is given to the hydrological forcing, the details of which are key to realizing simulations with stable meridional overturning circulations.
The atmospheric forcing from Large, W., Yeager, S., 2004. Diurnal to decadal global forcing for ocean and sea-ice models: the data sets and flux climatologies. NCAR Technical Note: NCAR/TN-460+STR. CGD Division of the National Center for Atmospheric Research was developed for coupled-ocean and sea ice models. We found it to be suitable for our purposes, even though its evaluation originally focussed more on the ocean than on the sea-ice. Simulations with this atmospheric forcing are presented from seven global ocean-ice models using the CORE-I design (repeating annual cycle of atmospheric forcing for 500 years). These simulations test the hypothesis that global ocean-ice models run under the same atmospheric state produce qualitatively similar simulations. The validity of this hypothesis is shown to depend on the chosen diagnostic. The CORE simulations provide feedback to the fidelity of the atmospheric forcing and model configuration, with identification of biases promoting avenues for forcing dataset and/or model development.
Despite slow rates of ocean mixing, observational and modeling studies suggest that buoyancy is redistributed to all depths of the ocean on surprisingly short interannual to decadal time scales. The ...mechanisms responsible for this redistribution remain poorly understood. This work uses an Earth system model to evaluate the global steady-state ocean buoyancy (and related steric sea level) budget, its interannual variability, and its transient response to a doubling of CO₂ over 70 years, with a focus on the deep ocean. At steady state, the simple view of vertical advective–diffusive balance for the deep ocean holds at low to midlatitudes. At higher latitudes, the balance depends on a myriad of additional terms, namely mesoscale and submesoscale advection, convection and overflows from marginal seas, and terms related to the nonlinear equation of state. These high-latitude processes rapidly communicate anomalies in surface buoyancy forcing to the deep ocean locally; the deep, high-latitude changes then influence the large-scale advection of buoyancy to create transient deep buoyancy anomalies at lower latitudes. Following a doubling of atmospheric carbon dioxide concentrations, the high-latitude buoyancy sinks are suppressed by a slowdown in convection and reduced dense water formation. This change is accompanied by a slowing of both upper and lower cells of the global meridional overturning circulation, reducing the supply of dense water to low latitudes beneath the pycnocline and the commensurate flow of light waters to high latitudes above the pycnocline. By this mechanism, changes in high-latitude buoyancy are communicated to the global deep ocean on relatively fast advective time scales.
•We focus on ACC and Southern Ocean MOC during 1958–2007 in 17 CORE-II forced models.•Most CORE-II simulations are close to eddy saturation.•Most CORE-II simulations are far from showing signs of ...eddy compensation.•Constant in time or space k results in poor representation of mesoscale eddy effects.•MOC has larger sensitivity than ACC transport even in eddy saturated state.
In the framework of the second phase of the Coordinated Ocean-ice Reference Experiments (CORE-II), we present an analysis of the representation of the Antarctic Circumpolar Current (ACC) and Southern Ocean meridional overturning circulation (MOC) in a suite of seventeen global ocean–sea ice models. We focus on the mean, variability and trends of both the ACC and MOC over the 1958–2007 period, and discuss their relationship with the surface forcing. We aim to quantify the degree of eddy saturation and eddy compensation in the models participating in CORE-II, and compare our results with available observations, previous fine-resolution numerical studies and theoretical constraints. Most models show weak ACC transport sensitivity to changes in forcing during the past five decades, and they can be considered to be in an eddy saturated regime. Larger contrasts arise when considering MOC trends, with a majority of models exhibiting significant strengthening of the MOC during the late 20th and early 21st century. Only a few models show a relatively small sensitivity to forcing changes, responding with an intensified eddy-induced circulation that provides some degree of eddy compensation, while still showing considerable decadal trends. Both ACC and MOC interannual variabilities are largely controlled by the Southern Annular Mode (SAM). Based on these results, models are clustered into two groups. Models with constant or two-dimensional (horizontal) specification of the eddy-induced advection coefficient κ show larger ocean interior decadal trends, larger ACC transport decadal trends and no eddy compensation in the MOC. Eddy-permitting models or models with a three-dimensional time varying κ show smaller changes in isopycnal slopes and associated ACC trends, and partial eddy compensation. As previously argued, a constant in time or space κ is responsible for a poor representation of mesoscale eddy effects and cannot properly simulate the sensitivity of the ACC and MOC to changing surface forcing. Evidence is given for a larger sensitivity of the MOC as compared to the ACC transport, even when approaching eddy saturation. Future process studies designed for disentangling the role of momentum and buoyancy forcing in driving the ACC and MOC are proposed.
Oceanic heat uptake (OHU) is a significant source of uncertainty in both the transient and equilibrium responses to increasing the planetary radiative forcing. OHU differs among climate models and is ...related in part to their representation of vertical and lateral mixing. This study examines the role of ocean model formulation—specifically the choice of the vertical coordinate and the strength of the background diapycnal diffusivity Kd
—in the millennial-scale near-equilibrium climate response to a quadrupling of atmospheric CO2. Using two fully coupled Earth system models (ESMs) with nearly identical atmosphere, land, sea ice, and biogeochemical components, it is possible to independently configure their ocean model components with different formulations and produce similar near-equilibrium climate responses. The SST responses are similar between the two models (r² = 0.75, global average ∼4.3°C) despite their initial preindustrial climate mean states differing by 0.4°C globally. The surface and interior responses of temperature and salinity are also similar between the two models. However, the Atlantic meridional overturning circulation (AMOC) responses are different between the two models, and the associated differences in ventilation and deep-water formation have an impact on the accumulation of dissolved inorganic carbon in the ocean interior. A parameter sensitivity analysis demonstrates that increasing the amount of Kd
produces very different near-equilibrium climate responses within a given model. These results suggest that the impact of the ocean vertical coordinate on the climate response is small relative to the representation of subgrid-scale mixing.