<p><strong>Abstract.</strong> Subglacial water influences the dynamics of ice masses. The state of subglacial pore water, whether liquid or frozen, is associated with differences in electrical resistivity that span several orders of magnitude, hence liquid water can be inferred from electrical resistivity depth profiles. Such profiles can be obtained from inversions of time domain electromagnetics (TEM) soundings, but these are often non-unique. Here, we adapt an existing Bayesian

more »

... al algorithm (<q>MuLTI</q>) to the inversion of TEM data constrained by independent depth constraints, to provide statistical properties and uncertainty analysis of the resistivity profile with depth. The method was applied to ground-based TEM data acquired on the terminus of the Norwegian glacier Midtdalsbreen, with depth constraints provided by co-located ground penetrating radar data. Our inversion shows that the glacier bed is directly underlain by material of resistivity 10²&amp;thinsp;Ωm ±&amp;thinsp;100&amp;thinsp;%, with thickness 5&amp;ndash;40&amp;thinsp;m, in turn underlain by a highly conductive basement (10⁰&amp;thinsp;Ωm ±&amp;thinsp;15&amp;thinsp;%). High resistivity material, 5&amp;thinsp;&amp;times;&amp;thinsp;10⁴&amp;thinsp;Ωm ±&amp;thinsp;25&amp;thinsp;%, exists at the front of the glacier. All uncertainties are defined by the interquartile range of the posterior resistivity distribution. Combining these resistivity profiles with co-located seismic shear-wave velocity inversions to further reduce ambiguity in the hydro-geological interpretation of the subsurface, we propose a new 3D interpretation of the Midtdalsbreen subglacial material partitioned into partially frozen sediment, frozen sediment/permafrost and weathered/fractured bedrock with saline water.</p>