Anisotropic backscatter in ice-penetrating radar data: potential mechanisms and implications

Airborne and ground-based radar have been used extensively in the past to measure ice thickness and to investigate the internal structure of ice sheets in terms of layering. The main reflection mechanisms for internal reflections are changes in density, conductivity, and crystal orientation fabric, which alter thepermittivity of the ice. Linking the different mechanisms to the individual reflection horizons enables thededuction of glaciological parameters like accumulation rates or age-depth estimates. If no sample material from snow pits or ice-cores are available, multi-frequency and multi-polarization measurements must be applied to distinguish between the different reflection mechanisms. The backscattered power of horizons caused by changes in conductivity varies with the center frequency whereas in the case of horizons originating from changing crystal orientation the backscattered power is dependent on the polarization plane of the carrier signal.In this study we examine a sample data set near the German summer station Kohnen (drill site for theEPICA-EDML ice core) on the Antarctic plateau. The data were acquired with an airplane sliding on ground, producing varying incident polarization with a circular profile and several cross profiles with different headings. We find that the backscattered power changes with varying antenna orientation (i.e. polarization). In the upper third of the ice column the backscatter has two maxima with a 180° symmetry. The maxima align with the direction of minimal surface strain. At approximately 900 m depth the anisotropy is shifted by 90° in heading azimuth, with the maxima now being parallel to the maximum in surface strain. This dataset is unique, as airborne systems (primarily designed for the sounding of ice thickness) are usually not used for ground-based applications. The observed anisotropy appears clearly and is intriguing as the reason for it is entirely unknown. As primary suspects we consider the role of changing crystal orientation and ellipsoidal shaped air bubbles. The effect is visible from 200 1400 m. It appears distributed along the entire interval, and not restricted to individual layers. It seems that the polarization dependence becomes visible by a changing background level of the acquired signal, which is otherwise largely dominated by layer-like, polarization independent reflections. Hence we apply a (semi-analytical) volume scattering model in order to understand the different reflection mechanisms better. From ice-core measurements it is known that the crystals in the upper hundred meters are only weakly aligned (if at all), and it is unclear how the crystal orientation changes overshort depth intervals (~10 m). The rotation of the anisotropy coincides with the clathrate transition in the ice core and thus we first focus on the effect of anisotropic air bubbles. In an in-coherent approach we treat the ice matrix as a random medium and use the vector radiative transfer theory to incorporate boundary conditions. In a second step we model the effect of crystal orientation to estimate both, the degree of alignment and the statistical variance in the permittivity tensor needed to generate the observed pattern in backscatter. Doing so, we eventually aim at pinning down the mechanisms for the anisotropy in the upper interval, lower interval and the interrelation of the two by a shift of 90°.Anisotropic air bubbles as well as aligned crystal orientation allow to deduce stress and strain rates and a potential change thereof along depth. So far it is largely unclear, how surface strain rates relate with strain rates within the ice. If one of the two suspected mechanisms can be excluded or confirmed, this study may serve as a case study for future polarimetric surveys with low-frequency radars, in order to supply ice-sheet modelling with adequate boundary conditions - including changes in the internal structure of ice sheets along depth

The trabecular architecture of the superior articular process of the lumbar spine (L2-S1)

The role of the facet joint in low back pain has gained public attention lately. The objective of our study was to investigate whether there is any difference in the adaptation of the cancellous bone in the superior articular process depending on the specific stress condition in different levels of the spine. Therefore, the trabecular structure of the superior articular processes of L2 and S1 of 15 cadavers (aged 63-100years) were studied using μCT (micro-computer tomography). Each sample was divided into five sections, each of which containing 20% of the slices. The following structure parameters were compared between L2 and S1 as well as within each process; bone-volume-fraction (BV/TV), trabecular number (Tb.N), trabecular thickness (Tb.Th), structure-model-index (SMI) and degree of anisotropy (DA). Statistically significant differences were observed between L2 and S1 for the BV/TV, SMI, Tb.Th and Tb.N in superior 2 sections. BV/TV, Tb.Th and Tb.N were higher in S1 than in L2. The SMI is lower, and even negative in S1 compared to L2, showing a more plate-like structure. Within the articular process all structure parameters show a similar distribution in L2 and S1. BV/TV, Tb.N and DA decreased from cranial to caudal while Tb.Th was highest in the most cranial and caudal sections, with the lowest value in the middle. The SMI, on the other hand, increased from cranial to caudal displaying more rod-like structures. These results can be explained by the different stress the processes of the different spinal levels are exposed to as well as the specific motion patterns of the facet joint. The processes of the os sacrum are exposed to a higher axial and ventral load due to their location and the lumbosacral flexion. In addition the upper sections of each process experience higher stress peaks than the lower ones. Therefore, this study shows the material distribution within the cancellous bone adapts to these specific stress conditions the facet joints are exposed t

Towards interpretation of the radio-stratigraphy of Antarctic ice shelves from modeling and observations: A case study for the Roi Baudouin Ice Shelf, East Antarctica

Ice shelves surrounding the Antarctic perimeter buttress ice flow from the continent towards the ocean, and their disintegration leads to an increase in ice discharge and sea level rise. The evolution and integrity of ice shelves is governed by surface accumulation, basal melting, and ice dynamics. We find history of these processes imprinted in the ice-shelf stratigraphy, which is mapped using isochrones imaged with radar. As an observational archive, the radar obtained stratigraphy combined with ice flow modeling has high potential to assist model calibration and reduce uncertainties in projections for the ice-sheet evolution. In this study we use a simplistic and observationally driven ice-dynamic forward model to predict the ice-shelf stratigraphy. We validate this approach with the full Stokes ice-flow model Elmer/Ice, and present a test-case for the Roi Baudouin Ice Shelf (East Antarctica) - where our model predictions agree well with radar obtained observations. The presented method enables us to investigate whether ice shelves are in steady-state, as well as to map spatial variations of how much of the ice-shelf volume is determined by its local surface mass balance. In the case of Roi Baudouin, we find the ice-shelf volume in the western part to be dominated by ice inflowing from the ice sheet, while the eastern part of the ice shelf is dominated by ice locally accumulated on the shelf. Such analysis serves as a metric for the susceptibility of ice shelves to climate change. We further apply our approach to other ice shelves in Antarctica

Constraining variable density of ice shelves using wide-angle radar measurements

International audienc

Predicting the steady-state isochronal stratigraphy of ice shelves using observations and modeling

Ice shelves surrounding the Antarctic perimeter moderate ice discharge towards the ocean through buttressing. Ice-shelf evolution and integrity depend on the local surface accumulation, basal melting and on the spatially variable ice-shelf viscosity. These components of ice-shelf mass balance are often poorly constrained by observations and introduce uncertainties in ice-sheet projections. Isochronal radar stratigraphy is an observational archive for the atmospheric, oceanographic and ice-flow history of ice shelves. Here, we predict the stratigraphy of locally accumulated ice on ice shelves with a kinematic forward model for a given atmospheric and oceanographic scenario. This delineates the boundary between local meteoric ice (LMI) and continental meteoric ice (CMI). A large LMI to CMI ratio hereby marks ice shelves whose buttressing strength is more sensitive to changes in atmospheric precipitation patterns. A mismatch between the steady-state predictions of the kinematic forward model and observations from radar can highlight inconsistencies in the atmospheric and oceanographic input data or be an indicator for a transient ice-shelf history not accounted for in the model. We discuss pitfalls in numerical diffusion when calculating the age field and validate the kinematic model with the full Stokes ice-flow model Elmer/Ice. The Roi Baudouin Ice Shelf (East Antarctica) serves as a test case for this approach. There, we find a significant east–west gradient in the LMI / CMI ratio. The steady-state predictions concur with observations on larger spatial scales (>10 km), but deviations on smaller scales are significant, e.g., because local surface accumulation patterns near the grounding zone are underestimated in Antarctic-wide estimates. Future studies can use these mismatches to optimize the input data or to pinpoint transient signatures in the ice-shelf history using the ever growing archive of radar observations of internal ice stratigraphy

Settings, current mass balance, and long-term evolution of ice rises in the Fimbul ice shelf, Western Dronning Maud Land

Ice rises, local ice caps on elevated ocean beds surrounded by ice shelves, impact the stability of the surrounding ice shelf and nearby outlet glaciers which, in turn, control the mass balance of the Antarctic Ice Sheet. In 2011, we started to investigate three ice rises in the vicinity of the Fimbul Ice Shelf, western Dronning Maud Land (DML), East Antarctica. Fimbul Ice Shelf is fed by Jutulstraumen outlet glacier that accounts for about 10% of the total outgoing mass from the DML sector of the Antarctic Ice Sheet. This mass flux is as large as that of Shirase and West Ragnhild Glaciers, two other most significant glaciers in DML. Our major findings include 1) low elevations of the ice-rise beds (100-350 m b.s.l.), 2) presence of englacial stacked upward arches in (Raymond Bumps) near the current ice-flow divides, 3) varying ice shelf thickness and surface velocity at the grounding zones on opposite sides of the ice rises, 4) orthographically influenced surface mass balance patterns, and 5) asymmetric surface flow speeds over the ice-flow divides. We use Results #4 and #5 to investigate recent (decadal) mass balance of the ice rises. Also, present-day thickness changes and ice rheology will be determined with repeated FM-CW radar measurements in the vicinity of the flow divides. These rheology estimates and Results #2, 3, 4, and 5 will be used to constrain thermo-mechanically coupled ice-flow models to examine timing of divide-flow onset, thinning/thickening rates (i.e. mass balance) of the ice rises and effects of changes of ice-shelf thicknesses at the end of the local flow fields on the ice-flow divide positions.第4回極域科学シンポジウム個別セッション：[OM] 気水圏11月15日（金） 統計数理研究所 ３階セミナー室１（D305

Quantifying the effect of ocean bed properties on ice sheet geometry over 40 000 years with a full-Stokes model

Simulations of ice sheet evolution over glacial cycles require integration of observational constraints using ensemble studies with fast ice sheet models. These include physical parameterisations with uncertainties, for example, relating to grounding-line migration. More complete ice dynamic models are slow and have thus far only be applied for  50 % under almost equal forcing. Grounding-line positions differ by up to 49 km, show significant hysteresis, and migrate non-steadily in both scenarios with long quiescent phases disrupted by leaps of rapid migration. The simulations quantify the evolution of two different ice sheet geometries (namely thick and slow vs. thin and fast), triggered by the variable grounding-line migration over the differing ocean beds. Our study extends the timescales of 3D full-Stokes by an order of magnitude compared to previous studies with the help of parallelisation. The extended time frame for full-Stokes models is a first step towards better understanding other processes such as erosion and sediment redistribution in the ice shelf cavity impacting the entire catchment geometry

High-resolution distributed vertical strain and velocity from repeat borehole logging by optical televiewer:Derwael Ice Rise, Antarctica

Abstract Direct measurements of spatially distributed vertical strain within ice masses are scientifically valuable but challenging to acquire. We use manual marker tracking and automatic cross correlation between two repeat optical televiewer (OPTV) images of an ~100 m-long borehole at Derwael Ice Rise (DIR), Antarctica, to reconstruct discretised, vertical strain rate and velocity at millimetre resolution. The resulting profiles decay with depth, from −0.07 a −1 at the surface to ~−0.002 a −1 towards the base in strain and from −1.3 m a −1 at the surface to ~−0.5 m a −1 towards the base in velocity. Both profiles also show substantial local variability. Three coffee-can markers installed at different depths into adjacent boreholes record consistent strain rates and velocities, although averaged over longer depth ranges and subject to greater uncertainty. Measured strain-rate profiles generally compare closely with output from a 2-D ice-flow model, while the former additionally reveal substantial high-resolution variability. We conclude that repeat OPTV borehole logging represents an effective means of measuring distributed vertical strain at millimetre scale, revealing high-resolution variability along the uppermost ~100 m of DIR, Antarctica.info:eu-repo/semantics/publishe