skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Global Ocean Cooling of 2.3°C During the Last Glacial Maximum
Abstract Quantitative constraints on past mean ocean temperature (MOT) critically inform our historical understanding of Earth's energy balance. A recently developed MOT proxy based on paleoatmospheric Xe, Kr, and N2ratios in ice core air bubbles is a promising tool rooted in the temperature dependences of gas solubilities. However, these inert gases are systematically undersaturated in the modern ocean interior, and it remains unclear how air‐sea disequilibrium may have changed in the past. Here, we carry out 30 tracer‐enabled model simulations under varying circulation, sea ice cover, and wind stress regimes to evaluate air‐sea disequilibrium in the Last Glacial Maximum (LGM) ocean. We find that undersaturation of all three gases was likely reduced, primarily due to strengthened high‐latitude winds, biasing reconstructed MOT by −0.38 ± 0.37°C (1σ). Accounting for air‐sea disequilibrium, paleoatmospheric inert gases indicate that LGM MOT was 2.27 ± 0.46°C (1σ) colder than the pre‐industrial era.  more » « less
Award ID(s):
2049359 2048926
PAR ID:
10506115
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
51
Issue:
9
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Proxy‐based reconstructions suggest that equilibrium changes in global mean sea surface temperature (ΔGMSST) are nearly equivalent to changes in mean ocean temperature (ΔMOT) on glacial‐interglacial timescales over the past 900,000 years. However, the underlying mechanisms responsible for this relationship remain poorly understood. Here we use simulations from Paleoclimate Modeling Intercomparison Project Phase 3 and 4 (PMIP3/4) to investigate equilibrium ΔMOT and its linkage to sea surface temperature changes between the Last Glacial Maximum (LGM, ∼21,000 years ago) and pre‐Industrial. Results show that PMIP3/4 simulations generally underestimate proxy‐based ΔMOT. Regression analysis reveals that LGM MOT is strongly modulated by mid‐latitude SST cooling, with the Southern Ocean having a greater influence compared to other oceanic regions, thus helping explain why models with similar ΔGMSSTs exhibit significantly different ΔMOTs. Additionally, we find a strong relationship between simulated Antarctic sea‐ice coverage and Southern Ocean SST changes, with implications for constraining sea‐ice reconstructions. 
    more » « less
  2. Abstract Southern Ocean sea ice plays a central role in the oceanic meridional overturning circulation, transforming globally prevalent watermasses through surface buoyancy loss and gain. Buoyancy loss due to surface cooling and sea ice growth promotes the formation of bottom water that flows into the Atlantic, Indian, and Pacific basins, while buoyancy gain due to sea ice melt helps transform the returning deep flow into intermediate and mode waters. Because northward expansion of Southern Ocean sea ice during the Last Glacial Maximum (LGM; 19–23 kyr BP) may have enhanced deep ocean stratification and contributed to lower atmospheric CO2levels, reconstructions of sea ice extent are critical to understanding the LGM climate state. Here, we present a new sea ice proxy based on the18O/16O ratio of foraminifera (δ18Oc). In the seasonal sea ice zone, sea ice formation during austral winter creates a cold surface mixed layer that persists in the sub‐surface during spring and summer. The cold sub‐surface layer, known as winter water, sits above relatively warm deep water, creating an inverted temperature profile. The unique surface‐to‐deep temperature contrast is reflected in estimates of equilibrium δ18Oc, implying that paired analysis of planktonic and benthic foraminifera can be used to infer sea ice extent. To demonstrate the feasibility of the δ18Ocmethod, we present a compilation ofN. pachydermaandCibicidoidesspp. results from the Atlantic sector that yields an estimate of winter sea ice extent consistent with modern observations. 
    more » « less
  3. Abstract Since the last glacial period, North America has experienced dramatic changes in regional climate, including the collapse of ice sheets and changes in precipitation. We use clumped isotope (∆47) thermometry and carbonate δ18O measurements of glacial and deglacial pedogenic carbonates from the Palouse Loess to provide constraints on hydroclimate changes in the Pacific Northwest. We also employ analysis of climate model simulations to help us further provide constraints on the hydroclimate changes in the Pacific Northwest. The coldest clumped isotope soil temperaturesT(47) (13.5 ± 1.9°C to 17.1 ± 1.7°C) occurred ∼34,000–23,000 years ago. Using a soil‐to‐air temperature transfer function, we estimate Last Glacial Maximum (LGM) mean annual air temperatures of ∼−5.5°C and warmest average monthly temperatures (i.e., mean summer air temperatures) of ∼4.4°C. These data indicate a regional warming of 16.4 ± 2.6°C from the LGM to the modern temperatures of 10.9°C, which was about 2.5–3 times the global average. Proxy data provide locality constraints on the boundary of the cooler anticyclone induced by LGM ice sheets, and the warmer cyclone in the Eastern Pacific Ocean. Climate model analysis suggests regional amplification of temperature anomalies is due to the proximal location of the study area to the Laurentide Ice Sheet margin and the impact of the glacial anticyclone on the region, as well as local albedo. Isotope‐enabled model experiments indicate variations in water δ18O largely reflect atmospheric circulation changes and enhanced rainout upstream that brings more depleted vapor to the region during the LGM. 
    more » « less
  4. null (Ed.)
    Abstract. Equilibrium climate sensitivity (ECS) has been directly estimated using reconstructions of past climates that are different than today's. A challenge to this approach is that temperature proxies integrate over the timescales of the fast feedback processes (e.g., changes in water vapor, snow, and clouds) that are captured in ECS as well as the slower feedback processes (e.g., changes in ice sheets and ocean circulation) that are not. A way around this issue is to treat the slow feedbacks as climate forcings and independently account for their impact on global temperature. Here we conduct a suite of Last Glacial Maximum (LGM) simulations using the Community Earth System Model version 1.2 (CESM1.2) to quantify the forcingand efficacy of land ice sheets (LISs) and greenhouse gases (GHGs) in order to estimate ECS. Our forcing and efficacy quantification adopts the effective radiative forcing (ERF) and adjustment framework and provides a complete accounting for the radiative, topographic, and dynamical impacts of LIS on surface temperatures. ERF and efficacy of LGM LIS are −3.2 W m−2 and 1.1, respectively. The larger-than-unity efficacy is caused by the temperature changes over land and the Northern Hemisphere subtropical oceans which are relatively larger than those in response to a doubling of atmospheric CO2. The subtropical sea-surface temperature (SST) response is linked to LIS-induced wind changes and feedbacks in ocean–atmosphere coupling and clouds. ERF and efficacy of LGM GHG are −2.8 W m−2 and 0.9, respectively. The lower efficacy is primarily attributed to a smaller cloud feedback at colder temperatures. Our simulations further demonstrate that the direct ECS calculation using the forcing, efficacy, and temperature response in CESM1.2 overestimates the true value in the model by approximately 25 % due to the neglect of slow ocean dynamical feedback. This is supported by the greater cooling (6.8 ∘C) in a fully coupled LGM simulation than that (5.3 ∘C) in a slab ocean model simulation with ocean dynamics disabled. The majority (67 %) of the ocean dynamical feedback is attributed to dynamical changes in the Southern Ocean, where interactions between upper-ocean stratification, heat transport, and sea-ice cover are found to amplify the LGM cooling. Our study demonstrates the value of climate models in the quantification of climate forcings and the ocean dynamical feedback, which is necessary for an accurate direct ECS estimation. 
    more » « less
  5. Carré, Matthieu (Ed.)
    Despite their importance for Earth’s climate and paleoceanography, the cycles of carbon (C) and its isotope13C in the ocean are not well understood. Models typically do not decompose C and13C storage caused by different physical, biological, and chemical processes, which makes interpreting results difficult. Consequently, basic observed features, such as the decreased carbon isotopic signature (δ13CDIC) of the glacial ocean remain unexplained. Here, we review recent progress in decomposing Dissolved Inorganic Carbon (DIC) into preformed and regenerated components, extend a precise and complete decomposition to δ13CDIC, and apply it to data-constrained model simulations of the Preindustrial (PI) and Last Glacial Maximum (LGM) oceans. Regenerated components, from respired soft-tissue organic matter and dissolved biogenic calcium carbonate, are reduced in the LGM, indicating a decrease in the active part of the biological pump. Preformed components increase carbon storage and decrease δ13CDICby 0.55 ‰ in the LGM. We separate preformed into saturation and disequilibrium components, each of which have biological and physical contributions. Whereas the physical disequilibrium in the PI is negative for both DIC and δ13CDIC, and changes little between climate states, the biological disequilibrium is positive for DIC but negative for δ13CDIC, a pattern that is magnified in the LGM. The biological disequilibrium is the dominant driver of the increase in glacial ocean C and the decrease in δ13CDIC, indicating a reduced sink of biological carbon. Overall, in the LGM, biological processes increase the ocean’s DIC inventory by 355 Pg more than in the PI, reduce its mean δ13CDICby an additional 0.52 ‰, and contribute 60 ppm to the lowering of atmospheric CO2. Spatial distributions of the δ13CDICcomponents are presented. Commonly used approximations based on apparent oxygen utilization and phosphate are evaluated and shown to have large errors. 
    more » « less