ALBERT

Description: Radio-echo sounding of the Antarctic and Greenlandic ice sheets often reveals a layer in the lowest hundreds of meters above bedrock more or less free of radio echoes, known as the echo-free zone (EFZ). The cause of this feature is unclear, so far lacking direct evidence for its origin. We compare echoes around the EPICA drill site in Dronning Maud Land, Antarctica, with the dielectric properties, crystal orientation fabrics and optical stratigraphy of the EPICA-DML ice core. We find that echoes disappear in the depth range where the dielectric contrast is blurred, and where the coherency of the layers in the ice core is lost due to disturbances caused by the ice flow. At the drill site, the EFZ onset at ~2100 m marks a boundary, below which the ice core may have experienced flow induced disturbances on various scales. The onset may indicate changing rheology which needs to be accounted for in the modeling of ice sheet dynamics.

Description: Radio-echo sounding of the Antarctic and Greenlandic ice sheets often reveals a layer in the lowest hundreds of meters above bedrock more or less free of radio echoes, known as the echo-free zone (EFZ). The cause of this feature is unclear, so far lacking direct evidence for its origin. We compare echoes around the EPICA drill site in Dronning Maud Land, Antarctica, with the microstructural and dielectrical properties of the EPICA-DML ice core. We find that echoes disappear in the depth range, where the coherency of the layers is lost due to disturbances caused by the ice flow. At the drill site, the EFZ onset at ~2100 m marks a boundary, below which the ice core may have experienced flow induced disturbances on various scales. The dating of the climate record becomes increasingly difficult below 1900 m, until correlation with the Dome C record is lost below 2417 m depth. The onset also indicates changing rheology which needs to be accounted for in the modeling of ice sheet dynamics.

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: 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 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 parametrisations have changed model behaviour. This study presents a detailed evaluation of model performance against a comprehensive multi-year 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:

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: The origin of a strong continuous radar reflector observed with airborne radio-echo sounding (RES) at the EPICA deep-drilling site in Dronning Maud Land, Antarctica, is identified as a transition in crystal fabric orientation from a vertical girdle- to increased single-pole orientation seen along the ice core. The reflector is observed with a 60 ns and 600 ns long pulse at a frequency of 150 MHz, spans one pulse length, is continuous over 5 km, and occurs at a depth of about 2020–2030 m at the drill site. Changes in conductivity as reflector origin are excluded by investigating the ice-core profile and synthetic RES data. Our observations allow to extrapolate the crystal orientation feature along the reflector in space, with implications for ice-sheet dynamics. As the conductivity profile of the EPICA shows no distinctive peak at this depths, we exclude changes in conductivity as the reflector origin. This is supported by application of numerical forward modelling of electromagnetic wave propagation, based on the conductivity profile, which is able to reproduce nearby reflections, but fails to reproduce this one. Because of background noise, the permittivity profile based on dielectric does not show prominent signals at these depths. We therefore interpret the observed reflector to originate from this change in crystal fabric.

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 origin of a strong continuous radar reflector observed with airborne radio-echo sounding (RES) at the EPICA deep-drilling site in Dronning Maud Land, Antarctica, is identified as a transition in crystal fabric orientation from a vertical girdle to an increased single-pole orientation seen along the ice core. The reflector is observed with a 60 ns and 600 ns long pulse at a frequency of 150 MHz, spans one pulse length, is continuous over 5 km, and occurs at a depth of about 2025–2045 m at the drill site. Changes in conductivity as reflector origin are excluded by investigating the ice-core profile, synthetic RES data, and a RES profile with different electromagnetic polarisation azimuths. The reflector's magnitude shows maximum values for polarisation parallel to the nearby ice divide and disappears for polarisation perpendicular to it, identifying the orientation of the girdle to lie in the vertical plane parallel to the ice divide. Observations allow us to extrapolate the crystal orientation feature along the reflector in space, with implications for ice-sheet dynamics and modeling.