Ageostrophic baroclinic instabilities develop within the surface mixed layer of the ocean at horizontal fronts and efficiently restratify the upper ocean. In this paper a parameterization for the restratification driven by finite-amplitude baroclinic instabilities of the mixed layer is proposed in terms of an overturning streamfunction that tilts isopycnals from the vertical to the horizontal. The streamfunction is proportional to the product of the horizontal density gradient, the mixed layer depth squared, and the inertial period. Hence restratification proceeds faster at strong fronts in deep mixed layers with a weak latitude dependence. In this paper the parameterization is theoretically motivated, confirmed to perform well for a wide range of mixed layer depths, rotation rates, and vertical and horizontal stratifications. It is shown to be superior to alternative extant parameterizations of baroclinic instability for the problem of mixed layer restratification. Two companion papers discuss the numerical implementation and the climate impacts of this parameterization.

No abstract

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.

The Modeling Eddies in the Southern Ocean (MESO) project uses numerical sensitivity studies to examine the role played by Southern Ocean winds and eddies in determining the density structure of the global ocean and the magnitude and structure of the global overturning circulation. A hemispheric isopycnal-coordinate ocean model (which avoids numerical diapycnal diffusion) with realistic geometry is run with idealized forcing at a range of resolutions from coarse (2°) to eddy-permitting ( 1 ⁄6°). A comparison of coarse resolutions with fine resolutions indicates that explicit eddies affect both the structure of the overturning and the response of the overturning to wind stress changes. While the presence of resolved eddies does not greatly affect the prevailing qualitative picture of the ocean circulation, it alters the overturning cells involving the Southern Ocean transformation of dense deep waters and light waters of subtropical origin into intermediate waters. With resolved eddies, the surface-to-intermediate water cell extends farther southward by hundreds of kilometers and the deep-to-intermediate cell draws on comparatively lighter deep waters. The overturning response to changes in the winds is also sensitive to the presence of eddies. In noneddying simulations, changing the Ekman transport produces comparable changes in the overturning, much of it involving transformation of deep waters and resembling the mean circulation. In the eddypermitting simulations, a significant fraction of the Ekman transport changes are compensated by eddyinduced transport drawing from lighter waters than does the mean overturning. This significant difference calls into question the ability of coarse-resolution ocean models to accurately capture the impact of changes in the Southern Ocean on the global ocean circulation.

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 2 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.

No abstract

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.

Abstract. The Ocean Model Intercomparison Project (OMIP) is an endorsed project in the Coupled Model Intercomparison Project Phase 6 (CMIP6). OMIP addresses CMIP6 science questions, investigating the origins and consequences of systematic model biases. It does so by providing a framework for evaluating (including assessment of systematic biases), understanding, and improving ocean, sea-ice, tracer, and biogeochemical components of climate and earth system models contributing to CMIP6. Among the WCRP Grand Challenges in climate science (GCs), OMIP primarily contributes to the regional sea level change and near-term (climate/decadal) prediction GCs.OMIP provides (a) an experimental protocol for global ocean/sea-ice models run with a prescribed atmospheric forcing; and (b) a protocol for ocean diagnostics to be saved as part of CMIP6. We focus here on the physical component of OMIP, with a companion paper (Orr et al., 2016) detailing methods for the inert chemistry and interactive biogeochemistry. The physical portion of the OMIP experimental protocol follows the interannual Coordinated Ocean-ice Reference Experiments (CORE-II). Since 2009, CORE-I (Normal Year Forcing) and CORE-II (Interannual Forcing) have become the standard methods to evaluate global ocean/sea-ice simulations and to examine mechanisms for forced ocean climate variability. The OMIP diagnostic protocol is relevant for any ocean model component of CMIP6, including the DECK (Diagnostic, Evaluation and Characterization of Klima experiments), historical simulations, FAFMIP (Flux Anomaly Forced MIP), C4MIP (Coupled Carbon Cycle Climate MIP), DAMIP (Detection and Attribution MIP), DCPP (Decadal Climate Prediction Project), ScenarioMIP, HighResMIP (High Resolution MIP), as well as the ocean/sea-ice OMIP simulations.