Application to the Mid-Norwegian Continental Margin

Total Page:16

File Type:pdf, Size:1020Kb

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. Gold Open Access: This article is published under the terms of the CC-BY 3.0 license. The concept of isostasy has been applied to several Continental margins may be especially well regional studies worldwide, helping to estimate lith- suited for isostatic studies since the crust is thinned ospheric strength and explain the shape and origin and its deeper structure and Moho are more easily of topography (e.g. Watts 2001; Kaban et al. imaged by seismic methods than, for example, 2003). Local studies, however, which investigate those of cratonic lithosphere. Seismic mapping is the crustal structure in more detail by employing especially extensive (and accessible) along refraction and deep reflection seismic data as well hydrocarbon-bearing continental margins such as as gravity modelling, often ignore the additional the mid-Norwegian one. Furthermore, thinned insight that could be gained from isostatic consider- crust and lithosphere are generally weaker, implying ations. Positive exceptions to this are, for example, that lateral density variations and Moho undulations Watts et al. (1985), Zeyen & Ferna`ndez (1994), will have a larger expression in the bathymetry or Mat¸enco et al. (1997), Hirsch et al. (2009) and basin shape. Thermal weakening during rifting and Fullea et al. (2010), where isostasy has entered the breakup as well as rigidification during subsequent model building and interpretation. Methods for cooling additionally affect the flexure of the margin. joint analysis and inversion of gravity and isostasy The enormous amount of geophysical and geo- have, for example, been developed by Zeyen & Fer- logical data that exists for the region of the NE na`ndez (1994), Stewart et al. (2000), Braitenberg Atlantic has been gathered, analysed and compiled et al. (2004) and Salem et al. (2014). in recent years into a comprehensive tectonostrati- Only a few of the studies involving isostasy graphic atlas (the NAG-TEC Atlas, Hopper et al. cover the mid-Norwegian margin (Kimbell et al. 2014). Based on these data, a new seismic-derived 2004; Ebbing 2007) and we show in this paper that Moho grid has been compiled (Funck et al. 2016) simple isostatic considerations can also improve the and subsequently further improved using gravity understanding of the crustal structure. inversion (Haase et al. 2016). The latter study also From:Pe´ron-Pinvidic, G., Hopper, J. R., Stoker, M. S., Gaina, C., Doornenbal, J. C., Funck,T.&A´ rting,U.E. (eds) The NE Atlantic Region: A Reappraisal of Crustal Structure, Tectonostratigraphy and Magmatic Evolution. Geological Society, London, Special Publications, 447, https://doi.org/10.1144/SP447.13 # 2017 The Author(s). Published by The Geological Society of London. Publishing disclaimer: http://www.geolsoc.org.uk/pub_ethics Downloaded from http://sp.lyellcollection.org/ by guest on September 28, 2021 S. GRADMANN ET AL. constructed a first-order density model over the formation in the early Eocene (c. 54 Ma, Talwani entire NE-Atlantic, which we utilize here to forward & Eldholm 1977; Skogseid & Eldholm 1987). calculate an isostatically consistent bathymetry that This volcanic activity was part of the North Atlantic reflects all internal loading of the model. The result- Igneous Province (Eldholm & Grue 1994; Meyer ing differences from the true bathymetry enable us et al. 2007) and is reflected along the present-day to detect consistencies and discrepancies between mid-Norwegian outer margin by extensive regions the datasets mapping the deep subsurface structure. with seaward-dipping reflectors and massive vol- Since gravity and isostatically balanced elevation canic flows. Magmatic activity also led to intrusions are independent observables but both depend on into the thinned continental crust and underplating the underlying density model, consistencies and dis- in the continent–ocean transition zone. Oceanic crepancies can point towards poorly understood or spreading along the Aegir and Mohn’s ridges led erroneously interpreted areas. The obtained eleva- to the complete separation of the Greenland and tion residuals are then discussed with respect to Eurasian plates. Continued, contemporaneous existing seismic profiles and structural interpreta- Cenozoic rifting in East Greenland and associated tions along the margin. By considering residuals breakup in late Oligocene (c. 25–24 Ma) led to of shorter wavelengths only, we focus on the crustal the formation of a second spreading axis and to structure, which is much better mapped and dis- the separation of the Jan Mayen microcontinent cussed than the underlying lithospheric mantle. from Greenland (Talwani & Eldholm 1977; Vogt Previous studies that investigated regional isostatic et al. 1980; Nunns 1982; Gudlaugsson et al. 1988; and gravity modelling in the NE Atlantic region Gaina et al. 2009; Gernigon et al. 2015). Today, (Kimbell et al. 2004; Reynisson et al. 2010) had the Kolbeinsey Ridge is the active spreading centre less complete datasets at hand and focussed on linking the Reykjanes Ridge in the south to the somewhat different questions. This study merely Mohn’s Ridge in the north. tests an existing subsurface model and does not A large number of seismic refraction profiles attempt to improve it. cover the mid-Norwegian margin but become more sparse to the north (Fig. 1). Some profiles are several decades old and large receiver spacing Geological history of the mid-Norwegian (10–20 km) often prevents the detailed mapping margin of smaller-scale features. Nevertheless, this exten- sive grid predominantly shaped our present-day The mid-Norwegian margin, like any passive rif- understanding of the deep margin structures. ted margin, consists of numerous rifted blocks The 300 km wide Møre Margin is characterized and interjacent sedimentary basins and sub-basins by a narrow platform and a rapid thinning of the (Fig. 1). The margin is subdivided into the three continental crust. The main area of crustal thinning major segments of the Møre, Vøring and Lofoten is confined to a narrow zone with a width of only a margins. The large shelf areas of the North Sea few tens of kilometres (Olafsson et al. 1992; Raum and Barents Sea extend to the south and north, 2000; Mjelde et al. 2008; Osmundsen & Ebbing respectively. The structure of the mid-Norwegian 2008; Kvarven et al. 2014; Nirrengarten et al. margin reflects a long period of extension following 2014). West of this major necking zone lies the the orogenic collapse of the Caledonian orogeny. nearly 200 km wide and long Cretaceous Møre Earliest extension during the post-orogenic collapse Basin, which reaches depths of more than 15 km. initiated in Devonian–Carboniferous times, form- This basin is underlain by highly thinned crust (in ing large detachment faults and basins that are places less than 10 km thick) and high-velocity, now partly exposed in onshore areas (e.g. Seranne high-density lower crustal bodies (LCBs) have & Seguret 1987; Osmundsen & Andersen 2001, been mapped both outbound and inbound of the and references therein) and possibly buried by youn- necking zone but not directly underneath it (Olafs- ger sediments in the offshore domain. Carbonifer- son et al. 1992; Kvarven et al. 2014; Nirrengarten ous to early Mesozoic extension created the early, et al. 2014). The outer limit of the margin compri- proximal rift basins (e.g. the Trøndelag Platform; ses the Møre Marginal High which is overlain by a Brekke 2000). Distal rift basins formed during Mid- series of volcanic flows. Whereas some authors dle Jurassic to Early/middle-Cretaceous (e.g. the interpret the basement of the Møre Marginal High Lofoten, Vøring and Møre basins; Brekke 2000). as thick continental crust (e.g. Reynisson et al. Late Cretaceous to Paleocene extension shaped the
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]
  • 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]
  • Toolbox for Analysis of Flexural Isostasy (TAFI)—A MATLAB Toolbox for Modeling Flexural Deformation of the Lithosphere GEOSPHERE; V
    Software Contribution GEOSPHERE Toolbox for Analysis of Flexural Isostasy (TAFI)—A MATLAB toolbox for modeling flexural deformation of the lithosphere GEOSPHERE; v. 13, no. 5 Sumant Jha, Dennis L. Harry, and Derek L. Schutt Department of Geosciences, Colorado State University, 1482 Campus Delivery, Fort Collins, Colorado 80523, USA doi:10.1130/GES01421.1 7 figures; 5 tables; 2 supplemental files ABSTRACT typically involving subsidence in the vicinity of the load (assuming the load represents an addition of weight to the lithosphere) and uplift in surrounding CORRESPONDENCE: sumant .jha@colostate .edu This paper describes a new interactive toolbox for the MATLAB computing areas (the peripheral “bulge”) (Fig. 1). The magnitude and wavelength of the platform called Toolbox for Analysis of Flexural Isostasy (TAFI). TAFI supports flexural deflection depends on the magnitude of the load and the strength of CITATION: Jha, S., Harry, D.L., and Schutt, D.L., two- and three-dimensional (2-D and 3-D) modeling of flexural subsidence the lithosphere, with greater strength resulting in more broadly distributed 2017, Toolbox for Analysis of Flexural Isostasy (TAFI)—A MATLAB toolbox for modeling flexural and uplift of the lithosphere in response to vertical tectonic loading. Flexural (larger-wavelength) deformation. Flexure results in variations in the regional defor ma tion of the lithosphere: Geosphere, v. 13, deformation is approximated as bending of a thin elastic plate overlying an gravity field due to deformation of density interfaces within the lithosphere. no. 5, p. 1555–1565, doi:10.1130/GES01421.1. inviscid fluid asthenosphere. The associated gravity anomaly is calculated by Such density interfaces may include the crust-mantle boundary, the interface summing the anomalies produced by flexure of each density interface within between the sedimentary basin fill and underlying rocks, the seafloor in marine Received 26 August 2016 the lithosphere, using Parker’s algorithm.
    [Show full text]