Holocene carbon cycle dynamics, links to model files

We are investigating the late Holocene rise in CO2 by performing four experiments with the climate-carbon-cycle model CLIMBER2-LPJ. Apart from the deep sea sediments, important carbon cycle processes considered are carbon uptake or release by the vegetation, carbon uptake by peatlands, and CO 2 release due to shallow water sedimentation of CaCO3. Ice core data of atmospheric CO2 between 8 ka BP and preindustrial climate can only be reproduced if CO2 outgassing due to shallow water sedimentation of CaCO3 is considered. In this case the model displays an increase of nearly 20 ppmv CO2 between 8 ka BP and present day. Model configurations that do not contain this forcing show a slight decrease in atmospheric CO2. We can therefore explain the late Holocene rise in CO2 by invoking natural forcing factors only, and anthropogenic forcing is not required to understand preindustrial CO2 dynamics.

Marine Isotope Stage (MIS) 11 results of model data with Community Climate System Model version 3 including the Dynamic Global Vegetation Model (CCSM3-DGVM) in NetCDF format at global and regional scale

The climate of Marine Isotope Stage (MIS) 11, the interglacial roughly 400,000 years ago, is investigated for four time slices, 416, 410, 400, and 394 ka. The overall picture is that MIS 11 was a relatively warm interglacial in comparison to preindustrial, with Northern Hemisphere (NH) summer temperatures early in MIS 11 (416-410 ka) warmer than preindustrial, though winters were cooler. Later in MIS 11, especially around 400 ka, conditions were cooler in the NH summer, mainly in the high latitudes. Climate changes simulated by the models were mainly driven by insolation changes, with the exception of two local feedbacks that amplify climate changes. Here, the NH high latitudes, where reductions in sea ice cover lead to a winter warming early in MIS 11, as well as the tropics, where monsoon changes lead to stronger climate variations than one would expect on the basis of latitudinal mean insolation change alone, are especially prominent. The results support a northward expansion of trees at the expense of grasses in the high northern latitudes early during MIS 11, especially in northern Asia and North America.

Results of modelled Arctic river discharge trends since 7 ka BP

The Arctic hydrological cycle throughout the Holocene is analyzed based on the results of transient simulations with the coupled atmosphere-ocean circulation model ECHO-G. The results suggest a ~ 2 % increase of mid-Holocene to preindustrial Arctic river discharges for the Eurasian continent. However, rivers of the North America Arctic realm show a moderate runoff decline of approximately 4 to 5 % for the same period. The total river discharge into the Arctic Ocean has remained at an approximately constant preindustrial level since the mid Holocene. The positive discharge trend within Eurasia is caused by a more rapid decrease in local net evaporation compared to a smaller decline in advected moisture and hence precipitation. This effect is neither recognized within the North American Arctic domain nor in the far eastern part of the Eurasian Arctic realm. A detailed comparison of these model findings with a variety of proxy studies is conducted. The collected proxy records show trends of continental surface temperatures and precipitation rates that are consistent with the simulations. A continuation of the transient Holocene runs for the 19th and 20th century with increased greenhouse gases indicates an increase of the total river influx into the Arctic Ocean of up to 7.6 %. The Eurasian river discharges increase by 7.5 %, the North American discharges by up to 8.4 %. The most rapid increases have been detected since the beginning of the 20th century. These results are corroborated by the observed rising of Arctic river discharges during the last century which is attributed to anthropogenic warming. The acceleration of the Arctic hydrological cycle in the 20th century is without precedence in the Holocene.

Geochemistry on Holocene sediments from the Mediterranean

The enhanced accumulation of organic matter in Eastern Mediterranean sapropels and their unusually low d15N values have been attributed to either enhanced nutrient availability which led to elevated primary production and carbon sequestration or to enhanced organic matter preservation under anoxic conditions. In order to evaluate these two hypothesis we have determined Ba/Al ratios, amino acid composition, N and organic C concentrations and d15N in sinking particles, surface sediments, eight spatially distributed core records of the youngest sapropel S1 (10-6 ka) and older sapropels (S5, S6) from two locations. These data suggest that (i) temporal and spatial variations in d15N of sedimentary N are driven by different degrees of diagenesis at different sites rather than by changes in N-sources or primary productivity and (ii) present day TOC export production would suffice to create a sapropel like S1 under conditions of deep-water anoxia. This implies that both enhanced TOC accumulation and d15N depletion in sapropels were due to the absence of oxygen in deep waters. Thus preservation plays a major role for the accumulation of organic-rich sediments casting doubt on the need of enhanced primary production for sapropel formation.

Palynology and radiocarbon dating on sediment profile PG1756 from Lake Billyakh, NE Siberia, Russia

In this study a radiocarbon-dated pollen record from Lake Billyakh (65°17'N, 126°47'E; 340 m a.s.l.) in the Verkhoyansk Mountains was used to reconstruct vegetation and climate change since about 15 kyr BP (1 kyr=1000 cal. yr). The pollen record and pollen-based biome reconstruction suggest that open cool steppe and grass and sedge tundra communities with Poaceae, Cyperaceae, Artemisia, Chenopodiaceae, Caryophyllaceae and Selaginella rupestris dominated the area from 15 to 13.5 kyr BP. On the other hand, the constant presence of Larix pollen in quantities comparable to today's values points to the constant presence of boreal deciduous conifer trees in the regional vegetation during the last glaciation. A major spread of shrub tundra communities, including birch (Betula sect. Nanae), alder (Duschekia fruticosa) and willow (Salix) species, is dated to 13.5-12.7 kyr BP, indicating a noticeable increase in precipitation toward the end of the last glaciation, particularly during the Allerød Interstadial. Between 12.7 and 11.4 kyr BP pollen percentages of herbaceous taxa rapidly increased, whereas shrub taxa percentages decreased, suggesting strengthening of the steppe communities associated with the relatively cold and dry Younger Dryas Stadial. However, the pollen data in hand indicate that Younger Dryas climate was less severe than the climate during the earlier interval from 15 to 13.5 kyr BP. The onset of the Holocene is marked in the pollen record by the highest values of shrub and lowest values of herbaceous taxa, suggesting a return of warmer and wetter conditions after 11.4 kyr BP. Percentages of tree taxa increase gradually and reach maximum values after 7 kyr BP, reflecting the spread of boreal cold deciduous and taiga forests in the region. An interval between 7 and 2 kyr BP is noticeable for the highest percentages of Scots spine (Pinus subgen. Diploxylon), spruce (Picea) and fir (Abies) pollen, indicating mid-Holocene spread of boreal forest communities in response to climate amelioration and degradation of the permafrost layer.

Alkenone-derived sea surface temperature records of three sediment cores from the southern Benguela upwelling system

We present alkenone-derived Sea Surface Temperature (SST) records from three marine cores collected within the southern Benguela Upwelling System (BUS) spanning the last 3 ka. The SST evolution over the last 3 millennia is marked by aperiodic millennial-scale oscillations that broadly correspond to climatic anomalies identified over the North Atlantic region. The BUS SST data further suggest cooling and warming trends opposite to the temperature evolution in the Moroccan upwelling region and in Antarctica. In contrast, the last decades are marked by a cooling of unprecedented magnitude in both the Benguela and Moroccan upwelling systems, which is not observed in the Antarctic record. These contrasted responses in Atlantic upwelling systems triggered by natural and anthropogenic forcings shed light on how different climatic mechanisms are mediated by ocean-atmosphere interactions and transmitted to the geological records of past and present climate changes.

Sea surface temperatures during Marine Isotope Stage (MIS) 11

The Marine Isotope Stage (MIS) 11 (424-374 ka) was characterized by a protracted deglaciation and an unusually long climatic optimum. It remains unclear to what degree the climate development during this interglacial reflects the unusually weak orbital forcing or greenhouse gas trends. Previously, arguments about the duration and timing of the MIS11 climatic optimum and about the pace of the deglacial warming were based on a small number of key records, which appear to show regional differences. In order to obtain a global signal of climate evolution during MIS11, we compiled a database of 78 sea surface temperature (SST) records from 57 sites spanning MIS11, aligned these individually on the basis of benthic (N = 28) or planktonic (N = 31) stable oxygen isotope curves to a common time frame and subjected 48 of them to an empirical orthogonal function (EOF) analysis. The analysis revealed a high commonality among all records, with the principal SST trend explaining almost 49% of the variability. This trend indicates that on the global scale, the surface ocean underwent rapid deglacial warming during Termination V, in pace with carbon dioxide rise, followed by a broad SST optimum centered at ~410 kyr. The second EOF, which explained ~18% of the variability, revealed the existence of a different SST trend, characterized by a delayed onset of the temperature optimum during MIS11 at ~398 kyr, followed by a prolonged warm period lasting beyond 380 kyr. This trend is most consistently manifested in the mid-latitude North Atlantic and Mediterranean Sea and is here attributed to the strength of the Atlantic meridional overturning circulation. A sensitivity analysis indicates that these results are robust to record selection and to age-model uncertainties of up to 3-6 kyr, but more sensitive to SST seasonal attribution and SST uncertainties >1 °C. In order to validate the CCSM3 (Community Climate System Model, version 3) predictive potential, the annual and seasonal SST anomalies recorded in a total of 74 proxy records were compared with runs for three time slices representing orbital configuration extremes during the peak interglacial of MIS11. The modeled SST anomalies are characterized by a significantly lower variance compared to the reconstructions. Nevertheless, significant correlations between proxy and model data are found in comparisons on the seasonal basis, indicating that the model captures part of the long-term variability induced by astronomical forcing, which appears to have left a detectable signature in SST trends.

A reconstruction of atmospheric carbon dioxide and its stable carbon isotopic composition from the penultimate glacial maximum to the glacial inception

The reconstruction of the stable carbon isotope evolution in atmospheric CO2 (d13Catm ), as archived in Antarctic ice cores, bears the potential to disentangle the contributions of the different carbon cycle fluxes causing past CO2 variations. Here we present a new record of d13Catm before, during and after the Marine Isotope Stage 5.5 (155 000 to 105 000 years BP). The record was derived with a well established sublimation method using ice from the EPICA Dome C (EDC) and the Talos Dome ice cores in East Antarctica. We find a 0.4 permil shift to heavier values between the mean d13Catm level in the Penultimate (~ 140 000 years BP) and Last Glacial Maximum (~ 22 000 years BP), which can be explained by either (i) changes in the isotopic composition or (ii) intensity of the carbon input fluxes to the combined ocean/atmosphere carbon reservoir or (iii) by long-term peat buildup. Our isotopic data suggest that the carbon cycle evolution along Termination II and the subsequent interglacial was controlled by essentially the same processes as during the last 24 000 years, but with different phasing and magnitudes. Furthermore, a 5000 years lag in the CO2 decline relative to EDC temperatures is confirmed during the glacial inception at the end of MIS 5.5 (120 000 years BP). Based on our isotopic data this lag can be explained by terrestrial carbon release and carbonate compensation.

Age determination, planktic foraminifera and stable istopes of sediment core MSM05/5_712-2

A sediment core from the West Spitsbergen continental margin was studied to reconstruct climate and paleoceanographic variability during the last ~9 ka in the eastern Fram Strait. Our multiproxy evidence suggests that the establishment of the modern oceanographic configuration in the eastern Fram Strait occurred stepwise, in response to the postglacial sea-level rise and the related onset of modern sea-ice production on the shallow Siberian shelves. The late Early and Mid Holocene interval (9 to 5 ka) was generally characterized by relatively unstable conditions. High abundance of the subpolar planktic foraminifer species Turborotalita quinqueloba implies strong intensity of Atlantic Water (AW) inflow with high productivity and/or high AW temperatures, resulting in a strong heat flux to the Arctic. A series of short-lived cooling events (8.2, 6.9. and 6.1 ka) occurred superimposed on the warm late Early and Mid Holocene conditions. Our proxy data imply that simultaneous to the complete postglacial flooding of Arctic shallow shelves and the initiation of modern sea-ice production, strong advance of polar waters initiated modern oceanographic conditions in the eastern Fram Strait at ~5.2 ka. The Late Holocene was marked by the dominance of the polar planktic foraminifer species Neogloboquadrina pachyderma, a significant expansion of sea ice/icebergs, and strong stratification of the water column. Although planktic foraminiferal assemblages as well as sea surface and subsurface temperatures suggest a return of slightly strengthened advection of subsurface Atlantic Water after 3 ka, a relatively stable cold-water layer prevailed at the sea surface and the study site was probably located within the seasonally fluctuating marginal ice zone during the Neoglacial period.