Abstract Astronomical cycles are strongly expressed in marine geological records, providing important insights into Earth system dynamics and an invaluable means of constructing age models. However, how various astronomical periods are filtered by the Earth system and the mechanisms by which carbon reservoirs and climate components respond, particularly in absence of dynamic ice sheets, is unclear. Using an Earth system model that includes feedbacks between climate, ocean circulation, and inorganic (carbonate) carbon cycling relevant to geological timescales, we systematically explore the impact of astronomically‐modulated insolation forcing and its expression in model variables most comparable to key paleoceanographic proxies (temperature, the δ13C of inorganic carbon, and sedimentary carbonate content). Temperature predominately responds to obliquity and is little influenced by the modeled carbon cycle feedbacks. In contrast, the cycling of nutrients and carbon in the ocean generates significant precession power in atmospheric CO2, benthic ocean δ13C, and sedimentary wt% CaCO3, while inclusion of marine sedimentary and weathering processes shifts power to the long eccentricity period. Our simulations produce reducedpCO2and dissolved inorganic carbon δ13C at long eccentricity maxima and, contrary to early Cenozoic marine records, CaCO3preservation in the model is enhanced during eccentricity modulated warmth. Additionally, the magnitude of δ13C variability simulated in our model underestimates marine proxy records. These model‐data discrepancies hint at the possibility that the Paleogene silicate weathering feedback was weaker than modeled here and that additional organic carbon cycle feedbacks are necessary to explain the full response of the Earth system to astronomical forcing.
more »
« less
This content will become publicly available on April 25, 2026
Obliquity disruption and Antarctic ice sheet dynamics over a 2.4-Myr astronomical grand cycle
Marine δ18O data reveal astronomical forcing of the climate and cryosphere during the Miocene, when atmosphericPco2was on par with emissions scenarios over the next century. This inspired hypotheses for how Milankovitch cycles, ice-ocean interactions, and greenhouse gases influence ice volume. Mass balance controls for marine and terrestrial ice sheets differ, and proxy data collected far from Antarctica provide valuable but limited insight into regional processes. We evaluate clast abundance data from Antarctic marine sedimentary records, observing a strong signal of eccentricity and precession coincident with a terrestrial ice sheet and a clear obliquity signal at the margins of a marine ice sheet. These analyses are integrated with a synthesis of proxy data, and we argue that high variance in obliquity forcing (mediated and enhanced by the ocean and atmosphere) can inhibit ice sheet growth, even when insolation forcing is conducive to glaciation. This “obliquity disruption” explains cryosphere variability before the existence of large northern hemisphere ice sheets.
more »
« less
- Award ID(s):
- 1947646
- PAR ID:
- 10628133
- Publisher / Repository:
- AAAS
- Date Published:
- Journal Name:
- Science Advances
- Volume:
- 11
- Issue:
- 17
- ISSN:
- 2375-2548
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Astronomical cycles are strongly expressed in marine geological records, providing important insights into Earth system dynamics and an invaluable means of constructing age models. However, how various astronomical periods are filtered by the Earth system and the mechanisms by which carbon reservoirs and climate components respond, particularly in absence of dynamic ice sheets, is unclear. Using an Earth system model that includes feedbacks between climate, ocean circulation, and inorganic (carbonate) carbon cycling relevant to geological timescales, we systematically explore the impact of astronomically modulated insolation forcing and its expression in model variables most comparable to key paleoceanographic proxies (temperature, the δ13C of inorganic carbon, and sedimentary carbonate content). Temperature predominately responds to short and long eccentricity and is little influenced by the modeled carbon cycle feedbacks. In contrast, the cycling of nutrients and carbon in the ocean generates significant precession power in atmospheric CO2, benthic ocean δ13C, and sedimentary wt% CaCO3, while inclusion of marine sedimentary and weathering processes shifts power to the long eccentricity period. Our simulations produce reducedpCO2and dissolved inorganic carbon (DIC) δ13C at long eccentricity maxima and, contrary to early Cenozoic marine records, CaCO3preservation in the model is enhanced during eccentricity‐modulated warmth. Additionally, the magnitude of δ13C variability simulated in our model underestimates marine proxy records. These model‐data discrepancies hint at the possibility that the Paleogene silicate weathering feedback was weaker than modeled here and that additional organic carbon cycle feedbacks are necessary to explain the full response of the Earth system to astronomical forcing.more » « less
-
Abstract Speleothem CaCO3δ18O is a commonly employed paleomonsoon proxy. However, inferring local rainfall amount from speleothem δ18O can be complicated due to changing source water δ18O, temperature effects, and rainout over the moisture transport path. These complications are addressed using δ18O of planktonic foraminiferal CaCO3, offshore from the Yangtze River Valley (YRV). The advantage is that the effects of global seawater δ18O and local temperature changes can be quantitatively removed, yielding a record of local seawater δ18O, a proxy that responds primarily to dilution by local precipitation and runoff. Whereas YRV speleothem δ18O is dominated by precession-band (23 ky) cyclicity, local seawater δ18O is dominated by eccentricity (100 ky) and obliquity (41 ky) cycles, with almost no precession-scale variance. These results, consistent with records outside the YRV, suggest that East Asian monsoon rainfall is more sensitive to greenhouse gas and high-latitude ice sheet forcing than to direct insolation forcing.more » « less
-
Changes in global mean sea level (GMSL) during the late Cenozoic remain uncertain. We use a reconstruction of changes in δ18O of seawater to reconstruct GMSL since 4.5 million years ago (Ma) that accounts for temperature-driven changes in the δ18O of global ice sheets. Between 4.5 and 3 Ma, sea level highstands remained up to 20 m above present whereas the first lowstands below present suggest onset of Northern Hemisphere glaciation at 4 Ma. Intensification of global glaciation occurred from 3 Ma to 2.5 Ma, culminating in lowstands similar to the Last Glacial Maximum lowstand at 21,000 years ago and that reoccurred throughout much of the Pleistocene. We attribute the middle Pleistocene transition in ice sheet variability (1.2 Ma to 0.62 Ma) to modulation of 41-thousand-year (kyr) obliquity forcing by an increase in ~100-kyr CO2variability.more » « less
-
Abstract Over the last 3.3 million years, the Antarctic Ice Sheet (AIS) has undergone phases of ice sheet growth and decay, impacting sea level and climate globally. Presently, the largely marine‐terminating AIS loses mass primarily by iceberg calving and basal melt of ice shelves. Quantifying past rates and timing of AIS melt is vital to understanding future cryosphere and sea level changes. One proxy for past ice sheet instabilities is iceberg rafted debris (IRD) fluxes. However, traditional methods of IRD quantification are labor‐intensive. Here, we present a new method of identifying IRD grains in sediment core X‐ray images using a convolutional neural network machine learning algorithm. We present a 3.3‐million‐year record of AIS IRD melt events using sediment cores from International Ocean Discovery Program Sites U1536, U1537, and U1538 in the Southern Ocean's “Iceberg Alley.” We identify two increases in the IRD fluxes throughout this period, at ∼1.8 and 0.43 Ma. We propose that after 1.8 Ma, the AIS expanded and transitioned from a primarily terrestrial‐terminating to a primarily marine‐terminating ice sheet. Therefore, after 1.8 Ma, glacial terminations and AIS iceberg discharge are associated with variations in global ice volume, presumably through the mechanism of sea level and, therefore, grounding line change. The second AIS regime change occurs during the Mid‐Brunhes Event (∼0.43 Ma). After this time, there are heightened and continuous IRD fluxes at each glacial termination, indicating increased AIS size and instability after this time.more » « less
An official website of the United States government
