A high-resolution multiproxy study of sediment core GeoB7186-3 from SE Pacific off southern South America

A high-resolution multiproxy study performed on a marine record from SE Pacific off southern South America was used to reconstruct past regional environmental changes and their relation to global climate, particularly to El Niño/Southern Oscillation (ENSO) phenomenon during the last 2200 years. Our results suggest a sustained northward shift in the position of the zonal systems, i.e. the Southern Westerly Wind belt and the Antarctic Circumpolar Current, which occurred between 1300 and 750 yr BP. The synchrony of the latitudinal shift with cooling in Antarctica and reduced ENSO activity observed in several marine and terrestrial archives across South America suggests a causal link between ENSO and the proposed displacement of the zonal systems. This shift might have acted as a positive feedback to more La Niña-like conditions between 1300 and 750 yr BP by steepening the hemispheric and tropical Pacific zonal sea surface temperature gradient. This scenario further suggests different boundary conditions for ENSO before 1300 and after 750 yr BP.

Fly-through movie for the reconstructed folds in the onset region of the Petermann Glacier, North Greenland

The increasing catalogue of high-quality ice-penetrating radar data provides a unique insight in the internal layering architecture of the Greenland ice sheet. The stratigraphy, an indicator of past deformation, highlights irregularities in ice flow and reveals large perturbations without obvious links to bedrock shape. In this work, to establish a new conceptual model for the formation process, we analysed the radar data at the onset of the Petermann Glacier, North Greenland, and created a three-dimensional model of several distinct stratigraphic layers. We demonstrate that the dominant structures are cylindrical folds sub-parallel to the ice flow. By numerical modelling, we show that these folds can be formed by lateral compression of mechanically anisotropic ice, while a general viscosity contrast between layers would not lead to folding for the same boundary conditions. We conclude that the folds primarily form by converging flow as the mechanically anisotropic ice is channelled towards the glacier.

Stable carbon and oxygen isotope ratios of foraminifera from Holocene sediments of the North Atlantic

The Denmark Strait Overflow (DSO) today compensates for the northward flowing Norwegian and Irminger branches of the North Atlantic Current that drive the Nordic heat pump. During the Last Glacial Maximum (LGM), ice sheets constricted the Denmark Strait aperture in addition to ice eustatic/isostatic effects which reduced its depth (today ~630 m) by ~130 m. These factors, combined with a reduced north-south density gradient of the water-masses, are expected to have restricted or even reversed the LGM DSO intensity. To better constrain these boundary conditions, we present a first reconstruction of the glacial DSO, using four new and four published epibenthic and planktic stable-isotope records from sites to the north and south of the Denmark Strait. The spatial and temporal distribution of epibenthic delta18O and delta13C maxima reveals a north-south density gradient at intermediate water depths from sigma0 ~28.7 to 28.4/28.1 and suggests that dense and highly ventilated water was convected in the Nordic Seas during the LGM. However, extremely high epibenthic delta13C values on top of the Mid-Atlantic Ridge document a further convection cell of Glacial North Atlantic Intermediate Water to the south of Iceland, which, however, was marked by much lower density (sigma0 ~28.1). The north-south gradient of water density possibly implied that the glacial DSO was directed to the south like today and fed Glacial North Atlantic Deep Water that has underthrusted the Glacial North Atlantic Intermediate Water in the Irminger Basin.

Geochemistry of sediments cores in the Mediterranean Sea

We investigated five time-equivalent core sections (180-110 kyr BP) from the Balearic Sea (Menorca Rise), the easternmost Levantine Basin and southwest, south, and southeast of Crete to reconstruct spatial patterns of productivity during deposition of sapropels S5 and S6 in the Mediterranean Sea. Our indicators are Ba, total organic carbon and carbonate contents. We found no indications of Ba remobilization within the investigated core intervals, and used the accumulation rate of biogenic Ba to compute paleoproductivity. Maximum surface water productivity (up to 350 g C/m2/yr) was found during deposition of S5 (isotope stage 5e) but pronounced spatial variability is evident. Coeval sediment intervals in the Balearic Sea show very little productivity change, suggesting that chemical and biological environments in the eastern and western Mediterranean basins were decoupled in this interval. We interpret the spatial variability as the result of two different modes of nutrient delivery to the photic zone: riverderived nutrient input and shoaling of the pycnocline/nutricline to the photic zone. The productivity increase during the formation of S6 was moderate compared to S5 and had a less marked spatial variability within the study area of the eastern Mediterranean Sea. Given that S6 formed during a glacial interval, glacial boundary conditions such as high wind stress and/or cooler surface water temperatures apparently favored lateral and vertical mixing and prevented the development of the spatial gradients within the Eastern Mediterranean Sea (EMS) observed for S5. A non-sapropel sediment interval with elevated Ba content and depleted 18O/16O ratios in planktonic foraminifer calcite was detected between S6 and S5 that corresponds to the weak northern hemisphere insolation maximum at 150 kyr. At this time, productivity apparently increased up to five times over surrounding intervals, but abundant benthic fauna show that the deep water remained oxic. Following our interpretation, the interval denotes a failed sapropel, when a weaker monsoon did not force the EMS into permanent stratification. The comparison of interglacial and glacial sapropels illustrates the relevance of climatic boundary conditions in the northern catchment in determining the facies and spatial variability of sapropels within the EMS.

Hydroclimate records and radiocarbon dating of sediment cores SO189/2 _39KL, SO189/2 _119KL, and SO189/2 _144KL from the eastern tropical Indian Ocean

The response of the tropical climate in the Indian Ocean realm to abrupt climate change events in the North Atlantic Ocean is contentious. Repositioning of the intertropical convergence zone is thought to have been responsible for changes in tropical hydroclimate during North Atlantic cold spells1, 2, 3, 4, 5, but the dearth of high-resolution records outside the monsoon realm in the Indian Ocean precludes a full understanding of this remote relationship and its underlying mechanisms. Here we show that slowdowns of the Atlantic meridional overturning circulation during Heinrich stadials and the Younger Dryas stadial affected the tropical Indian Ocean hydroclimate through changes to the Hadley circulation including a southward shift in the rising branch (the intertropical convergence zone) and an overall weakening over the southern Indian Ocean. Our results are based on new, high-resolution sea surface temperature and seawater oxygen isotope records of well-dated sedimentary archives from the tropical eastern Indian Ocean for the past 45,000 years, combined with climate model simulations of Atlantic circulation slowdown under Marine Isotope Stages 2 and 3 boundary conditions. Similar conditions in the east and west of the basin rule out a zonal dipole structure as the dominant forcing of the tropical Indian Ocean hydroclimate of millennial-scale events. Results from our simulations and proxy data suggest dry conditions in the northern Indian Ocean realm and wet and warm conditions in the southern realm during North Atlantic cold spells.

Stable carbon and oxygen isotope record of IODP Site 306-U1313, MIS 23-19

Since the seminal work by Hays et al. (1976), a plethora of studies has demonstrated a correlation between orbital variations and climatic change. However, information on how changes in orbital boundary conditions affected the frequency and amplitude of millennial-scale climate variability is still fragmentary. The Marine Isotope Stage (MIS) 19, an interglacial centred at around 785 ka, provides an opportunity to pursue this question and test the hypothesis that the long-term processes set up the boundary conditions within which the short-term processes operate. Similarly to the current interglacial, MIS 19 is characterised by a minimum of the 400-kyr eccentricity cycle, subdued amplitude of precessional changes, and small amplitude variations in insolation. Here we examine the record of climatic conditions during MIS 19 using high-resolution stable isotope records from benthic and planktonic foraminifera from a sedimentary sequence in the North Atlantic (Integrated Ocean Drilling Program Expedition 306, Site U1313) in order to assess the stability and duration of this interglacial, and evaluate the climate system's response in the millennial band to known orbitally induced insolation changes. Benthic and planktonic foraminiferal d18O values indicate relatively stable conditions during the peak warmth of MIS 19, but sea-surface and deep-water reconstructions start diverging during the transition towards the glacial MIS 18, when large, cold excursions disrupt the surface waters whereas low amplitude millennial scale fluctuations persist in the deep waters as recorded by the oxygen isotope signal. The glacial inception occurred at ~779 ka, in agreement with an increased abundance of tetra-unsaturated alkenones, reflecting the influence of icebergs and associated meltwater pulses and high-latitude waters at the study site. After having combined the new results with previous data from the same site, and using a variety of time series analysis techniques, we evaluate the evolution of millennial climate variability in response to changing orbital boundary conditions during the Early-Middle Pleistocene. Suborbital variability in both surface- and deep-water records is mainly concentrated at a period of ~11 kyr and, additionally, at ~5.8 and ~3.9 kyr in the deep ocean; these periods are equal to harmonics of precession band oscillations. The fact that the response at the 11 kyr period increased over the same interval during which the amplitude of the response to the precessional cycle increased supports the notion that most of the variance in the 11 kyr band in the sedimentary record is nonlinearly transferred from precession band oscillations. Considering that these periodicities are important features in the equatorial and intertropical insolation, these observations are in line with the view that the low-latitude regions play an important role in the response of the climate system to the astronomical forcing. We conclude that the effect of the orbitally induced insolation is of fundamental importance in regulating the timing and amplitude of millennial scale climate variability.

Effective stress, porosity, p-wave velocity and mineral composition of ODP Hole 174A-1073A sediments

Porosity, permeability, and compressional (P-wave) velocity were measured as a function of stress on sediments from Ocean Drilling Program Site 1073, U.S. Mid-Atlantic continental slope. Thin sections, scanning electron microscopy, and X-ray diffraction analyses provided mineralogical characteristics of the samples. Uniaxial strain boundary conditions were imposed on the samples during consolidation tests with the maximum effective axial stress reaching 13 MPa. The maximum effective radial stress necessary to maintain uniaxial strain was 7.6 MPa. Over an effective axial stress interval of 0 to 5.2 MPa, Sample 174A-1073A-26X-2, 82-89 cm (226.65 meters below seafloor [mbsf]), exhibited the largest decrease in porosity (51% to 41%), whereas Sample 71X-1, 2-8 cm (644.70 mbsf), exhibited the smallest decrease in porosity (48% to 45%). All samples showed negligible porosity increases during unloading. The permeability (on the order of 1 x 10-17 m**2) of Sample 174A-1073A-71X-1, 2-8 cm, was twice that measured on Sample 8H-1, 23-26 cm (63.75 mbsf), even though the former was considerably deeper and older. The differences in porosity-stress behavior and permeability between shallow and deep samples is related to lithologic, mineralogic, and diagenetic differences between the sediments above and below the Pliocene-Pleistocene to Miocene unconformity. P-wave velocity for Samples 174A-1073A-41X-5, 97-103 cm (372.35 mbsf), and 71X-1, 2-8 cm, increased with decreasing porosity, but did not change significantly during unloading.

The influence of seawater pH on U/ Ca ratios in the scleractinian cold-water coral Lophelia pertusa

The increasing pCO2 in seawater is a serious threat for marine calcifiers and alters the biogeochemistry of the ocean. Therefore, the reconstruction of past-seawater properties and their impact on marine ecosystems is an important way to investigate the underlying mechanisms and to better constrain the effects of possible changes in the future ocean. Cold-water coral (CWC) ecosystems are biodiversity hotspots. Living close to aragonite undersaturation, these corals serve as living laboratories as well as archives to reconstruct the boundary conditions of their calcification under the carbonate system of the ocean. We investigated the reef-building CWC Lophelia pertusa as a recorder of intermediate ocean seawater pH. This species-specific field calibration is based on a unique sample set of live in situ collected L. pertusa and corresponding seawater samples. These data demonstrate that uranium speciation and skeletal incorporation for azooxanthellate scleractinian CWCs is pH dependent and can be reconstructed with an uncertainty of ±0.15. Our Lophelia U / Ca-pH calibration appears to be controlled by the high pH values and thus highlighting the need for future coral and seawater sampling to refine this relationship. However, this study recommends L. pertusa as a new archive for the reconstruction of intermediate water mass pH and hence may help to constrain tipping points for ecosystem dynamics and evolutionary characteristics in a changing ocean.

A 300 kyr record of upwelling off Oman during the Late Quaternary: evidence of the Asian southwest monsoon

In the northwest Arabian Sea upwelling occurs each summer, driven by the strong SW monsoon winds. Upwelling results in high biological productivity and a distinctive assemblage of plankton species in the surface waters off Oman that are preserved in the sediments along the Oman continental margin, creating a geologic record of monsoon-driven upwelling. Sediments recovered from the Oman continental margin during Ocean Drilling Program leg 117 provide an opportunity to examine how upwelling has varied during the late Quaternary, spanning a longer interval than piston cores recovered prior to the ODP cruise. Variations in foraminifer shell accumulation and in the relative abundance of Globigerina bulloides indicate dominant cycles of variation at 1/100 kyr, the dominant frequency of glacial-interglacial variations, and at 1/23 kyr, the frequency of precessionally driven cycles in seasonal insolation. The strongest monsoon winds (indicated by increased upwelling) occurred during interglacial times when perihelion was aligned with the summer solstice, an orbital change that increased the insolation received during summer in the northern hemisphere. During glacial times upwelling was reduced, and although the precessional cycles were still present their amplitude was smaller. At both frequencies the upwelling cycles are in phase with minimum ice volume, evidence that glacial-interglacial climate changes also include changes to the climate system that influence the low-latitude monsoon. We attribute the decrease in the monsoon winds observed during glacial times to changes in bare land albedo over Asia and/or to changes in the areal extent and seasonal cycle in Asian snow cover that decrease the summer land-sea temperature contrast. Other mechanisms may also be involved. These new upwelling time series differ substantially from previous results, however the previous work relied on cores located farther offshore where upwelling is less intense and other physical mechanisms become important. Our results support the observations derived from atmospheric general circulation models of the atmosphere that indicate that both glacial boundary conditions, and the strength of summer insolation are important variables contributing to cycles in the monsoon winds during the late Quaternary.

Data compilation of ESOP

The work in this sub-project of ESOP focuses on the advective and convective transforma-tion of water masses in the Greenland Sea and its neighbouring areas. It includes observational work on the sub-mesoscale and analysis of hydrographic data up to the gyre-scale. Observations of active convective plumes were made with a towed chain equipped with up to 80 CTD sensors, giving a horizontal and vertical resolution of the hydrographic fields of a few metres. The observed scales of the penetrative convective plumes compare well with those given by theory. On the mesoscale the structure of homogeneous eddies formed as a result of deep convection was observed and the associated mixing and renewal of the intermediate layers quantified. The relative importance and efficiency of thermal and haline penetrative convection in relation to the surface boundary conditions (heat and salt fluxes and ice cover) and the ambient stratification are studied using the multi year time series of hydro-graphic data in the central Greenland Sea. The modification of the water column of the Greenland Sea gyre through advection from and mixing with water at its rim is assessed on longer time scales. The relative contributions are quantified using modern water mass analysis methods based on inverse techniques. Likewise the convective renewal and the spreading of the Arctic Intermediate Water from its formation area is quantified. The aim is to budget the heat and salt content of the water column, in particular of the low salinity surface layer, and to relate its seasonal and interannual variability to the lateral fluxes and the fluxes at the air-sea-ice interface. This will allow to estimate residence times for the different layers of the Greenland Sea gyre, a quantity important for the description of the Polar Ocean carbon cycle.

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.

(Table 1) Main constituents of coccolith size fractions expressed as carbonate contribution at Bass River during the PETM

Coccoliths, calcite plates produced by the marine phytoplankton coccolithophores, have previously shown a large array of carbon and oxygen stable isotope fractionations (termed "vital effects"), correlated to cell size and hypothesized to reflect the varying importance of active carbon acquisition strategies. Culture studies show a reduced range of vital effects between large and small coccolithophores under high CO2, consistent with previous observations of a smaller range of interspecific vital effects in Paleocene coccoliths. We present new fossil data examining coccolithophore vital effects over three key Cenozoic intervals reflecting changing climate and atmospheric partial pressure of CO2 (pCO2). Oxygen and carbon stable isotopes of size-separated coccolith fractions dominated by different species from well preserved Paleocene-Eocene thermal maximum (PETM, ~56 Ma) samples show reduced interspecific differences within the greenhouse boundary conditions of the PETM. Conversely, isotope data from the Plio-Pleistocene transition (PPT; 3.5-2 Ma) and the last glacial maximum (LGM; ~22 ka) show persistent vital effects of ~2 per mil. PPT and LGM data show a clear positive trend between coccolith (cell) size and isotopic enrichment in coccolith carbonate, as seen in laboratory cultures. On geological timescales, the degree of expression of vital effects in coccoliths appears to be insensitive topCO2 changes over the range ~350 ppm (Pliocene) to ~180 ppm (LGM). The modern array of coccolith vital effects arose after the PETM but before the late Pliocene and may reflect the operation of more diverse carbon acquisition strategies in coccolithophores in response to decreasing Cenozoic pCO2.

Ba/Ca ratios in benthic foraminifera from recent core tops

The carbon isotope ratio (delta13C) and cadmium content (Cd/Ca) of benthic foraminifera shells have been used to reconstruct deep-water circulation patterns of the glacial oceans. These tracers co-vary with phosphorus in the modern ocean because they are nearly quantitatively regenerated from sinking biological debris in the upper water column. Hence they can be used to reconstruct the distribution of labile nutrients in glacial water masses. Independent constraints on glacial deep-ocean circulation patterns could be provided by a tracer of the distribution of silica and alkalinity, the deeply regenerated constituents of planktonic hard parts. Barium shares key aspects of its behaviour with these refractory nutrients because it is removed from solution in surface waters and incorporated into sinking particles which slowly dissolve deep in the water column and in the sediments. The fractionation of Ba between deep-water masses of the major ocean basins is largely controlled by thermohaline circulation patterns, so Ba conforms to different boundary conditions from Cd and 13C. As Ba substitutes into trigonal carbonates, it is a potential palaeoceano-graphic tracer if the Ba content of foraminifera shells reflects ambient dissolved Ba concentrations. Here we present data from Recent core-top benthic foraminifera which indicate that the Ba content of some recent calcitic benthic foraminifera does co-vary with bottom-water Ba.

Results of calculated surface air temperature anomaly of COSMOS Pliocene experiment 2 in NetCDF format

In this manuscript we describe the experimental procedure employed at the Alfred Wegener Institute in Germany in the preparation of the simulations for the Pliocene Model Intercomparison Project (PlioMIP). We present a description of the utilized Community Earth System Models (COSMOS, version: COSMOS-landveg r2413, 2009) and document the procedures that we applied to transfer the Pliocene Research, Interpretation and Synoptic Mapping (PRISM) Project mid-Pliocene reconstruction into model forcing fields. The model setup and spin-up procedure are described for both the paleo- and preindustrial (PI) time slices of PlioMIP experiments 1 and 2, and general results that depict the performance of our model setup for mid-Pliocene conditions are presented. The mid-Pliocene, as simulated with our COSMOS setup and PRISM boundary conditions, is both warmer and wetter in the global mean than the PI. The globally averaged annual mean surface air temperature in the mid-Pliocene standalone atmosphere (fully coupled atmosphere-ocean) simulation is 17.35 °C (17.82 °C), which implies a warming of 2.23 °C (3.40 °C) relative to the respective PI control simulation.