Alan C. Mix

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.

Global warming preceded by increasing carbon dioxide concentrations during the last deglaciation

The covariation of carbon dioxide (CO2) concentration and temperature in Antarctic ice-core records suggests a close link between CO2 and climate during the Pleistocene ice ages. The role and relative importance of CO2 in producing these climate changes remains unclear, however, in part because the ice-core deuterium record reflects local rather than global temperature. Here we construct a record of global surface temperature from 80 proxy records and show that temperature is correlated with and generally lags CO2 during the last (that is, the most recent) deglaciation. Differences between the respective temperature changes of the Northern Hemisphere and Southern Hemisphere parallel variations in the strength of the Atlantic meridional overturning circulation recorded in marine sediments. These observations, together with transient global climate model simulations, support the conclusion that an antiphased hemispheric temperature response to ocean circulation changes superimposed on globally in-phase warming driven by increasing CO2 concentrations is an explanation for much of the temperature change at the end of the most recent ice age.

A Reconstruction of Regional and Global Temperature for the Past 11,300 Years

Exceptional Now The climate has been warming since the industrial revolution, but how warm is climate now compared with the rest of the Holocene? Marcott et al. (p. 1198) constructed a record of global mean surface temperature for more than the last 11,000 years, using a variety of land- and marine-based proxy data from all around the world. The pattern of temperatures shows a rise as the world emerged from the last deglaciation, warm conditions until the middle of the Holocene, and a cooling trend over the next 5000 years that culminated around 200 years ago in the Little Ice Age. Temperatures have risen steadily since then, leaving us now with a global temperature higher than those during 90% of the entire Holocene. Current global average surface air temperature is warmer than that for all but a small fraction of the past 11,300 years. Surface temperature reconstructions of the past 1500 years suggest that recent warming is unprecedented in that time. Here we provide a broader perspective by reconstructing regional and global temperature anomalies for the past 11,300 years from 73 globally distributed records. Early Holocene (10,000 to 5000 years ago) warmth is followed by ~0.7°C cooling through the middle to late Holocene (<5000 years ago), culminating in the coolest temperatures of the Holocene during the Little Ice Age, about 200 years ago. This cooling is largely associated with ~2°C change in the North Atlantic. Current global temperatures of the past decade have not yet exceeded peak interglacial values but are warmer than during ~75% of the Holocene temperature history. Intergovernmental Panel on Climate Change model projections for 2100 exceed the full distribution of Holocene temperature under all plausible greenhouse gas emission scenarios.

Ice sheets and sea level of the Last Glacial Maximum

This paper outlines the general issues regarding ice sheets and sea level of the Last Glacial Maximum (LGM), which formed the basis of an EPILOG Project workshop. Papers in this special issue of Quaternary Science Reviews provide a comprehensive assessment of these issues from the perspective of geological reconstructions of ice sheet extent, records of sea-level change, ice sheet modelling, geophysical models of glacial isostatic adjustment, and geochemical proxies of ice volume. This new assessment has substantially narrowed the uncertainties in the total changes in ice sheets and sea level and their proxies, suggesting a net decrease in eustatic sea level at the LGM ranging from 120 to 135 m. r 2001 Elsevier Science Ltd. All rights reserved.

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.

Oceanic Response to Orbital Forcing in the Late Quaternary: Observational and Experimental Strategies

Observations on deep-sea cores demonstrate that late Pleistocene climate is dominated by three broad-band cycles centered near periods of 23 ky, 41 ky, and 100 ky. These cycles permeate the global system, and include changes in the atmosphere, cryosphere, surface ocean, and deep ocean. The periods of these climatic cycles match orbital cycles of precession, obliquity, and eccentricity; and each orbit-climate pair is significantly correlated (coherent). These observations may be explained in different ways. We review various models and conclude that the climatic cycles can be explained as an interaction between orbitally forced and internally driven oscillations of the climate system. Depending on the cycle and model, the external forcing may influence climate either as part of a driving mechanism which determines both the amplitude and phase of the cycle, or as a pacing mechanism which sets the phase of an internal oscillation. Our goal is to search climatic data for clues about the mechanisms which operate within the climate system on orbital time scales. Our strategy is patterned after previous investigations which partition the climatic record into cyclic components, record the phase of system responses in each cycle, and examine these phase sequences for clues about the chain of causal mechanisms.

Global climate evolution during the last deglaciation

Deciphering the evolution of global climate from the end of the Last Glacial Maximum approximately 19 ka to the early Holocene 11 ka presents an outstanding opportunity for understanding the transient response of Earth’s climate system to external and internal forcings. During this interval of global warming, the decay of ice sheets caused global mean sea level to rise by approximately 80 m; terrestrial and marine ecosystems experienced large disturbances and range shifts; perturbations to the carbon cycle resulted in a net release of the greenhouse gases CO2 and CH4 to the atmosphere; and changes in atmosphere and ocean circulation affected the global distribution and fluxes of water and heat. Here we summarize a major effort by the paleoclimate research community to characterize these changes through the development of well-dated, high-resolution records of the deep and intermediate ocean as well as surface climate. Our synthesis indicates that the superposition of two modes explains much of the variability in regional and global climate during the last deglaciation, with a strong association between the first mode and variations in greenhouse gases, and between the second mode and variations in the Atlantic meridional overturning circulation.

Photosynthetic fractionation of 13C and concentrations of dissolved CO2 in the central equatorial Pacific during the last 255,000 years

Carbon isotopically based estimates of CO2 levels have been generated from a record of the photosynthetic fractionation of 13C [is equivalent to epsilon(p)] in a central equatorial Pacific sediment core that spans the last approximately 255 ka. Contents of 13C in phytoplanktonic biomass were determined by analysis of C37 alkadienones. These compounds are exclusive products of Prymnesiophyte algae which at present grow most abundantly at depths of 70-90 m in the central equatorial Pacific. A record of the isotopic composition of dissolved CO2 was constructed from isotopic analyses of the planktonic foraminifera Neogloboquadrina dutertrei, which calcifies at 70-90 m in the same region. Values of epsilon(p), derived by comparison of the organic and inorganic delta values, were transformed to yield concentrations of dissolved CO2 [is equivalent to c(e)] based on a new, site-specific calibration of the relationship between epsilon(p) and c(e). The calibration was based on reassessment of existing epsilon(p) versus c(e) data, which support a physiologically based model in which epsilon(p) is inversely related to c(e). Values of PCO2, the partial pressure of CO2 that would be in equilibrium with the estimated concentrations of dissolved CO2, were calculated using Henrys law and the temperature determined from the alkenone-unsaturation index U(K/37). Uncertainties in these values arise mainly from uncertainties about the appropriateness (particularly over time) of the site-specific relationship between epsilon(p) and 1/c(e). These are discussed in detail and it is concluded that the observed record of epsilon(p) most probably reflects significant variations in delta pCO2, the ocean-atmosphere disequilibrium, which appears to have ranged from approximately 110 microatmospheres during glacial intervals (ocean > atmosphere) to approximately 60 microatmospheres during interglacials. Fluxes of CO2 to the atmosphere would thus have been significantly larger during glacial intervals. If this were characteristic of large areas of the equatorial Pacific, then greater glacial sinks for the equatorially evaded CO2 must have existed elsewhere. Statistical analysis of air-sea pCO2 differences and other parameters revealed significant (p<0.01) inverse correlations of delta pCO2 with sea surface temperature and with the mass accumulation rate of opal. The former suggests response to the strength of upwelling, the latter may indicate either drawdown of CO2 by siliceous phytoplankton or variation of [CO2]/[Si(OH)4] ratios in upwelling waters.

Climate Sensitivity Estimated from Temperature Reconstructions of the Last Glacial Maximum

Last Glacial Maximum temperature reconstructions and model simulations can constrain the equilibrium climate sensitivity. Assessing the impact of future anthropogenic carbon emissions is currently impeded by uncertainties in our knowledge of equilibrium climate sensitivity to atmospheric carbon dioxide doubling. Previous studies suggest 3 kelvin (K) as the best estimate, 2 to 4.5 K as the 66% probability range, and nonzero probabilities for much higher values, the latter implying a small chance of high-impact climate changes that would be difficult to avoid. Here, combining extensive sea and land surface temperature reconstructions from the Last Glacial Maximum with climate model simulations, we estimate a lower median (2.3 K) and reduced uncertainty (1.7 to 2.6 K as the 66% probability range, which can be widened using alternate assumptions or data subsets). Assuming that paleoclimatic constraints apply to the future, as predicted by our model, these results imply a lower probability of imminent extreme climatic change than previously thought.