Glacial Rebound and Sea-Level Change in the British Isles

Total Page:16

File Type:pdf, Size:1020Kb

Glacial Rebound and Sea-Level Change in the British Isles TEmARESEARCH Glacial rebound and sea-level change in the British Isles Kurt Lambeck Research School of Earth Sciences, The Australian National University, Canberra ACT 2601, Australia ABSTRACT dominates over the rise in sea-level from the meltwater of the distant ice Observations of sea levels around the coastline of the British Isles for sheets, whereas in early Holocene time the past 10,000-15,000 years exhibit a major regional variation and the contributions from the meltwater of provide an important data base for testing models of glacial rebound as the distant ice sheets begins to dom- well as models of the Late Devensian ice sheet. A high-resolution inate. In the late Holocene time, after about 6000 yr BP, it is the crustal rebound rebound model has been developed which is consistent with both the that again is the dominant cause of sea- spatial and temporal patterns of sea-level change and which level change. Further away from the demonstrates that the observations are the result of (i) the glacio- centre of rebound, in NW England (Fig. isostatic crustal rebound in response to the unloading of the ice sheet 2a), the balance between the two con- over Britain and, to a lesser degree, of the ice sheet over Fennoscandia, tributions is primarily in favour of the and (ii) the rise in sea-level from the melting Late Pleistocene ice second contribution except in Latest De- vensian time when the rebound compo- sheets, including the response of the crust to the water loading (the nent dominates. The observations from hydro-isostatic effect). The agreement between model and sites in southern England (Fig. 2b-d) observations is such that there is no need to invoke vertical crustal exhibit less spatial variability and the movements for Great Britain and Ireland of other than glacio-hydro- sea-level change here is controlled isostatic origin. The rebound contributions are important throughout mainly, but not wholly, by the melt- the region and nowhere is it sufficiently small for the sea-level change water from the major ice sheets. to approximate the eustatic sea-level rise. The observational data There has been little quantitative modelling of the rebound of Britain and distribution around the periphery as well as from sites near the centre the associated sea-level fluctuations, of the former ice sheet is sufficient to permit constraints to be despite a number of geophysical, glaci- established on both earth model parameters specifying the mantle ological and geomorphological reasons viscosity and lithospheric thickness and the extent and volume of the for doing so. The rebound models can ice sheet at the time of the last glaciation. Preliminary solutions are provide constraints on (i) the Earth's presented which indicate an upper mantle viscosity of (3-5)1eoPa s, a response to surface loading; (5) the lithospheric thickness of about 100 km or less, and an ice model that ice sheet dimensions and the timing was not confluent with the Scandinavian ice sheet during the last of the retreat of the ice sheets; and (iii) the ages of some of the erosional fea- glaciation and whose maximum thickness over Scotland is unlikely to tures, such as the rock platforms of have exceeded about 1500 m. western Scotland (e.g. Dawson, 1984) or the older elevated beaches of eastern Terra Nova, 3,379-389 and northeastern Scotland (e.g. Firth, 1989), where conventional dating methods have failed. such model- I NTRO D UCTl ON Also, Europe, North America and Antarctica ling is desirable in order to establish Complex temporal and spatial patterns (e.g. Donner, 1970; Cullingford et al., whether sea-level observations from of sea-level change have been recorded 1981). The Scottish results, from sites different sites can be combined to form throughout Great Britain and Ireland near the centre of rebound, exhibit this representative sea-level curves for a during Late Devensian and Holocene spatial variability even within restricted given region. This question, in particu- time and the observations illustrated in areas such as the Forth, Tay, and Moray lar, is addressed in this paper. Fig. 1 are characteristic of this regional Firths of eastern Scotland and this is The solution of the glacial rebound variability. This pattern has been qual- perhaps best illustrated by the shoreline and sea-level problem for the Great Brit- itatively understood in terms of the gradient diagrams of Sissons (1983a) ain region requires a knowledge of the crustal rebound in response to the re- and others. The characteristic temporal Earths rheology (E) and of this history treat of the Late Devensian ice sheet patterns of the Scottish observations are of the global ice sheets (I). This ice sheet over Britain and to the rise in sea-level the combined result of these two con- includes the small but important British resulting from the concomitant melting tributions: in Late Devensian time the ice sheet ZB, as well as the more volumin- of distant major ice sheetsover northern glaao-isostatic rebound contribution ous Fennoscandian ice sheet lF and 379 KURT LAMBECK 40, . I of the various parameters is possible Arnprlor Beaulg Flrlh here (Lambeck, 1991) and that, in con- Uppsr Farth Valley - m Northeastern Scotland . sequence, the observation of sea-level p.. change in the recent past can contribute :I I, to several areas of earth science. 0 FORMU LATlON -10 .*--_.,. .___.-- I The formulation of Nakada and Lam- -LO ' J beck (1987, 1989) is used and the sea- IS 10 J 0 a level equation is written as Al(cp,A:t) = A&(t)+Ap(cp, A:t) Lochpll head +Af(cp,A:t)+A{r-'(cp,A:t) Western &oUand . (1) A{(cp,A:t) is the position of mean sea- level at latitude Q, longitude A and at time t relative to the present position at the same location. The term A&(t) is the equivalent sea-level function (or eu- . ..._._.- -.-_. .- .lo static sea-level) defined as I .a A&(t) = (change in uolume of meltwater)/ IS 10 J 0 10 J 0 C (surface area ofocenns) It includes the time dependence of the Brldge d EPrn Ardyne totality of all ice sheets. The three re- Firth of Tag . Firth of Clyde . maining terms on the right-hand-side of (1) include the contributions, for each of the three ice models IB, IF and f-f separately, from (i) the loading-unload- ._._.--._ ing of the crust by the ice sheets, (ii) the loading-unloading of the crust by the 0 ...~-.- meltwater deposited into or removed '. .. ,- I -10 ,._, from the oceans, and (iii) the gravita- IS 10 I 0 e tional self-attraction of the water and ice that ensures that the ocean surface re- Fig. 1. Observed sea-level heights (with error bars where available, and reduced to mean sea-level) mains an equipotential surface at all and predicted lmls at a number of sites in Scotland. The predicted lmls are based on the Scotland times. These 'corrective terms' are func- rebound model (Table 1)and the dashed lines represent the range of values predicted by the upper and tions of position and time. Their evalua- lower limits for the muntleviscosity. (See Fig. 5 for location map.)(a) Arnprior, upper Forth Valley. tion requires a knowledge of the ice Thedataarefrom Sissonsand Brooks (1971);(b) BridgeofEarth, FirthofTay, nearPerth ineastern models I(cp,A:t), of the earth model E, Scotland. Theobservations arefrom Cullingford et al. (2980) and Paterson et al. (7980)for the two and of the mean geometry (which will oldest data points (c) Beauly Firth, northeastern Scotland. The data are from Firth and Haggart itself be a function of time). The earth (2989)for the Holocene timeandfrom Firth (1989)forLate Devensian time; fd)Lochgilphead, Loch model parameters include the effective Fyne, western Scotland. Data are from Peacock et a1.(1978);(e) Ardyne, Firth of Clyde, western thickness of the lithosphere He as well Scotland. Data arefrom Pencock et a1 (2978); ff)Northern shoreofthe Solway Firth, southwestern as the effective viscosity of the mantle q, Scotland, between Dumfries and Kirkcudbright. Data are from Iardine (1982). where 7 will be a function of earth radius. Furthermore, both Heand qmay exhibit lateral variation. The observation model is the far-field ice sheets 1-1 for Laurentia, whether, in view of the uncertainties as- Antarctica and the BarentsKara Seas sociated with both the € and I models region. In glaciological studies it has and with incomplete data sets for past generally been assumed that E is known positions of sea-level, it is indeed pos- and ice models are constructed using a sible to establish independent improve- formulation relating the forces acting ments on both the € and 1 model on and within the ice sheet and the ice parameters as well as resolve some of rheology. In most geophysical studies it the ambiguities in the observational has been assumed that I is known and data. Preliminary modelling of the re- where AS, = observed sea-level (re- earth models are inferred from the sea- bound of Great Britain has, how- duced to mean sea-level) at location level observations. It may be asked ever, indicated that some separation (cp,A) and at time t; E~ = Observation 380 GLACIAL REBOUND AND SEA-LEVEL CHANGE 10 , I (iii) the lithosphere, of thickness Hp, is 0 0 of high (165 Pa s) viscosity and overlies an upper mantle (extend- -10 . ing to a depth of 650 km) with a viscosity qvmand a lower mantle of -m . viscosity qfm. .XI . The earth model parameters are there- Moreearnbe Bay Avonmoutb Norlhwest England Upper Brlstol Channel fore effective parameters only that de- 40 5 0 0 scribe the average regional response of a b lo the Earth to surface loading on time scales of ld-lo4 years.
Recommended publications
  • Geophysical Abstracts 156-159 January-December 1954
    Geophysical Abstracts 156-159 January-December 1954 GEOLOGICAL SURVEY BULLETIN 1022 Abstracts of current literature pertaining to the physics of the solid earth and geophysicq,l exploration UNITED STATES GOVERNMENT PRINTING OFFICE, WASHINGTON : 1955 UNITED STATESlDEPARTMENT OF THE INTERIOR Douglas McKay, Secretary GEOLOGICAL SURVEY W. ~· Wrather, Director CONTENTS [The letters in parentheses are those used to designate the chapters for separate publication] Page (A) Geophysical Abstracts 156, January-March------------------------ 1 (B) Geophysical Abstracts 157, April-June---------------------------- 71 (C) Geophysical Abstracts 158, July-September________________________ 135 (D) Geophysical Abstracts 159, October-December_____________________ 205 Under department orders, Geophysical Abstracts have been published at different times by the Bureau of Mines or the Geological Survey as noted below: 1-86, May 1929-June· 1936, Bureau of Mines Information Circulars. [Mimeo- graphed] 87, July-December 1936, Geological Survey Bulletin 887. 88-91, January-December 1937, Geological Survey Bulletin 895. 92-95, January-December 1938, Geological Survey Bulletin 909. 96-99, January-December 1939, Geological Survey Bulletin 915. 100-103, January-December 1940, Geological Survey Bulletin 925. 104-107, January-December 1941, Geological Survey Bulletin 932. 108-111, January-December 1942, Geological Survey Bulletin 939. 112-127, January 1943-December 1946, Bureau of Mines Information Circulars. [Mimeographed] 128-131, January-December 1947, Geological Survey Bulletin 957. 132-135, January-December 1948, Geological Survey Bulletin 959. 136-139, January-December 1949, Geological Survey Bulletin 966. 140-143, January-December 1950, Geological Survey Bulletin 976. 144-147, January-December 1951, Geological Survey Bulletin 981. 148-151, January-December 1952, Geological Survey Bulletin 991. 152-155, January-December 1953, Geological Survey Bulletin 1002.
    [Show full text]
  • The Nature of Lunar Isostasy
    45th Lunar and Planetary Science Conference (2014) 1630.pdf THE NATURE OF LUNAR ISOSTASY. Michael M. Sori1 and Maria T. Zuber1. 1Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, MA 02139, USA ([email protected]). Introduction: One way planetary topography can A key assumption in investigating the role of Pratt be supported is isostatic compensation, in which isostasy by looking at the relationship between crustal overburden pressure of rock is balanced at some depth. density and topography is that the density one observes The regions of the Moon that are not associated with near the surface is representative of the underlying maria or basins are generally isostatically compensated crustal column. We justify that assumption here by [1], an observation that was made when the first noting that calculation of the effective density of the detailed lunar gravity maps were constructed [2] and lunar crust as a function of spherical harmonic degree has held with each subsequently more precise data set considered results in a linear trend [14], supporting the [e.g., 3]. notion of a single-layer crust. There are two models of isostasy commonly Results: We make scatter plots of crustal density considered. In the Airy isostasy model [4], crustal as a function of elevation. One such scatter plot, for thickness is varied such that overburden pressures are the South Pole-Aitken basin, is shown in Figure 1. equal at some depth of compensation. The crust is a Points are sampled in a grid every ~8 km. For each layer of uniform density overlaying a mantle of higher scatter plot, we make a least-squared fit to the data and uniform density.
    [Show full text]
  • Glacial Isostatic Adjustment and Sea-Level Change – State of the Art Report Technical Report TR-09-11
    Glacial isostatic adjustment and sea-level change – State of the art report Glacial isostatic adjustment and sea-level change Technical Report TR-09-11 Glacial isostatic adjustment and sea-level change State of the art report Pippa Whitehouse, Durham University April 2009 Svensk Kärnbränslehantering AB Swedish Nuclear Fuel and Waste Management Co Box 250, SE-101 24 Stockholm Phone +46 8 459 84 00 TR-09-11 ISSN 1404-0344 CM Gruppen AB, Bromma, 2009 Tänd ett lager: P, R eller TR. Glacial isostatic adjustment and sea-level change State of the art report Pippa Whitehouse, Durham University April 2009 This report concerns a study which was conducted for SKB. The conclusions and viewpoints presented in the report are those of the author and do not necessarily coincide with those of the client. A pdf version of this document can be downloaded from www.skb.se. Preface This document contains information on the process of glacial isostatic adjustment (GIA) and how this affects sea-level and shore-line displacement, and the methods which are employed by researchers to study and understand these processes. The information will be used in e.g. the report “Climate and climate-related issues for the safety assessment SR-Site”. Stockholm, April 2009 Jens-Ove Näslund Person in charge of the SKB climate programme Contents 1 Introduction 7 1.1 Structure and purpose of this report 7 1.2 A brief introduction to GIA 7 1.2.1 Overview/general description 7 1.2.2 Governing factors 8 1.2.3 Observations of glacial isostatic adjustment 9 1.2.4 Time scales 9 2 Glacial
    [Show full text]
  • Isostasy and Flexure of the Lithosphere
    Isostasy and Flexure of the Lithosphere A. B. WATTS Department of Earth Sciences, Oxford University Oxford, United Kingdom PUBLISHED BY THE PRESS SYNDICATE OF THE UNIVERSITY OF CAMBRIDGE The Pitt Building, Trumpington Street, Cambridge, United Kingdom CAMBRIDGE UNIVERSITY PRESS The Edinburgh Building, Cambridge, United Kingdom http://www.cup.cam.ac.uk 40 West 20th Street, New York, NY 10011–4211, USA http://www.cup.org 10 Stamford Road, Oakleigh, Melbourne 3166, Australia Ruiz de Alarco´ n 13, 28014 Madrid, Spain # Cambridge University Press 2001 This book is in copyright. Subject to statutory exception and to the provisions of relevant collective licensing agreements, no reproduction of any part may take place without the written permission of Cambridge University Press. First published 2001 Printed in the United States of America Typeface Times Ten PS 10/12.5 System 3B2 [KW] A catalogue record for this book is available from the British Library Library of Congress Cataloging-in-Publication Data Watts, A. B. (Anthony Brian), 1945– Isostasy and flexure of the lithosphere / A. B. Watts. p. cm. Includes bibliographical references and index. ISBN 0-521-62272 – ISBN 0-521-00600-7 (pb) 1. Isostasy. 2. Earth – Crust. I. Title. QE511.W38 2001 551 – dc21 00-065146 ISBN 0 521 62272 7 hardback ISBN 0 521 00600 7 paperback Contents Preface page xi Acknowledgments xiii Notation xvii 1 The Development of the Concept of Isostasy 1 1.1 Introduction 1 1.2 First Isostatic Ideas 2 1.3 The Deflection of the Vertical in India 9 1.4 Isostasy According to Airy
    [Show full text]
  • Isostasy, Dynamic Topography, and the Elevation of the Apennines of Italy ∗ Claudio Faccenna A, , Thorsten W
    Earth and Planetary Science Letters 407 (2014) 163–174 Contents lists available at ScienceDirect Earth and Planetary Science Letters www.elsevier.com/locate/epsl Isostasy, dynamic topography, and the elevation of the Apennines of Italy ∗ Claudio Faccenna a, , Thorsten W. Becker b, Meghan S. Miller b, Enrico Serpelloni c, Sean D. Willett d a Laboratory Experimental Tectonics, Università Roma TRE, Roma, Italy b Department of Earth Sciences, University of Southern California, Los Angeles, CA, USA c Istituto Nazionale Geofisica e Vulcanologia, Centro Nazionale Terremoti, Bologna, Italy d Department of Earth Sciences, Swiss Federal Institute of Technology, Zurich, Switzerland a r t i c l e i n f o a b s t r a c t Article history: The elevation of an orogenic belt is commonly related to crustal/lithosphere thickening. Here, we discuss Received 19 March 2014 the Apennines as an example to show that topography at a plate margin may be controlled not only by Received in revised form 28 August 2014 isostatic adjustment but also by dynamic, mantle-driven processes. Using recent structural constraints for Accepted 14 September 2014 the crust and mantle we find that the expected crustal isostatic component explains only a fraction of Available online xxxx the topography of the belt, indicating positive residual topography in the central Apennines and negative Editor: Y. Ricard residual topography in the northern Apennines and Calabria. The trend of the residual topography Keywords: matches the mantle flow induced dynamic topography estimated from regional tomography models. Apennines We infer that a large fraction of the Apennines topography is related to mantle dynamics, producing elevation relative upwellings in the central Apennines and downwellings in the northern Apennines and Calabria subduction where subduction is still ongoing.
    [Show full text]
  • Isostasy and Its Meaning
    SVENSKA GEOFYSISKA FORENINGEN VOLUME I, NUMBER 3 Te I L u s AUGUST 1949 A QUARTERLY JOURNAL OF GEOPHYSICS Isostasy and its Meaning By B. GUTENBERG, California Institute of Technology1 (Manuscript received June 25, 1949) Abstract The theory of isostasy supposes that, in regions which have not been disturbed recently, each vertical column of the earth's crust with a certain minimum radius and extending to a depth of about IOO km has approximately the same mass. To find the deviation from this approximation in a given region, the density must be assumed as a function of depth. Such assumptions used at present for calculations are discussed critically. The resulting errors are greater than it is normally beleaved; errors in the calculated isostatic gravity anomalies exceeding ten milligals must be expected in certain regions. Systematic errors result from the usual assumption in routine calculations that the mean density in the earth's crustal layers under the bottom of the Pacific and in the continental areas is the same, and that in both the difference between the density of the layers above about 30 km and the layers below this depth is 0.6. The processes producing and maintaining isostatic equilibrium are discussed. In the theory of isostasy it is assumed that available, there are two ways to approach this in regions that have not been disturbed re- uestion. The first is, to use observations of cently each vertical column of the earth's crust tx e density and distribution of rocks in the with a given radius (at least, say, 10 km) and earth's crust and of gravity at the earth's surface extending down to a sufficient specified depth and to calculate the residuals (the so-called (apparently at least 60 km), has approximately gravity anomalies) against an assumed equili- the same mass, regardless of the surface condi- brium condition; the second is to consider the tion (continental or oceanic) or of the surface processes involved in establishing isostasy.
    [Show full text]
  • Application to the Mid-Norwegian Continental Margin
    Downloaded from http://sp.lyellcollection.org/ by guest on September 28, 2021 Isostasy as a tool to validate interpretations of regional geophysical datasets – application to the mid-Norwegian continental margin SOFIE GRADMANN1*, CLAUDIA HAASE1 &JO¨ RG EBBING1,2 1Geological Survey of Norway, Leiv Eirikssons vei 39, 7040 Trondheim, Norway 2Present address: Department of Geosciences, Kiel University, 24118 Kiel, Germany *Correspondence: sofi[email protected] Abstract: Isostasy is a well understood concept, yet rarely applied to its full capacity in regional interpretations of crustal structures. In this study, we utilize a recent density model for the entire NE Atlantic, based on refraction seismic data and gravity inversion, to calculate isostatically balanced bathymetry along the mid-Norwegian margin. Since gravity and isostatically balanced elevation are independent observables but both depend on the underlying density model, consis- tencies and discrepancies point towards model deficits, erroneously interpreted or poorly understood areas. Four areas of large isostatic residuals are identified. Along the outer Vøring Margin, a mass deficit points to more extensive high-density bodies or a shallower Moho than currently mapped. Farther seaward, along the Vøring Marginal High, a mass excess indicates inaccurate mapping of the continent–ocean boundary and surrounding structures. A number of eclogitic bodies along the proximal mid-Norwegian margin have been described in recent publications and their presence is now also confirmed by isostatic calculations. Major elevation and gravity residuals along the transition between the Vøring and Møre margins signify that the structure of this region is poorly understood and modifications to the mapped continent–ocean transition may be required.
    [Show full text]
  • Concepts of Isostasy
    Unit 3 Interior of the Earth: Structure and Composition .................................................................................................................................................................................................................................... UNIT 4 CONCEPTS OF ISOSTASY Structure 4.1 Introduction 4.4 Views on Isostasy Expected Learning Outcomes Airy’s Theory 4.2 Basic Principles Pratt’s Theory Hayford and Bowie Principle of Buoyancy Joly Density Heiskanen Lithospheric Flexure Vening-Meinesz or Flexural The Concept of Equilibrium Isostasy The Isostatic Adjustments Isostasy v/s No Isostasy 4.3 Development of the Concept of 4.5 Summary Isostasy 4.6 Terminal Questions Determination of Latitudes Isostasy as the Earth’s Balance 4.7 Answers Isostatic Effects 4.8 References/Further Reading 4.1 INTRODUCTION In the previous units of this course, you have learnt about the origin of the Earth, its internal structure and composition as well as the life forms it supports. You must have realized that it is a very complex system. Now we turn our attention to understand how and which features of the Earth helps in achieving this delicate balance. In particular, this unit is devoted to the discussion of isostasy, which gives us information about the mechanical stability between the upstanding features and low lying basins of the Earth. In this unit, we will trace the evolution of the concept of isostasy, which originated in a survey initiated to determine the shape of the Earth. 65 Block 1 GEO-Tectonics .................................................................................................................................................................................................................................... You might have visited different places in India as well as abroad and observed some relief features present on the Earth’s surface, such as mountains, plains, plateaus, valleys, lakes and oceans etc. You might have also noticed that such features are of different shapes and sizes.
    [Show full text]
  • Recent Global Sea Levels and Land Levels
    3 Recent global sea levels and land levels D.G. Aubrey and K.O. Emery Abstract Tide-gauge records from around the world ambiguously document the rise and fall of relative sea levels during the past century. Analyses of these records, non-uniformly distributed in time and space, reveal that land levels as well as ocean levels are changing, complicating the estimation of the eustatic component of sea level change. Because most tide gauges are concentrated in the Nonhero Hemisphere instead of the Southern Hemisphere (where most ocean area is located), tide-gauge data reflect major continental motions due to glacio-isostasy and neo-tectonism, which suppress and mask the signal from eustatic sea level change. Estimation of the magnitude of sea level change is important for many reasons, including the interpretation of possible effects of global climate change resulting from carbon dioxide and other trace gas loading of the atmosphere. Proper interpretation of past sea level changes is valuable for the calibration and evaluation of certain global climate, ocean and geophysical Earth models. This interpretation can be accomplished using tide-gauge records and other records of relative sea levels only if land level changes are sufficiently distinguished from actual sea level changes. Introduction Past and future climate changes drive various responses in the oceans. One such response is variation in sea levels as the oceans wann or cool, and as water is removed from or returned to the ocean basins and glaciers. During the past 15.000--18,000 years following the last glacial maximum, the Earth has experienced an average warming of about 4-5 ·c.
    [Show full text]
  • Dynamic Surface of the Earth
    MODULE - 2 Changing face of the Earth 3 Notes DYNAMIC SURFACE OF THE EARTH In the previous lesson, we have learnt that the interior of the earth is very hot. Earthquakes and volcanoes are concentrated along a few narrow belts. The type and density of rocks of the crust are variable. The surface features are dynamic in character. This dynamism is due to two forces — endogenetic and exogenetic. Endogenetic forces are those which are caused from below the surface. Due to this, an area may get elevated or gets submerged. These forces try to make the surface irregular while exogenetic froces are those which operate from above the surface. They try to eliminate the irregulari- ties of the surface through the process of denudation about which we will be reading in lesson. In this lesson we will be studying about the endogenetic forces. OBJECTIVES After studying this lesson, you will be able to : define isostasy; describe the variation in relief features on the earth’s surface; explain the isostatic adjustment by various experiments; explain the views of Airy and Pratt and distinguish between the ideas of both; explain the concept of continental drift; enumerate the evidences of continental drift; explain the concept of plate tectonics; identify and locate different plates on the world map; GEOGRAPHY 45 MODULE - 2 Dynamic Surface of the Earth Changing face of the Earth explain the mechanism of plate movement; identify various plate boundaries and associated features; explain the distribution of land and water on the globe and associate earthquakes and volcanoes with plate boundaries. Notes 3.1 CONCEPT OF ISOSTASY The term “Isostasy” is derived from “Isostasios”, a word of Greek language meaning the state of being in balance.
    [Show full text]
  • Crustal Uplift and Subsidence Due to the Interaction Between Tectonic and Surface Processes – an Integrated 3D Numerical Approach for Spatial Quantification
    Crustal uplift and subsidence due to the interaction between tectonic and surface processes – an integrated 3D numerical approach for spatial quantification Zur Erlangung des akademischen Grades eines DOKTORS DER NATURWISSENSCHAFTEN von der Fakult¨at fur¨ Physik der Universit¨at Karlsruhe (TH) genehmigte DISSERTATION von Diplom-Physiker Daniel Kurfeß aus Karlsruhe Tag der mundlichen¨ Prufung:¨ 5. Dezember 2008 Referent: Prof. Dr. Friedemann Wenzel Korreferent: Prof. Dr. Karl Fuchs Crustal uplift and subsidence due to the interaction between tectonic and surface processes – an integrated 3D numerical approach for spatial quantification by Daniel Kurfeß (E-mail: [email protected]) Geophysikalisches Institut, Universit¨at Karlsruhe (TH) 5 Contents Abstract 9 1 Introduction 11 1.1 Motivation.............................. 11 1.2 Objectivesofthisthesis . 13 2 Theoretical foundations 17 2.1 Upliftandsubsidence. 17 2.1.1 Feedback between tectonic and surface processes . 18 2.1.2 Orders of crustal displacement rates . 19 2.1.3 Definition of “uplift” and “subsidence” . 20 2.2 Isostasyanderosion......................... 21 2.2.1 Simpleerosionmodels . 21 2.2.2 Localisostasy ........................ 22 2.2.3 Regional isostasy . 23 2.3 Numerical landscape evolution modeling . 25 2.3.1 FE modeling with ABAQUS/StandardTM ........ 27 2.3.2 Surfaceprocessesmodeling. 30 6 CONTENTS 3 Coupled simulation integrating surface processes modeling and geomechanical FE modeling 35 3.1 Technique .............................. 35 3.2 Benchmarkmodels ......................... 39 3.2.1 ModelA(flexure). 39 3.2.2 ModelB(feedback). 42 3.3 Discussion.............................. 46 3.3.1 Limitations ......................... 47 3.3.2 Capabilities ......................... 47 4 Feedback process study 49 4.1 Impact of isostatic tectonic response . 49 4.2 Modelsetup............................. 50 4.3 Resultsanddiscussion ......................
    [Show full text]
  • Isostasy and Large Scale Gravity Chap. 9 Homework Answers (Dec
    Isostasy and Large Scale Gravity Chap. 9 Homework Answers (Dec. 2009) 1. Why does it seem a contradiction that the mantle can both transmit shear-waves (S- waves) and rebound isostatically? Explain on both an atomic and a bulk scale. The false paradox here is that how can the mantle both undergo ‘solid-state’ creep (i.e., flow like a liquid) due the application of loads (e.g., ice sheets, erosion) over millions of years AND transmit earthquake induced shear waves with wave periods of <100 s. The logic goes, if the earth behaves as a ‘liquid’ to long time scale load, and hence is ‘like a liquid’ internally, then since (ideal) liquids have no shear strength, then no earthquake shear waves should propagate through the earth. This was a big debate after earthquakes shear waves were finally measured in 1900! At the atomic scale (10-9 m), earthquake induced shear wave oscillations travel through the earth as ‘elastic waves’. For these short wave periods waves (<100 s), this means that the atoms can be well modeled as mass points with springs connecting the atoms as dictated by the crystal lattice arrangement of atoms. This means that shear waves can propagate through the atom ‘springs’ (atomic bonds) without much lose (dissipation) of the kinetic and strain energy associate that IS the seismic wave. Thus, the waves from an earthquake larger than magnitude 5 will be detectable by a seismometer on the other side of the planet. But, at the atomic scale, the addition or subtraction of mass loads (e.g., ice sheets, filling sediment basin) occurs very long time scales (i.e., million years) much longer than the time- scales for earthquakes waves (<100 s wave period).
    [Show full text]