Model simulation output data (CCSM3) related to uplift of African topography and the intensification of upwelling in the Benguela region

The Benguela Current, located off the west coast of southern Africa, is tied to a highly productive upwelling system**1. Over the past 12 million years, the current has cooled, and upwelling has intensified**2, 3, 4. These changes have been variously linked to atmospheric and oceanic changes associated with the glaciation of Antarctica and global cooling**5, the closure of the Central American Seaway**1, 6 or the further restriction of the Indonesian Seaway**3. The upwelling intensification also occurred during a period of substantial uplift of the African continent**7, 8. Here we use a coupled ocean-atmosphere general circulation model to test the effect of African uplift on Benguela upwelling. In our simulations, uplift in the East African Rift system and in southern and southwestern Africa induces an intensification of coastal low-level winds, which leads to increased oceanic upwelling of cool subsurface waters. We compare the effect of African uplift with the simulated impact of the Central American Seaway closure9, Indonesian Throughflow restriction10 and Antarctic glaciation**11, and find that African uplift has at least an equally strong influence as each of the three other factors. We therefore conclude that African uplift was an important factor in driving the cooling and strengthening of the Benguela Current and coastal upwelling during the late Miocene and Pliocene epochs.

Macrobenthic abundance, biomass, productivity and production in the deep Arctic ocean

The few existing studies on macrobenthic communities of the deep Arctic Ocean report low standing stocks, and confirm a gradient with declining biomass from the slopes down to the basins as commonly reported for deep-sea benthos. In this study we have further investigated the relationship of faunal abundance (N), biomass (B) as well as community production (P) with water depth, geographical latitude and sea ice concentration. The underlying dataset combines legacy data from the past 20 years, as well as recent field studies selected according to standardized quality control procedures. Community P/B and production were estimated using the multi-parameter ANN model developed by Brey (2012). We could confirm the previously described negative relationship of water depth and macrofauna standing stock in the Arctic deep-sea. Furthermore, the sea-ice cover increasing with high latitudes, correlated with decreasing abundances of down to < 200 individuals/m**2, biomasses of < 65 mg C/m**2 and P of < 75 mg C/m**2/y. Stations under influence of the seasonal ice zone (SIZ) showed much higher standing stock and P means between 400 - 1400 mg C/m**2/y; even at depths up to 3700 m. We conclude that particle flux is the key factor structuring benthic communities in the deep Arctic ocean, explaining both the low values in the ice-covered Arctic basins and the high values along the SIZ.

Stable isotope record and age determination of ODP Site 121-758 in the northeastern Indian Ocean

Numerous studies have shown that delta18O records from benthic and planktonic foraminifera, primarily a proxy of global ice volume variations, reflect Milankovitch periodicities. To study climatic response to orbital forcing at Ocean Drilling Program site 758, we have generated continuous delta18O and delta13C records from a single benthic foraminiferal species Cibicides wuellerstorfi for the last 3.6 m.y. and extended the planktonic foraminiferal isotope records of Farrell and Janecek (1991, doi:10.2973/odp.proc.sr.121.124.1991) (0-2.5 Ma, based on Globigerinoides sacculifer) to 3.6 Ma (Chen, 1994). We then constructed an age model by matching, correlating and tuning the benthic delta18O record to a model simulation of ice volume (Imbrie and Imbrie, 1980, doi:10.1126/science.207.4434.943). The filtered 41- and 23-kyr signals based on the resultant astronomically tuned age model are highly correlated to obliquity (r=0.83) and precession (r=0.75), respectively. Although derived with methodology different from Shackleton et al. (1990) and Hilgen (1991, doi:10.1016/0012-821X(91)90206-W, 1991, doi:10.1016/0012-821X(91)90082-S), our results generally agree with their published astronomical timescales for the time interval from 0 to 3.0 Ma, providing additional support for the newly emerging chronology based on orbital tuning. Slight discrepancies exist in the time interval from 3.0 to 3.6 Ma, suggesting several possibilities, including differences in the approaches of orbital tuning and the relatively low amplitude of delta18O variations in our record. However, even if the discrepancies are due to the relatively low amplitude of the isotope signals in our record at 3.0-3.6 Ma, our resultant timescale as a whole does not adversely affect our evaluation of the paleoclimatology and paleoceanography of the Indian Ocean, such as the evolution of the 100-, 41- and 23-kyr cycles, and variation of global ice volume and deepwater temperature during the past 3.6 m.y.

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.

Age model, stable isotopes, major element proportions and endmember unmixing results from four western tropical Atlantic cores

We investigate changes in the delivery and oceanic transport of Amazon sediments related to terrestrial climate variations over the last 250 ka. We present high-resolution geochemical records from four marine sediment cores located between 5 and 12° N along the northern South American margin. The Amazon River is the sole source of terrigenous material for sites at 5 and 9° N, while the core at 12° N receives a mixture of Amazon and Orinoco detrital particles. Using an endmember unmixing model, we estimated the relative proportions of Amazon Andean material ("%-Andes", at 5 and 9° N) and of Amazon material ("%-Amazon", at 12° N) within the terrigenous fraction. The %-Andes and %-Amazon records exhibit significant precessional variations over the last 250 ka that are more pronounced during interglacials in comparison to glacial periods. High %-Andes values observed during periods of high austral summer insolation reflect the increased delivery of suspended sediments by Andean tributaries and enhanced Amazonian precipitation, in agreement with western Amazonian speleothem records. Increased Amazonian rainfall reflects the intensification of the South American monsoon in response to enhanced land-ocean thermal gradient and moisture convergence. However, low %-Amazon values obtained at 12° N during the same periods seem to contradict the increased delivery of Amazon sediments. We propose that reorganizations in surface ocean currents modulate the northwestward transport of Amazon material. In agreement with published records, the seasonal North Brazil Current retroflection is intensified (or prolonged in duration) during cold substages of the last 250 ka (which correspond to intervals of high DJF or low JJA insolation) and deflects eastward the Amazon sediment and freshwater plume.

Age model of ODP Site 162-983

At ODP Site 983, relative geomagnetic paleointensity and planktic and benthic delta18O records have been acquired for the last 350 kyr. The mean sedimentation rate in this interval is 11.3 cm/kyr. Magnetic properties and hysteresis ratios indicate that pseudo-single domain magnetite is the remanence carrier. Volume susceptibility (kappa), anhysteretic (ARM) and isothermal (IRM) remanence values vary by a factor of 3-4, well within the criteria usually cited for paleointensity studies. Natural remanent magnetization (NRM) is normalized by ARM and IRM to acquire the paleointensity proxy. Arithmetic means of NRM/ARM and NRM/IRM, calculated for five demagnetization steps in the 25-45 mT range, constitute the relative paleointensity estimates. Some paleointensity lows (particularly those at ~40, ~120 and ~188 ka) are associated with directional excursions of the field, especially the event at ~188 ka (referred to here as the Iceland Basin Event) that constitutes a short-lived polarity reversal. For the last 200 kyr, the records can be correlated with other high-resolution paleointensity records such as those from the Labrador Sea, Mediterranean/Somali Basin and Sulu Sea, implying that the millennial scale features are globally synchronous. A labeling system for paleointensity features is proposed that ties prominent highs and lows to oxygen isotope stages.

Determination of millennial-scale climate variability based on alkenones in the western Atlantic Ocean

High-resolution records of alkenone-derived sea surface temperatures and elemental Ti/Ca ratios from a sediment core retrieved off northeastern Brazil (4° S) reveal short-term climate variability throughout the past 63,000 a. Large pulses of terrigenous sediment discharge, caused by increased precipitation in the Brazilian hinterland, coincide with Heinrich events and the Younger Dryas period. Terrigenous input maxima related to Heinrich events H6-H2 are characterized by rapid cooling of surface water ranging between 0.5 and 2° C. This signature is consistent with a climate model experiment where a reduction of the Atlantic meridional overturning circulation (AMOC) and related North Atlantic cooling causes intensification of NE trade winds and a southward movement of the Intertropical Convergence Zone, resulting in enhanced precipitation off northeastern Brazil. During deglaciation the surface temperature evolution at the core site predominantly followed the Antarctic warming trend, including a cooling, prior to the Younger Dryas period. An abrupt temperature rise preceding the onset of the Bølling/Allerød transition agrees with model experiments suggesting a Southern Hemisphere origin for the abrupt resumption of the AMOC during deglaciation caused by Southern Ocean warming and associated with northward flow anomalies of the South Atlantic western boundary current.

Transient Holocene experiments under orbital forcing using the comprehensive global climate model CCSM3 (Community Climate System Model 3)

The Southern Hemisphere Westerly Winds (SWW) have been suggested to exert a critical influence on global climate through wind-driven upwelling of deep water in the Southern Ocean and the potentially resulting atmospheric CO2 variations. The investigation of the temporal and spatial evolution of the SWW along with forcings and feedbacks remains a significant challenge in climate research. In this study, the evolution of the SWW under orbital forcing from the early Holocene (9 kyr BP) to pre-industrial modern times is examined with transient experiments using the comprehensive coupled global climate model CCSM3. Analyses of the model results suggest that the annual and seasonal mean SWW were subject to an overall strengthening and poleward shifting trend during the course of the early-to-late Holocene under the influence of orbital forcing, except for the austral spring season, where the SWW exhibited an opposite trend of shifting towards the equator.

Isotope ratios of smectites, calcites and whole rocks from DSDP Legs 51 and 52, Atlantic Ocean

Delta180 and 87Sr/86Sr isotopic data from smectites, calcites, and whole rocks, together with published isotopic age determinations, alkali element concentration data and petrographic observations suggest a sequential model of ocean floor alteration. The early stage lasts about 3 m.y. and is characterized by palagonite and smectite formation, and solutions with a large basaltic component, increasing with temperature which varies from 15° to 80° C at DSDP site 418A. Most carbonates are deposited after this stage from solutions with a negligible basaltic Sr component and temperatures of 15° to 40° C. Water of seawater Sr and O isotopic composition is shown to percolate to at least 500 m into the basaltic basement. No evidence was found for continuing exchange of strontium or oxygen after 3 m.y.

Temporal and spatial LGM variability of biomarker records in the Arctic Ocean

Using the sea ice proxy IP25 and phytoplankton-derived biomarkers (brassicasterol and dinosterol) Arctic sea-ice conditions were reconstructed for Marine Isotope Stage (MIS) 3 to 1 in sediment cores from the north of Barents Sea continental margin across the Central Arctic to the Southern Mendeleev Ridge. Our results suggest more extensive sea-ice cover than present-day during MIS 3, increasing sea-ice growth during MIS 2 and decreased sea-ice cover during the last deglacial. The summer ice edge sustained north of the Barents Sea even during extremely cold (i.e., Last Glacial Maximum (LGM)) as well as warm periods (i.e., Bølling-Allerød). During the LGM, the western Svalbard margin and the northern Barents Sea margin areas were characterized by high concentrations of both IP25 and phytoplankton biomarkers, interpreted as a productive ice-edge situation, caused by the inflow of warm Atlantic Water. In contrast, the LGM high Arctic proper (north of 84°N) was covered by thick permanent sea ice throughout the year with rare break up, indicated by zero or near-zero biomarker concentrations. The spring/summer sea-ice margin significantly extended southwards to the southern Lomonosov Ridge and Mendeleev Ridge during the LGM. Our proxy reconstructions are very consistent with published model results based on the North Atlantic/Arctic Ocean Sea Ice Model (NAOSIM).

Benthic foraminiferal stable isotope stratigraphy of ODP Site 138-846 in the tropical Pacific Ocean

A stable-isotope stratigraphy at Site 846 (tropical Pacific, 3°06'S, 90°49'W, 3307 m water depth), based on the benthic foraminifers Cibicides wuellerstorfi and Uvigerina peregrina, yields a high-resolution record of deep-sea delta18O and delta13C over the past 1.8 Ma, with an average sampling interval of 3 k.y. Variance in the delta18O and delta13C records is concentrated in the well-known orbital periods of 100, 41, and 23 k.y. In the 100-k.y. band, both isotopic signals grow from relatively low amplitudes prior to 1.2 Ma, to high amplitudes in the late Quaternary since 0.7 Ma. The amplitude of delta18O and especially of delta13C decreases in the 41-k.y. band as it grows in the 100-k.y. band, consistent with a transfer of energy into an orbitally-paced internal oscillation. A weak 30-k.y. rhythm, present in both delta18O and delta13C, may reflect nonlinear interaction between the 41-k.y. and 100-k.y. bands in the evolving climate system. In the 23-k.y. and 19-k.y. bands associated with orbital precession, delta18O and delta13C are not coherent with each other on long time scales, and do not evolve like the 100-k.y. and 41-k.y. bands. This suggests that the source of the growing 100-k.y. oscillation is not a nonlinear response to precession, in contrast to predictions of some climate models. Sedimentation rates at this site also vary with a strong 100-k.y. cycle. Unlike the isotope records, the amplitude of 100-k.y. variations in sedimentation rate is relatively constant over the past 1.8 Ma, ranging from about 15 to 70 m/m.y. Prior to 0.9 Ma, sedimentation rates co-vary with orbital eccentricity, rather than with global climate as reflected by delta18O or delta13C. A source of this 100-k.y. cycle of sedimentation rate in the absence of similar ice volume fluctuations may be precessional heating of equatorial land masses, which in an energy balance climate model drives variations of monsoonal climates with a 100-k.y. rhythm. For the interval younger than 0.9 Ma, high sedimentation rates in the 100-k.y. band are consistently associated with glacial stages. This change of pattern suggests that when the amplitude of glacial cycles become large enough, their global effects overpower a local monsoon-driven variation in sedimentation rate at Site 846.

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.