Stable isotope record and lithic fragments distribution in sediments of MIS9-14 of the North Atlantic

Stable isotope and ice-rafted debris records from three core sites in the mid-latitude North Atlantic (IODP Site U1313, MD01-2446, MD03-2699) are combined with records of ODP Sites 1056/1058 and 980 to reconstruct hydrographic conditions during the middle Pleistocene spanning Marine Isotope Stages (MIS) 9-14 (300-540 ka). Core MD03-2699 is the first high-resolution mid-Brunhes record from the North Atlantic's eastern boundary upwelling system covering the complete MIS 11c interval and MIS 13. The array of sites reflect western and eastern basin boundary current as well as north to south transect sampling of subpolar and transitional water masses and allow the reconstruction of transport pathways in the upper limb of the North Atlantic's circulation. Hydrographic conditions in the surface and deep ocean during peak interglacial MIS 9 and 11 were similar among all the sites with relative stable conditions and confirm prolonged warmth during MIS 11c also for the mid-latitudes. Sea surface temperature (SST) reconstructions further reveal that in the mid-latitude North Atlantic MIS 11c is associated with two plateaus, the younger one of which is slightly warmer. Enhanced subsurface northward heat transport in the eastern boundary current system, especially during early MIS 11c, is denoted by the presence of tropical planktic foraminifer species and raises the question how strongly it impacted the Portuguese upwelling system. Deep water ventilation at the onset of MIS 11c significantly preceded surface water ventilation. Although MIS 13 was generally colder and more variable than the younger interglacials the surface water circulation scheme was the same. The greatest differences between the sites existed during the glacial inceptions and glacials. Then a north - south trending hydrographic front separated the nearshore and offshore waters off Portugal. While offshore waters originated from the North Atlantic Current as indicated by the similarities between the records of IODP Site U1313, ODP Site 980 and MD01-2446, nearshore waters as recorded in core MD03-2699 derived from the Azores Current and thus the subtropical gyre. Except for MIS 12, Azores Current influence seems to be related to eastern boundary system dynamics and not to changes in the Atlantic overturning circulation.

Surface sediment analyses, bivalve and gastropod taxa of marine sediments from the Banc d'Arguin, Mauritania

High-nutrient tropical carbonate systems are known to produce sediments that, in terms of skeletal composition, are reminiscent of their extra-tropical counterparts. Such carbonate systems and associated carbonate grain assemblages in the tropics are rare in the present-day world. Nonetheless, it is crucial to gain a better understanding of those ecosystems, including their drivers and players because such settings potentially represent models for ancient depositional systems as well as for predicted future environmental conditions. One of the modern occurrences of eutrophic tropical carbonate systems is the northern Mauritanian Shelf. The marine environment is characterized by an eastern boundary upwelling system that pushes cool and nutrient-rich intermediate waters onto a wide epicontinental platform (Golfe d'Arguin) where the waters warm up to tropical temperatures. The resulting facies is mixed carbonate-siliciclastic with a dominant foramol association grading into bimol and barnamol grain assemblages in the shallowest areas forming the Banc d'Arguin. Besides this cool water-related heterozoan association, the carbonate sediment is characterized by tropical molluskan species, while chlorozoan biota (e.g., corals and algal symbiont-bearing foraminifers) are entirely absent. We here present a first comprehensive facies analysis of this model example of eutrophic tropical carbonates. Furthermore, we reconstruct the loci of carbonate production and provide a conclusive depositional model of the Banc d'Arguin that received little attention to date due to its poorly accessible nature.

Sea surface temperature record of ODP Site 202-1233

We present the first high-resolution alkenone-derived sea surface temperature (SST) reconstruction in the southeast Pacific (Ocean Drilling Program Site 1233) covering the major part of the last glacial period and the Holocene. The record shows a clear millennial-scale pattern that is very similar to climate fluctuations observed in Antarctic ice cores, suggesting that the Southern Hemisphere high-latitude climate changes extended into the midlatitudes, involving simultaneous changes in air temperatures over Antarctica, sea ice extent, extension of the Antarctic Circumpolar Current, and westerly atmospheric circulation. A comparison to other midlatitude surface ocean records suggests that this "Antarctic" millennial-scale pattern was probably a hemisphere-wide phenomenon. In addition, we performed SST gradient reconstructions over the complete latitudinal range of the Pacific Eastern Boundary Current System for different time intervals during the last 70 kyr. The main results suggest an equatorward displaced subtropical gyre circulation during marine isotope stages 2 and 4.

Stable oxygen isotope record and age determination of sediments from the North Atlantic Ocean

Sea surface temperature and salinity estimates reconstructed using planktonic foraminiferal abundance and delta18O records from core SU90-03 (40°N, 32°W, 2475 m water depth) reveal large climatic fluctuations linked to major instabilities in Northern Hemisphere ice sheets over the last 150 000 years. Episodes of enhanced ice rafted detritus (IRD) input were accompanied by discrete temperature minima, representing coolings of between 4 and 8°C, and reductions in surface salinity of up to 2.5-3.5 per mil. Several additional cooling episodes of a similar magnitude were documented during intervals of low IRD input that appear to be synchronous, within the limits of dating, with ice rafting events spatially confined to higher latitudes. Accelerator mass spectrometer 14C dates for Heinrich events (H1 - 14.2 ka, H2 - 21.4 ka, H3 - 26.7 ka, H4 - 34.8 ka, H5 - 47.2 ka) obtained from core SU90-03 agree well with other published age estimates and suggest a contemporaneous pattern of climate change throughout the North Atlantic during the last glacial period. This interpretation is supported by a comparison of IRD and palaeotemperature records from DSDP site 609 and core SU90-03, which clearly shows that the major climatic fluctuations identified at high latitudes were transmitted toward the subtropics. However, 14C dates suggest that ice rafting episodes may be diachronous to some extent. The northward migration of the polar front after the H1 event at 40°N in the mid-Atlantic occurred at 14 ka, approximately 500 years earlier than along the Portuguese margin, where the southerly advection of polar waters persisted within eastern boundary current system.

Barium in sediments off northwest Africa

High-resolution down-core analyses of the solid phase content of total barium (Batot) and total organic carbon (TOC) back to 25 kyr B.P. were performed on a gravity core from the upper continental slope off Cape Yubi (Morocco). The observed discrepancy between the two potential paleoproductivity proxies, Batot and TOC, initiated supplementary examinations of the pore water, the geochemistry of the clay fraction, X-ray diffraction analyses, and the application of a sequential Ba extraction method of selected samples. Additionally, we analyzed down-core samples of the planktonic foraminifera Turborotalita quinqueloba and Globorotalia inflata for their Ba/Ca ratios. These analyses, which were performed for the first time on these species, were used to reconstruct past oceanic Ba concentrations. We suggest that in the study area, which is characterized by high accumulation rates, the preserved TOC content is a valuable proxy for past primary productivity, whereas the solid phase Batot contents appear to be affected by other mechanisms and factors. Peaks of total barium content in the clay fraction and of Ba/Ca ratios in the planktonic foraminifera shells found during the Younger Dryas and the Heinrich 1 event are likely to result from increased meltwater influx into the northern North Atlantic. We suggest that Ba-enriched meltwater was transmitted by the eastern boundary current system from higher latitudes to the region of the Canary Islands. Total barium contents of the clay fraction (Batot,clay) and Ba/Ca in planktonic foraminifera shells seem to be reliable proxies for this discharge of meltwater.

Stable isotope record of planktonic foraminifera from the western Iberian margin during marine isotope stage 2

The surface water hydrography along the western Iberian margin, as part of the North Atlantic's eastern boundary upwelling system, consists of a complex, seasonally variable system of equatorward and poleward surface and subsurface currents and seasonal upwelling. Not much information exists to ascertain if the modern current and productivity patterns subsisted under glacial climate conditions, such as during marine isotope stage (MIS) 2, and how North Atlantic meltwater events, especially Heinrich events, affected them. To help answer these questions we are combining stable isotope records of surface to subsurface dwelling planktonic foraminifer species with sea surface temperature and export productivity data for four cores distributed along the western and southwestern Iberian margin (MD95-2040, MD95-2041, MD99-2336, and MD99-2339). The records reveals that with the exception of the Heinrich events and Greenland Stadial (GS) 4 hydrographic conditions along the western Iberian margin were not much different from the present. During the Last Glacial Maximum (LGM), subtropical surface and subsurface waters penetrated poleward to at least 40.6°N (site MD95-2040). Export productivity was, in general, high on the western margin during the LGM and low in the central Gulf of Cadiz, in agreement with the modern situation. During the Heinrich events and GS 4, on the other hand, productivity was high in the Gulf of Cadiz and suppressed in the upwelling regions along the western margin where a strong halocline inhibited upwelling. Heinrich event 1 had the strongest impact on the hydrography and productivity off Iberia and was the only period when subarctic surface waters were recorded in the central Gulf of Cadiz. South of Lisbon (39°N), the impact of the other Heinrich events was diminished, and not all of them led to a significant cooling in the surface waters. Thus, climatic impacts of Heinrich events highly varied with latitude and the prevailing hydrographic conditions in this region.

Benthic foraminifera in Southwest Indian Ocean surface sediments

The distribution of deep-sea benthonic foraminifera in core top samples from the southwest Indian Ocean is examined. Principal component analysis reveals two major assemblages. One assemblages between 3600 and 4800-m water depth is dominated by Episominella umbonifera and is associated with cold (Theta = -0.3 to 0.8°C), low salinity (34.66 to 34.72 * 10**-3) Antarctic Bottom Water in the Crozet Basin, in fracture zones, and on the flanks of the Southwest Indian Ridge. A second assemblage, dominated by Planulina wuellerstorfi, Globocassidulina subglobasa, Astrononion echolsi and Pullenia bulloides, is between 1600 and 3800 m on the Crozet Plateau, Madagascar Ridge, Central Indian Ridge, and Southwest Indian Ridge and is associated with relatively warm (Theta = 0.8 to 2.6°C), high salinity (34.72 to 34.76 * 10**-3) North Atlantic Deep Water. The third principal component divides the P. wuellerstorfi assemblage into two subgroups. One is dominated by Epistominella exigua, P. bulloides, P. wuellerstorfi, and A. echolsi and a second is dominated by G. subglobosa. The distribution of the E. umbonifera assemblage and previous hydrographic studies suggest that AABW flows as a western boundary contour current in the Crozet Basin and penetrates fracture zones in the Southwest Indian Ridge between 55 and 57°E and near 66°E as it travels northward into the Madagascar and Mascarene basins. The faunal-water mass associations from the southeast Indian Ocean are compared; the most notable faunal difference is the absence of Uvigerina as a dominant taxon in the southwest Indian Ocean. A comparison of dissolved oxygen and Uvigerina data shows that oxygen is not a major influence upon the distribution of Uvigerina. A correlation analysis of the faunal data and water depth, potential temperature, in situ temperature, salinity, dissolved oxygen, and 1 - Omega, an index of calcium carbonate undersaturation, was carried out to determine the relationships between fauna and hydrography. The second principal component has a significant positive correlation at the 99.9% level with temperature and negative correlations with water depth and 1 - Omega. A general faunal-water mass correlation exists, but it is not possible to determine which variable controls the faunal distributions.

Modeling the sensitivity of precipitation oxygen isotopes to the Last Glacial Maximum boundary conditions: A study using Community Atmosphere Model (CAM3.0)

To understand the validity of d18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (d18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and d18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and d18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in d18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the d18Oprecip distribution among the simulations.

Carbonate and barite content accross the Eocene/Oligocene boundary of ODP Leg 198 sites, Shatsky Rise

The barite and CaCO3 content (in weight percent) of marine sediments can be used to determine spatial and temporal changes in export production (organic and carbonate carbon flux) and/or CaCO3 preservation (inorganic carbon burial). Here we report barite and CaCO3 content in Eocene/Oligocene (E/O) boundary sediments from locations drilled on Shatsky Rise during Ocean Drilling Program Leg 198. Records of these indexes may be used along with other data to determine how the major E/O boundary climatic transition (initiation of Antarctic glaciation and resultant ocean-climate system changes) affected marine export production/preservation at Shatsky Rise. Such data are necessary to elucidate the timing and phasing of changes in the carbon cycle relative to fluctuations in oceanographic conditions across this climatically important interval.

Age determination of sediment cores of the California Current System

Alkenone sea surface temperature (SST) records were generated from the Ocean Drilling Program's (ODP) Sites 1014 and 1016 to examine the response of the California Current System to global climate change during the last 136 ka. The temperature differences between these sites (Delta SST(NEP)=SST(ODP1014)-SST(ODP1016)) reflected the intensity of the California Current and varied between 0.4 and 6.1 °C. A high Delta SST(NEP) (weaker California Current) was found for late marine isotope stage (MIS) 2 and early MIS 5e, while a low Delta SST(NEP) (stronger California Current) was detected for mid-MIS 5e and MIS 1. Spectral analysis indicated that this variation pattern dominated 23- (precession) and 30-ka periods. Comparison of the Delta SST(NEP) and SST based on data from core MD01-2421 at the Japan margin revealed anti-phase variation; the high Delta SST(NEP) (weakening of the California Current) corresponded to the low SST at the Japan margin (the southward displacement of the NW Pacific subarctic boundary), and vice versa. This variation was synchronous with a model prediction of the tropical El Niño-Southern Oscillation behavior. These findings suggest that the intensity of the North Pacific High varied in response to precessional forcing, and also that the response has been linked with the changes of tropical ocean-atmosphere interactions.

Seawater carbonate chemistry and concentration boundary layers around complex assemblages of macroalgae in a laboratory experiment

Metabolic processes have the potential to modulate the effects of ocean acidification (OA) in nearshore macroalgal beds. We investigated whether natural mixed assemblages of the articulate coralline macroalgae Arthrocardia corymbosa and understory crustose coralline algae (CCA) altered pH and O2 concentrations within and immediately above their canopies. In a unidirectional flume, we tested the effect of water velocity (0-0.1 m/s), bulk seawater pH (ambient pH 8.05, and pH 7.65), and irradiance (photosynthetically saturating light and darkness) on pH and O2 concentration gradients, and the derived concentration boundary layer (CBL) thickness. At bulk seawater pH 7.65 and slow velocities (0 and 0.015 m/s), pH at the CCA surface increased to 7.90-8.00 in the light. Although these manipulations were short term, this indicates a potential daytime buffering capacity that could alleviate the effects of OA. Photosynthetic activity also increased O2 concentrations at the surface of the CCA. However, this moderating capacity was flow dependent; the CBL thickness decreased from an average of 26.8 mm from the CCA surface at 0.015 m/s to 4.1 mm at 0.04 m/s. The reverse trends occurred in the dark, with respiration causing pH and O2 concentrations to decrease at the CCA surface. At all flow velocities the CBL thicknesses (up to 68 mm) were much greater than those previously published, indicating that the presence of canopies can alter the CBL substantially. In situ, the height of macroalgal canopies can be an order of magnitude larger than those used here, indicating that the degree of buffering to OA will be context dependent.

Diffusion boundary layers ameliorate the negative effects of ocean acidification on the temperate coralline macroalga Arthrocardia corymbosa

Anthropogenically-modulated reductions in pH, termed ocean acidification, could pose a major threat to the physiological performance, stocks, and biodiversity of calcifiers and may devalue their ecosystem services. Recent debate has focussed on the need to develop approaches to arrest the potential negative impacts of ocean acidification on ecosystems dominated by calcareous organisms. In this study, we demonstrate the role of a discrete (i.e. diffusion) boundary layer (DBL), formed at the surface of some calcifying species under slow flows, in buffering them from the corrosive effects of low pH seawater. The coralline macroalga Arthrocardia corymbosa was grown in a multifactorial experiment with two mean pH levels (8.05 'ambient' and 7.65 a worst case 'ocean acidification' scenario projected for 2100), each with two levels of seawater flow (fast and slow, i.e. DBL thin or thick). Coralline algae grown under slow flows with thick DBLs (i.e., unstirred with regular replenishment of seawater to their surface) maintained net growth and calcification at pH 7.65 whereas those in higher flows with thin DBLs had net dissolution. Growth under ambient seawater pH (8.05) was not significantly different in thin and thick DBL treatments. No other measured diagnostic (recruit sizes and numbers, photosynthetic metrics, %C, %N, %MgCO3) responded to the effects of reduced seawater pH. Thus, flow conditions that promote the formation of thick DBLs, may enhance the subsistence of calcifiers by creating localised hydrodynamic conditions where metabolic activity ameliorates the negative impacts of ocean acidification.

Calcium isotope ratios from marine barite and carbonate over Paleocene-Eocene Boundary (ODP Sites 1212 and 1221)

Carbonates are invaluable archives of the past, and have been used extensively to reconstruct paleoclimate and paleoceanographic conditions over geologic time scales. Such archives are susceptible to diagenetic alteration via dissolution, recrystallization and secondary precipitation, particularly during ocean acidification events when intense dissolution can occur. Despite the importance of diagenesis on proxy fidelity, the effects of diagenesis on the calcium isotopic composition (d44Ca) of carbonates are unclear. Accordingly, bulk carbonate d44Ca was measured at high resolution in two Pacific deep sea sediment cores (ODP Sites 1212 and 1221) with considerably different dissolution histories over the Paleocene-Eocene Thermal Maximum (PETM, ~55 Ma). The d44Ca of marine barite was also measured at the deeper Site 1221, which experienced severe carbonate dissolution during the PETM. Large (~0.8 per mil) variations in bulk carbonate d44Ca occur in the deeper site near the peak carbon isotope excursion, and are correlated with a large drop in carbonate weight percent. Such an effect is seen in neither the 1221 barite record nor the bulk carbonate record at the shallower, less dissolved Site 1212. We contend that ocean chemical changes associated with the abrupt and massive carbon release into the ocean-atmosphere system and subsequent ocean acidification at the PETM affected the bulk carbonate d44Ca record via diagenesis in the sedimentary column. Such changes are considerable, and need to be taken into account when interpreting and modeling Ca isotope data over extreme climatic events associated with ocean chemical evolution.

Planktonic foraminiferal faunas and stable isotope ratios of sediments from the Benguela Upwelling System

Two sediment cores retrieved from the continental slope in the Benguela Upwelling System, GeoB 1706 (19°33.7'S 11°10.5'E) and GeoB 1711 (23°18.9'S, 12°22.6'E), reveal striking variations in planktonic foraminiferal abundances during the last 160,000 years. These fluctuations are investigated to assess changes in the intensity and position of the upwelling centres off Namibia. Four species make up over 95% of the variation within the core, and enable the record to be divided into episodes characterized by particular planktonic foraminiferal assemblages. The fossil assemblages have meaningful ecological significance when compared to those of the modern day and the relationship to their environment. The cold-water planktonic foraminifer, Neogloboquadrina pachyderma sinistral [N. pachyderma (s)], dominates the modern-day, coastal upwelling centres, and Neogloboquadrina pachyderma dextral and Globigerina bulloides characterize the fringes of the upwelling cells. Globorotalia inflata is representative of the offshore boundary between newly upwelled waters and the transitional, reduced nutrient levels of the subtropical waters. In the fossil record, episodes of high N. pachyderma (s) abundances are interpreted as evidence of increased upwelling intensity, and the associated increase in nutrients. The N. pachyderma (s) record suggests temporal shifts in the intensity of upwelling, and corresponding trophic domains, that do not follow the typical glacial-interglacial pattern. Periods of high N. pachyderma (s) abundance describe rapid, discrete events dominating isotope stages 3 and 2. The timing of these events correlates to the temporal shifts of the Angola-Benguela Front (Jansen et al., 1997) situated to the north of the Walvis Ridge. Absence of high abundances of N. pachyderma (s) from the continental slope of the southern Cape Basin indicates that Southern Ocean surface water advection has not exerted a major influence on the Benguela Current System. The coincidence of increased upwelling intensity with the movement of the Angola-Benguela Front can be interpreted mainly by changes in strength and zonality of the trade wind system.

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.