Paleoceanography of the eastern equatorial Pacific during the Neogene: Synthesis of Leg 138 drilling results

The primary objective of Leg 138 was to provide detailed information about the ocean\u27s response to global climate change during the Neogene. Two north-south transects were drilled (95° and 110°W) within the region of equatorial divergence-driven upwelling (and thus high accumulation rates and resolution) and spanning the major equatorial ocean current boundaries (and thus recording a high-amplitude signal of the response of the sediment to climatically and/or tectonically driven changes in ocean circulation). The Neogene is marked by a number of well-known climatic and tectonic events (the closing of the Isthmus of Panama, the onset of North Atlantic Deep Water (NADW), the rapid uplift of the Himalayas, the major intensification of Northern Hemisphere glaciation), and the response of the ocean before and after these events was a key focus of Leg 138 drilling. To address these objectives at the highest resolution possible, the Leg 138 scientific staff developed a number of new shipboard strategies and analytical procedures. These included the real-time analysis of the near-continuous gamma ray attenuation porosity evaluator (GRAPE) and susceptibility profiles produced by the multisensor track (MST) on unsplit cores to monitor core recovery and, if necessary, to modify the drilling strategy to ensure proper offset of coring gaps; the collection of near-continuous color reflectance data on split cores; the logging of the first hole drilled at each site to optimize drilling and sampling strategies for subsequent holes; and the use of multiple continuous records to unambiguously construct complete composite sections for each site. The complete, continuous records provided by the GRAPE (with a temporal resolution of often yr), in conjunction with an excellent microfossil stratigraphy and often excellent magnetostratigraphy, allowed for astronomical tuning of the stratigraphic record and resulted in a set of internally consistent, high-resolution age models that provide a secure, absolute time scale for the past 6 m.y. For the period before 6 m.y., the absolute time calibration is less secure, but it is still better than any previously offered. The high-resolution stratigraphic framework of Leg 138 provided new insight into the previously ambiguous tectonic history of the region. By assuming that maximum sedimentation rates along the north-south transect would be expected at the equator, the Leg 138 stratigraphy supports the 1985 work of Cox and Engerbretson, which calls for two different poles of rotation of the Pacific Plate during the interval 0-20 Ma. The Leg 138 plate reconstructions also support several previously hypothesized ridge crest jumps and a slowing of the absolute motion of the Nazca Plate at about 5 Ma. Although Leg 138 data that predates about 13 Ma is limited, the impression that one can gain from these data is that the eastern equatorial Pacific was characterized by relatively high carbonate concentrations and accumulation rates before about 11 Ma. This pattern was interrupted occasionally by rapid massive outpourings of near-monospecific laminated diatom oozes that probably represent the formation of massive mats along strong surface-water fronts. The laminated diatom oozes (LDO) continue to be present in the Leg 138 record (many of them being expressed as seismic reflections) until about 4.4 Ma. Carbonate accumulation rates begin to decline slowly between 11 and 9.8 Ma, when, at about 9.5 Ma, a near-complete loss of carbonate (the carbonate crash ) takes place everywhere in the Leg 138 region (and beyond), except at the westernmost sites close to the equator. The carbonate crash was a time of fundamental change for the eastern equatorial Pacific, and perhaps for most of the ocean basins. Unlike many of the carbonate variations that precede and postdate it, this crash represents a major dissolution event whose effects can be traced seismically in the central and western Pacific. The changes in bottom-water chemistry associated with this event (or series of events) appear to be related to the early phases of the closing of the Panama Gateway. The role of NADW initiation and intensification for controlling carbonate accumulation in the eastern equatorial Pacific is still not resolved; however, ocean modeling demonstrates that the closing of the Panama Gateway may also have a direct influence on NADW production. Therefore, the effects of changes in the Panama Gateway sill depth and the production of NADW may be manifested in the history of eastern equatorial Pacific sedimentation. The carbonate crash was followed by a recovery of the carbonate system (except in the Guatemala and Peru basins, which never recovered) that led up to the late Miocene/early Pliocene sedimentation rate maxima, during which equatorial sedimentation rates are as much as five times greater than those of the late Pliocene or Pleistocene. Examination of modern productivity/preser vation relationships implies that the sedimentation rate maximum was the result of enhanced productivity. The distribution of eolian sediments and isotopic gradients, along with an analysis of the modes of variance in carbonate deposition over the last 6 m.y., suggest a more northerly position of the Intertropical Convergence Zone (ITCZ), a stronger north-south gradient across the equator, and a more zonal circulation focused along the equator during the time of maximum sedimentation. The mechanisms suggested for these changes in circulation patterns include the response of the eastern equatorial Pacific to the closing of the Isthmus of Panama, as well as a global increase in the flux of Ca and Si into the oceans, a possible response to evolution of the Himalayas and the Tibetan Plateau. In an effort to understand the response of the climate system to external (orbital) forcing, 6-m.y.-long, continuous records of carbonate (derived from GRAPE), δ 1 8 and insolation were analyzed and compared. Evolutionary spectral calculations of the variance and coherence among these records indicate that the insolation record is dominated by precessional frequencies, but that the relative importance of the two precessional frequencies has changed significantly over the last 6 m.y. In general, precessional forcing is not found in the carbonate or isotopic records. In the tilt band, however, a linear response is present between solar forcing and the carbonate and isotope records over some intervals. The carbonate record appears to be tightly coupled to the tilt component of insolation before about 1.9 Ma; however, the isotope record does not begin to show sensitivity to orbital tilt until about 4.5 Ma, the time of significant changes in sedimentation patterns in the eastern equatorial Pacific. Only during the last 500,000 yr do all frequencies respond in a similar manner; we also see a marked increase in the response of the isotopic record to orbital forcing (including 100,000- and 400,000-yr periods)

Model of late Pleistocene-Holocene variations in rate of sediment accumulation : Panama Basin, eastern equatorial Pacific

The assumption of constant quartz accumulation for the deep-sea sediment core Y69-106P, taken in the Panama Basin, has been used to date the core and construct a sedimentation rate versus time curve for it. Stratigraphic control for the calculated time scale includes three carbon-14 measurements, the extinction of the radiolarian Stylatractus universus, and correlation with oxygen isotope curves from other dated cores. The model sedimentation rates and sediment composition data allow the determination of sediment accumulation rates for calcium carbonate, opaline silica, and remaining "detritus". Fluctuations in calcium carbonate accumulation rates correlate closely with variations in oxygen isotope ratios in biogenous carbonate from two other equatorial Pacific sediment cores. High oxygen isotope ratios correlate with high calcium carbonate accumulation rates. Opaline silica accumulation rates reflect changes in the dominance of radiolarian fossil groups which can be related to surface circulation in the Panama Basin. The successive maxima in opal and calcium accumulation rates and the oxygen isotope ratios in Y69-106P suggest that at the end of the last glacial period, the rate of calcium carbonate dissolution increased in the Panama Basin first, to be followed by a decrease in the intensity of the eastern equatorial Pacific upwelling, and finally the volume of continental glaciers decreased. Spectral analysis of the rate of accumulation of calcium carbonate in Y69-106P and of the oxygen isotope record of core V28-238 indicate the presence of a 23,000 year periodicity. Analysis of opal accumulation rate suggest a 100,000 year periodicity. These two periodicities are found in fluctuations in the earth's orbital parameters. The eccentricity of the orbit fluctuates with an average period of 98,000 years and the precession of the equinoxes changes with a 22,000 year period

1995, Spatial and temporal variability of late Neogene equatorial Pacific carbonate

High-resolution, continuous records of GRAPE wet bulk density (a carbonate proxy) from Ocean Drilling Program Leg 138 provide one the opportunity for a detailed study of eastern equatorial Pacific Ocean carbonate sedimentation during the last 6 m.y. The transect of sites drilled spans both latitude and longitude in the eastern equatorial Pacific from 90° to 110°W and from 5°S to 10°N. Two modes of variability are resolved through the use of Empirical Orthogonal Function (EOF) analysis. In the presence of large tectonic and climatic boundary condition changes over the last 6 m.y., the dominant mode of spatial variability in carbonate sedimentation is remarkably constant. The first mode accounts for over 50% of the variance in the data, and is consistent with forcing by equatorial divergence. This mode characterizes both carbonate concentration and carbonate mass accumulation rate time series. Variability in the first mode is highly coherent with insolation, indicating a strong linear relationship between equatorial Pacific car bonate sedimentation and Milankovitch variability. Frequency domain analysis indicates that the coupling to equatorial divergence in carbonate sedimentation is strongest in the precession band (19-23 k.y.) and weakest though present at lower frequencies. The second mode of variability has a consistent spatial pattern of east-west asymmetry over the past 4 m.y. only; prior to 4 Ma, a different mode of spatial variability may have been present, possibly suggesting influence by closure of the Isthmus of Panama or other tectonic changes. The second mode of variability may indicate influence by CaCO3 dissolution. The second mode of variability is not highly coherent with insolation. Comparison of the modes of carbonate variability to a 4 m.y. record of benthic δ 1 8 indicates that although overall correlation between carbonate and δ 1 8 is low, both modes of variability in carbonate sedimentation are coherent with δ 1 8 changes at some frequencies. The first mode of carbonate variability is coherent with Sites 846/849 δ 1 8 at the dominant insolation periods, and the second mode is coherent at 100 k.y. during the last 2 m.y. The coherence between carbonate sedimentation and δ 1 8 in both EOF modes suggests that multiple uncorrelated modes of variability operated within the climate system during the late Neogene

Downhole Logging as a Paeoceanographic Tool on Ocean Drilling Program Leg 138: Interface Between High-Resolution Stratigraphy and Regional Syntheses

On Ocean Drilling Program (ODP) Leg 138, standard shipboard procedures were modified to allow for the real-time monitoring of several laboratory core-scanning systems that provide centimeter-scale measurements of saturated bulk density, magnetic susceptibility and digital color reflectance. These continuous, high-resolution data sets were used to ensure the proper offset of multiple holes and to splice together complete sedimentary sections. Typically, the spliced, continuous sediment sections were found to be about 10% longer than the section drilled, as measured by the length of the drill string. While the source of this elongation is not yet fully understood, it must be compensated for in order to property determine sediment fluxes and mass accumulation rates. Downhole logging, in conjunction with inverse correlation techniques provided a means to determine where the distortion occurred and to correct back to true in situ depths. Downhole logging also provides a means, through the generation of synthetic seismograms, of precisely relating the paleoceanographic events found in the core record to the high-resolution seismic record. Once correlated to the seismic record, the spatial and temporal extent of paleoceanographic events can be traced well beyond the borehole. Most seismic events in the equatorial Pacific are related to rapid changes in carbonate content that, in turn, are related to both productivity events (often expressed as monospecific laminated diatom oozes) and times of enhanced dissolution. While many of these events may have oceanwide extent, others, like the absence of carbonate in the late-Miocene to Recent in the Guatemala Basin have been shown to be regional and confined to only the deeper portions of the Guatemala Basin. As we identify and trace specific paleoceanographic events in the seismic record, we can begin to explore the response of the ocean through gradients of latitude, productivity, and depth

Downhole logging as a paeoceanographic tool on ocean drilling program leg 138: Interface between high-resolution stratigraphy and regional syntheses

On Ocean Drilling Program (ODP) Leg 138, standard shipboard procedures were modified to allow for the real-time monitoring of several laboratory core-scanning systems that provide centimeter-scale measurements of saturated bulk density, magnetic susceptibility and digital color reflectance. These continuous, high-resolution data sets were used to ensure the proper offset of multiple holes and to splice together complete sedimentary sections. Typically, the spliced, continuousediment sections were found to be about 10% longer than the section drilled, as measured by the length of the drill string. While the source of this elongation is not yet fully understood, it must be compensated for in order to property determine sediment fluxes and mass accumulation rates. Downhole logging, in conjunction with inverse correlation techniques provided a means to determine where the distortion occurred and to correct back to true in sire depths. Downhole logging also provides a means, through the generation of synthetic seismograms, of precisely relating the paleoceanographic events found in the core record to the high-resolution seismic record. Once correlated to the seismic record, the spatial and temporal extent of paleoceanographic events can be traced well beyond the borehole. Most seismic events in the equatorial Pacific are related to rapid changes in carbonate contenthat, in turn, are related to both productivity events (often expressed as monospecific laminated diatom oozes) and times of enhanced dissolution. While many of these events may have oceanwide extent, others, like the absence of carbonate in the late-Miocene to Recent in the Guatemala Basin have been shown to be regional and confined to only the deeper portions of the Guatemala Basin. As we identify and trace specific paleoceanographic events in the seismic record, we can begin to explore the response of the ocean through gradients of latitude, productivity, and depth

Patterns of CaCO₃deposition in the eastern tropical Pacific Ocean for the last 150 kyr: Evidence for a southeast Pacific depositional spike during marine isotope stage (MIS) 2

We constructed biogenic mass accumulation rate (MAR) time series for eastern Pacific core transects across the equator at ~105˚and ~85˚W and along the equator from 80˚to 140˚W. We used empirical orthogonal function (EOF) analysis to extract spatially coherent patterns of CaCO₃deposition for the last 150 kyr. EOF mode 1 (51% variance) is a CaCO₃ MAR spike centered in marine oxygen isotope stage 2 (MIS 2) found under the South Equatorial Current. EOF mode 2 (19% of variance) is high north of the equator. EOF mode 3 (9% of variance) is an east-west mode centered along the North Equatorial Counter Current. The MIS 2 CaCO₃spike is the largest event in the eastern Pacific for the last 150 kyr: CaCO₃MARs are 2–3 times higher at 18 ka than elsewhere in the record, including MIS 6. It is caused by high CaCO₃ production rather than minimal dissolution. EOF 2, while it resembles deep water flow patterns, nevertheless, shows coherence to Corg deposition and is probably also driven by CaCO₃production

Erosion by rivers and transport pathways in the ocean: A provenance tool using 40Ar--39Ar incremental heating on fine-grained sediment

We use Ar-40-Ar-39 incremental heating to fingerprint bulk fluvial sediment entering the northeast Pacific Ocean, with the long-term intent of tracking sediment source and transport changes from the terrestrial system to the marine environment through time. We show reproducible age spectra from individual rivers accounting for the majority of sediment delivered to the Pacific margin. Two tests are performed to confirm the validity of the bulk sediment Ar-40-Ar-39 incremental heating measurements and to address why polymineralic sediment might yield concordant age steps. The first model tests, in light of bulk mineralogy and diffusion of Ar from silicates, whether measured K/Ca spectra (measured from Ar-39 and Ar-37, respectively) are consistent with typical values for K- and Ca-bearing minerals. Calculations show that the bulk mineralogy is reflected in the outgassing K/Ca spectra and identify plagioclase as the dominant mineral contributing to the plateau-defining portion of the age spectra. A second model predicts bulk sediment ages from integrated bedrock cooling age-area estimates in order to examine whether bulk sediment plateau ages are representative of the average cooling age of rocks from a given river basin. Calculated and observed ages are notably similar in three river basins when topographic and lithologic effects are accounted for. Overall, this technique shows considerable promise, not only in tracking individual terrigenous sources in the marine realm but also for understanding processes such as erosion and sediment transport in terrestrial systems

Late Pleistocene paleoclimatology of the central equatorial Pacific: A quantitative record of eolian and carbonate deposition

Detailed records of [delta]18O, [delta]13C, percentage and mass accumulation rate of CaCO3, and eolian percentage, mass accumulation rate, and grainsize generated for core RC11-210 from the equatorial Pacific reveal the timing of paleoclimatic events over the past 950,000 yr. The CaCO3 percentage record shows the standard Pacific correlation of high CaCO3 content with glacial periods, but displays a marked change of character about 490,000 yr ago with older stages showing much less variability. The carbonate mass flux record, however, does not show such a noticeable change. Sedimentation rates vary from about 0.5 to 3.0 cm/1000 yr and, during the past 490,000 yr, sections with enhanced sedimentation rates correspond to periods of high CaCO3 percentage. Eolian mass accumulation rates, an indication of the aridity of the source region, are usually higher during glacial times. Eolian grainsize, an indication of the intensity of atmospheric circulation, generally fluctuates at a higher frequency than the 100,000-yr glacial cycle. The mid-Brunhes climatic event centered at 300,000 yr ago appears as a 50,000-yr interval of low intensity and reduced variability of atmospheric circulation. Furthermore, the nature of this entire record changes then, with the younger portion indicating less variation in wind intensity than the older part of the record. The late Matuyama increase in amplitude of paleoclimatic signals begins 875,000 yr ago in the eolian record, 25,000 yr before the [delta]18O and CaCO3 percentage amplitude increases about 850,000 yr ago.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/26521/1/0000059.pd

Modes of Global Climate Variability during Marine Isotope Stage 3 (60–26 ka)

Recent analysis of 38 globally distributed paleoclimatic records covering Marine Isotope Stage 3 (MIS 3) 60–26 ka demonstrated that the two leading empirical orthogonal functions (EOFs) explaining the data are the Greenland ice-core signal (“northern” signal) and the Antarctic ice-core signal (“southern” signal). Here singular spectral analysis (SSA) is used to show that millennial-scale variability of each of these two leading EOFs is characterized by two independent modes. The two modes of each EOF share similar relative distributions of variance, identical spectra, and, where each mode has spectral power, coherency spectra, which are significantly above the null hypothesis level at 95% confidence. The only difference between the modes of the northern and southern signals is that they are phase shifted. The phasing and long response time of the low-frequency mode, combined with its relationship to atmospheric CO2 and sea level, are consistent with coupled changes in the ocean, ice sheets, atmosphere, and carbon cycle, whereas the phasing and short response time of the high-frequency mode are consistent with an atmospheric transmission likely induced by changes in hemispheric sea ice distributions and attendant feedbacks.Keywords: Empirical orthogonal functions, Climate variabilit

Multivariate statistical analysis and partitioning of sedimentary geochemical data sets: General principles and specific MATLAB scripts

Multivariate statistical treatments of large data sets in sedimentary geochemical and other fields are rapidly becoming more popular as analytical and computational capabilities expand. Because geochemical data sets present a unique set of conditions (e.g., the closed array), application of generic off-the-shelf applications is not straightforward and can yield misleading results. We present here annotated MATLAB scripts (and specific guidelines for their use) for Q-mode factor analysis, a constrained least squares multiple linear regression technique, and a total inversion protocol, that are based on the well-known approaches taken by Dymond (1981), Leinen and Pisias (1984), Kyte et al. (1993), and their predecessors. Although these techniques have been used by investigators for the past decades, their application has been neither consistent nor transparent, as their code has remained in-house or in formats not commonly used by many of today’s researchers (e.g., FORTRAN). In addition to providing the annotated scripts and instructions for use, we discuss general principles to be considered when performing multivariate statistical treatments of large geochemical data sets, provide a brief contextual history of each approach, explain their similarities and differences, and include a sample data set for the user to test their own manipulation of the scripts.Keywords: geochemistry, statistics, multivariate statistic