Understanding Glacier-Climate Interaction with Simple Models
|
|
- Lionel Morrison
- 5 years ago
- Views:
Transcription
1 Understanding Glacier-Climate Interaction with Simple Models L / FS 203 Rhonegletscher, Switzerland, (VAW/ETHZ) A proglacial lake is forming in the terminus area of Rhonegletscher (Furkapass, Wallis). Well visible is the white tarp covering the touristic ice cave at the far side of the glacier. Martin Lüthi Versuchsanstalt für Wasserbau, Hydrologie und Glaziologie, VAW, ETH Zürich
2
3 Chapter Review of Glacier Dynamics This section is a short review on some essential aspects of glacier dynamics, which is mostly taken from the Physics of Glaciers I script. Good books on the topic include Cuey and Paterson (200) and Greve and Blatter (2009).. Field Equations To calculate velocities and stresses in a glacier we have to solve eld equations. For a mechanical problem (e.g. glacier ow) we need the continuity of mass and the force balance equations. Conservation of Mass The mass continuity equation for a compressible material of density is (in @z + r (v) = 0 = 0 (.a) (.b) If the density is i = 0) and constant (incompressible 0) we get, in dierent, equivalent notations Balance of Forces tr _" = _" ii = 0 r v = v i;i = 0 _" xx + _" yy + _" zz = 0 (.2a) (.2b) (.2c) (.2d) The force balance equation describes that all forces acting on a volume of ice, including the body force b = g (where g is gravity), need to be balanced by forces 3
4 Chapter Review of Glacier Dynamics acting on the sides of the volume. In compact tensor notation they read r + b = 0 ; (.3a) The same equations rewritten in index notation (summation convention) and in full, unabridged notation ij;j + b i j + b i = 0 xy yy + g x = 0 + g y = 0 (.3c) + g z = 0 These three equations describe how the body forces and boundary stresses are balanced by the stress gradients throughout the body. Rheology of Polycrystalline Ice Glacier ice usually considered as a nonlinear viscous uid, or more specically a power law uid. The most widely used ow relation for glacier ice is (Glen, 952; Nye, 957) _" ij = A n (d) ij : (.4) with n 3. The rate factor A = A(T ) depends on temperature and other parameters like water content, impurity content and crystal size. The quantity is the second invariant of the deviatoric stress tensor and dened by 2 := 2 (d) ij (d) ij = 2 (d)2 xx + (d)2 yy + (d)2 zz.2 Geometry Evolution + 2 (d)2 xy + 2 (d)2 xz + 2 (d)2 yz : (.5) If the surface of a glacier is described by a function F = 0, the following + v rf = _ b (.6) The net balance rate _ b describes how much volume is added or removed from the surface, and has units m/a. 4
5 Understanding Glacier-Climate Interaction with Simple Models If the surface is written as a function z s of the horizontal coordinates F (x; y; z; t) := z s (x; y; t) z = 0, this + s + s v z = _ b : (.7).3 Shallow Ice Approximation The vertical velocity distribution in a parallel sided slab n+ is u(z) = v x (z) = 2A n + (g sin )n H n+ (H z) {z } deformation velocity + {z} u b sliding velocity (.8) This is known as shallow ice equation, since it can be shown by rigorous scaling arguments that the longitudinal stress i i are negligible compared to the shear stress for shallow ice geometries such as the inland parts of ice sheets (except for the domes). 5
6
7 Chapter 2 Macroscopic Glacier Models Deep insights into the working of glaciers can obtained by omitting details of glacier dynamics, and concentrating on the main aspects of glacier response to climate. To this aim we will consider dierent Macroscopic Glacier Models. The idea is to simplify almost all aspects of glacier geometry and dynamics to elucidate the main processes. 2. Assumptions The volume V of a glacier depends on length L, width W and thickness H, and a shape factor f describing the fraction of the box W HL occupied by ice. It can be written V = f W HL : (2.) To treat glaciers as macroscopic entities means that one assumes the following Assumption : Details about ice thickness distribution are ignored. Assumption 2: Glaciers of dierent sizes are similar in shape. Þ Þ ¼ Þ Ä Î À Ä Î ¼ Ä Ü Figure 2.: Schematic geometry of a macroscopic glacier, and some important quantities. Assumption describes that one is only interested in total ice volume. Details of its distribution are lumped into the variables H (thickness) and f (shape). Assumption 7
8 Chapter 2 Macroscopic Glacier Models 2 describes an observation made in nature and from numerical models, that shapes of glaciers in steady state are often quite similar. Most glaciers have a characteristic longitudinal prole at the terminus. This similarity is (partially) captured in the shape factor f. The quantity H is a characteristic ice thickness scale which characterizes the whole glacier Elevation (m) Elevation (m) Distance (km) Distance (km) Figure 2.2: Plots illustrating the similarity in shape for dierent glaciers on a simple bedrock. Left: Surface geometries for dierent equilibrium line altitudes between Z = 00 : : : 600 m (indicated below the glacier terminus). Right: Surface geometries for Z = 400 m and dierent mass balance gradients (indicated below the glacier terminus). The dashed line shows the glacier geometry for Z = 500 m and _g = 0:006 a for comparison. Symbols on horizontal lines: Solid dots indicate the horizontal location of the equilibrium line, upward pointing triangles the maximum ice thickness, and downward pointing triangles maximum ow velocities (from Lüthi, 2009). Further simplications of glacier geometry are Assumption 3: Glacier width W is constant. Assumption 4: Bed slope s := tan is constant. Assumption 5: The volumetric balance rate _ b(t) is given by a function that depends only on position x or on elevation z. (The volumetric balance rate has units units m 3 =( m 2 a) = m a.) balance rate B _ is simply the integral of b B(t) = Z L 0 _b(x; t) dx or _B(t) = Z L The volume change rate is given by the total mass balance rate 0 Glacier net _b(z s (x); t) dx : (2.2) dv dt = _ B: (2.3) 8
9 Understanding Glacier-Climate Interaction with Simple Models Assumption 6: Glacier length L(t) changes according to an evolution equation. There are two very dierent cases: either length L is proportional to V (i.e. immediate reaction), or an independent evolution equation is prescribed dl dt = ( HW dv dt a (L L a ) L proportional to V relaxation equation (2.4) The rst equation describes immediate reaction with a proportionality between L and V. The second equation is a relaxation relation with a as the relaxation time scale (the limiting case a! 0 is equivalent to immediate reaction). The climateadjusted length L a can be calculated by either using a volume-area scaling relation, such that L a = L a (V ) (e.g. Lüthi, 2009), or by linking it directly to climate through a climate sensitivity factor c L a = cz (e.g. Klok and Oerlemans, 2003). To proceed, one has to make an assumption about macroscopic rheology. This essentially amounts to prescribing how ice thickness H depends on the state of the glacier. Assumption 7: Characteristic ice thickness H is due to a macroscopic rheology. Three approaches have been used that are discussed in more detail below H =8 > < >: H 0 constant thickness; 0 =(gs) perfect plasticity; a L scaling relation (self-similarity): (2.5a) The constant thickness case is the simplest of all, and has been employed by Harrison (203) to investigate large changes. Perfect plasticity and scaling have been investigated by Messié (20). The models by Harrison (203) (constant thickness) and Klok and Oerlemans (2003) implicitly belong to the last case (with = and = 2, respectively). 2.2 Mass balance / Volume change Glaciers change their volumes through mass balance processes at the surface (mainly, there are more processes within and under the ice). Alpine type glaciers have a mostly elevation-dependent mass balance, due to the atmospheric temperature lapse rate. In what follows we assume a linear increase of net mass balance rate with elevation. Net balance at the equilibrium line altitude (ELA) is zero (by denition). Assumption 5b: Net balance rate is given by _ b = _g(z s (x) z ELA ) 9
10 Chapter 2 Macroscopic Glacier Models The mass balance gradient is dened as _g := d _ b=dz. For a constant mass balance gradient _g we can write total balance rate as _B(t) = W = W = _g W Z L Z0 L Z 0 L 0 _b(z s (x); t) dx _g (z s (x) z ELA ) dx using z s (x) = z b (x) + h(x) (z b (x) + h(x) z ELA ) dx (2.6) To simplify this expression we dene Z := z 0 z ELA and use the denition V = W Z L 0 h(x) dx = flhw to arrive at _B = _g [V + (z b z ELA ) LW ] : (2.7) In a steady state the total mass balance rate is _ B = 0 and thus (with bars denoting averages over the glacier length) z s = z ELA and V = (z ELA z b ) LW = (z s z b ) LW ; (2.8) and therefore H = z s z b (obviously). To further simplify things, we assume a constant bedrock slope s Assumption 4: Bed slope s := tan is constant. The bed elevation therefore is z b (x) = z 0 sx and mean bedrock elevation for constant slope is z b = z 0 sl=2, which leads to V = Wz ELA z 0 + s 2 L L = W s 2 L2 ZL ; (2.9) where Z := z 0 z ELA has been introduced for convenience (actually, this amounts to redening the vertical coordinate, cf. Fig. 2.). This Equation shows that for the assumed simplications, glacier volume depends non-linearly upon glacier length, and that mean ice thickness is given by (Oerlemans, 200) H = fh = Now we can also write total balance in the simplied form _B = _g W(Z + fh)l V W L = s 2 L Z : (2.0) 0 L2 s : (2.) 2
11 Understanding Glacier-Climate Interaction with Simple Models 2.3 V-model (block model) First, we discuss the V-model (Harrison, 203). This model has the characteristics of immediate reaction of length to changes in volume (assumption 6a), and constant ice thickness H 0 (assumption 7a, thus f = ). The geometry is just a block of ice, which gives rise to the name block model. Since length and volume are related through V = W LH 0 ; (2.2) the volume evolution can be written by replacing L in Equation (2.) dv dt = B _ = _g (Z + H 0 ) V s V : (2.3) H 0 2W H 0 2! This equation can be simplied by scaling the variables. We use the denitions P := + Z H 0 V b := 2W H 2 0 s t b := _g equilibrium line parameter volume scale, i.e. V? = V=V b time scale, i.e. t? = t=t b Denoting scaled variables with asterisk (?), Equation (2.3) reads dv? dt? = P V? V?2 = V? (P V? ) : (2.4) Equation (2.4) is the logistic equation describing resource-limited growth (also called the Verhulst equation). For a constant climate (i.e. P = const) and an initial volume V it has the analytic solution?0 V? (t? ) = + P V?0 This function is also known as the logistic function. P e P t? : (2.5)
12 Chapter 2 Macroscopic Glacier Models 2.4 Volume-Length Scaling The V-model has a major shortcoming: the thickness does not vary with changing glacier length. A straightforward way to remedy this is to use a more physical approximation for ice mechanics. The V-model is reminiscent of perfect plasticity, i.e. the ice either deforms without limit if the shear stress at the base b exceeds a certain threshold 0, or it does not deform at all. Glacier ice behaves dierently, it is well described as a viscous uid, with a stress-dependent viscosity (or alternatively: a strain-rate dependent viscosity). It belongs to the class of power-law uids, often encountered in geophysics; for details refer to the Physics of Glaciers I script, or Cuey and Paterson (200). The only formula we make use of is the ice ux through a vertical section of an inclined, parallel-sided slab. With a thickness H and surface inclination s s = dzs dx = s + dh dx, the shallow ice approximation (i.e. neglecting longitudinal stress gradients; Hutter (983)) yields gss Q(H) = 2 _" 0 (n + 2) 0 n H n+2 = _" 0 s n s H n+2 ; (2.6) where g and denote gravitational acceleration and ice density, and = 2: m 2. Equation (2.6) is the integral of Equation (.8). Using standard values for temperate ice n = 3, A = 25 MPa 3 a (Paterson, 999) and a reference stress 0 = MPa = 0 5 Pa leads to a reference strain rate _" 0 = 25 a. Total mass balance rate in the accumulation area depends on the volume, which we denote by V A = f A W GH, where the shape factor f A describes by how much the accumulation area deviates from a block shape. Similar to Equation (2.) we get _B A = _gw(z + f A H)) G + G2 s 2 = _gw ( f A )ZG + f A 2 sg 2 (2.7a) ; (2.7b) where use was made of the denition of the ice thickness H at the equilibrium line Z = sg + H =) H = sg Z : (2.8) In a steady state the ice ux through a vertical section at the equilibrium line equals the total mass balance of the accumulation area. Therefore, combining Equations (2.6) and (2.7b) evaluated at x = G and for n = 3, yields (Lüthi, 2009) H 5 = _g _" 0 s 3 s G ( f A )Z + f A 2 sg : (2.9) Position G and ice thickness at the equilibrium line H can now be determined numerically from the system of equations (2.8) and (2.9) for a given value of Z. 2
13 Understanding Glacier-Climate Interaction with Simple Models For the special case s > 0, Z = 0 (ELA at highest elevation of the bed), this system of equations has the unique solution s H = sg = _g f A s _" 0 s 2 3 : (2.20) For a general analytic approximation, we assume s > 0 and Z > 0 which corresponds to a mountain glacier geometry. To derive a volumelength scaling relation one has to make assumptions about the surface slope and the mean ice thickness of the accumulation area. Results from Full-Stokes-model runs (FS-model) show that the mean ice thickness in the accumulation area is similar to the ice thickness at the equilibrium line, i.e. f A 0:95. Under the additional assumption that the surface slope at the equilibrium line is equal to the bed slope (s s s), Equation ' (2.9) yields _g 5 2 H G 2 _" 0 s 2 5 : (2.2) Since the accumulation area ratio (AAR; the fraction of accumulation area to the total area) is almost constant at r 0:53 for all steady FS-model glaciers, G can be replaced by rl, leading to ' _g H 2 _" 0 f _g V ' 2 _" 0 5 r s2 5 L 2 5 ; (2.22) 5 r s2 5 L 7 5 = al : (2.23) The rst expression in parentheses depends on the parameters for mass balance and ice deformation, the second on geometry. The scaling exponent = :4 is similar to.36 obtained by Bahr et al. (997), although for a dierent geometric setting. Moreover the dependence of the scaling factor a on mass balance gradient _g, ice ow rate factor _" 0, and bedrock slope s is explicitly given in Equation (2.23). 3
14 Chapter 2 Macroscopic Glacier Models 2.5 LV-model The LV-model is a simplied representation of glacier dynamics as a two-variable dynamical system in the variables length L and volume V (Lüthi, 2009). The dynamical system is formulated for unit width, and reproduces the essential inuence of mass balance and ice dynamics on glacier geometry on a macroscopic scale. Figure 2. illustrates the building blocks of the dynamical system: two reservoirs of volumes V A and V B which are linked by a ux element located at horizontal coordinate G. The resulting dynamical system is (Lüthi, 2009, Eq. 40) dv hv dt = _g " + ZL dl dt = V a a L2i s # ; (2.24a) 2 L : (2.24b) Equation (2.24a) is the same as Equation (2.3) in the V-model. Equation (2.24b) is a relaxation equation for the current glacier length L with time constant a (in years; a should not be confounded with the volume time scale v ). The steady state length for the current volume V is determined by the volume-length scaling relation given in Equation (2.23). This dynamical system (2.24) is driven by a forcing in the term Z(t) = z 0 z ELA (t). 4
15 Understanding Glacier-Climate Interaction with Simple Models List of Symbols _b volumetric mass balance rate m a _g vertical gradient of mass balance rate _g b _ _B total mass balance rate m 3 a s bed slope s = tan H glacier thickness m L glacier length m V glacier volume m 3 W glacier width m Z elevation span of accumulation area Z = z 0 z ELA m Z? elevation span of ablation area Z? = z ELA z L m x horizontal along-ow coordinate axis m y horizontal across-ow coordinate axis m z vertical coordinate axis m 5
16
17 Bibliography Bahr, D. B., Meier, M. F., and Peckham, S. D. (997). The physical basis of glacier volume-area scaling. Journal of Geophysical Research, 02(B9): Cuey, K. and Paterson, W. (200). The Physics of Glaciers. Elsevier. ISBN Glen, J. W. (952). Experiments on the deformation of ice. Journal of Glaciology, 2(2):4. Greve, R. and Blatter, H. (2009). Dynamics of Ice Sheets and Glaciers. Springer- Verlag Heidelberg. ISBN Harrison, W. (203). How do glaciers respond to climate? submitted. Journal of Glaciology. Hutter, K. (983). Theoretical glaciology; material science of ice and the mechanics of glaciers and ice sheets. D. Reidel Publishing Company/Tokyo, Terra Scientic Publishing Company. Klok, E. and Oerlemans, J. (2003). Deriving historical equilibrium-line altitudes from a glacier length record by linear inverse modelling. The Holocene, 3(3): Lüthi, M. P. (2009). Transient response of idealized glaciers to climate variations. Journal of Glaciology, 55(93): Messié, C. (20). Macroscopic response of idealized glaciers to a change in mass balance. Master's thesis, VAW, ETH Zürich. Nye, J. F. (957). The distribution of stress and velocity in glaciers and ice-sheets. Proceedings of the Royal Society of London, Ser. A, 239(26):333. Oerlemans, J. (200). Glaciers and climate change. Balkema Publishers. Paterson, W. S. B. (999). The Physics of Glaciers. Butterworth-Heinemann, third edition. 7
Dynamics of Glaciers
Dynamics of Glaciers McCarthy Summer School 01 Andy Aschwanden Arctic Region Supercomputing Center University of Alaska Fairbanks, USA June 01 Note: This script is largely based on the Physics of Glaciers
More informationIce in the climate system. Summary so far. Today. The Cryosphere. 1. Climate history of the Earth. 2. Paleo observations (1)
Ice in the climate system 1. Climate history of the Earth 2. Paleo observations (1) 3. Paleo observations (2) 4. Ice ages 5. Climate sensitivity 6. Ice in the climate system Summary so far Radiation (Milankovitch
More informationIn this section, mathematical description of the motion of fluid elements moving in a flow field is
Jun. 05, 015 Chapter 6. Differential Analysis of Fluid Flow 6.1 Fluid Element Kinematics In this section, mathematical description of the motion of fluid elements moving in a flow field is given. A small
More informationContinuum mechanism: Stress and strain
Continuum mechanics deals with the relation between forces (stress, σ) and deformation (strain, ε), or deformation rate (strain rate, ε). Solid materials, rigid, usually deform elastically, that is the
More informationHow do glaciers respond to climate? Perspectives from the simplest models
Journal of Glaciology, Vol. 59, No. 217, 2013 doi: 10.3189/2013JoG13J048 949 How do glaciers respond to climate? Perspectives from the simplest models W.D. HARRISON Geophysical Institute, University of
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION DOI: 1.138/NGEO1218 Supplementary information Ice speed of a calving glacier modulated by small fluctuations in basal water pressure Shin Sugiyama 1, Pedro Skvarca 2, Nozomu Naito
More informationChapter 1 Fluid Characteristics
Chapter 1 Fluid Characteristics 1.1 Introduction 1.1.1 Phases Solid increasing increasing spacing and intermolecular liquid latitude of cohesive Fluid gas (vapor) molecular force plasma motion 1.1.2 Fluidity
More informationThe effect of bottom boundary conditions in the ice-sheet to ice-shelf transition zone problem
Journal of Glaciology, Vol. 53, No. 182, 2007 363 The effect of bottom boundary conditions in the ice-sheet to ice-shelf transition zone problem Alexander V. WILCHINSKY Centre for Polar Observation and
More informationESS314. Basics of Geophysical Fluid Dynamics by John Booker and Gerard Roe. Conservation Laws
ESS314 Basics of Geophysical Fluid Dynamics by John Booker and Gerard Roe Conservation Laws The big differences between fluids and other forms of matter are that they are continuous and they deform internally
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More information2 A: The Shallow Water Equations
2 A: The Shallow Water Equations 2.1 Surface motions on shallow water Consider two-dimensional (x-z) motions on a nonrotating, shallow body of water, of uniform density, as shown in Fig. 1 below. The ow
More information2 GOVERNING EQUATIONS
2 GOVERNING EQUATIONS 9 2 GOVERNING EQUATIONS For completeness we will take a brief moment to review the governing equations for a turbulent uid. We will present them both in physical space coordinates
More informationA posteriori estimator for the accuracy of the shallow shelf approximation
Math Geosci manuscript No. will be inserted by the editor A posteriori estimator for the accuracy of the shallow shelf approximation Guillaume Jouvet Marco Picasso Received: date / Accepted: date Abstract
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationESS 431 Principles of Glaciology ESS 505 The Cryosphere
MID-TERM November 9, 2015 ESS 431 Principles of Glaciology ESS 505 The Cryosphere Instructions: Please answer the following 5 questions. [The actual 5 questions will be selected from these 12 questions
More informationModelling of dispersed, multicomponent, multiphase flows in resource industries. Section 3: Examples of analyses conducted for Newtonian fluids
Modelling of dispersed, multicomponent, multiphase flows in resource industries Section 3: Examples of analyses conducted for Newtonian fluids Globex Julmester 017 Lecture # 04 July 017 Agenda Lecture
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationContinuum Model of Avalanches in Granular Media
Continuum Model of Avalanches in Granular Media David Chen May 13, 2010 Abstract A continuum description of avalanches in granular systems is presented. The model is based on hydrodynamic equations coupled
More informationChapter 3. Shallow Water Equations and the Ocean. 3.1 Derivation of shallow water equations
Chapter 3 Shallow Water Equations and the Ocean Over most of the globe the ocean has a rather distinctive vertical structure, with an upper layer ranging from 20 m to 200 m in thickness, consisting of
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationEvolution of Rhonegletscher, Switzerland, over the past 125 years and in the future: application of an improved flowline model
268 Annals of Glaciology 46 2007 Evolution of Rhonegletscher, Switzerland, over the past 125 years and in the future: application of an improved flowline model Shin SUGIYAMA, 1,2 Andreas BAUDER, 2 Conradin
More informationStress, Strain, Mohr s Circle
Stress, Strain, Mohr s Circle The fundamental quantities in solid mechanics are stresses and strains. In accordance with the continuum mechanics assumption, the molecular structure of materials is neglected
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationEART162: PLANETARY INTERIORS
EART162: PLANETARY INTERIORS Francis Nimmo Last Week Global gravity variations arise due to MoI difference (J 2 ) We can also determine C, the moment of inertia, either by observation (precession) or by
More informationThermal flow in glaciers: Application to the Lys Glacier (Italian Western Alps)
Thermal flow in glaciers: Application to the Lys Glacier (Italian Western Alps) A. Deponti, L. De Biase & V. Pennati Environmental Science Department, University of Milano Bicocca, Italy Abstract In the
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationAE/ME 339. Computational Fluid Dynamics (CFD) K. M. Isaac. Momentum equation. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.
AE/ME 339 Computational Fluid Dynamics (CFD) 9//005 Topic7_NS_ F0 1 Momentum equation 9//005 Topic7_NS_ F0 1 Consider the moving fluid element model shown in Figure.b Basis is Newton s nd Law which says
More information10 Shallow Water Models
10 Shallow Water Models So far, we have studied the effects due to rotation and stratification in isolation. We then looked at the effects of rotation in a barotropic model, but what about if we add stratification
More informationNDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.
CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo
More informationChapter 2: Basic Governing Equations
-1 Reynolds Transport Theorem (RTT) - Continuity Equation -3 The Linear Momentum Equation -4 The First Law of Thermodynamics -5 General Equation in Conservative Form -6 General Equation in Non-Conservative
More informationRock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth
Rock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth References: Turcotte and Schubert, Geodynamics, Sections 2.1,-2.4, 2.7, 3.1-3.8, 6.1, 6.2, 6.8, 7.1-7.4. Jaeger and Cook, Fundamentals of
More informationConvection Heat Transfer
Convection Heat Transfer Department of Chemical Eng., Isfahan University of Technology, Isfahan, Iran Seyed Gholamreza Etemad Winter 2013 Heat convection: Introduction Difference between the temperature
More information12. Stresses and Strains
12. Stresses and Strains Finite Element Method Differential Equation Weak Formulation Approximating Functions Weighted Residuals FEM - Formulation Classification of Problems Scalar Vector 1-D T(x) u(x)
More informationNumerical Modelling in Geosciences. Lecture 6 Deformation
Numerical Modelling in Geosciences Lecture 6 Deformation Tensor Second-rank tensor stress ), strain ), strain rate ) Invariants quantities independent of the coordinate system): - First invariant trace:!!
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
More informationKINEMATICS OF CONTINUA
KINEMATICS OF CONTINUA Introduction Deformation of a continuum Configurations of a continuum Deformation mapping Descriptions of motion Material time derivative Velocity and acceleration Transformation
More informationStress, Strain, and Viscosity. San Andreas Fault Palmdale
Stress, Strain, and Viscosity San Andreas Fault Palmdale Solids and Liquids Solid Behavior: Liquid Behavior: - elastic - fluid - rebound - no rebound - retain original shape - shape changes - small deformations
More information1 Stress and Strain. Introduction
1 Stress and Strain Introduction This book is concerned with the mechanical behavior of materials. The term mechanical behavior refers to the response of materials to forces. Under load, a material may
More informationViscous Fluids. Amanda Meier. December 14th, 2011
Viscous Fluids Amanda Meier December 14th, 2011 Abstract Fluids are represented by continuous media described by mass density, velocity and pressure. An Eulerian description of uids focuses on the transport
More informationAE/ME 339. K. M. Isaac Professor of Aerospace Engineering. 12/21/01 topic7_ns_equations 1
AE/ME 339 Professor of Aerospace Engineering 12/21/01 topic7_ns_equations 1 Continuity equation Governing equation summary Non-conservation form D Dt. V 0.(2.29) Conservation form ( V ) 0...(2.33) t 12/21/01
More informationSummer School in Glaciology, Fairbanks/McCarthy, Exercises: Glacial geology
Bob Anderson Summer School in Glaciology, Fairbanks/McCarthy, 2010 Exercises: Glacial geology 1. Glacier thickness. We wish to estimate the local thickness of a glacier given only a topographic map of
More informationGlacier Thermodynamics: Ice Temperature and Heat Transfer Processes
Glacier Thermodynamics: Ice Temperature and Heat Transfer Processes ESS431: Principles of Glaciology ESS505: The Cryosphere Wednesday, 10/24 Ben Hills Today s Objectives: Why do we care about ice temperature?
More informationNotes 4: Differential Form of the Conservation Equations
Low Speed Aerodynamics Notes 4: Differential Form of the Conservation Equations Deriving Conservation Equations From the Laws of Physics Physical Laws Fluids, being matter, must obey the laws of Physics.
More informationUQ Benchmark Problems for Multiphysics Modeling
SIAMUQ UQ Challenge Benchmarks UQ Benchmark Problems for Multiphysics Modeling Maarten Arnst March 31, 2014 SIAMUQ UQ Challenge Benchmarks 1 / 25 Motivation Previous presentation at USNCCM2013 UQ Challenge
More informationCAN GLACIER IN ICE CAVE CUT U-SHAPED VALLEY - A NUMERICAL ANALYSIS
CAN GLACIER IN ICE CAVE CUT U-SHAPED VALLEY - A NUMERICAL ANALYSIS Shaohua Yang No.19A Yuquan Road Beijing 100049, China, yangshaohua09@sina.com Yaolin Shi No.19A Yuquan Road Beijing 100049, China, shyl@ucas.ac.cn
More information3D Elasticity Theory
3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.
More informationP = 1 3 (σ xx + σ yy + σ zz ) = F A. It is created by the bombardment of the surface by molecules of fluid.
CEE 3310 Thermodynamic Properties, Aug. 27, 2010 11 1.4 Review A fluid is a substance that can not support a shear stress. Liquids differ from gasses in that liquids that do not completely fill a container
More informationUnit IV State of stress in Three Dimensions
Unit IV State of stress in Three Dimensions State of stress in Three Dimensions References Punmia B.C.,"Theory of Structures" (SMTS) Vol II, Laxmi Publishing Pvt Ltd, New Delhi 2004. Rattan.S.S., "Strength
More informationGlacier Dynamics. Glaciers 617. Andy Aschwanden. Geophysical Institute University of Alaska Fairbanks, USA. October 2011
Glacier Dynamics Glaciers 617 Andy Aschwanden Geophysical Institute University of Alaska Fairbanks, USA October 2011 1 / 81 The tradition of glacier studies that we inherit draws upon two great legacies
More informationNumerical simulation of the evolution of glacial valley cross sections
Numerical simulation of the evolution of glacial valley cross sections arxiv:0901.1177v1 [physics.geo-ph] 9 Jan 2009 Hakime Seddik Ralf Greve Shin Sugiyama Institute of Low Temperature Science, Hokkaido
More informationChapter 4: Fluid Kinematics
Overview Fluid kinematics deals with the motion of fluids without considering the forces and moments which create the motion. Items discussed in this Chapter. Material derivative and its relationship to
More informationGeophysical Ice Flows: Analytical and Numerical Approaches
Geophysical Ice Flows: Analytical and Numerical Approaches Will Mitchell University of Alaska - Fairbanks July 23, 2012 Supported by NASA grant NNX09AJ38G Ice: an awesome problem...velocity, pressure,
More information1 Glaciers and climate
1 Glaciers and climate 1.1 Equilibrium length For a glacier to be in equilibrium with climate the (integrated) mass balance must be equal to zero. For a glacier of a uniform width and length l this implies
More informationThe Shallow Water Equations
The Shallow Water Equations Clint Dawson and Christopher M. Mirabito Institute for Computational Engineering and Sciences University of Texas at Austin clint@ices.utexas.edu September 29, 2008 The Shallow
More informationOCN/ATM/ESS 587. The wind-driven ocean circulation. Friction and stress. The Ekman layer, top and bottom. Ekman pumping, Ekman suction
OCN/ATM/ESS 587 The wind-driven ocean circulation. Friction and stress The Ekman layer, top and bottom Ekman pumping, Ekman suction Westward intensification The wind-driven ocean. The major ocean gyres
More informationAdvanced Course in Theoretical Glaciology
Advanced Course in Theoretical Glaciology Ralf Greve Institute of Low Temperature Science Hokkaido University Lecture Notes Sapporo 2015 These lecture notes are largely based on the textbook Greve, R.
More informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More informationModule 2: Governing Equations and Hypersonic Relations
Module 2: Governing Equations and Hypersonic Relations Lecture -2: Mass Conservation Equation 2.1 The Differential Equation for mass conservation: Let consider an infinitely small elemental control volume
More informationCE 240 Soil Mechanics & Foundations Lecture 7.1. in situ Stresses I (Das, Ch. 8)
CE 240 Soil Mechanics & Foundations Lecture 7.1 in situ Stresses I (Das, Ch. 8) Class Outline Stress tensor, stress units Effective stress, Stresses in saturated soil without seepage Stresses in saturated
More informationn i,j+1/2 q i,j * qi+1,j * S i+1/2,j
Helsinki University of Technology CFD-group/ The Laboratory of Applied Thermodynamics MEMO No CFD/TERMO-5-97 DATE: December 9,997 TITLE A comparison of complete vs. simplied viscous terms in boundary layer
More informationSimulating a glacier: application to Engabreen
Utrecht University Institute for Marine and Atmospheric Research Utrecht Bachelor thesis Simulating a glacier: application to Engabreen Supervisors: Carleen Tijm-Reijmer Hans Oerlemans Author: Jesse Reusen
More informationSec. 1.1: Basics of Vectors
Sec. 1.1: Basics of Vectors Notation for Euclidean space R n : all points (x 1, x 2,..., x n ) in n-dimensional space. Examples: 1. R 1 : all points on the real number line. 2. R 2 : all points (x 1, x
More informationDynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet
Dynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet or all equations you will probably ever need Definitions 1. Coordinate system. x,y,z or x 1,x
More informationDarcy s law in 3-D. K * xx K * yy K * zz
PART 7 Equations of flow Darcy s law in 3-D Specific discarge (vector) is calculated by multiplying te ydraulic conductivity (second-order tensor) by te ydraulic gradient (vector). We obtain a general
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Surface and body forces Tensors, Mohr circles. Theoretical strength of materials Defects Stress concentrations Griffith failure
More informationChapter 5. Sound Waves and Vortices. 5.1 Sound waves
Chapter 5 Sound Waves and Vortices In this chapter we explore a set of characteristic solutions to the uid equations with the goal of familiarizing the reader with typical behaviors in uid dynamics. Sound
More informationStability of Shear Flow
Stability of Shear Flow notes by Zhan Wang and Sam Potter Revised by FW WHOI GFD Lecture 3 June, 011 A look at energy stability, valid for all amplitudes, and linear stability for shear flows. 1 Nonlinear
More informationStudy of Glacier flow. Glacier Flow and Ice Profiles 10/8/09. Downslope transport of ice mass Accumula'on vs abla'on
Glacier Flow and Ice Profiles What can we learn about glacier movement and subglacial processes? Study of Glacier flow Based on three lines of observamon: 1. Field evidence 2. Laboratory experiments 3.
More informationM E 320 Professor John M. Cimbala Lecture 10
M E 320 Professor John M. Cimbala Lecture 10 Today, we will: Finish our example problem rates of motion and deformation of fluid particles Discuss the Reynolds Transport Theorem (RTT) Show how the RTT
More informationModelling ice shelf basal melt with Glimmer-CISM coupled to a meltwater plume model
Modelling ice shelf basal melt with Glimmer-CISM coupled to a meltwater plume model Carl Gladish NYU CIMS February 17, 2010 Carl Gladish (NYU CIMS) Glimmer-CISM + Plume February 17, 2010 1 / 24 Acknowledgements
More informationThe Nature of Kinematic Waves in Glaciers and their Application to Understanding the Nisqually Glacier, Mt. Rainier, Washington
Portland State University PDXScholar University Honors Theses University Honors College 2016 The Nature of Kinematic Waves in Glaciers and their Application to Understanding the Nisqually Glacier, Mt.
More informationAssessing the ability of numerical ice sheet models to simulate grounding line migration
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 110,, doi:10.1029/2004jf000202, 2005 Assessing the ability of numerical ice sheet models to simulate grounding line migration A. Vieli and A. J. Payne Centre for Polar
More informationComputer Applications in Engineering and Construction Programming Assignment #9 Principle Stresses and Flow Nets in Geotechnical Design
CVEN 302-501 Computer Applications in Engineering and Construction Programming Assignment #9 Principle Stresses and Flow Nets in Geotechnical Design Date distributed : 12/2/2015 Date due : 12/9/2015 at
More informationLoading σ Stress. Strain
hapter 2 Material Non-linearity In this chapter an overview of material non-linearity with regard to solid mechanics is presented. Initially, a general description of the constitutive relationships associated
More informationPAPER 333 FLUID DYNAMICS OF CLIMATE
MATHEMATICAL TRIPOS Part III Wednesday, 1 June, 2016 1:30 pm to 4:30 pm Draft 21 June, 2016 PAPER 333 FLUID DYNAMICS OF CLIMATE Attempt no more than THREE questions. There are FOUR questions in total.
More informationNumerical simulations of glacier evolution performed using flow-line models of varying complexity
Numerical simulations of glacier evolution performed using flow-line models of varying complexity Antonija Rimac 1, 2, Sharon van Geffen 1, 2, and Johannes Oerlemans 1 1 Institute for Marine and Atmospheric
More informationBasic Equations of Elasticity
A Basic Equations of Elasticity A.1 STRESS The state of stress at any point in a loaded bo is defined completely in terms of the nine components of stress: σ xx,σ yy,σ zz,σ xy,σ yx,σ yz,σ zy,σ zx,andσ
More informationMECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso
MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Preliminary Math Concept of Stress Stress Components Equilibrium
More informationLecture 8: Tissue Mechanics
Computational Biology Group (CoBi), D-BSSE, ETHZ Lecture 8: Tissue Mechanics Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015/16 7. Mai 2016 2 / 57 Contents 1 Introduction to Elastic Materials
More informationP = ρ{ g a } + µ 2 V II. FLUID STATICS
II. FLUID STATICS From a force analysis on a triangular fluid element at rest, the following three concepts are easily developed: For a continuous, hydrostatic, shear free fluid: 1. Pressure is constant
More informationGlaciology (as opposed to Glacial Geology) Why important? What are glaciers? How do they work?
Glaciology (as opposed to Glacial Geology) Why important? What are glaciers? How do they work? Glaciers are important because of their role in creating glacial landscapes (erosional and depositional features).
More information6 Temperatures in glaciers and ice sheets
Chapter 6 Temperatures in glaciers and ice sheets Glaciers are divided into three categories, depending on their thermal structure Cold The temperature of the ice is below the pressure melting temperature
More informationCHAPTER 2 Pressure and Head
FLUID MECHANICS Gaza, Sep. 2012 CHAPTER 2 Pressure and Head Dr. Khalil Mahmoud ALASTAL Objectives of this Chapter: Introduce the concept of pressure. Prove it has a unique value at any particular elevation.
More informationT. Perron Glaciers 1. Glaciers
T. Perron 12.001 Glaciers 1 Glaciers I. Why study glaciers? [PPT: Perito Moreno glacier, Argentina] Role in freshwater budget o Fraction of earth s water that is fresh (non-saline): 3% o Fraction of earth
More informationConstitutive Equations
Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module
More informationContinuum Mechanics. Continuum Mechanics and Constitutive Equations
Continuum Mechanics Continuum Mechanics and Constitutive Equations Continuum mechanics pertains to the description of mechanical behavior of materials under the assumption that the material is a uniform
More informationDecay of the Greenland Ice Sheet due to surface-meltwater-induced acceleration of basal sliding
Decay of the Greenland Ice Sheet due to surface-meltwater-induced acceleration of basal sliding arxiv:0905.07v [physics.geo-ph] May 009 Ralf Greve Shin Sugiyama Institute of Low Temperature Science, Hokkaido
More informationStress equilibrium in southern California from Maxwell stress function models fit to both earthquake data and a quasi-static dynamic simulation
Stress equilibrium in southern California from Maxwell stress function models fit to both earthquake data and a quasi-static dynamic simulation Peter Bird Dept. of Earth, Planetary, and Space Sciences
More informationAdvanced Heat and Mass Transfer by Amir Faghri, Yuwen Zhang, and John R. Howell
Laminar external natural convection on vertical and horizontal flat plates, over horizontal and vertical cylinders and sphere, as well as plumes, wakes and other types of free flow will be discussed in
More informationMath Review: Vectors and Tensors for Rheological Applications
Math Review: Vectors and Tensors for Rheological Applications Presented by Randy H. Ewoldt University of Illinois at Urbana-Champaign U. of Minnesota Rheological Measurements Short Course June 2016 MathReview-1
More informationParameterization of lateral drag in flowline models of glacier dynamics
Journal of Glaciology, Vol. 58, No. 212, 2012 doi: 10.3189/2012JoG12J018 1119 Parameterization of lateral drag in flowline models of glacier dynamics Surendra ADHIKARI, Shawn J. MARSHALL Department of
More informationMAE 323: Lecture 1. Review
This review is divided into two parts. The first part is a mini-review of statics and solid mechanics. The second part is a review of matrix/vector fundamentals. The first part is given as an refresher
More informationMultiple Integrals and Vector Calculus: Synopsis
Multiple Integrals and Vector Calculus: Synopsis Hilary Term 28: 14 lectures. Steve Rawlings. 1. Vectors - recap of basic principles. Things which are (and are not) vectors. Differentiation and integration
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 12.510 Introduction to Seismology Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Stress and Strain
More informationq v = - K h = kg/ν units of velocity Darcy's Law: K = kρg/µ HYDRAULIC CONDUCTIVITY, K Proportionality constant in Darcy's Law
Darcy's Law: q v - K h HYDRAULIC CONDUCTIVITY, K m/s K kρg/µ kg/ν units of velocity Proportionality constant in Darcy's Law Property of both fluid and medium see D&S, p. 62 HYDRAULIC POTENTIAL (Φ): Φ g
More informationDynamics of Ice Sheets and Glaciers
Dynamics of Ice Sheets and Glaciers Ralf Greve Institute of Low Temperature Science Hokkaido University Lecture Notes Sapporo 2004/2005 Literature Ice dynamics Paterson, W. S. B. 1994. The Physics of
More informationChapter 3 Permeability
3.2 Darcy s Law In 1856, Darcy investigated the flow of water through sand filters for water purification. His experimental apparatus is shown in Figure 3.11. By empirical observation Figure 3.11 Schematic
More informationMulti-Modal Flow in a Thermocoupled Model of the Antarctic Ice Sheet, with Verification
Multi-Modal Flow in a Thermocoupled Model of the Antarctic Ice Sheet, with Verification Craig Lingle 1 Jed Brown 2 Ed Bueler 2 1 Geophysical Institute University of Alaska Fairbanks, USA 2 Department of
More information