Research Collection
Doctoral Thesis
Erosion and weathering of the Northern Apennines with implications for the tectonics and kinematics of the orogen
Author(s): Erlanger, Erica
Publication Date: 2020
Permanent Link: https://doi.org/10.3929/ethz-b-000393261
Rights / License: In Copyright - Non-Commercial Use Permitted
This page was generated automatically upon download from the ETH Zurich Research Collection. For more information please consult the Terms of use.
ETH Library Diss. ETH No. 26370
Erosion and weathering of the Northern Apennines with implications for the tectonics and kinematics of the orogen
Erica Danielle Erlanger Cover artwork by Reed Olsen DISS. ETH NO. 26370
Erosion and weathering of the Northern Apennines with implications for the tectonics and kinematics of the orogen
A thesis submitted to attain the degree of
DOCTOR OF SCIENCES of ETH ZURICH
(Dr. sc. ETH Zurich)
presented by
ERICA DANIELLE ERLANGER
Master of Science, Purdue University
born on 22.03.1986
citizen of France and the United States of America
accepted on the recommendation of
Prof. Dr. Sean D. Willett Prof. Vincenzo Picotti Prof. Sean F. Gallen Prof. Dr. Frank J. Pazzaglia
2020
2019 Abstract
Mountainous landscapes reflect the competition between denudation, uplift, and climate, which produce, modify, and destroy relief and topography. Bedrock rivers are dynamic topographic features and a critical link between these processes, as they record and convey changes in tectonics, climate, and sea level across the landscape. River incision models, such as the stream power model, are often used to quantify the relationship between topography and rock motion in the context of landscapes at steady state. At steady state, the stream power model predicts higher denudation rates for steeper river channels, while accounting for only the vertical motion of rock due to rock uplift or denudation. However, natural landscapes often have more complicated histories, particularly in convergent orogens with asymmetric topography, where steady state requires that denudation must balance both vertical and horizontal rock motion.
This thesis addresses this central issue by comparing the spatial and temporal pattern of denudation with metrics of topographic steepness in the Northern Apennine Mountains of Italy, a young and active orogen with asymmetric topography. New and existing catchment-averaged denudation rates from cosmogenic 10Be concentrations demonstrate that the steeper flank of the Northern Apennines is eroding more slowly than the gentler flank. Long-term denudation rates inverted from low-temperature thermochronometers show that this pattern of denudation across the orogen is long-lived, since at the least 3—5 Ma, and that denudation rates have decreased on the Ligurian side through time. The apparent decoupling between denudation rates and topography is resolved with a kinematic model of the orogenic wedge that accounts for the full vertical and horizontal rock velocity field. This model reconciles the 10Be concentrations, geomorphic observations, and geodetic rates of rock motion with the topography of the Northern Apennines, and provides new estimates for slab retreat rates consistent with recent estimates from tomography, surface geology, and morphology.
This thesis also explores the partitioning of denudation into physical erosion and chemical weathering in the Northern Apennines. Chemical weathering in particular is an important control on landscape evolution and the global CO2 budget. Most studies have focused on weathering in orogens comprised of silicate-rich lithologies, which can remove CO2 from the atmosphere over geologic timescales, whereas carbonate weathering is generally considered to be CO2 neutral. However, even in silicate-rich landscapes, carbonate weathering dominates total solute fluxes. Recently uplifted orogens in particular are often characterized by carbonate-rich, marine sedimentary sequences, so the global weathering flux of carbon and calcium to the oceans should be more strongly influenced by these orogens. However, the partitioning of denudation fluxes remains largely unexplored in mixed lithology orogens, so, it is unclear whether the same processes that control erosion and weathering apply to both silicate-rich and mixed lithology settings. Here, denudation fluxes from the Northern Apennines are partitioned into carbonate and silicate chemical weathering and physical erosion fluxes. These fluxes demonstrate that denudation is dominated by physical erosion of both silicate and carbonate rocks; carbonate physical erosion is controlled by lithology; weathering fluxes are dominated by carbonate dissolution; and denudation is negatively correlated with runoff. Finally, denudation fluxes from the Northern Apennines are similar to other temperature mountain ranges (e.g. Southern Alps of New Zealand), although total weathering fluxes from this study are generally higher, due to greater carbonate weathering fluxes. The results from this thesis challenge current interpretations regarding denudation rates through space and time and contribute to a broader understanding of surface and crustal processes in the Northern Apennines.
2 Zusammenfassung
Gebirgslandschaften spiegeln den Wettstreit zwischen Denudation, Hebung und Klimaeinflüssen wider, welche Relief und Topographie hervorbringen, verändern und zerstören. Erosive Flüsse sind dynamische topographische Merkmale dieser Landschaften und stellen eine wichtige Verbindung zwischen den genannten Prozessen dar, da sie Veränderungen der Tektonik, des Klimas und des Meeresspiegels aufzeichnen und diese Veränderungen auf die Landschaft übertragen. Flusserosionsmodelle wie das Stream-Power-Modell werden häufig verwendet, um den Zusammenhang zwischen Topographie und Hebungsrate im Gleichgewicht zu quantifizieren. Im Gleichgewicht sagt das Stream-Power-Modell höhere Denudationsraten für steile Flusslängsprofile voraus, wobei aber nur die vertikale Bewegung des Gesteins wegen Hebung oder Denudation berücksichtigt wird. Jedoch weisen natürliche Landschaften meist eine kompliziertere Entwicklung auf, insbesondere bei konvergenten Orogenen mit asymmetrischer Topographie, bei denen der Gleichgewichtszustand erfordert, dass die Denudation sowohl die vertikale als auch die horizontale Gesteinsbewegung ausgleicht.
Diese Dissertation befasst sich mit diesem zentralen Problem durch den Vergleich des räumlichen und zeitlichen Denudationsmusters mit der Steilheit der Topographie im nördlichen Apennin Italiens, einer jungen und aktiven Gebirgskette mit asymmetrischer Topographie. Schätzungen der Denudationsraten anhand von neuen und bestehenden, über das Einzugsgebiet gemittelten Konzentrationsmessungen von kosmogenem 10Be in Flusssedimenten zeigen, dass die steilere Flanke des nördlichen Apennins langsamer erodiert als die flachere Flanke. Langfristige Denudationsraten bestimmt durch Niedertemperatur-Thermochronometrie zeigen, dass dieses Denudationsmuster seit mindestens 3-5 Ma besteht und dass die Denudationsraten auf der ligurischen Seite im Laufe der Zeit abgenommen haben. Die anscheinende Entkoppelung zwischen Denudationsraten und Topographie kann mithilfe eines kinematischen Modells des Gebirgskeils, dass das vollständige vertikale und horizontale Gesteinsgeschwindigkeitsfeld berücksichtigt, erklärt werden. Dieses Modell stimmt die 10Be-Konzentrationen, geomorphologischen Beobachtungen und geodätischen Geschwindigkeiten der Gesteinsbewegung mit der Topographie des nördlichen Apennins ab und liefert neue Schätzungen der Rückzugsraten der Lithosphärenzunge im Einklang mit Schätzungen basierend auf Tomographie, Oberflächengeologie und Morphologie.
In dieser Dissertation wird außerdem die Aufteilung der Denudation in chemische Verwitterung und physikalische Erosion im nördlichen Apennin untersucht. Insbesondere die chemische Verwitterung spielt eine wichtige Rolle für die Landschaftsentwicklung und das globale CO2-Budget. Bisherige Forschungsarbeiten konzentrieren sich hauptsächlich auf Verwitterung in Orogenen bestehend aus silikatreichen Gesteinen, da diese der Atmosphäre
CO2 entziehen. Karbonat-Verwitterung hingegen ist CO2-neutral. Allerdings dominiert auch in Silikat-reichen Landschaften die Karbonat-Verwitterung den Gesamtfluss gelöster Stoffe. Vor allem kürzlich emporgehobene Orogene sind häufig durch karbonatreiche, marine Sedimentsequenzen gekennzeichnet, weshalb der globale Verwitterungsfluss von Kohlenstoff und Kalzium in die Ozeane stärker von solchen Orogenen beeinflusst werden sollte. Trotzdem ist die Aufteilung der Denudation in Orogenen mit gemischten Lithologien noch weitgehend unerforscht. Daher ist unklar, ob die gleichen Prozesse, die Erosion und Verwitterung steuern, sowohl auf silikatreiche als auch auf gemischte Lithologien zutreffen. Hier werden die Denudationsflüsse des nördlichen Apennins in Karbonat- und Silikatverwitterungflüsse und in physikalische Erosionsflüsse aufgeteilt. Diese Flüsse zeigen, dass die Denudation negativ mit dem Abfluss korreliert, die Denudation von physikalischer Erosion sowohl von Silikat- als auch von Karbonatgesteinen dominiert wird, die physikalische Erosion des
3 Karbonatgesteins durch die Lithologie beherrscht wird und die Verwitterungsmengen von der Lösung von Karbonatgestein dominiert werden. Schließlich sind die Denudationsraten im nördlichen Apennin vergleichbar mit denen in anderen Gebirgsketten in gemäßigten Klimazonen (z.B. Neuseeländische Alpen), obwohl die Verwitterungsraten dieser Studie im Vergleich höher sind. Die Ergebnisse dieser Dissertation stellen aktuelle Interpretationen der Denudationsraten durch Raum und Zeit in Frage und tragen zu einem breiteren Verständnis der Oberflächen- und Krustenprozesse im nördlichen Apennin bei.
4 Acknowledgements
I am indebted to the many people who have contributed to and shaped my research at ETH and life in Zürich. Pursuing a Doctoral Degree at the ETH Zürich has been a privilege, and my time in Switzerland has been a wonderful chapter in my life, filled with many new faces who I hope will continue to play a role in my future endeavors and adventures.
I must first thank my advisor, Sean Willett, or “Big Sean”. Thank you for taking me on as your student, for challenging me and supporting me in my PhD research, and for helping me to obtain the postdoc that awaits me in Potsdam. You are also a man of many other interests and talents, so I must also thank you for the occasional plant care tips and for the few times we went rappelling down your waterfalls in Bonassola. I would also like to thank the rest of my doctoral committee. Sean Gallen, your teaching style and enthusiasm for geology have been an inspiration to me during my PhD, and I’m so thankful for your support and willingness to brainstorm and execute new ideas. Thank you for your contributions to my research, and I hope we’ll continue to work together in the future. Vincenzo, you’ve taught me so much about the Apennines and about Italian food, I’m not even sure which is more important at this point. I know that to find the strath of a river terrace, you should follow the trail of mint plants or water that seeps from the terrace. I can also safely say that I know a Minestrone Genovese is only ready when the spoon sticks straight out of the soup, and that the best brodo is made with an old chicken. Your thirst for geology and Italian food, in all their glorious detail, have made my time at ETH both enjoyable and delicious. I would also like to express my gratitude to Frank Pazzaglia for flying out to Switzerland and acting as the external member of my examining committee. Finally, my thanks go to Whitney Behr, who agreed to act as the chairperson for my defense. I must also express my sincere gratitude to Giuditta. It’s hard for me to describe how important you’ve been to me and my research at ETH. You trained me to prep samples, took the time to date my samples, and worked closely with me to make sure I understood the ins and outs of thermochronometry. Thank you for your commitment to my work, and for occasional pizza dates.
My collaborators have also helped me in no small measure. My gratitude goes to Aaron Bufe, without whom my weathering chapter and SNF proposal would’ve suffered from a fundamental knowledge gap. Thank you for your hundreds of poignant, intelligent, and insightful comments on my writing. I am thankful for your help, and I am honored that I will be able to work with you for the next two years at the GFZ. I would also like to thank Jeremy Rugenstein. You’ve played a crucial role in fostering my interest in weathering and geochemistry. Thank you for your help in the field, for measuring my water chemistry samples, for always having an open ear to brainstorm, and for your eternally positive attitude. Thank you to Maarten Lupker for your advice on all things cosmogenic nuclides, and to Negar Hagipour, for all your help with processing and measuring my difficult 10Be samples. Thank you to the group of Frederic Hermann for hosting me at UNIL and allowing me to use their OSL labs and equipment.
I would also like to extends my thanks to my colleagues and friends for the professional and personal support they have given me. To the ESD Group: Thank you to former and current members of the ESD group: Richard Ott, Yanyan Wang, Odin Marc, Chia-Yu Chen, Malwina San Jose, Loraine Gourbet, and Kosuke Ueda for stimulating discussions, feedback, hikes, and trips to the Limmat. I must also thank our group secretary, Katharin Fehr. I can’t thank you enough for your help with the admissions process and with all the logistics of moving to another continent, and finally for the wonderful John baker fruit breads you brought to our institute coffee. Thank you to Alina Fiedrich and Richard Ott for editing my terrible German writing and to Fabian 5 Kuhn for moral support. A special thanks goes to Sascha Winterberg: my office mate, Alpine expert, and friend. Thank you for the many hikes, trips, vacations, family holidays, and conversations about the Greater Alpine area that we’ve shared. I hope we stay in each other’s lives for a long time. Special thanks also go to Julia Krawielicki. We’ve had so many wonderful moments together, and I don’t know how I’ll manage to see action movies without you. Thank you for the countless times you helped to pin me into my latest fabric creation and for also including me in your family Christmas celebrations.
Finally, I must thank my family. To my mom and dad, Jackie and Robert, thank you for being understanding of my choice to live and study on another continent, for traveling to see me, and editing my scientific writing. Thank you to my Swiss family: Michael, Marianne, Simon, Annica, Max, Lucian, Andrea, and Coco. I met all of you only four years ago when I arrived in Switzerland, but you quickly felt like family. Thank you for allowing me to be a part of your lives, for sharing brunches with me, and making Switzerland feel more like a home.
6 Chapter 1 Introduction
7 1 Introduction
Active orogens have the highest rates of uplift and denudation on Earth, and these processes compete to build, modify, and destroy mountain topography (Chamberlin, 1899; Raymo et al., 1988; Molnar and England, 1990; Gaillardet et al., 1999; Willett et al., 2001; West et al., 2005). Classical, conceptual models of landscape evolution have long recognized that the growth of topography is driven by tectonic uplift, while the degradation and lowering of topography are driven by denudation processes (Davis, 1899; Penck, 1953; Hack, 1975).
Rivers provide a fundamental link between denudation processes, uplift, and the evolution of topography. In particular, bedrock rivers are archives of changes in tectonics, climate, and sea level and transmit these changes over the landscape (Whipple and Tucker, 1999). For example, preserved fluvial strath terraces can be dated to estimate paleo-erosion rates or the rate of river incision, a proxy for rock uplift (Pazzaglia and Brandon, 2001; Picotti and Pazzaglia, 2008; Cyr and Granger, 2008; Erlanger et al., 2012). As rivers are the major transporting agents of sediment and solutes from the continent to the ocean, and they integrate large portions of the land surface, we can use the wealth of information from rivers to estimate denudation processes on a regional scale.
There is an important distinction between denudation, which describes the overall removal of rock from the Earth’s surface, and its components: physical erosion and chemical weathering. Physical erosion is the mechanical breakdown of rock, and chemical weathering is the dissolution of rocks at the surface. Each of these processes can be measured in the landscape using various methods that integrate over different spatial and temporal scales. Physical erosion has historically been quantified with sediment yield measurements from river stream gauges and reflect short-term rates (1-10 years) (Clayton and Megahan, 1986; Granger and Riebe, 2013), and chemical weathering has been estimated from water solute fluxes (Gaillardet et al., 1997; Jacobson et al., 2003; West et al., 2005; Gaillardet et al., 2018). During the last couple of decades, cosmogenic nuclide dating has emerged as a standard way of measuring denudation, chemical weathering, and physical erosion in landscapes (Granger and Riebe, 2013). Chemical weathering and physical erosion rates can be quantified together from the mass balance of elements in soils and bedrock, (Riebe et al., 2001b, 2004; Granger and Riebe, 2013), and denudation can be quantified from bedrock landforms (e.g. inselbergs) (Bierman and Caffee, 2002), unconsolidated soils (Riebe et al., 2001a), or from alluvial sediment and watersheds (Granger et al., 1996, 1997). Estimating geologically recent (102 – 105 yrs), catchment-wide denudation rates with cosmogenic 10Be from bulk river sediment has become a well-established method in geomorphology (Granger et al., 1996; Bierman and Steig, 1996; von Blanckenburg, 2005). The concentration of cosmogenic 10Be reflects the residence time of a quartz mineral near the surface, as it is primarily produced within the upper meter below the ground surface. As rivers channels integrate hillslope and fluvial processes, they are a natural mixing agent for sediment, so that the 10Be concentration in river sediment reflects the spatially-averaged denudation rate of the entire upstream area that contributes sediment.
Over long timescales, orogenic systems are suggested to approach an equilibrium between rates of denudation and rock uplift, known as “steady-state”. At steady state, empirical scaling relationships between denudation and river channel steepness predict that bedrock river channels are steeper in more rapidly eroding landscapes. These relationships can be expressed by the stream power model, which describes the denudation rate in terms of area (A), slope (S), rock erodibility (K), and nondimensional exponents (m and n):
E = K⋅Am⋅Sn (1) 8 where rock erodibility includes substrate, hydrology, and climate processes. (Wobus et al., 2006; Kirby and Whipple, 2012). We can then solve for this equation in terms of slope:
S = (E/K)(1/n) ⋅A(m/n) (2)
Empirical observations of river profile geometries (Hack, 1957; Flint, 1974) first described slope interms of the power-law relationship between local channel slope (ks), contributing drainage area (A), and channel concavity (-θ):
-θ S = ks⋅A (3)
Combining equations (2) and (3), we can link the observed topographic form of the river channel with the stream power incision model, and solve for denudation rate in terms of channel steepness normalized for drainage area (ksn), rock erodibility (K), and the nondimensional exponent (n):
n E = K⋅ksn (4)
Most emperical studies have observed monotonic increases in denudation rates with respect to channel steepness, consistent with the function form of equation 1 (Safran et al., 2005; Harkins et al., 2007; DiBiase et al., 2010; Miller et al., 2013). In most of these cases, the relationship is non-linear and consistent with a power law exponent greater than 1 (Ouimet et al., 2009; Cyr et al., 2014; Hilley and Young, 2018) The Apennines Mountains in particular host some of the highest measured denudation rates, which are associated with some of the lowest channel steepness values (Figure 1). In this study, I explore relationships between denudation, weathering, erosion, and topography in the Apennines Mountains, in order to understand the dominant controls (e.g. tectonics, climate, lithology) on landscape evolution.
9 A) 800 E. Tibet, China (Ouimet et al., 2009) Appalachians, USA (Miller et al., 2013) Andes, Bolivia (Safran et al., 2005) Apennine Mtns, Italy (Cyr et al., 2010) NE Tibet, China (Hawkins et al., 2007) 600 San Gabriel Mtns, USA (DiBiase et al., 2010) E and S Alps, Europe (Norton et al., 2011) Rwenzori Mtns, E Africa (Roller et al., 2012) )
0.9 400 (m sn K
200
0 0 500 1000 1500 2000 2500 10Be Denudation Rate (m/Myr)
B) 1000
100 ) 0.9 (m sn K
E. Tibet, China (Ouimet et al., 2009) 10 Appalachians, USA (Miller et al., 2013) Andes, Bolivia (Safran et al., 2005) Apennine Mtns, Italy (Cyr et al., 2010) NE Tibet, China (Hawkins et al., 2007) San Gabriel Mtns, USA (DiBiase et al., 2010) E and S Alps, Europe (Norton et al., 2011) Rwenzori Mtns, E Africa (Roller et al., 2012) 1 1 10 100 1000 10000 10Be Denudation Rate (m/Myr) Figure 1. Comparison of basin-averaged denudation rate (m/Ma) with normalized channel steepness for orogens around the world on A) linear axes and B) log-log scale axes. Original data sources are: Bolivian Andes (Safran et al., 2005), Northeastern Tibet (Harkins et al., 2007), Eastern Tibet (Ouimet et al., 2009), San Gabriel Mountains (DiBiase et al., 2010), Apennines Mountains (Cyr et al., 2010), Eastern and Southern Alps (Norton et al., 2011), Rwenzori Mountains (Roller et al., 2012), Appalachian Mountains (Miller et al., 2013), and Santa Lucia Mountains (Hilley and Young, 2018).
10 1.1 Northern Apennines
One of the best-studied, active mountain chains in the world is the Apennine Mountains, which form the backbone of peninsular Italy (Figure 2). The Apennines are a young orogen, broadly formed due to the south- westward subduction of the Adriatic microplate beneath Eurasia since 30 Ma (Ricci Lucchi, 1986; Dewey et al., 1989). The Apennines developed as a subaqueous accretionary wedge, fed from the uplifting Central Alps until the Miocene, (Ricci Lucchi, 1975; Roveri et al., 2001; Gandolfi et al., 2007). From the Miocene to Plio- Pleistocene, rapid exhumation resulted in the sub-aerial exposure of the Apenninic wedge and the creation of topographic relief (Balestrieri et al., 1996; Boccaletti and Sani, 1998; Abbate et al., 1999; Fellin et al., 2007).
Unlike the classical model for convergent orogens, the Apennines evolved in response to rollback of the lower plate. Accretion of material into the upper plate coupled with retreat of the subduction system has produced contemporaneous crustal shortening in the external part of the orogen, and extension in the internal part of the orogen (Vai and Martini, 2001). Today, the Northern Apennines are characterized by relatively low average elevations of 400 m, mixed siliciclastic and carbonate lithologies, and a primary drainage divide that separates rivers draining to the Adriatic Sea (Adriatic side) from rivers draining to the Ligurian Sea (Ligurian side) (Figure 3).
ALPS N
Po Plain DINARIDES Adriatic Sea
Ligurian Sea APENNINES
Rome Ligurian Provençal Basin Tyrrhenian Sea
0 200 Km Figure 2. Colored hillshade map of Italy and ocean bathymetry, showing locations of mountain ranges in the greater Alpine area and adjacent ocean basins. Location of the study area is given by the white box. 11 9 E 10 E 11 E 12 E N 5 4
Parma
Bologna
Genoa
La Spezia N 4 4 Ligurian Sea 0 20 40 60 80 m Florence Figure 3. Topography of study area outlined in white on Figure 2. Black line illustrates the drainage divide that separates rivers draining to the Ligurian Sea and rivers that drain to the Adriatic Sea.
The topography of the Northern Apennines is asymmetric: at the 10-km scale, the Ligurian side is shorter and steeper, whereas the Adriatic side is longer and has gentler slopes. Numerical models suggest that horizontal rock motion is responsible for asymmetric topography in orogens, such as the Southern Alps of New Zealand, the Olympic Mountains of the USA, and the Central Range of Taiwan (Willett et al., 2001). These model results are consistent with the distribution of exhumation ages across the orogen (Fuller et al., 2006; Willett and Brandon, 2002; Willett et al., 2003), which reflects the unroofing history of rock as it approaches the surface, and is attributable to tectonic or denudational processes. Exhumation ages are determined with radioisotopic thermochronometers, which reflect the cooling history of rocks over timescales of 105-107 years and depths between 2-10 km (Reiners and Brandon, 2006). Using multiple thermochronometers within the same location can also constrain material paths in an orogen, and the relative importance of horizontal versus vertical rock motions. In the Northern Apennines, the distribution of exhumation ages from multiple thermochronometers suggests that material paths have an important component of horizontal rock motion (Thomson et al., 2010). Regionally, exhumation ages are youngest on the Adriatic side relative to the Ligurian side, which translates to an order of magnitude higher denudation rates on the Adriatic side between 3-5 Ma (Thomson et al., 2010). One of the main questions that remains is whether the evolution of the Northern Apennines and its asymmetric topography are also reflected in modern denudation rates?
Currently, modern denudation estimates exist only for the Adriatic side, and rates for major rivers vary between ~0.2–0.8 mm/yr (Cyr and Granger, 2008; Cyr et al., 2014; Wittmann et al., 2016). Previous studies comparing 10Be denudation rates, river incision rates, and paleo-erosion rates suggest that the Northern Apennines have been in a state of dynamic equilibrium over the past 900 ka (Cyr and Granger, 2008). This state of equilibrium suggests that the steeper Ligurian side should be eroding faster than the gentler Adriatic side, based on equation
1. However, low denudation rates coupled with high ksn values (Figure 1) have been observed along one Adriatic River, and this pattern was attributed to tributaries draining different lithologies with variable erodibility (K) values (Cyr et al., 2014). This study and others address a current debate in geomorphology on the role of lithology in modulating channel steepness. Some find a strong lithologic control on steepness from modeling (Stock and Montgomery, 1999) and field studies (Gallen, 2018), whereas others have found that denudation
12 rates and ksn are insensitive to rock type (Miller et al., 2013; Hilley and Young, 2018).
Although the role of lithology on denudation is debated, early studies on chemical weathering have shown that lithology is a primary control, as it dictates the minerals available at the surface for weathering (Stallard and Edmond, 1983; Meybeck, 1987; Galy and France-Lanord, 1999). Young orogens, such as the Northern Apennines, are typically characterized by marine sedimentary sequences that contain important volumes of carbonate, which can enhance chemical weathering, as carbonate weathers a factor of 3 times faster than silicates (Meybeck, 1987; Gaillardet et al., 2013). However, previous studies have estimated chemical weathering rates in orogens containing large volumes of silicate-rich minerals such as the Andes Mountains (Gaillardet et al., 1997), Sierra Nevada Mountains of California (Riebe et al., 2001b), Himalaya Mountains (West et al., 2005), and New Zealand Southern Alps (Jacobson et al., 2003; Jacobson and Blum, 2003), so chemical weathering rates are lacking from landscapes with more typical, mixed lithologic assemblages. As such, the Northern Apennine Mountains present an excellent opportunity to resolve denudation, physical erosion, and chemical weathering budgets in a young orogen.
Over the last several decades, numerous studies have addressed rates of exhumation (Balestrieri et al., 1996; Ventura et al., 2001; Zattin et al., 2002; Balestrieri et al., 2003; Fellin et al., 2007; Thomson et al., 2010; Malusà and Balestrieri, 2012; Carlini et al., 2013; Balestrieri et al., 2018), paleo-denudation and modern denudation rates on the Adriatic side of the Northern Apennines (Cyr and Granger, 2008; Cyr et al., 2014; Wittmann et al., 2016), geodetic uplift rates (D’Anastasio et al., 2006; Ferranti et al., 2006; Serpelloni et al., 2013), horizontal GPS rates (Bennett et al., 2012; Devoti et al., 2017) and orogen kinematics (Thomson et al., 2010). Despite numerous studies on the evolution of the Northern Apennines, and studies on the rates of denudation and surface deformation through time and space, a number of knowledge gaps still exist. (1) We have no estimates of catchment-averaged denudation rates or long-term denudation rates from rivers on the Ligurian side of the Northern Apennines to compare with existing estimates from the Adriatic side. (2) We lack information on chemical weathering and erosional fluxes for young orogens comprised of mixed siliciclastic-carbonate lithologies.
Another important point to consider is the fundamental assumption of a vertical reference frame for denudation rates and ksn estimates, which has implications for conditions of steady state. Steady state is commonly assumed to represent a balance between the vertical motions associated with denudation and rock uplift (England and Molnar, 1990), although previous modeling and field-based studies on active orogens have indicated that steady state, in fact, reflects a balance between both the horizontal and vertical motion of rock (Willett et al., 2001; Miller et al., 2007). Catchment-averaged denudation rates and ksn estimates also assume a vertical reference frame, whereas the pattern of exhumation rates in the Northern Apennines has shown an important horizontal component of rock motion (Thomson et al., 2010). Hence, one of the reasons the relationship between erodibility, denudation, and channel steepness in the Northern Apennines shows different patterns from other active orogens (Figure 1) could also be due to incorrect assumptions about the nature of denudation rates. As such, understanding the evolution and present topography of the Northern Apennines requires a holistic approach to quantifying and assessing the relationships between surface and crustal-scale processes. We address these knowledge gaps by taking advantage of the wealth of information provided by the geometry of river catchments, longitudinal river profiles, riverine water chemistry, and river sediments. The implications of this research are broad and will contribute new, quantitative data on the relationship between denudation rates through space and time, and the tectonics and kinematics of the Northern Apennines orogen. 13 1.2 Thesis Objectives The goal of this research is broadly to understand the spatial and temporal variability in surface and crustal- scale processes in the Northern Apennines. This thesis is divided into three chapters. Chapter 2 focuses on interpreting the kinematics of the Northern Apennines orogen, through a compilation of new and existing catchment-wide denudation rates derived from concentrations of cosmogenic 10Be (Cyr and Granger, 2008; Cyr et al., 2014; Wittmann et al., 2016) and geodetic rates of vertical and horizontal rock motion (D’Anastasio et al., 2006; Bennett et al., 2012; Devoti et al., 2017). Chapter 3 reviews existing exhumation data from low-temperature thermochronometers, contributes new detrital apatite fission-track ages, models long-term denudation rates, and proposes a new model for the spatial and temporal development of denudation through time across the Northern Apennines orogenic wedge. Chapter 4 focuses on calculating denudation, weathering, and erosional fluxes in the Northern Apennines from10 Be concentrations, water solutes, and the carbonate sand in the catchment. These fluxes are then partitioned between silicate and carbonate lithologies and compared with results from silicate-rich endmembers locations in other orogens. Finally, Chapter 5 provides a synthesis and concluding remarks on this thesis, as well as future work directions.
References Abbate, E., Balestrieri, M.L., Bigazzi, G., Ventura, B., Zattin, M., and Zuffa, G.G., 1999, An extensive apatite fission- track study throughout the Northern Apennines Nappe belt: Radiation Measurements, v. 31, p. 673–676, doi: Doi 10.1016/S1350-4487(99)00168-7. Balestrieri, M.L., Abbate, E., and Bigazzi, G., 1996, Insights on the thermal evolution of the Ligurian Apennines (Italy) through fission-track analysis: Journal of the Geological Society, v. 153, p. 419–425, doi: 10.1144/gsjgs.153.3.0419. Balestrieri, M.L., Benvenuti, M., and Catanzariti, R., 2018, Unravelling basin shoulder dynamics through detrital apatite fission-track signature: the case of the Quaternary Mugello Basin, Italy: Geological Magazine, v. 155, p. 1413– 1426, doi: 10.1017/s0016756817000073. Balestrieri, M.L., Bernet, M., Brandon, M.T., Picotti, V., Reiners, P.W., and Zattin, M., 2003, Pliocene and Pleistocene exhumation and uplift of two key areas of the Northern Apennines: Quaternary International, v. 101, p. 67– 73, http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&cmd=Retrieve&dopt=AbstractPlus&list_ uids=13036896733521009365related:1Vr4toFc7LQJ. Bennett, R.A., Serpelloni, E., Hreinsdóttir, S., Brandon, M.T., Buble, G., Basic, T., Casale, G., Cavaliere, A., Anzidei, M., Marjonovic, M., Minelli, G., Molli, G., and Montanari, A., 2012, Syn-convergent extension observed using the RETREAT GPS network, northern Apennines, Italy: Journal of Geophysical Research: Solid Earth, v. 117, doi: 10.1029/2011JB008744. Bierman, P.R., and Caffee, M., 2002, Cosmogenic exposure and erosion history of Australian bedrock landforms: Bulletin of the Geological Society of America, v. 114, p. 787–803, doi: 10.1130/0016-7606(2002)114<0787:CEAEHO>2.0.CO;2. Bierman, P.R., and Steig, E.J., 1996, Estimating rates of denudation using cosmogenic isotope abundances in sediment: Earth Surface Processes and Landforms, v. 21, p. 125–139, doi: 10.1002/(SICI)1096-9837(199602)21:2<125::AID- ESP511>3.0.CO;2-8. von Blanckenburg, F., 2005, The control mechanisms of erosion and weathering at basin scale from cosmogenic nuclides in river sediment: Earth and Planetary Science Letters, v. 237, p. 462–479, doi: 10.1016/j.epsl.2005.06.030. Boccaletti, M., and Sani, F., 1998, Cover thrust reactivations related to internal basement involvement during Neogene- Quaternary evolution of the northern Apennines: Tectonics, v. 17, p. 112–130, doi: 10.1029/97TC02067. Carlini, M., Artoni, A., Aldega, L., Balestrieri, M.L., Corrado, S., Vescovi, P., Bernini, M., and Torelli, L., 2013, Exhumation and reshaping of far-travelled/allochthonous tectonic units in mountain belts. New insights for the relationships between shortening and coeval extension in the western Northern Apennines (Italy): Tectonophysics, v. 608, p. 267–287, doi: 10.1016/j.tecto.2013.09.029. Chamberlin, T.C., 1899, An attempt to frame a working hypothesis of the cause of glacial periods on an atmospheric basis: The Journal of Geology, v. 7, p. 545–584. Clayton, J.L., and Megahan, W.F., 1986, Erosional and chemical denudation rates in the Southwestern Idaho batholith: Earth Surface Processes and Landforms, v. 11, p. 389–400, doi: 10.1002/esp.3290110405. Cyr, A.J., and Granger, D.E., 2008, Dynamic equilibrium among erosion, river incision, and coastal uplift in the northern 14 and central Apennines, Italy: Geology, v. 36, p. 103–106, doi: 10.1130/G24003A.1. Cyr, A.J., Granger, D.E., Olivetti, V., and Molin, P., 2014, Distinguishing between tectonic and lithologic controls on bedrock channel longitudinal profiles using cosmogenic10Be erosion rates and channel steepness index: Geomorphology, v. 209, p. 27–38, doi: 10.1016/j.geomorph.2013.12.010. Cyr, A.J., Granger, D.E., Olivetti, V., and Molin, P., 2010, Quantifying rock uplift rates using channel steepness and cosmogenic nuclide-determined erosion rates: Examples from northern and southern Italy: Lithosphere, v. 2, p. 188–198, doi: 10.1130/L96.1. D’Anastasio, E., De Martini, P.M., Selvaggi, G., Pantosti, D., Marchioni, A., and Maseroli, R., 2006, Short-term vertical velocity field in the Apennines (Italy) revealed by geodetic levelling data: Tectonophysics, v. 418, p. 219–234, doi: 10.1016/j.tecto.2006.02.008. Davis, W.M., 1899, The geographical cycle: The Geographical Journal, v. 14, p. 481–504. Devoti, R., D’Agostino, N., Serpelloni, E., Pietrantonio, G., Riguzzi, F., Avallone, A., Cavaliere, A., Cheloni, D., Cecere, G., D’Ambrosio, C., Falco, L., Selvaggi, G., Métois, M., Esposito, A., et al., 2017, A combined velocity field of the mediterranean region: Annals of Geophysics, v. 60, doi: 10.4401/ag-7059. Dewey, J.F., Helman, M.L., Knott, S.D., Turco, E., and Hutton, D.H.W., 1989, Kinematics of the western Mediterranean: Geological Society, London, Special Publications, v. 45, p. 265–283, doi: 10.1144/GSL.SP.1989.045.01.15. DiBiase, R.A., Whipple, K.X., Heimsath, A.M., and Ouimet, W.B., 2010, Landscape form and millennial erosion rates in the San Gabriel Mountains, CA: Earth and Planetary Science Letters, v. 289, p. 134–144, doi: 10.1016/j. epsl.2009.10.036. England, P., and Molnar, P., 1990, Surface uplift, uplift of rocks, and exhumation of rocks: Geology, v. 18, p. 1173–1177. Erlanger, E.D., Granger, D.E., and Gibbon, R.J., 2012, Rock uplift rates in South Africa from isochron burial dating of fl uvial and marine terraces: Geology, v. 40, p. 1019–1022, doi: 10.1130/G33172.1. Fellin, M.G., Reiners, P.W., Brandon, M.T., Wüthrich, E., Balestrieri, M.L., and Molli, G., 2007, Thermochronologic evidence for the exhumational history of the Alpi Apuane metamorphic core complex, northern Apennines, Italy: Tectonics, v. 26, doi: 10.1029/2006TC002085. Ferranti, L., Antonioli, F., Mauz, B., Amorosi, A., Dai Pra, G., Mastronuzzi, G., Monaco, C., Orrù, P., Pappalardo, M., Radtke, U., Renda, P., Romano, P., Sansò, P., and Verrubbi, V., 2006, Markers of the last interglacial sea-level high stand along the coast of Italy: Tectonic implications: Quaternary International, v. 145–146, p. 30–54, doi: 10.1016/j. quaint.2005.07.009. Flint, J.J., 1974, Stream gradient as a function of order, magnitude, and discharge: Water Resources Research, v. 10, p. 969–973, doi: 10.1029/WR010i005p00969. Fuller, C.W., Willett, S.D., Fisher, D., and Lu, C.Y., 2006, A thermomechanical wedge model of Taiwan constrained by fission-track thermochronometry: Tectonophysics, v. 425, p. 1–24, doi: 10.1016/j.tecto.2006.05.018. Gaillardet, J., Calmels, D., Romero-Mujalli, G., Zakharova, E., and Hartmann, J., 2018, Global climate control on carbonate weathering intensity: Chemical Geology, p. 1–11, doi: 10.1016/j.chemgeo.2018.05.009. Gaillardet, J., Dupre, B., Allegre, C.J., and Négrel, P., 1997, Chemical and physical denudation in the Amazon River Basin: Chemical geology, v. 142, p. 141–173. Gaillardet, J., Dupré, B., and Louvat, P., 1999, ScienceDirect - Chemical Geology : Global silicate weathering and CO2 consumption rates deduced from the chemistry of large rivers: Chemical Geology, http://www.sciencedirect. com/science/article/pii/S0009254199000315%5Cnpapers2://publication/uuid/BF9B0411-7B1F-44AB-A416- 1B7D0528E13B. Gaillardet, J., Viers, J., and Dupré, B., 2013, Trace Elements in River Waters: v. 7, 195–235 p., doi: 10.1016/B978-0- 08-095975-7.00507-6. Gallen, S.F., 2018, Lithologic controls on landscape dynamics and aquatic species evolution in post-orogenic mountains: Earth and Planetary Science Letters, v. 493, p. 150–160, doi: 10.1016/j.epsl.2018.04.029. Galy, A., and France-Lanord, C., 1999, Weathering processes in the Ganges-Brahmaputra basin and the riverine alkalinity budget: Chemical Geology, v. 159, p. 31–60. Gandolfi, G., Paganelli, L., and Cavazza, W., 2007, Heavy-mineral associations as tracers of limited compositional mixing during turbiditic sedimentation of the Marnoso-Arenacea Formation (Miocene, Northern Apennines, Italy): Developments in Sedimentology, v. 58, p. 621–645. Granger, D.E., Kirchner, J.W., and Finkel, R., 1996, Spatially Averaged Long-Term Erosion Rates Measured from in Situ-Produced Cosmogenic Nuclides in Alluvial Sediment: The Journal of Geology, v. 104, p. 249–257, doi: 10.1086/629823. Granger, D.E., Kirchner, J.W., and Riebe, C.S., 1997, Inferring exhumation rates and processes from cosmogenic nuclides in sediments: Geological Society of America, 1997 annual meeting Abstracts with Programs - Geological Society of America,. Granger, D.E., and Riebe, C.S., 2013, Cosmogenic Nuclides in Weathering and Erosion: Elsevier Ltd., v. 7, 401–436 p., doi: 10.1016/B978-0-08-095975-7.00514-3. 15 Hack, J.T., 1975, Dynamic equilibrium and landscape evolution: Theories of landform development, v. 1, p. 87–102. Hack, J.T., 1957, Studies of longitudinal stream profiles in Virginia and Maryland: USGS Professional Paper 249, p. 97. Harkins, N., Kirby, E., Heimsath, A., Robinson, R., and Reiser, U., 2007, Transient fluvial incision in the headwaters of the Yellow River, northeastern Tibet, China: Journal of Geophysical Research: Earth Surface, v. 112, p. 1–21, doi: 10.1029/2006JF000570. Hilley, G.E., and Young, H.H., 2018, Millennial-scale denudation rates of the Santa Lucia Mountains, California: Implications for landscape evolution in steep, high-relief, coastal mountain ranges: GSA Bulletin, v. 130, p. 1809– 1824, doi: 10.1130/B31907.1. Jacobson, A.D., and Blum, J.D., 2003, Relationship between mechanical erosion and atmospheric CO 2 consumption in the New Zealand Southern Alps: , p. 865–868. Jacobson, A.D., Blum, J.D., Chamberlain, C.P., Craw, D., and Koons, P.O., 2003, Climatic and tectonic controls on chemical weathering in the New Zealand Southern Alps: Geochimica et Cosmochimica Acta, v. 67, p. 29–46. Kirby, E., and Whipple, K.X., 2012, Expression of active tectonics in erosional landscapes: Journal of Structural Geology, v. 44, p. 54–75, doi: 10.1016/j.jsg.2012.07.009. Malusà, M.G., and Balestrieri, M.L., 2012, Burial and exhumation across the Alps-Apennines junction zone constrained by fission-track analysis on modern river sands: Terra Nova, v. 24, p. 221–226, doi: 10.1111/j.1365- 3121.2011.01057.x. Meybeck, M., 1987, Global chemical weathering of surficial rocks estimated from river dissolved loads: American journal of science, v. 287, p. 401–428. Miller, S.R., Sak, P.B., Kirby, E., and Bierman, P.R., 2013, Neogene rejuvenation of central Appalachian topography: Evidence for differential rock uplift from stream profiles and erosion rates: Earth and Planetary Science Letters, v. 369–370, p. 1–12, doi: 10.1016/j.epsl.2013.04.007. Miller, S.R., Slingerland, R.L., and Kirby, E., 2007, Characteristics of steady state fluvial topography above fault-bend folds: Journal of Geophysical Research: Earth Surface, v. 112, p. 1–21, doi: 10.1029/2007JF000772. Molnar, P., and England, P., 1990, Late Cenozoic uplift of mountain ranges and global climate change: chicken or egg? Nature, v. 346, p. 29. Norton, K.P., von Blanckenburg, F., DiBiase, R., Schlunegger, F., and Kubik, P.W., 2011, Cosmogenic 10Be-derived denudation rates of the Eastern and Southern European Alps: International Journal of Earth Sciences, v. 100, p. 1163–1179, doi: 10.1007/s00531-010-0626-y. Ouimet, W.B., Whipple, K.X., and Granger, D.E., 2009, Beyond threshold hillslopes: Channel adjustment to base-level fall in tectonically active mountain ranges: Geology, v. 37, p. 579–582, doi: 10.1130/G30013A.1. Pazzaglia, F.J., and Brandon, M.T., 2001, A fluvial record of long-term steady-state uplift and erosion across the Cascadia forearc high, western Washington State: American Journal of Science, v. 301, p. 385–431, doi: 10.2475/ ajs.301.4-5.385. Penck, W., 1953, Morphological analysis of landforms: a contribution to physical geology, New York, St. Martin’s Press. Picotti, V., and Pazzaglia, F.J., 2008, A new active tectonic model for the construction of the Northern Apennines mountain front near Bologna (Italy): Journal of Geophysical Research: Solid Earth, v. 113, p. 1–24, doi: 10.1029/2007JB005307. Raymo, M.E., Ruddiman, W.F., and Froelich, P.N., 1988, Influence of late Cenozoic mountain building on ocean geochemical cycles: Geology, v. 16, p. 649–653. Reiners, P.W., and Brandon, M.T., 2006, Using Thermochronology To Understand Orogenic Erosion: Annual Review of Earth and Planetary Sciences, v. 34, p. 419–466, doi: 10.1146/annurev.earth.34.031405.125202. Ricci Lucchi, F., 1975, Miocene paleogeography and basin analysis in the Periadriatic Apennines: Petroleum Exploration Society of Libya: Ricci Lucchi, F., 1986, The Oligocene to Recent foreland basins of the northern Appennines: Foreland Basins, p. 105– 140. Riebe, C.S., Kirchner, J.W., and Finkel, R.C., 2004, Erosional and climatic effects on long-term chemical weathering rates in granitic landscapes spanning diverse climate regimes: Earth and Planetary Science Letters, v. 224, p. 547– 562, doi: 10.1016/j.epsl.2004.05.019. Riebe, C.S., Kirchner, J.W., Granger, D.E., and Finkel, R.C., 2001a, Minimal climatic control on erosion rates in the Sierra Nevada , California: , p. 447–450. Riebe, C.S., Kirchner, J.W., Granger, D.E., and Finkel, R.C., 2001b, Strong tectonic and weak climatic control of long term chemical weathering rates: GSA Bulletin, p. 511–514, doi: 10.1130/0091-7613(2001)029<0511:STAWCC>2.0.CO;2. Roller, S., Wittmann, H., Kastowski, M., and Hinderer, M., 2012, Erosion of the Rwenzori Mountains, East African Rift, from in situ-produced cosmogenic 10Be: Journal of Geophysical Research: Earth Surface, v. 117, doi: 10.1029/2011JF002117. Roveri, M., Bassetti, M.A., and Ricci Lucchi, F., 2001, The mediterranean Messinian salinity crisis: An Apennine 16 foredeep perspective: Sedimentary Geology, v. 140, p. 201–214, doi: 10.1016/S0037-0738(00)00183-4. Safran, E.B., Bierman, P.R., Aalto, R., Dunne, T., Whipple, K.X., and Caffee, M., 2005, Erosion rates driven by channel network incision in the Bolivian Andes: Earth Surface Processes and Landforms, v. 30, p. 1007–1024, doi: 10.1002/ esp.1259. Serpelloni, E., Faccenna, C., Spada, G., Dong, D., and Williams, S.D.P., 2013, Vertical GPS ground motion rates in the Euro-Mediterranean region: New evidence of velocity gradients at different spatial scales along the Nubia-Eurasia plate boundary: Journal of Geophysical Research: Solid Earth, v. 118, p. 6003–6024, doi: 10.1002/2013JB010102. Stallard, R.F., and Edmond, J.M., 1983, Oltman , Recent studies Meade et al sediment: Journal of Geophysical Research, v. 88, p. 9671–9688. Stock, J.D., and Montgomery, D.R., 1999, Geologic constraints on bedrock river incision using the stream power law: Journal of Geophysical Research: Solid Earth, v. 104, p. 4983–4993, doi: 10.1029/98JB02139. Thomson, S.N., Brandon, M.T., Reiners, P.W., Zattin, M., Isaacson, P.J., and Balestrieri, M.L., 2010, Thermochronologic evidence for orogen-parallel variability in wedge kinematics during extending convergent orogenesis of the northern Apennines, Italy: Bulletin of the Geological Society of America, v. 122, p. 1160–1179, doi: 10.1130/B26573.1. Vai, F., and Martini, I.P., 2001, Anatomy of an orogen: the Apennines and adjacent Mediterranean basins: Springer. Ventura, B., Pini, G.A., and Zuffa, G.G., 2001, Thermal history and exhumation of the Northern Apennines ( Italy ): evidence from combined apatitefissio-track and vitrinite reflectance data from foreland basin sediments.: Basin Research, p. 435–448. West, A.J., Galy, A., and Bickle, M., 2005, Tectonic and climatic controls on silicate weathering: Earth and Planetary Science Letters, v. 235, p. 211–228, doi: 10.1016/j.epsl.2005.03.020. Willett, S.D., and Brandon, M.T., 2002, On steady states in mountain belts: Geology, v. 30, p. 175–178, doi: 10.1130/0091-7613(2002)030<0175:OSSIMB>2.0.CO;2. Willett, S.D., Fisher, D., Fuller, C., Yeh, E.C., and Lu, C.Y., 2003, Erosion rates and orogenic-wedge kinematics in Taiwan inferred from fission-track thermochronometry: Geology, v. 31, p. 945–948, doi: 10.1130/G19702.1. Willett, S.D., Slingerland, R., and Hovius, N., 2001, Uplift, shortenning, and steady state topography in active mountain belts: American Journal of Science, v. 301, p. 455–485. Wittmann, H., Malusà, M.G., Resentini, A., Garzanti, E., and Niedermann, S., 2016, The cosmogenic record of mountain erosion transmitted across a foreland basin: Source-to-sink analysis of in situ 10Be, 26Al and 21Ne in sediment of the Po river catchment: Earth and Planetary Science Letters, v. 452, p. 258–271, doi: 10.1016/j.epsl.2016.07.017. Wobus, C., Whipple, K.X., Kirby, E., Snyder, N., Johnson, J., Spyropolou, K., Crosby, B., and Sheehan, D., 2006, Tectonics from topography: procedurses, promise, and pitfalls: Geological Society of America Special Paper, v. 398, p. 55–74, doi: 10.1130/2006.2398(04). Zattin, M., Picotti, V., and Zuffa, G.G., 2002, Fission-Track Reconstruction of the Front of the Northern Apennine Thrust Wedge and Overlying Ligurian Unit: v. 302, p. 346–379.
17 CHAPTER 2 The paradox of topography and erosion rates in the active orogen of the Northern Apennines
Erica Erlanger, Yanyan Wang, Sean D. Willett, Sean F. Gallen, Vincenzo Picotti
Submitted to Geology Journal
18 The paradox of topography and erosion rates in the active orogen of the Northern Apennines
Erica D. Erlanger1, Yanyan Wang1, Sean D. Willett1, Sean F. Gallen2, and Vincenzo Picotti1
1Department of Earth Sciences, ETH Zürich, Zürich, 8092, Switzerland
2Department of Geosciences, Colorado State University, Fort Collins 80523, CO
ABSTRACT The Northern Apennines are an active orogenic wedge with tectonic deformation driven by subduction of the Adriatic plate in the east and backarc spreading in the Ligurian Sea to the west. The orogen is characterized by: 1) Late Cenozoic subduction rates of 5-10 mm/yr; 2) geodetic extension rates of 2–4 mm/yr across the range, nearly balanced by contraction across the Adriatic mountain front; 3) geodetic uplift rates of 0.5–2.5 mm/yr in the NE, but stability or subsidence in the SW; 4) an asymmetric topography with the main divide offset to the west resulting in a steeper southwest-facing mountain flank; and 5) high erosion rates in the NE and low erosion rates in the SW based on detrital cosmogenic 10Be concentrations presented in this paper. These last two observations imply a negative correlation between erosion rate and metrics of landscape steepness. In this paper, we demonstrate that this apparent paradox of decoupling between topographic form and erosion rate can be resolved by a kinematic model of slab retreat, in which erosional flux is described by a vector with horizontal and vertical components. This model reconciles geodetic data, 10Be concentrations, and geomorphic observations.
INTRODUCTION Subduction of the Adriatic plate under the Northern Apennines has been continuous over the last 30 Ma, and nearly equal in rate to the rollback of the slab and trench with respect to stable Europe (Faccenna et al., 2014; Rosenbaum and Piana Agostinetti, 2015), offering an opportunity to study the crustal kinematics of a retreating subduction system. Deformation in the Northern Apennines is dominated by mid-crustal thrusting at the mountain front, and upper crustal normal faulting focused in the center of the range (Malinverno and Ryan, 1986; Cavinato and De Celles, 1999) Horizontal GPS measurements (Bennett et al., 2012; Devoti et al., 2017) illustrate NE motion with increasing velocity from SW to NE, indicating that the Ligurian and part of the Adriatic mountain flanks are under extension, while contraction occurs primarily at the mountain front and buried thrusts underlying the Po Plain sediments. Geodetic releveling data (D’Anastasio et al., 2006) indicate uplift focused on the Adriatic side, with a maximum rate of 2.5 mm/yr. The topography of the Northern Apennines illustrates a ridge crest that is offset to the west (Fig. 1A), resulting in a regionally steeper southwestern flank and a gentler northeastern flank, an observation that we quantify in this paper. The topographic asymmetry of the range is inconsistent with the pattern of geodetic uplift, as higher uplift rates normally imply steeper topography. In this study, we present new catchment- averaged erosion rates from major drainage basins in the Northern Apennines, supplementing published data, to evaluate how erosion rates compare to rock uplift rates and topographic metrics. Our results illustrate a disconnect between the present distribution of topography, geodetic uplift and apparent erosion rates, and we 19 show that these contradictory observations can be reconciled with a self-consistent kinematic model of the orogen.
TOPOGRAPHIC CHARACTERISTICS To demonstrate that the asymmetry of the Northern Apennines extends to the scale of geomorphic processes, we conducted two topographic analyses (Figs. 1B-1C), calculating the normalized length of the river, χ (Perron and Royden, 2013; Willett et al., 2014), and the channel steepness normalized for upstream drainage area (ksn) (Wobus et al., 2006). The Ligurian side is shorter, in normalized length, with lower values of χ at the divide (Fig. 1B), and steeper, with higher ksn values, particularly in the upper reaches. In contrast, the Adriatic side is longer, with larger χ values at the divide and lower average ksn values (Fig. 1C, Fig. DR1).
The contrast in χ at the water divide is large. Consistent with the contrast in ksn across the divide, these metrics imply that the water divide should move NE for long-term equilibrium to be obtained. Differential uplift rate, precipitation or rock erodibility can alter that expectation, but uplift rates are an order of magnitude higher on the Adriatic side (D’Anastasio et al., 2006) and precipitation is higher on the Ligurian side (Crespi et al., 2018), each of which should produce asymmetry with the opposite sense. Rock erodibility could partially explain the asymmetry and is discussed later. Numerous wind gaps attesting to river capture or basin beheading also support motion of the main Apennines divide to the NE (Ferraris et al., 2012; Spagnolo and Firpo, 2007).
EROSION RATES To test how modern erosion rates reflect the geomorphic asymmetry, we measuredin situ cosmogenic 10Be concentrations from fluvial quartz sediment, obtained from 11 catchments across the Northern Apennines. These supplement published data (Cyr and Granger, 2008; Wittmann et al., 2016)river incision, and uplift rates in the northern and central Apennines, Italy, since 0.9 Ma, are determined from new cosmogenic nuclide data. Beryllium-10 concentrations in modern and middle Pleistocene sediments indicate erosion rates from 0.20 to 0.58 mm/yr. These rates are similar to estimates of sediment yield (0.12-0.44 mm/yr and provide good spatial coverage of the range (Fig. 1D). Erosion rates (Methods, Data Repository) vary from 0.254 ± 0.044 mm/yr to 0.643 ± 0.230 mm/yr on the Adriatic side and from 0.105 ± 0.016 mm/yr to 0.284 ± 0.036 mm/yr on the Ligurian side (Fig. 1D, Table 1). Erosion rates for Ligurian catchments are a factor of 2-6 times lower than their Adriatic counterparts (Figs. 1 & 2).
PARADOX OF EROSION RATES AND GEOMORPHIC ASYMMETRY In general, geomorphic process laws predict a positive correlation between metrics of steepness, be they relief, slope or channel steepness (Wobus et al., 2006; Kirby and Whipple, 2012) and both rock uplift and erosion rates. In the Northern Apennines, the regional steepness pattern from ksn, χ plot slope, and 10 km-scale topographic slope show the opposite trend, with the basins of the Ligurian side exhibiting denudation rates systematically lower than denudation rates from the gentler basins on the Adriatic side (Fig. 2). Geodetic and geomorphic measures of rock uplift rate show the same discrepancy, with high uplift rates on the Adriatic side, and low uplift rates on the Ligurian side. Furthermore, the prediction of χ mapping is that the main divide is moving NE, a process normally associated with higher erosion rates in the aggressor basins, in this case those draining to the SW. Lithology is variable, exposing different sedimentary units across the range and meta-sedimentary mélanges and ophiolite units that are more prevalent on the Ligurian side, suggesting the harder lithologies could be associated with the steeper catchments. However, many of the upper catchments on the Adriatic side expose the same, harder lithologies, with no systematic increase in steepness. Spatial changes in ksn coincide 20 with the drainage divide, not major lithologic changes. Although we recognize that rock erodibility plays a role, it cannot be the exclusive control on denudation. Cosmogenic 10Be is accumulated over 103 to 104 yrs, so it is possible these data record only a short-term transient state, but the rates obtained here are consistent with paleo-erosion rates from Early-Middle Pleistocene marine terraces in the Romagna Apennines (Cyr and Granger, 2008) and thermochronometric ages from both sides of the Apennines (Thomson et al., 2010). These contradictory observations thus constitute a paradox that fails to explain the present distribution of topography, kinematics and geomorphic process rates in the Northern Apennines.
HORIZONTAL MOTION We propose that the key to reconciling the difference between erosion rates, uplift rates, and topographic steepness is to recognize that 10Be concentrations can be interpreted in a more general manner to account for tectonic kinematics. Assuming secular equilibrium, such that the total 10Be produced in a catchment is equal to the total 10Be exported, the concentration of 10Be is the ratio of the total 10Be produced and the total volume of quartz converted to sediment by erosion. The volume of quartz released per unit time is a mass flux and can be expressed as the product of the catchment area projected onto a horizontal surface, as is implicitly done in every DEM analysis, and the vertical rock velocity with respect to the surface, which is the erosion rate. However, rock motion with respect to the Earth’s surface is not always vertical, particularly in tectonically active settings where horizontal motion and deformation are important (Willett et al., 2001; Miller et al., 2007). In this case, the mass flux of quartz from the Earth should be calculated using a velocity with a non-zero horizontal component, so the catchment surface should be projected onto a plane orthogonal to the direction of the rock motion (Fig. DR2). As an end-member, we can consider pure horizontal motion of the rock relative to the surface, where the appropriate area for calculating the flux is the basin surface projected onto a vertical place. The horizontal velocity in this formulation is the average velocity of the Earth’s surface with respect to the rock, but does not imply uniform, horizontal motion of the surface at each point within the catchment. A typical catchment, comprised of channels and hillslopes at various orientations, can be regarded as eroding downward, but at differential rates along the catchment, such that the net result is apparent horizontal motion of the land surface at the kilometer scale (Willett et al., 2018). Cosmogenic 10Be production is thus unaffected by the assumption of horizontal motion, but the quartz or rock velocity implied by a given flux is affected. We demonstrate this concept by calculating an apparent velocity consistent with our 10Be concentrations, assuming a purely horizontal flux, where the catchment surface is projected onto a vertical plane with a strike direction parallel to the main divide (Methods, Data Repository). A positive flux requires that the velocity be directed away from the main divide, so taking a reference frame as positive to the NE implies that apparent velocities on the Ligurian side will be negative. We find apparent horizontal velocities ranging from -1.6 to -5.1 mm/yr for the Ligurian side and 5.3 to 21.1 mm/yr for the Adriatic side (Table 1). The more general case with vertical and horizontal velocity components is discussed below.
KINEMATIC MODEL These concepts and data can be unified in a kinematic model of an orogenic wedge abovea subducting, retreating slab (Fig. 3), where the prowedge (Adriatic side) and retrowedge (Ligurian side) are the accreting and non-accreting sides of the orogen, respectively, sensu Willett et al. (1993) (Methods, Data
Repository). Subduction is driven entirely by slab retreat at a velocity EVH (Fig. 3A). The relief formed by orogenesis sits above the slab hinge zone, with a topographic surface that is fixed in the vertical and has a horizontal velocity EVS relative to Eurasia. The surface velocity, EVS, must be equal to the retreat velocity, EVH, or the topography would be progressively offset from the subduction zone. Accretion of material into the wedge 21 implies motion of rock within the wedge with respect to the surface, SVW, so there is a non-zero horizontal component to the erosional flux. Horizontal GPS velocities across the Apennines provide an estimate of the wedge velocity, EVW, indicating an average of 2 to 4 mm/yr towards the NE, and this velocity must always be less than the hinge retreat velocity, EVH. These velocities can be expressed with respect to one another, recognizing the vector addition relationship:
EVH = EVS = EVW - SVW = EVW + WVS (1) The interpretation of detrital 10Be concentrations is consistent with the vertical and horizontal kinematics of Figure 3A when the non-vertical flux theory described above is applied. This is illustrated by showing the erosional flux as a function of the horizontal and vertical velocities (Fig. 3B). Each point on the y-axis represents a 10Be concentration converted to an apparent vertical erosion rate, as in conventional analysis (Table 1); the same measurement converted to an apparent horizontal velocity (Table 1) plots as a single point on the x-axis. Because velocity is defined to be positive to the NE, retrowedge points plot on the negative x-axis and prowedge points on the positive x-axis. The solid green and gray lines (Fig. 3B) are fit to the means of all 10Be measurements for the retrowedge and prowedge, respectively, and points on these lines represent combinations of horizontal and vertical velocities consistent with the 10Be measurements. The horizonal velocity SVW will always be negative, but given the opposite slopes of the wedge surface, SVW decreases the contribution of horizontal motion to the prowedge denudational flux and increases the contribution of horizontal motion to the retrowedge flux. Given our assumption of a common horizontal velocity across the wedge, but independent vertical velocities for the prowedge UP and retrowedge UR (gray and green bands), the intersection of the gray and green lines with any vertical line defines a kinematic model consistent with the 10 Be data. A horizontal velocity, SVW, between approximately -2 and -5 mm/yr, thus provides an acceptable fit to both geodetic uplift data and cosmogenic 10Be data. Horizontal surface motion of 2–5 mm/yr is sufficient to fully account for the 10Be concentrations in catchments of the retrowedge, suggesting the erosional flux of the retrowedge is driven entirely by the horizontal velocity. For the prowedge, the model suggests a combination of vertical and horizontal motion to account for both the growth of the topography and the erosional flux.
The estimate of SVW can also be combined with the GPS-based estimate of EVW to obtain the hinge retreat rate, EVH, using equation (1), giving an estimate of ~4 to 9 mm/yr. The late Cenozoic retreat rate is estimated independently as 6 to 10 mm/yr, dependent on latitude (Faccenna et al., 2014; Rosenbaum and Piana Agostinetti, 2015), so although uncertainties are large, these estimates are essentially equivalent, suggesting that geodetic data and erosional flux data are consistent with a steady rate of retreat over the late Cenozoic. The asymmetric topography of the Northern Apennines is consistent with a large component of horizontal motion. Long-wavelength asymmetry is a consequence of horizontal advection of topography balanced by steeper rivers (Willett et al., 2001). Miller et al. (2007) showed that the Siwalik Hills are similarly characterized by spatially variable vertical rock uplift rates, where higher channel steepness and concavity in distal (retrowedge) rivers were attributed to horizontal advection in the direction of streamflow. The degree of asymmetry is difficult to estimate as it depends on assumptions regarding erosion processes. For example, Willett et al. (2018) showed that the importance of topographic advection depends strongly on the slope exponent, n, in a stream-power incision law. In the Northern Apennines, the steeper Tyrhennian topography is consistent with horizontal advection driven by tectonic crustal accretion towards the SW (Fig. DR3, wedge schematics, Fig. 3B). Erosion on the retrowedge counters the divide motion, and the 10Be concentrations are consistent with a balance between the horizontal tectonic velocity and the velocity of the surface with respect to the underlying rock. The relative motion between the surface and the underlying rock gives the appearance of divide motion to the 22 NE, consistent with geomorphic observations, when, in fact, the divide distance to baselevel at either mountain front is steady in time. Like most paradoxes, the disconnect between the topography, the uplift rates, and the erosion rates in the Northern Apennines is simply an issue of the conceptual framing of the problem. In landscapes characterized by asymmetric topography, the direction of horizontal advection and balance between vertical and horizontal rock motion are manifest in the morphology of river profiles, the pattern of asymmetric topography, and as we have shown here, our interpretation of detrital 10Be concentrations in terms of a rock flux that has non-vertical components.
ACKNOWLEDGMENTS
This research is funded by the SINERGIA Swiss Alp Array project (SNF number 154434).
REFERENCES CITED Bennett, R.A., Serpelloni, E., Hreinsdóttir, S., Brandon, M.T., Buble, G., Basic, T., Casale, G., Cavaliere, A., Anzidei, M., Marjonovic, M., Minelli, G., Molli, G., and Montanari, A., 2012, Syn-convergent extension observed using the RETREAT GPS network, northern Apennines, Italy: Journal of Geophysical Research: Solid Earth, v. 117, doi: 10.1029/2011JB008744. Cavinato, G.P., and De Celles, P.G., 1999, Extensional basins in the tectonically bimodal central Apennines fold-thrust belt, Italy: Response to corner flow above a subducting slab in retrograde motion: Geology, v. 27, p. 955–958, doi: 10.1130/0091-7613(1999)027<0955:EBITTB>2.3.CO;2. Crespi, A., Brunetti, M., Lentini, G., and Maugeri, M., 2018, 1961–1990 high-resolution monthly precipitation climatologies for Italy: International Journal of Climatology, v. 38, p. 878–895, doi: 10.1002/joc.5217. Cyr, A.J., and Granger, D.E., 2008, Dynamic equilibrium among erosion, river incision, and coastal uplift in the northern and central Apennines, Italy: Geology, v. 36, p. 103–106, doi: 10.1130/G24003A.1. D’Anastasio, E., De Martini, P.M., Selvaggi, G., Pantosti, D., Marchioni, A., and Maseroli, R., 2006, Short- term vertical velocity field in the Apennines (Italy) revealed by geodetic levelling data: Tectonophysics, v. 418, p. 219–234, doi: 10.1016/j.tecto.2006.02.008. Devoti, R., D’Agostino, N., Serpelloni, E., Pietrantonio, G., Riguzzi, F., Avallone, A., Cavaliere, A., Cheloni, D., Cecere, G., D’Ambrosio, C., Falco, L., Selvaggi, G., Métois, M., Esposito, A., et al., 2017, A combined velocity field of the mediterranean region: Annals of Geophysics, v. 60, doi: 10.4401/ag- 7059. Faccenna, C., Becker, T.W., Miller, M.S., Serpelloni, E., and Willett, S.D., 2014, Isostasy, dynamic topography, and the elevation of the Apennines of Italy: Earth and Planetary Science Letters, v. 407, p. 163–174, doi: 10.1016/j.epsl.2014.09.027. Ferraris, F., Firpo, M., and Pazzaglia, F.J., 2012, DEM analyses and morphotectonic interpretation: The Plio- Quaternary evolution of the eastern Ligurian Alps, Italy: Geomorphology, v. 149–150, p. 27–40, doi: 10.1016/j.geomorph.2012.01.009. Kirby, E., and Whipple, K.X., 2012, Expression of active tectonics in erosional landscapes: Journal of Structural Geology, v. 44, p. 54–75, doi: 10.1016/j.jsg.2012.07.009. Malinverno, A., and Ryan, W.B.F., 1986, Extension in the Ligurian Sea and shortening in the Apennines as result of arc migration driven by sinking of the lithosphere: Tectonics, v. 5, p. 227–245. Miller, S.R., Slingerland, R.L., and Kirby, E., 2007, Characteristics of steady state fluvial topography
23 above fault-bend folds: Journal of Geophysical Research: Earth Surface, v. 112, p. 1–21, doi: 10.1029/2007JF000772. Perron, J.T., and Royden, L., 2013, An integral approach to bedrock river profile analysis: Earth Surface Processes and Landforms, v. 38, p. 570–576, doi: 10.1002/esp.3302. Rosenbaum, G., and Piana Agostinetti, N., 2015, Crustal and upper mantle responses to lithospheric segmentation in the northern Apennines: Tectonics, v. 34, p. 648–661, doi: 10.1002/2013TC003498. Spagnolo, M., and Firpo, M., 2007, Geomorphic evolution of the seaward escarpment in the NE Ligurian Alps (Italy): Zeitschrift für Geomorphologie, v. 51, p. 115–134, doi: 10.1127/0372-8854/2007/0051-0115. Thomson, S.N., Brandon, M.T., Reiners, P.W., Zattin, M., Isaacson, P.J., and Balestrieri, M.L., 2010, Thermochronologic evidence for orogen-parallel variability in wedge kinematics during extending convergent orogenesis of the northern Apennines, Italy: Bulletin of the Geological Society of America, v. 122, p. 1160–1179, doi: 10.1130/B26573.1. Willett, S., Beaumont, C., and Fullsack, P., 1993, Mechanical model for the tectonics of doubly vergent compressional orogens: Geology, v. 21, p. 371–374, doi: 10.1130/0091-7613(1993)021<0371:MMFTTO>2.3.CO. Willett, S.D., McCoy, S.W., and Beeson, H.W., 2018, Transience of the North American High Plains landscape and its impact on surface water: Nature, p. 1, doi: 10.1038/s41586-018-0532-1. Willett, S.D., McCoy, S.W., Taylor Perron, J., Goren, L., and Chen, C.Y., 2014, Dynamic reorganization of River Basins: Science, v. 343, doi: 10.1126/science.1248765. Willett, S.D., Slingerland, R., and Hovius, N., 2001, Uplift, shortenning, and steady state topography in active mountain belts: American Journal of Science, v. 301, p. 455–485. Wittmann, H., Malusà, M.G., Resentini, A., Garzanti, E., and Niedermann, S., 2016, The cosmogenic record of mountain erosion transmitted across a foreland basin: Source-to-sink analysis of in situ 10Be, 26Al and 21Ne in sediment of the Po river catchment: Earth and Planetary Science Letters, v. 452, p. 258–271, doi: 10.1016/j.epsl.2016.07.017. Wobus, C., Whipple, K.X., Kirby, E., Snyder, N., Johnson, J., Spyropolou, K., Crosby, B., and Sheehan, D., 2006, Tectonics from topography: procedurses, promise, and pitfalls: Geological Society of America Special Paper, v. 398, p. 55–74, doi: 10.1130/2006.2398(04).
24 N N N N 5 4 5 4 4 4 4 4 E 2 E 1 2 ** 1 21 a a Adriatic Adriatic n n g g ** ** o o l l 20 19 o o B B Ligurian e e ** c c 9 E 18 n n Erosion Rate (mm/yr) Rate (mm/yr) Erosion E 1 e e 1 1 r 0.10 0.60 r 1 o o l l F F 17 8 a a m m χ (m) # a r a r 16 7 P P # 6 15 14 0 22 E E 13 0 0 5 1 1 # 12 m 4 m 11 a a 0 0 l l 5 5 o o # s s 3 10 a s a s 5 5 a 2 2 n n a 1 # o o o o E E 2 n B B n 9 9 e e B D 0 0 G G E 2 E 1 2 1 a a n n g g ) o o 0.9 l l o o (m B B sn K Allochthonous Jurassic to e e c c E n n E 1 e e 0 279 1 1 r r uscan metamorphic rocks 1 Marnoso Arenacea Unit- turbidites foredeep Miocene Macigno/Cervarola Units- turbidites foredeep to Mid-Miocene Ligurian Unit - sedimentary deepwater Cenozoic early and ophiolites mélange - Epiligurian deposits and continental shallow marine T o o l l F F a a D m m a r a r P P C E E 0 0 1 1 m B m a 0 0 l 5 5 a l o s o s a s A 5 5 a s 2 2 n a a n o o o E E o B n n 9 9 B e e C A 0 0 G G N N N N 5 4 5 4 4 4 4 4
25 Figure 1. A: Simplified geologic map of the Northern Apennines and 15 km-wide swath profiles (black rectangles). Letters (A-D) on swath profiles correspond with the arrangement of the swath profiles in Figure 2. The Ligurian Unit is divided into the Internal Ligurian (IL) Unit (dark green color) and External Ligurian (EL) Unit (light green color).The primary drainage divide (solid black line) and thrust front (sawtooth line) are indicated. B: Normalized length parameter, χ, mapped over catchments using a base level of 50 m. Location of all figures (red box), Adriatic Sea, and Ligurian Sea are shown in the inset image.C: ksn map of catchments using a 1 km smoothing window. D: Compilation of catchment-averaged erosion rates from this study, Cyr and Granger (2008) (numbers with ** symbol), and Wittmann et al. (2016) (numbers with # symbol).
A 2 SW NE 3 Erosion Rate (mm/yr) Reference Velocity 2 (2 mm/yr) 1
Elevation (km) 1
0 0 B Uplift/Erosion Rate (mm/yr) 2 3
2
1
1 Elevation (km)
0 0 C 2 3 Erosion Rate (mm/yr)
2 1 1 Elevation (km)
0 0
D Uplift/Erosion Rate (mm/yr) 2 3
2 Reno 1 1 Elevation (km)
0 0 50 100 Distance (km) 26 § § 6.5 9.3 6.6 9.1 6.5 1.66 1.98 3.43 2.31 5.14 3.06 1.84 5.29 9.56 21.1 12.0 16.4 12.2 13.7 N.D N.D ------velocity velocity (mm/yr) Horizontal Horizontal ) ) vert § § 2 5.4 7.0 8.6 8.1 30.9 12.2 19.3 29.5 59.6 20.8 64.1 19.5 11.0 31.3 42.4 35.7 20.0 38.8 29.0 N.D N.D (km Projected catchment area (A § § 16° 16° 35° 35° 35° 35° 35° (°) 215° 215° 215° 196° 196° 196° 215° 215° 215° 215° 215° 215° N.D N.D vector Velocity Velocity orientation † § § § N.D N.D N.D (mm/yr) Erosion rate Erosion 0.277 0.109 ± 0.330 0.051 ± 0.382 0.072 ± 0.435 0.076 ±
0.119 0.014 ± 0.104 0.017 ± 0.217 0.031 ± 0.199 0.027 ± 0.293 0.041 ± 0.194 0.023 ± 0.138 0.021 ±
0.229 0.042 ± 0.415 0.206 ± 0.598 0.120 ± 0.519 0.129 ± 0.694 0.242 ± 0.372 0.061 ± 0.234 0.041 ± § § § *† mf N.D N.D N.D 0.047 0.049 0.047 0.048 0.046 0.044 0.045 0.049 0.048 0.048 0.045 0.044 0.046 0.046 0.047 0.048 0.046 0.044 P
§ § § *† ms N.D N.D N.D 0.018 0.019 0.017 0.018 0.016 0.015 0.016 0.019 0.018 0.018 0.016 0.015 0.017 0.017 0.017 0.018 0.017 0.015 (at/g/yr) P § § § *† n N.D N.D N.D 8.469 9.858 7.999 8.836 7.436 6.362 7.098 9.802 8.680 8.965 6.664 6.205 7.875 7.516 8.209 9.113 7.534 6.519 P
(mm/yr) Erosion rate rate Erosion 0.274 0.109 ± 0.298 0.046 ± 0.401 0.073 ± 0.430 0.076 ± 0.268 0.048 ± 0.442 0.224 ± 0.603 0.123 ± 0.497 0.124 ± 0.643 0.230 ± 0.358 0.058 ± 0.479 0.100 ± 0.311 0.037 ± 0.583 0.112 ± 0.108 0.013 ± 0.105 0.016 ± 0.216 0.032 ± 0.193 0.026 ± 0.284 0.036 ± 0.190 0.022 ± 0.136 0.021 ± 0.254 0.044 ± ) ) 2 plan area 122.0 680.5 308.3 264.2 614.3 917.6 480.5 992.1 989.9 132.5 179.2 190.1 296.8 555.4 946.8 251.7 316.2 149.4 264.2 Basin (km (A 1248.9 1160.5 840 877 848 883 (m) 685 843 716 743 819 665 598 627 591 558 521 619 985 825 885 559 641 Mean Elevation * mf 0.047 0.048 0.047 0.048 0.046 0.047 0.046 0.046 0.047 0.045 0.045 0.045 0.045 0.044 0.044 0.045 0.049 0.047 0.048 0.044 0.045 P
* ). ms mf 0.018 0.018 0.018 0.018 0.016 0.018 0.017 0.017 0.017 0.016 0.016 0.016 0.016 0.016 0.015 0.016 0.019 0.018 0.018 0.016 0.016 (at/g/yr) P * n 8.355 8.800 8.462 8.719 P 7.390 8.451 7.594 7.815 8.369 7.229 6.733 6.921 6.750 6.679 6.403 7.088 9.481 8.403 8.767 6.546 7.163 Be] (P muons fast and ) at/g) 3 ms 10 [ 9.0 1.5 ± 9.2 2.8 ± 8.4 1.3 ± (10 21.6 7.5 ± 20.8 2.4 ± 14.9 2.3 ± 14.3 2.0 ± 13.5 5.7 ± 11.2 2.4 ± 14.5 1.8 ± 10.2 1.9 ± 16.1 0.8 ± 44.7 2.5 ± 44.6 4.8 ± 23.6 2.4 ± 34.3 2.9 ± 20.9 2.0 ± 32.5 2.0 ± 35.0 3.7 ± 20.3 2.7 ± 19.8 2.8 ± (°E) 9.6118° 8.8598° 9.5849° 9.3616° 9.8569° 9.9258° 9.0576° 10.1692° 10.9240° 10.2381° 10.0913° 10.4083° 10.7571° 11.2573° 11.6213° 11.6865° 11.8849° 10.3762° 10.5060° 10.5539° 11.1258° Longitude ), slow muons (P muons slow ), n BE CONCENTRATIONS, EROSION RATES, AND HORIZONTAL VELOCITY CALCULATIONS CALCULATIONS VELOCITY AND HORIZONTAL RATES, EROSION CONCENTRATIONS, BE 10 (°N) Latitude 44.8327° 44.6322° 44.4201° 44.5716° 44.7256° 44.9061° 44.6953° 44.6222° 44.5322° 44.3923° 44.2211° 44.1698° 44.1187° 44.3518° 44.1950° 44.1867° 44.1360° 43.9299° 43.9993° 43.9278° 44.8926° ** ** ** ** Be production rates for neutrons (P neutrons for rates production Be 10 # # # # # Nure Baganza Panaro Parma River/Location Scrivia Trebbia Taro Enza Secchia Basin) (Whole Reno Senio at Casola Valsenio Lamone at San Eufemia Davadola at Montone Entella Vara Magra Serchio at Filicaia at Serchio Serchio at Piaggone at Serchio Lima Bisenzio Staffora 2 3 4 5 6 7 8 9 1 Catchment sampled by Cyr and Granger (2008). Granger and by Cyr sampled Catchment 11 12 13 16 17 19 20 21 Basin-averaged in catchment. distribution quartz uneven for corrected rates erosion and rates Production not determined. = N.D. 14 15 18 10 Catchment sampled by Wittmann et al. (2016). (2016). al. et by Wittmann sampled Catchment TABLE 1. CATCHMENT METRICS, TABLE * † § ** II. LiguriaII. and Tuscany III. Emilia Romagna III. I. Piemonte #
27 Figure 2. Swath profiles extending near to or from the Ligurian coastline to the Po Plain. The shaded boxes show erosion rates for the Ligurian side (blue) and Adriatic side (red); the height of the box indicates the erosion rate with uncertainties, and the width of box shows the distance along the profile where the rate is applicable. Dashed lines show geodetic uplift rates (D’Anastasio et al., 2006), and arrows above the profile show horizontal velocities from GPS measurements (Bennett et al., 2012). Outlined arrows indicate GPS stations located on the profile, and hatched arrows indicate the station and velocity vector were projected onto the profile from a maximum distance of 15 km.
A
SW EVS EVS NE EVW
V SVE S W Corsica SVE Zone
Material Spatial Hinge EVS Point Feature +V E W E Zero Velocity
Velocity
Frame Reference Velocity S Subduction Hinge Eurasia
B V SVW SVW = 0 S W Retro Pro Retro Pro Retro Pro 1
0 9
0 8
Acceptable 0 Range of sVw
0 Retrowedge Average UP 0
0
0 Uplift (mm/yr) Prowedge 0
0 1
0 10 8 0 8 10 1 1 1 18 0 1 Average UR