Glacial Isostatic Adjustment Modelling of the Coast Mountains of British Columbia and Southeastern Alaska

The Coast Mountains in British Columbia and southeastern Alaska contain around 9040 km² of glaciers and ice fields at present. While these glaciers have followed an overall trend of mass loss since the Little Ice Age (or LIA around 300 years before present), the past decade has seen a significant increase in melting rate that is likely to continue due to the effects of climate change. The region is home to a complex tectonic setting, having proximity to the Queen Charlotte-Fairweather transform plate boundary in the northern region and the Cascadia subduction zone (CSZ) in the southern region, which has an associated active volcanic arc underlying the glaciated area. Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) glacier melt data collected between 2000 and 2019 represent a melt rate that is averaged between periods of relatively low mass loss (2000-2009) and high mass loss (2010-2019). As a preliminary test, this average melt rate was assumed to be constant back to the LIA. A history of gridded ice thicknesses was calculated to create an ice loading model for input to a series of forward modelling calculations to determine the crustal response. Predictions of vertical crustal motion are compared to available Global Navigation Satellite System (GNSS) measurements of uplift rate to constrain Earth rheology. The results using this simplified loading model favour a thin lithosphere (around 20-40 km thick) and asthenospheric viscosities on the order of 10¹⁹ Pa s. These values are significantly lower than those of rheological profiles used in extant global GIA models, but are in general agreement with previous GIA modelling of the forearc region of the CSZ. To improve the glacial history model, the Open Global Glacier Model (OGGM), driven by historic climate data and statistically downscaled climate projections, is being employed to create a more accurate loading model and refine our estimates of Earth rheology and regional crustal motion. The best-fitting models will be employed to separate GIA and tectonic components of crustal motion and to generate improved regional sea-level projections.