Deriving rock uplift histories from data-driven inversion of river profiles
|
|
- Leslie Parker
- 5 years ago
- Views:
Transcription
1 GSA DATA REPOSITORY Deriving rock uplift histories from data-driven inversion of river profiles Christoph Glotzbach 1 1 Institute of Geology, Leibniz University Hannover, Callinstr. 30, Hannover, Germany, glotzbach@geowi.uni-hannover.de DATA REPOSITORY DR1: Methodological details The rate of change in elevation of the river profile depends on the rock uplift U and erosion rate E e.g., Whipple and Tucker, 1999), and thus the evolution equation can be written as: z 1) t = U + E x,t x,t The lower cases in equation 1) indicate that the uplift and erosion rate can vary with time t) and location x). Erosion counteracts uplift and degrades uplifted river profiles. Based on geomorphological investigations, fluvial erosion laws have been implemented, such as the transport-limited or detachment-limited model e.g., Whipple and Tucker, 1999). Here the so-called stream power model is used, which assumes that the rate of incision is primarily controlled by discharge and channel slope, which is used as a proxy for stream power e.g. Whipple and Tucker 1999). The erosional coefficient K) is laterally and temporally constant in the presented applications, but the program also allows to vary K laterally. Temporal variations of discharge and K induced by orbital cycles have little effect on the presented river profile modelling, because the length and timescales of the model is much larger than the climate variations e.g. Paul et al. 2014). The concavity m r /n r ) and exponent of the river gradient n r ) are free parameters in all models and allowed to vary from and from 2/3 to 5/3, respectively e.g., Whipple, 2001). All models start with flat river profiles, which are set to the present-day outlet elevation of individual trunk streams. The latter is invariant during modelling and used as model boundary condition. The presented model enables to take into account the flexural isostatic uplift as a result of changes in topography. The isostatic uplift is modelled by solving the biharmonic equation similar to the approach in PECUBE Braun et al., 2012): D 4 Δu x 4 + D 4 Δu y 4 + 2D 4 Δu x 2 y 2 = ΔρgΔu + ρ 0gΔe 2) where D is the lateral constant flexural rigidity given by: ET e 3 D = 121 v 2 ) T e is the equivalent elastic thickness of the plate, E is the Young s modulus, v is the Poisson s ratio, g is the acceleration due to gravity and Δe and Δu are the increment in
2 surface topography and isostatic uplift over a time step Δt of computation. An efficient Fourier-transform approach was used to solve the so-called thin elastic plate equation Nunn and Aires, 1988). The isostatic uplift is added to the input rock uplift and thus either increase or decrease the overall rock uplift. The complex interaction between the transient topographic evolution and the input uplift history result in complex uplift paths that show distinct differences compared to the step-like uplift paths modelled without isostasy cf. Fig. DR3D and Fig. DR4D). The misfit of modelled and observed data is regulating the inversions, and is calculated with the log-likelihood lnl) function for Gaussian-distributed errors: n 2 ln2π) o lnl = + lnσ 2 i i p i 4) i=1 σ i where p i is the predicted value elevation, thermochronological age, nuclide concentration), and o i and σ i is the observed value and 1σ error of observation i and n is the total number of observations. The lnl of each observation class elevation of river nodes, cosmogenic nuclide concentrations, thermochronological ages and fission track length distributions) is normalised divided) by the number of observations and the sum of normalised lnls regulates the inversion, avoiding preferential fitting of the largest dataset. Analytical uncertainties of thermochronological and cosmogenic nuclide data are well known, ranging from approximately 5 to 20%. As an example the mean errors of the apatite fission track ages and 10 Be concentrations from the Sila massif are 10 and 6%, respectively. The vertical absolute height error is less than 16 m for 90% of the global shuttle radar topography mission SRTM) data e.g. Rabus et al. 2003), but river profiles extracted from steeply incised streams will be associated with much larger errors. I used an arbitrary height error of 75 m, which corresponds to a mean relative error of 10%. The neighbourhood algorithm inversion was used to efficiently explore the mulidimensional parameter space and extract statistical parameters Sambridge 1999a,b). This study used the parameterisation suggested by Glotzbach et al. 2011). Three main parameters control the convergence behaviour of the inversion: i) the number of iterations N i ), ii) the number of models generated at each iteration N m ) and iii) the number of Voronoi cells re-sampled at each iteration N r ). In this study N i is 200 inversion test) and 400 Sila massif), N m is 100 and N r is 90. The rather high resampling ratio 0.9) led to a rather slow convergence and ensures that the full parameter space is explored. Bayesian estimates such as 1D marginal PDFs are first derived separately for each dataset and finally they are combined with: PC D 1,D 2,D 3 ) =a n =3 n =1 PC D n ) PC) 2 5) where PC) is the prior probability of parameter C and PC D n ) are the posterior probabilities of independent datasets, here the river profile, the cosmogenic and thermochronological data. River long profiles were extracted from shuttle radar topography mission SRTM) digital elevation models DEM) with a vertical accuracy of ±16 m. Only those stream
3 segments with an accumulation area >1 km 2 were extracted from the DEM, and smoothed with a moving average window of 1 km. The resulting river profile is resampled to 1 km distance between successive river nodes. The Sila massif model consists of 129 river nodes, and each river node is connected to a 500 m long hillslope discretised into 10 nodes 50 m spacing). The computation time of a single forward model of the Sila massif is 16.4 CPU seconds on a 2.66 GHz Quad-Core Intel Xeon), whereas the solution of the riverhillslope model takes 7 seconds, the calculation of 10 Be concentration at all river and hillslope nodes takes 8 seconds and the 1D thermal modelling takes 1.4 seconds. Increasing the number of river nodes increases only sligthly the computational time of the river-hillslope model Fig. DR1); the calculation of 261 river nodes takes 10 seconds compared to 7 seconds for 129 river nodes CPU seconds y = 0.025x Number of river nodes Fig. DR1: Solution time of a the river-hillslope model for a single forward model as a function of the number of river nodes.
4 DR2: Testing the performance of the model Two synthetic datasets river profile, cosmogenic and thermochronological data) were generated with 1.) increasing uplift rates and 2.) decreasing uplift rates. In both examples, river nodes were extracted from the Acher catchment in the northern Black Forest in southwestern Germany. 1.) The synthetic data is modelled with an uplift rate of 0.1 km/ma from 50 to 1 Ma, followed by 0.5 km/ma uplift until present. The river profile was modeled assuming that river incision is proportional to stream power and the erosional coefficiency K) and exponents of discharge m) and slope n) are 1*10-6 yr -1, 0.5 and 1, respectively. The resulting river profile is characterised by a downstream increase in steepness with a distinct knickzone at 6-10 km upstream of the bounding node Fig. DR2A). Recorded erosion rates and cooling paths were used to calculate at all model nodes 63 river nodes and 630 hillslope nodes) the concentration of cosmogenic 10 Be and at three locations apatite U-Th)/He ages and apatite fission track ages with corresponding length distributions Fig. DR2A-C and DR3A-C). The synthetically generated analytical data was used for inverse modeling of seven parameters: P1) the time of change in uplift rate, P2-3) the uplift rate before and after P1, P4) the erosional coefficient K), P5) the exponent of slope n), P6) the concavity index m/n) and P7) the geothermal gradient. Convergence of the likelihood function is achieved after running a few thousands models in a few hours. Inversion results obtained by inverting only the river profile without analytical data reveal that single parameters, such as the exponent n and the concavity index m/n) can be successfully retrieved Fig. DR2D). The complete uplift history, especially the timing of the increase in uplift rate, however, can only be reconstructed by the integrated inversion of all datasets including the thermochronological and cosmogenic nuclide data Fig. DR2D).
5 O+.PF.CDB G3.-45,1C/H9I * % ) L ) /:4903. ;1:91/+,-.+/:4903. ;1:91/+,<=./:4903.?.:5@2,5/+,-.+/0+12,3. ABC58.5,D/+,-.+/0+12,3. E,33:310./>/+,<=./D+.:5 61<.3/C1<. +,-.+/0+12,3./45//64 F03,25 %,C,5,43/0+12,3./45//64.+1:,1C > % J L=./H64I O+.PF.CDB J,:54CD./2+19/K1FC<4+B/H9I ) * N N N N N N N N N N N N N % / :4903./C49. 7;M/H /45N>=I N N N N N N ) % * ;1:91/:4903./C49. % 7,9./HZ4I N) N N N N QRN QRN QRN QRN N N QRN QRN N N AJRN AJRN AJRN AJRN* N N N AJRN) N AJRN N QR*N QRN QRN% QRNS@ N QRN N QRN* QRN) N AJRN N AJRN AJRN AJRNS@ N N N AJRN N N AJRN AJRN N N N N N N N N % N N N N% N N N N% % N) N% N* N N N N N N N) 7,9./H64IT/U V03,25/H99>4IT/U V03,25/H99>4IT/U IT/U GW01C.C5/CT/U ;1CDN/,C<.W/9>CT/U XX/HY;>Z9IT/U) G+1:,1C/+45./H99>4I F+./HY;I N N N N N N N N N N %? ; 7,9./H64I Fig. DR2: Results of inverse modeled synthetic river profile, thermochronological and cosmogenic 10 Be data. Synthetic data was modeled with an erosional coefficiency K) and exponents of discharge m) and slope n) of 1*10-6 yr -1, 0.5 and 1, and 0.1 km/ma uplift from 50 to 1 Ma, followed by 0.5 km/ma until present. A) Synthetic river profile, apatite fission track age and length distribution, apatite U-Th)/He age and cosmogenic 10 Be nuclide concentrations of sample localities. The best-fit modeled data is shown in comparison to the synthetic data. B) Erosion rate evolution of modeled river nodes of the trunk stream for the last 1.2 Ma. C) Modeled cooling curves of three sample localities along the trunk stream used to calculate thermochronological data. D) 1D marginal distributions of forward models with their mean value and standard deviations. The black distributions are derived from an inversion based on all data river profile plus analytical data) and the red distributions are from an inversion based solely on the river profile. 2.) The synthetic data is modelled with an uplift rate of 1.0 km/ma from 20 to 1 Ma, followed by 0.1 km/ma uplift until present. The river profile was modeled assuming an erosional coefficiency K) of 1*10-7 yr -1 and exponents of discharge m) slope n) of 0.5 and 1.25, respectively. The resulting river profile is characterised by a upstream increase in steepness Fig. DR3A). Erosion rates and cooling paths were extracted from similar sample locations as in first example and used to calculate the concentration of cosmogenic 10 Be and apatite U-Th)/He ages and apatite fission track ages with corresponding length distributions Fig. DR3A-C). The synthetically generated analytical data was used for inverse modeling of seven parameters: P1) the time of change in uplift rate, P2-3) the uplift rate before and after P1, P4) the erosional coefficient K), P5) the exponent of slope n), P6) the concavity index m/n) and P7) the geothermal gradient. Convergence of the likelihood function is achieved after running a few thousands models. Extracted 1D marginal distributions of the
6 posterior probability distribution predict most of the free parameter correctly, except for the erosional coefficient k) and the exponent of slope n r ). The uplift history, however, is correctly retrieved by the inverse modelling Fig *7*89 K : : : : : : NQ,-./01 S4,TO,EB9 F : : : : )P,4+5- A0+<7,-E0+, Ø=1.13 SD=0.37 )GR-. -08?Q1 : F % L G5A+5-A0+<7,-E0+, : : : : : : F : %? Ø=1.01 SD=0.14 4*M,4-<45=*7,-08--/0 % K*A80EB,-=45+-65OEH L *E*8*07-<45=*7,-08--/0 : N Ø=0.15 ),+<,408O4,-.JG1 L % F : : : : U : : : G % )*+,-.>01 % )*+,-./01 Ø=2.8*10-6 SD=1.8* *7*89 % )*+,-./012-3 Ø=0.80 SD=0.07 : : ;<7*=8-.>+?/012-3 Ø=0.45 SD=0.03 : ;<7*=8-.>+?/012-3 % Ø=28.3 SD=4.4 ;<7*=8-.>+?/01 D 12-3 : : : *7*89-.C D : : F G5EB-*EH,C-+ 4?E : : % II-.JG?>+12-3F )*+,-./01 Fig. DR3: Results of inverse modeled synthetic river profile, thermochronological and cosmogenic 10 Be data. Synthetic data was modeled with an erosional coefficiency K) and exponents of discharge m) and slope n) of 1*10-7 yr -1, 0.5 and 1.25, and 1.0 km/ma uplift from 20 to 1 Ma, followed by 0.1 km/ma until present. A) Synthetic river profile, apatite fission track age and length distribution, apatite U-Th)He age and cosmogenic 10 Be nuclide concentrations of sample localities. The best-fit modeled data is shown in comparison to the synthetic data. B) Erosion rate evolution of modeled river nodes of the trunk stream for the last 1.2 Ma. C) Modeled cooling curves of three sample localities along the trunk stream used to calculate thermochronological data. D) 1D marginal distributions of forward models with their mean value and standard deviations.
7 DR3: Inversion results of the Sila massif taking into account isostasy D9/E,E.).-F D9/E,E.).-F < )*+,-./01234 AG A A % %.0.-.,)19.+*91>9/8.)* G % % C H % C H B,3>)*1BAAA B,3>)*1BAAA = IV*93/7V9/0/)/UF16,3>)* P/63/U* *916,3>)* WE6*9+*Q19.+*91>9/8.)* X.-V1).-V/)/UF =*6-O *91>9/8.)*?@A B5@A% AG A A B% 5.6-,07*189/318,:)-12;34?@A B5@A% BOH B BC B A A A A% B C =*17/07*0-9,-./0 2 % 1,-AMU4 %?@A% B5@A% B,3>)*1I9./0-/AAA C E <YI1,U*12J,4?@AC B5@A?@A B5@A%G G G A A A A AG A% A A A I.3*12J,4K1D I.3*12J,4K1D L>) M,4K1DG L>) M,4K1D% L>) M,4K1D A A A% A?@AH A?@A%C A?@AG B5@A A A A B5@A B5@AC A A A A A A A A A A A A A A% A% A% A%?@AT O A A% A% A A B5@AT O % A A A% A% A A AC G G % % 9/6A17/*88A12N O 4K1D N>/0* K1DC P/07A1.0Q*N13 9 M0 9 K1D RR12SPM;34K1DH I.3*12J,4 Fig. DR4: Results of inverse modelled Trionto river profile, thermochronological and cosmogenic 10 Be data, taking into account uplift caused by flexural isostasy with the following parameters: elastic thickness of 5 km, density of the crust and mantle 2700 and 3200 kg/m 3, Young s modulus 1*10 11 Pa, Poisson ratio of A) Observed and best-fit river profile, apatite fission track age, mean apatite fission track length MTL) and cosmogenic 10 Be nuclide concentrations of sample localities. The observed river profile is color coded according to lithology green: metamorphic rocks, blue: Mesozoic-Eocene sediments, red: granite). The best-fit model yield a normalized misfit of 19 for the following parameters P1-P9): 16.3 Ma, 0.86 Ma, 1.90 mm/a, 0.07 mm/a, 1.45 mm/a, 4.12*10-6, 0.69, 0.63, 21.3 C/km. B) 1D marginal distributions of the posterior probability distribution with their mean values and standard deviations. The lower right diagram shows the 2D marginal distribution of the uplift history derived by combining parameters P1-P5. The thick black and red line are the total rock uplift input uplift plus uplift caused by changes in topography) for the best-fit model for the lowest red) and higher black) river node. C) 2D marginal distributions of the posterior probability distribution of parameter pairs that show obvious tradeoffs. A% AG A A E P P/07A1.0Q*N13 9 M0 9 K1D L>) M,4K1DG L>) M,4K1D A A A% A A A % I.3*12J,4K1D A A A% A A D9/E,E.).-FD9/E,E.).-F A A A A A A A% L>).8-12;3MJ,4 A A A I.3*12J,4K1D A A N>/0* K1DC A A A A A A% A D9/E,E.).-F Probability *10-3 )
8 DR4: River profile evolution of the Trionto river, Sila massif, Italy % +,-./ ) * ) % 607/28-9:14;*<081,9=,/3>45 Fig. DR5: Modeled river profile at 20, 15, 10, 5, 1, 0.5, 0 and in 2 Ma of the Trionto river for the best-fit model with following parameters: i) an erosional coefficient k) of 4.36*10-6, ii) an exponent n of 1.25, iii) a concavity index m/n) of 0.40, iv) a geothermal gradient of 25.5 C/km and v) an initial uplift of 1.5 km/ma until 17.5 Ma, followed by 0.08 km/ma uplift until 0.9 Ma and finally 0.78 km/ma uplift. References cited in the Data Repository Braun, J., van der Beek, P.A., Valla, P., Robert, X., Herman, F., Glotzbach, C., Pedersen, V., Perry, C., Simon-Labric, T., and Prigent, C., 2012, Quantifying rates of landscape evolution and tectonic processes by thermochronology and numerical modeling of crustal heat transport using PECUBE. Tectonophysics, v , Glotzbach, C., van der Beek, P.A., and Spiegel, C., 2011, Episodic exhumation and relief growth in the Mont Blanc massif, Western Alps from numerical modelling of thermochrology data. Earth Planet. Sci. Lett., v. 304, Nunn, J., and Aires, J., 1988, Gravity anomalies and flexure of the lithosphere at the Middle Amazon Basin, Brazil. J. Geophys. Res., v. 93, Paul, J.D., Roberts, G.G., and White, N., 2014, The African landscape through space and time. Tectonics, v. 33, doi: /2013tc Rabus, B., Eineder, M., Roth, A., and Bamler, R., 2003, The shuttle radar topography mission a new class of digital elevation models acquired by spaceborne radar. ISPRS J. Photogramm. Remote Sens., v. 57, Sambridge, M., 1999a, Geophysical inversion with a neighbourhood algorithm I. Searching a parameter space. Geophys. J. Int., v. 138, Sambridge, M., 1999b, Geophysical inversion with a neighbourhood algorithm II.
9 Appraising the ensemble. Geophys. J. Int., v. 138, Whipple, K.X., and Tucker, G.E., 1999, Dynamics of the stream-power river incision model: Implications for height limits of mountain ranges, landscape response timescales, and research needs. J. Geophys. Res., v. 140, p Whipple, K.X., 2001, Fluvial landscape response time: how plausible is steady-state denudation? Am. J. Sci., v. 301, p
GSA DATA REPOSITORY Sternai et al. 1. Algorithm Flow Chart
GSA DATA REPOSITORY 2012311 Sternai et al. 1. Algorithm Flow Chart Figure DR1: Flow chart of the algorithm to further clarify the calculation scheme. 2. Reconstruction of the Pre-Glacial Alpine Topography
More informationDown-stream process transition (f (q s ) = 1)
Down-stream process transition (f (q s ) = 1) Detachment Limited S d >> S t Transport Limited Channel Gradient (m/m) 10-1 Stochastic Variation { Detachment Limited Equilibrium Slope S d = k sd A -θ d S
More informationEffects of transient topography and drainage basin evolution on detrital thermochronometer data
UNIVERSITY OF MICHIGAN Effects of transient topography and drainage basin evolution on detrital thermochronometer data Contents Acknowledgments...3 Abstract...4 1. Introduction...5 2. Model setup...6 2.1
More informationConstraints on the rate of post-orogenic erosional decay from low-temperature thermochronological data: application to the Dabie Shan, China
Earth Surface Processes and Landforms Post-orogenic Earth Surf. Process. erosional Landforms decay 30, 1203 1225 (2005) 1203 Published online in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/esp.1271
More informationEarth and Planetary Science Letters
Earth and Planetary Science Letters 295 (2010) 511 522 Contents lists available at ScienceDirect Earth and Planetary Science Letters journal homepage: www.elsevier.com/locate/epsl Inversion of thermochronological
More informationTopographic metrics and bedrock channels Outline of this lecture
Topographic metrics and bedrock channels Outline of this lecture Topographic metrics Fluvial scaling and slope-area relationships Channel steepness sensitivity to rock uplift Advancing understanding 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 informationZeumann and Hampel, 2017, Impact of Cocos Ridge (Central America) subduction on the forearc drainage system: Geology, doi: /g
GSA Data Repository 2017296 Zeumann and Hampel, 2017, Impact of Cocos Ridge (Central America) subduction on the forearc drainage system: Geology, doi:10.1130/g39251.1. DESCRIPTION OF CASQUS To implement
More informationLong-term landscape evolution: Linking tectonics and surface processes. Paul Bishop Department of Geographical and Earth Sciences
Long-term landscape evolution: Linking tectonics and surface processes Paul Bishop Department of Geographical and Earth Sciences Collaboration PhDs: Daniel Campanile, Miguel Castillo, Tibi Codilean, Kate
More informationConstraining the stream power law: a novel approach combining a landscape evolution model and an inversion method
doi:10.5194/esurf-2-155-2014 Author(s) 2014. CC Attribution 3.0 License. Earth Surface Dynamics Open Access Constraining the stream power law: a novel approach combining a landscape evolution model and
More informationBasin-scale analysis of long-term sediment-generation rates derived from 10 Be in river sediment:
Basin-scale analysis of long-term sediment-generation rates derived from 10 Be in river sediment: The Susquehanna River basin and beyond Joanna M. Reuter M.S. Proposal April 30, 2003 Paul Bierman, advisor
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 informationPreservation or Piracy: Diagnosing low relief, high elevation surface formation mechanisms. Supplemental Methods: Additional methodological details.
GSA Data Repository 2017023 Preservation or Piracy: Diagnosing low relief, high elevation surface formation mechanisms Kelin X. Whipple 1, Roman A. DiBiase 2, William B. Ouimet 3, Adam M. Forte 1 Contents
More informationData Repository. Spatiotemporal trends in erosion rates across a pronounced rainfall gradient: examples from the south central Andes
Data Repository Spatiotemporal trends in erosion rates across a pronounced rainfall gradient: examples from the south central Andes Bodo Bookhagen 1 and Manfred R. Strecker 2 1 Geography Department, Ellison
More informationDrainage Basin Geomorphology. Nick Odoni s Slope Profile Model
Drainage Basin Geomorphology Nick Odoni s Slope Profile Model Odoni s Slope Profile Model This model is based on solving the mass balance (sediment budget) equation for a hillslope profile This is achieved
More informationWeek 3 (Feb 12): Erosion and Sediment Transport Discussion Leader: Ariel Deutsch
Week 3 (Feb 12): Erosion and Sediment Transport Discussion Leader: Ariel Deutsch The papers this week explore the topics of erosion and sediment transport, with a major theme revolving around climate-tectonic
More informationMass Wasting and Landscape Evolution
Mass Wasting and Landscape Evolution 11-8-06 Uplift is a tectonic process Three types of uplift: 1. Collisional uplift 2. isostatic uplift 3. Extensional uplif. A physical experiment in isostasy: [crust
More informationSupplement of The influence of upper-plate advance and erosion on overriding plate deformation in orogen syntaxes
Supplement of Solid Earth, 9, 127 1224, 218 https://doi.org/1.5194/se-9-127-218-supplement Author(s) 218. This work is distributed under the Creative Commons Attribution 4. License. Supplement of The influence
More informationSensitivity analysis and parameter space exploration
GSA Data Repository 2018013 Gray et al., 2018, Off-fault deformation rate along the southern San Andreas fault at Mecca Hills, southern California, inferred from landscape modeling of curved drainages:
More informationGeodynamics Lecture 8 Thermal processes in the lithosphere
Geodynamics Lecture 8 Thermal processes in the lithosphere Lecturer: David Whipp david.whipp@helsinki.fi 25.9.2014 Geodynamics www.helsinki.fi/yliopisto 2 Goals of this lecture Introduce time dependence
More informationAn integral approach to bedrock river profile analysis
1 An integral approach to bedrock river profile analysis 2 3 J. Taylor Perron and Leigh Royden 4 5 6 Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, Cambridge
More informationDetermination of uplift rates of fluvial terraces across the Siwaliks Hills, Himalayas of central Nepal
Determination of uplift rates of fluvial terraces across the Siwaliks Hills, Himalayas of central Nepal Martina Böhme Institute of Geology, University of Mining and Technology, Freiberg, Germany Abstract.
More informationIntro to Quantitative Geology
Introduction to Quantitative Geology Lesson 13.2 Low-temperature thermochronology Lecturer: David Whipp david.whipp@helsinki.fi 4.12.17 3 Goals of this lecture Define low-temperature thermochronology Introduce
More informationKarl Mueller QUATERNARY RIFT FLANK UPLIFT OF THE PENINSULAR RANGES IN BAJA AND SOUTHERN CALIFORNIA BY REMOVAL OF MANTLE LITHOSPHERE
QUATERNARY RIFT FLANK UPLIFT OF THE PENINSULAR RANGES IN BAJA AND SOUTHERN CALIFORNIA BY REMOVAL OF MANTLE LITHOSPHERE August 2013 Quaternary Rift Flank Uplift of the Peninsular Ranges in Baja and Southern
More informationAmerican Journal of Science
[American Journal of Science, Vol. 301, April/May, 2001, P.313 325] American Journal of Science APRIL/MAY 2001 FLUVIAL LANDSCAPE RESPONSE TIME: HOW PLAUSIBLE IS STEADY-STATE DENUDATION? KELIN X. WHIPPLE
More informationQUANTITATIVE THERMOCHRONOLOGY
QUANTITATIVE THERMOCHRONOLOGY Thermochronology, the study of the thermal history of rocks, provides an important record of the vertical motions of bodies of rock over geological timescales, enabling us
More informationAccelerated extension of Tibet linked to the northward underthrusting of Indian crust
Accelerated extension of Tibet linked to the northward underthrusting of Indian crust Richard Styron*, Michael Taylor, and Kurt Sundell richard.h.styron@gmail.com DOI: 1.138/NGEO2336 Methods summary The
More informationStratigraphic modelling using CHILD
5 Stratigraphic modelling using CHILD 5.1 Triangular irregular network Surface process models are widely used in geomorphology and geology, and the developments in the field follow each other rapidly.
More informationmountain rivers fixed channel boundaries (bedrock banks and bed) high transport capacity low storage input output
mountain rivers fixed channel boundaries (bedrock banks and bed) high transport capacity low storage input output strong interaction between streams & hillslopes Sediment Budgets for Mountain Rivers Little
More informationSelected Presentation from the INSTAAR Monday Noon Seminar Series.
Selected Presentation from the INSTAAR Monday Noon Seminar Series. Institute of Arctic and Alpine Research, University of Colorado at Boulder. http://instaar.colorado.edu http://instaar.colorado.edu/other/seminar_mon_presentations
More informationPIERRE G. VALLA Lycée du Parc, Lyon, France Classes for high schools competitions (Biology, Chemistry/Physics and Earth Sciences
PIERRE G. VALLA DOCTOR IN EARTH SCIENCES Earth Surface Dynamics - Geological Institute Swiss Federal Institute of Technology - ETH Zurich Sonneggstrasse 5-8092 Zurich - Switzerland Phone: +41 44 632 36
More informationERS Track 98 SAR Data and InSAR Pairs Used in the Analysis
ERS Track 98 SAR Data and InSAR Pairs Used in the Analysis Date 1 Date 2 Date 1 Date 2 Date 1 Date 2 Date 1 Date 2 7/17/1992 6/19/2000 7/17/1992 7/2/1993 9/10/1993 10/28/1996 9/3/1995 10/18/1999 9/25/1992
More informationFRACTAL RIVER BASINS
FRACTAL RIVER BASINS CHANCE AND SELF-ORGANIZATION Ignacio Rodriguez-Iturbe Texas A & M University Andrea Rinaldo University of Padua, Italy CAMBRIDGE UNIVERSITY PRESS Contents Foreword Preface page xiii
More informationMacmillan Publishers Limited. All rights reserved
LETTER doi:.8/nature8 Lifespan of mountain ranges scaled by feedbacks between landsliding and erosion by rivers David L. Egholm, Mads F. Knudsen & Mike Sandiford An important challenge in geomorphology
More informationGeomorphology: Mechanics and Evolution of Landscapes (GES )
Geomorphology: Mechanics and Evolution of Landscapes (GES 206.2.5401) Instructor: Itai Haviv, haviv@bgu.ac.il, room 331, building 58, office hours: Tuesday 11:15-13:15 Teaching assistant: May-Tal Sadeh,
More informationModeling fluvial incision and transient landscape evolution: Influence of dynamic channel adjustment
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 113,, doi:10.1029/2007jf000893, 2008 Modeling fluvial incision and transient landscape evolution: Influence of dynamic channel adjustment M. Attal, 1 G. E. Tucker,
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 informationWeathering, Landscape Evolution, Supergene Enrichment, and Ore Genesis: Implications to Mineral Exploration. Paulo Vasconcelos
Weathering, Landscape Evolution, Supergene Enrichment, and Ore Genesis: Implications to Mineral Exploration Paulo Vasconcelos 1 Lateritic/Oxide Au, Fe, Mn, Al, Ni, Cu, Zn, 2 Lateritic/Oxide Au, Fe, Mn,
More informationUser Guide to TopoToolbox - Plotting ksn-values
User Guide to TopoToolbox - Plotting ksn-values Background... 1 Load and generate required data... 2 Calculating the concavity index... 4 Calculate and plot ksn values... 4 Aggregate ksn values in stream
More informationMarine Science and Oceanography
Marine Science and Oceanography Marine geology- study of the ocean floor Physical oceanography- study of waves, currents, and tides Marine biology study of nature and distribution of marine organisms Chemical
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 informationSupplementary material
GSA DATA REPOSITORY2014088 Supplementary material Dominance of tectonics over climate in Himalayan denudation Vincent Godard, Didier L. Bourlès, Françoise Spinabella, Douglas W. Burbank, Bodo Bookhagen,
More informationPredictions of steady state and transient landscape morphology using sediment-flux-dependent river incision models
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 112,, doi:10.1029/2006jf000567, 2007 Predictions of steady state and transient landscape morphology using sediment-flux-dependent river
More informationKlaus Gessner, Chris Wijns, Louis Moresi, Fabio Boschetti and Alison Ord
Flow partitioning in the lithosphere during core complex formation: An interactive evolutionary computation approach using particle-in-cell finite elements Klaus Gessner, Chris Wijns, Louis Moresi, Fabio
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 informationLandscape Development
CHAPTER 22 Landscape Development Chapter Summary Landscapes are described in terms of their topography: elevation, the altitude of the surface of the Earth above sea level; relief, the difference between
More informationIntro to Geomorphology Key Concepts
Intro to Geomorphology Key Concepts Geomorphology Literally means the study of earth landforms - Geo = Earth - Morph=Form - Logos= Study of Involves understanding of - Mineralogy - Structure - Tectonics
More informationUNRAVELING THE HISTORY OF A LANDSCAPE: Using geomorphology, tephrochronology, and stratigraphy. Photo by: Josh Roering
UNRAVELING THE HISTORY OF A LANDSCAPE: Using geomorphology, tephrochronology, and stratigraphy Photo by: Josh Roering Photo: Eric Bilderback Photo by: Josh Roering Goal 1. Reconstruct the paleo-landscape
More informationImplications of the Saltation Abrasion Bedrock Incision Model for Steady-State River Longitudinal Profile Relief and Concavity
Earth Surface Processes and Landforms Steady-State Earth Surf. Process. Bedrock Landforms River Longitudinal 33, 1129 1151 Profile (2008) Relief and Concavity 1129 Published online in Wiley InterScience
More informationEdinburgh Research Explorer
Edinburgh Research Explorer GEOMORPHOLOGY Rivers split as mountains grow Citation for published version: Attal, M 2009, 'GEOMORPHOLOGY Rivers split as mountains grow' Nature Geoscience, vol. 2, no. 11,
More informationSediments and bedrock erosion
Eroding landscapes: fluvial processes Sediments and bedrock erosion Mikaël ATTAL Marsyandi valley, Himalayas, Nepal Acknowledgements: Jérôme Lavé, Peter van der Beek and other scientists from LGCA (Grenoble)
More informationJOURNAL OF GEOPHYSICAL RESEARCH, VOL. 117, F01004, doi: /2011jf002043, 2012
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 117,, doi:10.1029/2011jf002043, 2012 Late Neogene exhumation and relief development of the Aar and Aiguilles Rouges massifs (Swiss Alps) from low-temperature thermochronology
More information!" &#'(&) %*!+,*" -./0"1$ 1% % % - % 8 99:; < % % % % = 1. % % 2 /0 2 8 $ ' 99!; & %% % 2,A 1% %,1 % % % 2 3 %3 % / % / "1 % ; /0 % 2% % % %36
!" #$ &#'(&) *!+,*" - /0"1$ 1 1/0/// 0/02 /04"1 /0//,1$ 5/ - ( 6/027/ ///0 (/0// // - /002220(2 8 99:; < (/ = 1 2 /0$17899; 2 /0 2 8 $ 99?6 @ ' 99!; & 2,A 1,1 2 / / "1 -,14/02- ; /0 2 6,; B,1$ 2"1/0
More informationPlaattektoniek en Mickey Mouse: de bewegingen van de Aarde en de geologie van Marokko. G. Bertotti - TUDelft
Plaattektoniek en Mickey Mouse: de bewegingen van de Aarde en de geologie van Marokko G. Bertotti - TUDelft Moving continents Continent with matching boundaries Same fauna in different continents Similar
More informationInvestigating Controls on Bedrock River Incision Using Natural and Laboratory Experiments. Alexander C. Whittaker M.A. (cantab) MSci.
Investigating Controls on Bedrock River Incision Using Natural and Laboratory Experiments Alexander C. Whittaker M.A. (cantab) MSci. Thesis submitted for the degree of Doctor of Philosophy University of
More informationSUPPLEMENTARY INFORMATION
The major uncertainties in our model predictions arise from the input parameters, which include mantle density models (i.e. seismic tomography and choices about scaling velocities to temperature), crustal
More informationSurface process models and the links between tectonics and topography
University of Wollongong Research Online Faculty of Science, Medicine and Health - Papers Faculty of Science, Medicine and Health 2006 Surface process models and the links between tectonics and topography
More informationNumerical modeling of the Cenozoic geomorphic evolution of the southern Sierra Nevada, California
Earth and Planetary Science Letters 259 (2007) 85 96 www.elsevier.com/locate/epsl Numerical modeling of the Cenozoic geomorphic evolution of the southern Sierra Nevada, California Jon D. Pelletier Department
More informationCorrelating River Steepness with Erosion Rate along the Sri Lankan Escarpment. Matthew Potako
1 Correlating River Steepness with Erosion Rate along the Sri Lankan Escarpment Matthew Potako ABSTRACT A recent body of work has utilized cosmogenic isotopes as a measure of erosion rate in efforts to
More informationA Catchment Based Assessment Of The 3-Arc Second SRTM Digital Elevation Data Set
A Catchment Based Assessment Of The 3-Arc Second SRTM Digital Elevation Data Set a Martinez, C., a G.R. Hancock, a,b K.G. Evans, and b D.R. Moliere a School of Environmental and Life Sciences, The University
More informationCritical form and feedbacks in mountain-belt dynamics: the role of. rheology as a tectonic governor.
Critical form and feedbacks in mountain-belt dynamics: the role of rheology as a tectonic governor. Gerard H. Roe Department of Earth and Space Sciences, University of Washington, Seattle, WA. Mark T.
More informationSetting the Stage 8/28/09 1. Review. Landscape history of WY. Climate, geology and vegetation patterns in Wyoming
Setting the Stage Climate, geology and vegetation patterns in Wyoming 8/28/09 1 Review Vegetation ecology has its roots in phytogeography: Observation and description of PATTERNS of plant species distributions
More informationDiscovery of Crustal Movements areas in Lebanon: A study of Geographical information systems
Bulletin of Environment, Pharmacology and Life Sciences Bull. Env. Pharmacol. Life Sci., Vol 6[5] April 2017: 48-52 2017 Academy for Environment and Life Sciences, India Online ISSN 2277-1808 Journal s
More informationGeology 101 Lab 7: Isostasy Laboratory
Name TA Name Geology 101 Lab 7: Isostasy Laboratory To what thickness can sediments accumulate on the earth s surface, and what controls this thickness? How deeply can the continents be eroded? These are
More informationDistribution of active rock uplift along the eastern margin of the Tibetan Plateau: Inferences from bedrock channel longitudinal profiles
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 108, NO. B4, 2217, doi:10.1029/2001jb000861, 2003 Distribution of active rock uplift along the eastern margin of the Tibetan Plateau: Inferences from bedrock channel
More informationSubglacial Control on Glacier Flow in Northern Greenland
Subglacial Control on Glacier Flow in Northern Greenland Beáta Csathó (University at Buffalo, SUNY, Buffalo, NY), C.J. van der Veen (U. of Kansas, Lawrence, KS) Ralph van Frese and Tim Leftwich (The Ohio
More informationContinental Landscapes
Continental Landscapes Landscape influenced by tectonics, climate & differential weathering Most landforms developed within the last 2 million years System moves toward an equilibrium Continental Landscapes
More informationErosion Rate is a Function of Erodibility and Excess Shear Stress = k ( o - c ) From Relation between Shear Stress and Erosion We Calculate c and
Equilibrium, Shear Stress, Stream Power and Trends of Vertical Adjustment Andrew Simon USDA-ARS, Oxford, MS asimon@msa-oxford.ars.usda.gov Non-Cohesive versus Cohesive Materials Non-cohesive: sands and
More informationIsostasy, Bathymetry and the Physiography of the Ocean Floor
Isostasy, Bathymetry and the Physiography of the Ocean Floor EPSS 15 Spring 2017 Ad art for NBC Radio, 1939 Archimedes (c. 287 BCE 212 BCE) Greek mathematician, physicist and engineer Archimedes Principle:
More informationTectonics from topography: Procedures, promise, and pitfalls
Geological Society of merica Special Paper 398 26 Tectonics from topography: Procedures, promise, and pitfalls Cameron Wobus Kelin X. Whipple Department of Earth, tmospheric and Planetary Sciences, Massachusetts
More informationA generalized power law approximation for fluvial incision of bedrock channels
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 116,, doi:1.129/29jf1655, 211 A generalized power law approximation for fluvial incision of bedrock channels Nicole M. Gasparini 1 and Mark T. Brandon 2 Received 2
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 informationGeography of the world s oceans and major current systems. Lecture 2
Geography of the world s oceans and major current systems Lecture 2 WHY is the GEOMORPHOLOGY OF THE OCEAN FLOOR important? (in the context of Oceanography) WHY is the GEOMORPHOLOGY OF THE OCEAN FLOOR important?
More informationUNIT 3 GEOLOGY VOCABULARY FLASHCARDS THESE KEY VOCABULARY WORDS AND PHRASES APPEAR ON THE UNIT 3 CBA
UNIT 3 GEOLOGY VOCABULARY FLASHCARDS THESE KEY VOCABULARY WORDS AND PHRASES APPEAR ON THE UNIT 3 CBA A map that shows Earth s Topographic Map surface topography, which is Earth s shape and features Contour
More informationRemote Sensing and GIS Applications for Hilly Watersheds SUBASHISA DUTTA DEPARTMENT OF CIVIL ENGINEERING IIT GUWAHATI
Remote Sensing and GIS Applications for Hilly Watersheds SUBASHISA DUTTA DEPARTMENT OF CIVIL ENGINEERING IIT GUWAHATI Deciding Alternative Land Use Options in a Watershed Using GIS Source: Anita Prakash
More informationIntro to Quantitative Geology
Introduction to Quantitative Geology Lesson 13.1 Basic concepts of thermochronology Lecturer: David Whipp david.whipp@helsinki.fi 4.12.17 3 Goals of this lecture Introduce the basic concepts of thermochronology
More informationWhen one of things that you don t know is the number of things that you don t know
When one of things that you don t know is the number of things that you don t know Malcolm Sambridge Research School of Earth Sciences, Australian National University Canberra, Australia. Sambridge. M.,
More informationAN EVALUATION ON THE DATA QUALITY OF SRTM DEM AT THE ALPINE AND PLATEAU AREA, NORTH-WESTERN OF CHINA
AN EVALUATION ON THE DATA QUALITY OF SRTM DEM AT THE ALPINE AND PLATEAU AREA, NORTH-WESTERN OF CHINA Y. Liu School of Earth and Environmental Sciences, Lanzhou University, Lanzhou, 730000 -liuy@lzu.edu.cn
More information1. From raw topographic data to a river long-profile analysis. At the Matlab prompt, type. DEM = GRIDobj('srtm_bigtujunga30m_utm11.
Earth and Planetary Surface Processes Winter 2017 - Lab 4. River channel long profiles. Wieboldt 310C, 10:30a-11:20a Grades are not assigned for lab, but attendance is required. If you are unable to make
More informationFOSS LANDFORMS 8/1/2012. Content Professional Development Tools Contributor: Ben Fackler Adams. Assessing with Learning Progressions in Science
8/1/2012 Assessing with Learning Progressions in Science FOSS LANDFORMS Photo by Joanne Johnson Content Professional Development Tools Contributor: Ben Fackler Adams Learning Progression Landforms Prerequisite
More informationEdinburgh Research Explorer
Edinburgh Research Explorer A statistical framework to quantify spatial variation in channel gradients using the integral method of channel profile analysis Citation for published version: Mudd, SM, Attal,
More informationMaximum Extent of Pleistocene Glaciation - 1/3 of land surface Most recent glacial maximum peaked 18,000 years ago and is considered to have ended
Glaciers Maximum Extent of Pleistocene Glaciation - 1/3 of land surface Most recent glacial maximum peaked 18,000 years ago and is considered to have ended 10,000 B.P. Current Extent of Glaciation - about
More informationTerrain Units PALEOGEOGRAPHY: LANDFORM CREATION. Present Geology of NYS. Detailed Geologic Map of NYS
NYS TOPOGRAPHY Why so? PALEOGEOGRAPHY: LANDFORM CREATION Prof. Anthony Grande AFG 014 Present Geology of NYS Detailed Geologic Map of NYS Generalized Geology Detailed Geology Hot links to the fold out
More informationThe Main Points. The View from the Surface. Geology of Mars. Lecture #20: Reading:
Surface of Mars Lecture #20: Geology and Geologic Processes View from the Surface History/Evolution of the surface Reading: Chapter 9.4 The Main Points Mars has had a geologically active past that has
More informationLab 8: Gravity and Isostasy (35 points)
it's not the most important thing in your life right now. But what is important is gravity. Arnold Schwarzenegger as Colonel John Matrix, Commando (Check out this classic of American cinema!) Lab 8: Gravity
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 informationTopographic outcomes predicted by stream erosion models: Sensitivity analysis and intermodel comparison
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 07, NO. B9, 279, doi:0.029/200jb00062, 2002 Topographic outcomes predicted by stream erosion models: Sensitivity analysis and intermodel comparison G. E. Tucker School
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 informationMountain Rivers. Gutta cavat lapidem. (Dripping water hollows out a stone) -Ovid, Epistulae Ex Ponto, Book 3, no. 10, 1. 5
Mountain Rivers Gutta cavat lapidem (Dripping water hollows out a stone) -Ovid, Epistulae Ex Ponto, Book 3, no. 10, 1. 5 Mountain Rivers Fixed channel boundaries (bedrock banks and bed) High transport
More informationUsing Earthscope and B4 LiDAR data to analyze Southern California s active faults
Using Earthscope and B4 LiDAR data to analyze Southern California s active faults Exercise 8: Simple landscape morphometry and stream network delineation Introduction This exercise covers sample activities
More informationCover effect in bedrock abrasion: A new derivation and its implications for the modeling of bedrock channel morphology
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 112,, doi:10.1029/2006jf000697, 2007 Cover effect in bedrock abrasion: A new derivation and its implications for the modeling of bedrock channel morphology Jens M.
More informationPrashant P. Magar 1 and Nirupama P. Magar 2 1
Journal of Indian Geomorphology Volume 4, 2016 ISSN 2320-0731 Indian Institute of Geomorphologists (IGI) Application of Hack s Stream Gradient Index (SL Index) to Longitudinal Profiles of the Rivers Flowing
More informationFINAL EXAM Crustal Deformation CONVERGE DIVERGENT PLATES MANTLE PLUMES FLUX BASALTIC GRANITIC
Crustal Deformation Reading: Chapter 10 Pages 283-304 Review Questions 4, 6, 7, 10, 12, 15, 18, 20 FINAL EXAM NOON TO 2 PM, TUESDAY DEC. 5 HERE: Natural Science 101 BRING A SCAN TRON TURN IN YOUR REVIEW
More informationGIS feature extraction tools in diverse landscapes
CE 394K.3 GIS in Water Resources GIS feature extraction tools in diverse landscapes Final Project Anna G. Kladzyk M.S. Candidate, Expected 2015 Department of Environmental and Water Resources Engineering
More informationCOSMOGENIC 10 BE EROSION HISTORY OF THE BLUE RIDGE ESCARPMENT A LONG-LIVED FEATURE OF THE SOUTHERN APPALACHAINS
COSMOGENIC 10 BE EROSION HISTORY OF THE BLUE RIDGE ESCARPMENT A LONG-LIVED FEATURE OF THE SOUTHERN APPALACHAINS A Thesis Progress Report Presented by Colleen L. Sullivan to The Faculty of the Geology Department
More information44031 Crustal Compensation Rate In Central Anatolian Region And Comparison Of Seismic Activity
44031 Crustal Compensation Rate In Central Anatolian Region And Comparison Of Seismic Activity B. Oruç* (Kocaeli University), T. Sönmez (Kocaeli University), E. Ulutaş (Kocaeli University), H.H. Selim
More informationEarth Systems Science Chapter 7. Earth Systems Science Chapter 7 11/11/2010. Seismology: study of earthquakes and related phenomena
Earth Systems Science Chapter 7 I. Structure of the Earth II. Plate Tectonics The solid part of the earth system includes processes, just like the atmosphere and oceans. However, the time scales for processes
More informationMagnus-Rex and the new refraction Moho Map for southern Norway
Magnus-Rex and the new refraction Moho Map for southern Norway W. Stratford and H. Thybo Department of Geography and Geology, University of Copenhagen When did the present mountains form? Topography has
More informationControls on clastic systems in the Angoche basin, Mozambique: tectonics, contourites and petroleum systems
P2-2-13 Controls on clastic systems in the Angoche basin, Mozambique: tectonics, contourites and petroleum systems Eva Hollebeek, Olivia Osicki, Duplo Kornpihl Schlumberger, London, UK Introduction Offshore
More informationErosion Surface Water. moving, transporting, and depositing sediment.
+ Erosion Surface Water moving, transporting, and depositing sediment. + Surface Water 2 Water from rainfall can hit Earth s surface and do a number of things: Slowly soak into the ground: Infiltration
More information