A GCM study of Antarctic glaciation

An atmospheric general circulation model, the NCAR CCM, has been used to investigate the possible effects of two specific tectonic mechanisms on Antarctic glaciation. These two mechanisms are: (1) closing the Drake Passage (connecting South America with Antarctica), which is assumed to effectively represent an increased meridional heat transport by the ocean; and (2) changing the elevation of Antarctica. Perpetual season (summer and winter) and seasonal cycle simulations with warmer sea-surface temperatures and no sea ice prescribed for mid- to high-latitude southern oceans have been made with both present-day (high) Antarctic elevations and with low Antarctic elevations (all points 200 m). The results suggest a relatively minor role for oceanic heat transport in the formation/elimination of Antarctic glaciation. That is, under the warmer conditions inferred to have prevailed prior to the opening of the Drake Passage, conditions would still have been favorable for the maintenance of an Antarctic ice-sheet. If anything, a moderate ocean warming would promote glaciation, by increasing snowfall. Lowering the elevation of Antarctica has a larger effect on the model simulations, reducing the likelihood of glacial conditions. In the absence of snowcover, summer temperatures over Antarctica can warm considerably, leading to a monsoon-like circulation. However, it may be difficult to achieve such snow-free conditions, even with greatly increased atmospheric carbon dioxide. A tundra-like climate is the closest the model has come to representing a non-glacial climate, even when both seasurface temperatures and elevations are maximally varied.     

A coupled climate model simulation of the Last Glacial Maximum,Part 2: approach to equilibrium

The climate of the last glacial maximum (LGM) is simulated with a coupled climate model. The simulated climate undergoes a rapid adjustment during the first several decades after imposition of LGM boundary conditions, as described in Part 1, and then evolves toward equilibrium over 900 model years. The climate simulated by the coupled model at this period is compared with observationally-based LGM reconstructions and with LGM results obtained with an atmosphere-mixed layer (slab) ocean version of the model in order to investigate the role of ocean dynamics in the LGM climate. Global mean surface air temperature and sea surface temperature (SST) decrease by about 10 °C and 5.6 °C in the coupled model which includes ocean dynamics, compared to decreases of 6.3 and 3.8 °C in slab ocean case. The coupled model simulates a cooling of about 6.5 °C over the tropics, which is larger than that of the CLIMAP reconstruction (1.7 °C) and larger than that of the slab ocean simulation (3.3 °C), but which is in reasonable agreement with some recent proxy estimates. The ocean dynamics of the coupled model captures features found in the CLIMAP reconstructions such as a relative maximum of ocean cooling over the tropical Pacific associated with a mean La Niña-like response and lead to a more realistic SST pattern than in the slab model case. The reduction in global mean precipitation simulated in the coupled model is larger (15%) than that simulated with the slab ocean model (~10%) in conjunction with the enhanced cooling. Some regions, such as the USA and the Mediterranean region, experience increased precipitation in accord with proxy paleoclimate evidence. The overall much drier climate over the ocean leads to higher sea surface salinity (SSS) in most ocean basins except for the North Atlantic where SSS is considerably lower due to an increase in the supply of fresh water from the Mississippi and Amazon rivers and presumably a decrease in salt transport by the weakened North Atlantic overturning circulation. The North Atlantic overturning stream function weakens to less than half of the control run value. The overturning is limited to a shallower depth (less than 1000 m) and its outflow is confined to the Northern Hemisphere. In the Southern Ocean, convection is much stronger than in the control run leading to a stronger overturning stream function associated with enhanced Antarctic Bottom Water formation. As a result, Southern Ocean water masses fill the entire deep ocean. The Antarctic Circumpolar Current (ACC) transport through the Drake Passage increases by about 25%. The ACC transport, despite weaker zonal winds, is enhanced due to changes in bottom pressure torque. The weakening of the overturning circulation in the North Atlantic and the accompanying 30% decrease in the poleward ocean heat transport contrasts with the strengthening of the overturning circulation in the Southern Ocean and a 40% increase in heat transport. As a result, sea ice coverage and thickness are affected in opposite senses in the two hemispheres. The LGM climate simulated by the coupled model is in reasonable agreement with paleoclimate proxy evidence. The dynamical response of the ocean in the coupled model plays an important role in determining the simulated, and undoubtedly, the actual, LGM climate.     

Impacts of ocean gateway and basin width on Tertiary tropical climate variability in a prototype model

Using a simple tropical climate model, we investigated possible impacts of changes in oceanic seaways (Panama and Tethys) and ocean basin sizes (great Pacific and narrow Atlantic) on tropical climate variability during Tertiary. Our model showed that the opening of seaways had little influence on climate variability in the tropical Pacific because the climate variability in the Pacific Ocean’s large basins were internally generated, regardless of the variation in the tropical Atlantic Ocean. Conversely, the climate variability in the tropical Atlantic Ocean was highly dependent on the tropical Pacific Ocean; thus, an opening seaway, particularly the Panama seaway, was crucial in generating the interannual variability in the tropical Atlantic Ocean. We also found that in the Pacific Ocean, basin size strongly modified the period and amplitude of the interannual variability of both the Pacific and Atlantic Oceans due to ocean wave dynamics.     

一个压力坐标海洋环流模式(PCOM1.0)的基本模拟性能评估

参照Griffies et al.(2009)提出的海洋—海冰耦合模式参考试验(Coordinated Ocean-ice Reference Experiments,COREs),设计了一个800年积分的数值试验,对一个质量严格守恒的压力坐标海洋环流模式(Pressure Coordinate Ocean Model,PCOM1.0)的基本模拟性能进行了评估,并与观测资料和再分析资料进行了对比。结果表明,PCOM1.0模拟的温盐场和基本流场与COREs模式的模拟水平基本接近。其中,模拟的大西洋经向翻转流在45°N附近达到18 Sv(1 Sv=106 m3 s-1),与观测估计值接近;对海表面温度的模拟误差主要集中在北太平洋黑潮区和北大西洋湾流区等中高纬度急流区;模拟的热带太平洋温跃层过于深厚;模拟的经德雷克海峡的体积输送达130 Sv,比大部分COREs模式及再分析资料都更接近于观测估计值。     

The Role of Ocean Dynamics in the Cross-equatorial Energy Transport under a Thermal Forcing in the Southern Ocean

Under external heating forcing in the Southern Ocean, climate models project anomalous northward atmosphere heat transport (AHT) across the equator, accompanied by a southward shift of the intertropical convergence zone (ITCZ). Comparison between a fully coupled and a slab ocean model shows that the inclusion of active ocean dynamics tends to partition the cross-equatorial energy transport and significantly reduce the ITCZ shift response by a factor of 10, a finding which supports previous studies. To understand how ocean dynamics damps the ITCZ's response to an imposed thermal heating in the Southern Ocean, we examine the ocean heat transport (OHT) and ocean circulation responses in a set of fully coupled experiments. Results show that both the Indo-Pacific and the Atlantic contribute to transport energy across the equator mainly through its Eulerian-mean component. However, different from previous studies that linked the changes in OHT to the changes in the wind-driven subtropical cells or the Atlantic meridional overturning circulation (AMOC), our results show that the cross-equatorial OHT anomaly is due to a broad clockwise overturning circulation anomaly below the subtropical cells (approximately bounded by the 5℃ to 20℃ isotherms and 50°S to 10°N). Further elimination of the wind-driven component, conducted by prescribing the climatological wind stress in the Southern Ocean heat perturbation experiments, leads to little change in OHT, suggesting that the OHT response is predominantly thermohaline-driven by air-sea thermal interactions.     

Meridional overturning circulation: stability and ocean feedbacks in a box model

A box model of the inter-hemispheric Atlantic meridional overturning circulation is developed, including a variable pycnocline depth for the tropical and subtropical regions. The circulation is forced by winds over a periodic channel in the south and by freshwater forcing at the surface. The model is aimed at investigating the ocean feedbacks related to perturbations in freshwater forcing from the atmosphere, and to changes in freshwater transport in the ocean. These feedbacks are closely connected with the stability properties of the meridional overturning circulation, in particular in response to freshwater perturbations. A separate box is used for representing the region north of the Antarctic circumpolar current in the Atlantic sector. The density difference between this region and the north of the basin is then used for scaling the downwelling in the north. These choices are essential for reproducing the sensitivity of the meridional overturning circulation observed in general circulation models, and therefore suggest that the southernmost part of the Atlantic Ocean north of the Drake Passage is of fundamental importance for the stability of the meridional overturning circulation. With this configuration, the magnitude of the freshwater transport by the southern subtropical gyre strongly affects the response of the meridional overturning circulation to external forcing. The role of the freshwater transport by the overturning circulation (M _ov ) as a stability indicator is discussed. It is investigated under which conditions its sign at the latitude of the southern tip of Africa can provide information on the existence of a second, permanently shut down, state of the overturning circulation in the box model. M _ov will be an adequate indicator of the existence of multiple equilibria only if salt-advection feedback dominates over other processes in determining the response of the circulation to freshwater anomalies. M _ov is a perfect indicator if feedbacks other than salt-advection are negligible.     

The role of meridional density differences for a wind-driven overturning circulation

Experiments with the coupled climate model CLIMBER-3α, which contains an oceanic general circulation model, show deep upwelling in the Southern Ocean to be proportional to the surface wind stress in the latitudinal band of Drake Passage. At the same time, the distribution of the Southern Ocean upwelling onto the oceanic basins is controlled by buoyancy distribution; the inflow into each basin being proportional to the respective meridional density difference. We observe approximately the same constant of proportionality for all basins, and demonstrate that it can be directly related to the flow geometry. For increased wind stress in the Southern Ocean, the overturning increases both in the Atlantic and the Indo-Pacific basin. For strongly reduced wind stress, the circulation enters a regime where Atlantic overturning is maintained through Pacific upwelling, in order to satisfy the transports set by the density differences. Previous results on surface buoyancy and wind stress forcing, obtained with different models, are reproduced within one model in order to distill a consistent picture. We propose that both Southern Ocean upwelling and meridional density differences set up a system of conditions that determine the global meridional overturning circulation.     

Atmospheric water vapor flux,bifurcation of the thermohaline circulation,and climate change

Latitudinal heat transport in the ocean and atmosphere represents a fundamental process of the Earth's climate system. The ocean component of heat transport is effected by the thermohaline circulation. Changes in this circulation, and hence latitudinal heat transport, would have a significant effect on global climate. Paleoclimate evidence from the Greenland ice cores and deep sea sediment cores suggests that during much of glacial time the climate system oscillated between two different states. Bimodal equilibrium states of the thermohaline circulation have been demonstrated in climate models. We address the question of the role of the atmospheric hydrological cycle on the global thermohaline circulation and the feedback to the climate system through changes in the ocean's latitudinal heat transport, with a simple coupled ocean-atmosphere energy-salt balance model. Two components of the atmospheric hydrological cycle, i.e., latitudinal water vapor transport and the net flux of water vapor from the Atlantic to the Pacific Ocean appear to play separate roles. If the inter-basin transport is sufficiently large, small changes in water vapor transport over the North Atlantic can effect bifurcation or a rapid transition between two different equilibria in the global thermohaline circulation; maximum difference between the modes occurs in the North Atlantic. If the inter-basin transport is from the Pacific to the Atlantic and sufficiently large, latitudinal vapor transport in the North Pacific controls the bifurcations, with maximum changes occurring in the North Pacific. For intermediate values of inter-basin transport, no rapid transitions occur in either basin. In the regime with vapor flux from the Atlantic to the Pacific, the on mode has strong production of deep water in the North Atlantic and a large flux of heat to the atmosphere from the high latitude North Atlantic. The off mode has strong deep water production in the Southern Ocean and weak production in the North Pacific. Heat transport into the high latitude North Atlantic by the ocean is reduced to about 20% of the on mode value. For estimated values of water vapor transport for the present climate the model asserts that while water vapor transport from the Atlantic to the Pacific Ocean is sufficiently large to make the North Atlantic the dominant region for deep water production, latitudinal water vapor transport is sufficiently low that the thermohaline circulation appears stable, i.e., far from a bifurcation point. This conclusion is supported to some extent by the fact that the high latitude temperature of the atmosphere as recorded in the Greenland ice cores has changed little over the last 9000 years.     

厄尔尼诺事件和拉尼娜事件的成因与预测

通过对南极气温资料、环南极海冰资料、臭氧变化资料、太平洋海温资料、地球自转速度变化资料、厄尔尼诺和拉尼娜资料的综合验证,发现了构造运动与厄尔尼诺因果关系。大气、海洋与岩石圈的角动量交换在南半球和北半球有不同的形式,这是由陆海分布的差异决定的。南极上空臭氧变化和环南极海冰变化是赤道海温和全球气候准两年振荡的原因。其中,德雷克海峡的海冰变化起主要作用。这个结论给出了作者提出的“海洋锅炉效应”、“臭氧洞漏能效应”、“德雷克海冰气候开关效应”和“大洋地壳跷跷板运动”的相互关系,证明构造运动对厄尔尼诺的重要影响。强潮汐准4a周期的发现,表明南极海冰变化、东太平洋海温变化、地球自转变化和厄尔尼诺都具有4a准周期变化的原因。海温和海冰开关的准2a周期和日食-厄尔尼诺系数理论有较好的预测效果。     

Impact of a realistic river routing in coupled ocean–atmosphere simulations of the Last Glacial Maximum climate

The presence of large ice sheets over North America and North Europe at the Last Glacial Maximum (LGM) strongly impacted Northern hemisphere river pathways. Despite the fact that such changes may significantly alter the freshwater input to the ocean, modified surface hydrology has never been accounted for in coupled ocean–atmosphere general circulation model simulations of the LGM climate. To reconstruct the LGM river routing, we use the ICE-5G LGM topography. Because of the uncertainties in the extent of the Fennoscandian ice sheet in the Eastern part of the Kara Sea, we consider two more realistic river routing scenarios. The first scenario is characterised by the presence of an ice dammed lake south of the Fennoscandian ice sheet, and corresponds to the ICE-5G topography. This lake is fed by the Ob and Yenisei rivers. In the second scenario, both these rivers flow directly into the Arctic Ocean, which is more consistent with the latest QUEEN ice sheet margin reconstructions. We study the impact of these changes on the LGM climate as simulated by the IPSL_CM4 model and focus on the overturning thermohaline circulation. A comparison with a classical LGM simulation performed using the same model and modern river basins as designed in the PMIP2 exercise leads to the following conclusions: (1) The discharge into the North Atlantic Ocean is increased by 2,000 m³/s between 38° and 54°N in both simulations that contain LGM river routing, compared to the classical LGM experiment. (2) The ice dammed lake is shown to have a weak impact, relative to the classical simulation, both in terms of climate and ocean circulation. (3) In contrast, the North Atlantic deep convection and meridional overturning are weaker than during the classical LGM run if the Ob and Yenisei rivers flow directly into the Arctic Ocean. The total discharge into the Arctic Ocean is increased by 31,000 m³/s, relative to the classical LGM simulation. Consequentially, northward ocean heat transport is weaker, and sea ice more extensive, in better agreement with existing proxy data.     

Freshwater transports in HadCM3

The hydrological cycle can influence climate through a great variety of processes. A good representation of the hydrological cycle in climate models is therefore crucial. Attempts to analyse the global hydrological cycle are hampered by a deficiency of suitable observations, particularly over the oceans. Fully coupled general circulation models are potentially powerful tools in interpreting the limited observational data in the context of large-scale freshwater exchanges. We have looked at large-scale aspects of the global freshwater budget in a simulation, of over 1000 years, by the Hadley Centre coupled climate model (HadCM3). Many aspects of the global hydrological cycle are well represented, but the model hydrological cycle appears to be too strong, with overly large precipitation and evaporation components in comparison with the observational datasets we have used. We show that the ocean basin-scale meridional transports of freshwater come into near balance with the surface freshwater fluxes on a time scale of about 400 years, with the major change being a relative increase of freshwater transport from the Southern Ocean into the Atlantic Ocean. Comparison with observations, supported by sensitivity tests, suggests that the major cause of a drift to more saline condition in the model Atlantic is an overestimate of evaporation, although other freshwater budget components may also play a role. The increase in ocean freshwater transport into the Atlantic during the simulation, primarily coming from the overturning circulation component, which changes from divergent to convergent, acts to balance this freshwater budget deficit. The stability of the thermohaline circulation in HadCM3 may be affected by these freshwater transport changes and this question is examined in the context of an existing conceptual model.     

Energy transports by ocean and atmosphere based on an entropy extremum principle. Part II: Two-dimensional transports

Summary The maximum entropy production (MEP) principle used in Part J has been extended to separate the two-dimensional required energy transports determined from Nimbus 7 satellite net radiation measurements into atmospheric and oceanic components. In terms of the meridional component of the ocean transport vectors, results show northward ocean transports throughout the entire Atlantic ocean from southern hemisphere high latitudes to northern hemisphere polar regions, southward transports throughout the entire Indian Ocean, and poleward transports separated at approximately 10°S in the Pacific Ocean. The ocean transport patterns are consistent with well-known features concerning heat transport within the three ocean basins. However, uncertainty remains in the magnitudes of the transports. Because of the large remaining discrepancies between published estimates based on direct measurements and indirect estimates derived from energy budget methods, assessing the accuracy of the magnitudes is difficult, although there is evidence that the limited model resolution leads to synergistic biases in the North Atlantic and North Pacific. In terms of the crossmeridional energy transport component, results suggest that most of the net energy transfer in the tropics takes place within the ocean. In the southern hemisphere high latitudes, the Pacific and Indian Oceans export heat cross-meridionally to the Atlantic Ocean through the passages below Cape Horn and the Cape of Good Hope, although the magnitudes of these inter-ocean heat exchanges are small. Another important aspect of the southern hemisphere results is that poleward transports are dominated by the atmospheric component with strong zonal asymmetry. By contrast, in the northern hemisphere, atmospheric transports over the ocean are generally weaker than the corresponding southern hemisphere terms, indicating that the northern hemisphere oceans are relatively more effective in transferring heat poleward. Finally, poleward atmospheric transports over the continental areas exceed those over the ocean at equivalent latitudes as a result of the generally greater energy deficits over the land areas.With 7 Figures     

Effects of variability of sea ice transport through the Fram Strait on the intensity of the Atlantic deep circulation

Results from an ice-ocean coupled model are used to investigate the impact of long-term variability in sea ice transport at the Fram Strait on the intensity of the Atlantic deep circulation. An increase (or decrease) in sea ice transport through the Fram Strait leads to a stronger (or weaker) deep circulation in the Atlantic. Change in the sea ice transport is accompanied by a salinity anomaly in the surface layer of the Arctic Ocean. Such an anomaly could inversely affect the Atlantic circulation once it reaches deep water formation regions. If the Canadian Archipelago is closed, the anomaly is subsequently transported through the Fram Strait, and counters the initial changes in the Atlantic deep circulation. On the other hand, if the Canadian Archipelago is open, some of the anomaly is transported to the Canadian Archipelago, and the initial change in the Atlantic deep circulation persists. In the Arctic Ocean basin, the time scale and path of the salinity anomalys propagation depends on the large-scale flow at the surface of the Arctic Ocean. Our results suggest that the salinity anomaly transport and its propagation pathway out of the Arctic Ocean are important determinants of the role of sea ice transport variability through the Fram Strait in controlling the intensity of the Atlantic deep circulation.     

Closure of the Panama Seaway during the Pliocene: implications for climate and Northern Hemisphere glaciation

The “Panama Hypothesis” states that the gradual closure of the Panama Seaway, between 13 million years ago (13 Ma) and 2.6 Ma, led to decreased mixing of Atlantic and Pacific water Masses, the formation of North Atlantic Deep water and strengthening of the Atlantic thermohaline circulation, increased temperatures and evaporation in the North Atlantic, increased precipitation in Northern Hemisphere (NH) high latitudes, culminating in the intensification of Northern Hemisphere Glaciation (NHG) during the Pliocene, 3.2–2.7 Ma. Here we test this hypothesis using a fully coupled, fully dynamic ocean-atmosphere general circulation model (GCM) with boundary conditions specific to the Pliocene, and a high resolution dynamic ice sheet model. We carry out two GCM simulations with “closed” and “open” Panama Seaways, and use the simulated climatologies to force the ice sheet model. We find that the models support the “Panama Hypothesis” in as much as the closure of the seaway results in a more intense Atlantic thermohaline circulation, enhanced precipitation over Greenland and North America, and ultimately larger ice sheets. However, the volume difference between the ice sheets in the “closed” and “open” configurations is small, equivalent to about 5 cm of sea level. We conclude that although the closure of the Panama Seaway may have slightly enhanced or advanced the onset of NHG, it was not a major forcing mechanism. Future work must fully couple the ice sheet model and GCM, and investigate the role of orbital and CO₂ effects in controlling NHG.     

A multi-model analysis of the role of the ocean on the African and Indian monsoon during the mid-Holocene

We investigate the role of the ocean feedback on the climate in response to insolation forcing during the mid-Holocene (6,000 year BP) using results from seven coupled ocean–atmosphere general circulation models. We examine how the dipole in late summer sea-surface temperature (SST) anomalies in the tropical Atlantic increases the length of the African monsoon, how this dipole structure is created and maintained, and how the late summer SST warming in the northwest Indian Ocean affects the monsoon retreat in this sector. Similar mechanisms are found in all of the models, including a strong wind evaporation feedback and changes in the mixed layer depth that enhance the insolation forcing, as well as increased Ekman transport in the Atlantic that sharpens the Atlantic dipole pattern. We also consider changes in interannual variability over West Africa and the Indian Ocean. The teleconnection between variations in SST and Sahelian precipitation favor a larger impact of the Atlantic dipole mode in this region. In the Indian Ocean, the strengthening of the Indian dipole structure in autumn has a damping effect on the Indian dipole mode at the interannual time scale.     

Kinetic energy analysis of the response of the Atlantic meridional overturning circulation to CO2-forced climate change

Atmosphere?Cocean general circulation models (AOGCMs) predict a weakening of the Atlantic meridional overturning circulation (AMOC) in response to anthropogenic forcing of climate, but there is a large model uncertainty in the magnitude of the predicted change. The weakening of the AMOC is generally understood to be the result of increased buoyancy input to the north Atlantic in a warmer climate, leading to reduced convection and deep water formation. Consistent with this idea, model analyses have shown empirical relationships between the AMOC and the meridional density gradient, but this link is not direct because the large-scale ocean circulation is essentially geostrophic, making currents and pressure gradients orthogonal. Analysis of the budget of kinetic energy (KE) instead of momentum has the advantage of excluding the dominant geostrophic balance. Diagnosis of the KE balance of the HadCM3 AOGCM and its low-resolution version FAMOUS shows that KE is supplied to the ocean by the wind and dissipated by viscous forces in the global mean of the steady-state control climate, and the circulation does work against the pressure-gradient force, mainly in the Southern Ocean. In the Atlantic Ocean, however, the pressure-gradient force does work on the circulation, especially in the high-latitude regions of deep water formation. During CO₂-forced climate change, we demonstrate a very good temporal correlation between the AMOC strength and the rate of KE generation by the pressure-gradient force in 50?C70°N of the Atlantic Ocean in each of nine contemporary AOGCMs, supporting a buoyancy-driven interpretation of AMOC changes. To account for this, we describe a conceptual model, which offers an explanation of why AOGCMs with stronger overturning in the control climate tend to have a larger weakening under CO₂ increase.     

The ocean’s role in setting the mean position of the Inter-Tropical Convergence Zone

Through study of observations and coupled climate simulations, it is argued that the mean position of the Inter-Tropical Convergence Zone (ITCZ) north of the equator is a consequence of a northwards heat transport across the equator by ocean circulation. Observations suggest that the hemispheric net radiative forcing of climate at the top of the atmosphere is almost perfectly symmetric about the equator, and so the total (atmosphere plus ocean) heat transport across the equator is small (order 0.2 PW northwards). Due to the Atlantic ocean’s meridional overturning circulation, however, the ocean carries significantly more heat northwards across the equator (order 0.4 PW) than does the coupled system. There are two primary consequences. First, atmospheric heat transport is southwards across the equator to compensate (0.2 PW southwards), resulting in the ITCZ being displaced north of the equator. Second, the atmosphere, and indeed the ocean, is slightly warmer (by perhaps 2 °C) in the northern hemisphere than in the southern hemisphere. This leads to the northern hemisphere emitting slightly more outgoing longwave radiation than the southern hemisphere by virtue of its relative warmth, supporting the small northward heat transport by the coupled system across the equator. To conclude, the coupled nature of the problem is illustrated through study of atmosphere–ocean–ice simulations in the idealized setting of an aquaplanet, resolving the key processes at work.     

The last glacial cycle: transient simulations with an AOGCM

A number of transient climate runs simulating the last 120?kyr have been carried out using FAMOUS, a fast atmosphere–ocean general circulation model (AOGCM). This is the first time such experiments have been done with a full AOGCM, providing a three-dimensional simulation of both atmosphere and ocean over this period. Our simulation thus includes internally generated temporal variability over periods from days to millennia, and physical, detailed representations of important processes such as clouds and precipitation. Although the model is fast, computational restrictions mean that the rate of change of the forcings has been increased by a factor of 10, making each experiment 12?kyr long. Atmospheric greenhouse gases (GHGs), northern hemisphere ice sheets and variations in solar radiation arising from changes in the Earth’s orbit are treated as forcing factors, and are applied either separately or combined in different experiments. The long-term temperature changes on Antarctica match well with reconstructions derived from ice-core data, as does variability on timescales longer than 10 kyr. Last Glacial Maximum (LGM) cooling on Greenland is reasonably well simulated, although our simulations, which lack ice-sheet meltwater forcing, do not reproduce the abrupt, millennial scale climate shifts seen in northern hemisphere climate proxies or their slower southern hemisphere counterparts. The spatial pattern of sea surface cooling at the LGM matches proxy reconstructions reasonably well. There is significant anti-correlated variability in the strengths of the Atlantic meridional overturning circulation (AMOC) and the Antarctic Circumpolar Current (ACC) on timescales greater than 10?kyr in our experiments. We find that GHG forcing weakens the AMOC and strengthens the ACC, whilst the presence of northern hemisphere ice-sheets strengthens the AMOC and weakens the ACC. The structure of the AMOC at the LGM is found to be sensitive to the details of the ice-sheet reconstruction used. The precessional component of the orbital forcing induces ~20?kyr oscillations in the AMOC and ACC, whose amplitude is mediated by changes in the eccentricity of the Earth’s orbit. These forcing influences combine, to first order, in a linear fashion to produce the mean climate and ocean variability seen in the run with all forcings.     

Seawater density variations in the North Atlantic and the Atlantic meridional overturning circulation

Seawater property changes in the North Atlantic Ocean affect the Atlantic meridional overturning circulation (AMOC), which transports warm water northward from the upper ocean and contributes to the temperate climate of Europe, as well as influences climate globally. Previous observational studies have focused on salinity and freshwater variability in the sinking region of the North Atlantic, since it is believed that a freshening North Atlantic basin can slow down or halt the flow of the AMOC. Here we use available data to show the importance of how density patterns over the upper ocean of the North Atlantic affect the strength of the AMOC. For the long-term trend, the upper ocean of the subpolar North Atlantic is becoming cooler and fresher, whereas the subtropical North Atlantic is becoming warmer and saltier. On a multidecadal timescale, the upper ocean of the North Atlantic has generally been warmer and saltier since 1995. The heat and salt content in the subpolar North Atlantic lags that in the subtropical North Atlantic by about 8–9 years, suggesting a lower latitude origin for the temperature and salinity anomalies. Because of the opposite effects of temperature and salinity on density for both long-term trend and multidecadal timescales, these variations do not result in a density reduction in the subpolar North Atlantic for slowing down the AMOC. Indeed, the variations in the meridional density gradient between the subpolar and subtropical North Atlantic Ocean suggest that the AMOC has become stronger over the past five decades. These observed results are supported by and consistent with some oceanic reanalysis products.     

What causes the spread of model projections of ocean dynamic sea-level change in response to greenhouse gas forcing?

Sea levels of different atmosphere–ocean general circulation models (AOGCMs) respond to climate change forcing in different ways, representing a crucial uncertainty in climate change research. We isolate the role of the ocean dynamics in setting the spatial pattern of dynamic sea-level (ζ) change by forcing several AOGCMs with prescribed identical heat, momentum (wind) and freshwater flux perturbations. This method produces a ζ projection spread comparable in magnitude to the spread that results from greenhouse gas forcing, indicating that the differences in ocean model formulation are the cause, rather than diversity in surface flux change. The heat flux change drives most of the global pattern of ζ change, while the momentum and water flux changes cause locally confined features. North Atlantic heat uptake causes large temperature and salinity driven density changes, altering local ocean transport and ζ. The spread between AOGCMs here is caused largely by differences in their regional transport adjustment, which redistributes heat that was already in the ocean prior to perturbation. The geographic details of the ζ change in the North Atlantic are diverse across models, but the underlying dynamic change is similar. In contrast, the heat absorbed by the Southern Ocean does not strongly alter the vertically coherent circulation. The Arctic ζ change is dissimilar across models, owing to differences in passive heat uptake and circulation change. Only the Arctic is strongly affected by nonlinear interactions between the three air-sea flux changes, and these are model specific.