Geophysical Research Abstracts Vol. 12, EGU2010-6973, 2010 EGU General Assembly 2010 Author(s) 2010 Three-dimensional modelling of calving processes on Johnsons Glacier, Livingston Island, Antarctica Jaime Otero (1), Francisco J Navarro (1), Carlos Martín (2), M Luisa Cuadrado (1), and M Isabel Corcuera (1) (1) Dept. Matemática Aplicada, ETSI detelecomunicación, Universidad Politécnica de Madrid, Spain (jotero@mat.upm.es), (2) British Antarctic Survey, High Cross, Madingley Road, Cambridge CB3 0ET, United Kingdom (cama@bas.ac.uk) Iceberg calving is an important mass loss mechanism from ice shelves and tidewater glaciers for many midand high-latitude glaciers and ice caps, yet the process is not well represented in prognostic models of ice dynamics. Benn and others (2007) proposed a calving criterion appropriate for both grounded and floating glacier tongues or ice shelves. This criterion assumes that the calving is triggered by the downward propagation of transverse surface crevasses, near the calving front, as a result of the extensional stress regime. The crevasse depth is calculated following Nye (1957), assuming that the base of a field of closely spaced crevasses lies at a depth where the longitudinal tensile strain rate tending to open the crevasse equals the creep closure resulting from the ice overburden pressure. Crevasses partially or totally filled with water will penetrate deeper, because of the contribution of water pressure to the opening of the crevasse. This criterion is readily incorporated into glacier and ice sheet models, but has not been fully validated with observations. We apply a three-dimensional extension of Benn and others (2007) criterion, incorporated into a full-stokes model of glacier dynamics, to estimate the current position of the calving front of Johnsons Glacier, Antarctica. We develop four experiments: (i) an straightforward three-dimensional extension of Benn and other s (2007) model; (2) an improvement to the latter that computes the tensile deviatoric stress opening the crevasse using the full-stress solution; (iii) a further improvement based on finding the depth at which the model-computed tensile deviatoric stress, considered as a function of depth, equals the ice overburden closure pressure; (iv) an experiment that adds, to the above, the effect of a threshold strain rate required for crevasses initiation. We found that the improvements considered in experiments (ii) and (iii) were necessary to reproduce accurately the observed calving front. Our modelling results also suggest that Johnsons Glacier has a polythermal structure, in contrast with the temperate structure suggested by earlier studies. REFERENCES: Benn, D.I., R.J. Hulton and R.H. Mottram. 2007a. Calving laws, sliding laws and the stability of tidewater glaciers. Ann. Glaciol., 46, 126-130. Nye, J.F. 1957. The distribution of stress and velocity in glaciers and ice-sheets. Proc. Roy. Soc., Ser. A, 239(1216), 113-133.
Three-dimensional modelling of calving processes on Johnsons glacier, Livingston Island, Antarctica Otero J. (1), Navarro F.J. (1), Martin C. (2), Cuadrado M.L. (1) and Corcuera M.I. (1) (1) Universidad Politecnica de Madrid (2) British Antartic Survey May 6, 2010 EGU General Assembly Vienna, Austria
Contents 1 Motivation The calving problem 2 Model of crevasse formation 3 A 3D calving model 3D extension of crevasse formation model Dynamical model Boundary conditions 4 Application to Johnsons glacier 5 Conclusions and outlook
Motivation The calving problem Most used models: Water depth at front. Problems: Highly empirical - distinct for different glaciers. Does not describe the physical mechanism. Height above buoyancy. Problem: Does not allow floating ice tongues/shelves. Other models: Force imbalance at terminal ice cliffs. Undercutting by subaqueous melting. Torque arising from buoyant forces. A recent model: Model of crevasse formation (Benn et al., 2007). Advantages: Strongly based on physics, allowing its use as pronostic model. Allows the development of floating ice tongues.
Model of crevasse formation Basics. Benn et al. (2007) x = L for d(x) = h(x) [ d = 1 ( ) 1 ] ǫ n 2 + (ρw gd w ) ρ i g A tensile stress = ice overburden pressure (opening) (closing) d d w h d 0 = 2 ρ i g ( ) 1 ǫ n A
Model of crevasse formation Glacier front evolution Z h X Dw L0 Z Z X X L 1 L 0 L0 L 1
Model of crevasse formation Limitations and extension Present limitations of crevasse formation model Bidimensional So far applied together with simple dynamical models: Driving stress balanced by basal drag Driving stress balanced by lateral drag A combination of both Our contribution 3D extension Application using a full-stokes dynamical model Improved crevasse depth computation
3D calving model Key aspects L becomes a function of x,y. ǫ ( = ǫ xx = u x direction in 2D model ) now becomes strain rate along ice flow ǫ u 2 u 1 x 2 x 1. ǫ determined from velocity field solution of a 3D full-stokes model. Improved crevasse depth computation Experiments.
3D calving model Experiments Experiment 1. Benn s model d 0 = 2 ( ǫ ρ i g A Without water Experiment 2 ) 1 n Tensile deviatoric stress calculated directly from constitutive relation [ d = 1 ( ) 1 ] ǫ n 2 + (ρw gd w ) ρ i g A With water filling the crevasses d 0 = 1 ρ i g B ǫ ǫ 1 n 1 Experiment 3 Exp. 2 + Tensile deviatoric stress as a function of depth Experiment 4 Exp. 3 + Introduce a yield strain rate
Dynamical model equations σ ij x j + ρg i = 0 conservation of linear momentum u i x i = 0 conservation of mass ǫ ij = Aτ n 1 τ ij constitutive relation
Boundary conditions XZ ice divide u v 0 YZ 0 Z accumulation stress free n 0 ij ij ablation calving u b basal sliding p b K p q e X
Boundary conditions Basal sliding Present implementation (1st approach) u b = K (ρgh)p p 1 h h x (ρgh ρ w gh w ) q, v b = K (ρgh)p p 1 h h y (ρgh ρ w gh w ) q, b w b = u b x + v b b y 120 100 60 80 d front,point d front,ela ELA d d ela,point ela,summit SUMMIT Hw ab = C d front,point + abs(min(0,b)) d front,ela H ac w = C(1 d ela,point d ela,summit ) Altitude a.s.l. (m) 300 200 100 ELA 0 200 400 600 800 1000 1200 1400 1600
Application: Johnsons glacier Location of Livingston Island and Hurd Peninsula, and surface map of Johnsons and Hurd glaciers (South Shetland Islands, Antartica). Limitations: Short time series of front position (5 yr.). Flat slope of seabed in proglacial area. Resulting in nearly constant front position. Application restricted to estimating model-predicted front position.
200 Surface topography and ELA 3050000 3049500 3049000 3048500 160 200 140 EJ29 120 100 60 80 EJ06 EJ05 EJ11 EJ15 EJ17 EJ14 140 EJ18 EJ16 160 EJ23 EJ10 EJ09 EJ26 EJ27 EJ28 3048000 240 EJ04 EJ22 240 260 280 260 280 EJ03 200 EJ21 3047500 300 EJ24 240 260 3047000 320 3046500 633500 634000 634500 635000 635500 636000 636500 637000 637500
Radar profiles, bathymetry and subglacial relief map 3050000 3049500 100 20 3049000 60 60 100 3048500 140 20 3048000 60 60 100 140 100 3047500 140 140 3047000 3046500 633500 634000 634500 635000 635500 636000 636500 637000 637500
200 Model parameters tuning. K uniform 0.32 a -1 E (m a ) b 1/3 B (Mpa a ) 0.3 0.28 0.26 0.24 0.22 0.2 9.2 9 8.8 8.6 8.4 8.2 8 7.8 7.6 7.4 7.2 7 6.8 6.6 6.4 6.2 6 5.8 5.6 5.4 5.2 5 Distance (km) 240 260 280 160 300 140 200 320 120 100 60 80 140 160 200 240 260 240 260 280 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2-1 -1 K (ma Pa ) Distance (km) Not in good agreement with observations near the glacier front.
200 Model parameters tuning. Nonuniform K K weighted according to position (2 zones, nearly coincident with accumulation & ablation). Aimed at closer agreement between computed and observed velocities. 0.28 0.27 a -1 E (m a ) b 7 1/3 B (Mpa a ) 0.26 0.25 0.24 0.23 6.8 6.6 6.4 6.2 6 5.8 5.6 5.4 5.2 5 4.8 Distance (km) 240 260 280 160 140 200 120 100 60 80 140 160 200 240 260 280 0.22 4.6 4.4 4.2 300 240 260 0.21 4 320 0.2 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2-1 -1 K (ma Pa ) Distance (km) Much better agreement with observations near the glacier front.
Predicted front position. Experiment 1 0 200 400 Water-free crevasses Water-filled crevasses (d w = H/2)
Predicted front position. Experiment 2 0 200 400 Water-free crevasses
Predicted front position. Experiment 3 0 200 400 Water-free crevasses Water-filled crevasses (d w = H/10)
Predicted front position. Experiment 4 0 200 400 Water-free crevasses Water-filled crevasses (d w = H/6)
Conclusions and outlook Conclusions Our three-dimensional extension of Benn s calving criterion, with Nye s formula, does not accurately reproduce the observed front position unless a large amount of water filling the near-front crevasses is hypothesized.
Conclusions and outlook Conclusions Our three-dimensional extension of Benn s calving criterion, with Nye s formula, does not accurately reproduce the observed front position unless a large amount of water filling the near-front crevasses is hypothesized. The modified criterion for crevasse depth, which computes the deviatoric longitudinal stress opening the crevasse using the full-stress solution, substantially improves the results. No water in crevasses is necessary. Crevasse depth slightly overestimated.
Conclusions and outlook Conclusions Our three-dimensional extension of Benn s calving criterion, with Nye s formula, does not accurately reproduce the observed front position unless a large amount of water filling the near-front crevasses is hypothesized. The modified criterion for crevasse depth, which computes the deviatoric longitudinal stress opening the crevasse using the full-stress solution, substantially improves the results. No water in crevasses is necessary. Crevasse depth slightly overestimated. The model that considers the tensile deviatoric stress opening the crevasse as a function of depth provides the best fit to observations with a small amount of water filling the crevasses.
Conclusions and outlook Conclusions Our three-dimensional extension of Benn s calving criterion, with Nye s formula, does not accurately reproduce the observed front position unless a large amount of water filling the near-front crevasses is hypothesized. The modified criterion for crevasse depth, which computes the deviatoric longitudinal stress opening the crevasse using the full-stress solution, substantially improves the results. No water in crevasses is necessary. Crevasse depth slightly overestimated. The model that considers the tensile deviatoric stress opening the crevasse as a function of depth provides the best fit to observations with a small amount of water filling the crevasses. Introducing a yield strain rate does not improve the fit to observations, unless a slightly larger amount of water is assumed. Such a yield strain rate is a physically plausible mechanism to overcome the fracture toughness of the ice.
Conclusions and outlook Outlook Improving the representation of basal sliding. Some alternatives: Improving the estimates of τb and H w. Application of the calving model to a glacier with a good record of front positions (e.g. Hansbreen), allowing use of transient model in prognostic mode.