Abstract

The role of oceanic background diapycnal diffusion for the equilibrium climate state is investigated in the global coupled climate model CM2G. Special emphasis is put on the oceanic meridional overturning and heat transport. Six runs with the model, differing only by their value of the background diffusivity, are run to steady state and the statistically steady integrations are compared. The diffusivity changes have large-scale impacts on many aspects of the climate system. Two examples are the volume-mean potential temperature, which increases by 3.6°C between the least and most diffusive runs, and the Antarctic sea ice extent, which decreases rapidly as the diffusivity increases. The overturning scaling with diffusivity is found to agree rather well with classical theoretical results for the upper but not for the lower cell. An alternative empirical scaling with the mixing energy is found to give good results for both cells. The oceanic meridional heat transport increases strongly with the diffusivity, an increase that can only partly be explained by increases in the meridional overturning. The increasing poleward oceanic heat transport is accompanied by a decrease in its atmospheric counterpart, which keeps the increase in the planetary energy transport small compared to that in the ocean.

Abstract

An isopycnal model of the North Pacific is used to demonstrate that the seasonal cycle of heating and cooling and the resulting mixed layer depth entrainment and detrainment cycle play a role in the propagation of wind-driven Rossby waves. The model is forced by realistic winds and seasonal heat flux to examine the interaction of nearly annual wind-driven Rossby waves with the seasonal mixed layer cycle. Comparison among four model runs, one adiabatic (without diapycnal mixing or explicit mixed layer dynamics), one diabatic (with diapycnal mixing and explicit mixed layer dynamics), one with the seasonal cycle of heating only, and one with only variable winds suggests that mixed layer entrainment changes the structure of the response substantially, particularly at midlatitudes. Specifically, the mixed layer seasonal cycle works against Ekman pumping in the forcing of first-mode Rossby waves between 17° and 28°N. South of there the mixed layer seasonal cycle has little influence on the Rossby waves, while in the north, seasonal Rossby waves do not propagate. To examine the first baroclinic mode response in detail, a modal decomposition of the numerical model output is done. In addition, a comparison of the forcing by diapycnal pumping and Ekman pumping is done by a projection of Ekman pumping and diapycnal velocities on to the quasigeostrophic potential vorticity equation for each vertical mode. The first baroclinic mode's forcing is split between Ekman pumping and diapycnal velocity at midlatitudes, providing an explanation for the changes in the response when a seasonal mixed layer response is included. This is confirmed by doing a comparison of the modal decomposition in the four runs described above, and by calculation of the first baroclinic mode Rossby wave response using the one-dimensional Rossby wave equation.

Abstract

The vertical structure of ocean eddies is generally surface-intensified, commonly attributed to the dominant baroclinic modes arising from the boundary conditions (BCs). Conventional BC considerations mostly focus on either flat- or rough-bottom conditions. The impact of surface buoyancy anomalies—often represented by surface potential vorticity (PV) anomalies—has not been fully explored. Here, we study the role of the surface PV in setting the vertical distribution of eddy kinetic energy (EKE) in an idealized adiabatic ocean model driven by wind stress. The simulated EKE profile in the extratropical ocean tends to peak at the surface and have an e-folding depth typically smaller than half of the ocean depth. This vertical structure can be reasonably represented by a single surface quasigeostrophic (SQG) mode at the energy-containing scale resulting from the large-scale PV structure. Due to isopycnal outcropping and interior PV homogenization, the surface meridional PV gradient is substantially stronger than the interior PV gradient, yielding surface-trapped baroclinically unstable modes with horizontal scales comparable to or smaller than the deformation radius. These surface-trapped eddies then grow in size both horizontally and vertically through an inverse energy cascade up to the energy-containing scale, which dominates the vertical distribution of EKE. As for smaller horizontal scales, the EKE distribution decays faster with depth. Guided by this interpretation, an SQG-based scale-aware parameterization of the EKE profile is proposed. Preliminary offline diagnosis of a high-resolution simulation shows the proposed scheme successfully reproducing the dependence of the vertical structure of EKE on the horizontal grid resolution.

Abstract

Two comprehensive Earth system models (ESMs), identical apart from their oceanic components, are used to estimate the uncertainty in projections of twenty-first-century sea level rise due to representational choices in ocean physical formulation. Most prominent among the formulation differences is that one (ESM2M) uses a traditional z-coordinate ocean model, while the other (ESM2G) uses an isopycnal-coordinate ocean. As evidence of model fidelity, differences in twentieth-century global-mean steric sea level rise are not statistically significant between either model and observed trends. However, differences between the two models’ twenty-first-century projections are systematic and both statistically and climatically significant. By 2100, ESM2M exhibits 18% higher global steric sea level rise than ESM2G for all four radiative forcing scenarios (28–49 mm higher), despite having similar changes between the models in the near-surface ocean for several scenarios. These differences arise primarily from the vertical extent over which heat is taken up and the total heat uptake by the models (9% more in ESM2M than ESM2G). The fact that the spun-up control state of ESM2M is warmer than ESM2G also contributes by giving thermal expansion coefficients that are about 7% larger in ESM2M than ESM2G. The differences between these models provide a direct estimate of the sensitivity of twenty-first-century sea level rise to ocean model formulation, and, given the span of these models across the observed volume of the ventilated thermocline, may also approximate the sensitivities expected from uncertainties in the characterization of interior ocean physical processes.

Abstract

Climate models simulate a wide range of climate changes at high northern latitudes in response to increased CO₂. They also have substantial disagreement on projected changes of the Atlantic meridional overturning circulation (AMOC). Here, two pairs of closely related climate models are used, with each containing members with large and small AMOC declines to explore the influence of AMOC decline on the high-latitude response to increased CO₂. The models with larger AMOC decline have less high-latitude warming and sea ice decline than their small AMOC decline counterpart. By examining differences in the perturbation heat budget of the 40°–90°N region, it is shown that AMOC decline diminishes the warming by weakening poleward ocean heat transport and increasing the ocean heat uptake. The cooling impact of this AMOC-forced surface heat flux perturbation difference is enhanced by shortwave feedback and diminished by longwave feedback and atmospheric heat transport differences. The magnitude of the AMOC decline within model pairs is positively related to the magnitudes of control climate AMOC and Labrador and Nordic Seas convection. Because the 40°–90°N region accounts for up to 40% of the simulated global ocean heat uptake over 100 yr, the process described here influences the global heat uptake efficiency.

Abstract

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 K _d —in the millennial-scale near-equilibrium climate response to a quadrupling of atmospheric CO₂. 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 K _d 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.

Abstract

The influence of alternative ocean and atmosphere subcomponents on climate model simulation of transient sensitivities is examined by comparing three GFDL climate models used for phase 5 of the Coupled Model Intercomparison Project (CMIP5). The base model ESM2M is closely related to GFDL’s CMIP3 climate model version 2.1 (CM2.1), and makes use of a depth coordinate ocean component. The second model, ESM2G, is identical to ESM2M but makes use of an isopycnal coordinate ocean model. The authors compare the impact of this “ocean swap” with an “atmosphere swap” that produces the GFDL Climate Model version 3 (CM3) by replacing the AM2 atmospheric component with AM3 while retaining a depth coordinate ocean model. The atmosphere swap is found to have much larger influence on sensitivities of global surface temperature and Northern Hemisphere sea ice cover. The atmosphere swap also introduces a multidecadal response time scale through its indirect influence on heat uptake. Despite significant differences in their interior ocean mean states, the ESM2M and ESM2G simulations of these metrics of climate change are very similar, except for an enhanced high-latitude salinity response accompanied by temporarily advancing sea ice in ESM2G. In the ESM2G historical simulation this behavior results in the establishment of a strong halocline in the subpolar North Atlantic during the early twentieth century and an associated cooling, which are counter to observations in that region. The Atlantic meridional overturning declines comparably in all three models.

Oceanic overflows are bottom-trapped density currents originating in semienclosed basins, such as the Nordic seas, or on continental shelves, such as the Antarctic shelf. Overflows are the source of most of the abyssal waters, and therefore play an important role in the large-scale ocean circulation, forming a component of the sinking branch of the thermohaline circulation. As they descend the continental slope, overflows mix vigorously with the surrounding oceanic waters, changing their density and transport significantly. These mixing processes occur on spatial scales well below the resolution of ocean climate models, with the result that deep waters and deep western boundary currents are simulated poorly. The Gravity Current Entrainment Climate Process Team was established by the U.S. Climate Variability and Prediction (CLIVAR) Program to accelerate the development and implementation of improved representations of overflows within large-scale climate models, bringing together climate model developers with those conducting observational, numerical, and laboratory process studies of overflows. Here, the organization of the Climate Process Team is described, and a few of the successes and lessons learned during this collaboration are highlighted, with some emphasis on the well-observed Mediterranean overflow. The Climate Process Team has developed several different overflow parameterizations, which are examined in a hierarchy of ocean models, from comparatively well-resolved regional models to the largest-scale global climate models.

The recent retreat and speedup of outlet glaciers, as well as enhanced surface melting around the ice sheet margin, have increased Greenland's contribution to sea level rise to 0.6 ± 0.1 mm yr⁻¹ and its discharge of freshwater into the North Atlantic. The widespread, near-synchronous glacier retreat, and its coincidence with a period of oceanic and atmospheric warming, suggests a common climate driver. Evidence points to the marine margins of these glaciers as the region from which changes propagated inland. Yet, the forcings and mechanisms behind these dynamic responses are poorly understood and are either missing or crudely parameterized in climate and ice sheet models. Resulting projected sea level rise contributions from Greenland by 2100 remain highly uncertain.

This paper summarizes the current state of knowledge and highlights key physical aspects of Greenland's coupled ice sheet–ocean–atmosphere system. Three research thrusts are identified to yield fundamental insights into ice sheet, ocean, sea ice, and atmosphere interactions, their role in Earth's climate system, and probable trajectories of future changes: 1) focused process studies addressing critical glacier, ocean, atmosphere, and coupled dynamics; 2) sustained observations at key sites; and 3) inclusion of relevant dynamics in Earth system models.

Understanding the dynamic response of Greenland's glaciers to climate forcing constitutes both a scientific and technological frontier, given the challenges of obtaining the appropriate measurements from the glaciers' marine termini and the complexity of the dynamics involved, including the coupling of the ocean, atmosphere, glacier, and sea ice systems. Interdisciplinary and international cooperation are crucial to making progress on this novel and complex problem.

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.