Refine
Year of publication
Document Type
- Article (163)
- Other (6)
- Postprint (5)
- Conference Proceeding (2)
- Review (2)
Keywords
- Himalaya (5)
- Iran (5)
- erosion (5)
- Andes (4)
- exhumation (4)
- Alborz Mountains (3)
- Central Andes (3)
- Lithospheric structure (3)
- Pamir (3)
- TerraceM (3)
Institute
- Institut für Geowissenschaften (178) (remove)
The mechanisms by which climate and vegetation affect erosion rates over various time scales lie at the heart of understanding landscape response to climate change. Plot-scale field experiments show that increased vegetation cover slows erosion, implying that faster erosion should occur under low to moderate vegetation cover. However, demonstrating this concept over long time scales and across landscapes has proven to be difficult, especially in settings complicated by tectonic forcing and variable slopes. We investigate this problem by measuring cosmogenic Be-10-derived catchment-mean denudation rates across a range of climate zones and hillslope gradients in the Kenya Rift, and by comparing our results with those published from the Rwenzori Mountains of Uganda. We find that denudation rates from sparsely vegetated parts of the Kenya Rift are up to 0.13 mm/yr, while those from humid and more densely vegetated parts of the Kenya Rift flanks and the Rwenzori Mountains reach a maximum of 0.08 mm/yr, despite higher median hillslope gradients. While differences in lithology and recent land-use changes likely affect the denudation rates and vegetation cover values in some of our studied catchments, hillslope gradient and vegetation cover appear to explain most of the variation in denudation rates across the study area. Our results support the idea that changing vegetation cover can contribute to complex erosional responses to climate or land-use change and that vegetation cover can play an important role in determining the steady-state slopes of mountain belts through its stabilizing effects on the land surface.
Uplift in the broken Andean foreland of the Argentine Santa Bárbara System (SBS) is associated with the contractional reactivation of basement anisotropies, similar to those reported from the thick-skinned Cretaceous-Eocene Laramide province of North America. Fault scarps, deformed Quaternary deposits and landforms, disrupted drainage patterns, and medium-sized earthquakes within the SBS suggest that movement along these structures may be a recurring phenomenon, with yet to be defined repeat intervals and rupture lengths. In contrast to the Subandes thrust belt farther north, where eastward-migrating deformation has generated a well-defined thrust front, the SBS records spatiotemporally disparate deformation along structures that are only known to the first order. We present herein the results of geomorphic desktop analyses, structural field observations, and 2D electrical resistivity tomography and seismic-refraction tomography surveys and an interpretation of seismic reflection profiles across suspected fault scarps in the sedimentary basins adjacent to the Candelaria Range (CR) basement uplift, in the south-central part of the SBS. Our analysis in the CR piedmont areas reveals consistency between the results of near-surface electrical resistivity and seismic-refraction tomography surveys, the locations of prominent fault scarps, and structural geometries at greater depth imaged by seismic reflection data. We suggest that this deformation is driven by deep-seated blind thrusting beneath the CR and associated regional warping, while shortening involving Mesozoic and Cenozoic sedimentary strata in the adjacent basins was accommodated by layer-parallel folding and flexural-slip faults that cut through Quaternary landforms and deposits at the surface.
Uplift in the broken Andean foreland of the Argentine Santa Bárbara System (SBS) is associated with the contractional reactivation of basement anisotropies, similar to those reported from the thick-skinned Cretaceous-Eocene Laramide province of North America. Fault scarps, deformed Quaternary deposits and landforms, disrupted drainage patterns, and medium-sized earthquakes within the SBS suggest that movement along these structures may be a recurring phenomenon, with yet to be defined repeat intervals and rupture lengths. In contrast to the Subandes thrust belt farther north, where eastward-migrating deformation has generated a well-defined thrust front, the SBS records spatiotemporally disparate deformation along structures that are only known to the first order. We present herein the results of geomorphic desktop analyses, structural field observations, and 2D electrical resistivity tomography and seismic-refraction tomography surveys and an interpretation of seismic reflection profiles across suspected fault scarps in the sedimentary basins adjacent to the Candelaria Range (CR) basement uplift, in the south-central part of the SBS. Our analysis in the CR piedmont areas reveals consistency between the results of near-surface electrical resistivity and seismic-refraction tomography surveys, the locations of prominent fault scarps, and structural geometries at greater depth imaged by seismic reflection data. We suggest that this deformation is driven by deep-seated blind thrusting beneath the CR and associated regional warping, while shortening involving Mesozoic and Cenozoic sedimentary strata in the adjacent basins was accommodated by layer-parallel folding and flexural-slip faults that cut through Quaternary landforms and deposits at the surface.
The Frolikha Fan : a large Pleistocene glaciolacustrine outwash fan in northern Lake Baikal, Siberia
(1998)
The removal, redistribution, and transient storage of sediments in tectonically active mountain belts is thought to exert a first-order control on shallow crustal stresses, fault activity, and hence on the spatiotemporal pattern of regional deformation processes. Accordingly, sediment loading and unloading cycles in intermontane sedimentary basins may inhibit or promote intrabasinal faulting, respectively, but unambiguous evidence for this potential link has been elusive so far. Here we combine 2D numerical experiments that simulate contractional deformation in a broken-foreland setting (i.e., a foreland where shortening is diachronously absorbed by spatially disparate, reverse faults uplifting basement blocks) with field data from intermontane basins in the NW Argentine Andes. Our modeling results suggest that thicker sedimentary fills (>0.7-1.0 km) may suppress basinal faulting processes, while thinner fills (<0.7 km) tend to delay faulting. Conversely, the removal of sedimentary loads via fluvial incision and basin excavation promotes renewed intrabasinal faulting. These results help to better understand the tectono-sedimentary history of intermontane basins that straddle the eastern border of the Andean Plateau in northwestern Argentina. For example, the Santa Maria and the Humahuaca basins record intrabasinal deformation during or after sediment unloading, while the Quebrada del Toro Basin reflects the suppression of intrabasinal faulting due to loading by coarse conglomerates. We conclude that sedimentary loading and unloading cycles may exert a fundamental control on spatiotemporal deformation patterns in intermontane basins of tectonically active broken forelands. (C) 2018 Elsevier B.V. All rights reserved.
Sedimentary basins in the interior of orogenic plateaus can provide unique insights into the early history of plateau evolution and related geodynamic processes. The northern sectors of the Iranian Plateau of the Arabia-Eurasia collision zone offer the unique possibility to study middle-late Miocene terrestrial clastic and volcaniclastic sediments that allow assessing the nascent stages of collisional plateau formation. In particular, these sedimentary archives allow investigating several debated and poorly understood issues associated with the long-term evolution of the Iranian Plateau, including the regional spatio-temporal characteristics of sedimentation and deformation and the mechanisms of plateau growth. We document that middle-late Miocene crustal shortening and thickening processes led to the growth of a basement-cored range (Takab Range Complex) in the interior of the plateau. This triggered the development of a foreland-basin (Great Pari Basin) to the east between 16.5 and 10.7Ma. By 10.7Ma, a fast progradation of conglomerates over the foreland strata occurred, most likely during a decrease in flexural subsidence triggered by rock uplift along an intraforeland basement-cored range (Mahneshan Range Complex). This was in turn followed by the final incorporation of the foreland deposits into the orogenic system and ensuing compartmentalization of the formerly contiguous foreland into several intermontane basins. Overall, our data suggest that shortening and thickening processes led to the outward and vertical growth of the northern sectors of the Iranian Plateau starting from the middle Miocene. This implies that mantle-flow processes may have had a limited contribution toward building the Iranian Plateau in NW Iran.
The idea that climatically modulated erosion may impact orogenic processes has challenged geoscientists for decades. Although modeling studies and physical calculations have provided a solid theoretical basis supporting this interaction, to date, field-based work has produced inconclusive results. The central-western Alborz Mountains in the northern sectors of the Arabia-Eurasia collision zone constitute a promising area to explore these potential feedbacks. This region is characterized by asymmetric precipitation superimposed on an orogen with a history of spatiotemporal changes in exhumation rates, deformation patterns, and prolonged, km-scale base-level changes. Our analysis suggests that despite the existence of a strong climatic gradient at least since 17.5 Ma, the early orogenic evolution (from similar to 36 to 9-6 Ma) was characterized by decoupled orographic precipitation and tectonics. In particular, faster exhumation and sedimentation along the more arid southern orogenic flank point to a north-directed accretionary flux and underthrusting of Central Iran. Conversely, from 6 to 3 Ma, erosion rates along the northern orogenic flank became higher than those in the south, where they dropped to minimum values. This change occurred during a similar to 3-Myr-long, km-scale base-level lowering event in the Caspian Sea. We speculate that mass redistribution processes along the northern flank of the Alborz and presumably across all mountain belts adjacent to the South Caspian Basin and more stable areas of the Eurasian plate increased the sediment load in the basin and ultimately led to the underthrusting of the Caspian Basin beneath the Alborz Mountains. This underthrusting in turn triggered a new phase of northward orogenic expansion, transformed the wetter northern flank into a new pro-wedge, and led to the establishment of apparent steady-state conditions along the northern orogenic flank (i.e., rock uplift equal to erosion rates). Conversely, the southern mountain front became the retro-wedge and experienced limited tectonic activity. These observations overall raise the possibility that mass-distribution processes during a pronounced erosion phase driven by base-level changes may have contributed to the inferred regional plate-tectonic reorganization of the northern Arabia-Eurasia collision during the last similar to 5 Ma. (C) 2015 Elsevier B.V. All rights reserved.
The Central Pontides of N Turkey represents a mobile orogenic belt of the southern Eurasian margin that experienced several phases of exhumation associated with the consumption of different branches of the Neo-Tethys Ocean and the amalgamation of continental domains. Our new low-temperature thermochronology data help to constrain the timing of these episodes, providing new insights into associated geodynamic processes. In particular, our data suggest that exhumation occurred at (1) similar to 110 to 90Ma, most likely during tectonic accretion and exhumation of metamorphic rocks from the subduction zone; (2) from similar to 60 to 40Ma, during the collision of the Kirehir and Anatolide-Tauride microcontinental domains with the Eurasian margin; (3) from similar to 0 to 25Ma, either during the early stages of the Arabia-Eurasia collision (soft collision) when the Arabian passive margin reached the trench, implying 70 to 530km of subduction of the Arabian passive margin, or during a phase of trench advance predating hard collision at similar to 20Ma; and (4) similar to 11Ma to the present, during transpression associated with the westward motion of Anatolia. Our findings document the punctuated nature of fault-related exhumation, with episodes of fast cooling followed by periods of slow cooling or subsidence, the role of inverted normal faults in controlling the Paleogene exhumation pattern, and of the North Anatolian Fault in dictating the most recent pattern of exhumation.
The southern foreland basin of the Alborz Mountains of northern Iran is characterized by an approximately 7.3-km-thick sequence of Miocene sedimentary rocks, constituting three basin-wde coarsening-upward units spanning a period of 10(6)years. We assess available magnetostratigraphy, paleoclimatic reconstructions, stratal architecture, records of depositional environments, and sediment-provenance data to characterize the relationships between tectonically-generated accommodation space (A) and sediment supply (S). Our analysis allows an inversion of the stratigraphy for particular forcing mechanisms, documenting causal relationships, and providing a basis to decipher the relative contributions of tectonics and climate (inferred changes in precipitation) in controlling sediment supply to the foreland basin. Specifically, A/S>1, typical of each basal unit (17.5-16.0, 13.8-13.1 and 10.3-9.6Ma), is associated with sharp facies retrogradation and reflects substantial tectonic subsidence. Within these time intervals, arid climatic conditions, changes in sediment provenance, and accelerated exhumation in the orogen suggest that sediment supply was most likely driven by high uplift rates. Conversely, A/S<1 (13.8 and 13.8-11Ma, units 1, and 2) reflects facies progradation during a sharp decline in tectonic subsidence caused by localized intra-basinal uplift. During these time intervals, climate continued to be arid and exhumation active, suggesting that sediment supply was again controlled by tectonics. A/S<1, at 11-10.3Ma and 9-6-7.6Ma (and possibly 6.2; top of units 2 and 3), is also associated with two episodes of extensive progradation, but during wetter phases. The first episode appears to have been linked to a pulse in sediment supply driven by an increase in precipitation. The second episode reflects a balance between a climatically-induced increase in sediment supply and a reduction of subsidence through the incorporation of the proximal foreland into the orogenic wedge. This in turn caused an expansion of the catchment and a consequent further increase in sediment supply.
A poorly understood lag time of 15-20 m.y. exists between the initial Arabia-Eurasia continental collision in late Eocene to early Oligocene time and the acceleration of tectonic and sedimentary processes across the collision zone in the early to late Miocene. The late Eocene to Miocene-Pliocene clastic and shallow-marine sedimentary rocks of the Kond, Eyvanekey, and Semnan Basins in the Alborz Mountains (northern Iran) offer the possibility to track the evolution of this orogen in the framework of collision processes. A transition from volcaniclastic submarine deposits to shallow-marine evaporites and terrestrial sediments occurred shortly after 36 Ma in association with reversals in sediment provenance, strata tilting, and erosional unroofing. These events followed the termination of subduction arc magmatism and marked a changeover from an extensional to a contractional regime in response to initiation of continental collision with the subduction of stretched Arabian lithosphere. This early stage of collision produced topographic relief associated with shallow foreland basins, suggesting that shortening and tectonic loading occurred at low rates. Starting from the early Miocene (17.5 Ma), flexural subsidence in response to foreland basin initiation occurred. Fast sediment accumulation rates and erosional unroofing trends point to acceleration of shortening by the early Miocene. We suggest that the lag time between the initiation of continental collision (36 Ma) and the acceleration of regional deformation (20-17.5 Ma) reflects a two-stage collision process, involving the "soft" collision of stretched lithosphere at first and "hard" collision following the arrival of unstretched Arabian continental litho sphere in the subduction zone.
Tectonic and climatic control on evolution of rift lakes in the Central Kenya Rift, East Africa
(2009)
The long-term histories of the neighboring Nakuru-Elmenteita and Naivasha lake basins in the Central Kenya Rift illustrate the relative importance of tectonic versus climatic effects on rift-lake evolution and the formation of disparate sedimentary environments. Although modem climate conditions in the Central Kenya Rift are very similar for these basins, hydrology and hydrochemistry of present-day lakes Nakuru, Elmenteita and Naivasha contrast dramatically due to tectonically controlled differences in basin geometries, catchment size, and fluvial processes. In this study, we use eighteen C-14 and Ar-40/Ar-39 dated fluvio-lacustrine sedimentary sections to unravel the spatiotemporal evolution of the lake basins in response to tectonic and climatic influences. We reconstruct paleoclimatic and ecological trends recorded in these basins based on fossil diatom assemblages and geologic field mapping. Our study shows a tendency towards increasing alkalinity and shrinkage of water bodies in both lake basins during the last million years. Ongoing volcano-tectonic segmentation of the lake basins, as well as reorganization of upstream drainage networks have led to contrasting hydrologic regimes with adjacent alkaline and freshwater conditions. During extreme wet periods in the past, such as during the early Holocene climate optimum, lake levels were high and all basins evolved toward freshwater systems. During drier periods some of these lakes revert back to alkaline conditions, while others maintain freshwater characteristics. Our results have important implications for the use and interpretation of lake sediment as climate archives in tectonically active regions and emphasize the need to deconvolve lacustrine records with respect to tectonics versus climatic forcing mechanisms.
Terrigenous sediment supply, marine transport, and depositional processes along tectonically active margins are key to decoding turbidite successions as potential archives of climatic and seismic forcings. Sequence stratigraphic models predict coarse-grained sediment delivery to deep-marine sites mainly during sea-level fall and lowstand. Marine siliciclastic deposition during transgressions and highstands has been attributed to sustained connectivity between terrigenous sources and marine sinks facilitated by narrow shelves. To decipher the controls on Holocene highstand turbidite deposition, we analyzed 12 sediment cores from spatially discrete, coeval turbidite systems along the Chile margin (29 degrees-40 degrees S) with changing climatic and geomorphic characteristics but uniform changes in sea level. Sediment cores from intraslope basins in north-central Chile (29 degrees-33 degrees S) offshore a narrow to absent shelf record a shut-off of turbidite deposition during the Holocene due to postglacial aridification. In contrast, core sites in south-central Chile (36 degrees-40 degrees S) offshore a wide shelf record frequent turbidite deposition during highstand conditions. Two core sites are linked to the Biobio river-canyon system and receive sediment directly from the river mouth. However, intraslope basins are not connected via canyons to fluvial systems but yield even higher turbidite frequencies. High sediment supply combined with a wide shelf and an undercurrent moving sediment toward the shelf edge appear to control Holocene turbidite sedimentation and distribution. Shelf undercurrents may play an important role in lateral sediment transport and supply to the deep sea and need to be accounted for in sediment-mass balances.
Much progress has been made in estimating recurrence intervals of great and giant subduction earthquakes using terrestrial, lacustrine, and marine paleoseismic archives. Recent detailed records suggest these earthquakes may have variable recurrence periods and magnitudes forming supercycles. Understanding seismic supercycles requires long paleoseismic archives that record timing and magnitude of such events. Turbidite paleoseismic archives may potentially extend past earthquake records to the Pleistocene and can thus complement commonly shorter-term terrestrial archives. However, in order to unambiguously establish recurring seismicity as a trigger mechanism for turbidity currents, synchronous deposition of turbidites in widely spaced, isolated depocenters has to be ascertained. Furthermore, characteristics that predispose a seismically active continental margin to turbidite paleoseismology and the correct sample site selection have to be taken into account.
Here we analyze 8 marine sediment cores along 950 km of the Chile margin to test for the feasibility of compiling detailed and continuous paleoseismic records based on turbidites. Our results suggest that the deposition of areally widespread, synchronous turbidites triggered by seismicity is largely controlled by sediment supply and, hence, the climatic and geomorphic conditions of the adjacent subaerial setting. The feasibility of compiling a turbidite paleoseismic record depends on the delicate balance between sufficient sediment supply providing material to fail frequently during seismic shaking and sufficiently low sedimentation rates to allow for coeval accumulation of planktonic foraminifera for high-resolution radiocarbon dating.
We conclude that offshore northern central Chile (29-32.5 degrees S) Holocene turbidite paleoseismology is not feasible, because sediment supply from the semi-arid mainland is low and almost no Holocene turbidity-current deposits are found in the cores. In contrast, in the humid region between 36 and 38 degrees S frequent Holocene turbidite deposition may generally correspond to paleoseismic events. However, high terrigenous sedimentation rates prevent high-resolution radiocarbon dating. The climatic transition region between 32.5 and 36 degrees S appears to be best suited for turbidite paleoseismology. (C) 2015 Elsevier Ltd. All rights reserved.
Newly acquired high-resolution bathymetric data (with 5 m and 2 m grid sizes) from the continental shelf off Concepcion (Chile), in combination with seismic reflection profiles, reveal a distinctly different evolution for the Biobio submarine canyon compared to that of one of its tributaries. Both canyons are incised into the shelf of the active margin. Whereas the inner shelf appears to be mantled with unconsolidated sediment, the outer shelf shows the influence of strong bottom currents that form drifts of loose sediment and transport -material into the Biobio submarine canyon and onto the continental slope.
The main stem of the Biobio Canyon is connected to the mouth of the Biobio River and currently provides a conduit for terrestrial sediment from the continental shelf to the deep seafloor. In contrast, the head of its tributary closest to the coast is located similar to 24 km offshore of the present-day coastline at 120 m water depth, and it is subject to passive sedimentation. However, canyon activity within the study area is interpreted to be controlled not only by the direct input of fluvial sediments into the canyon head facilitated by the river-mouth to canyon-head connection, but also by input from southward-directed bottom currents and possibly longshore drift. In addition, about 24 km offshore of the present-day coastline, the main stem of the Biobio Canyon has steep canyon walls next to sites of active tectonic deformation that are prone to wall failure. Mass-failure events may also foster turbidity currents and contribute to canyon feeding. In contrast, the tributary has less steep canyon walls with limited evidence of canyon-wall failure and is located down-system of bottom currents from the Biobio Canyon. It consequently receives neither fluvial nor longshore sediments. Therefore, the canyon's connectivity to fluvial or longshore sediment delivery pathways is affected by the distance of the canyon head from the coastline and the orientation of the canyon axis relative to the direction of bottom currents.
The ability of a submarine canyon to act as an active conduit for large quantities of terrestrial sediment toward the deep sea during sea-level highstands may be controlled by several different conditions simultaneously. These include bottom current direction, structural deformation of the seafloor affecting canyon location and orientation as well as canyon-wall failure, shelf gradient and associated distance from the canyon head to the coast, and fluvial networks. The complex interplay between these factors may vary even within an individual canyon system, resulting in distinct levels of canyon activity on a regional scale.
Understanding how Earth-surface processes respond to past climatic perturbations is crucial for making informed predictions about future impacts of climate change on sediment "uxes. Sedimentary records provide the archives for inferring these processes, but their interpretation is compromised by our incomplete understanding of how sediment-routing systems respond to millennial-scale climate cycles. We analyzed seven sediment cores recovered from marine turbidite depositional sites along the Chile continental margin. The sites span a pronounced arid-to-humid gradient with variable relief and related sediment connectivity of terrestrial and marine environments. These sites allowed us to study event related depositional processes in different climatic and geomorphic settings from the Last Glacial Maximum to the present day. The three sites reveal a steep decline of turbidite deposition during deglaciation. High rates of sea-level rise postdate the decline in turbidite deposition. Comparison with paleoclimate proxies documents that the spatio-temporal sedimentary pattern rather mirrors the deglacial humidity decrease and concomitant warming with no resolvable lag times. Our results let us infer that declining deglacial humidity decreased "uvial sediment supply. This signal propagated rapidly through the highly connected systems into the marine sink in north-central Chile. In contrast, in south-central Chile, connectivity between the Andean erosional zone and the "uvial transfer zone probably decreased abruptly by sediment trapping in piedmont lakes related to deglaciation, resulting in a sudden decrease of sediment supply to the ocean. Additionally, reduced moisture supply may have contributed to the rapid decline of turbidite deposition. These different causes result in similar depositional patterns in the marine sinks. We conclude that turbiditic strata may constitute reliable recorders of climate change across a wide range of climatic zones and geomorphic conditions. However, the underlying causes for similar signal manifestations in the sinks may differ, ranging from maintained high system connectivity to abrupt connectivity loss. (C) 2017 Elsevier B.V. All rights reserved.
Major earthquakes ( M > 8) have repeatedly ruptured the Nazca-South America plate interface of south-central Chile involving meter scale land-level changes. Earthquake recurrence intervals, however, extending beyond limited historical records are virtually unknown, but would provide crucial data on the tectonic behavior of forearcs. We analyzed the spatiotemporal pattern of Holocene earthquakes on Santa Maria Island (SMI; 37 degrees S), located 20 km off the Chilean coast and approximately 70 km east of the trench. SMI hosts a minimum of 21 uplifted beach berms, of which a subset were dated to calculate a mean uplift rate of 2.3 +/- 0.2 m/ky and a tilting rate of 0.022 +/- 0.002 degrees/ky. The inferred recurrence interval of strandline-forming earthquakes is similar to 180 years. Combining coseismic uplift and aseismic subsidence during an earthquake cycle, the net gain in strandline elevation in this environment is similar to 0.4 m per event
Aluminum-26 and beryllium-10 surface exposure dating on cut-and-fill river-terrace surfaces from the lower Sutlej Valley (northwest Himalaya) documents the close link between Indian Summer Monsoon (ISM) oscillations and intervals of enhanced fluvial incision. During the early Holocene ISM optimum, precipitation was enhanced and reached far into the internal parts of the orogen. The amplified sediment flux from these usually dry but glaciated areas caused alluviation of downstream valleys up to 120 m above present grade at ca. 9.9 k.y. B.P. Terrace formation (i.e., incision) in the coarse deposits occurred during century-long weak ISM phases that resulted in reduced moisture availability and most likely in lower sediment flux. Here, we suggest that the lower sediment flux during weak ISM phases allowed rivers to incise episodically into the alluvial fill
The tectonic and climatic boundary conditions of the broken foreland and the orogen interior of the southern Central Andes of northwestern Argentina cause strong contrasts in elevation, rainfall, and surface-process regimes. The climatic gradient in this region ranges from the wet, windward eastern flanks (similar to 2 m/yr rainfall) to progressively drier western basins and ranges (similar to 0.1 m/yr) bordering the arid Altiplano-Puna Plateau. In this study, we analyze the impact of spatiotemporal climatic gradients on surface erosion: First, we present 41 new catchment-mean erosion rates derived from cosmogenic nuclide inventories to document spatial erosion patterns. Second, we re-evaluate paleoclimatic records from the Calchaquies basin (66 W, 26 S), a large intermontane basin bordered by high (> 4.5 km) mountain ranges, to demonstrate temporal variations in erosion rates associated with changing climatic boundary conditions during the late Pleistocene and Holocene. Three key observations in this region emphasize the importance of climatic parameters on the efficiency of surface processes in space and time: (1) First-order spatial patterns of erosion rates can be explained by a simple specific stream power (SSP) approach. We explicitly account for discharge by routing high-resolution, satellite derived rainfall. This is important as the steep climatic gradient results in a highly non-linear relation between drainage area and discharge. This relation indicates that erosion rates (ER) scale with ER similar to SSP1.4 on cosmogenic-nuclide time scales. (2) We identify an intrinsic channel-slope behavior in different climatic compartments. Channel slopes in dry areas (< 0.25 m/yr rainfall) are slightly steeper than in wet areas (> 0.75 m/yr) with equal drainage areas, thus compensating lower amounts of discharge with steeper slopes. (3) Erosion rates can vary by an order of magnitude between presently dry (similar to 0.05 mm/yr) and well-defined late Pleistocene humid (similar to 0.5 mm/yr) conditions within an intemontane basin. Overall, we document a strong climatic impact on erosion rates and channel slopes. We suggest that rainfall reaching areas with steeper channel slopes in the orogen interior during wetter climate periods results in intensified sediment mass transport, which is primarily responsible for maintaining the balance between surface uplift, erosion, sediment routing and transient storage in the orogen.
The interplay between topography and Indian summer monsoon circulation profoundly controls precipitation distribution, sediment transport, and river discharge along the Southern Himalayan Mountain Front (SHF). The Higher Himalayas form a major orographic barrier that separates humid sectors to the south and and regions to the north. During the Indian summer monsoon, vortices transport moisture from the Bay of Bengal, swirl along the SHF to the northwest, and cause heavy rainfall when colliding with the mountain front. In the eastern and central parts of the Himalaya, precipitation measurements derived from passive microwave analysis (SSM/I) show a strong gradient, with high values at medium elevations and extensive penetration of moisture along major river valleys into the orogen. The end of the monsoonal conveyer belt is near the Sutlej Valley in the NW Himalaya, where precipitation is lower and rainfall maxima move to lower elevations. This region thus comprises a climatic transition zone that is very sensitive to changes in Indian summer monsoon strength. To constrain magnitude, temporal, and spatial distribution of precipitation, we analyzed high-resolution passive microwave data from the last decade and identified an abnormal monsoon year (AMY) in 2002. During the 2002 AMY, violent rainstorms conquered orographic barriers and penetrated far into otherwise and regions in the northwest Himalaya at elevations in excess of 3 km asl. While precipitation in these regions was significantly increased and triggered extensive erosional processes (i.e., debris flows) on sparsely vegetated, steep hillslopes, mean rainfall along the low to medium elevations was not significantly greater in magnitude. This shift may thus play an important role in the overall sediment flux toward the Himalayan foreland. Using extended precipitation and sediment flux records for the last century, we show that these events have a decadal recurrence interval during the present-day monsoon circulation. Hence, episodically occurring AMYs control geomorphic processes primarily in the high-elevation and sectors of the orogen, while annual recurring monsoonal rainfall distribution dominates erosion in the low- to medium- elevation parts along the SHF. (C) 2004 Elsevier B.V. All rights reserved
Late Quaternary intensified monsoon phases control landscape evolution in the northwest Himalaya
(2005)
The intensity of the Asian summer-monsoon circulation varies over decadal to millennial time scales and is reflected in changes in surface processes, terrestrial environments, and marine sediment records. However, the mechanisms of long-lived (2-5 k.y.) intensified monsoon phases, the related changes in precipitation distribution, and their effect on landscape evolution and sedimentation rates are not yet well understood. The and high-elevation sectors of the orogen correspond to a climatically sensitive zone that currently receives rain only during abnormal (i.e., strengthened) monsoon seasons. Analogous to present-day rainfall anomalies, enhanced precipitation during an intensified monsoon phase is expected to have penetrated far into these geomorphic threshold regions where hillslopes are close to the angle of failure. We associate landslide triggering during intensified monsoon phases with enhanced precipitation, discharge, and sediment flux leading to an increase in pore-water pressure, lateral scouring of rivers, and over- steepening of hillslopes, eventually resulting in failure of slopes and exceptionally large mass movements. Here we use lacustrine deposits related to spatially and temporally clustered large landslides (>0.5 km(3)) in the Sutlej Valley region of the northwest Himalaya to calculate sedimentation rates and to infer rainfall patterns during late Pleistocene (29-24 ka) and Holocene (10-4 ka) intensified monsoon phases. Compared to present-day sediment-flux measurements, a fivefold increase in sediment-transport rates recorded by sediments in landslide-dammed lakes characterized these episodes of high climatic variability. These changes thus emphasize the pronounced imprint of millennial-scale climate change on surface processes and landscape evolution
[1] The Puna-Altiplano plateau in South America is a high-elevation, low internal relief landform that is characterized by internal drainage and hyperaridity. Thermochronologic and sedimentologic observations from the Sierra de Calalaste region in the southwestern Puna plateau, Argentina, place new constraints on early plateau evolution by resolving the timing of uplift of mountain ranges that bound present-day basins and the filling pattern of these basins during late Eocene-Miocene time. Paleocurrent indicators, sedimentary provenance analyses, and apatite fission track thermochronology indicate that the original paleodrainage setting was disrupted by exhumation and uplift of the Sierra de Calalaste range between 24 and 29 Ma. This event was responsible for basin reorganization and the disruption of the regional fluvial system that has ultimately led to the formation of internal drainage conditions, which, in the Salar de Antofalla, were established not later than late Miocene. Upper Eocene-Oligocene sedimentary rocks flanking the range contain features that suggest an arid environment existed prior to and during its uplift. Provenance data indicate a common similar source located to the west for both the southern Puna and the Altiplano of Bolivia during the late Eocene- Oligocene with sporadic local sources. This suggests the existence of an extensive, longitudinally oriented foreland basin along the central Andes during this time
Intramontane sedimentary basins along the margin of continental plateaus often preserve strata that contain fundamental information regarding the pattern of orogenic growth. The sedimentary record of the elastic Miocene-Pliocene sequence deposited in the Fiambala Basin, at the southern margin of the Puna Plateau (NW Argentina), documents the late Miocene paleodrainage evolution from headwaters to the west, towards headwaters in the ranges that constitute the border of the Puna Plateau to the north. Apatite Fission track (AFT) thermochronology of sedimentary and basement rocks show that the southern Puna Plateau was the source for the youngest, middle Miocene, detrital population detected in late Miocene rocks; and that the margin of the Puna Plateau expressed a high relief, possibly similar to or higher than at present, by late Miocene time. Cooling ages obtained from basement rocks at the southern Puna margin suggest that exhumation started in the Oligocene and continued until the middle Miocene. We interpret the basin reorganization and the creation of a high relief plateau margin to be the direct response of the source-basin system to a wholesale surface uplift event that may have occurred during the late Cenozoic in the Puna-Altiplano region. At this time coeval paleodrainage reorganization is observed not only in the Fiambala Basin, but also in different basins along the southern and eastern Puna margin, suggesting a genetic link between the last stage of plateau formation and basin response. However, this event did not cause sufficient exhumation of basin bounding ranges to be recorded by AFT thermochronology. Our new data thus document a decoupling between late Cenozoic surface uplift and exhumation in the southern Puna Plateau. High relief achieved at the Puna margin by late Miocene time is linked to Oligocene-Miocene exhumation; no significant erosion (< 3 km) has occurred since in this and highland.
Recent studies have shown that the 1976-77 global climate shift strongly affected the South American climate. In our study, we observed a link between this climate shift and river-discharge variability in the subtropical Southern Central Andes. We analyzed the daily river-discharge time series between 1940 and 1999 from small to medium mountain drainage basins (10(2)-10(4) km(2) ) across a steep climatic and topographic gradient. We document that the discharge frequency distribution changed significantly, with higher percentiles exhibiting more pronounced trends. A change point between 1971 and 1977 marked an intensification of the hydrological cycle, which resulted in increased river discharge. In the upper Rio Bermejo basin of the northernmost Argentine Andes, the mean annual discharge increased by 40% over 7 years. Our findings are important for flood risk management in areas impacted by the 1976-77 climate shift; discharge frequency distribution analysis provides important insights into the variability of the hydrological cycle in the Andean realm.
The eastern flanks of the Central Andes are characterized by deep convection, exposing them to hydrometeorological extreme events, often resulting in floods and a variety of mass movements. We assessed the spatiotemporal pattern of rainfall trends and the changes in the magnitude and frequency of extreme events (ae<yen>95th percentile) along an E-W traverse across the southern Central Andes using rain-gauge and high-resolution gridded datasets (CPC-uni and TRMM 3B42 V7). We generated different climate indices and made three key observations: (1) an increase of the annual rainfall has occurred at the transition between low (< 0.5 km) and intermediate (0.5-3 km) elevations between 1950 and 2014. Also, rainfall increases during the wet season and, to a lesser degree, decreases during the dry season. Increasing trends in annual total amounts characterize the period 1979-2014 in the arid, high-elevation southern Andean Plateau, whereas trend reversals with decreasing annual total amounts were found at low elevations. (2) For all analyzed periods, we observed small or no changes in the median values of the rainfall-frequency distribution, but significant trends with intensification or attenuation in the 95th percentile. (3) In the southern Andean Plateau, extreme rainfall events exhibit trends towards increasing magnitude and, to a lesser degree, frequency during the wet season, at least since 1979. Our analysis revealed that low (< 0.5 km), intermediate (0.5-3 km), and high-elevation (> 3 km) areas respond differently to changing climate conditions, and the transition zone between low and intermediate elevations is characterized by the most significant changes.
This study analyzes the discharge variability of small to medium drainage basins (10(2)-10(4) km(2)) in the southern Central Andes of NW Argentina. The Hilbert-Huang Transform (HHT) was applied to evaluate non-stationary oscillatory modes of variability and trends, based on four time series of monthly normalized discharge anomaly between 1940 and 2015. Statistically significant trends reveal increasing discharge during the past decades and document an intensification of the hydrological cycle during this period. An Ensemble Empirical Mode Decomposition (EEMD) analysis revealed that discharge variability in this region can be best described by five quasi-periodic statistically significant oscillatory modes, with mean periods varying from 1 to 20 y. Moreover, we show that discharge variability is most likely linked to the phases of the Pacific Decadal Oscillation (PDO) at multi-decadal timescales (similar to 20 y) and, to a lesser degree, to the Tropical South Atlantic SST anomaly (TSA) variability at shorter timescales (similar to 2-5 y). Previous studies highlighted a rapid increase in discharge in the southern Central Andes during the 1970s, inferred to have been associated with the global 1976-77 climate shift. Our results suggest that the rapid discharge increase in the NW Argentine Andes coincides with the periodic enhancement of discharge, which is mainly linked to a negative to positive transition of the PDO phase and TSA variability associated with a long-term increasing trend. We therefore suggest that variations in discharge in this region are largely driven by both natural variability and the effects of global climate change. We furthermore posit that the links between atmospheric and hydrologic processes result from a combination of forcings that operate on different spatiotemporal scales. (C) 2017 Elsevier B.V. All rights reserved.
The timing and pattern of surface uplift of Miocene marine sediments capping the southern margin of the Central Anatolian Plateau in southern Turkey provide a first-order constraint on possible mechanisms of regional uplift. Nannofossil, ostracod, and planktic foraminifera biostratigraphy of the Basyayla section (Mut-Ermenek Basin) within the Mut and Kfiselerli Formations suggests a Tortonian age for marine sediments unconformably capping basement rocks at 2 km elevation. The identification of biozone MMi 12a (7.81-8.35 Ma) from planktic foraminifera in the upper part of the section provides the tightest constraint on the age, which is further limited to 8.35-8.108 Ma as a result of the reverse polarity of the collected samples (chron 4r.1 r or 4r.2r). This provides a limiting age for the onset of surface uplift at the margin of one of the world's major orogenic plateaus, from which an average uplift rate of 0.24-0.25 mm/yr can be calculated.
Subhorizontal beds of the uppermost marine sediments exposed throughout the Mut-Ermenek Basin suggest minimal localized deformation, with just minor faulting at the basin margin and broad antiformal deformation across the basin. This implies that the post-8 Ma uplift mechanism must be rooted deep within the crust or in the upper mantle. Published Pn-wave velocity data for the region are compatible with topography compensated by asthenosphere across the southern margin of the plateau, showing a close match to the highest topography when elevations are filtered with a 100-km-wide smoothing window. Uplift along the southern margin of the Central Anatolian Plateau is also reflected by the pattern of Miocene marine sediments capping the margin, which form an asymmetric drape fold over the topography. These observations, together with tomographic evidence for slab steepening and break-off beneath the Eastern Anatolian Plateau, suggest that at least some of the 2 km of post-8 Ma uplift of the southern Central Anatolian Plateau margin is compensated by low-density asthenospheric mantle that upwelled following slab break-off.
The arid Puna plateau of the southern Central Andes is characterized by Cenozoic distributed shortening forming intramontane basins that are disconnected from the humid foreland because of the defeat of orogen-traversing channels. Thick Tertiary and Quaternary sedimentary fills in Puna basins have reduced topographic contrasts between the compressional basins and ranges, leading to a typical low-relief plateau morphology. Structurally identical basins that are still externally drained straddle the eastern border of the Puna and document the eastward propagation of orographic barriers and ensuing aridification. One of them, the Angastaco basin, is transitional between the highly compartmentalized Puna highlands and the undeformed Andean foreland. Sandstone petrography, structural and stratigraphic analysis, combined with detrital apatite fission-track thermochronology from a similar to 6200-m-thick Miocene to Pliocene stratigraphic section in the Angastaco basin, document the late Eocene to late Pliocene exhumation history of source regions along the eastern border of the Puna (Eastern Cordillera (EC)) as well as the construction of orographic barriers along the southeastern flank of the Central Andes. Onset of exhumation of a source in the EC in late Eocene time as well as a rapid exhumation of the Sierra de Luracatao (in the EC) at about 20 Ma are recorded in the detrital sediments of the Angastaco basin. Sediment accumulation in the basin began similar to 15 Ma, a time at which the EC had already built sufficient topography to prevent Puna sourced detritus from reaching the basin. After similar to 13 Ma, shortening shifted eastward, exhuming ranges that preserve an apatite fission-track partial annealing zone recording cooling during the late Cretaceous rifting event. Facies changes and fossil content suggest that after 9 Ma, the EC constituted an effective orographic barrier that prevented moisture penetration into the plateau. Between 3.4 and 2.4 Ma, another orographic barrier was uplifted to the east, leading to further aridification and pronounced precipitation gradients along the mountain front. This study emphasizes the important role of tectonics in the evolution of climate in this part of the Andes
[ 1] For the Puna Plateau and Eastern Cordillera of NW Argentina, the temporal and spatial pattern of deformation and surface uplift remain poorly constrained. Analysis of completely and partially reset apatite fission track samples collected from vertical profiles along an ESE trending transect extending from the plateau interior across the southern Eastern Cordillera at similar to 25 degrees S reveals important constraints on the deformation and exhumation history of this part of the Andes. The data constrain the Neogene Andean development of the Eastern Cordillera as well as rift-related exhumation for some of the sampled locations in the Late Jurassic/Early Cretaceous. An intervening Eocene-Oligocene exhumation episode in the southern Eastern Cordillera was probably related to crustal shortening. Subsequent reburial of the area by Andean foreland basin strata commenced between 30 and 25 Myr. Magnitude and duration of sedimentation, revealed by thermal modeling, differ between the sample locations, pointing to an eastward propagating basin system. In the southern Eastern Cordillera, Andean deformation commenced at 22.5 - 21 Myr, predating both the inferred formation of significant topography by 5 - 7.5 Myr and preservation of sediments in the adjacent Cenozoic basins by 6.5 - 8 Myr. Comparing the calculated structural depth of partially reset samples suggests that newly formed west dipping reverse faults along the former Salta Rift margin accommodated most of the Neogene tectonic movement. Late Cenozoic deformation at the southern Eastern Cordillera began earlier in the west and subsequently propagated eastward. The lateral growth of the orogen is coupled with a foreland basin system developing in front of the range and then becomes subsequently compartmentalized by later emergent topography.
In the Himalaya of Chamba, NW India, a major orographic barrier in front of the Greater Himalayan Range extracts a high proportion of the monsoonal rainfall along its southern slopes and effectively shields the orogen interior from moisture-bearing winds. Along a similar to 100-km-long orogen perpendicular transect, 28 new apatite fission track (AFT) and 30 new zircon (U-Th)/He (ZHe) cooling ages reveal marked variations in age distributions and long-term exhumation rates between the humid frontal range and the semi-arid orogen interior. On the southern topographic front, very young, elevation-invariant AFT ages of <4 Ma have been obtained that are concentrated in a similar to 30-km-wide zone; 1-D-thermal modeling suggests a Plio-Pleistocene mean erosion rate of 0.8-1.9 mm yr(-1). In contrast, AFT and ZHe ages within the orogen interior are older (4-9 and 7-18 Ma, respectively), are positively correlated with sample elevation, and yield lower mean erosion rates (0.3-0.9 mm yr(-1)). Protracted low exhumation rates within the orogen interior over the last similar to 15 Myr prevailed contemporaneously with overall humid conditions and an effective erosional regime within the southern Himalaya. This suggests that the frontal Dhauladar Range was sufficiently high during this time to form an orographic barrier, focusing climatically enhanced erosional processes and tectonic deformation there. Thrusting along the two frontal range-bounding thrust, the Main Central Thrust and the Main Boundary Thrusts, was initiated at least similar to 15 Ma ago and has remained localized since then. The lack of evidence for localized uplift farther north indicates either a rather flat decollement with no ramp or the absence of active duplex systems beneath the interior of Chamba. Exhumational variability within Chamba is best explained as the result of continuous thrusting along a major basal decollement, with a flat beneath the slowly exhuming internal compartments and a steep frontal ramp at the rapidly exhuming frontal range. The pattern in Chamba contrasts with what is observed elsewhere along the Himalaya, where exhumation is focused in a zone similar to 150 km north of the orogenic front. In the NW Himalaya, preserved High Himalayan Crystalline nappes and Lesser Himalayan windows alternate on a relatively small scale of <100 km; these alternations are closely correlated with the pattern of exhumation. Although the spatial distribution of high-exhumation zones varies considerably between individual Himalayan sectors, all of these zones are closely correlated with locally higher rock-uplift rates, sharp topographic discontinuities, and focused orographic precipitation, suggesting strong feedbacks between tectonically driven rock uplift, orographically enhanced precipitation, and erosional processes.
Deciphering the response of sediment routing systems to climatic forcing is fundamental for understanding the impacts of climate change on landscape evolution. In the Kangra Basin (northwest Sub-Himalaya, India), upper Pleistocene to Holocene alluvial fills and fluvial terraces record periodic fluctuations of sediment supply and transport capacity on timescales of 10(3) to 10(5) yr. To evaluate the potential influence of climate change on these fluctuations, we compare the timing of aggradation and incision phases recorded within remnant alluvial fans and terraces with climate archives. New surface-exposure dating of six terrace levels with in-situ cosmogenic Be-10 indicates the onset of incision phases. Two terrace surfaces from the highest level (T1) sculpted into the oldest preserved alluvial fan (AF1) date back to 53.4 +/- 3.2 ka and 43.0 +/- 2.7 ka (1 sigma). T2 surfaces sculpted into the remnants of AF1 have exposure ages of 18.6 +/- 1.2 ka and 15.3 +/- 0.9 ka, while terraces sculpted into the upper Pleistocene-Holocene fan (AF2) provide ages of 9.3 +/- 0.4 ka (T3), 7.1 +/- 0.4 ka (T4), 5.2 +/- 0.4 ka (T5) and 3.6 +/- 0.2 ka (T6). Together with previously published OSL ages yielding the timing of aggradation, we find a correlation between variations in sediment transport with oxygen-isotope records from regions affected by the Indian Summer Monsoon. During periods of increased monsoon intensity and post-Last Glacial Maximum glacial retreat, aggradation occurred in the Kangra Basin, likely due to high sediment flux, whereas periods of weakened monsoon intensity or lower sediment supply coincide with incision. (C) 2016 Elsevier B.V. All rights reserved.
The southernmost thrust of the Himalayan orogenic wedge that separates the foreland from the orogen, the Main Frontal Thrust, is thought to accommodate most of the ongoing crustal shortening in the Sub-Himalaya. Steepened longitudinal river profile segments, terrace offsets, and back-tilted fluvial terraces within the Kangra reentrant of the NW Sub-Himalaya suggest Holocene activity of the Jwalamukhi Thrust (JMT) and other thrust faults that may be associated with strain partitioning along the toe of the Himalayan wedge. To assess the shortening accommodated by the JMT, we combine morphometric terrain analyses with in situ Be-10-based surface-exposure dating of the deformed terraces. Incision into upper Pleistocene sediments within the Kangra Basin created two late Pleistocene terrace levels (T1 and T2). Subsequent early Holocene aggradation shortly before similar to 10ka was followed by episodic reincision, which created four cut-and-fill terrace levels, the oldest of which (T3) was formed at 10.10.9ka. A vertical offset of 445m of terrace T3 across the JMT indicates a shortening rate of 5.60.8 to 7.51.1mma(-1) over the last similar to 10ka. This result suggests that thrusting along the JMT accommodates 40-60% of the total Sub-Himalayan shortening in the Kangra reentrant over the Holocene. We speculate that this out-of-sequence shortening may have been triggered or at least enhanced by late Pleistocene and Holocene erosion of sediments from the Kangra Basin.
The northward movement and collision of the Arabian plate with Eurasia generates compressive stresses and resulting shortening in Iran. Within the Alborz Mountains, North Iran, a complex and not well understood system of strike-slip and thrust faults accomodates a fundamental part of the NNE-SSW oriented shortening. On 28th of May 2004 the Mw 6.3 Baladeh earthquake hit the north-central Alborz Mountains. It is one of the rare and large events in this region in modern time and thus a seldom chance to study earthquake mechanisms and the local ongoing deformation processes. It also demonstrated the high vulnerability of this densily populated region.
The M (w) 6.2 Baladeh earthquake occurred on 28 May 2004 in the Alborz Mountains, northern Iran. This earthquake was the first strong shock in this intracontinental orogen for which digital regional broadband data are available. The Baladeh event provides a rare opportunity to study fault geometry and ongoing deformation processes using modern seismological methods. A joint inversion for hypocentres and a velocity model plus a surface-wave group dispersion curve analysis were used to obtain an adapted velocity model, customised for mid- and long-period waveform modelling. Based on the new velocity model, regional waveform data of the mainshock and larger aftershocks (M (w) a parts per thousand yen3.3) were inverted for moment tensors. For the Baladeh mainshock, this included inversion for kinematic parameters. All analysed earthquakes show dominant thrust mechanisms at depths between 14 and 26 km, with NW-SE striking fault planes. The mainshock ruptured a 28A degrees south-dipping area of 24 x 21 km along a north-easterly direction. The rupture plane of the mainshock does not coincide with the aftershock distribution, neither in map view nor with respect to depth. The considered aftershocks form two main clusters. The eastern cluster is associated with the mainshock. The western cluster does not appear to be connected with the rupture plane of the mainshock but, instead, indicates a second activated fault plane dipping at 85A degrees towards the north.
Erweiterte Momententensorinversion und ihre seismotektonische Anwendung : Elbursgebirge, Nordiran
(2009)
Der Elburs im Norden Irans ist ein durch die Konvergenz der Arabischen und Eurasischen Platte verursachtes doppelt konvergentes Gebirge. Das komplexe System von Blattverschiebungen und Überschiebungen sowie die Aufnahme der Deformation im Elburs ist noch nicht sehr gut verstanden. Eine neu zu entwicklende Methode zur Inversion von seismischen Momententensoren, die unterschiedliche Beobachtungen verschiedener Stationstypen kombiniert invertiert, soll die bisher hauptsächlich strukturelle/geomorphologische Datengrundlage um Momententensoren auch kleinerer Magnituden (M < 4.5) erweitern. Dies ist die notwendige Grundlage für detaillierte seismotektonische Studien, die wiederum die Basis für seismische Gefährdungsanalysen bilden.
The complex system of strike-slip and thrust faults in the Alborz Mountains, Northern Iran, are not well understood yet. Mainly structural and geomorphic data are available so far. As a more extensive base for seismotectonic studies and seismic hazard analysis we plan to do a comprehensive seismic moment tensor study also from smaller magnitudes (M < 4.5) by developing a new algorithm. Here, we present first preliminary results.
The Gofa Province and the Chew Bahir Basin of southern Ethiopia constitute tectonically active regions, where the Southern Main Ethiopian Rift converges with the Northern Kenya Rift through a wide zone of extensional deformation with several north to northeast-trending, left-stepping en-e & PRIME;chelon basins. This sector of the Southern Main Ethiopian Rift is characterized by a semi-arid climate and a largely uniform lithology, and thus provides ideal conditions for studying the different parameters that define the tectonic and geomorphic features of this complex kinematic transfer zone. In this study, the degree of tectonic activity, spatiotemporal variations in extension, and the nature of kinematic linkage between different fault systems of the transfer zone are constrained by detailed quantitative geomorphic analysis of river catchments and focused field work. We analyzed fluvial and landscape morphometric characteristics in combination with structural, seismicity, and climatic data to better evaluate the tectono-geomorphic history of this transfer zone. Our data reveal significant north-south variations in the degree of extension from the Sawula Basin in the north (mature) to the Chew Bahir Basin in the south (juvenile). First, normalized channel-steepness indices and the spatial arrangement of knickpoints in footwall-draining streams suggest a gradual, southward shift in extensional deformation and recent tectonic activity. Second, based on 1-k(m) radius local relief and mean-hillslope maximum values that are consistent with ksn anomalies, we confirm strain localization within zones of fault interaction. Third, morphometric indices such as hypsometry, basin asymmetry factor, and valley floor width to valley height ratio also indicate a north to south gradient in tectonic activity, highlighting the importance of such a wide transfer zone with diffuse extension linking different rift segments during the break-up of continental crust.
Despite a large number of dated glacial landforms in the Himalaya, the ice extent during the global Last Glacial Maximum (LGM) from 19 to 23 ka is only known to first order. New cosmogenic Be-10 exposure ages from well-preserved glacially polished surfaces, combined with published data, and an improved production rate scaling model allow reconstruction of the LGM ice extent and subsequent deglaciation in the Chandra Valley of NW India. We show that a >1000 m thick valley glacier retreated >150 km within a few thousand years after the onset of LGM deglaciation. By comparing the recession of the Chandra Valley Glacier and other Himalayan glaciers with those of Northern and Southern Hemisphere glaciers, we demonstrate that post-LGM deglaciation was similar and nearly finished prior to the Bolling/Allerod interstadial. Our study supports the view that many Himalayan glaciers advanced during the LGM, likely in response to global variations in temperature.
Despite remarkable tectonostratigraphic similarities along the Himalayan arc, pronounced topographic and exhumational variability exists in different morphotectonic segments. The processes responsible for this segmentation are debated. Of particular interest is a 30- to 40-km-wide orogen-parallel belt of rapid exhumation that extends from central Nepal to the western Himalaya and its possible linkage to a midcrustal ramp in the basal decollement, and the related growth of Lesser Himalayan duplex structures. Here we present 26 new apatite fission track cooling ages from the Beas-Lahul region, at the transition from the Central to the Western Himalaya (77 degrees-78 degrees E) to investigate segmentation in the Himalayan arc from a thermochronologic perspective. Together with previously published data from this part of the orogen, we document significant lateral changes in exhumation between the Dhauladar Range to the west, the Beas-Lahul region, and the Sutlej area to the east of the study area. In contrast to the Himalayan front farther east, exhumation in the far western sectors is focused at the frontal parts of the mountain range and associated with the hanging wall of the Main Boundary Thrust fault ramp. Our results allow us to spatially correlate the termination of the rapid exhumation belt with a midcrustal ramp to the west. We suggest that a plunging anticline at the northwestern edge of the Larji-Kullu-Rampur window represents the termination of the Central Himalayan segment, which is related to the evolution of the Lesser Himalayan duplex. Key Points
Alpine glacial erosion exerts a first-order control on mountain topography and sediment production, but its mechanisms are poorly understood. Observational data capable of testing glacial erosion and transport laws in glacial models are mostly lacking. New insights, however, can be gained from detrital tracer thermochronology. Detrital tracer thermochronology works on the premise that thermochronometer bedrock ages vary systematically with elevation, and that detrital downstream samples can be used to infer the source elevation sectors of sediments. We analyze six new detrital samples of different grain sizes (sand and pebbles) from glacial deposits and the modern river channel integrated with data from 18 previously analyzed bedrock samples from an elevation transect in the Leones Valley, Northern Patagonian Icefield, Chile (46.7 degrees S). We present 622 new detrital zircon (U-Th)/He (ZHe) single-grain analyses and 22 new bedrock ZHe analyses for two of the bedrock samples to determine age reproducibility. Results suggest that glacial erosion was focused at and below the Last Glacial Maximum and neoglacial equilibrium line altitudes, supporting previous modeling studies. Furthermore, grain age distributions from different grain sizes (sand, pebbles) might indicate differences in erosion mechanisms, including mass movements at steep glacial valley walls. Finally, our results highlight complications and opportunities in assessing glacigenic environments, such as dynamics of sediment production, transport, transient storage, and final deposition, that arise from settings with large glacio-fluvial catchments.
The thick-skinned fold-and-thrust belt on the eastern flank of the Andean Plateau in northwestern Argentina (NWA) is a zone of active contractional deformation characterized by fault-bounded mountain ranges with no systematic spatiotemporal pattern of tectonic activity. In contrast, the thin-skinned Subandean fold-and-thrust belt of northern Argentina and southern Bolivia is characterized primarily by in-sequence (i.e., west to east) fault progression, with a narrow zone of Quaternary deformation focused at the front of the orogenic wedge. To better understand how recent deformation is accommodated across these mountain ranges and the Argentinian portion of the orogen in particular, estimating and comparing deformation rates and patterns across different timescales is essential. We present Late Pleistocene shortening rates for the central Calchaqui intermontane valley in NWA associated with at least three episodes of deformation. Global Positioning System data for the same region reveal a gradual decrease in horizontal surface velocities from the Eastern Cordillera toward the foreland, which contrasts with the rapid velocity gradient associated with a locked decollement in the Subandean Ranges of southern Bolivia. Our new results represent a small view of regional deformation that, when considered in combination with the shallow crustal seismicity and decadal-scale surface velocities, support the notion that strain release in NWA is associated with numerous slowly deforming structures that are distributed throughout the orogen.
Tectonically active coasts are dynamic environments characterized by the presence of multiple marine terraces formed by the combined effects of wave erosion, tectonic uplift, and sea-level oscillations at glacialcycle timescales. Well-preserved erosional terraces from the last interglacial sea-level highstand are ideal marker horizons for reconstructing past sea-level positions and calculating vertical displacement rates. We carried out an almost continuous mapping of the last interglacial marine terrace along similar to 5000 km of the western coast of South America between 1 degrees N and 40 degrees S. We used quantitatively replicable approaches constrained by published terrace-age estimates to ultimately compare elevations and patterns of uplifted terraces with tectonic and climatic parameters in order to evaluate the controlling mechanisms for the formation and preservation of marine terraces and crustal deformation. Uncertainties were estimated on the basis of measurement errors and the distance from referencing points. Overall, our results indicate a median elevation of 30.1 m, which would imply a median uplift rate of 0.22 m kyr(-1) averaged over the past similar to 125 kyr. The patterns of terrace elevation and uplift rate display high-amplitude (similar to 100-200 m) and long-wavelength (similar to 10(2) km) structures at the Manta Peninsula (Ecuador), the San Juan de Marcona area (central Peru), and the Arauco Peninsula (south-central Chile). Medium-wavelength structures occur at the Mejillones Peninsula and Topocalma in Chile, while short-wavelength (< 10 km) features are for instance located near Los Vilos, Valparaiso, and Carranza, Chile. We interpret the long-wavelength deformation to be controlled by deep-seated processes at the plate interface such as the subduction of major bathymetric anomalies like the Nazca and Carnegie ridges. In contrast, short-wavelength deformation may be primarily controlled by sources in the upper plate such as crustal faulting, which, however, may also be associated with the subduction of topographically less pronounced bathymetric anomalies. Latitudinal differences in climate additionally control the formation and preservation of marine terraces. Based on our synopsis we propose that increasing wave height and tidal range result in enhanced erosion and morphologically well-defined marine terraces in south-central Chile. Our study emphasizes the importance of using systematic measurements and uniform, quantitative methodologies to characterize and correctly interpret marine terraces at regional scales, especially if they are used to unravel the tectonic and climatic forcing mechanisms of their formation. This database is an integral part of the World Atlas of Last Interglacial Shorelines (WALIS), published online at https://doi.org/10.5281/zenodo.4309748 (Freisleben et al., 2020).
We analyzed the Lomas de Carabajal area in the intermontane Lerma valley of the Cordillera Oriental to assess the level of neotectonic activity in a densely populated region of northwestern Argentina. In this region, Plio-Pleistocene synorogenic conglomerates are deformed, locally associated with high-angle faults, and NNW-SSE oriented en-echelon folds characterized by wavelengths of < 1 km. The deformed Quaternary units follow the same pattern of deformation as observed in the underlying Neogene deposits; growth-strata geometries are observed near faults. This configuration is compatible with local left-lateral transpressional tectonism driven by ENE-WSW buttressing against the NW-oriented border of a Cretaceous extensional basin (Alemania sub-basin). Optically Stimulated Luminescence analysis of sandy-silty layers interbedded within the folded late Pleistocene conglomeratic sequence helps to determine uplift rates of 0.83-0.87 mm/a during the last 30-40 ka. Nearby the Lomas de Carabajal, a WNW-striking, 3-m-high fault scarp disrupts radiocarbon dated, 10-ka-old loessic deposits providing a Holocene mean uplift rate of 0.30 mm/a. Our data unambiguously show that shallow crustal deformation in the intermontane Lerma valley is ongoing; some of this deformation may be associated with seismicity. Our findings support the notion of temporally and spatially disparate deformation processes in the broken foreland of the northwestern Argentinean Andes.
A potential human footprint on Western Central African rainforests before the Common Era has become the focus of an ongoing controversy. Between 3,000 y ago and 2,000 y ago, regional pollen sequences indicate a replacement of mature rainforests by a forest-savannah mosaic including pioneer trees. Although some studies suggested an anthropogenic influence on this forest fragmentation, current interpretations based on pollen data attribute the "rainforest crisis" to climate change toward a drier, more seasonal climate. A rigorous test of this hypothesis, however, requires climate proxies independent of vegetation changes. Here we resolve this controversy through a continuous 10,500-y record of both vegetation and hydrological changes from Lake Barombi in Southwest Cameroon based on changes in carbon and hydrogen isotope compositions of plant waxes. delta C-13-inferred vegetation changes confirm a prominent and abrupt appearance of C-4 plants in the Lake Barombi catchment, at 2,600 calendar years before AD 1950 (cal y BP), followed by an equally sudden return to rainforest vegetation at 2,020 cal y BP. delta D values from the same plant wax compounds, however, show no simultaneous hydrological change. Based on the combination of these data with a comprehensive regional archaeological database we provide evidence that humans triggered the rainforest fragmentation 2,600 y ago. Our findings suggest that technological developments, including agricultural practices and iron metallurgy, possibly related to the large-scale Bantu expansion, significantly impacted the ecosystems before the Common Era.
The 'wet' early to mid-Holocene of tropical Africa, with its enhanced monsoon, ended with an abrupt shift toward drier conditions and was ultimately replaced by a drier climate that has persisted until the present day. The forcing mechanisms, the timing, and the spatial extent of this major climatic transition are not well understood and remain the subject of ongoing research. We have used a detailed palaeo-shoreline record from Lake Turkana (Kenya) to decipher and characterise this marked climatic transition in East Africa. We present a high-precision survey of well-preserved palaeo-shorelines, new radiocarbon ages from shoreline deposits, and oxygen-isotope measurements on freshwater mollusk shells to elucidate the Holocene moisture history from former lake water-levels in this climatically sensitive region. In combination with previously published data our study shows that during the early Holocene the water-level in Lake Turkana was high and the lake overflowed temporarily into the White Nile drainage system. During the mid-Holocene (similar to 5270 +/- 300 cal. yr BP), however, the lake water-level fell by similar to 50 m, coeval with major episodes of aridity on the African continent. A comparison between palaeo-hydrological and archaeological data from the Turkana Basin suggests that the mid-Holocene climatic transition was associated with fundamental changes in prehistoric cultures, highlighting the significance of natural climate variability and associated periods of protracted drought as major environmental stress factors affecting human occupation in the East African Rift System. (
The African Humid Period (AHP) between similar to 15 and 5.5 cal. kyr BP caused major environmental change in East Africa, including filling of the Suguta Valley in the northern Kenya Rift with an extensive (similar to 2150 km(2)), deep (similar to 300 m) lake. Interfingering fluvio-lacustrine deposits of the Baragoi paleo-delta provide insights into the lake-level history and how erosion rates changed during this time, as revealed by delta-volume estimates and the concentration of cosmogenic Be-10 in fluvial sand. Erosion rates derived from delta-volume estimates range from 0.019 to 0.03 mm yr(-1). Be-10-derived paleo-erosion rates at similar to 11.8 cal. kyr BP ranged from 0.035 to 0.086 mm yr(-1), and were 2.7 to 6.6 times faster than at present. In contrast, at similar to 8.7 cal. kyr BP, erosion rates were only 1.8 times faster than at present. Because Be-10-derived erosion rates integrate over several millennia; we modeled the erosion-rate history that best explains the 10Be data using established non-linear equations that describe in situ cosmogenic isotope production and decay. Two models with different temporal constraints (15-6.7 and 12-6.7 kyr) suggest erosion rates that were 25 to 300 times higher than the initial erosion rate (pre-delta formation). That pulse of high erosion rates was short (similar to 4 kyr or less) and must have been followed by a rapid decrease in rates while climate remained humid to reach the modern Be-10-based erosion rate of,similar to 0.013 mm yr(-1). Our simulations also flag the two highest Be-10-derived erosion rates at 11.8 kyr BP related to nonuniform catchment erosion. These changes in erosion rates and processes during the AHP may reflect a strong increase in precipitation, runoff, and erosivity at the arid-to-humid transition either at 15 or similar to 12 cal. kyr BP, before the landscape stabilized again, possibly due to increased soil production and denser vegetation.
The subduction of bathymetric anomalies at convergent margins can profoundly affect subduction dynamics, magmatism, and the structural and geomorphic evolution of the overriding plate. The Northern Patagonian Icefield (NPI) is located east of the Chile Triple Junction at similar to 47 degrees S, where the Chile Rise spreading center collides with South America. This region is characterized by an abrupt increase in summit elevations and relief that has been controversially debated in the context of geodynamic versus glacial erosion effects on topography. Here we present geomorphic, thermochronological, and structural data that document neotectonic activity along hitherto unrecognized faults along the flanks of the NPI. New apatite (U-Th)/He bedrock cooling ages suggest faulting since 2-3 Ma. We infer the northward translation of an similar to 140 km long fore-arc sliver-the NPI block-results from enhanced partitioning of oblique plate convergence due to the closely spaced collision of three successive segments of the Chile Rise. In this model, greater uplift occurs in the hanging wall of the Exploradores thrust at the northern leading edge of the NPI block, whereas the Cachet and Liquine-Ofqui dextral faults decouple the NPI block along its eastern and western flanks, respectively. Localized extension possibly occurs at its southern trailing edge along normal faults associated with margin-parallel extension, tectonic subsidence, and lower elevations along the Andean crest line. Our neotectonic model provides a novel explanation for the abrupt topographic variations inland of the Chile Triple Junction and emphasizes the fundamental effects of local tectonics on exhumation and topographic patterns in this glaciated landscape.
Abstract. The Sea of Marmara, in northwestern Turkey, is a transition zone where the dextral North Anatolian Fault zone (NAFZ) propagates westward from the Anatolian Plate to the Aegean Sea Plate. The area is of interest in the context of seismic hazard of Istanbul, a metropolitan area with about 15 million inhabitants. Geophysical observations indicate that the crust is heterogeneous beneath the Marmara basin, but a detailed characterization of the crustal heterogeneities is still missing. To assess if and how crustal heterogeneities are related to the NAFZ segmentation below the Sea of Marmara, we develop new crustal-scale 3-D density models which integrate geological and seismological data and that are additionally constrained by 3-D gravity modeling. For the latter, we use two different gravity datasets including global satellite data and local marine gravity observation. Considering the two different datasets and the general non-uniqueness in potential field modeling, we suggest three possible “end-member” solutions that are all consistent with the observed gravity field and illustrate the spectrum of possible solutions. These models indicate that the observed gravitational anomalies originate from significant density heterogeneities within the crust. Two layers of sediments, one syn-kinematic and one pre-kinematic with respect to the Sea of Marmara formation are underlain by a heterogeneous crystalline crust. A felsic upper crystalline crust (average density of 2720 kgm⁻³) and an intermediate to mafic lower crystalline crust (average density of 2890 kgm⁻³) appear to be cross-cut by two large, dome-shaped mafic highdensity bodies (density of 2890 to 3150 kgm⁻³) of considerable thickness above a rather uniform lithospheric mantle (3300 kgm⁻³). The spatial correlation between two major bends of the main Marmara fault and the location of the highdensity bodies suggests that the distribution of lithological heterogeneities within the crust controls the rheological behavior along the NAFZ and, consequently, maybe influences fault segmentation and thus the seismic hazard assessment in the region.
Abstract. The Sea of Marmara, in northwestern Turkey, is a transition zone where the dextral North Anatolian Fault zone (NAFZ) propagates westward from the Anatolian Plate to the Aegean Sea Plate. The area is of interest in the context of seismic hazard of Istanbul, a metropolitan area with about 15 million inhabitants. Geophysical observations indicate that the crust is heterogeneous beneath the Marmara basin, but a detailed characterization of the crustal heterogeneities is still missing. To assess if and how crustal heterogeneities are related to the NAFZ segmentation below the Sea of Marmara, we develop new crustal-scale 3-D density models which integrate geological and seismological data and that are additionally constrained by 3-D gravity modeling. For the latter, we use two different gravity datasets including global satellite data and local marine gravity observation. Considering the two different datasets and the general non-uniqueness in potential field modeling, we suggest three possible “end-member” solutions that are all consistent with the observed gravity field and illustrate the spectrum of possible solutions. These models indicate that the observed gravitational anomalies originate from significant density heterogeneities within the crust. Two layers of sediments, one syn-kinematic and one pre-kinematic with respect to the Sea of Marmara formation are underlain by a heterogeneous crystalline crust. A felsic upper crystalline crust (average density of 2720 kgm⁻³) and an intermediate to mafic lower crystalline crust (average density of 2890 kgm⁻³) appear to be cross-cut by two large, dome-shaped mafic highdensity bodies (density of 2890 to 3150 kgm⁻³) of considerable thickness above a rather uniform lithospheric mantle (3300 kgm⁻³). The spatial correlation between two major bends of the main Marmara fault and the location of the highdensity bodies suggests that the distribution of lithological heterogeneities within the crust controls the rheological behavior along the NAFZ and, consequently, maybe influences fault segmentation and thus the seismic hazard assessment in the region.
Abstract. The aim of this study is to investigate the shallow thermal field differences for two differently aged passive continental margins by analyzing regional variations in geothermal gradient and exploring the controlling factors for these variations. Hence, we analyzed two previously published 3-D conductive and lithospheric-scale thermal models of the Southwest African and the Norwegian passive margins. These 3-D models differentiate various sedimentary, crustal, and mantle units and integrate different geophysical data such as seismic observations and the gravity field. We extracted the temperature–depth distributions in 1 km intervals down to 6 km below the upper thermal boundary condition. The geothermal gradient was then calculated for these intervals between the upper thermal boundary condition and the respective depth levels (1, 2, 3, 4, 5, and 6 km below the upper thermal boundary condition). According to our results, the geothermal gradient decreases with increasing depth and shows varying lateral trends and values for these two different margins. We compare the 3-D geological structural models and the geothermal gradient variations for both thermal models and show how radiogenic heat production, sediment insulating effect, and thermal lithosphere–asthenosphere boundary (LAB) depth influence the shallow thermal field pattern. The results indicate an ongoing process of oceanic mantle cooling at the young Norwegian margin compared with the old SW African passive margin that seems to be thermally equilibrated in the present day.
Abstract. The aim of this study is to investigate the shallow thermal field differences for two differently aged passive continental margins by analyzing regional variations in geothermal gradient and exploring the controlling factors for these variations. Hence, we analyzed two previously published 3-D conductive and lithospheric-scale thermal models of the Southwest African and the Norwegian passive margins. These 3-D models differentiate various sedimentary, crustal, and mantle units and integrate different geophysical data such as seismic observations and the gravity field. We extracted the temperature–depth distributions in 1 km intervals down to 6 km below the upper thermal boundary condition. The geothermal gradient was then calculated for these intervals between the upper thermal boundary condition and the respective depth levels (1, 2, 3, 4, 5, and 6 km below the upper thermal boundary condition). According to our results, the geothermal gradient decreases with increasing depth and shows varying lateral trends and values for these two different margins. We compare the 3-D geological structural models and the geothermal gradient variations for both thermal models and show how radiogenic heat production, sediment insulating effect, and thermal lithosphere–asthenosphere boundary (LAB) depth influence the shallow thermal field pattern. The results indicate an ongoing process of oceanic mantle cooling at the young Norwegian margin compared with the old SW African passive margin that seems to be thermally equilibrated in the present day.
The Victoria microplate between the Eastern and Western Branches of the East African Rift System is one of the largest continental microplates on Earth. In striking contrast to its neighboring plates, Victoria rotates counterclockwise with respect to Nubia. The underlying cause of this distinctive rotation has remained elusive so far. Using 3D numerical models, we investigate the role of pre-existing lithospheric heterogeneities in continental microplate rotation. We find that Victoria's rotation is primarily controlled by the distribution of rheologically stronger zones that transmit the drag of the major plates to the microplate and of the mechanically weaker mobile belts surrounding Victoria that facilitate rotation. Our models reproduce Victoria's GPS-derived counterclockwise rotation as well as key complexities of the regional tectonic stress field. These results reconcile competing ideas on the opening of the rift system by highlighting differences in orientation of the far-field divergence, local extension, and the minimum horizontal stress. One of the largest continental microplates on Earth is situated in the center of the East African Rift System, and oddly, the Victoria microplate rotates counterclockwise with respect to the neighboring African tectonic plate. Here, the authors' modelling results suggest that Victoria microplate rotation is caused by edge-driven lithospheric processes related to the specific geometry of rheologically weak and strong regions.
The Vicuna Pampa volcanic complex, at the SE edge of the arid Puna Plateau of the Central Andes, records the interplay between volcanic construction and degra-dational processes. The low-sloping Vicuna Pampa volcanic complex, with a 1200-m-deep, southeastward-opening depression, was previously interpreted as a collapse caldera based on morphological considerations. However, characteristic features associated with collapse calderas do not exist, and close inspection instead suggests that the Vicuna Pampa volcanic complex is a strongly eroded, broad, massif-type composite volcano of mainly basaltic to trachyandesitic composition. Construction of the Vicuna Pampa volcanic complex occurred during two distinct cycles separated by the development of the depression. The first and main cycle took place at ca. 12 Ma and was dominated by lava flows and subordinate scoria cones and domes. The second cycle, possibly late Miocene in age, affected the SW portion of the depression with the emplacement of domes. We interpret the central depression as the result of a possible sector collapse and subsequent intense fluvial erosion during middle to late Miocene time, facilitated by faulting, steepened topography, and wetter climate conditions compared to today. We estimate that similar to 65% of the initial edifice of similar to 240 km(3) was degraded. The efficiency of degradation processes for removing mass from the Vicuna Pampa volcanic complex is surprising, considering that today the region is arid, and the stream channels within the complex are predominantly transport limited, forming a series of coalesced, aggraded alluvial fans and eolian infill. Hence, the Vicuna Pampa volcanic complex records the effects of past degradation efficiency that differs substantially from that of today.
The northwest Argentine Andes constitute a premier natural laboratory to assess the complex interactions between isolated uplifts, orographic precipitation gradients, and related erosion and sedimentation patterns. Here we present new stratigraphic observations and age information from intermontane basin sediments to elucidate the Neogene to Quaternary shortening history and associated sediment dynamics of the broken Salta foreland. This part of the Andean orogen, which comprises an array of basement-cored range uplifts, is located at similar to 25 degrees S and lies to the east of the arid intraorogenic Altiplano/Puna plateau. In the Salta foreland, spatially and temporally disparate range uplift along steeply dipping inherited faults has resulted in foreland compartmentalization with steep basin-tobasin precipitation gradients. Sediment architecture and facies associations record a three-phase (similar to 10, similar to 5, and <2 Ma), east directed, yet unsystematic evolution of shortening, foreland fragmentation, and ensuing changes in precipitation and sediment transport. The provenance signatures of these deposits reflect the trapping of sediments in the intermontane basins of the Andean hinterland, as well as the evolution of a severed fluvial network. Present-day moisture supply to the hinterland is determined by range relief and basin elevation. The conspiring effects of range uplift and low rainfall help the entrapment and long-term storage of sediments, ultimately raising basin elevation in the hinterland, which may amplify aridification in the orogen interior.
[1] We document late Pleistocene - Holocene aggradation and incision processes at the mountain front of the Qilian Shan, an active intracontinental fold-and-thrust belt accommodating a significant portion of the India-Asia convergence. The Shiyou River cuts through a NNE vergent fault propagation fold with Miocene red beds in the core and Pliocene - Quaternary growth strata on the northern forelimb. South of the anticline, Miocene strata dip 20 degrees SSW, suggesting a similar orientation for the basal decollement. After aggradation of an similar to 150-m-thick, late Pleistocene valley fill, the Shiyou River formed three terraces. The highest terrace, located 170 m above the river, constitutes the top of the fill. The other terraces are fill cut terraces: their treads are located 130 - 105 m and 37 m above the river, respectively. The 10 Be exposure dating of the terraces suggests that river incision accelerated from 0.8 +/- 0.2 mm yr(-1) to similar to 10 mm yr(-1) at 10 - 15 kyr. We interpret fast Holocene river incision as largely unrelated to tectonic forcing. The late Pleistocene incision rate of 0.8 +/- 0.2 mm yr(-1) places an upper limit of 2.2 +/- 0.5 mm yr(-1) on the horizontal shortening rate, assuming that incision is solely caused by rock uplift above a decollement dipping 20 degrees. However, the actual shortening rate may lie between similar to 2.2 mm yr(-1) and zero because deformation of the terraces and the valley fill cannot be unequivocally demonstrated. Our estimate is consistent with the bulk shortening rate of similar to 5 - 10 mm yr(-1) across several faults in NE Tibet derived from neotectonic and GPS data, although in case of the Shiyou River, Holocene deformation is barely discernible owing to intense climate-induced river incision.
We derive a slip rate for a thrust at the central Qilian Shan mountain front by combining structural investigations, satellite imagery, topographic profiling, luminescence dating, and Be-10 exposure dating. The seismically active Zhangye thrust transects late Pleistocene alluvial fan deposits and forms a prominent north facing scarp. The fault consists of two segments that differ in orientation, scarp height, and age. A series of loess-covered terraces records the uplift history of the western thrust segment. Loess accumulation on all terraces started at 8.5 +/- 1.5 kyr and postdates terrace formation. Gravels from the highest terrace yielded a Be-10 exposure age of 90 +/- 11 kyr, which dates the onset of faulting. With a displacement of 55-60 m derived from fault scarp profiles, this yields a vertical slip rate of 0.64 +/- 0.08 mm yr(-1). Along the eastern thrust segment, three Be-10 ages from the uplifted alluvial fan constrain that faulting started at similar to31 +/- 5 kyr. Together with a displacement of 25-30 m this leads to a vertical faulting rate of 0.88 +/- 0.16 mm yr(-1). A dip estimate of 40degrees to 60degrees for the fault plane combined with lower and upper limits of similar to0.6 and similar to0.9 mm yr(-1) for the vertical slip rate gives minimum and maximum horizontal shortening rates of 0.4 and 1.1 mm yr(-1) across the Zhangye thrust. Our results are consistent with geologic and GPS constraints, which suggest that NNE directed shortening across the northeastern Tibetan Plateau is distributed on several active faults with a total shortening rate of 4 to 10 mm yr(-1)
A fault scaling law suggests that, over eight orders of magnitude, fault length L is linearly related to maximum displacement D. Individual faults may therefore retain a constant ratio of D/L as they grow. If erosion is minor compared with tectonic uplift, the length and along-strike relief of young mountain ranges should thus reflect fault growth. Topographic profiles along the crests of mountain ranges in the actively deforming foreland of north-east Tibet exhibit a characteristic shape with maximum height near their centre and decreasing elevation toward the tips. We interpret the along-strike relief of these ranges to reflect the slip distribution on high-angle reverse faults. A geometric model illustrates that the lateral propagation rate of such mountain ranges may be deciphered if their length- to-height ratio has remained constant. As an application of the model, we reconstruct the growth of the Heli Shan using a long-term uplift rate of similar to1.3 mm yr(-1) derived from Ne-21 and Be-10 exposure dating
Basement-cored uplift provinces are often characterized by high-angle reverse faulting along preexisting crustal heterogeneities, which may greatly affect the mechanics of deformation and the coupling between erosion and orogenic structure. Herein we construct a coupled deformation-erosion model to understand the mechanics and erosion of mountain belts in which the spatial distribution of deformation is largely influenced by the presence of preexisting high-angle faults. In this case, deformation is accommodated along, and topography is built above, these structures. This topographic loading leads to increasing lithostatic stresses beneath these regions. As a result, active deformation may migrate to frictionally stronger structures in adjacent regions where lithostatic loading is lower. The migration of deformation to such nearby structures depends on the Hubbert-Rubey pore fluid pressure ratio of the crust (lambda), the orientations of the frictionally weaker and stronger preexisting faults (beta(1) and beta(2), respectively), the friction coefficients (mu(b1) and mu(b2)) and Hubbert-Rubey fluid-pressure ratios (lambda(b1) and lambda(b2)) of these faults, and the surface slope of the topography above the frictionally weaker structure (alpha), assuming zero surface slope above the frictionally stronger structure. In general, we found that for a given alpha and beta(1), as mu(b1) increases lambda=lambda(b1)=lambda(b2) increases, and beta(2) decreases, the value of mu(b2) required to force deformation to migrate increases. However, as erosional processes lead to decreasing surface slopes, deformation will be inhibited from migrating to frictionally stronger structures in adjacent regions. Our model results may help to explain some aspects of the deformation observed and the possible correlation between precipitation and the migration of deformation within these tectonic provinces
Intramontane basins may act as important sediment storage areas, serve as recorders of the history of deformation, record syntectonic deposition, and document the evolution of climatic conditions during deposition. We document the timing, cyclicity, and processes that led to the filling and reexcavation of the intramontane Quebrada del Toro basin in NW Argentina. Geomorphic and geologic observations indicate that the basin was filled with sediment that has been subsequently excavated at least two times in the last similar to 8 m.y. The last filling and excavation cycle occurred within the last 0.98 m.y. and has led to the deposition and removal of similar to 61.4 km(3) of material from the basin, leading to a basin-wide averaged minimum denudation rate of 0.16 mm/yr. Aggradation within the basin took place due to channel steepening of the downstream fluvial system that connects the intramontane basin to the foreland. This portion of the fluvial system is actively incising through an uplifting bedrock zone. We use observations within the Toro to test a quasiphysically based model of channel aggradation behind a rising base level that rises due to downstream channel steepening. Our work shows that the bedrock incision rate constant required to reproduce conditions observed within the Toro basin is consistent with values measured independently in similar rock types. Therefore, in intramontane basins that experience similar processes of filling and evacuation, this model may be used to assess the relative importance of tectonic rock uplift, bedrock resistance to fluvial incision, and climate in determining the geomorphic and sedimentologic history of these basins
[1] Orogenic structure appears to be partially controlled by the addition to and removal of material from the mountain belt by tectonic accretion and geomorphic erosion, respectively. We developed a coupled erosion-deformation model for orogenic wedges that are in erosional steady state and deform at their Coulomb failure limit. Erosional steady state is reached when all material introduced into the wedge is removed by erosion that is limited by the rate at which rivers erode through bedrock. We found that the ultimate form of a wedge is controlled by the wedge mechanical properties, sole-out depth of the basal decollement, erosional exponents, basin geometry, and the ratio of the added material flux to the erosional constant. As this latter ratio is increased, wedge width and surface slopes increase. We applied these models to the Taiwan and Himalayan orogenic wedges and found that despite a higher flux of material entering the former, the inferred ratio was larger for the latter. Calculated values for the erodibility of each wedge showed at least an order of magnitude lower value for the Himalaya relative to Taiwan. These values are consistent with the lower precipitation regime in the Himalaya relative to Taiwan and the exposure of crystalline rocks within the Himalayan orogenic wedge. Independently determined rock erodibility estimates are consistent with the accretionary wedge sediments and metasediments and the crystalline and high-grade metamorphic rocks exposed within Taiwan and the Himalaya, respectively. Therefore differences in rock type and climate apparently lead to key differences in the erosion and hence orogenic structure of these two mountain belts
[1] The development of topography within and erosional removal of material from an orogen exerts a primary control on its structure. We develop a model that describes the temporal development of a frontally accreting, critically growing Coulomb wedge whose topography is largely limited by bedrock fluvial incision. We present general results for arbitrary initial critical wedge geometries and investigate the temporal development of a critical wedge with no initial topography. Increasing rock erodibility and/or precipitation, decreasing mass flux accreting to the wedge front, increasing wedge sole-out depth, decreasing wedge and basal decollement overpressure, and increasing basal decollement friction lead to narrow wedges. Large power law exponent values cause the wedge geometry to quickly reach a condition in which all material accreted to the front of the wedge is removed by erosion. We apply our model to the Aconcagua fold-and-thrust belt in the central Andes of Argentina where wedge development over time is well constrained. We solve for the erosional coefficient K that is required to recreate the field-constrained wedge growth history, and these values are within the range of independently determined values in analogous rock types. Using qualitative observations of rock erodibilities within the wedge, we speculate that power law exponents of 1/3 less than or equal to m less than or equal to 0.4 and 2/3 less than or equal to n less than or equal to 1 characterize the erosional growth of the Aconcagua fold-and-thrust belt. This general model may be used to understand the development of mountain belts where orogenic wedges grow as they deform at their Coulomb failure limit
Synorogenic extension has been recognized as an integral structural constituent of mountain belts and high-elevation plateaus during their evolution. In the Himalaya, both orogen-parallel and orogen-normal extension has been recognized. However, the underlying driving forces for extension and their timing are still a matter of debate. Here we present new fault kinematic data based on systematic measurements of hundreds of outcrop-scale brittle fault planes in the NW Indian Himalaya. This new data set, as well as field observations including crosscutting relationships, mineral fibers on fault planes, and correlations with deformation structures in lake sediments, allows us to distinguish different deformation styles. The overall strain pattern derived from our data reflects the large regional contractional deformation pattern very well but also reveals significant extensional deformation in a region, which is dominated by shortening. In total, we were able to identify six deformation styles, most of which are temporally and spatially linked, representing protracted shortening. Our observations also furnish the basis for a detailed overview of the younger deformation history in the NW Himalaya, which has been characterized by extension overprinting previously generated structures related to shortening. The four dominant deformation styles are (1) shortening parallel to the regional convergence direction; (2) arc-normal extension; (3) arc-parallel extension; and finally, (4) E-W extension. This is the first data set where a succession of both arc-normal and E-W extension has been documented in the Himalaya. Importantly, our observations help differentiate E-W extension triggered by processes within the Tibetan Plateau from arc-parallel and arc-normal extension originating from the curvature of the Himalayan orogen.
Explaining the presence of normal faults in overall compressive settings is a challenging problem in understanding the tectonics of active mountain belts. The Himalayan-Tibetan orogenic system is an excellent setting to approach this problem because it preserves one of the most dramatic records of long-term, contemporaneous shortening and extension. Over the past decades, several studies have described extensional features, not only in the Tibetan Plateau, but also in the Himalaya. For a long time, the favored model explained the function of the Southern Tibetan detachment system, a major fault zone in the Himalaya, as a decoupling horizon between the regime of crustal shortening forming the Himalayan wedge to the south and the extensional regime of the Tibetan Plateau to the north. However, in recent years, increasing evidence has shown that N-S-trending normal faults in the Central Himalaya crosscut not only the Southern Tibetan detachment system, but also the Main Central thrust. Here, we present new structural data and geologic evidence collected within the NW Indian Himalaya and combine them with previously published seismicity data sets in order to document pervasive E-W extension accommodated along N-S-trending faults extending as far south as the footwall of the Main Central thrust. We conducted a kinematic analysis of fault striations on brittle faults, documented and mapped fault scarps in Quaternary sedimentary deposits using satellite imagery, and made field observations in the Greater Sutlej region (Spiti, Lahul, Kinnaur) and the Garhwal Himalaya. Studies of extensional features within the regionally NW- SE-trending NW Indian Himalaya provide the advantage that arc-parallel and E-W extension can be separated, in contrast to the Central Himalaya. Therefore, our observations of E-W extension in the Indian NW Himalaya are well suited to test the applicability of current tectonic models for the whole Himalaya. We favor the interpretation of E-W extension in the NW Indian Himalaya as a propagation of extension driven by collapse of the Tibetan Plateau.
Stable isotope altimetry is a useful tool for estimating paleoelevation in sedimentary records. Yet questions remain regarding how source moisture, climate, and local topography can influence these estimates. Here we present stable isotope altimetry results on late Quaternary pedogenic carbonates of known elevation on both flanks of the Andean orogen at 33 degrees S. We measured delta O-18 values of pedogenic carbonates and river water samples from small drainages at regular elevation increments within the Rio Aconcagua (Chile) and Rio Mendoza (Argentina) catchments. The delta O-18 values of river waters correlate well with elevation and show similar isotopic gradients between the Chilean (-3.7 parts per thousand/km) and Argentine (-4.8 parts per thousand/km) sides of the range. Uncertainties associated with scatter in the river water data and assumptions about the temperature of carbonate formation indicate that elevation estimates have 1 sigma errors of 350-450 m. We estimate the isotopic composition of soil water from pedogenic carbonates on both sides of the range by assuming mean annual temperatures based the modern temperature lapse rate from meteorological station data. Combined, our data show that stable isotope altimetry produces reasonable estimates of modern elevation, with the majority of our samples (60%) within the 1 sigma uncertainties and 77% within 2 sigma.
We present a new three-dimensional density model of the Central Andes characterizing the structure and composition of the lithosphere together with a geodynamic simulation subjected to continental intraplate shortening. The principal aim of this study is to assess the link between heterogeneities in the lithosphere and different deformation patterns and styles along the orogen-foreland system of the Central Andes. First, we performed a 3D integration of new geological and geophysical data with previous models through forward modelling of Bouguer anomalies. Subsequently, a geodynamic model was set-up and parametrized from the previously obtained 3D structure and composition. We do not find a unambigous correlation between the resulting density configuration and terrane boundaries proposed by other authors. Our models reproduce the observed Bouguer anomaly and deformation patterns in the foreland. We find that thin-skinned deformation in the Subandean fold-and thrust belt is controlled by a thick sedimentary layer and coeval underthrusting of thin crust of the foreland beneath the thick crust of the Andean Plateau. In the adjacent thick-skinned deformation province of the inverted Cretaceous extensional Santa Barbara System sedimentary strata are much thinner and crustal thickness transitions from greater values in the Andean to a more reduced thickness in the foreland. Our results show that deformation processes occur where the highest gradients of lithospheric strength are present between the orogen and the foreland, thus suggesting a spatial correlation between deformation and lithospheric strength.
Constraining the pre-Neogene history of the Puna plateau is crucial for establishing the initial conditions that attended the early stage evolution of the southern extent of the Andean plateau. We apply high-to low-temperature thermochronology data from plutonic rocks in northwestern Argentina to quantify the Paleozoic, Mesozoic and early Tertiary cooling history of the Andean crust. U-Pb crystallization ages of zircons indicate that pluton intrusion occurred during the early mid-Ordovician (490-470 Ma) and the late Jurassic (160-150 Ma). Lower-temperature cooling histories from Ar-40/Ar-39 analyses of K-feldspar vary substantially. Basement rocks underlying the western Puna resided at temperatures below 200 degrees C (<6 km depth) since the Devonian (similar to 400 Ma). In contrast, basement rocks underlying the southeastern Puna were hotter (similar to 200-300 degrees C) throughout the Paleozoic and Jurassic and cooled to temperatures of <200 degrees C by similar to 120 Ma. The southeastern Puna basement records a rapid cooling phase coeval with active extension of the Cretaceous Salta rift at similar to 160-100 Ma that we associate with tectonic faulting and lithospheric thinning. The northeastern Puna experienced protracted cooling until the late Cretaceous with temperatures <200 degrees C during the Paleocene. Higher cooling rates between 78 and 55 Ma are associated with thermal subsidence during the postrift stage of the Salta rift and/or shortening-related flexural subsidence. Accelerated cooling and deformation during the Eocene was focused within a narrow zone along the eastern Puna/Eastern Cordillera transition that coincides with Paleozoic/Mesozoic structural and thermal boundaries. Our results constrain regional erosion-induced cooling throughout the Cenozoic to have been less than similar to 150 degrees C, which implies total Cenozoic denudation of <6-4 km.
The segmentation of major fault systems in subduction zones controls earthquake magnitude and location, but the causes for the existence of segment boundaries and the relationships between long-term deformation and the extent of earthquake rupture, are poorly understood. We compare permanent and seismic-cycle deformation patterns along the rupture zone of the 2010 Maule earthquake (M8.8), which ruptured 500 km of the Chile subduction margin. We analyzed the morphology of MIS-5 marine terraces using LiDAR topography and established their chronology and coeval origin with twelve luminescence ages, stratigraphy and geomorphic correlation, obtaining a virtually continuous distribution of uplift rates along the entire rupture zone. The mean uplift rate for these terraces is 0.5 m/ka. This value is exceeded in three areas, which have experienced rapid emergence of up to 1.6 m/ka; they are located at the northern, central, and southern sectors of the rupture zone, referred to as Topocalma, Carranza and Arauco, respectively. The three sectors correlate with boundaries of eight great earthquakes dating back to 1730. The Topocalma and Arauco sectors, located at the boundaries of the 2010 rupture, consist of broad zones of crustal warping with wavelengths of 60 and 90 km, respectively. These two regions coincide with the axes of oroclinal bending of the entire Andean margin and correlate with changes in curvature of the plate interface. Rapid uplift at Carranza, in turn, is of shorter wavelength and associated with footwall flexure of three crustal-scale normal faults. The uplift rate at Carranza is inversely correlated with plate coupling as well as with coseismic slip, suggesting permanent deformation may accumulate interseismically. We propose that the zones of upwarping at Arauco and Topocalma reflect changes in frictional properties of the megathrust resulting in barriers to the propagation of great earthquakes. Slip during the 1960 (M9.5) and 2010 events overlapped with the similar to 90-km-long zone of rapid uplift at Arauco; similarly, slip in 2010 and 1906 extended across the similar to 60-km-long section of the megathrust at Topocalma, but this area was completely breached by the 1730 (M similar to 9) event, which propagated southward until Carranza. Both Arauco and Topocalma show evidence of sustained rapid uplift since at least the middle Pleistocene. These two sectors might thus constitute discrete seismotectonic boundaries restraining most, but not all great earthquake ruptures. Based on our observations, such barriers might be breached during multi-segment super-cycle events. (C) 2015 Elsevier Ltd. All rights reserved.
The morphology of marine and lacustrine terraces has been largely used to measure past sea- and lake-level positions and estimate vertical deformation in a wealth of studies focused on climate and tectonic processes. To obtain accurate morphometric assessments of terrace morphology we present TerraceM-2, an improved version of our MatlabR (R) graphic-user interface that provides new methodologies for morphometric analyses as well as landscape evolution and fault-dislocation modeling. The new version includes novel routines to map the elevation and spatial distribution of terraces, to model their formation and evolution, and to estimate fault-slip rates from terrace deformation patterns. TerraceM-2 has significantly improves its processing speed and mapping capabilities, and includes separate functions for developing customized workflows beyond the graphic-user interface. We illustrate these new mapping and modeling capabilities with three examples: mapping lacustrine shorelines in the Dead Sea to estimate deformation across the Dead Sea Fault, landscape evolution modeling to estimate a history of uplift rates in southern Peru, and dislocation modeling of deformed marine terraces in California. These examples also illustrate the need to use topographic data of different resolutions. The new modeling and mapping routines of TerraceM-2 highlight the advantages of an integrated joint mapping and modeling approach to improve the efficiency and precision of coastal terrace metrics in both marine and lacustrine environments.
High-resolution topographic data greatly facilitate the remote identification of geomorphic features, furnishing valuable information concerning surface processes and characterization of reference markers for quantifying tectonic deformation. Marine terraces have been used as long baseline geodetic markers of relative past sea-level positions, reflecting the interplay between vertical crustal movements and sea-level oscillations. Uplift rates may be determined from the terrace age and the elevation of its shoreline angle, a geomorphic feature that can be correlated with past sea-levels positions. A precise definition of the shoreline angle in time and space is essential to obtain reliable uplift rates with coherent spatial correlation. To improve our ability to rapidly assess and map shoreline angles at regional and local scales, we have developed TerraceM, a MATLAB (R) graphical user interface that allows the shoreline angle and its associated error to be estimated using high-resolution topography. TerraceM uses topographic swath profiles oriented orthogonally to the terrace riser. Four functions are included to analyze the swath profiles and extract the shoreline angle, from both staircase sequences of multiple terraces and rough coasts characterized by eroded remnants of emerged terrace surfaces. The former are measured by outlining the paleocliffs and paieo-platforms and finding their intersection by extrapolating linear regressions, whereas the latter are assessed by automatically detecting peaks of sea-stack tops and back-projecting them to the modern sea cliff. In the absence of rigorous absolute age determinations of marine terraces, their geomorphic age may be estimated using previously published diffusion models. Postprocessing functions are included to obtain first-order statistics of shoreline-angle elevations and their spatial distribution. TerraceM has the ability to process series of profiles from several sites in an efficient and structured workflow. Results may be exported in Google Earth and ESRI shapefile formats. The precision and accuracy of the method have been estimated from a case study at Santa Cruz, California, by comparing TerraceM results with published field measurements. The repeatability was evaluated using multiple measurements made by inexperienced users. TerraceM will improve the efficiency and precision of estimating shoreline-angle elevations in wave-cut terraces in both marine and lacustrine environments.
Most of the deformation associated with the seismic cycle in subduction zones occurs offshore and has been therefore difficult to quantify with direct observations at millennial timescales. Here we study millennial deformation associated with an active splay-fault system in the Arauco Bay area off south central Chile. We describe hitherto unrecognized drowned shorelines using high-resolution multibeam bathymetry, geomorphic, sedimentologic, and paleontologic observations and quantify uplift rates using a Landscape Evolution Model. Along a margin-normal profile, uplift rates are 1.3m/ka near the edge of the continental shelf, 1.5m/ka at the emerged Santa Maria Island, -0.1m/ka at the center of the Arauco Bay, and 0.3m/ka in the mainland. The bathymetry images a complex pattern of folds and faults representing the surface expression of the crustal-scale Santa Maria splay-fault system. We modeled surface deformation using two different structural scenarios: deep-reaching normal faults and deep-reaching reverse faults with shallow extensional structures. Our preferred model comprises a blind reverse fault extending from 3km depth down to the plate interface at 16km that slips at a rate between 3.0 and 3.7m/ka. If all the splay-fault slip occurs during every great megathrust earthquake, with a recurrence of similar to 150-200years, the fault would slip similar to 0.5m per event, equivalent to a magnitude similar to 6.4 earthquake. However, if the splay-fault slips only with a megathrust earthquake every similar to 1000years, the fault would slip similar to 3.7m per event, equivalent to a magnitude similar to 7.5 earthquake.
The first step towards assessing hazards in seismically active regions involves mapping capable faults and estimating their recurrence times. While the mapping of active faults is commonly based on distinct geologic and geomorphic features evident at the surface, mapping blind seismogenic faults is complicated by the absence of on-fault diagnostic features. Here we investigated the Pichilemu Fault in coastal Chile, unknown until it generated a Mw 7.0 earthquake in 2010. The lack of evident surface faulting suggests activity along a partly-hidden blind fault. We used off-fault deformed marine terraces to estimate a fault-slip rate of 0.52 ± 0.04 m/ka, which, when integrated with satellite geodesy suggests a 2.12 ± 0.2 ka recurrence time for Mw~7.0 normal-faulting earthquakes. We propose that extension in the Pichilemu region is associated with stress changes during megathrust earthquakes and accommodated by sporadic slip during upper-plate earthquakes, which has implications for assessing the seismic potential of cryptic faults along convergent margins and elsewhere.
The first step towards assessing hazards in seismically active regions involves mapping capable faults and estimating their recurrence times. While the mapping of active faults is commonly based on distinct geologic and geomorphic features evident at the surface, mapping blind seismogenic faults is complicated by the absence of on-fault diagnostic features. Here we investigated the Pichilemu Fault in coastal Chile, unknown until it generated a Mw 7.0 earthquake in 2010. The lack of evident surface faulting suggests activity along a partly-hidden blind fault. We used off-fault deformed marine terraces to estimate a fault-slip rate of 0.52 ± 0.04 m/ka, which, when integrated with satellite geodesy suggests a 2.12 ± 0.2 ka recurrence time for Mw~7.0 normal-faulting earthquakes. We propose that extension in the Pichilemu region is associated with stress changes during megathrust earthquakes and accommodated by sporadic slip during upper-plate earthquakes, which has implications for assessing the seismic potential of cryptic faults along convergent margins and elsewhere.
Intraplate earthquakes pose a significant seismic hazard in densely populated rift systems like the Lower Rhine Graben in Central Europe. While the locations of most faults in this region are well known, constraints on their seismogenic potential and earthquake recurrence are limited. In particular, the Holocene deformation history of active faults remains enigmatic. In an exposure excavated across the Schafberg fault in the southwestern Lower Rhine Graben, south of Untermaubach, in the epicentral region of the 1756 Duren earthquake (M (L) 6.2), we mapped a complex deformation zone in Holocene fluvial sediments. We document evidence for at least one paleoearthquake that resulted in vertical surface displacement of 1.2 +/- 0.2 m. The most recent earthquake is constrained to have occurred after 815 AD, and we have modeled three possible earthquake scenarios constraining the timing of the latest event. Coseismic deformation is characterized by vertical offset of sedimentary contacts distributed over a 10-m-wide central damage zone. Faults were identified where they fracture and offset pebbles in the vertically displaced gravel layers and fracture orientation is consistent with the orientation of the Schafberg fault. This study provides the first constraint on the most recent surface-rupturing earthquake on the Schafberg fault. We cannot rule out that this fault acted as the source of the 1756 Duren earthquake. Our study emphasizes the importance of, and the need for, paleoseismic studies in this and other intracontinental regions, in particular on faults with subtle geomorphic expression that would not typically be recognized as being potentially seismically active. Our study documents textural features in unconsolidated sediment that formed in response to coseismic rupturing of the underlying bedrock fault. We suggest that these features, e.g., abundant oriented transgranular fractures in their context, should be added to the list of criteria used to identify a fault as potentially active. Such information would result in an increase of the number of potentially active faults that contribute to seismic hazards of intracontinental regions.
The Lower Rhine Graben (Central Europe) is a prime example of a seismically active low-strain rift zone characterized by pronounced anthropogenic and climatic overprint of structures, and long recurrence intervals of large earthquakes. These factors render the identification of active faults and surface ruptures difficult. We investigated two fault scarps in the Lower Rhine Graben, to decipher their structural character, offset and potential seismogenic origin. Both scarps were modified by anthropogenic activity. The Hemmerich site lies c. 20 km SW of Cologne, along the Erft Fault. The Untermaubach site lies SW of Duren, where the Schafberg Fault projects into the Rur River valley. At the Hemmerich site, geomorphic and geophysical data, as well as exploratory coring reveal evidence of repeated normal faulting. Geophysical analysis and palaeoseismological excavation at the Untermaubach site reveal a complex fault zone in Holocene gravels characterized by subtle gravel deformation. Differentiation of tectonic and fluvial features was only possible with trenching, because fault structures and grain sizes of the sediments were below the resolution of the geophysical data. Despite these issues, our investigation demonstrates that valuable insight into past earthquakes and seismogenic deformation in a low-strain environment can be revealed using a multidisciplinary approach.
Neighbouring faults can interact, potentially link up and grow, and consequently increase the seismic and related natural hazards in their vicinity. Despite evidence of Quaternary faulting, the kinematic relationships between the neighbouring Mosha Fasham Fault (MFF) and the North Tehran Thrust (NTT) and their temporal evolution in the Alborz mountains are not well understood. The ENE-striking NTT is a frontal thrust that delimits the Alborz mountains to the south with a 2000 m topographic front with respect to the proximal Tehran plain. However, no large instrumentally recorded earthquakes have been attributed to that fault. In contrast, the sigmoidally shaped MFF is a major strike-slip fault, located within the Alborz Mountains. Sinistral motion along the eastern part of the MFF is corroborated by microseismicity and fault kinematic analysis, which documents recent transtensional deformation associated with NNE-SSW oriented shortening. To better understand the activity of these faults on different timescales, we combined fault- kinematic analysis and geomorphic observations, to infer the kinematic history of these structures. Our fault kinematic study reveals an early dextral shear for the NTT and the central MFF, responsible for dextral strike-slip and oblique reverse faulting during NW-oriented shortening. This deformation regime was superseded by NE-oriented shortening, associated with sinistral-oblique thrusting along the NTT and the central-western MFF, sinistral strike-slip motion along subsidiary faults in the central MFF segment, and folding and tilting of Eocene to Miocene units in the MFF footwall. Continued thrusting along the NTT took place during the Quaternary. However, folding in the hanging wall and sinistral stream-offsets indicate a left-oblique component and Quaternary strike-slip reactivation of the eastern NTT- segment, close to its termination. This complex history of faulting under different stress directions has resulted in a composite landscape with inherited topographic signatures. Our study shows that the topography of the hanging wall of the NTT reflects a segmentation into sectors with semi-independent uplift histories. Areas of high topographic residuals and apparent high uplift underscore the fault kinematics. Combined, our data suggest an early mechanical linkage of the NTT and MFF fault systems during a former dextral transpressional stage, caused by NW-compression. During NE-oriented shortening, the NTT and MFF were reactivated and incorporated into a nascent transpressional duplex. The youngest manifestation of motion in this system is sinistral transtension. However, this deformation is not observed everywhere and has not yet resulted in topographic inversion.
The northern Tien Shan of Kyrgyzstan and Kazakhstan has been affected by a series of major earthquakes in the late 19th and early 20th centuries. To assess the significance of such a pulse of strain release in a continental interior, it is important to analyze and quantify strain release over multiple time scales. We have undertaken paleoseismological investigations at two geomorphically distinct sites (Panfilovkoe and Rot Front) near the Kyrgyz capital Bishkek. Although located near the historic epicenters, both sites were not affected by these earthquakes. Trenching was accompanied by dating stratigraphy and offset surfaces using luminescence, radiocarbon, and Be-10 terrestrial cosmogenic nuclide methods. At Rot Front, trenching of a small scarp did not reveal evidence for surface rupture during the last 5000 years. The scarp rather resembles an extensive debris-flow lobe. At Panfilovkoe, we estimate a Late Pleistocene minimum slip rate of 0.2 +/- 0.1 mm/a, averaged over at least two, probably three earthquake cycles. Dip-slip reverse motion along segmented, moderately steep faults resulted in hanging wall collapse scarps during different events. The most recent earthquake occurred around 3.6 +/- 1.3 kyr ago (1 sigma), with dip-slip offsets between 1.2 and 1.4 m. We calculate a probabilistic paleomagnitude to be between 6.7 and 7.2, which is in agreement with regional data from the Kyrgyz range. The morphotectonic signals in the northern Tien Shan are a prime example of deformation in a tectonically active intracontinental mountain belt and as such can help understand the longer-term coevolution of topography and seismogenic processes in similar structural settings worldwide.
The tectonically driven growth of mountains reflects the characteristics of the underlying fault systems and the applied tectonic forces. Over time, fault networks might be relatively static, but stress conditions could change and result in variations in fault slip orientation. Such a tectonic landscape would transition from a simple to a composite state: the topography of simple landscapes is correlated with a single set of tectonic boundary conditions, while composite landscapes contain inherited topography due to earlier deformation under different boundary conditions. We use fault interaction modeling to compare vertical displacement fields with topographic metrics to differentiate the two types of landscapes. By successively rotating the axis of maximum horizontal stress, we produce a suite of vertical displacement fields for comparison with real landscapes. We apply this model to a transpressional duplex in the south central Alborz Mountains of Iran, where NW oriented compression was superseded by neotectonic NE compression. The consistency between the modeled displacement field and real landforms indicates that the duplex topography is mostly compatible with the modern boundary conditions, but might include a small remnant from the earlier deformation phase. Our approach is applicable for various tectonic settings and represents an approach to identify the changing boundary conditions that produce composite landscapes. It may be particularly useful for identifying changes that occurred in regions where river profiles may no longer record a signal of the change or where the spatial pattern of uplift is complex.
Climatic changes are of major importance in landslide generation in the Argentine Andes. Increased humidity as a potential influential factor was inferred from the temporal clustering of landslide deposits during a period of significantly wetter climate, 30,000 years ago. A change in seasonality was tested by comparing past (inferred from annual-layered lake deposits, 30,000 years old) and modern (present-day observations) precipitation changes. Quantitative analysis of cross recurrence plots were developed to compare the influence of the El Nino/Southern Oscillation (ENSO) on present and past rainfall variations. This analysis has shown the stronger influence of NE trades in the location of landslide deposits in the intra-andean basin and valleys, what caused a higher contrast between summer and winter rainfall and an increasing of precipitation in La Nina years. This is believed to reduce thresholds for landslide generation in the arid to semiarid intra-andean basins and valleys.
Previous thermomechanical modeling studies indicated that variations in the temperature and strength of the crystalline crust might be responsible for the juxtaposition of domains with thin-skinned and thick-skinned crustal deformation along strike the foreland of the central Andes. However, there is no evidence supporting this hypothesis from data-integrative models. We aim to derive the density structure of the lithosphere by means of integrated 3-D density modeling, in order to provide a new basis for discussions of compositional variations within the crust and for future thermal and rheological modeling studies. Therefore, we utilize available geological and geophysical data to obtain a structural and density model of the uppermost 200km of the Earth. The derived model is consistent with the observed Bouguer gravity field. Our results indicate that the crystalline crust in northern Argentina can be represented by a lighter upper crust (2,800kg/m(3)) and a denser lower crust (3,100kg/m(3)). We find new evidence for high bulk crustal densities >3,000kg/m(3) in the northern Pampia terrane. These could originate from subducted Puncoviscana wackes or pelites that ponded to the base of the crystalline crust in the late Proterozoic or indicate increasing bulk content of mafic material. The precise composition of the northern foreland crust, whether mafic or felsic, has significant implications for further thermomechanical models and the rheological behavior of the lithosphere. A detailed sensitivity analysis of the input parameters indicates that the model results are robust with respect to the given uncertainties of the input data.
This work explores the control of fore-arc structure on segmentation of megathrust earthquake ruptures using coastal geomorphic markers. The Arauco-Nahuelbuta region at the south-central Chile margin constitutes an anomalous fore- arc sector in terms of topography, geology, and exhumation, located within the overlap between the Concepcion and Valdivia megathrust segments. This boundary, however, is only based on similar to 500 years of historical records. We integrate deformed marine terraces dated by cosmogenic nuclides, syntectonic sediments, published fission track data, seismic reflection profiles, and microseismicity to analyze this earthquake boundary over 10(2) -10(6) years. Rapid exhumation of Nahuelbuta's dome-like core started at 4 +/- 1.2 Ma, coeval with inversion of the adjacent Arauco basin resulting in emergence of the Arauco peninsula. Here, similarities between topography, spatiotemporal trends in fission track ages, Pliocene-Pleistocene growth strata, and folded marine terraces suggest that margin-parallel shortening has dominated since Pliocene time. This shortening likely results from translation of a fore-arc sliver or microplate, decoupled from South America by an intra-arc strike-slip fault. Microplate collision against a buttress leads to localized uplift at Arauco accrued by deep-seated reverse faults, as well as incipient oroclinal bending. The extent of the Valdivia segment, which ruptured last in 1960 with an M-w 9.5 event, equals the inferred microplate. We propose that mechanical homogeneity of the fore-arc microplate delimits the Valdivia segment and that a marked discontinuity in the continental basement at Arauco acts as an inhomogeneous barrier controlling nucleation and propagation of 1960-type ruptures. As microplate-related deformation occurs since the Pliocene, we propose that this earthquake boundary and the extent of the Valdivia segment are spatially stable seismotectonic features at million year scale.
Megathrust earthquakes are commonly accompanied by increased upper-plate seismicity and occasionally triggered fault slip. In Chile, crustal faults slipped during and after the 2010 Maule (M8.8) earthquake. We studied the El Yolki fault (EYOF), a transtensional structure midways the Maule rupture not triggered in 2010. We mapped a Holocene coastal plain using light detection and ranging, which did not reveal surface ruptures. However, the inner-edge and shoreline angles along the coastal plain as well as 4.3- to 4.0-ka intertidal sediments are back-tilted on the EYOF footwall block, documenting 10 m of vertical displacement. These deformed markers imply similar to 2-mm/year throw rate, and dislocation models a slip rate of 5.6 mm/year for the EYOF. In a 5-m-deep trench, the Holocene intertidal sediments onlap to five erosive steps, interpreted as staircase wave-cut landforms formed by discrete events of relative sea level drop. We tentatively associated these steps with coseismic uplift during EYOF earthquakes between 4.3 and 4.0 ka. The Maule earthquake rupture may be subdivided into three subsegments based on coseismic slip and gravity anomalies. Coulomb stress transfer models predict neutral stress changes at the EYOF during the Maule earthquake but positive changes for a synthetic slip distribution at the central subsegment. If EYOF earthquakes were triggered by megathrust events, their slip distribution was probably focused in the central subsegment. Our study highlights the millennial variability of crustal faulting and the megathrust earthquake cycle in Chile, with global implications for assessing the hazards posed by hidden but potentially seismogenic coastal faults along subduction zones.
Central Anatolia is a low-relief, high-elevation region where decadal-scale deformation rates estimated from space geodesy suggest low strain rates within a stiff microplate. However, numerous Quaternary faults have been mapped within this low-strain region and estimating their slip rate and seismic potential is important for hazard assessments in an area of increasing infrastructural development. Here we focus on the Sultanhani Fault (SF), which constitutes an integral part of the Eskisehir-Cihanbeyli Fault System, and use deformed maximum highstand shorelines of palaeo-lake Konya to estimate tectonic slip rates at millennial scale. Some of these shorelines were previously interpreted as fault scarps, but we provide conclusive evidence for their erosional origin. We found that shoreline-angle elevations estimated from differential GPS profiles record vertical displacements of 10.2 m across the SF. New radiocarbon ages of lacustrine molluscs suggest 22.4 m of relative lake-level fall between 22.1 +/- 0.3 and 21.7 +/- 0.4 cal. kaBP, constraining the timing of abrupt abandonment of the highstand shoreline. Models of lithospheric rebound associated with regressions of the Tuz Golu and Konya palaeolakes predict only similar to 1 m of regional-scale uplift across the Konya Basin. Dislocation models of displaced shorelines suggest fault-slip rates of 1.5 and 1.8 mm yr(-1) for planar and listric fault geometries, respectively, providing reasonable results for the latter. We found fault scarps in the Nasuhpinar mudflat that likely represent the most recent ground-breaking rupture of the SF, with an average vertical displacement of 1.2 +/- 0.5 m estimated from 54 topographic profiles, equivalent to a M similar to 6.5-6.9 earthquake based on empirical scaling laws. If such events were characteristic during the ultimate 21 ka, a relatively short recurrence time of similar to 800-900 yr would be needed to account for the millennial slip rate. Alternatively, the fault scarp at Nasuhpinar might represent a larger earthquake requiring more frequent smaller events to account for the millennial rate. The relatively fast slip rate of the SF over the past 21 ka is unlikely to have persisted over longer timescales and might reflect spatiotemporal variations in deformation rates within kinematically-linked fault systems within Central Anatolia, or a transient perturbation to the local stress field or fault strength. Such perturbation might have been related to climatically controlled changes in surface and near-surface loads and by interactions among the different tectonic processes that have been proposed to drive the overall slow uplift and associated extension in the Central Anatolian Plateau.
Thick sedimentary fills in intermontane valleys are common in formerly glaciated mountain ranges but difficult to quantify. Yet knowledge of the fill thickness distribution could help to estimate sediment budgets of mountain belts and to decipher the role of stored material in modulating sediment flux from the orogen to the foreland. Here we present a new approach to estimate valley fill thickness and bedrock topography based on the geometric properties of a landscape using artificial neural networks. We test the potential of this approach following a four-tiered procedure. First, experiments with synthetic, idealized landscapes show that increasing variability in surface slopes requires successively more complex network configurations. Second, in experiments with artificially filled natural landscapes, we find that fill volumes can be estimated with an error below 20%. Third, in natural examples with valley fill surfaces that have steeply inclined slopes, such as the Unteraar and the Rhone Glaciers in the Swiss Alps, for example, the average deviation of cross-sectional area between the measured and the modeled valley fill is 26% and 27%, respectively. Finally, application of the method to the Rhone Valley, an overdeepened glacial valley in the Swiss Alps, yields a total estimated sediment volume of 9711km(3) and an average deviation of cross-sectional area between measurements and model estimates of 21.5%. Our new method allows for rapid assessment of sediment volumes in intermontane valleys while eliminating most of the subjectivity that is typically inherent in other methods where bedrock reconstructions are based on digital elevation models.
Unravelling the spatiotemporal evolution of the Cenozoic Andean (Altiplano-Puna) plateau has been one of the most intriguing problems of South American geology. Despite a number of investigations, the early deformation and uplift history of this area remained largely enigmatic. This paper analyses the Paleogene tectono-sedimentary history of the Casa Grande Basin, in the present-day transition zone between the northern sector of the Puna Plateau and the northern part of the Argentine Eastern Cordillera. Our detailed mapping of synsedimentary structures records the onset of regional contractional deformation during the middle Eocene, revealing reactivation of Cretaceous extensional structures and the development of doubly vergent thrusts. This is in agreement with records from other southern parts of the Puna Plateau and the Eastern Cordillera. These observations indicate the existence of an Eocene broken foreland setting within the region, characterized by low-lying compressional basins and ranges with spatially disparate sectors of deformation, which was subsequently subjected to regional uplift resulting in the attainment of present-day elevations during the Neogene.