Thwaites Glacier Research Articles

Calving Dynamics and the Potential Impact of Mélange Buttressing at the Western Calving Front of Thwaites Glacier, West Antarctica

AbstractThe western region of the wide Thwaites Glacier terminus is characterized by a near‐vertical calving front. The grounding line at this western calving front (WCF) rests on a relatively high ridge, behind which exists a reverse‐sloping bed; retreat of the grounding line into this over‐deepening basin could therefore expose deep calving faces that may be subject to ice‐cliff failure. Here, we use the 3D Helsinki Discrete Element Model to identify the factors that control the calving dynamics in this location. We then focus on the ability of mélange to influence these dynamics given the wide embayment in which Thwaites Glacier terminates. We find that calving along the WCF is currently influenced by ice flow across the grounding line and consequent longitudinal tensile stress and rift formation. Calving is slowed in simulations that are initiated with a highly constricted mélange, with a thicker mélange suppressing calving entirely. We liken the constrained simulations to a scenario in which mélange piles behind a large grounded iceberg. In a future which may see calving become a more dominant control on the retreat of Thwaites Glacier, this type of blockage will be necessary for robust force chains to develop and transmit resistive forces to the terminus. The ability of the mélange to hinder calving at this location will be determined by the presence and rigidity of binding land‐fast sea ice and iceberg keel depths. Therefore, it is necessary to represent calving, mélange and sea ice in a single framework to predict the fate of Thwaites Glacier.

Acceleration of Antarctica glaciers at high subglacial heat flow

High subglacial heat flow and volcanic activity in West Antarctica contribute to instability and accelerated flow into the ocean of the West Antarctic ice sheet. In this case, a catastrophic rise in sea level by tens of centimeters – the first meters can occur in a very short geological time (years-decades) due to the rapid sliding of large masses of ice in West Antarctica into the ocean. If the Pine Island (50 cm sea level rise) or Thwaites (65 cm sea level rise) glaciers slide into the ocean, the West Antarctic Ice Sheet will lose support from these glaciers and may begin to collapse. In this case, the sea level will rise by a few meters. Based on Glen’s rheological law for a two-dimensional model of the movement of ice as a nonlinear viscous fluid, the flow velocities of a 3000 m thick glacier were calculated under conditions of adhesion to the bed (~20 m/year) and under conditions of sliding along the bedrock when the lower edge of the glacier melts due to increased heat flow from below (~3000 m/year). These velocities are in good agreement with the velocities of the Pine Island, Thwaites, Amery, Denman and Totten glaciers. The rapid movement of some outlet glaciers in East Antarctica is also likely caused by melting of their bases, suggesting increased subglacial heat flow in these areas of East Antarctica.

Coupled ice–ocean interactions during future retreat of West Antarctic ice streams in the Amundsen Sea sector

Abstract. The Amundsen Sea sector has some of the fastest-thinning ice shelves in Antarctica, caused by high, ocean-driven basal melt rates, which can lead to increased ice streamflow, causing increased sea level rise (SLR) contributions. In this study, we present the results of a new synchronously coupled ice-sheet–ocean model of the Amundsen Sea sector. We use the Wavelet-based, Adaptive-grid, Vertically Integrated ice sheet model (WAVI) to solve for ice velocities and the Massachusetts Institute of Technology general circulation model (MITgcm) to solve for ice thickness and three-dimensional ocean properties, allowing for full mass conservation in the coupled ice–ocean system. The coupled model is initialised in the present day and run forward under idealised warm and cold ocean conditions with a fixed ice front. We find that Thwaites Glacier dominates the future SLR from the Amundsen Sea sector, with a SLR that evolves approximately quadratically over time. The future evolution of Thwaites Glacier depends on the lifespan of small pinning points that form during the retreat. The rate of melting around these pinning points provides the link between future ocean conditions and the SLR from this sector and will be difficult to capture without a coupled ice–ocean model. Grounding-line retreat leads to a progressively larger Thwaites Ice Shelf cavity, leading to a positive trend in total melting, resulting from the increased ice basal surface area. Despite these important sensitivities, Thwaites Glacier retreats even in a scenario with zero ocean-driven melting. This demonstrates that a tipping point may have been passed in these simulations and some SLR from this sector is now committed.

Responses of the Pine Island and Thwaites glaciers to melt and sliding parameterizations

Abstract. The Pine Island and Thwaites glaciers are the two largest contributors to sea level rise from Antarctica. Here we examine the influence of basal friction and ice shelf basal melt in determining projected losses. We examine both Weertman and Coulomb friction laws with explicit weakening as the ice thins to flotation, which many friction laws include implicitly via the effective pressure. We find relatively small differences with the choice of friction law (Weertman or Coulomb) but find losses to be highly sensitive to the rate at which the basal traction is reduced as the area upstream of the grounding line thins. Consistent with earlier work on Pine Island Glacier, we find sea level contributions from both glaciers to vary linearly with the melt volume averaged over time and space, with little influence from the spatial or temporal distribution of melt. Based on recent estimates of melt from other studies, our simulations suggest that the combined melt-driven and sea level rise contribution from both glaciers may not exceed 10 cm by 2200, although the uncertainty in model parameters allows for larger increases. We do not include other factors, such as ice shelf breakup, that might increase loss, or factors such as increased accumulation and isostatic uplift that may mitigate loss.

Widespread seawater intrusions beneath the grounded ice of Thwaites Glacier, West Antarctica

Warm water from the Southern Ocean has a dominant impact on the evolution of Antarctic glaciers and in turn on their contribution to sea level rise. Using a continuous time series of daily-repeat satellite synthetic-aperture radar interferometry data from the ICEYE constellation collected in March-June 2023, we document an ice grounding zone, or region of tidally controlled migration of the transition boundary between grounded ice and ice afloat in the ocean, at the main trunk of Thwaites Glacier, West Antarctica, a strong contributor to sea level rise with an ice volume equivalent to a 0.6-m global sea level rise. The ice grounding zone is 6 km wide in the central part of Thwaites with shallow bed slopes, and 2 km wide along its flanks with steep basal slopes. We additionally detect irregular seawater intrusions, 5 to 10 cm in thickness, extending another 6 km upstream, at high tide, in a bed depression located beyond a bedrock ridge that impedes the glacier retreat. Seawater intrusions align well with regions predicted by the GlaDS subglacial water model to host a high-pressure distributed subglacial hydrology system in between lower-pressure subglacial channels. Pressurized seawater intrusions will induce vigorous melt of grounded ice over kilometers, making the glacier more vulnerable to ocean warming, and increasing the projections of ice mass loss. Kilometer-wide, widespread seawater intrusion beneath grounded ice may be the missing link between the rapid, past, and present changes in ice sheet mass and the slower changes replicated by ice sheet models.

Characterizing sub-glacial hydrology using radar simulations

Abstract. The structure and distribution of sub-glacial water directly influences Antarctic ice mass loss by reducing or enhancing basal shear stress and accelerating grounding line retreat. A common technique for detecting sub-glacial water involves analyzing the spatial variation in reflectivity from an airborne radar echo sounding (RES) survey. Basic RES analysis exploits the high dielectric contrast between water and most other substrate materials, where a reflectivity increase ≥ 15 dB is frequently correlated with the presence of sub-glacial water. There are surprisingly few additional tools to further characterize the size, shape, or extent of hydrological systems beneath large ice masses. We adapted an existing radar backscattering simulator to model RES reflections from sub-glacial water structures using the University of Texas Institute for Geophysics (UTIG) Multifrequency Airborne Radar Sounder with Full-phase Assessment (MARFA) instrument. Our series of hypothetical simulation cases modeled water structures from 5 to 50 m wide, surrounded by bed materials of varying roughness. We compared the relative reflectivity from rounded Röthlisberger channels and specular flat canals, showing both types of channels exhibit a positive correlation between size and reflectivity. Large (> 20 m), flat canals can increase reflectivity by more than 20 dB, while equivalent Röthlisberger channels show only modest reflectivity gains of 8–13 dB. Changes in substrate roughness may also alter observed reflectivity by 3–6 dB. All of these results indicate that a sophisticated approach to RES interpretation can be useful in constraining the size and shape of sub-glacial water features. However, a highly nuanced treatment of the geometric context is necessary. Finally, we compared simulated outputs to actual reflectivity from a single RES flight line collected over Thwaites Glacier in 2022. The flight line crosses a previously proposed Röthlisberger channel route, with an obvious bright bed reflection in the radargram. Through multiple simulations comparing various water system geometries, such as canals and sub-glacial lakes, we demonstrated the important role that topography and water geometry can play in observed RES reflectivity. From the scenarios that we tested, we concluded the bright reflector from our RES flight line cannot be a Röthlisberger channel but could be consistent with a series of flat canals or a sub-glacial lake. However, we note our simulations were not exhaustive of all possible sub-glacial water configurations. The approach outlined here has broad applicability for studying the basal environment of large glaciers. We expect to apply this technique when constraining the geometry and extent of many sub-glacial hydrologic structures in the future. Further research may also include comprehensive investigations of the impact of sub-glacial roughness, substrate heterogeneity, and computational efficiencies enabling more complex and complete simulations.

Frozen Bodies

This essay presents a visual journal that I created for a field season on Thwaites Glacier, Antarctica, as a way both to record states of dissociation/depersonalization and track my gender identity. Drawing on autoethnography as a method of inquiry, I created a pocket-sized journal of body outlines and used oil pastels to document my embodied and emotional responses to conducting fieldwork. I reflect on how my visual journal helped me understand how my body responded to living and working in extreme conditions, how my work and body relate to one another, and how the creation and completion of this journal gave me agency over my experiences of dissociation and queerness. Through this account, I aim to contribute to the growing literature on the personal and political dimensions of polar research—critical glaciology—as well as to demonstrate the potential of autoethnographic methods to facilitate embodied ways of knowing in deep field contexts.

Response of Onshore Oceanic Heat Supply to Yearly Changes in the Amundsen Sea Icescape (Antarctica)

AbstractThe heat transfer between the warm oceanic water and the floating portion of the Antarctic ice sheet (the ice shelves) occurs in a dynamic environment with year‐to‐year changes in the distribution of icebergs and fast‐ice (the “icescape”). Dramatic events such as the collapse of glacier tongues are apparent in satellite images but oceanographic observations are insufficient to capture the synoptic impact of such events on the supply of oceanic heat to ice shelves. This study uses a 3D numerical model and semi‐idealized experiments to examine whether the current high melting rates of ice shelves in the Amundsen Sea could be mitigated by certain icescape configurations. Specifically, the experiments quantify the impacts on oceanic heat supply of presence/absence of the Thwaites Glacier Tongue, Bear Ridge Iceberg Chain, tabular iceberg B22, and fast‐ice cover seaward of Pine Island Ice Shelf (PIS). The experiments reveal that future changes in the coastal icescape are unlikely to reverse the high ice shelf melting rates of the Amundsen Sea, and that icescape changes between 2011 and 2022 actually enhanced them slightly. Ice shelves such as Crosson and Thwaites are found to have multiple viable sources of oceanic heat whose relative importance may shift following icescape reconfigurations but the overall heat supply remains high. Similarly, the formation of a fast‐ice cover seaward of PIS slows down the cavity circulation (by 7%) but does not reduce its heat supply.

Synchronous retreat of Thwaites and Pine Island glaciers in response to external forcings in the presatellite era

Today, relatively warm Circumpolar Deep Water is melting Thwaites Glacier at the base of its ice shelf and at the grounding zone, contributing to significant ice retreat. Accelerating ice loss has been observed since the 1970s; however, it is unclear when this phase of significant melting initiated. We analyzed the marine sedimentary record to reconstruct Thwaites Glacier's history from the early Holocene to present. Marine geophysical surveys were carried out along the floating ice-shelf margin to identify core locations from various geomorphic settings. We use sedimentological data and physical properties to define sedimentary facies at seven core sites. Glaciomarine sediment deposits reveal that the grounded ice in the Amundsen Sea Embayment had already retreated to within ~45 km of the modern grounding zone prior to ca. 9,400 y ago. Sediments deposited within the past 100+ y record abrupt changes in environmental conditions. On seafloor highs, these shifts document ice-shelf thinning initiating at least as early as the 1940s. Sediments recovered from deep basins reflect a transition from ice proximal to slightly more distal conditions, suggesting ongoing grounding-zone retreat since the 1950s. The timing of ice-shelf unpinning from the seafloor for Thwaites Glacier coincides with similar records from neighboring Pine Island Glacier. Our work provides robust new evidence that glacier retreat in the Amundsen Sea was initiated in the mid-twentieth century, likely associated with climate variability.

Geological insights from the newly discovered granite of Sif Island between Thwaites and Pine Island glaciers

AbstractLarge-scale geological structures have controlled the long-term development of the bed and thus the flow of the West Antarctic Ice Sheet (WAIS). However, complete ice cover has obscured the age and exact positions of faults and geological boundaries beneath Thwaites Glacier and Pine Island Glacier, two major WAIS outlets in the Amundsen Sea sector. Here, we characterize the only rock outcrop between these two glaciers, which was exposed by the retreat of slow-flowing coastal ice in the early 2010s to form the new Sif Island. The island comprises granite, zircon U-Pb dated to ~177–174 Ma and characterized by initial ɛNd, 87Sr/86Sr and ɛHf isotope compositions of -2.3, 0.7061 and -1.3, respectively. These characteristics resemble Thurston Island/Antarctic Peninsula crustal block rocks, strongly suggesting that the Sif Island granite belongs to this province and placing the crustal block's boundary with the Marie Byrd Land province under Thwaites Glacier or its eastern shear margin. Low-temperature thermochronological data reveal that the granite underwent rapid cooling following emplacement, rapidly cooled again at ~100–90 Ma and then remained close to the Earth's surface until present. These data help date vertical displacement across the major tectonic structure beneath Pine Island Glacier to the Late Cretaceous.

Disentangling the drivers of future Antarctic ice loss with a historically calibrated ice-sheet model

Abstract. We use an observationally calibrated ice-sheet model to investigate the future trajectory of the Antarctic ice sheet related to uncertainties in the future balance between sub-shelf melting and ice discharge, on the one hand, and the surface mass balance, on the other. Our ensemble of simulations, forced by a panel of climate models from the sixth phase of the Coupled Model Intercomparison Project (CMIP6), suggests that the ocean will be the primary driver of short-term Antarctic mass loss, initiating ice loss in West Antarctica already during this century. The atmosphere initially plays a mitigating role through increased snowfall, leading to an Antarctic contribution to global mean sea-level rise by 2100 of 6 (−8 to 15) cm under a low-emission scenario and 5.5 (−10 to 16) cm under a very high-emission scenario. However, under the very high-emission pathway, the influence of the atmosphere shifts beyond the end of the century, becoming an amplifying driver of mass loss as the ice sheet's surface mass balance decreases. We show that this transition occurs when Antarctic near-surface warming exceeds a critical threshold of +7.5 ∘C, at which the increase in surface runoff outweighs the increase in snow accumulation, a signal that is amplified by the melt–elevation feedback. Therefore, under the very high-emission scenario, oceanic and atmospheric drivers are projected to result in a complete collapse of the West Antarctic ice sheet along with significant grounding-line retreat in the marine basins of the East Antarctic ice sheet, leading to a median global mean sea-level rise of 2.75 (6.95) m by 2300 (3000). Under a more sustainable socio-economic pathway, we find that the Antarctic ice sheet may still contribute to a median global mean sea-level rise of 0.62 (1.85) m by 2300 (3000). However, the rate of sea-level rise is significantly reduced as mass loss is likely to remain confined to the Amundsen Sea Embayment, where present-day climate conditions seem sufficient to commit to a continuous retreat of Thwaites Glacier.

What can radar-based measures of subglacial hydrology tell us about basal shear stress? A case study at Thwaites Glacier, West Antarctica

Abstract Ice sheet models use observations to infer basal shear stress, but the variety of methods and datasets available has resulted in a wide range of estimates. Radar-based metrics such as reflectivity and specularity content have been used to characterize subglacial hydrologic conditions that are linked to spatial variations in basal shear stress. We explore whether radar metrics can be used to inform models about basal shear stress. At Thwaites Glacier, West Antarctica, we sample basal shear stress inversions across a wide range of ice sheet models to see how the basal shear stress distribution changes in regions of varying relative reflectivity and specularity content. Our results reveal three key findings: (1) Regions of high specularity content exhibit lower mean basal shear stresses (2) Wet and bumpy regions, as characterized by high relative reflectivity and low specularity content, exhibit higher mean basal shear stresses (3) Models disagree about what basal shear stress should be at the onset of rapid ice flow and high basal melt where relative reflectivity and specularity content are low.

What is the primary culprit that causes the "Doomsday Glacier" in Antarctica?

Antarctica, situated in the southernmost part of the earth is surrounded by the Southern Ocean. It has an extremely important role in tuning global climate, both because of its geographic location, and its giant ice mass, that comprises ~90% of the world's ice. In the past decades, the Thwaites Glacier in the Amundsen Sea of West Antarctica is undergoing the fastest recession in the region. The ice loss in Thwaites Glacier is currently responsible for roughly four percent of the global sea-level rise, which has been attributed to climate change and ocean warming. Due to the continuous collapsing and melting of Thwaites Glacier and the severe threat to humans, scientists gave it a terrifying name "Doomsday Glacier". With increasingly geological and geophysical studies conducted in West Antarctica, geothermal heat flux has been discovered to play a vital role in icesheet retreating. The rapidly retreating Thwaites and Pope glaciers are underlain by areas of largely elevated geothermal heat flow, which relates to the tectonic and magmatic history of the West Antarctic Rift System in this region, suggesting that this area is coupled to the dynamics of the underlying lithosphere. The collapsing and melting of the West Antarctic glaciers are without doubt a realistic and complex issue. From the current geological surveys, the heat flux of the crust of West Antarctica appears to be accelerating, coupled with frequent earthquakes and volcanoes. Geothermal features such as hot water lakes, thermal rivers, and giant ice caves beneath the glaciers have been continuously discovered. Therefore, it is reasonable to speculate that geothermal effects play a significant role in modifying the vast ice masses, causing glacier sliding, cracking, collapsing, and ultimately melting, and create conditions for a warming climate to melt West Antarctic glaciers. Many studies also suggest that warm ocean water is intruding beneath the glaciers across the grounding line, leading to melting of glacier bottom, which has been generally considered to be associated with human-emitted greenhouse gases, but it is thought here not to be a primary factor for glacier melting, but most likely a secondary factor. It is believed that primary and secondary factors, Milankovitch orbits, black body radiation, solar activities, human activities, and other factors are all interconnected to form a feedback loop between the glacier base and the ocean, or even a positive feedback loop, which further accelerates the collapsing and melting of the west Antarctic glaciers

Insights into the vulnerability of Antarctic glaciers from the ISMIP6 ice sheet model ensemble and associated uncertainty

Abstract. The Antarctic Ice Sheet represents the largest source of uncertainty in future sea level rise projections, with a contribution to sea level by 2100 ranging from −5 to 43 cm of sea level equivalent under high carbon emission scenarios estimated by the recent Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). ISMIP6 highlighted the different behaviors of the East and West Antarctic ice sheets, as well as the possible role of increased surface mass balance in offsetting the dynamic ice loss in response to changing oceanic conditions in ice shelf cavities. However, the detailed contribution of individual glaciers, as well as the partitioning of uncertainty associated with this ensemble, have not yet been investigated. Here, we analyze the ISMIP6 results for high carbon emission scenarios, focusing on key glaciers around the Antarctic Ice Sheet, and we quantify their projected dynamic mass loss, defined here as mass loss through increased ice discharge into the ocean in response to changing oceanic conditions. We highlight glaciers contributing the most to sea level rise, as well as their vulnerability to changes in oceanic conditions. We then investigate the different sources of uncertainty and their relative role in projections, for the entire continent and for key individual glaciers. We show that, in addition to Thwaites and Pine Island glaciers in West Antarctica, Totten and Moscow University glaciers in East Antarctica present comparable future dynamic mass loss and high sensitivity to ice shelf basal melt. The overall uncertainty in additional dynamic mass loss in response to changing oceanic conditions, compared to a scenario with constant oceanic conditions, is dominated by the choice of ice sheet model, accounting for 52 % of the total uncertainty of the Antarctic dynamic mass loss in 2100. Its relative role for the most dynamic glaciers varies between 14 % for MacAyeal and Whillans ice streams and 56 % for Pine Island Glacier at the end of the century. The uncertainty associated with the choice of climate model increases over time and reaches 13 % of the uncertainty by 2100 for the Antarctic Ice Sheet but varies between 4 % for Thwaites Glacier and 53 % for Whillans Ice Stream. The uncertainty associated with the ice–climate interaction, which captures different treatments of oceanic forcings such as the choice of melt parameterization, its calibration, and simulated ice shelf geometries, accounts for 22 % of the uncertainty at the ice sheet scale but reaches 36 % and 39 % for Institute Ice Stream and Thwaites Glacier, respectively, by 2100. Overall, this study helps inform future research by highlighting the sectors of the ice sheet most vulnerable to oceanic warming over the 21st century and by quantifying the main sources of uncertainty.

Characterizing bed roughness on the Antarctic continental margin

Abstract Spatial variability in bed topography, characterized as bed roughness, impacts ice-sheet flow and organization and can be used to infer subglacial conditions and processes, yet is difficult to quantify due to sparse observations. Paleo-subglacial beds of formerly expanded glaciers found across the Antarctic continental shelf are well preserved, have relatively limited post-glacial sediment cover and contain glacial landforms that can be resolved at sub-meter vertical scales. We analyze high-resolution bathymetry offshore of Pine Island and Thwaites glaciers in the Amundsen Sea to explore spatial variability of bed roughness where streamlined subglacial landforms allow for the determination of ice-flow direction. We quantify bed roughness using std dev. and Fast Fourier Transform methods, each employed at local (100 km) and regional (101–2 km) scales and in along- and across-flow orientations to determine roughness expressions across spatial scales. We find that the magnitude of roughness is impacted by the parameters selected – which are often not sufficiently reported in studies – to quantify roughness. Important spatial patterns can be discerned from high-resolution bathymetry, highlighting both its usefulness in identifying patterns of streaming ice flow and underscores the need for a standardized way of characterizing topographic variability.

The stability of present-day Antarctic grounding lines – Part 2: Onset of irreversible retreat of Amundsen Sea glaciers under current climate on centennial timescales cannot be excluded

Abstract. Observations of ocean-driven grounding-line retreat in the Amundsen Sea Embayment in Antarctica raise the question of an imminent collapse of the West Antarctic Ice Sheet. Here we analyse the committed evolution of Antarctic grounding lines under the present-day climate. To this aim, we first calibrate a sub-shelf melt parameterization, which is derived from an ocean box model, with observed and modelled melt sensitivities to ocean temperature changes, making it suitable for present-day simulations and future sea level projections. Using the new calibration, we run an ensemble of historical simulations from 1850 to 2015 with a state-of-the-art ice sheet model to create model instances of possible present-day ice sheet configurations. Then, we extend the simulations for another 10 000 years to investigate their evolution under constant present-day climate forcing and bathymetry. We test for reversibility of grounding-line movement in the case that large-scale retreat occurs. In the Amundsen Sea Embayment we find irreversible retreat of the Thwaites Glacier for all our parameter combinations and irreversible retreat of the Pine Island Glacier for some admissible parameter combinations. Importantly, an irreversible collapse in the Amundsen Sea Embayment sector is initiated at the earliest between 300 and 500 years in our simulations and is not inevitable yet – as also shown in our companion paper (Part 1, Hill et al., 2023). In other words, the region has not tipped yet. With the assumption of constant present-day climate, the collapse evolves on millennial timescales, with a maximum rate of 0.9 mm a−1 sea-level-equivalent ice volume loss. The contribution to sea level by 2300 is limited to 8 cm with a maximum rate of 0.4 mm a−1 sea-level-equivalent ice volume loss. Furthermore, when allowing ice shelves to regrow to their present geometry, we find that large-scale grounding-line retreat into marine basins upstream of the Filchner–Ronne Ice Shelf and the western Siple Coast is reversible. Other grounding lines remain close to their current positions in all configurations under present-day climate.

Simulation of flexural-gravity wave propagation for elastic plates in shallow water using an energy-stable finite difference method with weakly enforced boundary and interface conditions

We introduce an energy stable, high-order-accurate finite difference approximation of the dynamic, pure bending Kirchhoff plate equations for complex geometries and spatially variable properties. We utilize the summation-by-parts (SBP) framework to discretize the biharmonic operator with variable coefficients, with attention given to free and clamped boundary conditions and corner conditions. Energy conservation is established by combining SBP boundary closures with weak enforcement of the boundary and interface conditions using a penalty (simultaneous approximation term, SAT) technique. Then we couple the plate equations to the shallow water equations to study flexural-gravity wave propagation, and prove that the semi-discrete system of equations is self-adjoint. We demonstrate the stability and accuracy properties of the method on curvilinear multiblock grids using the method of manufactured solutions. The method, which we provide in an open-source code, is then used to model ocean wave interactions with the Thwaites Glacier and Pine Island Ice Shelf in the Amundsen Sea off the coast of West Antarctica.

Unsupervised machine learning detection of iceberg populations within sea ice from dual-polarisation SAR imagery

Accurate quantification of iceberg populations is essential to inform estimates of Southern Ocean freshwater and heat balances as well as shipping hazards. The automated operational monitoring of icebergs remains challenging, largely due to a lack of generality in existing approaches. Previous efforts to map icebergs have often exploited synthetic aperture radar (SAR) data but the majority are designed for open water situations, require significant operator input, and are susceptible to the substantial spatial and temporal variability in backscatter that characterises SAR time-series. We propose an adaptive unsupervised classification procedure based on Sentinel 1 SAR data and a recursive Dirichlet Process implementation of Bayesian Gaussian Mixture Model. The approach is robust to inter-scene variability and can identify icebergs even within complex environments containing mixtures of open water, sea ice and icebergs of various sizes. For the study area in the Amundsen Sea Embayment, close to the calving front of Thwaites Glacier, our classifier achieved a mean pixel-wise F1 score against manual iceberg delineations from the SAR scenes of 0.960±0.018 with a corresponding object-level F1 score of 0.729±0.086. The method provides an excellent basis for estimation of total near-shore iceberg populations and has inherent potential for scalability that other approaches lack.

Towards modelling of corrugation ridges at ice-sheet grounding lines

Abstract. Improvements in the resolution of sea-floor mapping techniques have revealed extremely regular, sub-metre-scale ridge landforms produced by the tidal flexure of ice-shelf grounding lines as they retreated very rapidly (i.e. at rates of several kilometres per year). Guided by such novel sea-floor observations from Thwaites Glacier, West Antarctica, we present three mathematical models for the formation of these corrugation ridges at a tidally migrating grounding line (that is retreating at a constant rate), where each ridge is formed by either constant till flux to the grounding line, till extrusion from the grounding line, or the resuspension and transport of grains from the grounding-zone bed. We find that both till extrusion (squeezing out till like toothpaste as the ice sheet re-settles on the sea floor) and resuspension and transport of material can qualitatively reproduce regular, delicate ridges at a retreating grounding line, as described by sea-floor observations. By considering the known properties of subglacial sediments, we agree with existing schematic models that the most likely mechanism for ridge formation is till extrusion at each low-tide position, essentially preserving an imprint of the ice-sheet grounding line as it retreated. However, when realistic (shallow) bed slopes are used in the simulations, ridges start to overprint one another, suggesting that, to preserve the regular ridges that have been observed, grounding line retreat rates (driven by dynamic thinning?) may be even higher than previously thought.

Large Earthquakes in Subduction Zones around the Polar Regions as a Possible Reason for Rapid Climate Warming in the Arctic and Glacier Collapse in West Antarctica

A correlation is observed between changes in the level of Earth’s seismic activity and increments of the atmospheric methane concentration over the past 40 years. Trigger mechanisms are proposed for methane emissions and glacier collapse in polar regions. These mechanisms are due to deformation waves caused by large earthquakes in subduction zones located near the polar regions: the Aleutian and Kuril–Kamchatka subduction zones, closest to the Arctic, and the Antarctica–Chilean and Tonga–Kermadec–Macquarie subduction zones. Disturbances of the lithosphere are transmitted over the distances of 3000–4000 km and more at a speed of about 100 km/year. Additional associated stresses come to the Arctic and Antarctica several decades after the occurrence of large earthquakes. In the Arctic zone, additional stresses affect the low-permeability structure of gas bearing sedimentary strata, causing increased methane emission and climate warming. In West Antarctica, deformation waves could trigger the acceleration and intensive collapse of West Antarctic glaciers, which has been observed since the 1970s. These waves are also capable of activating dormant volcanoes located under the sheet glaciers of West Antarctica, leading to an increase in heat flux, to the melting of ice at the glaciers’ base, and to their accelerated sliding towards the ocean, as is happening with the Thwaites Glacier.