Premium
Joint Inversion for Surface Accumulation Rate and Geothermal Heat Flow From Ice‐Penetrating Radar Observations at Dome A, East Antarctica. Part I: Model Description, Data Constraints, and Inversion Results
Author(s) -
Wolovick M. J.,
Moore J. C.,
Zhao L.
Publication year - 2021
Publication title -
journal of geophysical research: earth surface
Language(s) - English
Resource type - Journals
eISSN - 2169-9011
pISSN - 2169-9003
DOI - 10.1029/2020jf005937
Subject(s) - geology , ice sheet , antarctic ice sheet , inversion (geology) , radar , ice stream , ice sheet model , geothermal gradient , greenland ice sheet , ground penetrating radar , geophysics , geothermal heating , climatology , geomorphology , structural basin , cryosphere , sea ice , computer science , geothermal energy , telecommunications
Ice‐penetrating radar data contain a wealth of information about the bed and internal structure of the ice sheet. While these data have long been used to diagnose the presence of basal water, infer attenuation rates, or explore the internal stratigraphy of the ice sheet, they have rarely been used jointly in a formal inverse model for the ice sheet temperature structure. Here, we invert a coupled thermomechanical ice sheet and basal hydrology model to infer both geothermal heat flow (GHF) and accumulation rate from multiple classes of radar observations in the area around Dome A, East Antarctica. Our forward model solves for a coupled steady state between the ice sheet flow field, temperature, and basal hydrology, including melt, water transport, and freeze‐on. We fit radar observations of basal water, freeze‐on, and internal layers, along with a GHF prior based on aeromagnetic observations. We minimize the combined misfit function by first using an evolutionary algorithm followed by localized perturbation tests. In addition to inferring the spatial distribution of GHF and accumulation rate, we are also able to estimate the uncertainty about our best‐fit answer, as well as quantify how our result depends on each individual data constraint. Our results demonstrate a new method for combining multiple glaciological constraints into a single inverse model of the ice sheet, and give us a more rigorous picture of the information content provided by each data set. In a companion paper we analyze and interpret the best‐fit model.