Search

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.

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 (r2 = 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.

Greenland Ice Sheet (GIS) might have lost a large amount of its volume during the last interglacial and may do so again in the future due to climate warming. In this study, we test whether the climate response to the glacial meltwater is sensitive to its discharging location. Two fully coupled atmosphere–ocean general circulation models, CM2G and CM2M, which have completely different ocean components are employed to do the test. In each experiment, a prescribed freshwater flux of 0.1 Sv is discharged from one of the four locations around Greenland—Petermann, 79 North, Jacobshavn and Helheim glaciers. The results from both models show that the AMOC weakens more when the freshwater is discharged from the northern GIS (Petermann and 79 North) than when it is discharged from the southern GIS (Jacobshavn and Helheim), by 15% (CM2G) and 31% (CM2M) averaged over model year 50–300 (CM2G) and 70–300 (CM2M), respectively. This is due to easier access of the freshwater from northern GIS to the deepwater formation site in the Nordic Seas. In the long term (> 300 year), however, the AMOC change is nearly the same for freshwater discharged from any location of the GIS. The East Greenland current accelerates with time and eventually becomes significantly faster when the freshwater is discharged from the north than from the south. Therefore, freshwater from the north is transported efficiently towards the south first and then circulates back to the Nordic Seas, making its impact to the deepwater formation there similar to the freshwater discharged from the south. The results indicate that the details of the location of meltwater discharge matter if the short-term ( 300 years) climate response is focused upon.

Model and observational studies have concluded that geothermal heating significantly alters the global overturning circulation and the properties of the widely distributed Antarctic Bottom Water. Here two distinct geothermal heat flux datasets are tested under different experimental designs in a fully coupled model that mimics the control run of a typical Coupled Model Intercomparison Project (CMIP) climate model. The regional analysis herein reveals that bottom temperature and transport changes, due to the inclusion of geothermal heating, are propagated throughout the water column, most prominently in the Southern Ocean, with the background density structure and major circulation pathways acting as drivers of these changes. While geothermal heating enhances Southern Ocean abyssal overturning circulation by 20%–50%, upwelling of warmer deep waters and cooling of upper ocean waters within the Antarctic Circumpolar Current (ACC) region decrease its transport by 3–5 Sv (1 Sv = 106 m3 s−1). The transient responses in regional bottom temperature increases exceed 0.1°C. The large-scale features that are shown to transport anomalies far from their geothermal source all exist in the Southern Ocean. Such features include steeply sloping isopycnals, weak abyssal stratification, voluminous southward flowing deep waters and exported bottom waters, the ACC, and the polar gyres. Recently the Southern Ocean has been identified as a prime region for deep ocean warming; geothermal heating should be included in climate models to ensure accurate representation of these abyssal temperature changes.

The Arctic Ocean simulated in 14 global ocean-sea ice models in the framework of the Coordinated Ocean-ice Reference Experiments, phase II (CORE-II) is analyzed in this study. The focus is on the Arctic liquid freshwater (FW) sources and freshwater content (FWC). The models agree on the interannual variability of liquid FW transport at the gateways where the ocean volume transport determines the FW transport variability. The variation of liquid FWC is induced by both the surface FW flux (associated with sea ice production) and lateral liquid FW transport, which are in phase when averaged on decadal time scales. The liquid FWC shows an increase starting from the mid-1990s, caused by the reduction of both sea ice formation and liquid FW export, with the former being more significant in most of the models. The mean state of the FW budget is less consistently simulated than the temporal variability. The model ensemble means of liquid FW transport through the Arctic gateways compare well with observations. On average, the models have too high mean FWC, weaker upward trends of FWC in the recent decade than the observation, and low consistency in the temporal variation of FWC spatial distribution, which needs to be further explored for the purpose of model development.

We characterise the representation of the Southern Ocean water mass structure and sea ice within a suite of 15 global ocean-ice models run with the Coordinated Ocean-ice Reference Experiment Phase II (CORE-II) protocol. The main focus is the representation of the present (1988-2007) mode and intermediate waters, thus framing an analysis of winter and summer mixed layer depths; temperature, salinity, and potential vorticity structure; and temporal variability of sea ice distributions. We also consider the interannual variability over the same 20 year period. Comparisons are made between models as well as to observation-based analyses where available. The CORE-II models exhibit several biases relative to Southern Ocean observations, including an underestimation of the model mean mixed layer depths of mode and intermediate water masses in March (associated with greater ocean surface heat gain), and an overestimation in September (associated with greater high latitude ocean heat loss and a more northward winter sea-ice extent). In addition, the models have cold and fresh/warm and salty water column biases centred near 50 degrees S. Over the 1933-2007 period, the CORE-II models consistently simulate spatially variable trends in sea-ice concentration, surface freshwater fluxes, mixed layer depths, and 200-700 in ocean heat content. In particular, sea-ice coverage around most of the Antarctic continental shelf is reduced, leading to a cooling and freshening of the near surface waters. The shoaling of the mixed layer is associated with increased surface buoyancy gain, except in the Pacific where sea ice is also influential. The models are in disagreement, despite the common CORE-II atmospheric state, in their spatial pattern of the 20-year trends in the mixed layer depth and sea-ice. (C) 2015 Elsevier Ltd. All rights reserved

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

The influence of changing ocean currents on climate change is evaluated by comparing an earth system model’s response to increased CO2 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.

Highlights: • We compare the simulated Arctic Ocean in 15 global ocean–sea ice models. • There is a large spread in temperature bias in the Arctic Ocean between the models. • Warm bias models have a strong temperature anomaly of inflow of Atlantic Water. • Dense outflows formed on Arctic shelves are not captured accurately in the models. In this paper we compare the simulated Arctic Ocean in 15 global ocean-sea ice models in the framework of the Coordinated Ocean-ice Reference Experiments, phase II (CORE-II). Most of these models are the ocean and sea-ice components of the coupled climate models used in the Coupled Model Intercomparison Project Phase 5 (CMIP5) experiments. We mainly focus on the hydrography of the Arctic interior, the state of Atlantic Water layer and heat and volume transports at the gateways of the Davis Strait, the Bering Strait, the Fram Strait and the Barents Sea Opening. We found that there is a large spread in temperature in the Arctic Ocean between the models, and generally large differences compared to the observed temperature at intermediate depths. Warm bias models have a strong temperature anomaly of inflow of the Atlantic Water entering the Arctic Ocean through the Fram Strait. Another process that is not represented accurately in the CORE-II models is the formation of cold and dense water, originating on the eastern shelves. In the cold bias models, excessive cold water forms in the Barents Sea and spreads into the Arctic Ocean through the St. Anna Through. There is a large spread in the simulated mean heat and volume transports through the Fram Strait and the Barents Sea Opening. The models agree more on the decadal variability, to a large degree dictated by the common atmospheric forcing. We conclude that the CORE-II model study helps us to understand the crucial biases in the Arctic Ocean. The current coarse resolution state-of-the-art ocean models need to be improved in accurate representation of the Atlantic Water inflow into the Arctic and density currents coming from the shelves.

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.

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.

A parameterization for the restratification by finite-amplitude, submesoscale, mixed layer eddies, formulated as an overturning streamfunction, has been recently proposed to approximate eddy fluxes of density and other tracers. Here, the technicalities of implementing the parameterization in the coarse-resolution ocean component of global climate models are made explicit, and the primary impacts on model solutions of implementing the parameterization are discussed. Three global ocean general circulation models including this parameterization are contrasted with control simulations lacking the parameterization. The MLE parameterization behaves as expected and fairly consistently in models differing in discretization, boundary layer mixing, resolution, and other parameterizations. The primary impact of the parameterization is a shoaling of the mixed layer, with the largest effect in polar winter regions. Secondary impacts include strengthening the Atlantic meridional overturning while reducing its variability, reducing CFC and tracer ventilation, modest changes to sea surface temperature and air–sea fluxes, and an apparent reduction of sea ice basal melting., National Science Foundation (U.S.) (Grant OCE-0612143), National Science Foundation (U.S.) (Grant OCE-0612059), National Science Foundation (U.S.) (Grant OCE-0825376), National Science Foundation (U.S.) (Grant DMS-0855010), National Science Foundation (U.S.) (Grant OCE-0934737)

The Arctic Ocean simulated in fourteen global ocean-sea ice models in the framework of the Coordinated Ocean-ice Reference Experiments, phase II (CORE II) is analyzed. The focus is on the Arctic sea ice extent, the solid freshwater (FW) sources and solid freshwater content (FWC). Available observations are used for model evaluation. The variability of sea ice extent and solid FW budget is more consistently reproduced than their mean state in the models. The descending trend of September sea ice extent is well simulated in terms of the model ensemble mean. Models overestimating sea ice thickness tend to underestimate the descending trend of September sea ice extent. The models underestimate the observed sea ice thinning trend by a factor of two. When averaged on decadal time scales, the variation of Arctic solid FWC is contributed by those of both sea ice production and sea ice transport, which are out of phase in time. The solid FWC decreased in the recent decades, caused mainly by the reduction in sea ice thickness. The models did not simulate the acceleration of sea ice thickness decline, leading to an underestimation of solid FWC trend after 2000. The common model behavior, including the tendency to underestimate the trend of sea ice thickness and March sea ice extent, remains to be improved.

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. (C) 2015 Elsevier Ltd. All rights reserved.

In many global ocean climate models, mesoscale eddies are parameterized as along isopycnal diffusion and eddyinduced advection (or equivalently skew-diffusion). The eddy-induced advection flattens isopycnals and acts as a sink of available potential energy, whereas the isopycnal diffusion mixes tracers along neutral directions. While much effort has gone into estimating diffusivities associated with this closure, less attention has been paid to the details of how this closure (which tries to flatten isopycnals) interacts with the mixed layer (in which vertical mixing tries to drive the isopycnals vertical). In order to maintain numerical stability, models often stipulate a maximum slope Smax which in combination with the thickness diffusivity Agm defines a maximum eddy-induced advective transport Agm*Smax. In this paper, we examine the impact of changing Smax within the GFDL global coupled climate model. We show that this parameter produces significant changes in wintertime mixed layer depth, with implications for wintertime temperatures in key regions, the distribution of precipitation, and the vertical structure of heat uptake. Smaller changes are seen in details of ventilation and currents, and even smaller changes as regards the overall hydrography. The results suggest that not only the value of the coefficient, but the details of the tapering scheme, need to be considered when comparing isopycnal mixing schemes in models.

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.

This paper summarizes the formulation of the ocean component to the Geophysical Fluid Dynamics Laboratory's (GFDL) climate model used for the 4th IPCC Assessment (AR4) of global climate change. In particular, it reviews the numerical schemes and physical parameterizations that make up an ocean climate model and how these schemes are pieced together for use in a state-of-the-art climate model. Features of the model described here include the following: (1) tripolar grid to resolve the Arctic Ocean without polar filtering, (2) partial bottom step representation of topography to better represent topographically influenced advective and wave processes, (3) more accurate equation of state, (4) three-dimensional flux limited tracer advection to reduce overshoots and undershoots, (5) incorporation of regional climatological variability in shortwave penetration, (6) neutral physics parameterization for representation of the pathways of tracer transport, (7) staggered time stepping for tracer conservation and numerical efficiency, (8) anisotropic horizontal viscosities for representation of equatorial currents, (9) parameterization of exchange with marginal seas, (10) incorporation of a free surface that accomodates a dynamic ice model and wave propagation, (11) transport of water across the ocean free surface to eliminate unphysical ``virtual tracer flux" methods, (12) parameterization of tidal mixing on continental shelves. We also present preliminary analyses of two particularly important sensitivities isolated during the development process, namely the details of how parameterized subgridscale eddies transport momentum and tracers.

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.

Highlights: • Global mean sea level simulated in interannual CORE simulations. • Regional sea level patterns simulated in interannual CORE simulations. • Theoretical foundation for analysis of global mean sea level and regional patterns. Abstract: We provide an assessment of sea level simulated in a suite of global ocean-sea ice models using the interannual CORE atmospheric state to determine surface ocean boundary buoyancy and momentum fluxes. These CORE-II simulations are compared amongst themselves as well as to observation-based estimates. We focus on the final 15 years of the simulations (1993–2007), as this is a period where the CORE-II atmospheric state is well sampled, and it allows us to compare sea level related fields to both satellite and in situ analyses. The ensemble mean of the CORE-II simulations broadly agree with various global and regional observation-based analyses during this period, though with the global mean thermosteric sea level rise biased low relative to observation-based analyses. The simulations reveal a positive trend in dynamic sea level in the west Pacific and negative trend in the east, with this trend arising from wind shifts and regional changes in upper 700 m ocean heat content. The models also exhibit a thermosteric sea level rise in the subpolar North Atlantic associated with a transition around 1995/1996 of the North Atlantic Oscillation to its negative phase, and the advection of warm subtropical waters into the subpolar gyre. Sea level trends are predominantly associated with steric trends, with thermosteric effects generally far larger than halosteric effects, except in the Arctic and North Atlantic. There is a general anti-correlation between thermosteric and halosteric effects for much of the World Ocean, associated with density compensated changes.

By convention, the ocean’s large-scale circulation is assumed to be a thermohaline overturning driven by the addition and extraction of buoyancy at the surface and vertical mixing in the interior. Previous work suggests that the overturning should die out as vertical mixing rates are reduced to zero. In this paper, a formal energy analysis is applied to a series of ocean general circulation models to evaluate changes in the large-scale circulation over a range of vertical mixing rates. Two different model configurations are used. One has an open zonal channel and an Antarctic Circumpolar Current (ACC). The other configuration does not. The authors find that a vigorous large-scale circulation persists at the limit of no mixing in the model with a wind-driven ACC. A wind-powered overturning circulation linked to the ACC can exist without vertical mixing and without much energy input from surface buoyancy forces.

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 Oceanice 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. (C) 2013 Elsevier Ltd. All rights reserved.

Brine rejection during the formation of Antarctic sea ice is known to enhance the salinity of dense shelf waters in the Weddell and Ross Seas. As these shelf waters flow off the shelves and descend to the bottom, they entrain ambient deep water to create new bottom water. It is not uncommon for ocean modelers to modify salinity boundary conditions around Antarctica in an attempt to include a “sea ice effect” in their models. However, the degree to which Antarctic salinities are enhanced is usually not quantified or defended. In this paper, studies of shelf hydrography and δ 18O are reviewed to assess the level of salinity enhancement appropriate for ocean general circulation models. The relevant quantities are 1) the salinity difference between the water masses modified on the shelves and the final offshelf flow and 2) the flux of salt (or freshwater) that gives rise to this salinity difference. Onshelf/offshelf salinity changes in the Weddell and Ross Seas appear to be fairly small, 0.15–0.20 sa...

[1] A simple model of the temperature-dependent biological decay of dissolved oil is embedded in an ocean-climate model and used to simulate underwater plumes of dissolved and suspended oil originating from a point source in the northern Gulf of Mexico, with an upper-bound supply rate estimated from the contemporary analysis of the Deepwater Horizon blowout. The behavior of plumes at different depths is found to be determined by the combination of sheared current strength and the vertical profile of decay rate. For all plume scenarios, toxic levels of dissolved oil remain confined to the northern Gulf of Mexico, and abate within weeks after the spill stops. An estimate of oxygen consumption due to microbial oxidation of hydrocarbons suggests that a deep plume of hydrocarbons could lead to localized regions of prolonged hypoxia near the source, but only when oxidation of methane is included.

International audience; 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.

[1] The impact of the penetration length scale of shortwave radiation into the surface ocean is investigated with a fully coupled ocean, atmosphere, land and ice model. Oceanic shortwave radiation penetration is assumed to depend on the chlorophyll concentration. As chlorophyll concentrations increase the distribution of shortwave heating becomes shallower. This change in heat distribution impacts mixed-layer depth. This study shows that removing all chlorophyll from the ocean results in a system that tends strongly towards an El Nino state—suggesting that chlorophyll is implicated in maintenance of the Pacific cold tongue. The regions most responsible for this response are located off-equator and correspond to the oligotrophic gyres. Results from a suite of surface chlorophyll perturbation experiments suggest a potential positive feedback between chlorophyll concentration and a non-local coupled response in the fully coupled ocean-atmosphere system.

[1] This paper considers the impact of the parameterization of subgridscale mixing on ocean heat transport in coarseresolution ocean models of the type used in coupled climate models. Increasing the vertical diffusion increases poleward heat transport in both hemispheres. Increasing lateral diffusion associated with transient eddies increases poleward heat transport in the southern hemisphere while decreasing it in the northern hemisphere. The results are interpreted in the context of a simple analytical model. INDEX TERMS: 4532 Oceanography: Physical: General circulation; 4568 Turbulence, diffusion, and mixing processes; 4203 Oceanography: General: Analytical modeling; 4255 Numerical modeling; 4279 Upwelling and convergences. Citation: Gnanadesikan, A., B. L. Samuels, and R. D. Slater, Sensitivity of water mass transformation and heat transport to subgridscale mixing in coarse-resolution ocean models, Geophys. Res. Lett., 30(18), 1967, doi:10.1029/2003GL018036, 2003.