ALBERT

Description: The aim of our study was to use the coupled climate-carbon cycle model HadCM3C to quantify climate impact of ecosystem changes over recent decades and under future scenarios, due to changes in both atmospheric CO2 and surface albedo. We use two future scenarios – the IPCC SRES A1B scenario, and a climate stabilisation scenario (2C20), allowing us to assess the impact of climate mitigation on results. We performed a pair of simulations under each scenario – one in which vegetation was fixed at the initial state and one in which vegetation changes dynamically in response to climate change, as determined by the interactive vegetation model within HadCM3C. In our simulations with interactive vegetation, relatively small changes in global vegetation coverage were found, mainly dominated by increases in scrub and needleleaf trees at high latitudes and losses of broadleaf trees and grasses across the Amazon. Globally this led to a loss of terrestrial carbon, mainly from the soil. Global changes in carbon storage were related to the regional losses from the Amazon and gains at high latitude. Regional differences in carbon storage between the two scenarios were largely driven by the balance between warming-enhanced decomposition and altered vegetation growth. Globally, interactive vegetation reduced albedo acting to enhance albedo changes due to climate change. This was mainly related to the darker land surface over high latitudes (due to vegetation expansion, particularly during winter and spring); small increases in albedo occurred over the Amazon. As a result, there was a relatively small impact of vegetation change on most global annual mean climate variables, which was generally greater under A1B than 2C20, with markedly stronger local-to-regional and seasonal impacts. Globally, vegetation change amplified future annual temperature increases by 0.24 and 0.15 K (under A1B and 2C20, respectively) and increased global precipitation, with reductions in precipitation over the Amazon and increases over high latitudes. In general, changes were stronger over land – for example, global temperature changes due to interactive vegetation of 0.43 and 0.28 K under A1B and 2C20, respectively. Regionally, the warming influence of future vegetation change in our simulations was driven by the balance between driving factors. For instance, reduced tree cover over the Amazon reduced evaporation (particularly during summer), outweighing the cooling influence of any small albedo changes. In contrast, at high latitudes the warming impact of reduced albedo (particularly during winter and spring) due to increased vegetation cover appears to have offset any cooling due to small evaporation increases. Climate mitigation generally reduced the impact of vegetation change on future global and regional climate in our simulations. Our study therefore suggests that there is a need to consider both biogeochemical and biophysical effects in climate adaptation and mitigation decision making.

Description: The oceans currently take up around a quarter of the carbon dioxide (CO2) emitted by human activity. While stored in the ocean, this CO2 is not influencing Earth's radiation budget; the ocean CO2 sink therefore plays an important role in mitigating global warming. CO2 uptake by the oceans is heterogeneous, with the subpolar North Atlantic being the strongest CO2 sink region. Observations over the last two decades have indicated that CO2 uptake by the subpolar North Atlantic sink can vary rapidly. Given the importance of this sink and its apparent variability, it is critical that we understand the mechanisms behind its operation. Here we explore subpolar North Atlantic CO2 uptake across a large ensemble of Earth System Model simulations, and find that models show a peak in sink strength around the middle of the century after which CO2 uptake begins to decline. We identify different drivers of change on interannual and multidecadal timescales. Short-term variability appears to be driven by fluctuations in regional seawater temperature and alkalinity, whereas the longer-term evolution throughout the coming century is largely occurring through a counterintuitive response to rising atmospheric CO2 concentrations. At high atmospheric CO2 concentrations the contrasting Ravelle factors between the subtropical and subpolar gyres, combined with the transport of surface waters from the subtropical to subpolar gyre, means that the subpolar CO2 uptake capacity is largely satisfied from its southern boundary rather than through air–sea CO2 flux. Our findings indicate that: (i) we can explain the mechanisms of subpolar North Atlantic CO2 uptake variability across a broad range of Earth System Models, (ii) a focus on understanding the mechanisms behind contemporary variability may not directly tell us about how the sink will change in the future, (iii) to identify long-term change in the North Atlantic CO2 sink we should focus observational resources on monitoring subtropical as well as the subpolar seawater CO2, (iv) recent observations of a weakening subpolar North Atlantic CO2 sink suggests that the sink strength is already in long-term decline.

Description: A common time scale for the EPICA ice cores from Dome C (EDC) and Dronning Maud Land (EDML) has been established. Since the EDML core was not drilled on a dome, the development of the EDML1 time scale for the EPICA ice core drilled in Dronning Maud Land was based on the creation of a detailed stratigraphic link between EDML and EDC, which was dated by a simpler 1D ice-flow model. The synchronisation between the two EPICA ice cores was done through the identification of several common volcanic signatures. This paper describes the rigorous method, using the signature of volcanic sulfate, which was employed for the last 52 kyr of the record. We estimated the discrepancies between the modelled EDC and EDML glaciological age scales during the studied period, by evaluating the ratio R of the apparent duration of temporal intervals between pairs of isochrones. On average R ranges between 0.8 and 1.2 corresponding to an uncertainty of up to 20% in the estimate of the time duration in at least one of the two ice cores. Significant deviations of R up to 1.4–1.5 are observed between 18 and 28 kyr before present (BP), where present is defined as 1950. At this stage our approach does not allow us unequivocally to find out which of the models is affected by errors, but assuming that the thinning function at both sites and accumulation history at Dome C (which was drilled on a dome) are correct, this anomaly can be ascribed to a complex spatial accumulation variability (which may be different in the past compared to the present day) upstream of the EDML core.

Description: The El Niño-Southern Oscillation (ENSO), an ocean–atmosphere coupled oscillation over the equatorial Pacific, is the most important source of global climate variability on inter-annual time scales. It has substantial environmental and socio-economic consequences such as devastation of South American fish populations and increased forest fires in Indonesia. The instrumental ENSO record is too short for analysing long-term trends and variability, hence proxy data is used to extend the record. However, different proxy sources have produced varying reconstructions of ENSO, with some evidence for a temperature–precipitation divergence in ENSO trends over the past millennium, in particular during the Mediaeval Climate Anomaly (MCA; AD 800–1300) and the Little Ice Age (LIA; AD 1400–1850). This throws into question the stability of the modern ENSO system and its links to the global climate, which has implications for future projections. Here we use a new statistical approach using EOF-based weighting to create two new large-scale ENSO reconstructions derived independently from precipitation proxies and temperature proxies respectively. The method is developed and validated using pseudoproxy experiments that address the effects of proxy dating error, resolution and noise to improve uncertainty estimations. The precipitation ENSO reconstruction displays a significantly more El Niño-like state during the LIA than the MCA, while the temperature reconstruction shows no significant difference. The trends shown in the precipitation ENSO reconstruction are relatively robust to variations in the precipitation EOF pattern. However, the temperature reconstruction suffers significantly from a lack of high-quality, favourably located proxy records, which limits its ability to capture the large-scale ENSO signal. Further expansion of the palaeo-database and improvements to instrumental, satellite and model representations of ENSO are needed to fully resolve the discrepancies found among proxy records.

Description: Ice core data from Antarctica provide detailed insights into the characteristics of past climate, atmospheric circulation, as well as changes in the aerosol load of the atmosphere. We present high-resolution records of soluble calcium (Ca2+), non-sea-salt soluble calcium (nssCa2+), and particulate mineral dust aerosol from the East Antarctic Plateau at a depth resolution of 1 cm, spanning the past 800 000 years. Despite the fact that all three parameters are largely dust-derived, the ratio of nssCa2+ to particulate dust is dependent on the particulate dust concentration itself. We used principal component analysis to extract the joint climatic signal and produce a common high-resolution record of dust flux. This new record is used to identify Antarctic warming events during the past eight glacial periods. The phasing of dust flux and CO2 changes during glacial-interglacial transitions reveals that iron fertilization of the Southern Ocean during the past nine glacial terminations was not the dominant factor in the deglacial rise of CO2 concentrations. Rapid changes in dust flux during glacial terminations and Antarctic warming events point to a rapid response of the southern westerly wind belt in the region of southern South American dust sources on changing climate conditions. The clear lead of these dust changes on temperature rise suggests that an atmospheric reorganization occurred in the Southern Hemisphere before the Southern Ocean warmed significantly.

Description: The many different proxy records from the European Project for Ice Coring in Antarctica (EPICA) Dome C ice core allow for the first time a comparison of nine glacial terminations in great detail. Despite the fact that all terminations cover the transition from a glacial maximum into an interglacial, there are large differences between single terminations. For some terminations, Antarctic temperature increased only moderately, while for others, the amplitude of change at the termination was much larger. For the different terminations, the rate of change in temperature is more similar than the magnitude or duration of change. These temperature changes were accompanied by vast changes in dust and sea salt deposition all over Antarctica. Here we investigate the phasing between a South American dust proxy (non-sea-salt calcium flux, nssCa2+), a sea ice proxy (sea salt sodium flux, ssNa+) and a proxy for Antarctic temperature (deuterium, δD). In particular, we look into whether a similar sequence of events applies to all terminations, despite their different characteristics. All proxies are derived from the EPICA Dome C ice core, resulting in a relative dating uncertainty between the proxies of less than 20 years. At the start of the terminations, the temperature (δD) increase and dust (nssCa2+ flux) decrease start synchronously. The sea ice proxy (ssNa+ flux), however, only changes once the temperature has reached a particular threshold, approximately 5°C below present day temperatures (corresponding to a δD value of −420‰). This reflects to a large extent the limited sensitivity of the sea ice proxy during very cold periods with large sea ice extent. At terminations where this threshold is not reached (TVI, TVIII), ssNa+ flux shows no changes. Above this threshold, the sea ice proxy is closely coupled to the Antarctic temperature, and interglacial levels are reached at the same time for both ssNa+ and δD. On the other hand, once another threshold at approximately 2°C below present day temperature is passed (corresponding to a δD value of −402‰), nssCa2+ flux has reached interglacial levels and does not change any more, despite further warming. This threshold behaviour most likely results from a combination of changes to the threshold friction velocity for dust entrainment and to the distribution of surface wind speeds in the dust source region.

Description: The many different proxy records from the European Project for Ice Coring in Antarctica (EPICA) Dome C ice core allow for the first time a comparison of nine glacial terminations in great detail. Despite the fact that all terminations cover the transition from a glacial maximum into an interglacial, there are large differences between single terminations. For some terminations, Antarctic temperature increased only moderately, while for others, the amplitude of change at the termination was much larger. For the different terminations, the rate of change in temperature is more similar than the magnitude or duration of change. These temperature changes were accompanied by vast changes in dust and sea salt deposition all over Antarctica. Here we investigate the phasing between a South American dust proxy (non-sea-salt calcium flux, nssCa), a sea ice proxy (sea salt sodium flux, ssNa) and a proxy for Antarctic temperature (deuterium, δD). In particular, we look into whether a similar sequence of events applies to all terminations, despite their different characteristics. All proxies are derived from the EPICA Dome C ice core, resulting in a relative dating uncertainty between the proxies of less than 20 years. At the start of the terminations, the temperature (δD) increase and dust (nssCa flux) decrease start synchronously. The sea ice proxy (ssNa flux), however, only changes once the temperature has reached a particular threshold, approximately 5°C below present day temperatures (corresponding to a δD value of –420‰). This reflects to a large extent the limited sensitivity of the sea ice proxy during very cold periods with large sea ice extent. At terminations where this threshold is not reached (TVI, TVIII), ssNa flux shows no changes. Above this threshold, the sea ice proxy is closely coupled to the Antarctic temperature, and interglacial levels are reached at the same time for both ssNa and δD. On the other hand, once another threshold at approximately 2°C below present day temperature is passed (corresponding to a δD value of –402‰), nssCa flux has reached interglacial levels and does not change any more, despite further warming. This threshold behaviour most likely results from a combination of changes to the threshold friction velocity for dust entrainment and to the distribution of surface wind speeds in the dust source region.