Last interglacial temperature evolution — a model inter-comparison

There is a growing number of proxy-based reconstructions detailing the climatic changes that occurred during the last interglacial period (LIG). This period is of special interest, because large parts of the globe were characterized by a warmer-than-present-day climate, making this period an interesting test bed for climate models in light of projected global warming. However, mainly because synchronizing the different palaeoclimatic records is difficult, there is no consensus on a global picture of LIG temperature changes. Here we present the first model inter-comparison of transient simulations covering the LIG period. By comparing the different simulations, we aim at investigating the common signal in the LIG temperature evolution, investigating the main driving forces behind it and at listing the climate feedbacks which cause the most apparent inter-model differences. The model inter-comparison shows a robust Northern Hemisphere July temperature evolution characterized by a maximum between 130–125 ka BP with temperatures 0.3 to 5.3 K above present day. A Southern Hemisphere July temperature maximum, −1.3 to 2.5 K at around 128 ka BP, is only found when changes in the greenhouse gas concentrations are included. The robustness of simulated January temperatures is large in the Southern Hemisphere and the mid-latitudes of the Northern Hemisphere. For these regions maximum January temperature anomalies of respectively −1 to 1.2 K and −0.8 to 2.1 K are simulated for the period after 121 ka BP. In both hemispheres these temperature maxima are in line with the maximum in local summer insolation. In a number of specific regions, a common temperature evolution is not found amongst the models. We show that this is related to feedbacks within the climate system which largely determine the simulated LIG temperature evolution in these regions. Firstly, in the Arctic region, changes in the summer sea-ice cover control the evolution of LIG winter temperatures. Secondly, for the Atlantic region, the Southern Ocean and the North Pacific, possible changes in the characteristics of the Atlantic meridional overturning circulation are crucial. Thirdly, the presence of remnant continental ice from the preceding glacial has shown to be important when determining the timing of maximum LIG warmth in the Northern Hemisphere. Finally, the results reveal that changes in the monsoon regime exert a strong control on the evolution of LIG temperatures over parts of Africa and India. By listing these inter-model differences, we provide a starting point for future proxy-data studies and the sensitivity experiments needed to constrain the climate simulations and to further enhance our understanding of the temperature evolution of the LIG period.

Oxygen and indicators of stress for marine life in multi-model global warming projections

Decadal-to-century scale trends for a range of marine environmental variables in the upper mesopelagic layer (UML, 100–600 m) are investigated using results from seven Earth System Models forced by a high greenhouse gas emission scenario. The models as a class represent the observation-based distribution of oxygen (O2) and carbon dioxide (CO2), albeit major mismatches between observation-based and simulated values remain for individual models. By year 2100 all models project an increase in SST between 2 °C and 3 °C, and a decrease in the pH and in the saturation state of water with respect to calcium carbonate minerals in the UML. A decrease in the total ocean inventory of dissolved oxygen by 2% to 4% is projected by the range of models. Projected O2 changes in the UML show a complex pattern with both increasing and decreasing trends reflecting the subtle balance of different competing factors such as circulation, production, remineralization, and temperature changes. Projected changes in the total volume of hypoxic and suboxic waters remain relatively small in all models. A widespread increase of CO2 in the UML is projected. The median of the CO2 distribution between 100 and 600m shifts from 0.1–0.2 mol m−3 in year 1990 to 0.2–0.4 mol m−3 in year 2100, primarily as a result of the invasion of anthropogenic carbon from the atmosphere. The co-occurrence of changes in a range of environmental variables indicates the need to further investigate their synergistic impacts on marine ecosystems and Earth System feedbacks.

Historical and idealized climate model experiments: an intercomparison of Earth system models of intermediate complexity

Both historical and idealized climate model experiments are performed with a variety of Earth system models of intermediate complexity (EMICs) as part of a community contribution to the Intergovernmental Panel on Climate Change Fifth Assessment Report. Historical simulations start at 850 CE and continue through to 2005. The standard simulations include changes in forcing from solar luminosity, Earth's orbital configuration, CO2, additional greenhouse gases, land use, and sulphate and volcanic aerosols. In spite of very different modelled pre-industrial global surface air temperatures, overall 20th century trends in surface air temperature and carbon uptake are reasonably well simulated when compared to observed trends. Land carbon fluxes show much more variation between models than ocean carbon fluxes, and recent land fluxes appear to be slightly underestimated. It is possible that recent modelled climate trends or climate–carbon feedbacks are overestimated resulting in too much land carbon loss or that carbon uptake due to CO2 and/or nitrogen fertilization is underestimated. Several one thousand year long, idealized, 2 × and 4 × CO2 experiments are used to quantify standard model characteristics, including transient and equilibrium climate sensitivities, and climate–carbon feedbacks. The values from EMICs generally fall within the range given by general circulation models. Seven additional historical simulations, each including a single specified forcing, are used to assess the contributions of different climate forcings to the overall climate and carbon cycle response. The response of surface air temperature is the linear sum of the individual forcings, while the carbon cycle response shows a non-linear interaction between land-use change and CO2 forcings for some models. Finally, the preindustrial portions of the last millennium simulations are used to assess historical model carbon-climate feedbacks. Given the specified forcing, there is a tendency for the EMICs to underestimate the drop in surface air temperature and CO2 between the Medieval Climate Anomaly and the Little Ice Age estimated from palaeoclimate reconstructions. This in turn could be a result of unforced variability within the climate system, uncertainty in the reconstructions of temperature and CO2, errors in the reconstructions of forcing used to drive the models, or the incomplete representation of certain processes within the models. Given the forcing datasets used in this study, the models calculate significant land-use emissions over the pre-industrial period. This implies that land-use emissions might need to be taken into account, when making estimates of climate–carbon feedbacks from palaeoclimate reconstructions.

Carbon dioxide and climate impulse response functions for the computation of greenhouse gas metrics: a multi-model analysis

The responses of carbon dioxide (CO2) and other climate variables to an emission pulse of CO2 into the atmosphere are often used to compute the Global Warming Potential (GWP) and Global Temperature change Potential (GTP), to characterize the response timescales of Earth System models, and to build reduced-form models. In this carbon cycle-climate model intercomparison project, which spans the full model hierarchy, we quantify responses to emission pulses of different magnitudes injected under different conditions. The CO2 response shows the known rapid decline in the first few decades followed by a millennium-scale tail. For a 100 Gt-C emission pulse added to a constant CO2 concentration of 389 ppm, 25 ± 9% is still found in the atmosphere after 1000 yr; the ocean has absorbed 59 ± 12% and the land the remainder (16 ± 14%). The response in global mean surface air temperature is an increase by 0.20 ± 0.12 °C within the first twenty years; thereafter and until year 1000, temperature decreases only slightly, whereas ocean heat content and sea level continue to rise. Our best estimate for the Absolute Global Warming Potential, given by the time-integrated response in CO2 at year 100 multiplied by its radiative efficiency, is 92.5 × 10−15 yr W m−2 per kg-CO2. This value very likely (5 to 95% confidence) lies within the range of (68 to 117) × 10−15 yr W m−2 per kg-CO2. Estimates for time-integrated response in CO2 published in the IPCC First, Second, and Fourth Assessment and our multi-model best estimate all agree within 15% during the first 100 yr. The integrated CO2 response, normalized by the pulse size, is lower for pre-industrial conditions, compared to present day, and lower for smaller pulses than larger pulses. In contrast, the response in temperature, sea level and ocean heat content is less sensitive to these choices. Although, choices in pulse size, background concentration, and model lead to uncertainties, the most important and subjective choice to determine AGWP of CO2 and GWP is the time horizon.

Amplified inception of European Little Ice Age by sea ice-ocean-atmosphere feedbacks

The inception of the Little Ice Age (~1400–1700 AD) is believed to have been driven by an interplay of external forcing and climate system internal variability. While the hemispheric signal seems to have been dominated by solar irradiance and volcanic eruptions, the understanding of mechanisms shaping the climate on a continental scale is less robust. In an ensemble of transient model simulations and a new type of sensitivity experiments with artificial sea ice growth, the authors identify a sea ice–ocean–atmosphere feedback mechanism that amplifies the Little Ice Age cooling in the North Atlantic–European region and produces the temperature pattern suggested by paleoclimatic reconstructions. Initiated by increasing negative forcing, the Arctic sea ice substantially expands at the beginning of the Little Ice Age. The excess of sea ice is exported to the subpolar North Atlantic, where it melts, thereby weakening convection of the ocean. Consequently, northward ocean heat transport is reduced, reinforcing the expansion of the sea ice and the cooling of the Northern Hemisphere. In the Nordic Seas, sea surface height anomalies cause the oceanic recirculation to strengthen at the expense of the warm Barents Sea inflow, thereby further reinforcing sea ice growth. The absent ocean–atmosphere heat flux in the Barents Sea results in an amplified cooling over Northern Europe. The positive nature of this feedback mechanism enables sea ice to remain in an expanded state for decades up to a century, favoring sustained cold periods over Europe such as the Little Ice Age. Support for the feedback mechanism comes from recent proxy reconstructions around the Nordic Seas.

A reconstruction of radiocarbon production and total solar irradiance from the Holocene ¹⁴C and CO₂ records: implications of data and model uncertainties

Radiocarbon production, solar activity, total solar irradiance (TSI) and solar-induced climate change are reconstructed for the Holocene (10 to 0 kyr BP), and TSI is predicted for the next centuries. The IntCal09/SHCal04 radiocarbon and ice core CO2 records, reconstructions of the geomagnetic dipole, and instrumental data of solar activity are applied in the Bern3D-LPJ, a fully featured Earth system model of intermediate complexity including a 3-D dynamic ocean, ocean sediments, and a dynamic vegetation model, and in formulations linking radiocarbon production, the solar modulation potential, and TSI. Uncertainties are assessed using Monte Carlo simulations and bounding scenarios. Transient climate simulations span the past 21 thousand years, thereby considering the time lags and uncertainties associated with the last glacial termination. Our carbon-cycle-based modern estimate of radiocarbon production of 1.7 atoms cm−2 s−1 is lower than previously reported for the cosmogenic nuclide production model by Masarik and Beer (2009) and is more in-line with Kovaltsov et al. (2012). In contrast to earlier studies, periods of high solar activity were quite common not only in recent millennia, but throughout the Holocene. Notable deviations compared to earlier reconstructions are also found on decadal to centennial timescales. We show that earlier Holocene reconstructions, not accounting for the interhemispheric gradients in radiocarbon, are biased low. Solar activity is during 28% of the time higher than the modern average (650 MeV), but the absolute values remain weakly constrained due to uncertainties in the normalisation of the solar modulation to instrumental data. A recently published solar activity–TSI relationship yields small changes in Holocene TSI of the order of 1 W m−2 with a Maunder Minimum irradiance reduction of 0.85 ± 0.16 W m−2. Related solar-induced variations in global mean surface air temperature are simulated to be within 0.1 K. Autoregressive modelling suggests a declining trend of solar activity in the 21st century towards average Holocene conditions.

Upper ocean vertical supply: A neglected primary factor controlling the distribution of neodymium concentrations of open ocean surface waters?

Neodymium (Nd) isotopes are an important geochemical tool to trace the present and past water mass mixing as well as continental inputs. The distribution of Nd concentrations in open ocean surface waters (0�100 m) is generally assumed to be controlled by lateral mixing of Nd from coastal surface currents and by removal through reversible particle scavenging. However, using 228Ra activity as an indicator of coastal water mass influence, surface water Nd concentration data available on key oceanic transects as a whole do not support the above scenario. From a global compilation of available data, we find that more stratified regions are generally associated with low surface Nd concentrations. This implies that upper ocean vertical supply may be an as yet neglected primary factor in determining the basin-scale variations of surface water Nd concentrations. Similar to the mechanism of nutrients supply, it is likely that stratification inhibits vertical supply of Nd from the subsurface thermocline waters and thus the magnitude of Nd flux to the surface layer. Consistently, the estimated required input flux of Nd to the surface layer to maintain the observed concentrations could be nearly two orders of magnitudes larger than riverine/dust flux, and also larger than the model-based estimation on shelf-derived coastal flux. In addition, preliminary results from modeling experiments reveal that the input from shallow boundary sources, riverine input, and release from dust are actually not the primary factors controlling Nd concentrations most notably in the Pacific and Southern Ocean surface waters.

Paleoclimate data-model comparison and the role of climate forcings over the past 1500 years

The past 1500 years provide a valuable opportunity to study the response of the climate system to external forcings. However, the integration of paleoclimate proxies with climate modeling is critical to improving the understanding of climate dynamics. In this paper, a climate system model and proxy records are therefore used to study the role of natural and anthropogenic forcings in driving the global climate. The inverse and forward approaches to paleoclimate data–model comparison are applied, and sources of uncertainty are identified and discussed. In the first of two case studies, the climate model simulations are compared with multiproxy temperature reconstructions. Robust solar and volcanic signals are detected in Southern Hemisphere temperatures, with a possible volcanic signal detected in the Northern Hemisphere. The anthropogenic signal dominates during the industrial period. It is also found that seasonal and geographical biases may cause multiproxy reconstructions to overestimate the magnitude of the long-term preindustrial cooling trend. In the second case study, the model simulations are compared with a coral δ18O record from the central Pacific Ocean. It is found that greenhouse gases, solar irradiance, and volcanic eruptions all influence the mean state of the central Pacific, but there is no evidence that natural or anthropogenic forcings have any systematic impact on El Niño–Southern Oscillation. The proxy climate relationship is found to change over time, challenging the assumption of stationarity that underlies the interpretation of paleoclimate proxies. These case studies demonstrate the value of paleoclimate data–model comparison but also highlight the limitations of current techniques and demonstrate the need to develop alternative approaches.

Global ocean storage of anthropogenic carbon

The global ocean is a significant sink for anthropogenic carbon (Cant), absorbing roughly a third of human CO2 emitted over the industrial period. Robust estimates of the magnitude and variability of the storage and distribution of Cant in the ocean are therefore important for understanding the human impact on climate. In this synthesis we review observational and model-based estimates of the storage and transport of Cant in the ocean. We pay particular attention to the uncertainties and potential biases inherent in different inference schemes. On a global scale, three data-based estimates of the distribution and inventory of Cant are now available. While the inventories are found to agree within their uncertainty, there are considerable differences in the spatial distribution. We also present a review of the progress made in the application of inverse and data assimilation techniques which combine ocean interior estimates of Cant with numerical ocean circulation models. Such methods are especially useful for estimating the air–sea flux and interior transport of Cant, quantities that are otherwise difficult to observe directly. However, the results are found to be highly dependent on modeled circulation, with the spread due to different ocean models at least as large as that from the different observational methods used to estimate Cant. Our review also highlights the importance of repeat measurements of hydrographic and biogeochemical parameters to estimate the storage of Cant on decadal timescales in the presence of the variability in circulation that is neglected by other approaches. Data-based Cant estimates provide important constraints on forward ocean models, which exhibit both broad similarities and regional errors relative to the observational fields. A compilation of inventories of Cant gives us a "best" estimate of the global ocean inventory of anthropogenic carbon in 2010 of 155 ± 31 PgC (±20% uncertainty). This estimate includes a broad range of values, suggesting that a combination of approaches is necessary in order to achieve a robust quantification of the ocean sink of anthropogenic CO2.

State of the Antarctic and Southern Ocean Climate System

This paper reviews developments in our understanding of the state of the Antarctic and Southern Ocean climate and its relation to the global climate system over the last few millennia. Climate over this and earlier periods has not been stable, as evidenced by the occurrence of abrupt changes in atmospheric circulation and temperature recorded in Antarctic ice core proxies for past climate. Two of the most prominent abrupt climate change events are characterized by intensification of the circumpolar westerlies (also known as the Southern Annular Mode) between similar to 6000 and 5000 years ago and since 1200-1000 years ago. Following the last of these is a period of major trans-Antarctic reorganization of atmospheric circulation and temperature between A. D. 1700 and 1850. The two earlier Antarctic abrupt climate change events appear linked to but predate by several centuries even more abrupt climate change in the North Atlantic, and the end of the more recent event is coincident with reorganization of atmospheric circulation in the North Pacific. Improved understanding of such events and of the associations between abrupt climate change events recorded in both hemispheres is critical to predicting the impact and timing of future abrupt climate change events potentially forced by anthropogenic changes in greenhouse gases and aerosols. Special attention is given to the climate of the past 200 years, which was recorded by a network of recently available shallow firn cores, and to that of the past 50 years, which was monitored by the continuous instrumental record. Significant regional climate changes have taken place in the Antarctic during the past 50 years. Atmospheric temperatures have increased markedly over the Antarctic Peninsula, linked to nearby ocean warming and intensification of the circumpolar westerlies. Glaciers are retreating on the peninsula, in Patagonia, on the sub-Antarctic islands, and in West Antarctica adjacent to the peninsula. The penetration of marine air masses has become more pronounced over parts of West Antarctica. Above the surface, the Antarctic troposphere has warmed during winter while the stratosphere has cooled year-round. The upper kilometer of the circumpolar Southern Ocean has warmed, Antarctic Bottom Water across a wide sector off East Antarctica has freshened, and the densest bottom water in the Weddell Sea has warmed. In contrast to these regional climate changes, over most of Antarctica, near-surface temperature and snowfall have not increased significantly during at least the past 50 years, and proxy data suggest that the atmospheric circulation over the interior has remained in a similar state for at least the past 200 years. Furthermore, the total sea ice cover around Antarctica has exhibited no significant overall change since reliable satellite monitoring began in the late 1970s, despite large but compensating regional changes. The inhomogeneity of Antarctic climate in space and time implies that recent Antarctic climate changes are due on the one hand to a combination of strong multidecadal variability and anthropogenic effects and, as demonstrated by the paleoclimate record, on the other hand to multidecadal to millennial scale and longer natural variability forced through changes in orbital insolation, greenhouse gases, solar variability, ice dynamics, and aerosols. Model projections suggest that over the 21st century the Antarctic interior will warm by 3.4 degrees +/- 1 degrees C, and sea ice extent will decrease by similar to 30%. Ice sheet models are not yet adequate enough to answer pressing questins about the effect of projected warming on mass balance and sea level. Considering the potentially major impacts of a warming climate on Antarctica, vigorous efforts are needed to better understand all aspects of the highly coupled Antarctic climate system as well as its influence on the Earth's climate and oceans.

Model-Based Estimates of Calanus Finmarchicus Abundance in the Gulf of Maine

Ocean observing systems and satellites routinely collect a wealth of information on physical conditions in the ocean. With few exceptions, such as chlorophyll concentrations, information on biological properties is harder to measure autonomously. Here, we present a system to produce estimates of the distribution and abundance of the copepod Calanus finmarchicus in the Gulf of Maine. Our system uses satellite-based measurements of sea surface temperature and chlorophyll concentration to determine the developmental and reproductive rates of C. finmarchicus. The rate information then drives a population dynamics model of C. finmarchicus that is embedded in a 2-dimensional circulation field. The first generation of this system produces realistic information on interannual variability in C. finmarchicus distribution and abundance during the winter and spring. The model can also be used to identify key drivers of interannual variability in C. finmarchicus. Experiments with the model suggest that changes in initial conditions are overwhelmed by variability in growth rates after approximately 50 d. Temperature has the largest effect on growth rate. Elevated chlorophyll during the late winter can lead to increased C. finmarchicus abundance during the spring, but the effect of variations in chlorophyll concentrations is secondary to the other inputs. Our system could be used to provide real-time estimates or even forecasts of C. finmarchicus distribution. These estimates could then be used to support management of copepod predators such as herring and right whales.

Modeling Responses of Diatom Productivity and Biogenic Silica Export to Iron Enrichment in the Equatorial Pacific Ocean

Using a three-dimensional physical-biogeochemical model, we have investigated the modeled responses of diatom productivity and biogenic silica export to iron enrichment in the equatorial Pacific, and compared the model simulation with in situ (IronEx II) iron fertilization results. In the eastern equatorial Pacific, an area of 540,000 km(2) was enhanced with iron by changing the photosynthetic efficiency and silicate and nitrogen uptake kinetics of phytoplankton in the model for a period of 20 days. The vertically integrated Chl a and primary production increased by about threefold 5 days after the start of the experiment, similar to that observed in the IronEx II experiment. Diatoms contribute to the initial increase of the total phytoplankton biomass, but decrease sharply after 10 days because of mesozooplankton grazing. The modeled surface nutrients (silicate and nitrate) and TCO(2) anomaly fields, obtained from the difference between the "iron addition'' and "ambient'' (without iron) concentrations, also agreed well with the IronEx II observations. The enriched patch is tracked with an inert tracer similar to the SF6 used in the IronEx II. The modeled depth-time distribution of sinking biogenic silica (BSi) indicates that it would take more than 30 days after iron injection to detect any significant BSi export out of the euphotic zone. Sensitivity studies were performed to establish the importance of fertilized patch size, duration of fertilization, and the role of mesozooplankton grazing. A larger size of the iron patch tends to produce a broader extent and longer-lasting phytoplankton blooms. Longer duration prolongs phytoplankton growth, but higher zooplankton grazing pressure prevents significant phytoplankton biomass accumulation. With the same treatment of iron fertilization in the model, lowering mesozooplankton grazing rate generates much stronger diatom bloom, but it is terminated by Si(OH)(4) limitation after the initial rapid increase. Increasing mesozooplankton grazing rate, the diatom increase due to iron addition stays at minimum level, but small phytoplankton tend to increase. The numerical model experiments demonstrate the value of ecosystem modeling for evaluating the detailed interaction between biogeochemical cycle and iron fertilization in the equatorial Pacific.

Decadal Variability in the Northeast Pacific in a Physical-Ecosystem Model: Role of Mixed Layer Depth and Trophic Interactions

A basin-wide interdecadal change in both the physical state and the ecology of the North Pacific occurred near the end of 1976. Here we use a physical-ecosystem model to examine whether changes in the physical environment associated with the 1976-1977 transition influenced the lower trophic levels of the food web and if so by what means. The physical component is an ocean general circulation model, while the biological component contains 10 compartments: two phytoplankton, two zooplankton, two detritus pools, nitrate, ammonium, silicate, and carbon dioxide. The model is forced with observed atmospheric fields during 1960-1999. During spring, there is a similar to 40% reduction in plankton biomass in all four plankton groups during 1977-1988 relative to 1970-1976 in the central Gulf of Alaska (GOA). The epoch difference in plankton appears to be controlled by the mixed layer depth. Enhanced Ekman pumping after 1976 caused the halocline to shoal, and thus the mixed layer depth, which extends to the top of the halocline in late winter, did not penetrate as deep in the central GOA. As a result, more phytoplankton remained in the euphotic zone, and phytoplankton biomass began to increase earlier in the year after the 1976 transition. Zooplankton biomass also increased, but then grazing pressure led to a strong decrease in phytoplankton by April followed by a drop in zooplankton by May: Essentially, the mean seasonal cycle of plankton biomass was shifted earlier in the year. As the seasonal cycle progressed, the difference in plankton concentrations between epochs reversed sign again, leading to slightly greater zooplankton biomass during summer in the later epoch.

Time of emergence of trends in ocean biogeochemistry

For the detection of climate change, not only the magnitude of a trend signal is of significance. An essential issue is the time period required by the trend to be detectable in the first place. An illustrative measure for this is time of emergence (ToE), that is, the point in time when a signal finally emerges from the background noise of natural variability. We investigate the ToE of trend signals in different biogeochemical and physical surface variables utilizing a multi-model ensemble comprising simulations of 17 Earth system models (ESMs). We find that signals in ocean biogeochemical variables emerge on much shorter timescales than the physical variable sea surface temperature (SST). The ToE patterns of pCO2 and pH are spatially very similar to DIC (dissolved inorganic carbon), yet the trends emerge much faster – after roughly 12 yr for the majority of the global ocean area, compared to between 10 and 30 yr for DIC. ToE of 45–90 yr are even larger for SST. In general, the background noise is of higher importance in determining ToE than the strength of the trend signal. In areas with high natural variability, even strong trends both in the physical climate and carbon cycle system are masked by variability over decadal timescales. In contrast to the trend, natural variability is affected by the seasonal cycle. This has important implications for observations, since it implies that intra-annual variability could question the representativeness of irregularly sampled seasonal measurements for the entire year and, thus, the interpretation of observed trends.

Ocean (De)oxygenation Across the Last Deglaciation: Insights for the Future

Anthropogenic warming is expected to drive oxygen out of the ocean as the water temperature rises and the rate of exchange between subsurface waters and the atmosphere slows due to enhanced upper ocean density stratification. Observations from recent decades are tantalizingly consistent with this prediction, though these changes remain subtle in the face of natural variability. Earth system model projections unanimously predict a long-term decrease in the global ocean oxygen inventory, but show regional discrepancies, particularly in the most oxygen-depleted waters, owing to the complex interplay between oxygen supply pathways and oxygen consumption. The geological record provides an orthogonal perspective, showing how the oceanic oxygen content varied in response to prior episodes of climate change. These past changes were much slower than the current, anthropogenic change, but can help to appraise sensitivities, and point toward potentially dominant mechanisms of change. Consistent with the model projections, marine sediments recorded an overall expansion of low-oxygen waters in the upper ocean as it warmed at the end of the last ice age. This expansion was not linearly related with temperature, though, but reached a deoxygenation extreme midway through the warming. Meanwhile, the deep ocean became better oxygenated, opposite the general expectation. These observations require that significant changes in apparent oxygen utilization occurred, suggesting that they will also be important in the future.