Search for: All records

Seismic travel time inversion code from: Geophysical evidence of tectonic structures beneath Thwaites Glacier, West Antarctica: influence on glacier dynamics. Submitted to Journal of Geophysical Research: Solid Earth 

Flexure modeling code for manuscript: Geophysical evidence of tectonic structures beneath Thwaites Glacier, West Antarctica: influence on glacier dynamics. Submitted to Journal of Geophysical Research: Solid Earth 

Results from seismic travel time inversion and potential fields modeling in: Geophysical evidence of tectonic structures beneath Thwaites Glacier, West Antarctica: influence on glacier dynamics 

Abstract. The contribution of the Greenland Ice Sheet (GIS) to sea level rise (SLR) is accelerating and there is an urgent need to improve predictions of when and from what parts of the ice sheet Greenland will contribute its first meter. Estimating the volume of Greenland ice that was lost during past warm periods offers a way to constrain the ice sheet’s response to future warming. Sub-ice sediment and bedrock, retrieved from deep ice core campaigns or targeted drilling efforts, yield critical and direct information about past ice-free conditions. However, it is challenging to scale the few available sub-ice point measurements to the geometry of the entire ice sheet. Here, we provide a framework for assessing sea-level potential, which we define as the amount the GIS has contributed to sea level when a particular location in Greenland is ice-free, from an ensemble of ice-sheet model simulations representing a wide range of plausible deglaciation scenarios. An assessment of dominant sources of uncertainty in our paleo ice sheet modelling, including climate forcing, ice-sheet initialization, and solid-Earth properties, reveals spatial patterns in the sensitivity of the ice sheet to these processes and related feedbacks. We find that the sea-level potential of central Greenland is most sensitive to lithospheric feedbacks and ice-sheet initialization, whereas the ice-sheet margins are most sensitive to climate forcing parameters. Our framework allows us to quantify the local and regional uncertainty in sea-level potential, which we use to evaluate the GIS bedrock according to the usefulness of information sub-ice sediments and bedrock provide about past ice-sheet geometry. Through our ensemble approach, we can assign a plausible range of GIS contributions to global sea level for deglaciated conditions at any site. Our results identify primarily areas in southwest Greenland, and secondarily north Greenland, as best-suited for subglacial access drilling that seeks to constrain the response of the ice sheet to past and future warming. 

Abstract Thwaites Glacier (TG) plays an important role in future sea-level rise (SLR) contribution from the West Antarctic Ice Sheet. Recent observations show that TG is losing mass, and its grounding zone is retreating. Previous modeling has produced a wide range of results concerning whether, when, and how rapidly further retreat will occur under continued warming. These differences arise at least in part from ill-constrained processes, including friction from the bed, and future atmosphere and ocean forcing affecting ice-shelf and grounding-zone buttressing. Here, we apply the Ice Sheet and Sea-level System Model (ISSM) with a range of specifications of basal sliding behavior in response to varying ocean forcing. We find that basin-wide bed character strongly affects TG's response to sub-shelf melt by modulating how changes in driving stress are balanced by the bed as the glacier responds to external forcing. Resulting differences in dynamic thinning patterns alter modeled grounding-line retreat across Thwaites' catchment, affecting both modeled rates and magnitudes of SLR contribution from this critical sector of the ice sheet. Bed character introduces large uncertainties in projections of TG under equal external forcing, pointing to this as a crucial constraint needed in predictive models of West Antarctica. 

Abstract Thwaites Glacier in West Antarctica has been identified as a route to destabilization of the whole West Antarctic Ice Sheet, potentially leading to several meters of sea‐level rise. However, future evolution of Thwaites Glacier remains uncertain due to a lack of detailed knowledge about its basal boundary that will affect how its retreat proceeds. Here we aim to improve understanding of the basal boundary in the lower part of Thwaites Glacier by modeling the crustal structures that are related to the bed‐type distribution and therefore influence the basal slip. We combine long‐offset seismic, and gravity‐ and magnetic‐anomaly data to model the crustal structures along two 120 km lines roughly parallel to ice flow. We find a sedimentary basin 40 km in length in the along‐flow direction, with a maximum thickness of 1.7 0.2 km, and two mafic intrusions at 5–10 km depth that vary in maximum thickness between 3.8 and 8.6 km. The sedimentary basin and major mafic intrusions we modeled are likely related to the multi‐stage tectonic evolution of the West Antarctic Rift System. Thwaites Glacier flows across a tectonic boundary within our study site, indicating it flows across tectonically formed structures. The varying geology and resulting variations in bed types demonstrate the influence of tectonics on Thwaites Glacier dynamics. 

Abstract The rapidly changing Thwaites Ice Shelf is crucial for understanding ice‐shelf dynamical processes and their implications for sea‐level rise from Antarctica. Fractures, particularly their vertical structure, are key to ice‐shelf structural integrity but remain poorly measured. To address this, we developed a fracture‐characterization workflow using ICESat‐2 ATL03 geolocated photon heights, producing the first time‐series vertical measurements of fractures across Thwaites from 2018 to 2024. We introduced the fracture depth/freeboard ratio as a normalized metric to quantify vertical fracture extent, serving as an indicator of structural damage. This metric enabled us to track fracture evolution in both the eastern ice shelf and western glacier tongue. In the eastern section, fracturing intensified along the northwestern shear zone and near the grounding line, in a positive feedback loop between enhanced fracturing and accelerated flow. The western section maintained an active rift formation zone about 15 km downstream of the historical grounding line. Flow velocity changes in this section were primarily confined to the unconstrained downstream portion, exhibiting an overall deceleration trend, while the upstream area remained stable. This contrast highlights the role of lateral margin conditions in governing ice‐shelf fracture and flow behavior. Changes in the eastern section showed some correspondence with warm winter air temperatures, reduced sea ice, and persistent warm ocean anomalies at shallower depths, suggesting that atmosphere‐sea ice‐ocean interactions influence ice‐shelf structural integrity through basal processes. Future research should integrate satellite‐derived fracture observations with numerical models of ice fracture and flow to better capture the dynamics of ice‐shelf weakening and retreat. 

Abstract Glacier-bed characteristics that are poorly known and modeled are important in projected sea-level rise from ice-sheet changes under strong warming, especially in the Thwaites Glacier drainage of West Antarctica. Ocean warming may induce ice-shelf thinning or loss, or thinning of ice in estuarine zones, reducing backstress on grounded ice. Models indicate that, in response, more-nearly-plastic beds favor faster ice loss by causing larger flow acceleration, but more-nearly-viscous beds favor localized near-coastal thinning that could speed grounding-zone retreat into interior basins where marine-ice-sheet instability or cliff instability could develop and cause very rapid ice loss. Interpretation of available data indicates that the bed is spatially mosaicked, with both viscous and plastic regions. Flow against bedrock topography removes plastic lubricating tills, exposing bedrock that is eroded on up-glacier sides of obstacles to form moats with exposed bedrock tails extending downglacier adjacent to lee-side soft-till bedforms. Flow against topography also generates high-ice-pressure zones that prevent inflow of lubricating water over distances that scale with the obstacle size. Extending existing observations to sufficiently large regions, and developing models assimilating such data at the appropriate scale, present large, important research challenges that must be met to reliably project future forced sea-level rise. 

Abstract. Direct observations of the size of the Greenland Ice Sheet during Quaternary interglaciations are sparse yet valuable for testing numerical models of ice-sheet history and sea level contribution. Recent measurements of cosmogenicnuclides in bedrock from beneath the Greenland Ice Sheet collected duringpast deep-drilling campaigns reveal that the ice sheet was significantlysmaller, and perhaps largely absent, sometime during the past 1.1 millionyears. These discoveries from decades-old basal samples motivate new,targeted sampling for cosmogenic-nuclide analysis beneath the ice sheet.Current drills available for retrieving bed material from the US IceDrilling Program require < 700 m ice thickness and a frozen bed,while quartz-bearing bedrock lithologies are required for measuring a largesuite of cosmogenic nuclides. We find that these and other requirementsyield only ∼ 3.4 % of the Greenland Ice Sheet bed as asuitable drilling target using presently available technology. Additionalfactors related to scientific questions of interest are the following: which areas of thepresent ice sheet are the most sensitive to warming, where would a retreating icesheet expose bare ground rather than leave a remnant ice cap, andwhich areas are most likely to remain frozen bedded throughout glacialcycles and thus best preserve cosmogenic nuclides? Here we identifylocations beneath the Greenland Ice Sheet that are best suited for potentialfuture drilling and analysis. These include sites bordering Inglefield Landin northwestern Greenland, near Victoria Fjord and Mylius-Erichsen Land innorthern Greenland, and inland from the alpine topography along the icemargin in eastern and northeastern Greenland. Results from cosmogenic-nuclide analysis in new sub-ice bedrock cores from these areas would help to constrain dimensions of the Greenland Ice Sheet in the past.