The impact of 21st century Arctic sea ice decline on the climate and mass balance of the Greenland ice sheet and Svalbard’s glaciers and ice caps

The Arctic is a region particularly susceptible to rapid climate change. General circulation models (GCMs) suggest a polar amplification of any global warming signal by a factor of about 1.5 due, in part, to sea ice feedbacks. The dramatic recent decline in multi-year sea ice cover lies outside the standard deviation of the CMIP3 ensemble GCM predictions. Sea ice acts as a barrier between cold air and warmer oceans during winter, as well as inhibiting evaporation from the ocean surface water during the summer. An ice free Arctic would likely have an altered hydrological cycle with more evaporation from the ocean surface leading to changes in precipitation distribution and amount. Using the U.K. Met Office Regional Climate Model (RCM), HadRM3, the atmospheric effects of the observed and projected reduction in Arctic sea ice are investigated. The RCM is driven by the atmospheric GCM HadAM3. Both models are forced with sea surface temperature and sea ice for the period 2061-2090 from the CMIP3 HadGEM1 experiments. Here we use an RCM at 50km resolution over the Arctic and 25km over Svalbard, which captures well the present-day pattern of precipitation and provides a detailed picture of the projected changes in the behaviour of the oceanic-atmosphere moisture fluxes and how they affect precipitation. These experiments show that the projected 21stCentury sea ice decline alone causes large impacts to the surface mass balance (SMB) on Svalbard. However Greenland’s SMB is not significantly affected by sea ice decline alone, but responds with a strongly negative shift in SMB when changes to SST are incorporated into the experiments. This is the first study to characterise the impact of changes in future sea ice to Arctic terrestrial cryosphere mass balance.

Changes in area and geodetic mass balance of small glaciers, polar Urals, Russia, 1950-2008

Changes in area of 30 small glaciers (mostly <1 km2) in the northern Polar Urals (67.5-68.25 °N) between 1953 and 2000 were assessed using historic aerial photography from 1953 and 1960, ASTER and panchromatic Landsat ETM+ imagery from 2000, and data from 1981 and 2008 terrestrial surveys. Changes in volume and geodetic mass balance of IGAN and Obruchev glaciers were calculated using data from terrestrial surveys in 1963 and 2008. In total, glacier area declined by 22.3 ± 3.9% in the 1953/60-2000 period. The areas of individual glaciers decreased by 4-46%. Surfaces of Obruchev and IGAN glaciers lowered by 22.5 ± 1.7 m and 14.9 ± 2.1 m. Over 45 years, geodetic mass balances of Obruchev and IGAN glaciers were -20.66 ± 2.91 and -13.54 ± 2.57 m w.e. respectively. Glacier shrinkage in the Polar Urals is related to a summer warming of 1 °C between 1953-81 and 1981-2008 and its rates are consistent with other regions of northern Asia but are higher than in Scandinavia. While glacier shrinkage intensified in the 1981-2000 period relative to 1953-81, increasing winter precipitation and shading effects slowed glacier wastage in 2000-08.

The REFLEX project: Comparing different algorithms and implementations for the inversion of a terrestrial ecosystem model against eddy covariance data

We describe a model-data fusion (MDF) inter-comparison project (REFLEX), which compared various algorithms for estimating carbon (C) model parameters consistent with both measured carbon fluxes and states and a simple C model. Participants were provided with the model and with both synthetic net ecosystem exchange (NEE) of CO2 and leaf area index (LAI) data, generated from the model with added noise, and observed NEE and LAI data from two eddy covariance sites. Participants endeavoured to estimate model parameters and states consistent with the model for all cases over the two years for which data were provided, and generate predictions for one additional year without observations. Nine participants contributed results using Metropolis algorithms, Kalman filters and a genetic algorithm. For the synthetic data case, parameter estimates compared well with the true values. The results of the analyses indicated that parameters linked directly to gross primary production (GPP) and ecosystem respiration, such as those related to foliage allocation and turnover, or temperature sensitivity of heterotrophic respiration, were best constrained and characterised. Poorly estimated parameters were those related to the allocation to and turnover of fine root/wood pools. Estimates of confidence intervals varied among algorithms, but several algorithms successfully located the true values of annual fluxes from synthetic experiments within relatively narrow 90% confidence intervals, achieving >80% success rate and mean NEE confidence intervals <110 gC m−2 year−1 for the synthetic case. Annual C flux estimates generated by participants generally agreed with gap-filling approaches using half-hourly data. The estimation of ecosystem respiration and GPP through MDF agreed well with outputs from partitioning studies using half-hourly data. Confidence limits on annual NEE increased by an average of 88% in the prediction year compared to the previous year, when data were available. Confidence intervals on annual NEE increased by 30% when observed data were used instead of synthetic data, reflecting and quantifying the addition of model error. Finally, our analyses indicated that incorporating additional constraints, using data on C pools (wood, soil and fine roots) would help to reduce uncertainties for model parameters poorly served by eddy covariance data.

Terrestrial biogeochemical feedbacks in the climate system

The terrestrial biosphere is a key regulator of atmospheric chemistry and climate. During past periods of climate change, vegetation cover and interactions between the terrestrial biosphere and atmosphere changed within decades. Modern observations show a similar responsiveness of terrestrial biogeochemistry to anthropogenically forced climate change and air pollution. Although interactions between the carbon cycle and climate have been a central focus, other biogeochemical feedbacks could be as important in modulating future climate change. Total positive radiative forcings resulting from feedbacks between the terrestrial biosphere and the atmosphere are estimated to reach up to 0.9 or 1.5 W m−2 K−1 towards the end of the twenty-first century, depending on the extent to which interactions with the nitrogen cycle stimulate or limit carbon sequestration. This substantially reduces and potentially even eliminates the cooling effect owing to carbon dioxide fertilization of the terrestrial biota. The overall magnitude of the biogeochemical feedbacks could potentially be similar to that of feedbacks in the physical climate system, but there are large uncertainties in the magnitude of individual estimates and in accounting for synergies between these effects.

How will organic carbon stocks in mineral soils evolve under future climate? Global projections using RothC for a range of climate change scenarios

We use a soil carbon (C) model (RothC), driven by a range of climate models for a range of climate scenarios to examine the impacts of future climate on global soil organic carbon (SOC) stocks. The results suggest an overall global increase in SOC stocks by 2100 under all scenarios, but with a different extent of increase among the climate model and emissions scenarios. The impacts of projected land use changes are also simulated, but have relatively minor impacts at the global scale. Whether soils gain or lose SOC depends upon the balance between C inputs and decomposition. Changes in net primary production (NPP) change C inputs to the soil, whilst decomposition usually increases under warmer temperatures, but can also be slowed by decreased soil moisture. Underlying the global trend of increasing SOC under future climate is a complex pattern of regional SOC change. SOC losses are projected to occur in northern latitudes where higher SOC decomposition rates due to higher temperatures are not balanced by increased NPP, whereas in tropical regions, NPP increases override losses due to higher SOC decomposition. The spatial heterogeneity in the response of SOC to changing climate shows how delicately balanced the competing gain and loss processes are, with subtle changes in temperature, moisture, soil type and land use, interacting to determine whether SOC increases or decreases in the future. Our results suggest that we should stop looking for a single answer regarding whether SOC stocks will increase or decrease under future climate, since there is no single answer. Instead, we should focus on improving our prediction of the factors that determine the size and direction of change, and the land management practices that can be implemented to protect and enhance SOC stocks.

A Possible Change in Mass Balance of Greenland and Antarctic Ice Sheets in the Coming Century

A high-resolution GCM is found to simulate precipitation and surface energy balance of high latitudes with high accuracy. This opens new possibilities to investigate the future mass balance of polar glaciers and its effect on sea level. The surface mass balance of the Greenland and the Antarctic ice sheets is simulated using the ECHAM3 GCM with TI06 horizontal resolution. With this model, two 5-year integrations for the present and doubled carbon dioxide conditions based on the boundary conditions provided by the ECHAM1/T21 transient experiment have been conducted. A comparison of the two experiments over Greenland and Antarctica shows to what extent the effect of climate change on the mass balance on the two largest glaciers of the world can differ. On Greenland one sees a slight decrease in accumulation and a substantial increase in melt, while on Antarctica a large increase in accumulation without melt is projected. Translating the mass balances into terms of sea-level equivalent. the Greenland discharge causes a sea level rise of 1.1 mm yr−1, while the accumulation on Antarctica tends to lower it by 0.9 mm yr−1. The change in the combined mass balance of the two continents is almost zero. The sea level change of the next century can be affected more effectively by the thermal expansion of seawater and the mass balance of smaller glaciers outside of Greenland and Antarctica.

Stable isotope and modelling evidence for CO2 as a driver of glacial–interglacial vegetation shifts in southern Africa

Atmospheric CO2 concentration is hypothesized to influence vegetation distribution via tree–grass competition, with higher CO2 concentrations favouring trees. The stable carbon isotope (δ13C) signature of vegetation is influenced by the relative importance of C4 plants (including most tropical grasses) and C3 plants (including nearly all trees), and the degree of stomatal closure – a response to aridity – in C3 plants. Compound-specific δ13C analyses of leaf-wax biomarkers in sediment cores of an offshore South Atlantic transect are used here as a record of vegetation changes in subequatorial Africa. These data suggest a large increase in C3 relative to C4 plant dominance after the Last Glacial Maximum. Using a process-based biogeography model that explicitly simulates 13C discrimination, it is shown that precipitation and temperature changes cannot explain the observed shift in δ13C values. The physiological effect of increasing CO2 concentration is decisive, altering the C3/C4 balance and bringing the simulated and observed δ13C values into line. It is concluded that CO2 concentration itself was a key agent of vegetation change in tropical southern Africa during the last glacial–interglacial transition. Two additional inferences follow. First, long-term variations in terrestrial δ13Cvalues are not simply a proxy for regional rainfall, as has sometimes been assumed. Although precipitation and temperature changes have had major effects on vegetation in many regions of the world during the period between the Last Glacial Maximum and recent times, CO2 effects must also be taken into account, especially when reconstructing changes in climate between glacial and interglacial states. Second, rising CO2 concentration today is likely to be influencing tree–grass competition in a similar way, and thus contributing to the "woody thickening" observed in savannas worldwide. This second inference points to the importance of experiments to determine how vegetation composition in savannas is likely to be influenced by the continuing rise of CO2 concentration.

Impact of brine-induced stratification on the glacial carbon cycle

During the cold period of the Last Glacial Maximum (LGM, about 21 000 years ago) atmospheric CO2 was around 190 ppm, much lower than the pre-industrial concentration of 280 ppm. The causes of this substantial drop remain partially unresolved, despite intense research. Understanding the origin of reduced atmospheric CO2 during glacial times is crucial to comprehend the evolution of the different carbon reservoirs within the Earth system (atmosphere, terrestrial biosphere and ocean). In this context, the ocean is believed to play a major role as it can store large amounts of carbon, especially in the abyss, which is a carbon reservoir that is thought to have expanded during glacial times. To create this larger reservoir, one possible mechanism is to produce very dense glacial waters, thereby stratifying the deep ocean and reducing the carbon exchange between the deep and upper ocean. The existence of such very dense waters has been inferred in the LGM deep Atlantic from sediment pore water salinity and δ18O inferred temperature. Based on these observations, we study the impact of a brine mechanism on the glacial carbon cycle. This mechanism relies on the formation and rapid sinking of brines, very salty water released during sea ice formation, which brings salty dense water down to the bottom of the ocean. It provides two major features: a direct link from the surface to the deep ocean along with an efficient way of setting a strong stratification. We show with the CLIMBER-2 carbon-climate model that such a brine mechanism can account for a significant decrease in atmospheric CO2 and contribute to the glacial-interglacial change. This mechanism can be amplified by low vertical diffusion resulting from the brine-induced stratification. The modeled glacial distribution of oceanic δ13C as well as the deep ocean salinity are substantially improved and better agree with reconstructions from sediment cores, suggesting that such a mechanism could have played an important role during glacial times.

Systematic study of the impact of fresh water fluxes on the glacial carbon cycle

During glacial periods, atmospheric CO2 concentration increases and decreases by around 15 ppm. At the same time, the climate changes gradually in Antarctica. Such climate changes can be simulated in models when the AMOC (Atlantic Meridional Oceanic Circulation) is weakened by adding fresh water to the North Atlantic. The impact on the carbon cycle is less straightforward, and previous studies give opposite results. Because the models and the fresh water fluxes were different in these studies, it prevents any direct comparison and hinders finding whether the discrepancies arise from using different models or different fresh water fluxes. In this study we use the CLIMBER-2 coupled climate carbon model to explore the impact of different fresh water fluxes. In both preindustrial and glacial states, the addition of fresh water and the resulting slow-down of the AMOC lead to an uptake of carbon by the ocean and a release by the terrestrial biosphere. The duration, shape and amplitude of the fresh water flux all have an impact on the change of atmospheric CO2 because they modulate the change of the AMOC. The maximum CO2 change linearly depends on the time integral of the AMOC change. The different duration, amplitude, and shape of the fresh water flux cannot explain the opposite evolution of ocean and vegetation carbon inventory in different models. The different CO2 evolution thus depends on the AMOC response to the addition of fresh water and the resulting climatic change, which are both model dependent. In CLIMBER-2, the rise of CO2 recorded in ice cores during abrupt events can be simulated under glacial conditions, especially when the sinking of brines in the Southern Ocean is taken into account. The addition of fresh water in the Southern Hemisphere leads to a decline of CO2, contrary to the addition of fresh water in the Northern Hemisphere.

Impact of oceanic processes on the carbon cycle during the last termination

During the last termination (from ~18 000 years ago to ~9000 years ago), the climate significantly warmed and the ice sheets melted. Simultaneously, atmospheric CO2 increased from ~190 ppm to ~260 ppm. Although this CO2 rise plays an important role in the deglacial warming, the reasons for its evolution are difficult to explain. Only box models have been used to run transient simulations of this carbon cycle transition, but by forcing the model with data constrained scenarios of the evolution of temperature, sea level, sea ice, NADW formation, Southern Ocean vertical mixing and biological carbon pump. More complex models (including GCMs) have investigated some of these mechanisms but they have only been used to try and explain LGM versus present day steady-state climates. In this study we use a coupled climate-carbon model of intermediate complexity to explore the role of three oceanic processes in transient simulations: the sinking of brines, stratification-dependent diffusion and iron fertilization. Carbonate compensation is accounted for in these simulations. We show that neither iron fertilization nor the sinking of brines alone can account for the evolution of CO2, and that only the combination of the sinking of brines and interactive diffusion can simultaneously simulate the increase in deep Southern Ocean δ13C. The scenario that agrees best with the data takes into account all mechanisms and favours a rapid cessation of the sinking of brines around 18 000 years ago, when the Antarctic ice sheet extent was at its maximum. In this scenario, we make the hypothesis that sea ice formation was then shifted to the open ocean where the salty water is quickly mixed with fresher water, which prevents deep sinking of salty water and therefore breaks down the deep stratification and releases carbon from the abyss. Based on this scenario, it is possible to simulate both the amplitude and timing of the long-term CO2 increase during the last termination in agreement with ice core data. The atmospheric δ13C appears to be highly sensitive to changes in the terrestrial biosphere, underlining the need to better constrain the vegetation evolution during the termination.

Improved evaporative flux partitioning and carbon flux in the land surface model JULES: impact on the simulation of land surface processes in temperate Europe

The primary role of land surface models embedded in climate models is to partition surface available energy into upwards, radiative, sensible and latent heat fluxes. Partitioning of evapotranspiration, ET, is of fundamental importance: as a major component of the total surface latent heat flux, ET affects the simulated surface water balance, and related energy balance, and consequently the feedbacks with the atmosphere. In this context it is also crucial to credibly represent the CO2 exchange between ecosystems and their environment. In this study, JULES, the land surface model used in UK weather and climate models, has been evaluated for temperate Europe. Compared to eddy covariance flux measurements, the CO2 uptake by the ecosystem is underestimated and the ET overestimated. In addition, the contribution to ET from soil and intercepted water evaporation far outweighs the contribution of plant transpiration. To alleviate these biases, adaptations have been implemented in JULES, based on key literature references. These adaptations have improved the simulation of the spatio-temporal variability of the fluxes and the accuracy of the simulated GPP and ET, including its partitioning. This resulted in a shift of the seasonal soil moisture cycle. These adaptations are expected to increase the fidelity of climate simulations over Europe. Finally, the extreme summer of 2003 was used as evaluation benchmark for the use of the model in climate change studies. The improved model captures the impact of the 2003 drought on the carbon assimilation and the water use efficiency of the plants. It, however, underestimates the 2003 GPP anomalies. The simulations showed that a reduction of evaporation from the interception and soil reservoirs, albeit not of transpiration, largely explained the good correlation between the carbon and the water fluxes anomalies that was observed during 2003. This demonstrates the importance of being able to discriminate the response of individual component of the ET flux to environmental forcing.

Highly contrasting effects of different climate forcing agents on terrestrial ecosystem services

Many atmospheric constituents besides carbon dioxide (CO2) contribute to global warming, and it is common to compare their influence on climate in terms of radiative forcing, which measures their impact on the planetary energy budget. A number of recent studies have shown that many radiatively active constituents also have important impacts on the physiological functioning of ecosystems, and thus the ‘ecosystem services’ that humankind relies upon. CO2 increases have most probably increased river runoff and had generally positive impacts on plant growth where nutrients are non-limiting, whereas increases in near-surface ozone (O3) are very detrimental to plant productivity. Atmospheric aerosols increase the fraction of surface diffuse light, which is beneficial for plant growth. To illustrate these differences, we present the impact on net primary productivity and runoff of higher CO2, higher near-surface O3, and lower sulphate aerosols, and for equivalent changes in radiative forcing.We compare this with the impact of climate change alone, arising, for example, from a physiologically inactive gas such as methane (CH4). For equivalent levels of change in radiative forcing, we show that the combined climate and physiological impacts of these individual agents vary markedly and in some cases actually differ in sign. This study highlights the need to develop more informative metrics of the impact of changing atmospheric constituents that go beyond simple radiative forcing.

Black carbon vertical profiles strongly affect its radiative forcing uncertainty

The impact of black carbon (BC) aerosols on the global radiation balance is not well constrained. Here twelve global aerosol models are used to show that at least 20% of the present uncertainty in modeled BC direct radiative forcing (RF) is due to diversity in the simulated vertical profile of BC mass. Results are from phases 1 and 2 of the global aerosol model intercomparison project (AeroCom). Additionally, a significant fraction of the variability is shown to come from high altitudes, as, globally, more than 40% of the total BC RF is exerted above 5 km. BC emission regions and areas with transported BC are found to have differing characteristics. These insights into the importance of the vertical profile of BC lead us to suggest that observational studies are needed to better characterize the global distribution of BC, including in the upper troposphere.

Deglacial changes of the southern margin of the southern westerly winds revealed by terrestrial records from SW Patagonia (52 degrees S)

Much of the ongoing discussion regarding synchrony or bipolar asynchrony of paleoclimate events has centered on the timing and structure of the last glacial termination in the southern mid- latitudes, in particular the southwestern Patagonian region (50�e55�S). Its location adjacent to the Drake Passage andnear the southern margin of the southern westerly winds (SWW) allows examining the postulated links between the Southern Oceane SWW coupled system and tmospheric CO2 variations through the last glacial termination. Results from two sites located in the Última Esperanza area (52�S) allow us to infer SWW-driven changes in hydrologic balance during this critical time interval. These findings indicate peatland development under temperate/wet conditions between 14,600 and 14,900 cal yr BP, followed by cooling and a lake transgressive phase that led to a shallow lake during the early part of the Antarctic Cold Reversal (ACR, 13,600-14,600 cal yr BP), followed in turn by a deeper lake and modest warming during Younger Dryas time (YD, 11,800-13,000 cal yr BP), superseded by terrestrialization and forest expansion at the beginning of the Holocene. We propose that the SWW (i) strengthened and shifted northward during ACR time causing a precipitation rise in northwestern and southwestern Patagonia coeval with mid- and high-latitude cooling and a halt in the deglacial atmospheric CO2 rise; (ii) shifted southward during YD time causing a precipitation decline/increase in NW/SW Patagonia, respectively, high-latitude warming, and invigorated CO2 release from the Southern Ocean; (iii) became weaker between 10,000 and 11,500 cal yr BP causing a precipitation decline throughout Patagonia, concurrent with peak mid- and high-latitude temperatures and atmospheric CO2 concentrations.

Serum carbon and nitrogen stable isotopes as potential biomarkers of dietary intake and their relation with incident type 2 diabetes: the EPIC-Norfolk study

Background: Stable-isotope ratios of carbon (13C/12C, expressed as δ13C) and nitrogen (15N/14N, or δ15N) have been proposed as potential nutritional biomarkers to distinguish between meat, fish, and plant-based foods. Objective: The objective was to investigate dietary correlates of δ13C and δ15N and examine the association of these biomarkers with incident type 2 diabetes in a prospective study. Design: Serum δ13C and δ15N (‰) were measured by using isotope ratio mass spectrometry in a case-cohort study (n = 476 diabetes cases; n = 718 subcohort) nested within the European Prospective Investigation into Cancer and Nutrition (EPIC)–Norfolk population-based cohort. We examined dietary (food-frequency questionnaire) correlates of δ13C and δ15N in the subcohort. HRs and 95% CIs were estimated by using Prentice-weighted Cox regression. Results: Mean (±SD) δ13C and δ15N were −22.8 ± 0.4‰ and 10.2 ± 0.4‰, respectively, and δ13C (r = 0.22) and δ15N (r = 0.20) were positively correlated (P < 0.001) with fish protein intake. Animal protein was not correlated with δ13C but was significantly correlated with δ15N (dairy protein: r = 0.11; meat protein: r = 0.09; terrestrial animal protein: r = 0.12, P ≤ 0.013). δ13C was inversely associated with diabetes in adjusted analyses (HR per tertile: 0.74; 95% CI: 0.65, 0.83; P-trend < 0.001], whereas δ15N was positively associated (HR: 1.23; 95% CI: 1.09, 1.38; P-trend = 0.001). Conclusions: The isotope ratios δ13C and δ15N may both serve as potential biomarkers of fish protein intake, whereas only δ15N may reflect broader animal-source protein intake in a European population. The inverse association of δ13C but a positive association of δ15N with incident diabetes should be interpreted in the light of knowledge of dietary intake and may assist in identifying dietary components that are associated with health risks and benefits.