Nicklas G. Pisias

Age dating and the orbital theory of the ice ages: Development of a high-resolution 0 to 300,000-year chronostratigraphy

Using the concept of “orbital tuning”, a continuous, high-resolution deep-sea chronostratigraphy has been developed spanning the last 300,000 yr. The chronology is developed using a stacked oxygen-isotope stratigraphy and four different orbital tuning approaches, each of which is based upon a different assumption concerning the response of the orbital signal recorded in the data. Each approach yields a separate chronology. The error measured by the standard deviation about the average of these four results (which represents the “best” chronology) has an average magnitude of only 2500 yr. This small value indicates that the chronology produced is insensitive to the specific orbital tuning technique used. Excellent convergence between chronologies developed using each of five different paleoclimatological indicators (from a single core) is also obtained. The resultant chronology is also insensitive to the specific indicator used. The error associated with each tuning approach is estimated independently and propagated through to the average result. The resulting error estimate is independent of that associated with the degree of convergence and has an average magnitude of 3500 yr, in excellent agreement with the 2500-yr estimate. Transfer of the final chronology to the stacked record leads to an estimated error of ±1500 yr. Thus the final chronology has an average error of ±5000 yr.

On the Structure and Origin of Major Glaciation Cycles 1. Linear Responses to Milankovitch Forcing

Time series of ocean properties provide a measure of global ice volume and monitor key features of the wind-driven and density-driven circulations over the past 400,000 years. Cycles with periods near 23,000, 41,000, and 100,000 years dominate this climatic narrative. When the narrative is examined in a geographic array of time series, the phase of each climatic oscillation is seen to progress through the system in essentially the same geographic sequence in all three cycles. We argue that the 23,000- and 41,000-year cycles of glaciation are continuous, linear responses to orbitally driven changes in the Arctic radiation budget; and we use the phase progression in each climatic cycle to identify the main pathways along which the initial, local responses to radiation are propagated by the atmosphere and ocean. Early in this progression, deep waters of the Southern Ocean appear to act as a carbon trap. To stimulate new observations and modeling efforts, we offer a process model that gives a synoptic view of climate at the four end-member states needed to describe the systems evolution, and we propose a dynamic system model that explains the phase progression along causal pathways by specifying inertial constants in a chain of four subsystems. “Solutions to problems involving systems of such complexity are not born full grown like Athena from the head of Zeus. Rather they evolve slowly, in stages, each of which requires a pause to examine data at great lengths in order to guarantee a sure footing and to properly choose the next step.” —Victor P. Starr

On the structure and origin of major glaciation cycles 2. The 100,000‐year cycle

Climate over the past million years has been dominated by glaciation cycles with periods near 23,000, 41,000, and 100,000 years. In a linear version of the Milankovitch theory, the two shorter cycles can be explained as responses to insolation cycles driven by precession and obliquity. But the 100,000-year radiation cycle (arising from eccentricity variation) is much too small in amplitude and too late in phase to produce the corresponding climate cycle by direct forcing. We present phase observations showing that the geographic progression of local responses over the 100,000-year cycle is similar to the progression in the other two cycles, implying that a similar set of internal climatic mechanisms operates in all three. But the phase sequence in the 100,000-year cycle requires a source of climatic inertia having a time constant (similar to 15,000 years) much larger than the other cycles (similar to 5,000 years). Our conceptual model identifies massive northern hemisphere ice sheets as this larger inertial source. When these ice sheets, forced by precession and obliquity, exceed a critical size, they cease responding as linear Milankovitch slaves and drive atmospheric and oceanic responses that mimic the externally forced responses. In our model, the coupled system acts as a nonlinear amplifier that is particularly sensitive to eccentricity-driven modulations in the 23,000-year sea level cycle. During an interval when sea level is forced upward from a major low stand by a Milankovitch response acting either alone or in combination with an internally driven, higher-frequency process, ice sheets grounded on continental shelves become unstable, mass wasting accelerates, and the resulting deglaciation sets the phase of one wave in the train of 100,000-year oscillations. Whether a glacier or ice sheet influences the climate depends very much on the scale....The interesting aspect is that an effect on the local climate can still make an ice mass grow larger and larger, thereby gradually increasing its radius of influence.

The role of the thermohaline circulation in abrupt climate change

The possibility of a reduced Atlantic thermohaline circulation in response to increases in greenhouse-gas concentrations has been demonstrated in a number of simulations with general circulation models of the coupled ocean–atmosphere system. But it remains difficult to assess the likelihood of future changes in the thermohaline circulation, mainly owing to poorly constrained model parameterizations and uncertainties in the response of the climate system to greenhouse warming. Analyses of past abrupt climate changes help to solve these problems. Data and models both suggest that abrupt climate change during the last glaciation originated through changes in the Atlantic thermohaline circulation in response to small changes in the hydrological cycle. Atmospheric and oceanic responses to these changes were then transmitted globally through a number of feedbacks. The palaeoclimate data and the model results also indicate that the stability of the thermohaline circulation depends on the mean climate state.

Foraminiferal faunal estimates of paleotemperature: Circumventing the No‐analog problem yields cool Ice Age tropics

The sensitivity of the tropics to climate change, particularly the amplitude of glacial-to-interglacial changes in sea surface temperature (SST), is one of the great controversies in paleoclimatology. Here we reassess faunal estimates of ice age SSTs, focusing on the problem of no-analog planktonic foraminiferal assemblages in the equatorial oceans that confounds both classical transfer function and modern analog methods. A new calibration strategy developed here, which uses past variability of species to define robust faunal assemblages, solves the no-analog problem and reveals ice age cooling of 5 o to 6oC in the equatorial current systems of the Atlantic and eastern Pacific Oceans. Classical transfer functions underestimated temperature changes in some areas of the tropical oceans because core-top assemblages misrepresented the ice age faunal assemblages. Our finding is consistent with some geochemical estimates and model predictions of greater ice age cooling in the tropics than was inferred by Climate: Long-Range Investigation, Mapping, and Prediction (CLIMAP) (1981) and thus may help to resolve a long-standing controversy. Our new foraminiferal transfer function suggests that such cooling was limited to the equatorial current systems, however, and supports CLIMAPs inference of stability of the subtropical gyre centers.

Ferromanganese nodules from MANOP Sites H, S, and R-Control of mineralogical and chemical composition by multiple accretionary processes

Abstract The chemical composition of ferromanganese nodules from the three nodule-bearing MANOP sites in the Pacific can be accounted for in a qualitative way by variable contributions of distinct accretionary processes. These accretionary modes are: 1. (1) hydrogenous, i.e., direct precipitation or accumulation of colloidal metal oxides in seawater, 2. (2) oxic diagenesis which refers to a variety of ferromanganese accretion processes occurring in oxic sediments; and 3. (3) suboxic diagenesis which results from reduction of Mn+4 by oxidation of organic matter in the sediments. Geochemical evidence suggests processes (1) and (2) occur at all three MANOP nodule-bearing sites, and process (3) occurs only at the hemipelagic site, H, which underlies the relatively productive waters of the eastern tropical Pacific. A normative model quantitatively accounts for the variability observed in nearly all elements. Zn and Na, however, are not well explained by the three end-member model, and we suggest that an additional accretionary process results in greater variability in the abundances of these elements. Variable contributions from the three accretionary processes result in distinct top-bottom compositional differences at the three sites. Nodule tops from H are enriched in Ni, Cu, and Zn, instead of the more typical enrichments of these elements in nodule bottoms. In addition, elemental correlations typical of most pelagic nodules are reversed at site H. The three accretionary processes result in distinct mineralogies. Hydrogenous precipitation produces δMnO2. Oxic diagenesis, however, produces Cu-Ni-rich todorokite, and suboxic diagenesis results in an unstable todorokite which transforms to a 7 A phase (“birnessite”) upon dehydration. The presence of Cu and Ni as charge-balancing cations influence the stability of the todorokite structure. In the bottoms of H nodules, which accrete dominantly by suboxic diagenesis, Na+ and possibly Mn+2 provide much of the charge balance for the todorokite structure. Limited growth rate data for H nodules suggest suboxic accretion is the fastest of the three processes, with rates at least 200 mm/106 yr. Oxic accretion is probably 10 times slower and hydrogenous 100 times slower. Since these rates predict more suboxic component in bulk nodules than is calculated by the normative analysis, we propose that suboxic accretion is a non-steady-state process. Variations in surface water productivity cause pulses of particulate flux to the sea floor which result in transient Mn reduction in the surface sediments and reprecipitation on nodule surfaces.

High resolution stratigraphic correlation of benthic oxygen isotopic records spanning the last 300,000 years

Abstract In order to compare the response of different oceanographic regions to global climate change, very detailed stratigraphic techniques are required. The global signal of ice volume changes recorded in the oxygen isotopic composition of foraminifera can provide the tool for developing the necessary high resolution stratigraphy. In order to evaluate the resolution of a stratigraphy based on detailed isotopic records, two techniques are used to correlate a set of benthic oxygen isotope records from seven piston cores taken in the North and South Atlantic, the Indian, and the equatorial and North Pacific oceans. The first technique was modified from the graphic correlation procedure of Shaw (1964). This procedure requires the identification of isotopic events that are correlated from core to core. Detailed correlations for intervals of the cores between events are provided by a series of straight line segments connecting all common events. The second technique developed by Martinson et al. (1982) uses inverse procedures to define a continuous non-linear mapping function that correlates the isotopic records. The mapping function maximizes the correlation coefficient between data sets being compared. The techniques are independent in that they rely on different criteria for correlating the data series. Stratigraphic correlations obtained by these procedures are in excellent agreement. The mean difference between the correlations is on the order of the sampling intervals of each core and, when corrected for sedimentation rates, suggests that benthic isotope records from the suite of seven cores can be correlated to a resolution of 2000 to 4000 yrs.

SPATIAL AND TEMPORAL OCEANOGRAPHIC VARIABILITY OF THE EASTERN EQUATORIAL PACIFIC DURING THE LATE PLEISTOCENE : EVIDENCE FROM RADIOLARIA MICROFOSSILS

Eight 150,000 year long records of sea surface temperatures combined with two additional records spanning 400,000 years constrain the spatial and temporal patterns of oceanographic change in the eastern equatorial Pacific and possible mechanisms of variability in the region. Empirical orthogonal function analysis shows two important modes of variability, one associated with the eastern boundary current and another associated with the North Equatorial Countercurrent. The two long time series located in the equatorial divergence and within the Peru Current have very different patterns of change. The spectrum for the time series from the Peru Current is dominated by orbital periods of 100, 41, and 23 kyr and is similar in variance distribution and phase to records from the Southern Ocean. In contrast, the equatorial divergence site has spectral concentrations at the orbital frequencies and also concentration of variance at the nonorbital 31,000 year period. The phase and amplitude spectra of these two sites support the importance of changes in eastern boundary advection and also document a nonlinear response of the equatorial Pacific to orbital changes. Finally, these data provide a new evaluation of the temperature change in the eastern equatorial Pacific during the last glacial maximum. Cooling in the Peru Current region is predicted to be about 4°C, and cooling in the equatorial divergence is estimated to be 3° to 5°C. The estimated cooling of the region is of the order of 2°C greater than the cooling predicted by Climate: Long-Range Investigation, Mapping, and Prediction (CLIMAP).

Paleoproductivity and carbon burial across the California Current: The multitracers transect, 42°N

The Multitracers Experiment studied a transect of water column, sediment trap, and sediment data taken across the California Current to develop quantitative methods for hindcasting paleoproductivity. The experiment used three sediment trap moorings located 120 km, 270 km, and 630 km from shore at the Oregon/California border in North America. We report here about the sedimentation and burial of particulate organic carbon (Corg) and CaCO3. In order to observe how the integrated CaCO3 and Corg burial across the transect has changed since the last glacial maximum, we have correlated core from the three sites using time scales constrained by both radiocarbon and oxygen isotopes. By comparing surface sediments to a two-and-a-half year sediment trap record, we have also defined the modern preservation rates for many of the labile sedimentary materials. Our analysis of the Corg data indicates that significant amounts (20–40%) of the total Corg being buried today in surface sediments is terrestrial. At the last glacial maximum, the terrestrial Corg fraction within 300 km of the coast was about twice as large. Such large fluxes of terrestrial Corg obscure the marine Corg record, which can be interpreted as productivity. When we corrected for the terrestrial organic matter, we found that the mass accumulation rate of marine Corg roughly doubled from the glacial maximum to the present. Because preservation rates of organic carbon are high in the high sedimentation rate cores, corrections for degradation are straightforward and we can be confident that organic carbon rain rate (new productivity) also doubled. As confirmation, the highest burial fluxes of other biogenic components (opal and Ba) also occur in the Holocene. Productivity off Oregon has thus increased dramatically since the last glacial maximum. CaCO3 fluxes also changed radically through the deglaciation; however, they are linked not to CaCO3 production but rather to changes in deepwater carbonate chemistry between 18 Ka and now.

Carbon 13 in Pacific Deep and Intermediate Waters, 0‐370 ka: Implications for Ocean Circulation and Pleistocene CO2

Stable isotopes in benthic foraminifera from Pacific sediments are used to assess hypotheses of systematic shifts in the depth distribution of oceanic nutrients and carbon during the ice ages. The carbon isotope differences between ∼1400 and ∼3200 m depth in the eastern Pacific are consistently greater in glacial than interglacial maxima over the last ∼370 kyr. This phenomenon of “bottom heavy” glacial nutrient distributions, which Boyle proposed as a cause of Pleistocene CO2 change, occurs primarily in the 1/100 and 1/41 kyr−1 “Milankovitch” orbital frequency bands but appears to lack a coherent 1/23 kyr−1 band related to orbital precession. Averaged over oxygen-isotope stages, glacial δ13C gradients from ∼1400 to ∼3200 m depth are 0.1‰ greater than interglacial gradients. The range of extreme shifts is somewhat larger, 0.2 to 0.5‰. In both cases, these changes in Pacific δ13C distributions are much smaller than observed in shorter records from the North Atlantic. This may be too small to be a dominant cause of atmospheric pCO2 change, unless current models underestimate the sensitivity of pCO2 to nutrient redistributions. This dampening of Pacific relative to Atlantic δ13C depth gradient favors a North Atlantic origin of the phenomenon, although local variations of Pacific intermediate water masses can not be excluded at present.