Your search

Crossing a key atmospheric CO2 threshold triggered a fundamental global climate reorganisation ~34 million years ago (Ma) establishing permanent Antarctic ice sheets. Curiously, a more dramatic CO2 decline (~800–400 ppm by the Early Oligocene(~27 Ma)), postdates initial ice sheet expansion but the mechanisms driving this later, rapid drop in atmospheric carbon during the early Oligocene remains elusive and controversial. Here we use marine seismic reflection and borehole data to reveal an unprecedented accumulation of early Oligocene strata (up to 2.2 km thick over 1500 × 500 km) with a major biogenic component in the Australian Southern Ocean. High-resolution ocean simulations demonstrate that a tectonically-driven, one-off reorganisation of ocean currents, caused a unique period where current instability coincided with high nutrient input from the Antarctic continent. This unrepeated and short-lived environment favoured extreme bioproductivity and enhanced sediment burial. The size and rapid accumulation of this sediment package potentially holds ~1.067 × 1015 kg of the ‘missing carbon’ sequestered during the decline from an Eocene high CO2-world to a mid-Oligocene medium CO2-world, highlighting the exceptional role of the Southern Ocean in modulating long-term climate.

The East Antarctic Ice Sheet stores a vast amount of freshwater, which makes it the single largest potential contributor to future sea-level rise. However, the lack of well-constrained geological records of past ice sheet changes impedes model validation, hampers mass balance estimates, and inhibits examination of ice loss mechanisms. Here we identify rapid ice-sheet thinning in coastal Dronning Maud Land from Early to Middle Holocene (9000–5000 years ago) using a deglacial chronology based on in situ cosmogenic nuclide surface exposure dates from central Dronning Maud Land, in concert with numerical simulations of regional and continental ice-sheet evolution. Regional sea-level changes reproduced from our refined ice-load history show a highstand at 9000–8000 years ago. We propose that sea-level rise and a concomitant influx of warmer Circumpolar Deep Water triggered ice shelf breakup via the marine ice sheet instability mechanism, which led to rapid thinning of upstream coastal ice sheet sectors.

We investigated mass balance changes over five ice rises in the last few decades near Fimbul and Nivl ice shelves in central Dronning Maud Land. We use the Input-Output Method constrained using field-based geophysical measurements conducted during the austral summers of 2012–14 over three ice rises near the Fimbul Ice Shelf. Further, we use satellite altimetry data (ICESat, ICESat-2, and CryoSat-2) to estimate geodetic mass balance over all five ice rises in recent decades. Both field- and satellite-based estimates show that until 2010, three out of five ice rises were thickening (0.4–0.2 m_ieq a⁻¹) while two were close to balance. However, over the last decade, the ice rises thickening previously started to thin (−0.2–−0.6 m_ieq a⁻¹) while the other two remained close to balance. Much of this variability is likely associated with regional surface mass balance trends, with each ice rise exhibiting its characteristics depending on its local glaciological settings.

The coastal Droning Maud Land in East Antarctica is characterized by small ice shelves with numbers of promontories and locally grounded isles, both called ice rises. These ice rises are typically dome-shaped and surface elevations are hundreds of meters above the surrounding ice shelves, which cause strong orographic effects on surface mass balance (SMB). We conducted shallow ice-penetrating radar sounding to visualize firn stratigraphy in the top 35 m over ~400 km of profiles across the Nivlisen Ice Shelf, and in a grid pattern over two adjacent ice rises (Djupranen and Leningradkollen). We tracked six reflectors (isochrones) and dated them using two ice cores taken at the ice rise summits, from which SMB over six periods in the past three decades was retrieved. The overall SMB pattern across the ice shelf remained similar for all periods; however, the eastwest contrast in SMB varies by a factor of 1.5–2 between the Leningradkollen and Djupranen grounding lines. The SMB patterns over the ice rises are more varied owing to complex interactions between topography, snowfall and wind. We use our results to evaluate the regional climate model RACMO2.3p2 in terms of the spatial SMB distribution and temporal changes over the ice shelf and ice rises at regional scale.

Changes in ocean-driven basal melting have a key influence on the stability of ice shelves, the mass loss from the ice sheet, ocean circulation, and global sea level rise. Coupled ice sheet–ocean models play a critical role in understanding future ice sheet evolution and examining the processes governing ice sheet responses to basal melting. However, as a new approach, coupled ice sheet–ocean systems come with new challenges, and the impacts of solutions implemented to date have not been investigated. An emergent feature in several contributing coupled models to the 1st Marine Ice Sheet–Ocean Model Intercomparison Project (MISOMIP1) was a time-varying oscillation in basal melt rates. Here, we use a recently developed coupling framework, FISOC (v1.1), to connect the modified ocean model ROMSIceShelf (v1.0) and ice sheet model Elmer/Ice (v9.0), to investigate the origin and implications of the feature and, more generally, the impact of coupled modeling strategies on the simulated basal melt in an idealized ice shelf cavity based on the MISOMIP setup. We found the spatial-averaged basal melt rates (3.56 m yr−1) oscillated with an amplitude ∼0.7 m yr−1 and approximate period of ∼6 years between year 30 and 100 depending on the experimental design. The melt oscillations emerged in the coupled system and the standalone ocean model using a prescribed change of cavity geometry. We found that the oscillation feature is closely related to the discretized ungrounding of the ice sheet, exposing new ocean, and is likely strengthened by a combination of positive buoyancy–melt feedback and/or melt–geometry feedback near the grounding line, and the frequent coupling of ice geometry and ocean evolution. Sensitivity tests demonstrate that the oscillation feature is always present, regardless of the choice of coupling interval, vertical resolution in the ocean model, tracer properties of cells ungrounded by the retreating ice sheet, or the dependency of friction velocities to the vertical resolution. However, the amplitude, phase, and sub-cycle variability of the oscillation varied significantly across the different configurations. We were unable to ultimately determine whether the feature arises purely due to numerical issues (related to discretization) or a compounding of multiple physical processes amplifying a numerical artifact. We suggest a pathway and choices of physical parameters to help other efforts understand the coupled ice sheet–ocean system using numerical models.