Relative abundance and ranges of selected diatoms from Pliocene-Pleistocene sections of ODP Leg 177, Atlantic sector of the Southern Ocean

During Ocean Drilling Program (ODP) Leg 177, seven sites were drilled aligned on a transect across the Antarctic Circumpolar Current in the Atlantic sector of the Southern Ocean. The primary scientific objective of Leg 177 was the study of the Cenozoic paleoceanographic and paleoclimatic history of the southern high latitudes and its relationship with the Antarctic cryosphere development. Of special emphasis was the recovery of Pliocene-Pleistocene sections, allowing paleoceanographic studies at millennial or higher time resolution, and the establishment of refined biostratigraphic zonations tied to the geomagnetic polarity record and stable isotope records. At most sites, multiple holes were drilled to ensure complete recovery of the section. A description of the recovered sections and the construction of a multihole splice for the establishment of a continuous composite is presented in the Leg 177 Initial Reports volume for each of the sites (Gersonde, Hodell, Blum, et al., 1999). Here we present the relative abundance pattern and the stratigraphic ranges of diatom taxa encountered from shore-based light microscope studies completed on the Pliocene-Pleistocene sequences from six of the drilled sites (Sites 1089-1094). No shore-based diatom studies have been conducted on the Pliocene-Pleistocene sediments obtained at Site 1088, located on the northern crest of the Agulhas Ridge, because of the scattered occurrence and poor preservation of diatoms in these sections (Shipboard Scientific Party, 1999b). The data included in our report present the baseline of a diatom biostratigraphic study of Zielinski and Gersonde (2002), which (1) includes a refinement of the southern high-latitude Pliocene-Pleistocene diatom zonation, in particular for the middle and late Pleistocene, and (2) presents a biostratigraphic framework for the establishment of age models of the recovered sediment sections. Zielinski and Gersonde (2002) correlated the diatom ranges with the geomagnetic polarity record established shipboard (Sites 1090 and 1092) (Shipboard Scientific Party, 1999c, 1999d) and on shore (Sites 1089, 1091, 1093, and 1094) by Channell and Stoner (2002). The Pliocene-Pleistocene diatom zonation proposed by Zielinski and Gersonde (2002) relies on a diatom zonation from Gersonde and Bárcena (1998) for the northern belt of the Southern Ocean. Because of latitudinal differentiation of sea-surface temperature, nutrients, and salinity between Antarctic and Subantarctic/subtropical water masses, the Pliocene-Pleistocene stratigraphic marker diatoms are not uniformly distributed in the Southern Ocean (Fenner, 1991; Gersonde and Bárcena, 1998). As a consequence, Zielinski and Gersonde (2002) propose two diatom zonations for application in the Antarctic Zone south of the Polar Front (Southern Zonation, Sites 1094 and 1093) and the area encompassing the Polar Front Zone (PFZ) and the Subantarctic Zone (Northern Zonation, Sites 1089-1092). This accounts especially for the Pleistocene zonation where Hemidiscus karstenii, whose first abundant occurrence datum and last occurrence datum defines the subzonation of the northern Thalassiosira lentiginosa Zone, occurs only sporadically in the cold-water realm south of the PFZ and thus is not applicable in sections from this area. However, newly established marker species assigned to the genus Rouxia (Rouxia leventerae and Rouxia constricta) are more related to cold-water environments and allow a refinement of the Pleistocene stratigraphic zonation for the southern cold areas. A study relying on quantitative counts of both Rouxia species confirms the utility of these stratigraphic markers for the identification of sequences attributed to marine isotope Stages 6 and 8 in the southern Southern Ocean (Zielinski et al., 2002).

Age models, color measurements, and calculated sea surface temperatures of sediment cores from the Bering Sea and the Sea of Okhotsk

Past changes in North Pacific sea surface temperatures and sea-ice conditions are proposed to play a crucial role in deglacial climate development and ocean circulation but are less well known than from the North Atlantic. Here, we present new alkenone-based sea surface temperature records from the subarctic northwest Pacific and its marginal seas (Bering Sea and Sea of Okhotsk) for the time interval of the last 15 kyr, indicating millennial-scale sea surface temperature fluctuations similar to short-term deglacial climate oscillations known from Greenland ice-core records. Past changes in sea-ice distribution are derived from relative percentage of specific diatom groups and qualitative assessment of the IP25 biomarker related to sea-ice diatoms. The deglacial variability in sea-ice extent matches the sea surface temperature fluctuations. These fluctuations suggest a linkage to deglacial variations in Atlantic meridional overturning circulation and a close atmospheric coupling between the North Pacific and North Atlantic. During the Holocene the subarctic North Pacific is marked by complex sea surface temperature trends, which do not support the hypothesis of a Holocene seesaw in temperature development between the North Atlantic and the North Pacific.

Age models of sediment cores from the Southern Ocean

Recent geochemical models invoke ocean alkalinity changes, particularly in the surface Southern Ocean, to explain glacial age pCO2 reduction. In such models, alkalinity increases in glacial periods are driven by reductions in North Atlantic Deep Water (NADW) supply, which lead to increases in deep-water nutrients and dissolution of carbonate sediments, and to increased alkalinity of Circumpolar Deep Water upwelling in the surface Southern Ocean. We use cores from the Southeast Indian Ridge and from the deep Cape Basin in the South Atlantic to show that carbonate dissolution was enhanced during glacial stages in areas now bathed by Circumpolar Deep Water. This suggests that deep Southern Ocean carbonate ion concentrations were lower in glacial stages than in interglacials, rather than higher as suggested by the polar alkalinity model [Broecker and Peng, 1989, doi:10.1029/GB001i001p00015]. Our observations show that changes in Southern Ocean CaCO3 preservation are coherent with changes in the relative flux of NADW, suggesting that Southern Ocean carbonate chemistry is closely linked to changes in deepwater circulation. The pattern of enhanced dissolution in glacials is consistent with a reduction in the supply of nutrient-depleted water (NADW) to the Southern Ocean and with an increase of nutrients in deep water masses. Carbonate mass accumulation rates on the Southeast Indian Ridge (3200-3800 m), and in relatively shallow cores (<3000 m) from the Kerguelen Plateau and the South Pacific were significantly reduced during glacial stages, by about 50%. The reduced carbonate mass accumulation rates and enhanced dissolution during glacials may be partly due to decreases in CaCO3:Corg flux ratios, acting as another mechanism which would raise the alkalinity of Southern Ocean surface waters. The polar alkalinity model assumes that the ratio of organic carbon to carbonate production on surface alkalinity is constant. Even if overall productivity in the Southern Ocean were held constant, a decrease in the CaCO3:Corg ratio would result in increased alkalinity and reduced pCO2 in Southern Ocean surface waters during glacials. This ecologically driven surface alkalinity change may enhance deepwater-mediated changes in alkalinity, and amplify rapid changes in pCO2.

Age models of ODP Sites 117-722 and 117-728

Pliocene and Pleistocene sediments of the Oman margin and Owen Ridge are characterized by continuous alternation of light and dark layers of nannofossil ooze and marly nannofossil ooze and cyclic variation of wet-bulk density. Origin of the wet-bulk density and color cycles was examined at Ocean Drilling Program Site 722 on the Owen Ridge and Site 728 on the Oman margin using 3.4-m.y.-long GRAPE (gamma ray attenuation) wet-bulk density records and records of sediment color represented as changes in gray level on black-and-white core photographs. At Sites 722 and 728 sediments display a weak correlation of decreasing wet-bulk density with increasing darkness of sediment color. Wet-bulk density is inversely related to organic carbon concentration and displays little relation to calcium carbonate concentration, which varies inversely with the abundance of terrigenous sediment components. Sediment color darkens with increasing terrigenous sediment abundance (decreasing carbonate content) and with increasing organic carbon concentration. Upper Pleistocene sediments at Site 722 display a regular pattern of dark colored intervals coinciding with glacial periods, whereas at Site 728 the pattern of color variation is more irregular. There is not a consistent relationship between the dark intervals and their relative wet-bulk density in the upper Pleistocene sections at Sites 722 and 728, suggesting that dominance of organic matter or terrigenous sediment as primary coloring agents varies. Spectra of wet-bulk density and optical density time series display concentration of variance at orbital periodicities of 100, 41, 23, and 19 k.y. A strong 41-k.y. periodicity characterizes wet-bulk density and optical density variation at both sites throughout most of the past 3.4 m.y. Cyclicity at the 41-k.y. periodicity is characterized by a lack of coherence between wet-bulk density and optical density suggesting that the bulk density and color cycles reflect the mixed influence of varying abundance of terrigenous sediments and organic matter. The 23-k.y. periodicity in wet-bulk density and sediment color cycles is generally characterized by significant coherence between wet-bulk density and optical density, which reflects an inverse relationship between these parameters. Varying organic matter abundance, associated with changes in productivity or preservation, is inferred to more strongly influence changes in wet-bulk density and sediment color at this periodicity.

Age models of upper Pliocene samples of the North Atlantic

Abundance variations of six Pliocene species of discoasters have been analyzed over the time interval from 1.89 to 2.95 Ma at five contrasting sites in the North Atlantic: Deep Sea Drilling Project Sites 552 (56°N) and 607 (41°N) and Ocean Drilling Program 658 (20°N), 659 (18°N), and 662 (1°S). A sampling interval equivalent to approximately 3 k.y. was used. Total Discoaster abundance showed a reduction with increasing latitude and from the effects of upwelling. This phenomenon is most obvious in Discoaster brouweri, the only species that survived over the entire time interval studied. Prior to 2.38 Ma, Discoaster pentaradiatus and Discoaster surculus are important components of the Discoaster assemblage: Discoaster pentaradiatus increases slightly with latitude up to 41°N, and its abundance relative to D. brouweri increases up to 56°N; D. surculus increases in abundance with latitude and with upwelling conditions relative to both D. brouweri and D. pentaradiatus and is dominant to the latter species at upwelling Site 658 and at the highest latitude sites. Discoaster asymmetricus and Discoaster tamalis appear to increase in abundance with latitude relative to D. brouweri. Many of the abundance changes observed appear to be connected with the initiation of glaciation in the North Atlantic at 2.4 Ma. The long-term trend of decreasing Discoaster abundance probably reflects the fall of sea-surface temperatures. This trend of cooling is overprinted by short-term variations that are probably associated with orbital forcing. Evidence for the astronomical elements of eccentricity and obliquity periodicities were found at all sites; however, only at Sites 607, 659, and 662 were precessional periodicities detected. Furthermore, only at Site 659 was precession found to be dominant to obliquity. Abundance peaks of individual species were found to cross-correlate between sites. The distinct abundance fluctuations observed especially in the tropics suggest that temperature is not the only factor responsible for this variation. This study reveals for the first time the importance of productivity pressure on the suppression of Discoaster abundance.

Helium isotope ratios and sedimentation rate models of ODP holes

A continuous age model for the brief climate excursion at the Paleocene-Eocene boundary has been constructed by assuming a constant flux of extraterrestrial 3He (3He[ET]) to the seafloor. 3He[ET] measurements from ODP Site 690 provide quantitative evidence for the rapid onset (models indicating extremely rapid release of isotopically light carbon, possibly from seafloor methane hydrate, as the proximal cause of the event. However, the 3He[ET] technique indicates a previously unrecognized and extreme increase in sedimentation rate coincident with the return of climate proxies to pre-event values. The 3He[ET]-based age model thus suggests a far more rapid recovery from the climatic perturbation than previously proposed or predicted on the basis of the modern carbon cycle, and so may indicate additional or accelerated mechanisms of carbon removal from the ocean-atmosphere system during this period. 3He[ET] was also measured at ODP Site 1051 to test the validity of the Site 690 chronology. Comparison of these data sets seems to require removal of several tens of kyr of sediment within the climatic excursion at Site 1051, an observation consistent with sediment structures and previous age modeling efforts. The Site 1051 age model shows a ~30 kyr period in which climate proxies return toward pre-event values, after which they remain invariant for ~80 kyr. If this rise represents the recovery interval identified at Site 690, then the 3HeET-based age models of the two sites are in good agreement. However, alternative interpretations are possible, and work on less disrupted sites is required to evaluate the reliability of the proposed new chronology of the climate excursion. Regardless of these details, this work shows that the 3HeET technique can provide useful independent evidence for the development and testing of astronomically calibrated age models.

Age models of sediment cores from the continental margin of northwest Africa

Diatom abundance and species composition were quantitatively studied in two latest Quaternary (~130 ka to the Present) sequences from the continental margin of northwest Africa. Off this region, coastal upwelling is well developed under the influence of the NE trade winds. Variations in diatom abundance in these cores are inferred to represent changes caused by varying degrees of the upwelling fertility. Times of high productivity are marked by high relative frequencies of Chaetoceros, while low productivity is marked by the dominance of Aulacoseira granulata. Upwelling increased during glacial episodes (isotopic stages 2-4 and 6) relative to isotopic stages 1 and 5. During the late Holocene, primary productivity levels are similar to those for Stage 5, but in the early Holocene upwelling intensities seem to have been weaker than today. The paleoproductivity reconstruction based on the diatom record is supported by paleoproductivity estimations based on the organic carbon content of the sediments (Sarnthein et al., 1987).

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

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

3D models created from ROV videos and their corresponding subtransect dimensions

Underwater video transects have become a common tool for quantitative analysis of the seafloor. However a major difficulty remains in the accurate determination of the area surveyed as underwater navigation can be unreliable and image scaling does not always compensate for distortions due to perspective and topography. Depending on the camera set-up and available instruments, different methods of surface measurement are applied, which make it difficult to compare data obtained by different vehicles. 3-D modelling of the seafloor based on 2-D video data and a reference scale can be used to compute subtransect dimensions. Focussing on the length of the subtransect, the data obtained from 3-D models created with the software PhotoModeler Scanner are compared with those determined from underwater acoustic positioning (ultra short baseline, USBL) and bottom tracking (Doppler velocity log, DVL). 3-D model building and scaling was successfully conducted on all three tested set-ups and the distortion of the reference scales due to substrate roughness was identified as the main source of imprecision. Acoustic positioning was generally inaccurate and bottom tracking unreliable on rough terrain. Subtransect lengths assessed with PhotoModeler were on average 20% longer than those derived from acoustic positioning due to the higher spatial resolution and the inclusion of slope. On a high relief wall bottom tracking and 3-D modelling yielded similar results. At present, 3-D modelling is the most powerful, albeit the most time-consuming, method for accurate determination of video subtransect dimensions.

Age models and sedimentation rates of sites in the tropical Atlantic

Late Pleistocene signals of calcium carbonate, organic carbon, and opaline silica concentration and accumulation are documented in a series of cores from a zonal/meridional/depth transect in the equatorial Atlantic Ocean to reconstruct the regional sedimentary history. Spectral analysis reveals that maxima and minima in biogenous sedimentation occur with glacial-interglacial cyclicity as a function of both (1) primary production at the sea surface modulated by orbitally forced variation in trade wind zonality and (2) destruction at the seafloor by variation in the chemical character of advected intermediate and deep water from high latitudes modulated by high-latitude ice volume. From these results a pattern emerges in which the relative proportion of signal variance from the productivity signal centered on the precessional (23 kyr) band decreases while that of the destruction signal centered on the obliquity (41 kyr) and eccentricity (100 kyr) periods increases below ~3600-m ocean depth.

Age models and stable isotope record of sediment cores from the Atlantic Ocean

The high-resolution delta18O and delta13C records of benthic foraminifera from a 150,000-year long core from the Caribbean Sea indicate that there was generally high delta13C during glaciations and low delta13C during interglaciations. Due to its 1800-m sill depth, the properties of deep water in the Caribbean Sea are similar to those of middepth tropical Atlantic water. During interglaciations, the water filling the deep Caribbean Sea is an admixture of low delta13C Upper Circumpolar Water (UCPW) and high delta13C Upper North Atlantic Deep Water (UNADW). By contrast, only high delta13C UNADW enters during glaciations. Deep ocean circulation changes can influence atmospheric CO2 levels (Broecker and Takahashi, 1985; Boyle, 1988 doi:10.1029/JC093iC12p15701; Keir, 1988 doi:10.1029/PA003i004p00413; Broecker and Peng, 1989 doi:10.1029/GB003i003p00215). By comparing delta13C records of benthic foraminifera from cores lying in Southern Ocean Water, the Caribbean Sea, and at several other Atlantic Ocean sites, the thermohaline state of the Atlantic Ocean (how close it was to a full glacial or full interglacial configuration) is characterized. A continuum of circulation patterns between the glacial and interglacial extremes appears to have existed in the past. Subtracting the deep Pacific (~mean ocean water) delta13C record from the Caribbean delta13C record yields a record which describes large changes in the Atlantic Ocean thermohaline circulation. The delta13C difference varies as the vertical nutrient distribution changes. This new proxy record bears a striking resemblance to the 150,000-year-long atmospheric CO2 record (Barnola et al., 1987 doi:10.1038/329408a0). This favorable comparison between the new proxy record and the atmospheric CO2 record is consistent with Boyle's (1988a) model that vertical nutrient redistribution has driven large atmospheric CO2 changes in the past. Changes in the relative contribution of NADW and Pacific outflow water to the Southern Ocean are also consistent with Broecker and Peng's (1989) recent model for atmospheric CO2 changes.