Postglacial Rebound Modeling with Power-Law Rheology
|
|
- Brice Sutton
- 5 years ago
- Views:
Transcription
1 Postglacial Rebound Modeling with Power-Law Rheology Patrick Wu Dept. of Geology & Geophysics, University of Calgary, Calgary, Alberta T2N-1N4, Canada, Keywords: rheology, power law creep, relative sea levels, ambient tectonic stress Abstract. Recent progress in postglacial rebound modeling with power-law creep is reviewed. Results of earlier 2D finite element models with simple nonlinear rheology are extended to 3D models with realistic ice histories. In the absence of ambient tectonic stress, mantle with power law rheology is shown to be characterized by a localized flow with distinct uplift patterns inside and outside the former ice margin. An important implication of this is that relative sea-level observations in and around the centre of rebound in Laurentia cannot be simultaneously reconciled unless nonlinear rheology is limited to lie within a thin zone ( 2 km) in the upper mantle below the lithosphere. Results of 2D modeling shows that ambient tectonic stress have an important effect on the postglacial rebound process if the tectonic stress level exceeds 1 MPa. In the presence of large ambient stress, power law mantle loses its uplift characteristics inside the ice margin and thus the flow appears to be linear. However, the uplift characteristic of a power law medium outside the ice margin is still retained even with ambient stress as high as 1 MPa. Thus, in the presence of ambient stress, nonlinear mantles continue to have difficulty explaining the relative sea level data outside former Laurentia ice margin. 1. Introduction An important question in the study of the dynamics of the earth, is whether the flow law of the mantle is linear (Newtonian) or nonlinear (power-law). Creep experiments on rocks suggest that mantle deformation follows a linear creep law (diffusion creep) when the stress level is low or when the grain size is small. However, at high stress level or large grain size, the dominant creep mechanism is dislocation creep (i.e. power-law creep). It has been suggested that nonlinear rheology prevails in the upper part of the mantle [4]. However, due to the uncertainty of the transition conditions between diffusion and dislocation creep, it is not clear if the rest of the mantle is nonlinear, and if there is a transition in the dominant creep mechanism, where that may occur. Recently, Karato & Li [5] suggested that the lower mantle may not be nonlinear. Karato & Wu [6] showed that rheological transition may occur in the upper mantle and suggested that this question of rheological transition may be resolved with seismic anisotropy observations and postglacial rebound modeling. The purpose of this chapter is to review recent progress in postglacial rebound modeling with
2 power-law creep. In the next section, the method of calculation is briefly reviewed. In section 3 & 4, simple disc load calculations are presented. The characteristics of power-law halfspaces are discussed in section 3, and in the section 4, relative sea level data is used to infer whether the whole mantle is nonlinear. The depth of the rheological transition zone is also addressed. In section 5, results for a realistic 3D ice deglaciation model is considered. In section 6, the effect of ambient tectonic stress is investigated and the results are summarized the section The Model In the early studies of postglacial rebound with nonlinear rheology, simplifying assumptions were made to keep the calculation tractable: Post and Griggs [9], and later Crough [2], assumed that the surface depression kept a constant shape during the rebound process; Brennen [1] assumed that the strain rate within the earth decreases exponentially with depth; Yokokura & Saito [2] assumed that the relaxation time only depends on time and not on space; Nakada [8] assumed a shallow channel flow. Unfortunately, all these simplifying assumptions have been shown to be invalid [16], calling into question the results of these earlier studies. In this paper, 2D and 3D Finite Element Models will be used to calculate the deformation of the earth due to glacial loading and unloading. The advantage of the FE method is that none of the above simplifying assumptions have to be made and the solution is much more rigorous. The earth models considered are isotropic, incompressible viscoelastic flat-earths with power-law rheology and no self-gravitation. The gravitational acceleration, g, and the density of the halfspace, ρ o, are taken to be constants. The flat-earth approximation has been demonstrated to be adequate in describing the rebound process for loads even as large as the Laurentian Ice sheet [19]. In this paper, it is assumed that the rebound process sees the steady state creep and not the transient creep. For a viscoelastic material, the deformation can be expressed as a sum of the elastic response and the creep response. The elastic strain is related to the stress by Hooke s law while the steady state creep law is given by: C ε ij =A * n-1 σ E σ ij (1) where A* is the creep parameter, σ' ij are the deviatoric stress components and σ' E is the equivalent deviatoric stress with σ' E = 1 σ' ij σ' ij (2) 2 For linear rheology, n equals 1. For nonlinear rheology, n normally lies between 2 and 6 for mantle rocks. For later discussion and comparison of results with other papers, it is useful to define the Table 1. Elastic parameters of halfspace earth models to be used with a load with size comparable to the Fennoscandian ice sheet. Density (kg/m 3 ) 338 Young's Moulus (Pa).35x Poisson's ratio. 5
3 effective viscosity as : η eff = 1 3A * (n - 1) σ E (3) Note that this differs from the usual definition [11, ] by a factor of 3/2, and the reason is discussed in Ranalli [1] (p.75-8). 3. Characteristic Deformation Pattern of Nonlinear Halfspaces The deformation pattern of a nonlinear mantle is quite distinct from that for a linear one. In this section, we shall review these differences in the deformation pattern as the stress exponent n increases from 1 to 4. We shall also investigate how the creep parameter A* affect this characteristic deformation pattern. The implication of the characteristic deformation pattern for a power-law medium on relative sea level (RSL) data is also discussed. For simplicity, the earth is considered to be a viscoelastic uniform halfspace with no ambient tectonic stress. The parameters of the earth model is listed in Table 1. The ice history is taken to be the Heaviside loading of a disc load with 9 km radius and a uniform magnitude of 2 MPa inside the load. Displ. (m) Vert t =. ka ka A = 3.33E n = 1 (a) t =. ka ka A=1.66E-28 n=2 (b) Displ. (m) Vert t =. ka ka A=8.33E-35 n=3 (c) t =. ka ka A=4.16E-41 n=4 (d) Radius (km) Radius (km) Figure 1. Vertical displacements at different times after the Heaviside loading of a uniform disc load on viscoelastic halfspaces with different stress exponent n. Elastic parameters of the models are given in Table 1.
4 The effects of stress exponent n on the vertical displacements is shown in Fig. 1. Here, the creep parameters, A*, are chosen so that the effective viscosity is 1 21 Pa-s. Comparing the deformation for the linear halfspace (Fig. 1a) with the nonlinear halfspaces, we see that power-law rheological models (n>1) are characterized by: i) Higher initial rate of deformation followed by much slower rate under the load [8]. Wu [16, 17, 18] have shown that this is due to the stress dependence of the effective viscosity for power-law medium (see equation 3). Thus, the initial flow rate is fast when the rebound stress is large, but slows down as the rebound stress relaxes. ii) Sharper cornered trough near the edge of the load. This is due to the high stress concentration near the edge of the load which results in low effective viscosity and fast relaxation of the shorter wavelengths of deformation. iii) Narrower peripheral trough and a larger amplitude, sharper cornered, peripheral bulge. This behavior is very similar to that for a linear channel and again can be explained by the ridge of high stress that extends beneath the load and peaks near the edge of the load (see Fig. 3 in Wu [17]). iv) Existence of a narrow viscously stationary zone (VSZ) or point. Within this narrow zone or point, the viscous displacement (i.e. displacement other than the elastic one at t=), almost remain stationary with time. From the centre of the load to this VSZ, land sinks continuously after the application of the load. Outside this VSZ, land rises continuously (see also Figs. 2b, c & d). On the other hand, a linear halfspace is characterized by a broad RSL transition zone outside the edge of the load, where land submergence is followed by emergence (see also Fig. 2a). The existence of a VSZ Displ. (m) Vert t= (a) A=3E n= ka t=.2.5 (b) A=1.66E n= ka Displ. (m) Vert t=.5 2. VSZ (c) ka Radius (km) A=3.33E-35 n= VSZ (d) A=1.66E-35 n= ka Radius (km) Figure 2. Vertical displacements at different times after the Heaviside loading of a uniform disc load on viscoelastic halfspaces with different combinations of A* and n. Elastic parameters of the models are given in Table 1.
5 VSZ and the lack of a broad RSL transition zone for power-law medium implies that nonlinear halfspaces cannot explain the sea-level data in the transition zone outside the former ice margin. As we shall see later, this existence of a VSZ or the lack of a broad RSL transition zone is characteristic of power-law medium even in the presence of ambient tectonic stress. Due to the importance of this VSZ on the interpretation of RSL data, we shall now investigate whether the creep parameter A* can affect the existence of the VSZ. In Fig. 2, the deformation near the edge of the load is plotted for nonlinear halfspaces (Figs. 2b-d) with n=3 and values of A* ranging from 1.66x1-35 to 1.66x1. For comparison, the deformation pattern for the linear halfspace is also shown (Fig. 2a). Fig. 2d shows that for small values of A*, the VSZ is located around 95 km. With the value of A* increased to 1.66x1, a VSZ is also found around 95 km for t < 2 ka, but as time progresses, this zone no longer remains stationary. Since the main difference between A*=1.66x1-35 and A*=1.66x1 is that the latter has a smaller effective viscosity than the former, this means that the VSZ in the former model is only stationary for the time under consideration. This zone is actually quasi-stationary - given enough time, VSZ does not actually exist. In order to see how the quasi stationary zone affect relative sea levels, vertical displacements relative to that at t= ka have been computed for the Heaviside unloading of the same uniform disc load, and plotted in Fig. 3 for four sites. Fig. 3a confirms that, at the centre of rebound (r= km), the nonlinear halfspaces give faster initial rate of uplift than the linear halfspace, but the present rates of uplift for the nonlinear halfspaces are much lower. The model with A*=1.66x1 has the lowest rate of uplift at the present because it has the smallest effective viscosity. Outside the load at r=92 km, the linear halfspace and the nonlinear model with A*= 1.66x1-35 show continuos land emergence(also shown in Fig. 2d), while the other nonlinear models show land emergence followed by submergence. Displ. (m) (a) R = km A=1.66E-35 n=3 A=8.33E-35 n=3 A=1.66E n=3 A=3.33E n=1 (b) R = 92 km Rel Displ. (m) (c) R = 948 km (d) R=165 km 25 Rel Time (ka BP) Time (ka BP) -5 Figure 3. Relative vertical displacement curves at different radius from the centre of a Fennoscandian size ice load.
6 However, the transition from land emergence to submergence takes place soon after the Heaviside unloading event around ka BP. At the VSZ around r=948 km, the model with A*= 1.66x1-35 show small land emergence followed by zero relative displacement after about 1 ka BP. Again, the other two nonlinear models show land emergence followed by submergence with very early time of transition. On the other hand, the linear halfspace predict a transition from land emergence to submergence around 6 ka BP. For r > 95 km, land submergence is predicted by all the nonlinear models except immediately after deglaciation at ka BP. However, the amount of submergence are very small in the last 8 ka. In contrast, the linear halfspace continue to predict land emergence followed by submergence around 8 ka BP, and the amount of submergence is about twice as much as that predicted by the nonlinear models. In summary, linear halfspace predicts a broad RSL transition zone, whereas nonlinear halfspace predicts a viscously quasi-stationary zone which result in a very narrow or nonexistent RSL transition zone. Even for those nonlinear models with a narrow RSL transition zone, the time of transition from emergence to submergence is much earlier than that predicted by the linear model and the amplitude of submergence is much smaller than that of the linear halfspace. Thus, RSL data near the former ice margin can be used to distinguish between linear and nonlinear halfspaces. 4. More complex earth models with parabolic disc load In this section, we are interested in finding earth models with nonlinear rheology and zero ambient tectonic stress that have the potential in explaining relative sea level data in Hudson Bay and along the east coast of Canada and the United States. For a load with size comparable to the Laurentia ice sheet, the deformation is sensitive to earth structure much deeper than the 67 km transition. Thus, the earth models considered here are stratified viscoelastic earths with a 22 km thick lithosphere (and no ambient tectonic stress). The elastic parameters of the earth models are listed in Table 2. Here, the Laurentian ice sheet is taken to be axially symmetric and the ice thickness at radius r is given by H ice (r) = H max α r where α =19 km is the radius of the disc load and H max is the maximum ice height at glacial maximum. Again, the load is assumed to be at isostatic equilibrium before being removed instantaneously from the surface of the earth at time T D. Note that H max and T D are not well determined, so the values for a linear uniform 1 22 Pa-s mantle ( H max =35 meter and T D = ka BP) will be used initially. However, in order to fit the sea level data with the nonlinear earth models, the timing of deglaciation T D is allowed to shift by 2 ka and H max can also be modified. Table 2. Elastic Parameters of stratified earth models Depth Density (kg/ m 3 ) Young's Modulus (Pa) Poisson Ratio Lithosphere Layer Layer Lower 1 2 Mantle - 22 km 34.87x km 34.87x km x1 below 67 km x
7 Model Name U22 U34 U35 U34B Lithospheric thickness (km) n in the mantle A* in the mantle Table 3. Uniform rheological earth models and ice models (Pa s - 1 ) 3.33x1 3.33x x1 3.33x1 Maximum ice thickness (m) Time of deglaciation (ka BP) 1 4.a Linear & Nonlinear Uniform Mantle Let us start with models that consists of a 22 km thick elastic lithospheres overlying stratified mantle with uniform rheology. The parameters for earth rheology and ice history are listed in Table 3. The letter U in the name of the models denote that rheology in the mantle is uniform, the numbers that follow are the negative of the exponent in the creep parameter. Fig. 4 shows the comparison of the observed sea level data and the predicted vertical displacements relative to the present. As shown in Wu & Johnson [19], relative vertical displacement curves give good approximations to relative sea level curves in the last 7 ka. Inspection of Fig. 4 shows that Model U22 (with linear rheology) is able to fit the sea level data simultaneously at the four sites - although a slightly thicker ice sheet will give a better fit to the data in Ottawa Island. On the other hand, all the nonlinear models predict land emergence instead of the observed submergence for sites outside the former ice margin (e.g. Boston Rel. Displ. (m) Ottawa Is Model U22 Model U35 Model U34 Model U34B NW NewFoundland Rel. Displ. (m) Boston Delaware Age (ka BP) Age (ka BP) -2 Figure 4. Comparing the predicted relative displacement curves with the observed sea level data (error bars) at 4 sites in Canada and along the east coast of the U.S. Axisymmetric (2D) parabolic ice models are used. Model parameters are given in Table 3.
8 and Delaware). This first order effect is due to the presence of VSZ in nonlinear mantle models and no modification of ice thickness nor time of deglaciation can change the predicted land emergence to submergence. Thus, nonlinear mantles are rejected by relative sea level data outside the former ice margin. Fig. 4 also illustrates the problems nonlinear mantles have in fitting the sea level data inside the ice margin. For H max =35 meter and T D = ka BP, Model U34 (short dashed line) is able to fit the RSL data in NW Newfoundland but under-predicts the height of the ancient beaches and the present-day rate of uplift in Ottawa Island. The RSL data in Ottawa Is. can be better explained by Model U34 if one delays the deglaciation time T D by 2 ka (Model U34B). However, Model U34B ( long dashed lines) now over predicts the height of the ancient beaches in NW Newfoundland, and for sites outside the former ice margin, land emergence is still predicted. Similarly, Model U35 can either explain the data in Ottawa Is. or in NW Newfoundland, but not both simultaneously. Thus, unlike the linear case, a simple ice profile cannot be used by nonlinear models to fit the RSL data in Ottawa Is. and Newfoundland simultaneously. 4.b Nonlinear Zone in the Upper Mantle with different thickness From the above, we saw that nonlinear rheology in the whole mantle is rejected by relative sea level data outside the former ice margin. In this subsection, nonlinear rheology is confined to a zone beneath the elastic lithosphere. This scenario may occur in the upper mantle because the viscosity for dislocation creep is generally smaller than the viscosity for diffusion near the base of the lithosphere, but the viscosity for diffusion becomes lower than that for dislocation for deeper depth [6]. Thus, the dominant creep mechanism may be dislocation near the base of the lithosphere, but may change to diffusion at deeper depth. The creep parameter in the NLZ is kept fixed at A*=3.33x1 (Pa -3 s -1 ), but the thickness of this nonlinear zone (NLZ) is varied (see Table 4). The purpose here is to find the thickest NLZ that can explain the RSL data in and around the centre of rebound. In the names of these models, the numbers following the label NLZ denote the thickness of the nonlinear zone in km. The predictions of these models are plotted in Fig. 5. For reference, the model U34 is also plotted. Table 4. Models with Nonlinear Zones of various thickness Model Name NLZ45 NLZ2 NLZ1 Lithospheric thickness (km) Thickness of NLZ in upper mantle (km) n in NLZ A* of NLZ (Pa s - 1 ) 3.33x1 3.33x1 3.33x1 n below NLZ A* in the mantle below the NLZ (Pa s - 1 ) 3.33x1 3.33x1 3.33x1 Maximum ice thickness (m) Time of deglaciation (ka BP)
9 Displ. (m) Rel Ottawa Is Model U34 Model NLZ45 Model NLZ2 Model NLZ1 NW NewFoundland Displ. (m) 2 1 Boston Delaware 2 1 Rel Age (ka BP) Age (ka BP) -2 Figure 5. Same as Fig. 4 except for the models in Table 4. Fig. 5 shows that with linear rheology in the lower mantle, land submergence is now predicted outside the former ice margin in Boston and Delaware. However, a nonlinear upper mantle above 67 km depth (Model NLZ45) under predicts the amount of land emergence in Ottawa Is. and also under predicts the amount of submergence outside the former ice margin. With the NLZ thickness reduced to 2 km or less, the RSL data both inside and immediately outside the former ice margin can be better explained, although a slightly thicker ice sheet would give even better fit to the data in Ottawa Is. Table 5. Models with Nonlinear Zones and various lower mantle viscosities Name NLZ45 Model NLZ45A NLZ45B Lithospheric thickness (km) Thickness of NLZ in upper mantle (km) n in NLZ A* of NLZ (Pa s - 1 ) 3.33x1 3.33x1 3.33x1 n in the lower mantle A* in the lower mantle (Pa s - 1 ) 3.33x1 1.11x x1 Maximum ice thickness (m) Time of deglaciation (ka BP)
10 Displ. (m) Rel Ottawa Is Model NLZ45 Model NLZ45A Model NLZ45B NW NewFoundland Displ. (m) Boston Delaware Rel Age (ka BP) Age (ka BP) Figure 6. Same as Fig. 4 except for the models in Table c Nonlinear Upper Mantle and Linear High Viscosity Lower Mantle For the previous models with nonlinear zones, the rheology of the lower mantle is linear with viscosity kept at 1x1 21 Pa-s. In this set of computation, we investigate the effect of increasing lower mantle viscosity to 3x1 21 Pa-s and 1x1 22 Pa-s (see table 5). The results, which are plotted in Fig. 6, show that increased lower mantle viscosity give land emergence instead of the observed submergence outside the former ice sheet. This is because increased lower mantle viscosity promotes channel flow which causes the forebulge to migrate outwards and result in land emergence outside the ice margin. Within the ice margin, the high viscosity in the lower mantle result in higher present-day rate of uplift, and thus the height of the ancient beaches are over-predicted. Thus, nonlinear upper mantle with high viscosity lower mantle is not compatible with sea level observations outside the former ice margin. Table 6. Earth models with Nonlinear Lower Mantle Model Name NLM21 NLM22 NLM23 Lithospheric thickness (km) n in upper mantle A* of upper mantle (Pa s - 1 ) 3.33x x x1 n in the lower mantle below 67 km A* in the lower mantle (Pa s - 1 ) 3.33x1 3.33x1 3.33x1 Maximum ice thickness (m) Time of deglaciation (ka BP)
11 4.d Nonlinear Lower Mantle and Linear Upper Mantle with Various Viscosities For the sake of completeness, let us consider the less likely case where nonlinear rheology is restricted to lie in the lower mantle and the viscosity in the upper mantle is varied from 1 2 Pa-s to 1 22 Pa-s (see table 6). In the name of the models, NLM stands for nonlinear lower mantle and the numbers that follow are the exponent of the creep parameter. The predicted relative displacement curves are compared with the observational data in Fig. 7. Inspection of Fig. 7 shows that a low viscosity (1 2 Pa-s) upper mantle on top of a nonlinear lower mantle (Model NLM21) cannot explain the RSL data outside the edge of the former ice margin because land emergence instead of the observed submergence is predicted. When the viscosity of the upper mantle is increased to 1 21 Pa-s (NLM22), a small amount of submergence is predicted outside the former ice margin, but the present-day rate of uplift outside the former ice margin is close to zero and the transition between land submergence and emergence is predicted to be as early as 11 ka BP. Furthermore, the height of the ancient beaches in Ottawa Is. are under-predicted by Model NLM22. When the viscosity of the upper mantle is further increased to 1 22 Pa-s (NLM23), the height of the 4-7 ka BP beaches in Ottawa Is. and the 8-1 ka BP beaches in NW Newfoundland can be explained, however, the amount of submergence outside the former ice margin is under predicted. To explain the amount of submergence outside the former ice margin, the thickness of the ice sheet has to be drastically increased, however, that would significantly over predict the height of the ancient beaches in Ottawa Is. and Newfoundland. Thus, a nonlinear lower mantle cannot satisfactorily explain all the sea level data in and around the centre of rebound simultaneously. In summary, of all the earth models investigated, the only type of earth model with nonlinear rheology and zero ambient tectonic stress that is consistent with sea level data inside and outside the former ice margin of the Laurentia is one that has a thin ( 2 km ) nonlinear zone. Rel. Displ. (m) a) Ottawa Is Model NLM21 Model NLM22 Model NLM23 b) NW NewFoundland Rel. Displ. (m) 4 c) Boston 2 d) Delaware Age (ka BP) Age (ka BP) Figure 7. Same as Fig. 4 except for the models in Table 6.
12 5. Nonlinear Mantle Rheology with Realistic ICE3G Model In the last two sections, only simple disc load with Heaviside loading/unloading history has been considered. In reality, the ice profile at glacial maximum was more complicated than a parabolic profile. In addition, the rate of deglaciation was not instantaneous but varied with time and location. Furthermore, the important effect of migration of the ice margin is not captured in those simple ice models. In this section, the realistic ICE3G model of Tushingham and Peltier [13, 14] with several saw-tooth glacial cycles that have slow buildup time of 9 ka but rapid deglacial time of 1 ka are used. The details of the ice model and the elastic parameters of the earth model can be found in Wu & Johnston [19]. In this section, only four uniform rheological models are considered - L22, L34, L35 & L36 (see Table 7). The letter L in their names denote that they all have a 15 km thick elastic lithosphere. The rheology of the underlying mantle is uniform and the numbers in their name represent the negative of the exponent in the creep parameter of the mantle. Except for Model L22, all the other models are nonlinear. The ambient tectonic stress is assumed to be zero here. In Fig. 8, the relative displacement curves for these models are compared with the observed RSL data in six sites in North America. Taking into account the problem of the earth s curvature for far-away sites (i.e. Brigantine and Southport) [19], Fig. 8 still show that the linear earth model has the best chance to explain all the RSL data simultaneously. However, in Churchill (Manitoba), neither the linear model nor the nonlinear models can predict the observed height of the ancient beaches - apparently, a thicker local ice or a delayed deglaciation is required by any of the models to give a better fit to the observations. In Ottawa Is., Model L34 significantly under predicts the height of the ancient beaches whereas in NW Newfoundland, Model L36 under predicts their height. In Boston, Model L36 predicts emergence rather than the observed submergence; Model L35 predicts zero rate of uplift but Model L34 is able to predict a small amount of submergence. In Brigantine (New Jersey), all three nonlinear models are able to predict submergence, but again the amount of submergence are too small. In Southport (North Carolina), all models, except Model L34, have no problem explaining the height of the submerged beaches. Comparing Figs. 4 & 8 shows that the effect of a realistic ice sheet is mainly in sites close to the ice margin at the last glacial maxima: With the inclusion of shrinking ice margin, land submergence are now predicted by nonlinear mantles with ICE3G. However, the amount of submergence is still not large enough to explain the observed RSL data. In summary, 3D models with realistic ice history show that nonlinear mantles are still rejected by RSL data just outside the ice margin at the last glacial maxima. However, more studies with realistic ice histories are needed to see if the conclusions for the other earth models in section 3 remain valid. Table 7. Uniform rheological earth models used with ICE3G Model Name L22 L34 L35 L36 Lithospheric thickness (km) n in the mantle A* of the mantle (Pa s - 1 ) 3.33x1 3.33x x x1
13 RSL (m) Ottawa Island Model L22 Model L36 Model L35 Model L34 Churchill RSL (m) NW Newfoundland Boston RSL (m) Brigantine, NJ Southport, NC Time (ka BP) Time (ka BP) -2 Figure 8. Same as Fig. 4 except for the earth models in Table 7 and realistic ice model ICE3G. 6. Effect of Ambient Tectonic Stress So far, all the rebound calculations in this paper have assumed that the ambient tectonic stress is zero. In the presence of ambient tectonic stress, σ T ij, the state of stress can be written as σ ij = σ ij T +σ ij R (4) where R σ ij is the deviatoric rebound stress. Substituting equation (4) into equations (1) to (2), and following the procedure outlined in Schmeling [11], one obtains: σ E = σ TE Σ i,j σ ij T σ ij R σ (5) where σ TE = 1 2 Σσ T T ij ij σ ij is the equivalent deviatoric tectonic stress. Also the creep law can
14 be written as: C T ε ij = εij +A n-1 σ TE σ R ij +A* n-1 σ 2 n-3 TE σ ijt Σ σ T R i,j ij σ ij (6) where the first term on the right side is the strain rate due to tectonic stress alone and is defined as: T ε ij =A * n-1 T σ TE σ ij The strain rate seen by the rebound process is thus: R ε n-1 ij =A*σ TE σ R ij +A* n-1 2 σ n-3 TE σ ijt Σσ T R ij σ ij i,j (8) Furthermore, equation (3) shows that the effective viscosity experienced by the rebound process in the presence of tectonic stress is: (7) η eff = 1+ n ηt eff Σ i,j σ ij T σ ij R σ TE 2 (9) where, η T eff = 1 n-1 3A*σ TE (1) is the effective viscosity seen by the tectonic process in the absence of rebound stress. Thus, equation (9) shows that the viscosity seen by postglacial rebound is different from that seen by mantle convection [] when the rheology is nonlinear. Some insight on how the interaction between rebound stress and tectonic stress affect creep in the mantle can be gained if we consider the following special case: Assume that the tectonic stress and the rebound stress components are orthogonal except for the kl-th component, then Σσ T i,j ij σ R ij = σ T R kl σ kl (11) Thus, from equation (8), the kl-th component of the strain rate is: R ε kl = 1 A* σ n-1 2 TE T 2 n-1 σ 2 R kl σ TE +2 σ kl () and the rebound stress component R σ kl sees a linear creep law but the viscosity is dependent on the tectonic stress distribution [11, 15]. If one further assumes that only the kl-th component of the tectonic stress contribute to the equivalent deviatoric tectonic stress, i.e. then equation () can be written as, σ TE = 1 2 σ kl T T σ kl (13) R ε n-1 kl =A*σ TE n-1 σ 4 TE +1 σr kl (14) The creep law for the other components (where rebound stress and tectonic stress are orthogonal) is: ε ij R =A*σ TE n-1 n-1 σ T ij T σ σ kl R R + σ ij kl (15)
15 If one further assumes that the ij-components of tectonic stress is zero, then: ε ij R =A*σ TE n-1 σ ij R (16) Again we have an effectively linear creep law for the ij-th component. Comparing equations (14) and (15) or (16), we see that the creep law for the ij-th component is different from that for the kl-th component - this means that the rebound process sees a linear but anisotropic creep law [11]. If the conclusions of this special case is applicable in general, then the rebound process will basically see a linear rheology even though the creep law is actually non-linear. This seems to imply that the RSL data inside and outside the former ice margin can be explained simultaneously by nonlinear mantle provided the right amount of ambient stress is applied. But are the conclusions of this special case applicable in general? Gasperini et al. [3] used the finite element method to study the vertical displacement history of a nonlinear upper mantle when the ambient tectonic stress is included in the equivalent stress (σ' E ) term only. With this simplified nonlinear creep law, they were able to show that the uplift history of linear models at the centre of the load can be approximated by nonlinear upper mantle models. So this does lend some support to the conclusions of the special case. Unfortunately, the uplift history in the RSL transition zone was not investigated by Gasperini et al [3]. Instantaneous loading/removal constant g of ice load o ICE Elastic Lithosphere ViscoElastic Power-law Mantle Fluid core Figure 9. Schematic diagram of the model that includes the interaction between ambient tectonic stress and rebound stress. The length of the arrows indicate the magnitude of the tectonic shear stress. After Wu [18]. Model Name A22 A36 A38 AUM36 Max. Ambient Stress (MPa) Lithospheric thickness (km) n in the upper mantle A* of upper mantle n in the lower mantle below 67km depth Table 8. Earth models with Ambient Tectonic Stress (Pa s - 1 ) 3.33x1 A* in the lower mantle below 67 km (Pa s x1 ) x x x x x1 3.33x1 Maximum ice thickness (m) Time of deglaciation (ka BP)
16 Vertical Displacement (m) A) Max. Shear =.1 MPa A=4E-35, n=3 t= B) Max. Shear = 1 MPa A=4E-35, n=3 t= kyr 16 kyr C) Max. Shear = 1 MPa A=4E-35, n=3 Ambient Stress (lines) vs No Ambient Stress (symbols) Ambient Stress (lines) vs No Ambient Stress (symbols) t = ka t = 1 ka t = 4 ka t =16 ka Distance from the centre of the load (km) Figure 1. Effect of ambient stress magnitude (.1, 1 and 1 MPa) on the deformation pattern at four time steps. Deformation calculated with the inclusion of ambient stress (lines) are compared with that calculated without ambient stress (symbols). After Wu [18]. Wu [18] also investigated this problem with the finite element method. In his study, the full interaction between tectonic stress and rebound stress in a power-law (n=3) medium was taken into account and no approximation was made in the creep law. However, to keep the problem simple, the plane strain assumption was employed, i.e. the Fennoscandian sized ice load was assumed to be a boxcar load with infinite extent in the y-direction perpendicular to the plane (see Fig. 9). The load was applied/removed instantaneously above the centre of a large convection cell, so that the tectonic stress seen were effectively horizontal shear stresses that varies quasi-linearly with depth and reached maxima at the top and bottom of the cell. The magnitude of the maximum shear was not fixed but was varied from.1 to 1 MPa in order to study the effect of ambient stress magnitude. His results showed that: (1) Ambient stress can be neglected in the calculation if its magnitude is much less than 1 MPa (Fig. 1a); (2) For ambient stress level of the order of 1 MPa, the effect of ambient stress is still small within the ice margin, but becomes important outside the former ice margin (Fig. 1b); (3) When the ambient stress level is much larger than 1 MPa, the large ambient stress results in a smaller but more time-independent effective viscosity (equation 3), thus the rate of defor-
17 Displ. (m) Rel A) Centre of Load Model A22 Model A36 Model A38 Model AUM Time (ka BP) B) Transition zone Time (ka BP) Figure 11. Effect of ambient stress on relative vertical displacement curves for the models in Table 8. Modified after Wu [18]. mation is higher inside the ice margin (Fig. 1c) but there is no rapid change in the rate of uplift as seen in nonlinear halfspaces (Fig. 3a). In this way, the deformation inside the ice margin appears to be linear. This is illustrated in Fig. 11a, where the relative displacement curves for the four earth models of Table 8 are shown. The letter A in the name of the models indicates that ambient tectonic stress is superimposed. The letters UM indicate that the upper mantle is nonlinear but the lower mantle is linear. The numbers that follow are the exponent of the creep parameter. Fig. 11a shows that at the centre of the ice load, models with nonlinear mantle do deform like the linear model, that is without the strong time dependent rate of deformation. The model with nonlinear rheology restricted to the upper mantle gives even better match to the linear model. On the other hand, Fig. 11b shows that RSL data outside the former ice margin can still be used to distinguish between a linear mantle from nonlinear mantle or nonlinear upper mantle because the nonlinear models produce too small rates of uplift at the present and the amount of submergence for the past 8 ka is under-predicted. The reason is again due to the rebound stress concentration near the edge of the ice sheet. Thus, the results of simple models show that large ambient tectonic stress (>> 1 MPa) can affect the rebound process in a power-law medium. In such case, the creep law may appear linear for the postglacial rebound within the ice margin, confirming the predictions of Weertman [15] and Schmeling [11]. On the other hand, for sites outside the ice margin, the RSL signature of a power-law mantle is still distinguishable from that due to deep flow in a linear mantle. This however needs to be verified with more realistic rebound modeling. 7. Conclusions Assuming that the rebound process sees the steady state creep and not the transient creep, it is shown that : 1) If the ambient tectonic stress level is low (<<1 MPa), nonlinear mantle rheology are characterized by a localized flow with strong time dependent uplift rate inside the ice margin and a viscously quasi-stationary zone outside the ice margin - thus the relative sea level observations outside the Laurentia ice margin cannot be explained by power-law mantle. 2) In the presence of large ambient tectonic stress (>> 1 MPa), land uplift in a power law mantle loses its strong time dependency for uplift rate inside the ice margin but maintains its characteristics outside the ice margin. 3) If the ambient stress level is low, a nonlinear lower mantle or a high viscosity lower mantle beneath
18 a nonlinear upper mantle are rejected by the sea level data outside the Laurentia ice margin. (However one has to be cautioned that sea level data may not be able to resolve the rheology at depth greater than about 15 km depth [7], thus lower mantle here actually means the upper part of the lower mantle.) 4) If the ambient stress level is low, a thin nonlinear rheology zone, with thickness around 2 km or less, embedded between a 22 km thick elastic lithosphere and a linear mantle is permitted by the sea level data in North America. Thus, the transition in dominant creep mechanism may occur around 3-4 km depth. This is consistent with the observations in seismic anisotropy [6]. 5) Preliminary 3D models, with realistic deglaciation histories but without ambient stress, confirm the characteristics of power law mantles found with 2D simple ice models. However, more 3D modeling with realistic ice histories are needed to confirm these finding for layered nonlinear rheology in the presence of ambient stress. References [1] C. Brennen, J. Geophys. Res., 79 (1974), p [2] S.T. Crough, Geophys.J.R.astr.Soc. 5 (1977), p [3] P. Gasperini, D.A. Yuen & R. Sabadini, Geophys. Res. Lett. 19 (1992), p [4] C. Goetze & D.L. Kohlstedt, J. Geophys. Res. 78 (1974), p [5] S. Karato & P. Li, Science 255 (1992), p [6] S. Karato & P. Wu, Science 26 (1993), p [7] J.X. Mitrovica & W.R. Peltier, Geophys.J.Int. 2 (1995), p [8] Nakada, M., J.Phys.Earth 31 (1983), p [9] R.L. Post & D.T. Griggs, Science 181 (1973), p [1] G. Ranalli, Rheology of the Earth: Deformation and Flow processes in Geophysics and Geodynamics, Allen & Unwin Inc. (1987). [11] H. Schmeling, H., Earth plant.sci Lett. 84 (1987), p [] D.L. Turcotte & G.Schubert, Geodynamics: applications of continuum physics to geological problems, Wiley & Sons, New York (1982). [13] A.M. Tushingham & W.R. Peltier, J.Geophys.Res. 96 (1991), p.4,497-4,523. [14] A.M. Tushingham & W.R. Peltier, J.Geophys.Res. 97 (1992), p.3,285-4,523. [15] J. Weertman, Phil.Trans.R.Soc.Lond.A. 288 (1978), p [16] P. Wu, Geophys.J.Int. 18 (1992), p [17] P. Wu, Geophys.J.Int. 114 (1993), p [18] P. Wu, Geophys. Res. Lett. 22 (1995), p [19] P. Wu & P. Johnston, this volume. [2] T. Yokokura & M. Saito, J.Phys.Earth 26 (1978), p
High-Harmonic Geoid Signatures due to Glacial Isostatic Adjustment, Subduction and Seismic Deformation
High-Harmonic Geoid Signatures due to Glacial Isostatic Adjustment, Subduction and Seismic Deformation L.L.A. Vermeersen (1), H. Schotman (1), M.-W. Jansen (1), R. Riva (1) and R. Sabadini (2) (1) DEOS,
More informationRheology of the Mantle and Plates (part 1): Deformation mechanisms and flow rules of mantle minerals
(part 1): Deformation mechanisms and flow rules of mantle minerals What is rheology? Rheology is the physical property that characterizes deformation behavior of a material (solid, fluid, etc) solid mechanics
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 informationHow important are elastic deflections in the Fennoscandian postglacial uplift?
How important are elastic deflections in the Fennoscandian postglacial uplift? WILLY FJELDSKAAR Fjeldskaar, W. How important are elastic deftections in the Fennoscandian postglacial uplift? Norsk Geologisk
More information7. Post Glacial Rebound. Ge 163 4/16/14-
7. Post Glacial Rebound Ge 163 4/16/14- Outline Overview Order of magnitude estimate of mantle viscosity Essentials of fluid mechanics Viscosity Stokes Flow Biharmonic equation Half-space model Channel
More informationSubduction II Fundamentals of Mantle Dynamics
Subduction II Fundamentals of Mantle Dynamics Thorsten W Becker University of Southern California Short course at Universita di Roma TRE April 18 20, 2011 Rheology Elasticity vs. viscous deformation η
More informationSupplementary information on the West African margin
Huismans and Beaumont 1 Data repository Supplementary information on the West African margin Interpreted seismic cross-sections of the north Angolan to south Gabon west African passive margins 1-3, including
More informationConstraints on lithosphere and mantle rheology from in-situ observations
Constraints on lithosphere and mantle rheology from in-situ observations Shijie Zhong Department of Physics University of Colorado at Boulder Collaborators: Archie Paulson and John Wahr on post-glacial
More informationDETAILS ABOUT THE TECHNIQUE. We use a global mantle convection model (Bunge et al., 1997) in conjunction with a
DETAILS ABOUT THE TECHNIQUE We use a global mantle convection model (Bunge et al., 1997) in conjunction with a global model of the lithosphere (Kong and Bird, 1995) to compute plate motions consistent
More informationSeismic and flexure constraints on lithospheric rheology and their dynamic implications
Seismic and flexure constraints on lithospheric rheology and their dynamic implications Shijie Zhong Dept. of Physics, University of Colorado Boulder, Colorado, USA Acknowledgement: A. B. Watts Dept. of
More informationConstraints on Mantle Structure from Surface Observables
MYRES I: Heat, Helium & Whole Mantle Convection Constraints on Mantle Structure from Surface Observables Magali Billen University of California, Davis Department of Geology The Goal Use observations of
More informationThree-dimensional finite-element modelling of Earth s viscoelastic deformation: effects of lateral variations in lithospheric thickness
Geophys. J. Int. (23) 1, 679 695 Three-dimensional finite-element modelling of Earth s viscoelastic deformation: effects of lateral variations in lithospheric thickness Shijie Zhong, Archie Paulson and
More informationThe importance of the South-American plate motion and the Nazca Ridge subduction on flat subduction below South Peru
Chapter 7 The importance of the South-American plate motion and the Nazca Ridge subduction on flat subduction below South Peru Abstract Flat subduction near Peru occurs only where the thickened crust of
More informationSurface changes caused by erosion and sedimentation were treated by solving: (2)
GSA DATA REPOSITORY 214279 GUY SIMPSON Model with dynamic faulting and surface processes The model used for the simulations reported in Figures 1-3 of the main text is based on two dimensional (plane strain)
More informationCourse Business. Today: isostasy and Earth rheology, paper discussion
Course Business Today: isostasy and Earth rheology, paper discussion Next week: sea level and glacial isostatic adjustment Email did you get my email today? Class notes, website Your presentations: November
More informationGravitational constraints
Gravitational constraints Reading: Fowler p172 187 Gravity anomalies Free-air anomaly: g F = g g( λ ) + δg obs F Corrected for expected variations due to the spheroid elevation above the spheroid Bouguer
More informationIntroduction to Geology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 12.001 Introduction to Geology Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. RHEOLOGICAL MODELS Rheology
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 informationLecture 5. Rheology. Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm
Lecture 5 Rheology Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm WW Norton; unless noted otherwise Rheology is... the study of deformation and flow of
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 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 informationSeismotectonics of intraplate oceanic regions. Thermal model Strength envelopes Plate forces Seismicity distributions
Seismotectonics of intraplate oceanic regions Thermal model Strength envelopes Plate forces Seismicity distributions Cooling of oceanic lithosphere also increases rock strength and seismic velocity. Thus
More informationISSN Scientific Technical Report STR 07/05. DOI: /GFZ.b GeoForschungsZentrum Potsdam
ISSN 1610-0956 Ingo Sasgen, Robert Mulvaney, Volker Klemann and Detlef Wolf Glacial-isostatic adjustment and sea-level change near Berkner Island, Antarctica Ingo Sasgen Dep. 1: Geodesy and Remote Sensing
More informationCan we see evidence of post-glacial geoidal adjustment in the current slowing rate of rotation of the Earth?
Can we see evidence of post-glacial geoidal adjustment in the current slowing rate of rotation of the Earth? BARRETO L., FORTIN M.-A., IREDALE A. In this simple analysis, we compare the historical record
More information(b) What is the amplitude at the altitude of a satellite of 400 km?
Practice final quiz 2015 Geodynamics 2015 1 a) Complete the following table. parameter symbol units temperature T C or K thermal conductivity heat capacity density coefficient of thermal expansion volumetric)
More informationPreviously Documented Basin-localized Extension on Mercury
GSA DATA REPOSITORY 2012323 Watters et al. Previously Documented Basin-localized Extension on Mercury The most widely distributed extensional landforms on Mercury seen in images from the Mariner 10 and
More informationRheology. What is rheology? From the root work rheo- Current: flow. Greek: rhein, to flow (river) Like rheostat flow of current
Rheology What is rheology? From the root work rheo- Current: flow Greek: rhein, to flow (river) Like rheostat flow of current Rheology What physical properties control deformation? - Rock type - Temperature
More informationGravity Tectonics Volcanism Atmosphere Water Winds Chemistry. Planetary Surfaces
Gravity Tectonics Volcanism Atmosphere Water Winds Chemistry Planetary Surfaces Gravity & Rotation Polar flattening caused by rotation is the largest deviation from a sphere for a planet sized object (as
More informationDynamic analysis. 1. Force and stress
Dynamic analysis 1. Force and stress Dynamics is the part of structural geology that involves energy, force, stress, and strength. It's very important to distinguish dynamic concepts from kinematic ones.
More informationSignificance of the fundamental mantle rotational relaxation mode in polar wander simulations
Geophys. J. Int. (1996) 127, F5-F9 FAST-TRACK PAPER Significance of the fundamental mantle rotational relaxation mode in polar wander simulations L. L. A. Vermeersen* and R. Sabadini Dipartimento di Scienze
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 informationInference of mantle viscosity from GRACE and relative sea level data
Geophys. J. Int. (7) 7, 97 58 doi:./j.365-6x.7.3556.x FAST TRACK PAPER Inference of mantle viscosity from GRACE and relative sea level data Archie Paulson, Shijie Zhong and John Wahr Department of Physics,
More informationAnswers: Internal Processes and Structures (Isostasy)
Answers: Internal Processes and Structures (Isostasy) 1. Analyse the adjustment of the crust to changes in loads associated with volcanism, mountain building, erosion, and glaciation by using the concept
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION Supplementary information: Our 1080 km x 360 km model setup includes a 7 km thick oceanic crust adjacent to a 60 km thick, 250 km wide orogenic crust. Both the oceanic and the
More informationWhat is the LAB Dynamically: Lithosphere and Asthenosphere Rheology from Post-loading Deformation
What is the LAB Dynamically: Lithosphere and Asthenosphere Rheology from Post-loading Deformation Roland Bürgmann, UC Berkeley With contributions by Pascal Audet, Daula Chandrasekhar, Georg Dresen, Andy
More informationRheology III. Ideal materials Laboratory tests Power-law creep The strength of the lithosphere The role of micromechanical defects in power-law creep
Rheology III Ideal materials Laboratory tests Power-law creep The strength of the lithosphere The role of micromechanical defects in power-law creep Ideal materials fall into one of the following categories:
More informationDeformation cycles of great subduction earthquakes in a viscoelastic Earth
Deformation cycles of great subduction earthquakes in a viscoelastic Earth Kelin Wang Pacific Geoscience Centre, Geological Survey of Canada School of Earth and Ocean Science, University of Victoria????
More informationCrust : wet quartzite Arc root : dry olivine mantle = 2840 kg/m km = 3300 kg/m km (Arc root thickness) 280 km (Arc width)
Crust : wet quartzite Arc root : dry olivine mantle = 2840 kg/m 3 41.5 km = 3300 kg/m 3 o 118.5 km Temperature (C) free surface z = 0 550 1350 160 km (Arc root thickness) 280 km (Arc width) sub-lithospheric
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 informationRheology and the Lithosphere
Rheology and the Lithosphere Processes in Structural Geology & Tectonics Ben van der Pluijm WW Norton+Authors, unless noted otherwise 3/8/2017 16:51 We Discuss Rheology and the Lithosphere What is rheology?
More informationThe numerical method used for experiments is based on an explicit finite element
Bialas 1 Model Supplementary Data The numerical method used for experiments is based on an explicit finite element technique similar to the Fast Lagrangian Analysis of Continua (FLAC) method (Cundall,
More informationby west-dipping subduction of the North American Plate
GSA Data Repository 2019133 Supplementary information for Kimberlite magmatism induced by west-dipping subduction of the North American Plate Wenbo Zhang 1*, Stephen T. Johnston 1, and Claire A. Currie
More informationBenchmarks for subduction zone models Subduction zone workshop, University of Michigan, July 2003
Introduction Benchmarks for subduction zone models Subduction zone workshop, University of Michigan, July 2003 In early October 2002 a group of researchers met at the University of Michigan at Ann Arbor,
More informationDefining the former elevation and shape of the lithosphere, in particular the elevation of the Earth s surface,
Isostasy in Move Defining the former elevation and shape of the lithosphere, in particular the elevation of the Earth s surface, is important in the restoration of a model as it aids in reducing uncertainty
More informationThe influence of short wavelength variations in viscosity on subduction dynamics
1 Introduction Deformation within the earth, driven by mantle convection due primarily to cooling and subduction of oceanic lithosphere, is expressed at every length scale in various geophysical observations.
More informationExercise: concepts from chapter 8
Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic
More informationA mechanical model of the San Andreas fault and SAFOD Pilot Hole stress measurements
GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L15S13, doi:10.1029/2004gl019521, 2004 A mechanical model of the San Andreas fault and SAFOD Pilot Hole stress measurements Jean Chéry Laboratoire Dynamique de la
More informationROTATING RING. Volume of small element = Rdθbt if weight density of ring = ρ weight of small element = ρrbtdθ. Figure 1 Rotating ring
ROTATIONAL STRESSES INTRODUCTION High centrifugal forces are developed in machine components rotating at a high angular speed of the order of 100 to 500 revolutions per second (rps). High centrifugal force
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION DOI: 10.1038/NGEO1887 Diverse calving patterns linked to glacier geometry J. N. Bassis and S. Jacobs 1. Supplementary Figures (a) (b) (c) Supplementary Figure S1 Schematic of
More informationtopography half is commonly represented in maps, and how it interacts with
Topography T. Perron 12.001 We ll spend a large fraction of the second half of the course discussing Earth s surface. Today we ll do two things: First, we ll discuss the ways topography is commonly represented
More informationThe anisotropic and rheological structure of the oceanic upper mantle from a simple model of plate shear
Geophys. J. Int. (24) 158, 287 296 doi: 1.1111/j.1365-246X.24.225.x The anisotropic and rheological structure of the oceanic upper mantle from a simple model of plate shear Noah S. Podolefsky, Shijie Zhong
More informationRheology: What is it?
Schedule Rheology basics Viscous, elastic and plastic Creep processes Flow laws Yielding mechanisms Deformation maps Yield strength envelopes Constraints on the rheology from the laboratory, geology, geophysics
More informationThe Frictional Regime
The Frictional Regime Processes in Structural Geology & Tectonics Ben van der Pluijm WW Norton+Authors, unless noted otherwise 1/25/2016 10:08 AM We Discuss The Frictional Regime Processes of Brittle Deformation
More informationDevelopment of Anisotropic Structure by Solid-State Convection in the Earth s Lower Mantle
1 Development of Anisotropic Structure by Solid-State Convection in the Earth s Lower Mantle Allen K. McNamara *, Peter E. van Keken, * & Shun-Ichiro Karato ** * Department of Geological Sciences, University
More informationD DAVID PUBLISHING. Deformation of Mild Steel Plate with Linear Cracks due to Horizontal Compression. 1. Introduction
Journal of Control Science and Engineering 1 (2015) 40-47 doi: 10.17265/2328-2231/2015.01.005 D DAVID PUBLISHING Deformation of Mild Steel Plate with Linear Cracks due to Horizontal Compression Mitsuru
More informationBrittle Deformation. Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm
Lecture 6 Brittle Deformation Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm WW Norton, unless noted otherwise Brittle deformation EarthStructure (2 nd
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 informationOriginally published as:
Originally published as: Lorenzo Martín, F., Wang, R., Roth, F. (2002): The effect of input parameters on visco-elastic models of crustal deformation. - Física de la Tierra, 14, 33-54 The effect of input
More informationNeogene Uplift of The Barents Sea
Neogene Uplift of The Barents Sea W. Fjeldskaar A. Amantov Tectonor/UiS, Stavanger, Norway FORCE seminar April 4, 2013 The project (2010-2012) Funding companies Flat Objective The objective of the work
More informationLimitations on the inversion for mantle viscosity from postglacial rebound
Geophys. J. Int. (27) 168, 119 129 doi: 1.1111/j.136-246X.26.3222.x Limitations on the inversion for mantle viscosity from postglacial rebound Archie Paulson, Shijie Zhong and John Wahr Department of Physics,
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationPlate tectonics. Chapter Slab pull. Dynamics of the Mantle and Lithosphere
Chapter 6 In a number of aspects the Earth s mantle can be regarded as a fluid that shows simple fluid-dynamical behavior. or instance, the mantle of the Earth is probably showing a form of convection.
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 informationDynamics 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 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 informationGeophysics Departmental Exam: 2004 Part 1
2004 Geophysics Departmental Exam: 2004 Part 1 This section is 90 minutes, closed book, and consists of questions designed to test your knowledge of facts and figures in the geosciences. The focus will
More informationRheology, or the study of the flow of matter. Panta rei. (Panta rhei)
Rheology, or the study of the flow of matter Panta rei (Panta rhei) Overview Basics of rheology Linear elasticity Linear viscosity Linear viscoelasticity To infinity... and beyond! Coming back to Earth
More informationJournal of Applied Mathematics and Computation (JAMC), 2018, 2(7),
Journal of Applied Mathematics and Computation (JAMC), 2018, 2(7), 266-270 http://www.hillpublisher.org/journal/jamc ISSN Online:2576-0645 ISSN Print:2576-0653 The Solution Of 2D Hydrodynamic Equations
More informationEarthquake and Volcano Deformation
Earthquake and Volcano Deformation Paul Segall Stanford University Draft Copy September, 2005 Last Updated Sept, 2008 COPYRIGHT NOTICE: To be published by Princeton University Press and copyrighted, c
More informationAsish Karmakar 1, Sanjay Sen 2 1 (Corresponding author, Assistant Teacher, Udairampur Pallisree Sikshayatan (H.S.), Udairampur, P.O.
IOSR Journal of Applied Geology and Geophysics (IOSR-JAGG) e-issn: 3 99, p-issn: 3 98.Volume 4, Issue 5 Ver. III (Sep. - Oct. 6), PP 39-58 www.iosrjournals.org A Sudden Movement across an Inclined Surface
More information4. Basins due to orogenic loading
4. Basins due to orogenic loading Orogenic loading and lithospheric flexure: peripheral vs. retroarc foreland basin Elastic vs. viscoelastic models Gravity anomalies, foreland flexure, and the development
More informationLECTURE 5 - Wave Equation Hrvoje Tkalčić " 2 # & 2 #
LECTURE 5 - Wave Equation Hrvoje Tkalčić " 2 # "t = ( $ + 2µ ) & 2 # 2 % " 2 (& ' u r ) = µ "t 2 % & 2 (& ' u r ) *** N.B. The material presented in these lectures is from the principal textbooks, other
More informationTheory at a Glance (for IES, GATE, PSU)
1. Stress and Strain Theory at a Glance (for IES, GATE, PSU) 1.1 Stress () When a material is subjected to an external force, a resisting force is set up within the component. The internal resistance force
More informationPhysics and Chemistry of the Earth and Terrestrial Planets
MIT OpenCourseWare http://ocw.mit.edu 12.002 Physics and Chemistry of the Earth and Terrestrial Planets Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationModeling the Thermal-Mechanical Behavior of Mid-Ocean Ridge Transform Faults
Excerpt from the Proceedings of the COMSOL Conference 2008 Boston Modeling the Thermal-Mechanical Behavior of Mid-Ocean Ridge Transform Faults Emily C Roland *1, Mark Behn,2 and Greg Hirth 3 1 MIT/WHOI
More informationVertical Motion from Satellite Altimetry and Tide gauges
Vertical Motion from Satellite Altimetry and Tide gauges Alexander Braun, C.K. Shum and C.-Y. Kuo Byrd Polar Research Center and Laboratory for Space Geodesy and Remote Sensing, The Ohio State University,
More informationModification of the lithospheric stress field by lateral variations in plate-mantle coupling
GEOPHYSICAL RESEARCH LETTERS, VOL. 36, L22307, doi:10.1029/2009gl040484, 2009 Modification of the lithospheric stress field by lateral variations in plate-mantle coupling J. B. Naliboff, 1 C. P. Conrad,
More information7.2.1 Seismic waves. Waves in a mass- spring system
7..1 Seismic waves Waves in a mass- spring system Acoustic waves in a liquid or gas Seismic waves in a solid Surface waves Wavefronts, rays and geometrical attenuation Amplitude and energy Waves in a mass-
More information**********************************************************************
Department of Civil and Environmental Engineering School of Mining and Petroleum Engineering 3-33 Markin/CNRL Natural Resources Engineering Facility www.engineering.ualberta.ca/civil Tel: 780.492.4235
More informationSEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by
SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e
More informationOn the Purpose of Toroidal Motion in a Convecting Mantle
GEOPHYSICAL RESEARCH LETTERS, VOL. 22, PAGES 3107-3110, 1995 On the Purpose of Toroidal Motion in a Convecting Mantle David Bercovici Department of Geology Geophysics, SOEST, University of Hawaii, Honolulu
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION doi: 10.1038/ngeo739 Supplementary Information to variability and distributed deformation in the Marmara Sea fault system Tobias Hergert 1 and Oliver Heidbach 1,* 1 Geophysical
More informationReference material Reference books: Y.C. Fung, "Foundations of Solid Mechanics", Prentice Hall R. Hill, "The mathematical theory of plasticity",
Reference material Reference books: Y.C. Fung, "Foundations of Solid Mechanics", Prentice Hall R. Hill, "The mathematical theory of plasticity", Oxford University Press, Oxford. J. Lubliner, "Plasticity
More informationElements of Rock Mechanics
Elements of Rock Mechanics Stress and strain Creep Constitutive equation Hooke's law Empirical relations Effects of porosity and fluids Anelasticity and viscoelasticity Reading: Shearer, 3 Stress Consider
More informationRHEOLOGY & LINEAR ELASTICITY. B Importance of fluids and fractures in deformation C Linear elasticity for homogeneous isotropic materials
GG303 Lecture 2 0 9/4/01 1 RHEOLOGY & LINEAR ELASTICITY I II Main Topics A Rheology: Macroscopic deformation behavior B Importance of fluids and fractures in deformation C Linear elasticity for homogeneous
More informationForces That Shape Earth. How do continents move? What forces can change rocks? How does plate motion affect the rock cycle?
Forces That Shape Earth How do continents move? What forces can change rocks? How does plate motion affect the rock cycle? Plate Motion Mountain ranges are produced by plate tectonics. The theory of plate
More informationLateral, radial, and temporal variations in upper mantle viscosity and rheology under Scandinavia
Article Volume 12, Number 1 29 January 2011 Q01007, doi:10.1029/2010gc003290 ISSN: 1525 2027 Lateral, radial, and temporal variations in upper mantle viscosity and rheology under Scandinavia Auke Barnhoorn
More informationEvaluation of Flexural Stiffness for RC Beams During Fire Events
3 rd International Structural Specialty Conference 3 ième conférence internationale spécialisée sur le génie des structures Edmonton, Alberta June 6-9, 202 / 6 au 9 juin 202 Evaluation of Flexural Stiffness
More informationFrictional rheologies have a wide range of applications in engineering
A liquid-crystal model for friction C. H. A. Cheng, L. H. Kellogg, S. Shkoller, and D. L. Turcotte Departments of Mathematics and Geology, University of California, Davis, CA 95616 ; Contributed by D.
More informationWhy does the Nazca plate slow down since the Neogene? Supplemental Information
GSA DATA REPOSITORY 2013004 Why does the Nazca plate slow down since the Neogene? Supplemental Information Javier Quinteros a,b,, Stephan V. Sobolev a,c a Deutsches GeoForschungsZentrum GFZ, Telegrafenberg,
More informationFinite element modelling of fault stress triggering due to hydraulic fracturing
Finite element modelling of fault stress triggering due to hydraulic fracturing Arsalan, Sattari and David, Eaton University of Calgary, Geoscience Department Summary In this study we aim to model fault
More informationEffect of Anisotropic Permeability on Band Growth
Effect of Anisotropic Permeability on Band Growth Jesse Taylor-West October 7, 24 Introduction The mechanics of partially molten regions of the mantle are not well understood- in part due to the inaccessibility
More informationA short review of continuum mechanics
A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material
More informationAnomalous secular sea-level acceleration in the Baltic Sea caused by isostatic adjustment
ANNALS OF GEOPHYSICS, 57, 4, 2014, S0432; doi:10.4401/ag-6548 Anomalous secular sea-level acceleration in the Baltic Sea caused by isostatic adjustment Giorgio Spada 1,*, Marco Olivieri 2, Gaia Galassi
More information12 Gravity and Topography
12. Gravity and Topography 126 12 Gravity and Topography In the limit of perfect hydrostatic equilibrium, there is an exact relationship between gravity and topography... and no new information emerges:
More informationGLY Geomorphology Notes
GLY 5705 - Geomorphology Notes Dr. Peter N. Adams Revised: Sept. 2012 10 Flexure of the Lithosphere Associated Readings: Anderson and Anderson (2010), pp. 86-92 Here we ask the question: What is the response
More informationFrequently Asked Questions
Frequently Asked Questions Why do we have to make the assumption that plane sections plane? How about bars with non-axis symmetric cross section? The formulae derived look very similar to beam and axial
More informationThe finite difference code (fully staggered grid) includes a marker-in-cell Lagrangian marker
GSA Data Repository 2018289 Ruh et al., 2018, Shale-related minibasins atop a massive olistostrome in an active accretionary wedge setting: Two-dimensional numerical modeling applied to the Iranian Makran:
More informationElizabeth H. Hearn modified from W. Behr
Reconciling postseismic and interseismic surface deformation around strike-slip faults: Earthquake-cycle models with finite ruptures and viscous shear zones Elizabeth H. Hearn hearn.liz@gmail.com modified
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 information