Crustal uplift rates of 4 GPS stations near Jakobshavn Isbrae (2006-2009)

We analyze 2006-2009 data from four continuous Global Positioning System (GPS) receivers located between 5 and 150 km from the glacier Jakobshavn Isbrae, West Greenland. The GPS stations were established on bedrock to determine the vertical crustal motion due to the unloading of ice from Jakobshavn Isbrae. All stations experienced uplift, but the uplift rate at Kangia North, only 5 km from the glacier front, was about 10 mm/yr larger than the rate at Ilulissat, located only ~45 km further away. This suggests that most of the uplift is due to the unloading of the Earth's surface as Jakobshavn thins and loses mass. Our estimate of Jakobshavn's contribution to uplift rates at Kangia North and Ilulissat are 14.6 ± 1.7 mm/yr and 4.9 ± 1.1 mm/yr, respectively. The observed rates are consistent with a glacier thinning model based on repeat altimeter surveys from NASA's Airborne Topographic Mapper (ATM), which shows that Jakobshavn lost mass at an average rate of 22 ± 2 km**3/yr between 2006 and 2009. At Kangia North and Ilulissat, the predicted uplift rates computed using thinning estimates from the ATM laser altimetry are 12.1 ± 0.9 mm/yr and 3.2 ± 0.3 mm/yr, respectively. The observed rates are slightly larger than the predicted rates. The fact that the GPS uplift rates are much larger closer to Jakobshavn than further away, and are consistent with rates inferred using the ATM-based glacier thinning model, shows that GPS measurements of crustal motion are a potentially useful method for assessing ice-mass change models.

Holocene precipitation change in different monsoon sub-regions (time-slices and transient data) simulated by different global climate models, with links to model results

The recently proposed global monsoon hypothesis interprets monsoon systems as part of one global-scale atmospheric overturning circulation, implying a connection between the regional monsoon systems and an in-phase behaviour of all northern hemispheric monsoons on annual timescales (Trenberth et al., 2000). Whether this concept can be applied to past climates and variability on longer timescales is still under debate, because the monsoon systems exhibit different regional characteristics such as different seasonality (i.e. onset, peak, and withdrawal). To investigate the interconnection of different monsoon systems during the pre-industrial Holocene, five transient global climate model simulations have been analysed with respect to the rainfall trend and variability in different sub-domains of the Afro-Asian monsoon region. Our analysis suggests that on millennial timescales with varying orbital forcing, the monsoons do not behave as a tightly connected global system. According to the models, the Indian and North African monsoons are coupled, showing similar rainfall trend and moderate correlation in rainfall variability in all models. The East Asian monsoon changes independently during the Holocene. The dissimilarities in the seasonality of the monsoon sub-systems lead to a stronger response of the North African and Indian monsoon systems to the Holocene insolation forcing than of the East Asian monsoon and affect the seasonal distribution of Holocene rainfall variations. Within the Indian and North African monsoon domain, precipitation solely changes during the summer months, showing a decreasing Holocene precipitation trend. In the East Asian monsoon region, the precipitation signal is determined by an increasing precipitation trend during spring and a decreasing precipitation change during summer, partly balancing each other. A synthesis of reconstructions and the model results do not reveal an impact of the different seasonality on the timing of the Holocene rainfall optimum in the different sub-monsoon systems. They rather indicate locally inhomogeneous rainfall changes and show, that single palaeo-records should not be used to characterise the rainfall change and monsoon evolution for entire monsoon sub-systems.

Age models of sediment cores from the Southern Ocean

Recent geochemical models invoke ocean alkalinity changes, particularly in the surface Southern Ocean, to explain glacial age pCO2 reduction. In such models, alkalinity increases in glacial periods are driven by reductions in North Atlantic Deep Water (NADW) supply, which lead to increases in deep-water nutrients and dissolution of carbonate sediments, and to increased alkalinity of Circumpolar Deep Water upwelling in the surface Southern Ocean. We use cores from the Southeast Indian Ridge and from the deep Cape Basin in the South Atlantic to show that carbonate dissolution was enhanced during glacial stages in areas now bathed by Circumpolar Deep Water. This suggests that deep Southern Ocean carbonate ion concentrations were lower in glacial stages than in interglacials, rather than higher as suggested by the polar alkalinity model [Broecker and Peng, 1989, doi:10.1029/GB001i001p00015]. Our observations show that changes in Southern Ocean CaCO3 preservation are coherent with changes in the relative flux of NADW, suggesting that Southern Ocean carbonate chemistry is closely linked to changes in deepwater circulation. The pattern of enhanced dissolution in glacials is consistent with a reduction in the supply of nutrient-depleted water (NADW) to the Southern Ocean and with an increase of nutrients in deep water masses. Carbonate mass accumulation rates on the Southeast Indian Ridge (3200-3800 m), and in relatively shallow cores (<3000 m) from the Kerguelen Plateau and the South Pacific were significantly reduced during glacial stages, by about 50%. The reduced carbonate mass accumulation rates and enhanced dissolution during glacials may be partly due to decreases in CaCO3:Corg flux ratios, acting as another mechanism which would raise the alkalinity of Southern Ocean surface waters. The polar alkalinity model assumes that the ratio of organic carbon to carbonate production on surface alkalinity is constant. Even if overall productivity in the Southern Ocean were held constant, a decrease in the CaCO3:Corg ratio would result in increased alkalinity and reduced pCO2 in Southern Ocean surface waters during glacials. This ecologically driven surface alkalinity change may enhance deepwater-mediated changes in alkalinity, and amplify rapid changes in pCO2.

Age models and stabile isotope ratios from SO35 cores

Stable isotope measurements on the planktonic foraminifer Globigerinoides ruber (white) have been carried out on a number of selected deep-seas sediment cores from the South Lau and Norlh Fiji Basins. The d18O-curves show good correlation with the inter-ocean oraphic correlation composite d18O-record of the standard reference section (Prell et al. 1986), which, in combination with the chronostratigraphic classifications of Herterich & Sarnthein (1984, modified) and Imbrie et al. 1984), allows a detailed dating of the sedimentary sequences. The deepest layers in core no. 119 (southern Lau Basin) could be assigned to Isotope Stage 24. Measurements made on bulk carbonate in two cores show a much higher glacial-interglacial amplitude, allowing the general identification of the conventional oxygen isotope stages. The d13C-values of the benthic foraminifer Cibicidoides wuellerstorfi show progressively lighter values northwards reflecting an increasing contribution of the isotopically lighter CO2 from the remineralisation of organic matter during the general northward movement of the deep water masses. Cyclicities in the sedimentation rates were observed in core nos. 117 and 119 (both southern Lau Basin) where the interglacials exhibit higher levels than the glacials. Calculated new or export paleoproductivity show that the glacials had higher productivity in the euphotic zone. From the oxygen isotope stratigraphy, the five ash layers in core nos. 117 and 119 could be dated as about 530 ka B.P. in Stage 14, 695 ka B.P. in Stage 18, 775 ka B.P. in Stage 21, 790 ka B.P. and 825 ka B.P. in Stage 22. Carbonate dissolution occurred during stages 5, 8 and 10 to 12.

Age models and stable isotope ratios from the Mid-Miocene

Mid-Miocene pelagic sedimentary sections can be correlated using intermediate and high resolution oxygen and carbon isotopic records of benthic foraminifera. Precision of a few tens of thousands of years is readily achievable at sites with high sedimentation rates, for example, Deep Sea Drilling Project sites 289 and 574. The mid-Miocene carbon isotope records are characterized by an interval of high d13C values between 17 and 13.5 Ma (the Monterey Excursion of Vincent and Berger 1985) upon which are superimposed a series of periodic or quasi-periodic fluctuations in d13C values. These fluctuations have a period of approximately 440 kyr, suggestive of the 413 kyr cycle predicted by Milankovitch theory. Vincent and Berger proposed that the Monterey Excursion was the result of increased organic carbon burial in continental margins sediments. The increased d13C values (called 13C maxima) superimposed on the generally high mid-Miocene signal coincide with increases in d18O values suggesting that periods of cooling and/or ice buildup were associated with exceptionally rapid burial of organic carbon and lowered atmospheric CO2 levels. It is likely that during the Monterey Excursion the ocean/atmosphere system became progressively more sensitive to small changes in insolation, ultimately leading to major cooling of deep water and expansion of continental ice. We have assigned an absolute chronology, based on biostratigraphic and magneto-biostratigraphic datum levels, to the isotope stratigraphy and have used that chronology to correlate unconformities, seismic reflectors, carbonate minima, and dissolution intervals. Intervals of sediment containing 13C maxima are usually better preserved than the overlying and underlying sediments, indicating that the d13C values of TCO2 in deep water and the corrosiveness of seawater are inversely correlated. This again suggests that the 13C maxima were associated with rapid burial of organic carbon and reduced levels of atmospheric CO2. The absolute chronology we have assigned to the isotopic record indicates that the major mid-Miocene deepwater cooling/ice volume expansion took 2 m.y. and was not abrupt as had been reported previously. The cooling appears abrupt at many sites because the interval is characterized by a number of dissolution intervals. The cooling was not monotonic, and the 2 m.y. interval included an episode of especially rapid cooling as well as a brief return to warmer conditions before the final phase of the cooling period. The increase in d18O values of benthic foraminifera between 14.9 and 12.9 Ma was greatest at deeper water sites and at sites closest to Antarctica. The data suggest that the d18O value of seawater increased by no more than about 1.1 per mil during this interval and that the remainder of the change in benthic d18O values resulted from cooling in Antarctic regions of deepwater formation. Equatorial planktonic foraminifera from sites 237 and 289 exhibit a series of 0.4 per mil steplike increases in d13C values. Only one of these increases in planktonic d13C is correlated with any of the features in the mid-Miocene benthic carbon isotope record.

Age models of five ice cores from the Dronning Maud Land, Antarctica

For the first time, annually resolved accumulation rates have been determined in central Antarctica by means of counting seasonal signals of ammonium, calcium, and sodium. All records, obtained from three intermediate depth ice cores from Dronning Maud Land, East Antarctica, show rather constant accumulation rates throughout the last 9 centuries with mean values of 63, 61, and 44 mm H2O yr**-1 and a typical year-to-year variation of about 30%. For the last few decades, no trend was detected accounting for the high natural variability of all records. A significant weak intersite correlation is apparent only between two cores when the high-frequency part with periods less than 30 years is removed. By analyzing the records in the frequency domain, no persistent periods were found. This suggests that the snow accumulation in this area is mainly influenced by local deposition patterns and may be additionally masked by redistribution of snow due to wind. By comparing accumulation rates over the last 2 millennia a distinct change in the layer thickness in one of the three cores was found, which might be attributed either to an area upstream of the drilling site with lower accumulation rates, or to deposition processes influenced by surface undulations. The missing of a clear correlation between the accumulation rate histories at the three locations is also important for the interpretation of small, short time variations of past precipitation records obtained from deep ice cores.

Elevation Models of Antarctica and Greenland derived from CryoSat-2 in the period 2011 to 2013

This study focuses on the present-day surface elevation of the Greenland and Antarctic ice sheets. Based on 3 years of CryoSat-2 data acquisition we derived new elevation models (DEMs) as well as elevation change maps and volume change estimates for both ice sheets. Here we present the new DEMs and their corresponding error maps. The accuracy of the derived DEMs for Greenland and Antarctica is similar to those of previous DEMs obtained by satellite-based laser and radar altimeters. Comparisons with ICESat data show that 80% of the CryoSat-2 DEMs have an uncertainty of less than 3 m ± 15 m. The surface elevation change rates between January 2011 and January 2014 are presented for both ice sheets. We compared our results to elevation change rates obtained from ICESat data covering the time period from 2003 to 2009. The comparison reveals that in West Antarctica the volume loss has increased by a factor of 3. It also shows an anomalous thickening in Dronning Maud Land, East Antarctica which represents a known large-scale accumulation event. This anomaly partly compensates for the observed increased volume loss of the Antarctic Peninsula and West Antarctica. For Greenland we find a volume loss increased by a factor of 2.5 compared to the ICESat period with large negative elevation changes concentrated at the west and southeast coasts. The combined volume change of Greenland and Antarctica for the observation period is estimated to be -503 ± 107 km**3/yr. Greenland contributes nearly 75% to the total volume change with -375 ± 24 km**3/yr.

Age models for sediment cores in the eastern tropical Pacific

We present new high-resolution N isotope records from the Gulf of Tehuantepec and the Nicaragua Basin spanning the last 50-70 ka. The Tehuantepec site is situated within the core of the north subtropical denitrification zone while the Nicaragua site is at the southern boundary. The d15N record from Nicaragua shows an 'Antarctic' timing similar to denitrification changes observed off Peru-Chile but is radically different from the northern records. We attribute this to the leakage of isotopically heavy nitrate from the South Pacific oxygen minimum zone (OMZ) into the Nicaragua Basin. The Nicaragua record leads the other eastern tropical North Pacific (ETNP) records by about 1000 years because denitrification peaks in the eastern tropical South Pacific (ETSP) before denitrification starts to increase in the Northern Hemisphere OMZ, i.e., during warming episodes in Antarctica. We find that the influence of the heavy nitrate leakage from the ETSP is still noticeable, although attenuated, in the Gulf of Tehuantepec record, particularly at the end of the Heinrich events, and tends to alter the recording of millennial timescale denitrification changes in the ETNP. This implies (1) that sedimentary d15N records from the southern parts of the ETNP cannot be used straightforwardly as a proxy for local denitrification and (2) that denitrification history in the ETNP, like in the Arabian Sea, is synchronous with Greenland temperature changes. These observations reinforce the conclusion that on millennial timescales during the last ice age, denitrification in the ETNP is strongly influenced by climatic variations that originated in the high-latitude North Atlantic region, while commensurate changes in Southern Ocean hydrography more directly, and slightly earlier, affected oxygen concentrations in the ETSP. Furthermore, the d15N records imply ongoing physical communication across the equator in the shallow subsurface continuously over the last 50-70 ka.

Three species distribution models for the marine diatom Fragilariopsis kerguelensis of the Southern Ocean, links to NetCDF files

Fragilariopsis kerguelensis, a dominant diatom species throughout the Antarctic Circumpolar Current, is coined to be one of the main drivers of the biological silicate pump. Here, we study the distribution of this important species and expected consequences of climate change upon it, using correlative species distribution modeling and publicly available presence-only data. As experience with SDM is scarce for marine phytoplankton, this also serves as a pilot study for this organism group. We used the maximum entropy method to calculate distribution models for the diatom F. kerguelensis based on yearly and monthly environmental data (sea surface temperature, salinity, nitrate and silicate concentrations). Observation data were harvested from GBIF and the Global Diatom Database, and for further analyses also from the Hustedt Diatom Collection (BRM). The models were projected on current yearly and seasonal environmental data to study current distribution and its seasonality. Furthermore, we projected the seasonal model on future environmental data obtained from climate models for the year 2100. Projected on current yearly averaged environmental data, all models showed similar distribution patterns for F. kerguelensis. The monthly model showed seasonality, for example, a shift of the southern distribution boundary toward the north in the winter. Projections on future scenarios resulted in a moderately to negligibly shrinking distribution area and a change in seasonality. We found a substantial bias in the publicly available observation datasets, which could be reduced by additional observation records we obtained from the Hustedt Diatom Collection. Present-day distribution patterns inferred from the models coincided well with background knowledge and previous reports about F. kerguelensis distribution, showing that maximum entropy-based distribution models are suitable to map distribution patterns for oceanic planktonic organisms. Our scenario projections indicate moderate effects of climate change upon the biogeography of F. kerguelensis.

The Gepatschferner from 1850 - 2006, with links to digital elevation models and glacier margins

This theses investigates changes at Gepatschferner in length, area and volume since the last glacier maximum in 1850. Changes are discussed for the following time periods: 1850-1922, 1922-1971, 1971-1997, 1997-2006. Digital elevation models were created for 1850 from geomorphological data and for 1922 and 1971 from historical maps. Existing DEMs for 1997 and 2006 were further analysed. Since 1850 Gepatschferner has retreated by 2 km in length and has lost 32% of its area and 36% of its volume. The rate of loss of volume is increasing faster than the rate of loss of area and losses in the upper regions of the glacier are becoming increasingly more important to overall losses. The largest losses per 50 m elevation increment occur at the tongue. These losses are greatest in the most recent time step studied, 1997-2006, and exceed previous values by 40% and more. The data base includes the glacier margins, elevations models as they have been compiled within the thesis (DEMs of 1997 and 2006 are part of the glacier inventories, length changes are part of the length change data base of the Austrian Alpine Club).

Age models of ODP Leg 175 holes

A multiproxy approach including the use of stable isotopes, magnetic characterization analyses, and organic geochemistry has been adopted to consider factors such as productivity and terrigenous input over the past 1.5 m.y. at two areas off the western coast of Africa. These factors can, in turn, be used to consider variability in ocean circulation and upwelling in addition to changes in climate on the African continent. In particular, studies focused on the influence of glacial-interglacial cycles and evidence for the mid-Pleistocene revolution (MPR), a complex change in climate that occurred at ~1 Ma. A comparison of the records from the two areas drilled during Ocean Drilling Program Leg 175, the Congo Basin, at a latitude of 5°S (Holes 1076A and 1077A), and the Walvis Ridge, at 17°S (Hole 1081A), demonstrates that these sites are affected by different localized factors. The sites in the Congo Basin are strongly influenced by freshwater and sediment from the Congo River, whereas the site at the Walvis Ridge is located in the center of oceanic upwelling and contains a more marine signal. Evidence also suggests that the two sites responded differently to both long- and short-term climatic variations. In particular, the response at the Walvis Ridge to the MPR occurred over an extended period, from 1.1 to 0.8 Ma, and was associated with a change in the dominant source of terrigenous input to the site in conjunction with a change in the productivity signal. In the Congo Basin, the response to the MPR was more rapid, occurring between 0.9 and 0.8 Ma. During this period, the influence of the Congo River became significant. However, productivity records only began to respond toward the end of this interval, at 0.8 Ma.

Aprubt climate change events in the Northeastern Atlantic Ocean

The western Iberian margin has been one of the key locations to study abrupt glacial climate change and associated interhemispheric linkages. The regional variability in the response to those events is being studied by combining a multitude of published and new records. Looking at the trend from Marine Isotope Stage (MIS) 10 to 2, the planktic foraminifer data, conform with the alkenone record of Martrat et al. [2007], shows that abrupt climate change events, especially the Heinrich events, became more frequent and their impacts in general stronger during the last glacial cycle. However, there were two older periods with strong impacts on the Atlantic meridional overturning circulation (AMOC): the Heinrich-type event associated with Termination (T) IV and the one occurring during MIS 8 (269 to 265 ka). During the Heinrich stadials of the last glacial cycle, the polar front reached the northern Iberian margin (ca. 41°N), while the arctic front was located in the vicinity of 39°N. During all the glacial periods studied, there existed a boundary at the latter latitude, either the arctic front during extreme cold events or the subarctic front during less strong coolings or warmer glacials. Along with these fronts sea surface temperatures (SST) increased southward by about 1°C per one degree of latitude leading to steep temperature gradients in the eastern North Atlantic and pointing to a close vicinity between subpolar and subtropical waters. The southern Iberian margin was always bathed by subtropical water masses - surface and/ or subsurface ones -, but there were periods when these waters also penetrated northward to 40.6°N. Glacial hydrographic conditions were similar during MIS 2 and 4, but much different during MIS 6. MIS 6 was a warmer glacial with the polar front being located further to the north allowing the subtropical surface and subsurface waters to reach at minimum as far north as 40.6°N and resulting in relative stable conditions on the southern margin. In the vertical structure, the Greenland-type climate oscillations during the last glacial cycle were recorded down to 2465 m during the Heinrich stadials, i.e. slightly deeper than in the western basin. This deeper boundary is related to the admixing of Mediterranean Outflow Water, which also explains the better ventilation of the intermediate-depth water column on the Iberian margin. This compilation revealed that latitudinal, longitudinal and vertical gradients existed in the waters along the Iberian margin, i.e. in a relative restricted area, but sufficient paleo-data exists now to validate regional climate models for abrupt climate change events in the northeastern North Atlantic Ocean.

Age models for ODPO Leg 130 sites

The Neogene carbonate stratigraphy of five sites drilled on Ontong Java Plateau during Leg 130 reveals a number of patterns which are unexpected, and which we refer to as loss paradox, equatorial insensitivity, and climate paradox. They denote the following unresolved questions. 1 The loss of carbonate at depth (as derived from differences in accumulation rates) is much greater than suggested by the change in carbonate percentages (calculated under the assumption that carbonate dissolution is the cause of loss). This indicates an important role for redeposition processes, such as winnowing (bottom currents), sifting (resuspension and catabatic flow) and episodic sloughing or solifluction (presumably stimulated by earthquakes). 2 Accumulation rates are not markedly increased at the time a site crosses the equator. There are several possible reasons. Equatorial upwelling may be unimportant in controlling sedimentation rates this far in the western Pacific, or its output may be spread over a considerable distance from the equator. Alternatively, increased supply below the equator is compensated for by increased removal (e.g. from resuspension by bioturbation, combined with catabatic flow). It is conceivable that errors in the timescale could also produce the effect seen. 3 There is an overall tendency for agreement between the stratigraphic patterns of carbonate content and of accumulation rates, but neither pattern is readily explained by reference to changes in climate (represented by benthic delta18O) or in sea-level (as derived from sequence stratigraphy).

Age models and stable isotope ratios of SO12 and SO35 sediment cores

Based on organic carbon accumulation rates, nine time slices of oceanic export paleoproductivity (Pnew) are presented which depict the variability of Pnew on a global scale through the last 30,000 years and document that the basic distribution patterns did not change through glacial and interglacial times. However, the glacial ocean shows an increased contrast of high- versus low-productivity zones. d13C values of near-surface-dwelling planktonic foraminifera Globigerinoides ruber suggest that the same contrast applies to the glacial nutrient inventories of the ambient surface waters, with a significant glacial transfer of PO4 from low- to high-productivity zones. In this way, glacial Pnew increased by a global average of about 2-4 Gt C/yr and led, via an enhanced CaCO3 dissolution and alkalinity in the deep ocean, to a significant extraction of CO2 from the surface water and the atrnosphere.

Gridded climate data from 5 Global Climate Models (GCM) of the Last Glacial Maximum (LGM) downscaled to 30 arc seconds for Europe, with links to NetCDF files

Studies on the impact of historical, current and future global change require very high-resolution climate data (less or equal 1km) as a basis for modelled responses, meaning that data from digital climate models generally require substantial rescaling. Another shortcoming of available datasets on past climate is that the effects of sea level rise and fall are not considered. Without such information, the study of glacial refugia or early Holocene plant and animal migration are incomplete if not impossible. Sea level at the last glacial maximum (LGM) was approximately 125m lower, creating substantial additional terrestrial area for which no current baseline data exist. Here, we introduce the development of a novel, gridded climate dataset for LGM that is both very high resolution (1km) and extends to the LGM sea and land mask. We developed two methods to extend current terrestrial precipitation and temperature data to areas between the current and LGM coastlines. The absolute interpolation error is less than 1°C and 0.5 °C for 98.9% and 87.8% of all pixels for the first two 1 arc degree distance zones. We use the change factor method with these newly assembled baseline data to downscale five global circulation models of LGM climate to a resolution of 1km for Europe. As additional variables we calculate 19 'bioclimatic' variables, which are often used in climate change impact studies on biological diversity. The new LGM climate maps are well suited for analysing refugia and migration during Holocene warming following the LGM.