Numerical modeling of marine radiocarbon reservoir ages during the last deglaciation

A critical problem in radiocarbon dating is the spatial and temporal variability of marine reservoir ages (MRAs). We assessed the MRA evolution during the last deglaciation by numerical modeling, applying a self-consistent iteration scheme in which an existing radiocarbon chronology (derived by Hughen et al., Quat. Sci. Rev., 25, pp. 3216-3227, 2006) was readjusted by transient, 3-D simulations of marine and atmospheric Delta14C. To estimate the uncertainties regarding the ocean ventilation during the last deglaciation, we considered various ocean overturning scenarios which are based on different climatic background states (PD: modern climate, GS: LGM climate conditions). Minimum and maximum MRAs are included in file 'MRAminmax_21-14kaBP.nc'. Three further files include MRAs according to equilibrium simulations of the preindustrial ocean (file 'C14age_preindustrial.nc'; this is an update of our results published in 2005) and of the glacial ocean (files 'C14age_spinupLGM_GS.nc' and 'C14age_spinupLGM_PD.nc').

Past circulation along the western Iberian margin

Fifteen Iberian margin sediment cores, distributed between 43°12'N and 35°53'N, have been used to reconstruct spatial and temporal (sub)surface circulation along the Iberian margin since the Last Glacial period. Time-slice maps of planktonic foraminiferal derived summer sea surface temperature (SST) and export productivity (Pexp) were established for specific time intervals within the last 35 ky: the Holocene (Recent and last 8 ky), Younger Dryas (YD), Heinrich Stadials (HS) 1, 2a, 2b, 3, and the Last Glacial Maximum (LGM). The SST during the Holocene shows the same latitudinal gradient along the western Iberian margin as present-day with cold but productive areas that reflect the influence of coastal upwelling centers. The LGM appears as a slightly less warm, but more productive period relative to the Holocene and present-day conditions, suggesting that sea-level minima forced a westward displacement of the coastal upwelling centers possibly accompanied by a strengthening of northward winds. During the YD, a longitudinal thermal front is depicted at 10°W, with cold polar waters offshore and warmer subtropical waters inshore, suggesting that the subtropical Paleo-Iberian Poleward Current more likely flowed at a more inshore location masking the local SST signal and amplitude of variation. A substantial cooling and drop in productivity is observed during all HS, in particular HS1 and HS3, reflecting the penetration of icebergs-derived meltwater. These most extreme southward extensions of very cold waters define a strong SST gradient that marks a possible Paleo-Azores Front. Higher production south of this front was likely fed by frontal nutrient advection.

Stable oxygen and carbon isotope record of benthic foraminifera in Pliocene sediments of the Californian margin

New benthic foraminiferal stable isotopic records of northeast Pacific intermediate water (ODP Site 1014, 1177 m) and mid-depth water (ODP Site 1018, 2476 m) were compared to isotopic records of deep water in the tropical Pacific (ODP Site 849, 3851 m) for the reconstruction of vertical profiles of nutrient and physical properties from the Early Pliocene to the Early Pleistocene (approx. 5-1.4 Ma). Our data indicate that, for the entire interval, there was enhanced north Pacific intermediate water ventilation relative to today, and a similar to modern circulation pattern with northward flowing Pacific Bottom Water (PBW) beneath its southward flowing return flow. However, the core of maximally aged return flow resided as deep as ~2500 m (as compared to ~1500 m today), probably due to the strengthened intermediate water flow. Less apparent aging of deep water along its path before 2.7 Ma indicates that thermohaline overturning may have been more rapid in the warm period of the Early Pliocene. In addition, prior to 2.7 Ma, foraminiferal oxygen isotopic values at mid-depth sites are higher than at deep sites (a reversed vertical gradient relative to today) in both the Atlantic and Pacific Oceans. We suggest that NADW was warmer and more saline than today and that it influenced mid-depth waters throughout the Atlantic and Pacific Oceans. Enhanced Pliocene formation of warmer/saltier intermediate water in the north Pacific, and deep water in the north Atlantic, may have been a result of higher than modern high/mid-latitude sea surface temperatures, evaporation, and salinity.

Model results of sensitivity experiments for marine nitrogen cycle in four NetCDF files

The marine nitrogen (N) inventory is thought to be stabilized by negative feedback mechanisms that reduce N inventory excursions relative to the more slowly overturning phosphorus inventory. Using a global biogeochemical ocean circulation model we show that negative feedbacks stabilizing the N inventory cannot persist if a close spatial association of N2 fixation and denitrification occurs. In our idealized model experiments, nitrogen deficient waters, generated by denitrification, stimulate local N2 fixation activity. But, because of stoichiometric constraints, the denitrification of newly fixed nitrogen leads to a net loss of N. This can enhance the N deficit, thereby triggering additional fixation in a vicious cycle, ultimately leading to a runaway N loss. To break this vicious cycle, and allow for stabilizing negative feedbacks to occur, inputs of new N need to be spatially decoupled from denitrification. Our idealized model experiments suggest that factors such as iron limitation or dissolved organic matter cycling can promote such decoupling and allow for negative feedbacks that stabilize the N inventory. Conversely, close spatial co-location of N2 fixation and denitrification could lead to net N loss.

Isopollen maps for 18,000 years B.P. of the offshore of Northwest Africa: Evidence for paleowind circulation

Distribution patterns of the most important pollen types from southern European and northwest African source areas for the 18,000 years B.P. time slice are reconstructed from pollen records of 14 well-dated deep-sea cores located between 37° and 9°N and compared with the modern pollen distribution in this area. It is concluded that the belt with maximum African Easterly Jet transport did not shift latitudinally during the last glacial-interglacial transition but remained at about 20°N. Furthermore, it is substantiated that the trade winds did not shift latitudinally during the last glacial-interglacial transition. This evidence is not compatible with an atmospheric circulation model that assumes a zone of surface westerlies in the northern part of northwest Africa. Trade winds during glacial episodes did, however, intensify, especially from about 36° to 24° N.

On the effects of circulation, sediment resuspension and biological incorporation by diatoms in an ocean model of aluminium, link to model data

The distribution of dissolved aluminium in the West Atlantic Ocean shows a mirror image with that of dissolved silicic acid, hinting at intricate interactions between the ocean cycling of Al and Si. The marine biogeochemistry of Al is of interest because of its potential impact on diatom opal remineralisation, hence Si availability. Furthermore, the dissolved Al concentration at the surface ocean has been used as a tracer for dust input, dust being the most important source of the bio-essential trace element iron to the ocean. Previously, the dissolved concentration of Al was simulated reasonably well with only a dust source, and scavenging by adsorption on settling biogenic debris as the only removal process. Here we explore the impacts of (i) a sediment source of Al in the Northern Hemisphere (especially north of ~ 40° N), (ii) the imposed velocity field, and (iii) biological incorporation of Al on the modelled Al distribution in the ocean. The sediment source clearly improves the model results, and using a different velocity field shows the importance of advection on the simulated Al distribution. Biological incorporation appears to be a potentially important removal process. However, conclusive independent data to constrain the Al / Si incorporation ratio by growing diatoms are missing. Therefore, this study does not provide a definitive answer to the question of the relative importance of Al removal by incorporation compared to removal by adsorptive scavenging.

Nd isotopes of sediments, foraminifera and fish debris of the western North Atlantic

Understanding changes in ocean circulation during the last deglaciation is crucial to unraveling the dynamics of glacial-interglacial and millennial climate shifts. We used neodymium isotope measurements on postdepositional iron-manganese oxide coatings precipitated on planktonic foraminifera to reconstruct changes in the bottom water source of the deep western North Atlantic at the Bermuda Rise. Comparison of our deep water source record with overturning strength proxies shows that both the deep water mass source and the overturning rate shifted rapidly and synchronously during the last deglacial transition. In contrast, any freshwater perturbation caused by Heinrich event 1 could have only affected shallow overturning. These findings show how changes in upper-ocean overturning associated with millennial-scale events differ from those associated with whole-ocean deglacial climate events.

Circulation, eddies, oxygen and nutrient changes in the eastern tropical South Pacific Ocean

A large, subsurface oxygen deficiency zone is located in the eastern tropical South Pacific Ocean (ETSP). The large-scale circulation in the eastern equatorial Pacific and off Peru in November/December 2012 shows the influence of the equatorial current system, the eastern boundary currents, and the northern reaches of the subtropical gyre. In November 2012 the Equatorial Undercurrent is centered at 250 m depth, deeper than in earlier observations. In December 2012 the equatorial water is transported southeastward near the shelf in the Peru-Chile Undercurrent with a mean transport of 1.6 Sv. In the oxygen minimum zone (OMZ) the flow is overlaid with strong eddy activity on the poleward side of the OMZ. Floats with parking depth at 400 m show fast westward flow in the mid-depth equatorial channel and sluggish flow in the OMZ. Floats with oxygen sensors clearly show the passage of eddies with oxygen anomalies. The long-term float observations in the upper ocean lead to a net community production estimate at about 18° S of up to 16.7 mmol C m?3 yr1 extrapolated to an annual rate and 7.7 mmol C m?3 yr?1 for the time period below the mixed layer. Oxygen differences between repeated ship sections are influenced by the Interdecadal Pacific Oscillation, by the phase of El Niño, by seasonal changes, and by eddies and hence have to be interpreted with care. At and south of the equator the decrease in oxygen in the upper ocean since 1976 is related to an increase in nitrate, phosphate, and in part in silicate.

On the role of circulation and mixing in the ventilation of oxygen minimum zones with a focus on the eastern tropical North Atlantic

Ocean observations carried out in the framework of the Collaborative Research Center 754 (SFB 754) "Climate-Biogeochemistry Interactions in the Tropical Ocean" are used to study (1) the structure of tropical oxygen minimum zones (OMZs), (2) the processes that contribute to the oxygen budget, and (3) long-term changes in the oxygen distribution. The OMZ of the eastern tropical North Atlantic (ETNA), located between the well-ventilated subtropical gyre and the equatorial oxygen maximum, is composed of a deep OMZ at about 400 m depth with its core region centred at about 20° W, 10° N and a shallow OMZ at about 100 m depth with lowest oxygen concentrations in proximity to the coastal upwelling region off Mauritania and Senegal. The oxygen budget of the deep OMZ is given by oxygen consumption mainly balanced by the oxygen supply due to meridional eddy fluxes (about 60%) and vertical mixing (about 20%, locally up to 30%). Advection by zonal jets is crucial for the establishment of the equatorial oxygen maximum. In the latitude range of the deep OMZ, it dominates the oxygen supply in the upper 300 to 400 m and generates the intermediate oxygen maximum between deep and shallow OMZs. Water mass ages from transient tracers indicate substantially older water masses in the core of the deep OMZ (about 120-180 years) compared to regions north and south of it. The deoxygenation of the ETNA OMZ during recent decades suggests a substantial imbalance in the oxygen budget: about 10% of the oxygen consumption during that period was not balanced by ventilation. Long-term oxygen observations show variability on interannual, decadal and multidecadal time scales that can partly be attributed to circulation changes. In comparison to the ETNA OMZ the eastern tropical South Pacific OMZ shows a similar structure including an equatorial oxygen maximum driven by zonal advection, but overall much lower oxygen concentrations approaching zero in extended regions. As the shape of the OMZs is set by ocean circulation, the widespread misrepresentation of the intermediate circulation in ocean circulation models substantially contributes to their oxygen bias, which might have significant impacts on predictions of future oxygen levels.

Carbonate preservation reconstructed from ODP Sites 154-927 and 154-929 at Ceará Rise

Enhanced Atlantic overturning during the Pliocene was first proposed almost 10 yrs ago. Evidence for this Pliocene super conveyor scenario has been collected using a number of proxies (e.g., benthic d13C, Nd isotopic composition of manganese crusts). The present study contributes to the existing evidences by using carbonate dissolution and current vigour history of early Pliocene sediments from the Ceará Rise (ODP Sites 927 and 929). In order to reveal carbonate dissolution history, a number of commonly used and newly established proxies were applied, i.e., sand and carbonate contents, foraminifer fragmentation index, Bulloides Dissolution Index and carbonate silt grain-size distributions. Terrigenous silt grain-size distributions were used to unravel variations in relative current strength and sediment input to the two sites. Overall good carbonate preservation at the shallow Site 927 (3314 m water depth) shows that this level was bathed in North Atlantic Deep Water throughout the early Pliocene. The contrastingly poor carbonate preservation record of the deeper Site 929 (4358 m water depth, at present exposed to Antarctic Bottom Water) is frequently interrupted by phases of good carbonate preservation. These results indicate that the depth of the calcite lysocline was mainly tied to present level (ab. 4200 m water depth), and sometimes even dropped to water depths greater than 4360 m due to even more enhanced circulation. Surprisingly the expansion of NADW is not clearly reflected by an increase in current speed as shown by continuously fine terrigenous grain size.

Radiogenic isotope record of Arctic Ocean circulation and weathering inputs of the past 15 million years

Lead (Pb), neodymium (Nd), and strontium (Sr) isotopic analyses were carried out on sediment leachates (reflecting the isotope composition of past seawater) and digests of the bulk residues (reflecting detrital continental inputs) of Integrated Ocean Drilling Program (IODP) Leg 302 and core PS2185 from the Lomonosov Ridge (Arctic Ocean). Our records are interpreted to reflect changes in continental erosion and oceanic circulation, driven predominantly by tectonic forcing on million-year timescales in the older (pre-2 Ma) part of the record and by climatic forcing of weathering and erosion of the Eurasian continental margin on thousand-year timescales in the younger (post-2 Ma) part. These data, covering the past ~15 Ma, show that continental inputs to the central Arctic Ocean have been more closely linked to glacial and hydrological processes occurring on the Eurasian margin than on continental North America and Greenland. The constancy of the detrital input signatures supports the early existence of an Arctic sea ice cover, whereas the major initiation of Northern Hemisphere glaciation at 2.7 Ma appears to have had little impact on the weathering regime of the Eurasian continental margin.

Productivity- and ocean circulation changes in the sub-Antarctic Atlantic during the last deglaciation and Marine Isotope Stage 3

Millennial-scale climate changes during the last glacial period and deglaciation were accompanied by rapid changes in atmospheric CO2 that remain unexplained. While the role of the Southern Ocean as a 'control valve' on ocean-atmosphere CO2 exchange has been emphasized, the exact nature of this role, in particular the relative contributions of physical (for example, ocean dynamics and air-sea gas exchange) versus biological processes (for example, export productivity), remains poorly constrained. Here we combine reconstructions of bottom-water [O2], export production and 14C ventilation ages in the sub-Antarctic Atlantic, and show that atmospheric CO2 pulses during the last glacial- and deglacial periods were consistently accompanied by decreases in the biological export of carbon and increases in deep-ocean ventilation via southern-sourced water masses. These findings demonstrate how the Southern Ocean's 'organic carbon pump' has exerted a tight control on atmospheric CO2, and thus global climate, specifically via a synergy of both physical and biological processes.

Estimation of sea surface temperatures and salinity of ODP Site 177-1089 (Appendix A)

ODP Site 1089 is optimally located in order to monitor the occurrence of maxima in Agulhas heat and salt spillage from the Indian to the Atlantic Ocean. Radiolarian-based paleotemperature transfer functions allowed to reconstruct the climatic history for the last 450 kyr at this location. A warm sea surface temperature anomaly during Marine Isotope Stage (MIS) 10 was recognized and traced to other oceanic records along the surface branch of the global thermohaline (THC) circulation system, and is particularly marked at locations where a strong interaction between oceanic and atmospheric overturning cells and fronts occurs. This anomaly is absent in the Vostok ice core deuterium, and in oceanic records from the Antarctic Zone. However, it is present in the deuterium excess record from the Vostok ice core, interpreted as reflecting the temperature at the moisture source site for the snow precipitated at Vostok Station. As atmospheric models predict a subtropical Indian source for such moisture, this provides the necessary teleconnection between East Antarctica and ODP Site 1089, as the subtropical Indian is also the source area of the Agulhas Current, the main climate agent at our study location. The presence of the MIS 10 anomaly in the delta13C foraminiferal records from the same core supports its connection to oceanic mechanisms, linking stronger Agulhas spillover intensity to increased productivity in the study area. We suggest, in analogy to modern oceanographic observations, this to be a consequence of a shallow nutricline, induced by eddy mixing and baroclinic tide generation, which are in turn connected to the flow geometry, and intensity, of the Agulhas Current as it flows past the Agulhas Bank. We interpret the intensified inflow of Agulhas Current to the South Atlantic as responding to the switch between lower and higher amplitude in the insolation forcing in the Agulhas Current source area. This would result in higher SSTs in the Cape Basin during the glacial MIS 10, due to the release into the South Atlantic of the heat previously accumulating in the subtropical and equatorial Indian and Pacific Ocean. If our explanation for the MIS 10 anomaly in terms of an insolation variability switch is correct, we might expect that a future Agulhas SSST anomaly event will further delay the onset of next glacial age. In fact, the insolation forcing conditions for the Holocene (the current interglacial) are very similar to those present during MIS 11 (the interglacial preceding MIS 10), as both periods are characterized by a low insolation variability for the Agulhas Current source area. Natural climatic variability will force the Earth system in the same direction as the anthropogenic global warming trend, and will thus lead to even warmer than expected global temperatures in the near future.

Physical oceanography during POLARSTERN cruise ARK-XIX/4b

Exchanges between the North Atlantic and the Arctic Ocean result in the most dramatic water mass conversions in the World Ocean: warm and saline Atlantic waters, flowing through the Nordic Seas into the Arctic Ocean, are modified by cooling, freezing and melting to become shallow fresh waters, ice and saline deep waters. The outflow from the Nordic Seas to the south provides the initial driving of the global thermohaline circulation cell. Knowledge of these fluxes and understanding of the modification processes is a major prerequisite for the quantification of the rate of overturning within the large circulation cells of the Arctic and the Atlantic Oceans, and is also a basic requirement for understanding the role of these ocean areas in climate variability on interannual to decadal time scales. The Fram Strait represents the only deep connection between the Arctic Ocean and the Nordic Seas. Just as the freshwater transport from the Arctic Ocean is of major influence on convection in the Nordic Seas and further south, the transport of warm and saline Atlantic water affects the water mass characteristics in the Arctic Ocean which has consequences for the internal circulation and possibly influences also ice and atmosphere. The West Spitsbergen Current carrying Atlantic Water northward. The East Greenland Current, carrying water from the Arctic Ocean southwards has a concentrated core above the continental slope. It is our aim to measure the oceanic fluxes through Fram Strait and to determine their variability in seasonal to decadal time scales. 53 CTD profiles were taken at 51 stations. Two CTD systems from Sea-Bird Electronics Inc SBE911+ were used. Mainly SN 561 with duplicate T and C sensors (temperature sensors SBE3, SN 2685 and 2678, conductivity sensors SBE4, SN 2325 and 2618 and pressure sensor Digiquartz 410K-105 SN 75659) was in service. For the control of the temperature sensors a SBE35 RT digital reversing thermometer, SN 27 was applied. The CTD was connected to a SBE32 Carousel Water Sampler, SN 273 (24 12-liter bottles). For 3 CTD-Stations (726-3, 727-1, 728-1) the Sea-Bird 911+ probe SN 485 was used with temperature sensor SBE3 SN 2460, conductivity sensor SBE4 SN 2054, pressure sensor Digiquartz 410K SN 68997 and the SBE32 Carousel Water Sampler SN 202. Additionally Benthos Altimeters Model 2110-2, SN 189 and SN 208 and Wetlabs C-Star Transmissiometers SN 403 and SN 267 were mounted on the carousels. During the cruise a total number of 184 water samples were analysed with a Guildline Autosal 8400B salinometer, and IAPSO standard seawater batch number P141, K=0.99993. 20 salinity samples were brought back to AWI for onshore analysis. The CTD sensors were calibrated before and after the cruise by Sea-Bird Electronics.

(Table 1) Age determination of sediment cores HU89038_8 and OCE326-GGC5

Models indicate that a complete shutdown of deep and intermediate water production is a possible consequence of extreme climate conditions in the northern North Atlantic, and the high ratio of 231Pa to 230Th on Bermuda Rise is evidence that this might have happened ?17 ka during Heinrich event 1 (H1). However, new radiocarbon data from bivalves that lived at ?4.6 km on the Bermuda Rise during H1 lead to a different conclusion. The bivalve data do indeed indicate ventilation of the deep western North Atlantic was suppressed during H1 but not as much as it was during the last glacial maximum. We propose that high diatom flux to the Bermuda Rise during H1 is at least in part responsible for increased 231Pa/230Th at that time. Although we cannot say for sure why opal production was so high in a gyre center location at that time, increased leakage of silica rich waters from the Southern Ocean to the North Atlantic is one possibility.