Localization of deep water formation: role of atmospheric moisture transport and geometrical constraints on ocean circulation

A series of coupled atmosphere–ocean–ice aquaplanet experiments is described in which topological constraints on ocean circulation are introduced to study the role of ocean circulation on the mean climate of the coupled system. It is imagined that the earth is completely covered by an ocean of uniform depth except for the presence or absence of narrow barriers that extend from the bottom of the ocean to the sea surface. The following four configurations are described: Aqua (no land), Ridge (one barrier extends from pole to pole), Drake (one barrier extends from the North Pole to 35°S), and DDrake (two such barriers are set 90° apart and join at the North Pole, separating the ocean into a large basin and a small basin, connected to the south). On moving from Aqua to Ridge to Drake to DDrake, the energy transports in the equilibrium solutions become increasingly “realistic,” culminating in DDrake, which has an uncanny resemblance to the present climate. Remarkably, the zonal-average climates of Drake and DDrake are strikingly similar, exhibiting almost identical heat and freshwater transports, and meridional overturning circulations. However, Drake and DDrake differ dramatically in their regional climates. The small and large basins of DDrake exhibit distinctive Atlantic-like and Pacific-like characteristics, respectively: the small basin is warmer, saltier, and denser at the surface than the large basin, and is the main site of deep water formation with a deep overturning circulation and strong northward ocean heat transport. A sensitivity experiment with DDrake demonstrates that the salinity contrast between the two basins, and hence the localization of deep convection, results from a deficit of precipitation, rather than an excess of evaporation, over the small basin. It is argued that the width of the small basin relative to the zonal fetch of atmospheric precipitation is the key to understanding this salinity contrast. Finally, it is argued that many gross features of the present climate are consequences of two topological asymmetries that have profound effects on ocean circulation: a meridional asymmetry (circumpolar flow in the Southern Hemisphere; blocked flow in the Northern Hemisphere) and a zonal asymmetry (a small basin and a large basin).

Ocean heat transport and water vapor greenhouse in a warm equable climate: a new look at the low gradient paradox

The authors study the role of ocean heat transport (OHT) in the maintenance of a warm, equable, ice-free climate. An ensemble of idealized aquaplanet GCM calculations is used to assess the equilibrium sensitivity of global mean surface temperature and its equator-to-pole gradient (ΔT) to variations in OHT, prescribed through a simple analytical formula representing export out of the tropics and poleward convergence. Low-latitude OHT warms the mid- to high latitudes without cooling the tropics; increases by 1°C and ΔT decreases by 2.6°C for every 0.5-PW increase in OHT across 30° latitude. This warming is relatively insensitive to the detailed meridional structure of OHT. It occurs in spite of near-perfect atmospheric compensation of large imposed variations in OHT: the total poleward heat transport is nearly fixed. The warming results from a convective adjustment of the extratropical troposphere. Increased OHT drives a shift from large-scale to convective precipitation in the midlatitude storm tracks. Warming arises primarily from enhanced greenhouse trapping associated with convective moistening of the upper troposphere. Warming extends to the poles by atmospheric processes even in the absence of high-latitude OHT. A new conceptual model for equable climates is proposed, in which OHT plays a key role by driving enhanced deep convection in the midlatitude storm tracks. In this view, the climatic impact of OHT depends on its effects on the greenhouse properties of the atmosphere, rather than its ability to increase the total poleward energy transport.

Controlling spurious diapycnal mixing in eddy-resolving height-coordinate ocean models: insights from virtual deliberate tracer release experiments

A perceived limitation of z-coordinate models associated with spurious diapycnal mixing in eddying, frontal flow, can be readily addressed through appropriate attention to the tracer advection schemes employed. It is demonstrated that tracer advection schemes developed by Prather and collaborators for application in the stratosphere, greatly improve the fidelity of eddying flows, reducing levels of spurious diapycnal mixing to below those directly measured in field experiments, ∼1 × 10−5 m2 s−1. This approach yields a model in which geostrophic eddies are quasi-adiabatic in the ocean interior, so that the residual-mean overturning circulation aligns almost perfectly with density contours. A reentrant channel configuration of the MIT General Circulation Model, that approximates the Antarctic Circumpolar Current, is used to examine these issues. Virtual analogs of ocean deliberate tracer release field experiments reinforce our conclusion, producing passive tracer solutions that parallel field experiments remarkably well.

Predictable climate impacts of the decadal changes in the ocean in the 1990s

During the 1990s there was a major change in the state of the world's oceans. In particular, the North Atlantic underwent a rapid warming, with sea surface temperatures (SSTs) in the subpolar gyre region increasing by 1°C in just a few years. Associated with the changes in SST patterns were changes in the surface climate, in particular, a tendency for warm and dry conditions over areas of North America in all seasons, and warm springs and wet summers over areas of Europe. Here, the extent to which a climate prediction system initialized using observations of the ocean state is able to capture the observed changes in seasonal mean surface climate is investigated. Rather than examining predictions of the mid-1990s North Atlantic warming event itself, this study compares hindcasts started before and after the warming, relative to hindcasts that do not assimilate information. It is demonstrated that the hindcasts capture many aspects of the observed changes in seasonal mean surface climate, especially in North, South, and Central America and in Europe. Furthermore, the prediction system retains skill beyond the first year. Finally, it is shown that, in addition to memory of Atlantic SSTs, successfully predicting Pacific SSTs was likely important for the hindcasts to predict surface climate over North America.

Mid-Holocene monsoons: a multi-model analysis of the inter-hemispheric differences in the responses to orbital forcing and ocean feedbacks

The response of monsoon circulation in the northern and southern hemisphere to 6 ka orbital forcing has been examined in 17 atmospheric general circulation models and 11 coupled ocean–atmosphere general circulation models. The atmospheric response to increased summer insolation at 6 ka in the northern subtropics strengthens the northern-hemisphere summer monsoons and leads to increased monsoonal precipitation in western North America, northern Africa and China; ocean feedbacks amplify this response and lead to further increase in monsoon precipitation in these three regions. The atmospheric response to reduced summer insolation at 6 ka in the southern subtropics weakens the southern-hemisphere summer monsoons and leads to decreased monsoonal precipitation in northern South America, southern Africa and northern Australia; ocean feedbacks weaken this response so that the decrease in rainfall is smaller than might otherwise be expected. The role of the ocean in monsoonal circulation in other regions is more complex. There is no discernable impact of orbital forcing in the monsoon region of North America in the atmosphere-only simulations but a strong increase in precipitation in the ocean–atmosphere simulations. In contrast, there is a strong atmospheric response to orbital forcing over northern India but ocean feedback reduces the strength of the change in the monsoon although it still remains stronger than today. Although there are differences in magnitude and exact location of regional precipitation changes from model to model, the same basic mechanisms are involved in the oceanic modulation of the response to orbital forcing and this gives rise to a robust ensemble response for each of the monsoon systems. Comparison of simulated and reconstructed changes in regional climate suggest that the coupled ocean–atmosphere simulations produce more realistic changes in the northern-hemisphere monsoons than atmosphere-only simulations, though they underestimate the observed changes in precipitation in all regions. Evaluation of the southern-hemisphere monsoons is limited by lack of quantitative reconstructions, but suggest that model skill in simulating these monsoons is limited.

Sensitivity analysis of an ocean carbon cycle model in the North Atlantic: an investigation of parameters affecting the air-sea CO2 flux, primary production and export of detritus

The sensitivity of the biological parameters in a nutrient-phytoplankton-zooplankton-detritus (NPZD) model in the calculation of the air-sea CO2 flux, primary production and detrital export is analysed. We explore the effect on these outputs of variation in the values of the twenty parameters that control ocean ecosystem growth in a 1-D formulation of the UK Met Office HadOCC NPZD model used in GCMs. We use and compare the results from one-at-a-time and all-at-a-time perturbations performed at three sites in the EuroSITES European Ocean Observatory Network: the Central Irminger Sea (60° N 40° W), the Porcupine Abyssal Plain (49° N 16° W) and the European Station for Time series in the Ocean Canary Islands (29° N 15° W). Reasonable changes to the values of key parameters are shown to have a large effect on the calculation of the air-sea CO2 flux, primary production, and export of biological detritus to the deep ocean. Changes in the values of key parameters have a greater effect in more productive regions than in less productive areas. The most sensitive parameters are generally found to be those controlling well-established ocean ecosystem parameterisations widely used in many NPZD-type models. The air-sea CO2 flux is most influenced by variation in the parameters that control phytoplankton growth, detrital sinking and carbonate production by phytoplankton (the rain ratio). Primary production is most sensitive to the parameters that define the shape of the photosynthesis-irradiance curve. Export production is most sensitive to the parameters that control the rate of detrital sinking and the remineralisation of detritus.

Modeling ocean primary production: sensitivity to spectral resolution of attenuation and absorption of light

Modeling the vertical penetration of photosynthetically active radiation (PAR) through the ocean, and its utilization by phytoplankton, is fundamental to simulating marine primary production. The variation of attenuation and absorption of light with wavelength suggests that photosynthesis should be modeled at high spectral resolution, but this is computationally expensive. To model primary production in global 3d models, a balance between computer time and accuracy is necessary. We investigate the effects of varying the spectral resolution of the underwater light field and the photosynthetic efficiency of phytoplankton (α∗), on primary production using a 1d coupled ecosystem ocean turbulence model. The model is applied at three sites in the Atlantic Ocean (CIS (∼60°N), PAP (∼50°N) and ESTOC (∼30°N)) to include the effect of different meteorological forcing and parameter sets. We also investigate three different methods for modeling α∗ – as a fixed constant, varying with both wavelength and chlorophyll concentration [Bricaud, A., Morel, A., Babin, M., Allali, K., Claustre, H., 1998. Variations of light absorption by suspended particles with chlorophyll a concentration in oceanic (case 1) waters. Analysis and implications for bio-optical models. J. Geophys. Res. 103, 31033–31044], and using a non-spectral parameterization [Anderson, T.R., 1993. A spectrally averaged model of light penetration and photosynthesis. Limnol. Oceanogr. 38, 1403–1419]. After selecting the appropriate ecosystem parameters for each of the three sites we vary the spectral resolution of light and α∗ from 1 to 61 wavebands and study the results in conjunction with the three different α∗ estimation methods. The results show modeled estimates of ocean primary productivity are highly sensitive to the degree of spectral resolution and α∗. For accurate simulations of primary production and chlorophyll distribution we recommend a spectral resolution of at least six wavebands if α∗ is a function of wavelength and chlorophyll, and three wavebands if α∗ is a fixed value.

On the role of the ocean in projected atmospheric stability changes in the Atlantic polar low region

The occurrence of destructive mesoscale ‘polar low’ cyclones in the subpolar North Atlantic is projected to decline under anthropogenic change, due to an increase in atmospheric static stability. This letter reports on the role of changes in ocean circulation in shaping the atmospheric stability. In particular, the Atlantic Meridional Overturning Circulation (AMOC) is projected to weaken in response to anthropogenic forcing, leading to a local minimum in warming in this region. The reduced warming is restricted to the lower troposphere, hence contributing to the increase in static stability. Linear correlation analysis of the CMIP3 climate model ensemble suggests that around half of the model uncertainty in the projected stability response arises from the varied response of the AMOC between models.

Impact of strong deep ocean stratification on the glacial carbon cycle

During the Last Glacial Maximum, the climate was substantially colder and the carbon cycle was clearly different from the late Holocene. According to proxy data deep oceanic δ13C was very low, and the atmospheric CO2 concentration also reduced. Several mechanisms have been proposed to explain these changes, but none can fully explain the data, especially the very low deep ocean δ13C values. Oceanic core data show that the deep ocean was very cold and salty, which would lead to enhanced deep ocean stratification. We show that such an enhanced stratification in the coupled climate model CLIMBER-2 helps get very low deep oceanic δ13C values. Indeed the simulated δ13C reaches values as low as −0.8‰ in line with proxy data evidences. Moreover it increases the oceanic carbon reservoir leading to a small, yet robust, atmospheric CO2 drop of approximately 10 ppm.

Quantifying the roles of ocean circulation and biogeochemistry in governing ocean carbon-13 and atmospheric carbon dioxide at the last glacial maximum

We use a state-of-the-art ocean general circulation and biogeochemistry model to examine the impact of changes in ocean circulation and biogeochemistry in governing the change in ocean carbon-13 and atmospheric CO2 at the last glacial maximum (LGM). We examine 5 different realisations of the ocean's overturning circulation produced by a fully coupled atmosphere-ocean model under LGM forcing and suggested changes in the atmospheric deposition of iron and phytoplankton physiology at the LGM. Measured changes in carbon-13 and carbon-14, as well as a qualitative reconstruction of the change in ocean carbon export are used to evaluate the results. Overall, we find that while a reduction in ocean ventilation at the LGM is necessary to reproduce carbon-13 and carbon-14 observations, this circulation results in a low net sink for atmospheric CO2. In contrast, while biogeochemical processes contribute little to carbon isotopes, we propose that most of the change in atmospheric CO2 was due to such factors. However, the lesser role for circulation means that when all plausible factors are accounted for, most of the necessary CO2 change remains to be explained. This presents a serious challenge to our understanding of the mechanisms behind changes in the global carbon cycle during the geologic past.

The effect of windstress change on future sea level change in the Southern Ocean

AOGCMs of the two latest phases (CMIP3 and CMIP5) of the Coupled Model Intercomparison Project, like earlier AOGCMs, predict large regional variations in future sea level change. The model-mean pattern of change in CMIP3 and CMIP5 is very similar, and its most prominent feature is a zonal dipole in the Southern Ocean: sea level rise is larger than the global mean north of 50°S and smaller than the global mean south of 50°S in most models. The individual models show widely varying patterns, although the inter-model spread in local sea level change is smaller in CMIP5 than in CMIP3. Here we investigate whether changes in windstress can explain the different patterns of projected sea level change, especially the Southern Ocean feature, using two AOGCMs forced by the changes in windstress from the CMIP3 and CMIP5 AOGCMs. We show that the strengthening and poleward shift of westerly windstress accounts for the most of the large spread among models in magnitude of this feature. In the Indian, North Pacific and Arctic Oceans, the windstress change is influential, but does not completely account for the projected sea level change.

The residual circulation of the Southern Ocean: which spatio-temporal scales are needed?

The Southern Ocean circulation consists of a complicated mixture of processes and phenomena that arise at different time and spatial scales which need to be parametrized in the state-of-the-art climate models. The temporal and spatial scales that give rise to the present-day residual mean circulation are here investigated by calculating the Meridional Overturning Circulation (MOC) in density coordinates from an eddy-permitting global model. The region sensitive to the temporal decomposition is located between 38°S and 63°S, associated with the eddy-induced transport. The ‘‘Bolus’’ component of the residual circulation corresponds to the eddy-induced transport. It is dominated by timescales between 1 month and 1 year. The temporal behavior of the transient eddies is examined in splitting the ‘‘Bolus’’ component into a ‘‘Seasonal’’, an ‘‘Eddy’’ and an ‘‘Inter-monthly’’ component, respectively representing the correlation between density and velocity fluctuations due to the average seasonal cycle, due to mesoscale eddies and due to large-scale motion on timescales longer than one month that is not due to the seasonal cycle. The ‘‘Seasonal’’ bolus cell is important at all latitudes near the surface. The ‘‘Eddy’’ bolus cell is dominant in the thermocline between 50°S and 35°S and over the whole ocean depth at the latitude of the Drake Passage. The ‘‘Inter-monthly’’ bolus cell is important in all density classes and is maximal in the Brazil–Malvinas Confluence and the Agulhas Return Current. The spatial decomposition indicates that a large part of the Eulerian mean circulation is recovered for spatial scales larger than 11.25°, implying that small-scale meanders in the Antarctic Circumpolar Current (ACC), near the Subantarctic and Polar Fronts, and near the Subtropical Front are important in the compensation of the Eulerian mean flow.

The vertical structure of ocean heat transport

One of the most important contributions the ocean makes to Earth's climate is through its poleward heat transport: about 1.5 PW or more than 30% of that accomplished by the ocean-atmosphere system (Trenberth and Caron, 2001). Recently, concern has arisen over whether global warming could affect this heat transport (Watson et al., 2001), for example, reducing high latitude convection and triggering a collapse of the deep overturning circulation (Rahmstorf, 1995). While the consequences of abrupt changes in oceanic circulation should be of concern, we argue that the attention devoted to deep circulations is disproportionate to their role in heat transport. For this purpose, we introduce a heat function which identifies the contribution to the heat transport by different components of the oceanic circulation. A new view of the ocean emerges in which a shallow surface intensified circulation dominates the poleward heat transport.