Description: The shape of ice shelf cavities are a major source of uncertainty in understanding ice‐ocean interactions. This limits assessments of the response of the Antarctic ice sheets to climate change. Here we use vibroseis seismic reflection surveys to map the bathymetry beneath the Ekström Ice Shelf, Dronning Maud Land. The new bathymetry reveals an inland‐sloping trough, reaching depths of 1,100 m below sea level, near the current grounding line, which we attribute to erosion by palaeo‐ice streams. The trough does not cross‐cut the outer parts of the continental shelf. Conductivity‐temperature‐depth profiles within the ice shelf cavity reveal the presence of cold water at shallower depths and tidal mixing at the ice shelf margins. It is unknown if warm water can access the trough. The new bathymetry is thought to be representative of many ice shelves in Dronning Maud Land, which together regulate the ice loss from a substantial area of East Antarctica.

Description: The onset and high upstream ice surface velocities of the North East Greenland Ice Stream (NEGIS) are not yet well reproducible in ice sheet models. A major uncertainty remains the understanding of basal sliding and a parameterization of basal conditions. In this study, we assess the slow-flowing part of the NEGIS in a systematic analysis of the basal conditions and investigate the increased ice flow. We analyze the spectral basal roughness in correlation with basal return power from an airborne radar survey with AWIs ultra-wideband radar system in 2018 and compare our results with current ice flow geometry and ice surface flow. We observe a roughness anisotropy where the ice stream widens, indicating a change from a smooth and soft bed to a harder bedrock as well as the evolution of elongated subglacial landforms. In addition, at the upstream part of the NEGIS we find a clear zoning of the bedrock return power, indicating an increased water content at the base of the ice stream. At the downstream part, we observe an increased bedrock return power throughout the entire width of the ice stream and outside its margins, indicating enhanced melting and the distribution of basal water beyond the shear zones.

Description: Outlet glaciers of the Greenland Ice Sheet transport ice from the interior to the ocean and contribute directly to sea level rise because because discharge and ablation often exceed the accumulation. To develop a better understanding of these fast flowing glaciers, we investigate the basal conditions of Store Glacier, a large outlet glacier flowing into Uummannaq Fjord in West Greenland. We use two crossing seismic profiles acquired near the centreline, 30 km upstream of the calving front, to interpret the physical nature of the ice and bed. We identify one notably englacial and two notably subglacial seismic reflections on both profiles. The englacial reflection represents a change in crystal orientation fabric, interpreted to be the Holocene–Wisconsin transition. From Amplitude Versus Angle (AVA) analysis we infer that the deepest ∼80 m of ice of the parallel-flow profile below this reflection is anisotropic with an enhancement of simple shear of ∼2. The ice is underlain by ∼45 m of unconsolidated sediments, below which there is a strong reflection caused by the transition to consolidated sediments. In the across-flow profile subglacial properties vary over small scale and the polarity of the ice–bed reflection switches from positive to negative. We interpret these as patches of different basal slipperiness associated with variable amounts of water. Our results illustrate variability in basal properties, and hence ice-bed coupling, at a spatial scale of ∼100 m, highlighting the need for direct observations of the bed to improve the basal boundary conditions in ice-dynamic models.

Description: To resolve the mechanisms behind the major climate reorganisation which occurred between 0.9 and 1.2Ma, the recovery of a suitable 1.5 million-year-old ice core is fundamental. The quest for such an Oldest Ice core requires a number of key boundary conditions, of which the poorly known basal geothermal heat flux (GHF) is lacking. We use a transient thermodynamical 1D vertical model that solves for the rate of change of temperature in the vertical, with surface temperature and modelled GHF as boundary conditions. For each point on the ice sheet, the model is forced with variations in atmospheric conditions over the last 2Ma, and modelled ice-thickness variations. The process is repeated for a range of GHF values to determine the value of GHF that marks the limit between frozen and melting conditions over the whole ice sheet, taking into account 2Ma of climate history. These threshold values of GHF are statistically compared to existing GHF data sets. The new probabilistic GHF fields obtained for the ice sheet thus provide the missing boundary conditions in the search for Oldest Ice. High spatial resolution radar data are examined locally in the Dome Fuji and Dome C regions, as these represent the ice core community's primary drilling sites. GHF, bedrock variability, ice thickness and other essential criteria combined highlight a dozen major potential Oldest Ice sites in the vicinity of Dome Fuji and Dome C, where GHF allows for Oldest Ice.

Description: We investigate the propagation of seismic waves in anisotropic ice. Two effects are important: (i) sudden changes in crystal orientation fabric (COF) lead to englacial reflections; (ii) the anisotropic fabric induces an angle dependency on the seismic velocities and, thus, recorded travel times. Velocities calculated from the polycrystal elasticity tensor derived for the anisotropic fabric from measured COF eigenvalues of the EDML ice core, Antarctica, show good agreement with the velocity trend determined from vertical seismic profiling. The agreement of the absolute velocity values, however, depends on the choice of the monocrystal elasticity tensor used for the calculation of the polycrystal properties. We make use of abrupt changes in COF as a common reflection mechanism for seismic and radar data below the firn–ice transition to determine COF-induced reflections in either data set by joint comparison with ice-core data. Our results highlight the possibility to complement regional radar surveys with local, surface-based seismic experiments to separate isochrones in radar data from other mechanisms. This is important for the reconnaissance of future ice-core drill sites, where accurate isochrone (i.e. non-COF) layer integrity allows for synchronization with other cores, as well as studies of ice dynamics considering non-homogeneous ice viscosity from preferred crystal orientations.

Description: The widely used detailed SNOWPACK model has undergone constant development over the years. A notable recent extension is the introduction of a Richards equation (RE) solver as an alternative for the bucket-type approach for describing water transport in the snow and soil layers. In addition, continuous updates of snow settling and new snow density parameterizations have changed model behavior. This study presents a detailed evaluation of model performance against a comprehensive multiyear data set from Weissfluhjoch near Davos, Switzerland. The data set is collected by automatic meteorological and snowpack measurements and manual snow profiles. During the main winter season, snow height (RMSE: 〈 4.2 cm), snow water equivalent (SWE, RMSE: 〈 40 mm w.e.), snow temperature distributions (typical deviation with measurements: 〈 1.0 °C) and snow density (typical deviation with observations: 〈 50 kg m−3) as well as their temporal evolution are well simulated in the model and the influence of the two water transport schemes is small. The RE approach reproduces internal differences over capillary barriers but fails to predict enough grain growth since the growth routines have been calibrated using the bucket scheme in the original SNOWPACK model. However, the agreement in both density and grain size is sufficient to parameterize the hydraulic properties successfully. In the melt season, a pronounced underestimation of typically 200 mm w.e. in SWE is found. The discrepancies between the simulations and the field data are generally larger than the differences between the two water transport schemes. Nevertheless, the detailed comparison of the internal snowpack structure shows that the timing of internal temperature and water dynamics is adequately and better represented with the new RE approach when compared to the conventional bucket scheme. On the contrary, the progress of the meltwater front in the snowpack as detected by radar and the temporal evolution of the vertical distribution of melt forms in manually observed snow profiles do not support this conclusion. This discrepancy suggests that the implementation of RE partly mimics preferential flow effects.

Description: We investigate the propagation of seismic waves in anisotropic ice. Two effects are important: (i) sudden changes in crystal orientation fabric (COF) lead to englacial reflections; (ii) the anisotropic fabric induces an angle dependency on the seismic velocities and, thus, recorded traveltimes. Velocities calculated from the polycrystal elasticity tensor derived for the anisotropic fabric from measured COF eigenvalues of the EDML ice core, Antarctica, show good agreement with the velocity trend determined from a vertical seismic profiling. The agreement of the absolute velocity values, however, depends on the choice of the monocrystal elasticity tensor used for the calculation of the polycrystal properties. With this validation of seismic velocities we make use of abrupt changes in COF as common reflection mechanism for seismic and radar data below the firn–ice transition to investigate their occurrence by comparison with ice-core data. Our results highlight the possibility to complement regional radar surveys with local, surface-based seismic deployment to separate isochrones in radar data from other mechanisms. This is important for the reconnaissance of future ice-core drill sites, where accurate isochrone (i.e. non-COF) layer integrity allows for synchronization with other cores, as well as studies of ice dynamics considering non-homogeneous viscosity from preferred crystal orientations.

Description: A preferred orientation of the anisotropic ice crystals influences the viscosity of the ice bulk and the dynamic behaviour of glaciers and ice sheets. Knowledge about the distribution of crystal anisotropy, to understand its contribution to ice dynamics, is mainly provided by crystal orientation fabric (COF) data from ice cores. However, the developed anisotropic fabric does not only influence the flow behaviour of ice, but also the propagation of seismic waves. Two effects are important: (i) sudden changes in COF lead to englacial reflections and (ii) the anisotropic fabric induces an angle dependency on the seismic velocities and, thus, also recorded traveltimes. A framework is presented here to connect COF data with the elasticity tensor to determine seismic velocities and reflection coefficients for cone and girdle fabrics from ice-core data. We connect the microscopic anisotropy of the crystals with the macroscopic anisotropy of the ice mass, observable with seismic methods. Elasticity tensors for different fabrics are calculated and used to investigate the influence of the anisotropic ice fabric on seismic velocities and reflection coefficients, englacially as well as for the ice-bed contact. Our work, therefore, provides a contribution to remotely determine the state of bulk ice anisotropy.

Description: A preferred orientation of the anisotropic ice crystals influences the viscosity of the ice bulk and the dynamic behaviour of glaciers and ice sheets. Knowledge about the distribution of crystal anisotropy is mainly provided by crystal orientation fabric (COF) data from ice cores. However, the developed anisotropic fabric influences not only the flow behaviour of ice but also the propagation of seismic waves. Two effects are important: (i) sudden changes in COF lead to englacial reflections, and (ii) the anisotropic fabric induces an angle dependency on the seismic velocities and, thus, recorded travel times. A framework is presented here to connect COF data from ice cores with the elasticity tensor to determine seismic velocities and reflection coefficients for cone and girdle fabrics.We connect the microscopic anisotropy of the crystals with the macroscopic anisotropy of the ice mass, observable with seismic methods. Elasticity tensors for different fabrics are calculated and used to investigate the influence of the anisotropic ice fabric on seismic velocities and reflection coefficients, englacially as well as for the ice–bed contact. Hence, it is possible to remotely determine the bulk ice anisotropy.

Description: The breakup of Gondwana is manifested by coeval early Jurassic Karoo magmatism in South Africa and East Antarctica. In South Africa, the large volumes of volcanic rocks of the adjoining Lebombo and Mwenetzi-Save monoclines represent a volcanic rift margin, and in East Antarctica, a corresponding feature, the Explora Wedge is buried below sediments and floating ice shelves on the continental margin of Dronning Maud Land. We use the seismic vibrator source to explore the sub-ice geology in Antarctica, and the new seismic reflection and available regional aeromagnetic data enable us to outline a dogleg landward extent of the Explora Wedge in Dronning Maud Land. The congruent inboard wedge geometries on the two continents define a high quality constraint, which facilitate for the first time, a geologically consistent and tight reconstruction of Africa relative to East Antarctica within Gondwana. The uncertainties in correlations of major geological features (mobile belts) from one continent to the other may now be of the order of ten’s of kilometers rather than hundreds of kilometers.