Major and trace element content of the 17 lithostratigraphic units recovered during DSDP Leg 49 (Table 1)

IPOD Leg 49 recovered basalts from 9 holes at 7 sites along 3 transects across the Mid-Atlantic Ridge: 63°N (Reykjanes), 45°N and 36°N (FAMOUS area). This has provided further information on the nature of mantle heterogeneity in the North Atlantic by enabling studies to be made of the variation of basalt composition with depth and with time near critical areas (Iceland and the Azores) where deep mantle plumes are thought to exist. Over 150 samples have been analysed for up to 40 major and trace elements and the results used to place constraints on the petrogenesis of the erupted basalts and hence on the geochemical nature of their source regions. It is apparent that few of the recovered basalts have the geochemical characteristics of typical "depleted" midocean ridge basalts (MORB). An unusually wide range of basalt compositions may be erupted at a single site: the range of rare earth patterns within the short section cored at Site 413, for instance, encompasses the total variation of REE patterns previously reported from the FAMOUS area. Nevertheless it is possible to account for most of the compositional variation at a single site by partial melting processes (including dynamic melting) and fractional crystallization. Partial melting mechanisms seem to be the dominant processes relating basalt compositions, particularly at 36°N and 45°N, suggesting that long-lived sub-axial magma chambers may not be a consistent feature of the slow-spreading Mid-Atlantic Ridge. Comparisons of basalts erupted at the same ridge segment for periods of the order of 35 m.y. (now lying along the same mantle flow line) do show some significant inter-site differences in Rb/Sr, Ce/Yb, 87Sr/86Sr, etc., which cannot be accounted for by fractionation mechanisms and which must reflect heterogeneities in the mantle source. However when hygromagmatophile (HYG) trace element levels and ratios are considered, it is the constancy or consistency of these HYG ratios which is the more remarkable, implying that the mantle source feeding a particular ridge segment was uniform with respect to these elements for periods of the order of 35 m.y. and probably since the opening of the Atlantic. Yet these HYG element ratios at 63°N are very different from those at 45°N and 36°N and significantly different from the values at 22°N and in "MORB". The observed variations are difficult to reconcile with current concepts of mantle plumes and binary mixing models. The mantle is certainly heterogeneous, but there is not simply an "enriched" and a "depleted" source, but rather a range of sources heterogeneous on different scales for different elements - to an extent and volume depending on previous depletion/enrichment events. HYG element ratios offer the best method of defining compositionally different mantle segments since they are little modified by the fractionation processes associated with basalt generation.

Bulk and molecular proxies of sediment core GeoB6518-1

A 20 kyr long sediment sequence from the Congo deep sea fan (core GeoB 6518-1), one of the world's largest deep sea river fans, has been analysed for bulk and molecular proxies in order to reconstruct the marine, soil and plant organic carbon (OC) contributions to these sediments since the last glacial maximum. The bulk proxies applied, C/N ratio and d13Corg, ranged from 10 to 12.5 and from -24.5 to -21 per mill VPDB, respectively. As molecular proxies, concentrations of marine derived alkenones and terrestrial derived odd-numbered n-alkanes were used, which varied between 0.2 and 4 µg/g dry weight sediment. In addition, the branched vs. isoprenoid tetraether (BIT) index, a proxy for soil organic matter input, was used, which varied from 0.3 to 0.5 in this core. Application of binary mixing models, based on the different individual proxies, showed estimates for terrestrial OC input varying by up to 50% due to the heterogeneous nature of the OC. Application of a three end-member mixing model using the d13Corg content, the C/N ratio and the BIT index, enabled the distinction of soil and plant organic matter as separate contributors to the sedimentary OC pool. The results show that marine OC accounts for 20% to 40% of the total OC present in the deep sea fan sediments over the last 20 kyr and that soil OC accounts for about half (45% on average) of the OC present. This suggests that soil OC represents the majority of the terrestrial OC delivered to the fan sediments. Accumulation rates of the plant and soil OC fractions over the last 20 kyr varied by a factor of up to 5, and are strongly related to sediment accumulation rates. They showed an increase starting at ca. 17 kyr BP, a decline during the Younger Dryas, peak values during the early Holocene and lower values in the late Holocene. This pattern matches with reconstructions of past central African humidity and Congo River discharge from the same core and revealed that central African precipitation patterns exert a dominant control on terrestrial OC deposition in the Congo deep sea fan. Marine OC accumulation rates are only weakly related to sediment accumulation rates and vary only little over time compared to the terrigenous fractions. These variations are likely a result of enhanced preservation during times of higher sedimentation rates and of relative small fluctuations in primary production due to wind-driven upwelling.

n-alkane carbon isotope ratios in plants and recent lake sediments from Cameroon

Trees and shrubs in tropical Africa use the C3 cycle as a carbon fixation pathway during photosynthesis, while grasses and sedges mostly use the C4 cycle. Leaf-wax lipids from sedimentary archives such as the long-chain n-alkanes (e.g., n-C27 to n-C33) inherit carbon isotope ratios that are representative of the carbon fixation pathway. Therefore, n-alkane d13C values are often used to reconstruct past C3/C4 composition of vegetation, assuming that the relative proportions of C3 and C4 leaf waxes reflect the relative proportions of C3 and C4 plants. We have compared the d13C values of n-alkanes from modern C3 and C4 plants with previously published values from recent lake sediments and provide a framework for estimating the fractional contribution (areal-based) of C3 vegetation cover (fC3) represented by these sedimentary archives. Samples were collected in Cameroon, across a latitudinal transect that accommodates a wide range of climate zones and vegetation types, as reflected in the progressive northward replacement of C3-dominated rain forest by C4-dominated savanna. The C3 plants analysed were characterised by substantially higher abundances of n-C29 alkanes and by substantially lower abundances of n-C33 alkanes than the C4 plants. Furthermore, the sedimentary d13C values of n-C29 and n-C31 alkanes from recent lake sediments in Cameroon (-37.4 per mil to -26.5 per mil) were generally within the range of d13C values for C3 plants, even when from sites where C4 plants dominated the catchment vegetation. In such cases simple linear mixing models fail to accurately reconstruct the relative proportions of C3 and C4 vegetation cover when using the d13C values of sedimentary n-alkanes, overestimating the proportion of C3 vegetation, likely as a consequence of the differences in plant wax production, preservation, transport, and/or deposition between C3 and C4 plants. We therefore tested a set of non-linear binary mixing models using d13C values from both C3 and C4 vegetation as end-members. The non-linear models included a sigmoid function (sine-squared) that describes small variations in the fC3 values as the minimum and maximum d13C values are approached, and a hyperbolic function that takes into account the differences between C3 and C4 plants discussed above. Model fitting and the estimation of uncertainties were completed using the Monte Carlo algorithm and can be improved by future data addition. Models that provided the best fit with the observed d13C values of sedimentary n-alkanes were either hyperbolic functions or a combination of hyperbolic and sine-squared functions. Such non-linear models may be used to convert d13C measurements on sedimentary n-alkanes directly into reconstructions of C3 vegetation cover.

Age models and foraminifer isotopes of four ODP Leg 202 sites

We present benthic isotope stratigraphies for Sites 1236, 1237, 1239, and 1241 that span the late Miocene-Pliocene time interval from 6 to 2.4 Ma. Orbitally tuned timescales were generated for Sites 1237 and 1241 by correlating the high-frequency variations in gamma ray attenuation density, percent sand of the carbonate fraction, and benthic d13C to variations in Earth's orbital parameters. The astronomical timescales for Sites 1237 and 1241 are in agreement with the one from Atlantic Site 925/926 (Ocean Drilling Program Leg 154). The comparison of benthic d18O and d13C records from the east Pacific sites and Atlantic Site 925/926 revealed a surprising clarity of the "41-k.y. signal" in d13C records and a remarkably good correlation between their d13C records. This suggests that the late Miocene-Pliocene amplitudes of obliquity-related d13C cycles reflect a magnitude of global response often larger than that provided by obliquity-related d18O cycles. At Site 1237, the orbitally derived ages of Pliocene magnetic reversal boundaries between the base of Réunion and the top of Thvera confirm astronomical datings of the generally accepted ATNTS2004 timescale, except for the top of Kaena and the base of Sidufjall. Our astronomical age for the top of Kaena is about one obliquity cycle older. The base of Sidufjall appears to be about one precession cycle younger. The age models of Sites 1236 and 1239 were established by correlating their benthic d18O and d13C records directly to the orbitally tuned isotope record of Site 1241.

Age models and stable isotope analysis on sediment core Site 199-1218 from the equatorial Pacific

A 13-million-year continuous record of Oligocene climate from the equatorial Pacific reveals a pronounced "heartbeat" in the global carbon cycle and periodicity of glaciations. This heartbeat consists of 405,000-, 127,000-, and 96,000-year eccentricity cycles and 1.2-million-year obliquity cycles in periodically recurring glacial and carbon cycle events. That climate system response to intricate orbital variations suggests a fundamental interaction of the carbon cycle, solar forcing, and glacial events. Box modeling shows that the interaction of the carbon cycle and solar forcing modulates deep ocean acidity as well as the production and burial of global biomass. The pronounced 405,000-year eccentricity cycle is amplified by the long residence time of carbon in the oceans.

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.

Age models, organic analyses and major-element analyses of sediment cores on a transect of the tropical African rainbelt

The distribution of rainfall in tropical Africa is controlled by the African rainbelt**1, which oscillates on a seasonal basis. The rainbelt has varied on centennial to millennial timescales along with changes in Northern Hemisphere high-latitude climate**2, 3, 4, 5, the Atlantic meridional overturning circulation**6 and low-latitude insolation**7 over the past glacial-interglacial cycle. However, the overall dynamics of the African rainbelt remain poorly constrained and are not always consistent with a latitudinal migration**2, 4, 5, 6, as has been proposed for other regions**8, 9. Here we use terrestrially derived organic and sedimentary markers from marine sediment cores to reconstruct the distribution of vegetation, and hence rainfall, in tropical Africa during extreme climate states over the past 23,000 years. Our data indicate that rather than migrating latitudinally, the rainbelt contracted and expanded symmetrically in both hemispheres in response to changes in climate. During the Last Glacial Maximum and Heinrich Stadial 1, the rainbelt contracted relative to the late Holocene, which we attribute to a latitudinal compression of atmospheric circulation associated with lower global mean temperatures**10. Conversely, during the mid-Holocene climatic optimum, the rainbelt expanded across tropical Africa. In light of our findings, it is not clear whether the tropical rainbelt has migrated latitudinally on a global scale, as has been suggested**8,9.

Surface velocity fields, digital elevation models and ice front positions derived from multi-mission SAR remote sensing data at Sjögren Inlet glaciers, Antarctic Peninsula

Substantial retreat or disintegration of numerous ice shelves have been observed on the Antarctic Peninsula. The ice shelf in the Prince Gustav Channel retreated gradually since the late 1980's and broke-up in 1995. Tributary glaciers reacted with speed-up, surface lowering and increased ice discharge, consequently contributing to sea level rise. We present a detailed long-term study (1993-2014) on the dynamic response of Sjögren Inlet glaciers to the disintegration of Prince Gustav Ice Shelf. We analyzed various remote sensing datasets to observe the reactions of the glaciers to the loss of the buttressing ice shelf. A strong increase in ice surface velocities was observed with maximum flow speeds reaching 2.82±0.48 m/d in 2007 and 1.50±0.32 m/d in 2004 at Sjögren and Boydell glaciers respectively. Subsequently, the flow velocities decelerated, however in late 2014, we still measured about two times the values of our first measurements in 1996. The tributary glaciers retreated 61.7±3.1 km² behind the former grounding line of the ice shelf. In regions below 1000 m a.s.l., a mean surface lowering of -68±10 m (-3.1 m/a) was observed in the period 1993-2014. The lowering rate decreased to -2.2 m/a in recent years. Based on the surface lowering rates, geodetic mass balances of the glaciers were derived for different time steps. High mass loss rate of -1.21±0.36 Gt/a was found in the earliest period (1993-2001). Due to the dynamic adjustments of the glaciers to the new boundary conditions the ice mass loss reduced to -0.59±0.11 Gt/a in the period 2012-2014, resulting in an average mass loss rate of -0.89±0.16 Gt/a (1993-2014). Including the retreat of the ice front and grounding line, a total mass change of -38.5±7.7 Gt and a contribution to sea level rise of 0.061±0.013 mm were computed. Analysis of the ice flux revealed that available bedrock elevation estimates at Sjögren Inlet are too shallow and are the major uncertainty in ice flux computations. This temporally dense time series analysis of Sjögren Inlet glaciers shows that the adjustments of tributary glaciers to ice shelf disintegration are still going on and provides detailed information of the changes in glacier dynamics.

Age models for sediment cores off Peru

Planktonic foraminiferal faunas of the southeast Pacific indicate that sea surface temperatures (SST) have varied by as much as 8-10°C in the Peru Current, and by ?5-7°C along the equator, over the past 150,000 years. Changes in SST at times such as the Last Glacial Maximum reflect incursion of high-latitude species Globorotalia inflata and Neogloboquadrina pachyderma into the eastern boundary current and as far north as the equator. A simple heat budget model of the equatorial Pacific shows that observed changes in Peru Current advection can account for about half of the total variability in equatorial SSTs. The remaining changes in equatorial SST, which are likely related to local changes in upwelling or pycnocline depth, precede changes in polar climates as recorded by d18O. This partitioning of processes in eastern equatorial Pacific SST reveals that net ice-age cooling here reflects first a rapid response of equatorial upwelling to insolation, followed by a later response to changes in the eastern boundary current associated with high-latitude climate (which closely resembles variations in atmospheric CO2 as recorded in the Vostok ice core). Although precise mechanisms responsible for the equatorial upwelling component of climate change remain uncertain, one likely candidate that may operate independently of the ice sheets is insolation-driven changes in El Niño/Southern Oscillation (ENSO) frequency. Early responses of equatorial SST detected both here and elsewhere highlight the sensitivity of tropical systems to small changes in seasonal insolation. The scale of tropical changes we have observed are substantially greater than model predictions, suggesting a need for further quantitative assessment of processes associated with long-term climate change.

Surface velocity fields, digital elevation models, ice front positions and grounding line derived from remote sensing data at Dinsmoor-Bombardier-Edgeworth glacier system, Antarctic Peninsula (1992-2014)

The northern Antarctic Peninsula is one of the fastest changing regions on Earth. The disintegration of the Larsen-A Ice Shelf in 1995 caused tributary glaciers to adjust by speeding up, surface lowering, and overall increased ice-mass discharge. In this study, we investigate the temporal variation of these changes at the Dinsmoor-Bombardier-Edgeworth glacier system by analyzing dense time series from various spaceborne and airborne Earth observation missions. Precollapse ice shelf conditions and subsequent adjustments through 2014 were covered. Our results show a response of the glacier system some months after the breakup, reaching maximum surface velocities at the glacier front of up to 8.8 m/d in 1999 and a subsequent decrease to ~1.5 m/d in 2014. Using a dense time series of interferometrically derived TanDEM-X digital elevation models and photogrammetric data, an exponential function was fitted for the decrease in surface elevation. Elevation changes in areas below 1000 m a.s.l. amounted to at least 130±15 m130±15 m between 1995 and 2014, with change rates of ~3.15 m/a between 2003 and 2008. Current change rates (2010-2014) are in the range of 1.7 m/a. Mass imbalances were computed with different scenarios of boundary conditions. The most plausible results amount to -40.7±3.9 Gt-40.7±3.9 Gt. The contribution to sea level rise was estimated to be 18.8±1.8 Gt18.8±1.8 Gt, corresponding to a 0.052±0.005 mm0.052±0.005 mm sea level equivalent, for the period 1995-2014. Our analysis and scenario considerations revealed that major uncertainties still exist due to insufficiently accurate ice-thickness information. The second largest uncertainty in the computations was the glacier surface mass balance, which is still poorly known. Our time series analysis facilitates an improved comparison with GRACE data and as input to modeling of glacio-isostatic uplift in this region. The study contributed to a better understanding of how glacier systems adjust to ice shelf disintegration.

Response of A. aurita ephyrae to three stressors: high-temperature, high-CO2 and low-oxygen

Global change is affecting marine ecosystems through a combination of different stressors such as warming, ocean acidification and oxygen depletion. Very little is known about the interactions among these factors, especially with respect to gelatinous zooplankton. Therefore, in this study we investigated the direct effects of pH, temperature and oxygen availability on the moon jellyfish Aurelia aurita, concentrating on the ephyral life stage. Starved one-day-old ephyrae were exposed to a range of pCO2 (400-4000 ppm) and three different dissolved oxygen levels (from saturated to hypoxic conditions), in two different temperatures (5 and 15 °C) for 7 days. Carbon content and swimming activity were analysed at the end of the incubation period, and mortality noted. General linearized models were fitted through the data, with the best fitting models including two- and three-way interactions between pCO2, temperature and oxygen concentration. The combined effect of the stressors was small but significant, with the clearest negative effect on growth caused by the combination of all three stressors present (high temperature, high CO2, low oxygen). We conclude that A. aurita ephyrae are robust and that they are not likely to suffer from these environmental stressors in a near future.

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.