Extraction of GRACE/GRACE-FO observed mass change patterns across Antarctica via independent component analysis (ICA).

Here we qualitatively analyse the mass change patterns across Antarctica via independent component analysis (ICA), a statistics-based blind source separation method to extract signals from complex data sets, in an attempt to reduce uncertainties in the glacial isostatic adjustment (GIA) effects and improve understanding of Antarctic Ice Sheet (AIS) mass-balance. We extract the six leading independent components from gravimetric data acquired during the Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow-On (GRACE-FO) missions. The results reveal that the observed continental-scale mass changes can be effectively separated into several spatial patterns that may be dominated by different physical processes. Although the hidden independent physical processes cannot be completely isolated, some significant signals, such as glacier melt, snow accumulation, periodic climatic signals, and GIA effects, can be determined without introducing any external information. We also observe that the time period of the analysed data set has a direct impact on the ICA results, as the impacts of extreme events, such as the anomalously large snowfall events in the late 2000s, may cause dramatic spatial and temporal changes in the ICA results. ICA provides a unique and informative approach to obtain a better understanding of both AIS-scale mass changes and specific regional-scale spatiotemporal signal variations.

A sea-level plateau preceding the Marine Isotope Stage 2 minima revealed by Australian sediments.

Further understanding of past climate requires a robust estimate of global ice volume fluctuations that in turn rely on accurate global sea-level reconstructions. An advantage of Marine Isotope Stage 2 (MIS 2) is the availability of suitable material for radiocarbon dating to allow comparison of sea-level data with other paleoclimatic proxies. However, the number and accuracy of sea-level records during MIS 2 is currently lacking. Here we present the history of MIS 2 eustatic sea-level change as recorded in the Bonaparte Gulf, northwestern Australia by reconstructing relative sea level and then modeling glacial isostatic adjustment. The isostatically-corrected global sea-level history indicates that sea-level plateaued from 25.9 to 20.4 cal kyr BP (modeled median probability) prior reaching its minimum (19.7 to 19.1 cal kyr BP). Following the plateau, we detect a 10-m global sea-level fall over ~1,000 years and a short duration of the Last Glacial Maximum (global sea-level minimum; 19.7 to 19.1 cal kyr BP). These large changes in ice volume over such a short time indicates that the continental ice sheets never reached their isostatic equilibrium during the Last Glacial Maximum.

Rapid glaciation and a two-step sea level plunge into the Last Glacial Maximum.

The approximately 10,000-year-long Last Glacial Maximum, before the termination of the last ice age, was the coldest period in Earth's recent climate history1. Relative to the Holocene epoch, atmospheric carbon dioxide was about 100 parts per million lower and tropical sea surface temperatures were about 3 to 5 degrees Celsius lower2,3. The Last Glacial Maximum began when global mean sea level (GMSL) abruptly dropped by about 40 metres around 31,000 years ago4 and was followed by about 10,000 years of rapid deglaciation into the Holocene1. The masses of the melting polar ice sheets and the change in ocean volume, and hence in GMSL, are primary constraints for climate models constructed to describe the transition between the Last Glacial Maximum and the Holocene, and future changes; but the rate, timing and magnitude of this transition remain uncertain. Here we show that sea level at the shelf edge of the Great Barrier Reef dropped by around 20 metres between 21,900 and 20,500 years ago, to -118 metres relative to the modern level. Our findings are based on recovered and radiometrically dated fossil corals and coralline algae assemblages, and represent relative sea level at the Great Barrier Reef, rather than GMSL. Subsequently, relative sea level rose at a rate of about 3.5 millimetres per year for around 4,000 years. The rise is consistent with the warming previously observed at 19,000 years ago1,5, but we now show that it occurred just after the 20-metre drop in relative sea level and the related increase in global ice volumes. The detailed structure of our record is robust because the Great Barrier Reef is remote from former ice sheets and tectonic activity. Relative sea level can be influenced by Earth's response to regional changes in ice and water loadings and may differ greatly from GMSL. Consequently, we used glacio-isostatic models to derive GMSL, and find that the Last Glacial Maximum culminated 20,500 years ago in a GMSL low of about -125 to -130 metres.

State dependence of climatic instability over the past 720,000 years from Antarctic ice cores and climate modeling.

Climatic variabilities on millennial and longer time scales with a bipolar seesaw pattern have been documented in paleoclimatic records, but their frequencies, relationships with mean climatic state, and mechanisms remain unclear. Understanding the processes and sensitivities that underlie these changes will underpin better understanding of the climate system and projections of its future change. We investigate the long-term characteristics of climatic variability using a new ice-core record from Dome Fuji, East Antarctica, combined with an existing long record from the Dome C ice core. Antarctic warming events over the past 720,000 years are most frequent when the Antarctic temperature is slightly below average on orbital time scales, equivalent to an intermediate climate during glacial periods, whereas interglacial and fully glaciated climates are unfavourable for a millennial-scale bipolar seesaw. Numerical experiments using a fully coupled atmosphere-ocean general circulation model with freshwater hosing in the northern North Atlantic showed that climate becomes most unstable in intermediate glacial conditions associated with large changes in sea ice and the Atlantic Meridional Overturning Circulation. Model sensitivity experiments suggest that the prerequisite for the most frequent climate instability with bipolar seesaw pattern during the late Pleistocene era is associated with reduced atmospheric CO2 concentration via global cooling and sea ice formation in the North Atlantic, in addition to extended Northern Hemisphere ice sheets.

Insolation-driven 100,000-year glacial cycles and hysteresis of ice-sheet volume.

The growth and reduction of Northern Hemisphere ice sheets over the past million years is dominated by an approximately 100,000-year periodicity and a sawtooth pattern (gradual growth and fast termination). Milankovitch theory proposes that summer insolation at high northern latitudes drives the glacial cycles, and statistical tests have demonstrated that the glacial cycles are indeed linked to eccentricity, obliquity and precession cycles. Yet insolation alone cannot explain the strong 100,000-year cycle, suggesting that internal climatic feedbacks may also be at work. Earlier conceptual models, for example, showed that glacial terminations are associated with the build-up of Northern Hemisphere 'excess ice', but the physical mechanisms underpinning the 100,000-year cycle remain unclear. Here we show, using comprehensive climate and ice-sheet models, that insolation and internal feedbacks between the climate, the ice sheets and the lithosphere-asthenosphere system explain the 100,000-year periodicity. The responses of equilibrium states of ice sheets to summer insolation show hysteresis, with the shape and position of the hysteresis loop playing a key part in determining the periodicities of glacial cycles. The hysteresis loop of the North American ice sheet is such that after inception of the ice sheet, its mass balance remains mostly positive through several precession cycles, whose amplitudes decrease towards an eccentricity minimum. The larger the ice sheet grows and extends towards lower latitudes, the smaller is the insolation required to make the mass balance negative. Therefore, once a large ice sheet is established, a moderate increase in insolation is sufficient to trigger a negative mass balance, leading to an almost complete retreat of the ice sheet within several thousand years. This fast retreat is governed mainly by rapid ablation due to the lowered surface elevation resulting from delayed isostatic rebound, which is the lithosphere-asthenosphere response. Carbon dioxide is involved, but is not determinative, in the evolution of the 100,000-year glacial cycles.

Ice-sheet collapse and sea-level rise at the Bølling warming 14,600 years ago.

Past sea-level records provide invaluable information about the response of ice sheets to climate forcing. Some such records suggest that the last deglaciation was punctuated by a dramatic period of sea-level rise, of about 20 metres, in less than 500 years. Controversy about the amplitude and timing of this meltwater pulse (MWP-1A) has, however, led to uncertainty about the source of the melt water and its temporal and causal relationships with the abrupt climate changes of the deglaciation. Here we show that MWP-1A started no earlier than 14,650 years ago and ended before 14,310 years ago, making it coeval with the Bølling warming. Our results, based on corals drilled offshore from Tahiti during Integrated Ocean Drilling Project Expedition 310, reveal that the increase in sea level at Tahiti was between 12 and 22 metres, with a most probable value between 14 and 18 metres, establishing a significant meltwater contribution from the Southern Hemisphere. This implies that the rate of eustatic sea-level rise exceeded 40 millimetres per year during MWP-1A.