Age model of firn core DML05C98_32 (B32)

In the framework of the European Project for Ice Coring in Antarctica (EPICA), a comprehensive glaciological pre-site survey has been carried out on Amundsenisen, Dronning Maud Land, East Antarctica, in the past decade. Within this survey, four intermediate-depth ice cores and 13 snow pits were analyzed for their ionic composition and interpreted with respect to the spatial and temporal variability of volcanic sulphate deposition. The comparison of the non-sea-salt (nss)-sulphate peaks that are related to the well-known eruptions of Pinatubo and Cerro Hudson in AD 1991 revealed sulphate depositions of comparable size (15.8 ± 3.4 kg/km**2) in 11 snow pits. There is a tendency to higher annual concentrations for smaller snow-accumulation rates. The combination of seasonal sodium and annually resolved nss-sulphate records allowed the establishment of a time-scale derived by annual-layer counting over the last 2000 years and thus a detailed chronology of annual volcanic sulphate deposition. Using a robust outlier detection algorithm, 49 volcanic eruptions were identified between AD 165 and 1997. The dating uncertainty is ±3 years between AD 1997 and 1601, around ±5 years between AD 1601 and 1257, and increasing to ±24 years at AD 165, improving the accuracy of the volcanic chronology during the penultimate millennium considerably.

Simulated Bromoform emission and concentration using the ocean biogeochemistry model MPIOM/HAMOCC forced by 6-hourly NCEP data

Bromoform (CHBr3) is one important precursor of atmospheric reactive bromine species that are involved in ozone depletion in the troposphere and stratosphere. In the open ocean bromoform production is linked to phytoplankton that contains the enzyme bromoperoxidase. Coastal sources of bromoform are higher than open ocean sources. However, open ocean emissions are important because the transfer of tracers into higher altitude in the air, i.e. into the ozone layer, strongly depends on the location of emissions. For example, emissions in the tropics are more rapidly transported into the upper atmosphere than emissions from higher latitudes. Global spatio-temporal features of bromoform emissions are poorly constrained. Here, a global three-dimensional ocean biogeochemistry model (MPIOM-HAMOCC) is used to simulate bromoform cycling in the ocean and emissions into the atmosphere using recently published data of global atmospheric concentrations (Ziska et al., 2013) as upper boundary conditions. Our simulated surface concentrations of CHBr3 match the observations well. Simulated global annual emissions based on monthly mean model output are lower than previous estimates, including the estimate by Ziska et al. (2013), because the gas exchange reverses when less bromoform is produced in non-blooming seasons. This is the case for higher latitudes, i.e. the polar regions and northern North Atlantic. Further model experiments show that future model studies may need to distinguish different bromoform-producing phytoplankton species and reveal that the transport of CHBr3 from the coast considerably alters open ocean bromoform concentrations, in particular in the northern sub-polar and polar regions.

R-script modified VS-Lite model

Drought is a key factor affecting forest ecosystem processes at different spatio-temporal scales. For accurately modeling tree functioning ? and thus for producing reliable simulations of forest dynamics ? the consideration of the variability in the timing and extent of drought effects on tree growth is essential, particularly in strongly seasonal climates such as in the Mediterranean area. Yet, most dynamic vegetation models (DVMs) do not include this intra-annual variability of drought effects on tree growth. We present a novel approach for linking tree-ring data to drought simulations in DVMs. A modified forward model of tree-ring width (VS-Lite) was used to estimate seasonal- and site-specific growth responses to drought of Scots pine (Pinus sylvestris L.), which were subsequently implemented in the DVM ForClim. Ring-width data from sixteen sites along a moisture gradient from Central Spain to the Swiss Alps, including the dry inner Alpine valleys, were used to calibrate the forward ring-width model, and inventory data from managed Scots pine stands were used to evaluate ForClim performance. The modified VS-Lite accurately estimated the year-to-year variability in ring-width indices and produced realistic intra-annual growth responses to soil drought, showing a stronger relationship between growth and drought in spring than in the other seasons and thus capturing the strategy of Scots pine to cope with drought. The ForClim version including seasonal variability in growth responses to drought showed improved predictions of stand basal area and stem number, indicating the need to consider intra-annual differences in climate-growth relationships in DVMs when simulating forest dynamics. Forward modeling of ring-width growth may be a powerful tool to calibrate growth functions in DVMs that aim to simulate forest properties in across multiple environments at large spatial scales.

Snow-pit measurements on Atka Bay landfast sea ice during ANT-Land 2012/2013 season

Up to now, snow cover on Antarctic sea ice and its impact on radar backscatter, particularly after the onset of freeze/thaw processes, are not well understood. Here we present a combined analysis of in situ observations of snow properties from the landfast sea ice in Atka Bay, Antarctica, and high-resolution TerraSAR-X backscatter data, for the transition from austral spring (November 2012) to summer (January 2013). The physical changes in the seasonal snow cover during that time are reflected in the evolution of TerraSAR-X backscatter. We are able to explain 76-93% of the spatio-temporal variability of the TerraSAR-X backscatter signal with up to four snowpack parameters with a root-mean-squared error of 0.87-1.62 dB, using a simple multiple linear model. Over the complete study, and especially after the onset of early-melt processes and freeze/thaw cycles, the majority of variability in the backscatter is influenced by changes in snow/ice interface temperature, snow depth and top-layer grain size. This suggests it may be possible to retrieve snow physical properties over Antarctic sea ice from X-band SAR backscatter.

Snow depths of volcanic events and accumulation rates of firn cores from the IPY traverse, Antarctica

Volcanic signatures in ice-core records provide an excellent means to date the cores and obtain information about accumulation rates. From several ice cores it is thus possible to extract a spatio-temporal accumulation pattern. We show records of electrical conductivity and sulfur from 13 firn cores from the Norwegian-USA scientific traverse during the International Polar Year 2007-2009 (IPY) through East Antarctica. Major volcanic eruptions are identified and used to assess century-scale accumulation changes. The largest changes seem to occur in the most recent decades with accumulation over the period 1963-2007/08 being up to 25% different from the long-term record. There is no clear overall trend, some sites show an increase in accumulation over the period 1963 to present while others show a decrease. Almost all of the sites above 3200 m above sea level (asl) suggest a decrease. These sites also show a significantly lower accumulation value than large-scale assessments both for the period 1963 to present and for the long-term mean at the respective drill sites. The spatial accumulation distribution is influenced mainly by elevation and distance to the ocean (continentality), as expected. Ground-penetrating radar data around the drill sites show a spatial variability within 10-20% over several tens of kilometers, indicating that our drill sites are well representative for the area around them. Our results are important for large-scale assessments of Antarctic mass balance and model validation.

Last Interglacial synthesis of high-latitude temperature: temperature anomalies and associated errors for 4 time slices

The Last Interglacial (LIG, 129-116 thousand of years BP, ka) represents a test bed for climate model feedbacks in warmer-than-present high latitude regions. However, mainly because aligning different palaeoclimatic archives and from different parts of the world is not trivial, a spatio-temporal picture of LIG temperature changes is difficult to obtain. Here, we have selected 47 polar ice core and sub-polar marine sediment records and developed a strategy to align them onto the recent AICC2012 ice core chronology. We provide the first compilation of high-latitude temperature changes across the LIG associated with a coherent temporal framework built between ice core and marine sediment records. Our new data synthesis highlights non-synchronous maximum temperature changes between the two hemispheres with the Southern Ocean and Antarctica records showing an early warming compared to North Atlantic records. We also observe warmer than present-day conditions that occur for a longer time period in southern high latitudes than in northern high latitudes. Finally, the amplitude of temperature changes at high northern latitudes is larger compared to high southern latitude temperature changes recorded at the onset and the demise of the LIG. We have also compiled four data-based time slices with temperature anomalies (compared to present-day conditions) at 115 ka, 120 ka, 125 ka and 130 ka and quantitatively estimated temperature uncertainties that include relative dating errors. This provides an improved benchmark for performing more robust model-data comparison. The surface temperature simulated by two General Circulation Models (CCSM3 and HadCM3) for 130 ka and 125 ka is compared to the corresponding time slice data synthesis. This comparison shows that the models predict warmer than present conditions earlier than documented in the North Atlantic, while neither model is able to produce the reconstructed early Southern Ocean and Antarctic warming. Our results highlight the importance of producing a sequence of time slices rather than one single time slice averaging the LIG climate conditions.

Age determination of seven sediment cores from the South Pacific Ocean

During the last deglaciation, the opposing patterns of atmospheric CO2 and radiocarbon activities (D14C) suggest the release of 14C-depleted CO2 from old carbon reservoirs. Although evidences point to the deep Pacific as a major reservoir of this 14C-depleted carbon, its extent and evolution still need to be constrained. Here we use sediment cores retrieved along a South Pacific transect to reconstruct the spatio-temporal evolution of D14C over the last 30,000 years. In ~2,500-3,600 m water depth, we find 14C-depleted deep waters with a maximum glacial offset to atmospheric 14C (DD14C = -1,000 per mil). Using a box model, we test the hypothesis that these low values might have been caused by an interaction of aging and hydrothermal CO2 influx. We observe a rejuvenation of circumpolar deep waters synchronous and potentially contributing to the initial deglacial rise in atmospheric CO2. These findings constrain parts of the glacial carbon pool to the deep South Pacific.

Lake Son Kol geochemistry, µXRF, biogeochemistry and pollen record

In general, a moderate drying trend is observed in mid-latitude arid Central Asia since the Mid-Holocene, attributed to the progressively weakening influence of the mid-latitude Westerlies on regional climate. However, as the spatio-temporal pattern of this development and the underlying climatic mechanisms are yet not fully understood, new high-resolution paleoclimate records from this region are needed. Within this study, a sediment core from Lake Son Kol (Central Kyrgyzstan) was investigated using sedimentological, (bio)geochemical, isotopic, and palynological analyses, aiming at reconstructing regional climate development during the last 6000 years. Biogeochemical data, mainly reflecting summer moisture conditions, indicate predominantly wet conditions until 4950 cal. yr BP, succeeded by a pronounced dry interval between 4950 and 3900 cal. yr BP. In the following, a return to wet conditions and a subsequent moderate drying trend until present times are observed. This is consistent with other regional paleoclimate records and likely reflects the gradual Late Holocene diminishment of the amount of summer moisture provided by the mid-latitude Westerlies. However, climate impact of the Westerlies was apparently not only restricted to the summer season but also significant during winter as indicated by recurrent episodes of enhanced allochthonous input through snowmelt, occurring before 6000 cal. yr BP and at 5100-4350, 3450-2850, and 1900-1500 cal. yr BP. The distinct ~1500-year periodicity of these episodes of increased winter precipitation in Central Kyrgyzstan resembles similar cyclicities observed in paleoclimate records around the North Atlantic, likely indicating a hemispheric-scale climatic teleconnection and an impact of North Atlantic Oscillation (NAO) variability in Central Asia.

Chert intervals in DSDP and ODP sites (Table 1)

Radiolarian cherts in the Tethyan realm of Jurassic age were recently interpreted as resulting from high biosiliceous productivity along upwelling zones in subequatorial paleolatitudes the locations of which were confirmed by revised paleomagnetic estimates. However, the widespread occurrence of cherts in the Eocene suggests that cherts may not always be reliable proxies of latitude and upwelling zones. In a new survey of the global spatio-temporal distribution of Cenozoic cherts in Deep Sea Drilling Project (DSDP) and Ocean Drilling Program (ODP) sediment cores, we found that cherts occur most frequently in the Paleocene and early Eocene, with a peak in occurrences at ~50 Ma that is coincident with the time of highest bottom water temperatures of the early Eocene climatic optimum (EECO) when the global ocean was presumably characterized by reduced upwelling efficiency and biosiliceous productivity. Cherts occur less commonly during the subsequent Eocene global cooling trend. Primary paleoclimatic factors rather than secondary diagenetic processes seem therefore to control chert formation. This timing of peak Eocene chert occurrence, which is supported by detailed stratigraphic correlations, contradicts currently accepted models that involve an initial loading of large amounts of dissolved silica from enhanced weathering and/or volcanism in a supposedly sluggish ocean of the EECO, followed during the subsequent middle Eocene global cooling by more vigorous oceanic circulation and consequent upwelling that made this silica reservoir available for enhanced biosilicification, with the formation of chert as a result of biosilica transformation during diagenesis. Instead, we suggest that basin-basin fractionation by deep-sea circulation could have raised the concentration of EECO dissolved silica especially in the North Atlantic, where an alternative mode of silica burial involving widespread direct precipitation and/or absorption of silica by clay minerals could have been operative in order to maintain balance between silica input and output during the upwelling-deficient conditions of the EECO. Cherts may therefore not always be proxies of biosiliceous productivity associated with latitudinally focused upwelling zones.

The distributions, C- and O-isotopic compositions, and frequency of giant-ooids

Early Triassic oceans were characterized by deposition of a number of "anachronistic facies", including microbialites, seafloor carbonate cement fans, and giant ooids. Giant ooids were particularly prevalent in Lower Triassic sections across South China and exhibit unusual features that may provide insights into marine environmental conditions following the end-Permian mass extinction. The section at Moyang (Guizhou Province) contains abundant giant ooids ranging in size between 2 and 6 mm (maximum 12 mm) and exhibiting various cortical structures, including regular, deformed, compound, regenerated and "domed". Preservation of ooid cortical structure is generally good as indicated by petrographic observations, and trace element and carbon isotope analyses suggest that diagenesis occurred in a closed diagenetic system. All ooids exhibit fine concentric laminae, frequently alternating between light-colored coarsely crystalline and dark-colored finely crystalline layers probably reflecting variation in organic content or original mineralogy. Under scanning electron microscope, biomineralized filaments or biofilms and tiny carbonate fluorapatite (CFA) crystals are commonly found in the finely crystalline layers. We infer that the precipitation of CFA was related to adsorption of P via microbial activity on the surfaces of ooids following episodic incursions of deep waters rich in carbon dioxide, hydrogen sulfide and phosphate into shallow-marine environments. Giant ooid precipitation may have been promoted in shallow ramp settings during these events by increased watermass agitation and supersaturation with respect to calcium carbonate, as well as reduced carbonate removal rates through biotic skeletal formation. Spatio-temporal distribution data reveal that giant ooids were widespread in the Tethyan region during the Early Triassic, and that they were most abundant immediately after the end-Permian crisis and disappeared gradually as metazoans repopulated marine environments.