Role of marine biology in glacial-interglacial CO2 cycles

It has been hypothesized that changes in the marine biological pump caused a major portion of the glacial reduction of atmospheric carbon dioxide by 80 to 100 parts per million through increased iron fertilization of marine plankton, increased ocean nutrient content or utilization, or shifts in dominant plankton types. We analyze sedimentary records of marine productivity at the peak and the middle of the last glacial cycle and show that neither changes in nutrient utilization in the Southern Ocean nor shifts in plankton dominance explain the CO2 drawdown. Iron fertilization and associated mechanisms can be responsible for no more than half the observed drawdown.

The drivers of projected North Atlantic sea level change

Sea level change predicted by the CMIP5 atmosphere–ocean general circulation models (AOGCMs) is not spatially homogeneous. In particular, the sea level change in the North Atlantic is usually characterised by a meridional dipole pattern with higher sea level rise north of 40°N and lower to the south. The spread among models is also high in that region. Here we evaluate the role of surface buoyancy fluxes by carrying out simulations with the FAMOUS low-resolution AOGCM forced by surface freshwater and heat flux changes from CO2-forced climate change experiments with CMIP5 AOGCMs, and by a standard idealised surface freshwater flux applied in the North Atlantic. Both kinds of buoyancy flux change lead to the formation of the sea level dipole pattern, although the effect of the heat flux has a greater magnitude, and is the main cause of the spread of results among the CMIP5 models. By using passive tracers in FAMOUS to distinguish between additional and redistributed buoyancy, we show that the enhanced sea level rise north of 40°N is mainly due to the direct steric effect (the reduction of sea water density) caused by adding heat or freshwater locally. The surface buoyancy forcing also causes a weakening of the Atlantic meridional overturning circulation, and the consequent reduction of the northward ocean heat transport imposes a negative tendency on sea level rise, producing the reduced rise south of 40°N. However, unlike previous authors, we find that this indirect effect of buoyancy forcing is generally less important than the direct one, except in a narrow band along the east coast of the US, where it plays a major role and leads to sea level rise, as found by previous authors.

On the relationship between decadal buoyancy anomalies and variability of the Atlantic meridional overturning circulation

Owing to the role of the Atlantic meridional overturning circulation (AMOC) in ocean heat transport, AMOC variability is thought to play a role in climate variability on a wide range of time scales. This paper focuses on the potential role of the AMOC in climate variability on decadal time scales. Coupled and ocean-only general circulation models run in idealized geometries are utilized to study the relationships between decadal AMOC and buoyancy variability and determine whether the AMOC plays an active role in setting sea surface temperature on decadal time scales.DecadalAMOC variability is related to changes in the buoyancy field along the western boundary according to the thermal wind relation. Buoyancy anomalies originate in the upper ocean of the subpolar gyre and travel westward as baroclinic Rossby waves. When the buoyancy anomalies strike the western boundary, they are advected southward by the deep western boundary current, leading to latitudinally coherent AMOC variability. The AMOC is observed to respond passively to decadal buoyancy anomalies: although variability of the AMOC leads to meridional ocean heat transport anomalies, these transports are not responsible for creating the buoyancy anomalies in the subpolar gyre that drive AMOC variability.

Climate at high-obliquity

The question of climate at high obliquity is raised in the context of both exoplanet studies (e.g. habitability) and paleoclimates studies (evidence for low-latitude glaciation during the Neoproterozoic and the ”Snowball Earth” hypothesis). States of high obliquity, φ, are distinctive in that, for φ ≥54◦, the poles receive more solar radiation in the annual mean than the Equator, opposite to the present day situation. In addition, the seasonal cycle of insolation is extreme, with the poles alternatively “facing” the sun and sheltering in the dark for months. The novelty of our approach is to consider the role of a dynamical ocean in controlling the surface climate at high obliquity, which in turn requires understanding of the surface winds patterns when temperature gradients are reversed. To address these questions, a coupled ocean-atmosphere-sea ice GCM configured on an aquaplanet is employed. Except for the absence of topography and modified obliquity, the set-up is Earth-like. Two large obliquities φ, 54◦ and 90◦, are compared to today’s Earth value, φ=23.5◦. Three key results emerge at high obliquity: 1) despite reversed temper- ature gradients, mid-latitudes surface winds are westerly and trade winds exist at the equator (as for φ=23.5◦) although the westerlies are confined to the summer hemisphere, 2) a habitable planet is possible with mid-latitude temperatures in the range 300-280 K and 3) a stable climate state with an ice cap limited to the equatorial region is unlikely. We clarify the dynamics behind these features (notably by an analysis of the potential vorticity structure and conditions for baroclinic instability of the atmosphere). Interestingly, we find that the absence of a stable partially glaciated state is critically linked to the absence of ocean heat transport during winter, a feature ultimately traced back to the high seasonality of baroclinic instability conditions in the atmosphere.

An event-based approach to understanding decadal fluctuations in the Atlantic meridional overturning circulation

Many previous studies have shown that unforced climate model simulations exhibit decadal-scale fluctuations in the Atlantic meridional overturning circulation (AMOC), and that this variability can have impacts on surface climate fields. However, the robustness of these surface fingerprints across different models is less clear. Furthermore, with the potential for coupled feedbacks that may amplify or damp the response, it is not known whether the associated climate signals are linearly related to the strength of the AMOC changes, or if the fluctuation events exhibit nonlinear behaviour with respect to their strength or polarity. To explore these questions, we introduce an objective and flexible method for identifying the largest natural AMOC fluctuation events in multicentennial/multimillennial simulations of a variety of coupled climate models. The characteristics of the events are explored, including their magnitude, meridional coherence and spatial structure, as well as links with ocean heat transport and the horizontal circulation. The surface fingerprints in ocean temperature and salinity are examined, and compared with the results of linear regression analysis. It is found that the regressions generally provide a good indication of the surface changes associated with the largest AMOC events. However, there are some exceptions, including a nonlinear change in the atmospheric pressure signal, particularly at high latitudes, in HadCM3. Some asymmetries are also found between the changes associated with positive and negative AMOC events in the same model. Composite analysis suggests that there are signals that are robust across the largest AMOC events in each model, which provides reassurance that the surface changes associated with one particular event will be similar to those expected from regression analysis. However, large differences are found between the AMOC fingerprints in different models, which may hinder the prediction and attribution of such events in reality.

Broad range of 2050 warming from an observationally constrained large climate model ensemble

Incomplete understanding of three aspects of the climate system—equilibrium climate sensitivity, rate of ocean heat uptake and historical aerosol forcing—and the physical processes underlying them lead to uncertainties in our assessment of the global-mean temperature evolution in the twenty-first century1,2. Explorations of these uncertainties have so far relied on scaling approaches3,4, large ensembles of simplified climate models1,2, or small ensembles of complex coupled atmosphere–ocean general circulation models5,6 which under-represent uncertainties in key climate system properties derived from independent sources7–9. Here we present results from a multi-thousand-member perturbed-physics ensemble of transient coupled atmosphere–ocean general circulation model simulations. We find that model versions that reproduce observed surface temperature changes over the past 50 years show global-mean temperature increases of 1.4–3 K by 2050, relative to 1961–1990, under a mid-range forcing scenario. This range of warming is broadly consistent with the expert assessment provided by the Intergovernmental Panel on Climate Change Fourth Assessment Report10, but extends towards larger warming than observed in ensemblesof-opportunity5 typically used for climate impact assessments. From our simulations, we conclude that warming by the middle of the twenty-first century that is stronger than earlier estimates is consistent with recent observed temperature changes and a mid-range ‘no mitigation’ scenario for greenhouse-gas emissions.

Climate and carbon cycle changes from 1850 to 2100 in MPI-ESM simulations for the Coupled Model Intercomparison Project phase 5

The new Max-Planck-Institute Earth System Model (MPI-ESM) is used in the Coupled Model Intercomparison Project phase 5 (CMIP5) in a series of climate change experiments for either idealized CO2-only forcing or forcings based on observations and the Representative Concentration Pathway (RCP) scenarios. The paper gives an overview of the model configurations, experiments related forcings, and initialization procedures and presents results for the simulated changes in climate and carbon cycle. It is found that the climate feedback depends on the global warming and possibly the forcing history. The global warming from climatological 1850 conditions to 2080–2100 ranges from 1.5°C under the RCP2.6 scenario to 4.4°C under the RCP8.5 scenario. Over this range, the patterns of temperature and precipitation change are nearly independent of the global warming. The model shows a tendency to reduce the ocean heat uptake efficiency toward a warmer climate, and hence acceleration in warming in the later years. The precipitation sensitivity can be as high as 2.5% K−1 if the CO2 concentration is constant, or as small as 1.6% K−1, if the CO2 concentration is increasing. The oceanic uptake of anthropogenic carbon increases over time in all scenarios, being smallest in the experiment forced by RCP2.6 and largest in that for RCP8.5. The land also serves as a net carbon sink in all scenarios, predominantly in boreal regions. The strong tropical carbon sources found in the RCP2.6 and RCP8.5 experiments are almost absent in the RCP4.5 experiment, which can be explained by reforestation in the RCP4.5 scenario.

The inconstancy of the transient climate response parameter under increasing CO₂

In the Coupled Model Intercomparison Project Phase 5 (CMIP5), the model-mean increase in global mean surface air temperature T under the 1pctCO2 scenario (atmospheric CO2 increasing at 1% yr−1) during the second doubling of CO2 is 40% larger than the transient climate response (TCR), i.e. the increase in T during the first doubling. We identify four possible contributory effects. First, the surface climate system loses heat less readily into the ocean beneath as the latter warms. The model spread in the thermal coupling between the upper and deep ocean largely explains the model spread in ocean heat uptake efficiency. Second, CO2 radiative forcing may rise more rapidly than logarithmically with CO2 concentration. Third, the climate feedback parameter may decline as the CO2 concentration rises. With CMIP5 data, we cannot distinguish the second and third possibilities. Fourth, the climate feedback parameter declines as time passes or T rises; in 1pctCO2, this effect is less important than the others. We find that T projected for the end of the twenty-first century correlates more highly with T at the time of quadrupled CO2 in 1pctCO2 than with the TCR, and we suggest that the TCR may be underestimated from observed climate change.

Recent progress in understanding and projecting regional and global mean sea-level change

Considerable progress has been made in understanding the present and future regional and global sea level in the 2 years since the publication of the Fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change. Here, we evaluate how the new results affect the AR5’s assessment of (i) historical sea level rise, including attribution of that rise and implications for the sea level budget, (ii) projections of the components and of total global mean sea level (GMSL), and (iii) projections of regional variability and emergence of the anthropogenic signal. In each of these cases, new work largely provides additional evidence in support of the AR5 assessment, providing greater confidence in those findings. Recent analyses confirm the twentieth century sea level rise, with some analyses showing a slightly smaller rate before 1990 and some a slightly larger value than reported in the AR5. There is now more evidence of an acceleration in the rate of rise. Ongoing ocean heat uptake and associated thermal expansion have continued since 2000, and are consistent with ocean thermal expansion reported in the AR5. A significant amount of heat is being stored deeper in the water column, with a larger rate of heat uptake since 2000 compared to the previous decades and with the largest storage in the Southern Ocean. The first formal detection studies for ocean thermal expansion and glacier mass loss since the AR5 have confirmed the AR5 finding of a significant anthropogenic contribution to sea level rise over the last 50 years. New projections of glacier loss from two regions suggest smaller contributions to GMSL rise from these regions than in studies assessed by the AR5; additional regional studies are required to further assess whether there are broader implications of these results. Mass loss from the Greenland Ice Sheet, primarily as a result of increased surface melting, and from the Antarctic Ice Sheet, primarily as a result of increased ice discharge, has accelerated. The largest estimates of acceleration in mass loss from the two ice sheets for 2003–2013 equal or exceed the acceleration of GMSL rise calculated from the satellite altimeter sea level record over the longer period of 1993–2014. However, when increased mass gain in land water storage and parts of East Antarctica, and decreased mass loss from glaciers in Alaska and some other regions are taken into account, the net acceleration in the ocean mass gain is consistent with the satellite altimeter record. New studies suggest that a marine ice sheet instability (MISI) may have been initiated in parts of the West Antarctic Ice Sheet (WAIS), but that it will affect only a limited number of ice streams in the twenty-first century. New projections of mass loss from the Greenland and Antarctic Ice Sheets by 2100, including a contribution from parts of WAIS undergoing unstable retreat, suggest a contribution that falls largely within the likely range (i.e., two thirds probability) of the AR5. These new results increase confidence in the AR5 likely range, indicating that there is a greater probability that sea level rise by 2100 will lie in this range with a corresponding decrease in the likelihood of an additional contribution of several tens of centimeters above the likely range. In view of the comparatively limited state of knowledge and understanding of rapid ice sheet dynamics, we continue to think that it is not yet possible to make reliable quantitative estimates of future GMSL rise outside the likely range. Projections of twenty-first century GMSL rise published since the AR5 depend on results from expert elicitation, but we have low confidence in conclusions based on these approaches. New work on regional projections and emergence of the anthropogenic signal suggests that the two commonly predicted features of future regional sea level change (the increasing tilt across the Antarctic Circumpolar Current and the dipole in the North Atlantic) are related to regional changes in wind stress and surface heat flux. Moreover, it is expected that sea level change in response to anthropogenic forcing, particularly in regions of relatively low unforced variability such as the low-latitude Atlantic, will be detectable over most of the ocean by 2040. The east-west contrast of sea level trends in the Pacific observed since the early 1990s cannot be satisfactorily accounted for by climate models, nor yet definitively attributed either to unforced variability or forced climate change.

Assimilation impacts on Arctic Ocean circulation, heat and freshwater

We investigate the Arctic basin circulation, freshwater content (FWC) and heat budget by using a high-resolution global coupled ice–ocean model implemented with a state-of-the-art data assimilation scheme. We demonstrate that, despite a very sparse dataset, by assimilating hydrographic data in and near the Arctic basin, the initial warm bias and drift in the control run is successfully corrected, reproducing a much more realistic vertical and horizontal structure to the cyclonic boundary current carrying the Atlantic Water (AW) along the Siberian shelves in the reanalysis run. The Beaufort Gyre structure and FWC and variability are also more accurately reproduced. Small but important changes in the strait exchange flows are found which lead to more balanced budgets in the reanalysis run. Assimilation fluxes dominate the basin budgets over the first 10 years (P1: 1987–1996) of the reanalysis for both heat and FWC, after which the drifting Arctic upper water properties have been restored to realistic values. For the later period (P2: 1997–2004), the Arctic heat budget is almost balanced without assimilation contributions, while the freshwater budget shows reduced assimilation contributions compensating largely for surface salinity damping, which was extremely strong in this run. A downward trend in freshwater export at the Canadian Straits and Fram Strait is found in period P2, associated with Beaufort Gyre recharge. A detailed comparison with observations and previous model studies at the individual Arctic straits is also included.

Global climate and ocean circulation on an aquaplanet ocean-atmosphere general circulation model

A low resolution coupled ocean-atmosphere general circulation model OAGCM is used to study the characteristics of the large scale ocean circulation and its climatic impacts in a series of global coupled aquaplanet experiments. Three configurations, designed to produce fundamentally different ocean circulation regimes, are considered. The first has no obstruction to zonal flow, the second contains a low barrier that blocks zonal flow in the ocean at all latitudes, creating a single enclosed basin, whilst the third contains a gap in the barrier to allow circumglobal flow at high southern latitudes. Warm greenhouse climates with a global average air surface temperature of around 27C result in all cases. Equator to pole temperature gradients are shallower than that of a current climate simulation. Whilst changes in the land configuration cause regional changes in temperature, winds and rainfall, heat transports within the system are little affected. Inhibition of all ocean transport on the aquaplanet leads to a reduction in global mean surface temperature of 8C, along with a sharpening of the meridional temperature gradient. This results from a reduction in global atmospheric water vapour content and an increase in tropical albedo, both of which act to reduce global surface temperatures. Fitting a simple radiative model to the atmospheric characteristics of the OAGCM solutions suggests that a simpler atmosphere model, with radiative parameters chosen a priori based on the changing surface configuration, would have produced qualitatively different results. This implies that studies with reduced complexity atmospheres need to be guided by more complex OAGCM results on a case by case basis.

Vertical heat transports in the ocean and their effect on time-dependent climate change

In response to increasing atmospheric con- centrations of greenhouse gases, the rate of time- dependent climate change is determined jointly by the strength of climate feedbacks and the e�ciency of pro- cesses which remove heat from the surface into the deep ocean. This work examines the vertical heat transport processes in the ocean of the HADCM2 atmosphere± ocean general circulation model (AOGCM) in experi- ments with CO2 held constant (control) and increasing at 1% per year (anomaly). The control experiment shows that global average heat exchanges between the upper and lower ocean are dominated by the Southern Ocean, where heat is pumped downwards by the wind- driven circulation and di�uses upwards along sloping isopycnals. This is the reverse of the low-latitude balance used in upwelling±di�usion ocean models, the global average upward di�usive transport being against the temperature gradient. In the anomaly experiment, weakened convection at high latitudes leads to reduced diffusive and convective heat loss from the deep ocean, and hence to net heat uptake, since the advective heat input is less a�ected. Reduction of deep water produc- tion at high latitudes results in reduced upwelling of cold water at low latitudes, giving a further contribution to net heat uptake. On the global average, high-latitude processes thus have a controlling in¯uence. The impor- tant role of di�usion highlights the need to ensure that the schemes employed in AOGCMs give an accurate representation of the relevant sub-grid-scale processes.

Atlantic meridional heat transports in two ocean reanalyes evaluated against the RAPID array

The Atlantic meridional overturning circulation in two versions of the NEMO ¼° global ocean model has been compared with the RAPID transport array at 26oN. Both model versions reproduce the mean MOC strength well although the Florida Straits flows differ because the pathway of the Gulf Stream is not strongly constrained at this resolution. Both models however have a mean meridional heat transport of 1.07PW, much lower than the 1.35PW from RAPID observations in Apr04-Oct07. Much of the heat transport discrepancy is due to lower transports in summer across the MidOcean (Bahamas-Africa) section, due to stronger southward geostrophic flows in the top 100m where the water is warmest. Seasonal thermocline changes increase temperature differences across the basin driving stronger geostrophic shear, but this effect is much weaker in the top 100m of the RAPID velocity data. The effect accounts for a reduction of 1.1Sv in MOC and 0.1PW in heat transports. The rest of the discrepancy comes from lower Ekman transports from using ERAInterim winds instead of QuikSCAT, a smaller zonally-varying “Eddy” heat transport component, estimated from repeat XBT sections in the observations, and the southward throughflow in the model. Other differences in depth structure of the model MOC and RAPID observations are described but have much less impact on heat transports.

Freshwater and heat transports from global ocean synthesis

[1] An eddy-permitting ¼° global ocean reanalysis based on the Operational Met Office FOAM data assimilation system has been run for 1989–2010 forced by ERA-Interim meteorology. Freshwater and heat transports are compared with published estimates globally and in each basin, with special focus on the Atlantic. The meridional transports agree with observations within errors at most locations, but where eddies are active the transports by the mean flow are nearly always in better agreement than the total transports. Eddy transports are down gradient and are enhanced relative to a free run. They may oppose or reinforce mean transports and provide 40–50% of the total transport near midlatitude fronts, where eddies with time scales <1 month provide up to 15%. Basin-scale freshwater convergences are calculated with the Arctic/Atlantic, Indian, and Pacific oceans north of 32°S, all implying net evaporation of 0.33 ± 0.04 Sv, 0.65 ± 0.07 Sv, and 0.09 ± 0.04 Sv, respectively, within the uncertainty of observations in the Atlantic and Pacific. The Indian is more evaporative and the Southern Ocean has more precipitation (1.07 Sv). Air-sea fluxes are modified by assimilation influencing turbulent heat fluxes and evaporation. Generally, surface and assimilation fluxes together match the meridional transports, indicating that the reanalysis is close to a steady state. Atlantic overturning and gyre transports are assessed with overturning freshwater transports southward at all latitudes. At 26°N eddy transports are negligible, overturning transport is 0.67 ± 0.19 Sv southward and gyre transport is 0.44 ± 0.17 Sv northward, with divergence between 26°N and the Bering Strait of 0.13 ± 0.23 Sv over 2004–2010.