Research

I use a number of tools to help improve understanding of ice-sheet dynamics: idealized mathematical models; laboratory-scale fluid-mechanical experiments; and analysis of geophysical data, such as ice-surface velocities and ice-penetrating radar.

 

Henry Ice Rise:

The Henry Ice Rise is a grounded area of ice found within the floating Ronne Ice Shelf, in the Weddell Sea Sector of West Antarctica. Ice rises act as pinning points for ice shelves, providing resistance to flow and buttressing the ice upstream of them.

Ground-based ice-penetrating radar data has revealed a collection of relic crevasses at the base of the Henry Ice Rise, which could not be formed in the current flow configuration. Their presence suggests that the ice rise has grown from a previously smaller configuration and crevasses formed at a series of former grounding line positions have since been buried and deformed.

See the details of this in my poster from the AGU 2017 Fall Meeting in New Orleans below.

AGU_poster_2017

AGU_poster_2017

 

Analysis of Ice Surface Velocities:

Insights into the flow and dynamics of ice shelves can be sought through analysis of the surface velocity field. Using ice surface velocities, derived from interferometric synthetic aperture radar (InSAR) (reference: MEaSURES:  Rignot, E., Mouginot, J., & Scheuchl, B. (2011). Ice Flow of the Antarctic Ice Sheet. Science) strain rate fields can be generated and analysed. Figure (1) shows plots for the Fimbul Ice Shelf, Antarctica. Ice flows out from a narrow glacier at the lower boundary of the images and spreads out into a wide ice shelf. The most prominent areas of deformation are in the margins of this tongue of fast flow, with large strain and shear rates. These areas correspond to the largest principal strain rates, which are orientated at approximate 45 degrees to the main flow. From the visual MODIS images (Haran, T., Bohlander, J., Scambos, T., Painter, T., & Fahnestock, M. (2005). MODIS Mosaic of Antarctica 2003-2004 (MOA2004) Image Map.) it is clear that these areas also contain extensive crevassing.

Fimbul_strain_shear_P1_P2_modis
Figure (1): Fimbul Ice Shelf: (a) MODIS visual imagery, (b) speed, (c) along-flow strain rate, (d) shear rate, (e) 1st Principal Strain Axis, (f) 2nd Principal Strain Axis. Reference: Wearing, M. G. (2016). The Flow Dynamics and Buttressing of Ice Shelves. University of Cambridge.

 

 

 

Idealized mathematical models:

A large part of my PhD work consisted of using mathematical models to improve understanding of the flow and buttressing of ice shelves. Ice shelves are the floating extensions of ice sheets, which flow out over the ocean. They are often confined within embayments and provide resistance to the flow of grounded ice upstream, which is referred to as ice-shelf buttressing.

I considered the flow of ice shelves in parallel and linearly diverging channels. In these cases the flow of the ice shelf is mainly controlled by the shear stress between the lateral side walls, until a region near the ice front (the end of the channel), where both extensional and shear stresses control the dynamics.

Figure (2) below shows the analytical solution for a steady-state ice shelf under the assumption that flow is controlled by lateral shear. This can be compared to a numerical solution for an ice shelf with both extension and shear dynamics in Figure (3).

Steady_state_n3_parallel_3figs
Figure (2): Steady state; thickness profile, speed and along-flow strain rate for a shear-dominated ice shelf with Glen’s Flow Law rheology (power law) n=3. Reference: Wearing, M. G. (2016). The Flow Dynamics and Buttressing of Ice Shelves. University of Cambridge.

 

n3_shear-and-ext_Hss
Figure (3): Steady state; thickness profile, speed, along-flow strain rate and centreline strain rate for an ice shelf with shear and extensional dynamics. Reference: Wearing, M. G. (2016). The Flow Dynamics and Buttressing of Ice Shelves. University of Cambridge.

 

I have also used geophysical data in the form of ice surface velocities, strain rates, shelf widths and ice thickness to confirm an analytical scaling relationship for the flow of ice at a laterally confined calving front, where the ice shelf contains ice of a uniform rheology. This scaling relationship takes the form:

u∝ w(εH)^(3/4),  where u is the centreline velocity, w is the calving front width, ε is the along flow strain rate and H is the ice thickness, all evaluated at the calving front.

Figure (4a) shows a linear regression in log-log space, confirming the relationship determined analytically. Data points represented by purple crosses all lie above the regression line and correspond to ice shelves that are not laterally confined, or that have weak margins leading to greater than expected velocities at the calving front.

15J116_figure_8
Figure (4a): Linear regression in log-log space, with slope 0.72. Purple crosses represent ice shelves that are not laterally confined at calving front or that have damaged margins. Reference: Wearing, M. G., Hindmarsh, R. C. A., & Worster, M. G. (2015). Assessment of ice flow dynamics in the zone close to the calving front of Antarctic ice shelves. Journal of Glaciology.
15J116_figure_9
Figure (4b): Regression relationship plotted in linear space u∝ w (εH)^(0.72). Reference: Wearing, M. G., Hindmarsh, R. C. A., & Worster, M. G. (2015). Assessment of ice flow dynamics in the zone close to the calving front of Antarctic ice shelves. Journal of Glaciology.

For more details see:

https://www.cambridge.org/core/journals/journal-of-glaciology/article/assessment-of-ice-flow-dynamics-in-the-zone-close-to-the-calving-front-of-antarctic-ice-shelves/153C41279EDF977D2627CAA213D3DC01

pre-print pdf

 

 

Lab-scale fluid-mechanical experiments:

During my PhD at the Department of Applied Mathematics and Theoretical Physics at the University of Cambridge I undertook a series of fluid-mechanical experiments designed to simulate ice shelves confined within an idealized embayment. In these experiments parameters such as, channel width, length, fluid rheology and flux rate could be varied and directly compared with numerical models. Schematics of the experimental set-up are shown in Figures (5a) & (5b).

Ice_shelf_in_parallel_channel_diagram
Figure (5a): Experimental set-up side view
Ice_shelf_in_parallel_channel_diagram_plan_view
Figure (5b): Experimental set-up plan view

Figure (6) shows a side view of the floating current, flowing from right to left. The surface of the viscous current is marked by a dashed orange line, while the ocean surface is marked with a red dashed line. White tracer particles are distributed on the surface of the current and are used to infer the surface velocity field from images captured from a plan-view camera.

thickness_profile_example_image_annotated
Figure 6: Side view of experiment. Red dashed line denoted ocean surface, orange dashed line denotes viscous fluid surface.

Below (Fig. (7)) are some results from a Newtonian model ice shelf, where the viscous current is made of Golden Syrup. The red dashed line corresponds to the location of the grounding line, with the current flowing from top to bottom. A clear shear profile can be seen across the channel with no-slip at the side walls and maximum speed in the centre of the channel. In the upper two thirds of the channel the dynamics are dominated by shear and the along-flow strain rate is low. In the final section there is an increase in extensional strain, as predicted by theoretical models. This can be seen in both the plan view strain rate and centreline velocity.

gaus_velocity_5int_mm_walls
Figure (7): Flow field for a Golden Syrup (Newtonian) Ice Shelf. From left to right: Low-pass Gaussian filter velocity field; Time-averaged velocity field; Along-flow strain rate field; Centreline velocity. Red dashed lines represent position of grounding line. Yellow dashed line marks end of channel. Reference: Wearing, M. G. (2016). The Flow Dynamics and Buttressing of Ice Shelves. University of Cambridge.