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)

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

Synorogenic evolution of large-scale drainage patterns: Isotope paleohydrology of sequential Laramide basins

In the past decade, we and others have compiled an extensive dataset of O, C and Sr isotope stratigraphies from sedimentary basins throughout the Paleogene North American Cordillera. In this study, we present new results from the Piceance Creek Basin of northwest Colorado, which record the evolving hydrology of the Eocene Green River Lake system. We then place the new data in the context of the broader Cordilleran dataset and summarize implications for understanding the synorogenic evolution of large-scale drainage patterns. The combined data reflect (1) a period of throughgoing foreland rivers heading in the Sevier fold-and-thrust belt and flowing east, (2) ponding of freshwater lakes in the foredeep as Laramide uplifts blocked drainage, (3) hydrologic closure that led to both intensive evaporation in the terminal sink of the Piceance Creek Basin and integration of catchments over length-scales \u3e1000 km, (4) infilling of basin accommodation by southward migrating magmatism in distal catchments, leading to the freshening and demise of intraforeland lakes that also stepped south over time

Terrigenous Fe input and biogenic sedimentation in the glacial and interglacial equatorial Pacific Ocean

Many ocean regions important to the global carbon budget, including the equatorial Pacific Ocean, have low chlorophyll concentrations despite high levels of conventional nutrients. Iron may instead be the limiting nutrient, and elevated input of terrigenous Fe during windy glacial episodes has been hypothesized to stimulate oceanic productivity through time and thus regulate the oceanic and atmospheric CO2 balance. To test whether particulate Fe input is related to the accumulation of biogenic matter in one important low chlorophyll‐high nutrient area, that is, the equatorial Pacific Ocean, we present results from a suite of sediment cores that collectively record biogenic deposition through the last six glacial‐interglacial cycles (∼600,000 years). Our data set includes new chemical data on total Fe, terrigenous, and biogenic components in three cores as well as previously published mineralogic records of eolian input to the region. Chemical, spectral, and stratigraphic analysis indicates that (1) terrigenous input to the region shows no consistent pattern of either glacial or interglacial maxima, (2) the accumulation of particulate Fe is closely related to the accumulation of terrigenous matter (linear r2 = 0.81–0.98), (3) there are no coherent spectral relationships between Fe input and glacial periodicity (i.e., δ18O) in any of the orbital frequency bands, (4) the linear and cross‐spectral correlations between Fe or eolian input and CaCO3 concentration are most commonly the strongest observed relationships between Fe and any biogenic component, yet indicate a largely inverse pattern, with higher Fe being associated with low CaCO3, (5) there is no consistent linear r2 correlation or spectral coherence between the accumulation of Fe and that of CaCO3, Corg, or opal. Thus in total there is no relationship between terrigenous Fe input and sedimentary sequestering of carbon. Additionally, although we cannot specifically address the potential for changes in solubility of the terrigenous fraction that may be driven by a terrigenous compositional change, the Fe/Ti ratio (which monitors first‐order mineralogic changes) records only slight variations that also are linearly and spectrally unrelated to glacial periodicity, the bulk Fe flux, and the accumulation of any biogenic component. Finally, we find that the paleoceanographic flux of Fe is several order‐of‐magnitudes larger than modern observations of eolian Fe input, suggesting that the long‐term importance of Fe input by dust storms (which deliver Fe on the order of the sedimentary burial) may be underestimated. The removal of particulate terrigenous Fe from the recently discovered source within the Equatorial Undercurrent, however, remains unquantified and may also prove significant

Increased dissolved terrestrial input to the coastal ocean during the last deglaciation

Here we present the first downcore results for a new paleoproxy, the Mn/Ca ratio of foraminiferal calcite, applied to sediment accumulated in the extreme Eastern Tropical North Pacific (ETNP) over the last 30,000 years. The Mn/Ca results are compared to oxygen isotopes and sea surface temperature calculated from Mg/Ca. We determined metal ratios using flow-through time-resolved analysis to minimize the effects of secondary mineralization. The foraminiferal species used for this study calcify at different depths. Core top ratios of these variant species change in proportion to the concentration of dissolved manganese in the water column at the depth of calcification. Since terrestrial input and oxidation reduction reactions control the levels of dissolved Mn in the oceans today, it therefore should be possible to use the Mn/Ca ratios of foraminifera as a proxy for these processes in the past. Mn/Ca of a mixed-layer species (G. ruber) suggest that dissolved terrestrial input to the surface waters of the ETNP during the last glacial maximum was lower than today but began to increase with initial sea level rise and reached a maximum at 15 ka B.P. before coming down to present-day levels at the end of sea level rise in the mid-Holocene (7–5 ka). Ratios of a deeper calcifying species (N. dutertrei) mimic those of G. ruber over this same time period, consistent with shoaling of the 18°C thermocline. Mn/Ca of a benthic species (U. peregrina) does not show a maximum at 15 ka, suggesting that Mn was efficiently remineralized in the water column during deglaciation. Assuming that the period from the last glacial until the mid-Holocene was a time of increased productivity, as elevated Mn might imply, the oxygen minimum zone (OMZ) was at least as well developed during deglaciation as it is today. Expansion of the OMZ may have contributed to the Mn/Ca trends we observe through time

The Atlantic Ocean at the last glacial maximum: 1. Objective mapping of the GLAMAP sea-surface conditions

Recent efforts of the German paleoceanographic community have resulted in a unique data set of reconstructed sea-surface temperature for the Atlantic Ocean during the Last Glacial Maximum, plus estimates for the extents of glacial sea ice. Unlike prior attempts, the contributing research groups based their data on a common definition of the Last Glacial Maximum chronozone and used the same modern reference data for calibrating the different transfer techniques. Furthermore, the number of processed sediment cores was vastly increased. Thus the new data is a significant advance not only with respect to quality, but also to quantity. We integrate these new data and provide monthly data sets of global sea-surface temperature and ice cover, objectively interpolated onto a regular 1°x1° grid, suitable for forcing or validating numerical ocean and atmosphere models. This set is compared to an existing subjective interpolation of the same base data, in part by employing an ocean circulation model. For the latter purpose, we reconstruct sea surface salinity from the new temperature data and the available oxygen isotope measurements

Reduced admixture of North Atlantic Deep Water to the deep central South Pacific during the last two glacial periods

Key Points: • Little deep water circulation changes in the past 240,000 years in the central South Pacific • Reduced North Atlantic Deep Water admixture during glacials to the Southern Ocean • South Pacific lithogenic material mainly sourced from SE Australia and South New Zealand The South Pacific is a sensitive location for the variability of the global oceanic thermohaline circulation given that deep waters from the Atlantic Ocean, the Southern Ocean, and the Pacific basin are exchanged. Here we reconstruct the deep-water circulation of the central South Pacific for the last two glacial cycles (from 240,000 years ago to the Holocene) based on radiogenic neodymium (Nd) and lead (Pb) isotope records complemented by benthic stable carbon data obtained from two sediment cores located on the flanks of the East Pacific Rise. The records show small but consistent glacial/interglacial changes in all three isotopic systems with interglacial average values of -5.8 and 18.757 for εNd and 206Pb/204Pb, respectively, whereas glacial averages are -5.3 and 18.744. Comparison of this variability of Circumpolar Deep Water (CDW) to previously published records along the pathway of the global thermohaline circulation is consistent with reduced admixture of North Atlantic Deep Water (NADW) to CDW during cold stages. The absolute values and amplitudes of the benthic δ13C variations are essentially indistinguishable from other records of the Southern Hemisphere and confirm that the low central South Pacific sedimentation rates did not result in a significant reduction of the amplitude of any of the measured proxies. In addition, the combined detrital Nd and strontium (87Sr/86Sr) isotope signatures imply that Australian and New Zealand dust has remained the principal contributor of lithogenic material to the central South Pacific

Increased dissolved terrestrial input to the coastal ocean during the last deglaciation

Here we present the first downcore results for a new paleoproxy, the Mn/Ca ratio of foraminiferal calcite, applied to sediment accumulated in the extreme Eastern Tropical North Pacific (ETNP) over the last 30,000 years. The Mn/Ca results are compared to oxygen isotopes and sea surface temperature calculated from Mg/Ca. We determined metal ratios using flow‐through time‐resolved analysis to minimize the effects of secondary mineralization. The foraminiferal species used for this study calcify at different depths. Core top ratios of these variant species change in proportion to the concentration of dissolved manganese in the water column at the depth of calcification. Since terrestrial input and oxidation reduction reactions control the levels of dissolved Mn in the oceans today, it therefore should be possible to use the Mn/Ca ratios of foraminifera as a proxy for these processes in the past. Mn/Ca of a mixed‐layer species (G. ruber) suggest that dissolved terrestrial input to the surface waters of the ETNP during the last glacial maximum was lower than today but began to increase with initial sea level rise and reached a maximum at 15 ka B.P. before coming down to present‐day levels at the end of sea level rise in the mid‐Holocene (7–5 ka). Ratios of a deeper calcifying species (N. dutertrei) mimic those of G. ruber over this same time period, consistent with shoaling of the 18°C thermocline. Mn/Ca of a benthic species (U. peregrina) does not show a maximum at 15 ka, suggesting that Mn was efficiently remineralized in the water column during deglaciation. Assuming that the period from the last glacial until the mid‐Holocene was a time of increased productivity, as elevated Mn might imply, the oxygen minimum zone (OMZ) was at least as well developed during deglaciation as it is today. Expansion of the OMZ may have contributed to the Mn/Ca trends we observe through time.Keywords: manganese, paleoprox

The Cenozoic Climatic and Topographic Evolution of the Western North American Cordillera

Herein we present oxygen isotope records from Cretaceous to Recent terrestrial sediments in the western North American Cordillera. The purpose of this analysis is to use oxygen isotope records to understand the coupled surface elevation and climate histories of this region through the Cenozoic. To do this we constructed δ¹⁸(O) maps of surface waters for time intervals that trace the development of topography of western North America. These maps are based on 4861 oxygen isotope analyses from both published (4478) and new (383) data. We determined the δ¹⁸(O) values of surface waters using paleotemperatures derived previously from floral assemblages and the appropriate isotope fractionation factors. These data suggest that in the late Cretaceous to early Eocene the Sevier hinterland formed a plateau of unknown height. Around 50 Ma, a topographic wave developed in British Columbia and eastern Washington that swept southward reaching northeastern Nevada at ~40 to 38 Ma, and southern Nevada ~23 Ma. This southward encroachment of an Eocene Plateau (SWEEP) caused reorganization of drainage patterns such that the intraforeland basins of Wyoming and Utah drainages extended deep within the Sevier hinterland as the wave swept southward. The landscape within the Sevier hinterland developed into a rugged and high mountain range with the hypsometric mean elevation of ~4 km and relief of ~1.5 km. This Eocene highland was bordered on the west by a high Sierra Nevada ramp and on the east by the intraforeland basins that captured water draining these growing highlands. The spatial and temporal evolution of this highland correlates with the timing of volcanism and extension. These observations support tectonic models that call for north to south removal of the Farallon slab or piecemeal removal of mantle lithosphere. The isotopic data show that prior to growth of this highland the North American Monsoon (NAM) penetrated much farther north in the Paleocene/Eocene than today. The combined effects of global cooling, increasing latitudinal temperature gradients, and the generation of the orographic barrier created by the growing north to south highland produced a southward migration of the NAM front. By the Oligocene the hydrologic regime that we observe today was in place. It has been modified since then as a result of Basin and Range extension and collapse of the highlands in the mid-Miocene. This collapse allowed the NAM to penetrate farther north into the Great Basin of Nevada and Utah