Refine
Has Fulltext
- no (173)
Year of publication
Document Type
- Article (173) (remove)
Language
- English (173) (remove)
Is part of the Bibliography
- yes (173)
Keywords
- Himalaya (5)
- erosion (5)
- thermochronology (5)
- Andes (4)
- Pamir (4)
- exhumation (4)
- Alborz Mountains (3)
- Central Andes (3)
- TerraceM (3)
- cosmogenic nuclides (3)
The seismicity of the Kenya rift is characterized by high-frequency low-magnitude events concentrated along the rift axis. Its seismic character is typical for magmatically active continental rifts, where igneous material at a shallow depth causes extensive grid faulting and geothermal activity. Thermal overprinting and dike intrusion prohibit the buildup of large elastic strains, therefore prohibiting the generation of large-magnitude earthquakes. On 6 January 1928, the M-S 6.9 Subukia earthquake occurred on the Laikipia-Marmanet fault, the eastern rift-bounding structure of the central Kenya rift. It is the largest instrumentally recorded seismic event in the Kenya rift, standing in contrast to the current model of the rift's seismic character in which large earthquakes are not anticipated. Furthermore, the proximity of the ruptured fault and the rift axis is intriguing: The rift-bounding structure that ruptured in 1928 remains seismically active, capable of generating large-magnitude earthquakes, even though thermally weakened crust and better oriented structures are present along the rift axis nearby, prohibiting any significant buildup of elastic strain. We excavated the surface rupture of the 1928 Subukia earthquake to find evidence for preceding ground-rupturing earthquakes. We also made a total station survey of the site topography and mapped the site geology. We show that the Laikipia-Marmanet fault was repeatedly activated during the late Quaternary. We found evidence for six ground-rupturing earthquakes, including the 1928 earthquake. The topographic survey around the trench site revealed a degraded fault scarp of approximate to 7.5 m in height, offsetting a small debris slide. Using scarp-diffusion modeling, we estimated an uplift rate of U = 0.09-0.15 mm/yr, constraining the scarp age to 50-85 ka. Assuming an average fault dip of 55 degrees-75 degrees, the preferred uplift rate (0.15 mm/yr) accommodates approximately 10%-20% of the recent rate of extension (0.5 mm/yr) across the Kenya rift.
Surface uplift at the northern margin of the Central Anatolian Plateau (CAP) is integrally tied to the evolution of the Central Pontides (CP), between the North Anatolian Fault (NAF) and the Black Sea. Our regional morphometric and plate kinematic analyses reveal topographic anomalies, steep channel gradients, and local high relief areas as indicators of ongoing differential surface uplift, which is higher in the western CP compared to the eastern CP and fault-normal components of geodetic slip vectors and the character of tectonic activity of the NAF suggest that stress is accumulated in its broad restraining bend. Seismic reflection and structural field data show evidence for a deep structural detachment horizon responsible for the formation of an actively northward growing orogenic wedge with a positive flower-structure geometry across the CP and the NAF. Taken together, the tectonic, plate kinematic, and geomorphic observations imply that the NAF is the main driving mechanism for wedge tectonics and uplift in the CP. In addition, the NAF Zone defines the boundary between the extensional CAP and the contractional CP. The syntectonic deposits within inverted intermontane basins and deeply incised gorges suggest that the formation of relief, changes in sedimentary dynamics, and > 1 km fluvial incision resulted from accelerated uplift starting in the early Pliocene. The Central Pontides thus provide an example of an accretionary wedge with surface-breaking faults that play a critical role in mountain building processes, sedimentary basin development, and ensuing lateral growth of a continental plateau since the end of the Miocene.
We document Quaternary fluvial incision driven by fault-controlled surface deformation in the inverted intermontane Gökirmak Basin in the Central Pontide mountains along the northern margin of the Central Anatolian Plateau. In-situ-produced Be-10, Ne-21, and Cl-36 concentrations from gravel-covered fluvial terraces and pediment surfaces along the trunk stream of the basin (the Gökirmak River) yield model exposure ages ranging from 71ka to 34645ka and average fluvial incision rates over the past similar to 350ka of 0.280.01mm a(-1). Similarities between river incision rates and coastal uplift rates at the Black Sea coast suggest that regional uplift is responsible for the river incision. Model exposure ages of deformed pediment surfaces along tributaries of the trunk stream range from 605ka to 110 +/- 10ka, demonstrating that the thrust faults responsible for pediment deformation were active after those times and were likely active earlier as well as explaining the topographic relief of the region. Together, our data demonstrate cumulative incision that is linked to active internal shortening and uplift of similar to 0.3mm a(-1) in the Central Pontide orogenic wedge, which may ultimately contribute to the lateral growth of the northern Anatolian Plateau.
Emerged marine terraces and paleoshorelines along plate margins are prominent geomorphic markers that can be used to quantify the rates and patterns of crustal deformation. The northern margin of the Central Anatolian Plateau has been interpreted as an actively deforming orogenic wedge between the North Anatolian Fault and the Black Sea. Here we use uplifted marine terraces across principal faults on the Sinop Peninsula at the central northern side of the Pontide orogenic wedge to unravel patterns of Quaternary faulting and orogenic wedge behavior. We leveled the present-day elevations of paleoshorelines and dated marine terrace deposits using optically stimulated luminescence (OSL) to determine coastal uplift. The elevations of the paleoshorelines vary between 4 +/- 0.2 and 67 +/- 1.4 m above sea level and OSL ages suggest terrace formation episodes during interglacial periods at ca 125, 190, 400 and 570 ka, corresponding to marine isotopic stages (MIS) 5e, 7a, 11 and 15. Mean apparent vertical displacement rates (without eustatic correction) deduced from these terraces range between 0.02 and 0.18 mm/a, with intermittent faster rates of up to 0.26 mm/a. We obtained higher rates at the eastern and southern parts of the peninsula, toward the hinterland, indicating non-uniform uplift across the different morphotectonic segments of the peninsula. Our data are consistent with active on- and offshore faulting across the Sinop Peninsula. When integrated with regional tectonic observations, the faulting pattern reflects shortening distributed over a broad region of the northern margin of the Central Anatolian Plateau during the Quaternary.
Timing and magnitude of surface uplift are key to understanding the impact of crustal deformation and topographic growth on atmospheric circulation, environmental conditions, and surface processes. Uplift of the East African Plateau is linked to mantle processes, but paleoaltimetry data are too scarce to constrain plateau evolution and subsequent vertical motions associated with rifting. Here, we assess the paleotopographic implications of a beaked whale fossil (Ziphiidae) from the Turkana region of Kenya found 740 km inland from the present-day coastline of the Indian Ocean at an elevation of 620 m. The specimen is similar to 17 My old and represents the oldest derived beaked whale known, consistent with molecular estimates of the emergence of modern straptoothed whales (Mesoplodon). The whale traveled from the Indian Ocean inland along an eastward-directed drainage system controlled by the Cretaceous Anza Graben and was stranded slightly above sea level. Surface uplift from near sea level coincides with paleoclimatic change from a humid environment to highly variable and much drier conditions, which altered biotic communities and drove evolution in east Africa, including that of primates.
Cenozoic uplift of the East African Plateau has been associated with fundamental climatic and environmental changes in East Africa and adjacent regions. While this influence is widely accepted, the timing and the magnitude of plateau uplift have remained unclear. This uncertainty stems from the lack of datable, geomorphically meaningful reference horizons that could record surface uplift. Here, we document the existence of significant relief along the East African Plateau prior to rifting, as inferred from modeling the emplacement history of one of the longest terrestrial lava flows, the similar to 300-km-long Yatta phonolite flow in Kenya. This 13.5 Ma lava flow originated on the present-day eastern Kenya Rift flank, and utilized a riverbed that once routed runoff from the eastern rim of the plateau. Combining an empirical viscosity model with subsequent cooling and using the Yatta lava flow geometry and underlying paleotopography (slope angle), we found that the prerift slope was at least 0.2 degrees, suggesting that the lava flow originated at a minimum elevation of 1400 m. Hence, high paleotopography in the Kenya Rift region must have existed by at least 13.5 Ma. We infer from this that middle Miocene uplift occurred, which coincides with the two-step expansion of grasslands, as well as important radiation and speciation events in tropical Africa.
The Subandean fold and thrust belt of Bolivia constitutes the easternmost part of the Andean orogen that reflects thin-skinned shortening and eastward propagation of the Andean deformation front. The exact interplay of tectonics, climate, and erosion in the deposition of up to 7.5 km of late Cenozoic strata exposed in the Subandes remains unclear. To better constrain these relationships, we use four W-E industry seismic reflection profiles, eight new zircon U-Pb ages from Mio-Pliocene sedimentary strata, and cross-section balancing to evaluate the rates of thrust propagation, shortening, and deposition pinch-out migration. Eastward thrusting arrived in the Subandean belt at similar to 12.4 +/- 0.5 Ma and propagated rapidly toward the foreland unit approximately 6 Ma. This was followed by out-of- sequence deformation from ca. 4 to 2.1 Ma and by renewed eastward propagation thereafter. Our results show that the thrust-front propagation- and deposition pinch-out migration rates mimic the sediment accumulation rate. The rates of deposition pinchout migration and thrust propagation increased three- and two fold, respectively (8 mm/a; 3.3 mm/a) at 86 Ma. The three-fold increase in deposition pinch-out migration rate at this time is an indication of enhanced erosional efficiency in the hinterland, probably coupled with flexural rebound of the basin. Following the pulse of pinch-out migration, the Subandean belt witnessed rapid similar to 80 km eastward propagation of thrusting to the La Vertiente structure at 6 Ma. As there is no evidence for this event of thrust front migration being linked to an increase in shortening rate, the enhanced frontal accretion suggests a shift to supercritical wedge taper conditions. We propose that the supercritical state was due to a drop in basal strength, caused by sediment loading and pore fluid overpressure. This scenario implies that climate-controlled variation in erosional efficiency was the driver of late Miocene mass redistribution, which induced flexural rebound of the Subandean thrust belt, spreading of a large clastic wedge across the basin, and subsequent thrust-front propagation.
Lake sediments in 10 Ethiopian, Kenyan, and Tanzanian rift basins suggest that there were three humid periods at 2.7 to 2.5 million years ago (Ma), 1.9 to 1.7 Ma, and 1.1 to 0.9 Ma, superimposed on the longer-term aridification of East Africa. These humid periods correlate with increased aridity in northwest and northeast Africa and with substantial global climate transitions. These episodes could have had important impacts on the speciation and dispersal of mammals and hominins, because a number of key events, such as the origin of the genus Homo and the evolution of the species Homo erectus, took place in this region during that time.
The late Cenozoic climate of East Africa is punctuated by episodes of short, alternating periods of extreme wetness and aridity, superimposed on a regime of subdued moisture availability exhibiting a long-term drying trend. These periods of extreme climate variability appear to correlate with maxima in the 400-thousand-year (kyr) component of the Earth's eccentricity cycle. Prior to 2.7 Ma the wet phases appear every 400 kyrs, whereas after 2.7 Ma, the wet phases appear every 800 kyrs, with periods of precessional-forced extreme climate variability at 2.7-2.5 Ma, 1.9-1.7 Ma, and 1.1-0.9 Ma before present. The last three major lake phases occur at the times of major global climatic transitions, such as the onset of Northern Hemisphere Glaciation (2.7-2.5 Ma), intensification of the Walker Circulation (1.9-1.7 Ma), and the Mid-Pleistocene Revolution (1.0-0.7 Ma). High-latitude forcing is required to compress the Intertropical Convergence Zone so that East Africa becomes locally sensitive to precessional forcing, resulting in rapid shifts from wet to dry conditions. These periods of extreme climate variability may have provided a catalyst for evolutionary change and driven key speciation and dispersal events amongst mammals and hominins in East Africa. (C) 2007 Elsevier Ltd. All rights reserved.
The development of rise Cenozoic East African Rift System (EARS) profoundly re-shaped the landscape and significantly increased the amplitude of short-term environmental response to climate variation. In particular, the development of amplifier lakes in rift basins after three million years ago significantly contributed to this exceptional sensitivity of East Africa to climate change compared to elsewhere on the African continent. Amplifier lakes are characterized by tectonically-formed graben morphologies in combination with an extreme contrast between high precipitation in the elevated parts of the catchment and high evaporation in the lake area. Such amplifier lakes respond rapidly to moderate, precessional-forced climate shifts, and as they do so apply dramatic environmental pressure to the biosphere. Rift basins, when either extremely dry or lake-filled, form important barriers for migration, mixing and competition of different populations of animals and hominins. Amplifier lakes link long-term, high-amplitude tectonic processes and short-term environmental fluctuations. East Africa may have become the place where early humans evolved as a consequence of this strong link between different time scales. (C) 2010 Elsevier Ltd. All rights reserved.
Variations in the temporal and spatial distribution of solar radiation caused by orbital changes provide a partial explanation for the observed long-term fluctuations in African lake levels. The understanding of such relationships is essential for designing climate-prediction models for the tropics. Our assessment of the nature and timing of East African climate change is based on lake-level fluctuations of Lake Naivasha in the Central Kenya Rift (0°55'S 36°20'E), inferred from sediment characteristics, diatoms, authigenic mineral assemblages and 17 single-crystal 40Ar/39Ar age determinations. Assuming that these fluctuations reflect climate changes, the Lake Naivasha record demonstrates that periods of increased humidity in East Africa mainly followed maximum equatorial solar radiation in March or September. Interestingly, the most dramatic change in the Naivasha Basin occurred as early as 146 kyr BP and the highest lake level was recorded at about 139 to 133 kyr BP. This is consistent with other well-dated low-latitude climate records, but does not correspond to peaks in Northern Hemisphere summer insolation as the trigger for the ice- age cycles. The Naivasha record therefore provides evidence for low-latitude forcing of the ice-age climate cycles.
Present erosion and sediment flux in the semi-arid intramontane Santa Maria Basin, NW Argentina are compared with conditions during a period of wetter and more variable climate at about 30,000 14C yrs ago. The results suggest that the influence of climate change on the overall erosional sediment budget is significant, mainly because of a change in the erosion regime coupled with an increase in mass movements. The most effective mechanism to increase landslide activity in this environment is a highly variable climate on inter-annual timescales. In contrast, Quaternary changes in erosional budgets due to variations in moisture regimes is small in the Santa Maria Basin. Since the magnitude of a potential increase in background erosion as well as enhanced landsliding is smaller than typical levels of uncertainty of erosional budgets for such large basins, it is not likely that climate-driven erosional unloading can influence tectonic style and rates in this semi-arid environment on time scales of several 103 to 104 years.
Multiple landslide clusters record quaternary climate changes in the northwestern Argentine andes
(2003)
The chronology of multiple landslide deposits and related lake sediments in the semi-arid eastern Argentine Cordillera suggests that major mass movements cluster in two time periods during the Quaternary, i.e. between 40 and 25 and after 5 14C kyr BP. These clusters may correspond to the Minchin (maximum at around 28-27 14C kyr BP) and Titicaca wet periods (after 3.9 14C kyr BP). The more humid conditions apparently caused enhanced landsliding in this environment. In contrast, no landslide-related damming and associated lake sediments occurred during the Coipasa (11.5- 10 14C yr BP) and Tauca wet periods (14.5-11 14C yr BP). The two clusters at 40-25 and after 5 14C kyr BP may correspond to periods where the El Niño-Southern Oscillation (ENSO) and Tropical Atlantic Sea Surface Temperature Variability (TAV) were active. This, however, was not the case during the Coipasa and Tauca wet periods. Lake-balance modelling of a landslide-dammed lake suggests a 10-15% increase in precipitation and a 3-4 ° C decrease in temperature at ~30 14C kyr BP as compared to the present. In addition, time-series analysis reveals a strong ENSO and TAV during that time. The landslide clusters in northwestern Argentina are therefore best explained by periods of more humid and more variable climates.
Fluvial fill terraces in intermontane basins are valuable geomorphic archives that can record tectonically and/or climatically driven changes of the Earth-surface process system. However, often the preservation of fill terrace sequences is incomplete and/or they may form far away from their source areas, complicating the identification of causal links between forcing mechanisms and landscape response, especially over multi-millennial timescales. The intermontane Toro Basin in the southern Central Andes exhibits at least five generations of fluvial terraces that have been sculpted into several-hundred-meter-thick Quaternary valley-fill conglomerates. New surface-exposure dating using nine cosmogenic Be-10 depth profiles reveals the successive abandonment of these terraces with a 100 kyr cyclicity between 75 +/- 7 and 487 +/- 34 ka. Depositional ages of the conglomerates, determined by four Al-26/Be-10 burial samples and U-Pb zircon ages of three intercalated volcanic ash beds, range from 18 +/- 141 to 936 +/- 170 ka, indicating that there were multiple cut-and-fill episodes. Although the initial onset of aggradation at similar to 1 Ma and the overall net incision since ca. 500 ka can be linked to tectonic processes at the narrow basin outlet, the superimposed 100 kyr cycles of aggradation and incision are best explained by eccentricity-driven climate change. Within these cycles, the onset of river incision can be correlated with global cold periods and enhanced humid phases recorded in paleoclimate archives on the adjacent Bolivian Altiplano, whereas deposition occurred mainly during more arid phases on the Altiplano and global interglacial periods. We suggest that enhanced runoff during global cold phases - due to increased regional precipitation rates, reduced evapotranspiration, or both - resulted in an increased sediment-transport capacity in the Toro Basin, which outweighed any possible increases in upstream sediment supply and thus triggered incision. Compared with two nearby basins that record precessional (21-kyr) and long-eccentricity (400-kyr) forcing within sedimentary and geomorphic archives, the recorded cyclicity scales with the square of the drainage basin length. (C) 2017 Elsevier B.V. All rights reserved.
Terrestrial cosmogenic nuclide (TCN) concentrations in fluvial sediment, from which denudation rates are commonly inferred, can be affected by hillslope processes. TCN concentrations in gravel and sand may differ if localized, deep-excavation processes (e.g. landslides, debris flows) affect the contributing catchment, whereas the TCN concentrations of sand and gravel tend to be more similar when diffusional processes like soil creep and sheetwash are dominant. To date, however, no study has systematically compared TCN concentrations in different detrital grain-size fractions with a detailed inventory of hillslope processes from the entire catchment. Here we compare concentrations of the TCN Be-10 in 20 detrital sand samples from the Quebrada del Toro (southern Central Andes, Argentina) to a hillslope-process inventory from each contributing catchment. Our comparison reveals a shift from low-slope gullying and scree production in slowly denuding, low-slope areas to steep-slope gullying and landsliding in fast-denuding, steep areas. To investigate whether the nature of hillslope processes (locally excavating or more uniformly denuding) may be reflected in a comparison of the Be-10 concentrations of sand and gravel, we define the normalized sand-gravel index (NSGI) as the Be-10-concentration difference between sand and gravel divided by their summed concentrations. We find a positive, linear relationship between the NSGI and median slope, such that our NSGI values broadly reflect the shift in hillslope processes from low-slope gullying and scree production to steep-slope gullying and landsliding. Higher NSGI values characterize regions affected by steep-slope gullying or landsliding. We relate the large scatter in the relationship, which is exhibited particularly in low-slope areas, to reduced hillslope-channel connectivity and associated transient sediment storage within those catchments. While high NSGI values in well-connected catchments are a reliable signal of deep-excavation processes, hillslope excavation processes may not be reliably recorded by NSGI values where sediment experiences transient storage. (c) 2018 John Wiley & Sons, Ltd.
The northward motion of the Pamir indenter with respect to Eurasia has resulted in coeval thrusting, strike-slip faulting, and normal faulting. The eastern Pamir is currently deformed by east-west oriented extension, accompanied by uplift and exhumation of the Kongur Shan (7719m) and Muztagh Ata (7546m) gneiss domes. Both domes are an integral part of the footwall of the Kongur Shan extensional fault system (KES), a 250 km long, north-south oriented graben. Why active normal faulting within the Pamir is primarily localized along the KES and not distributed more widely throughout the orogen has remained unclear. In addition, relatively little is known about how deformation has evolved throughout the Cenozoic, despite refined estimates on present-day crustal deformation rates and microseismicity, which indicate where crustal deformation is presently being accommodated. To better constrain the spatiotemporal evolution of faulting along the KES, we present 39 new apatite fission track, zircon U-Th-Sm/He, and Ar-40/Ar-39 cooling ages from a series of footwall transects along the KES graben shoulder. Combining these data with present-day topographic relief, 1-D thermokinematic and exhumational modeling documents successive stages, rather than synchronous deformation and gneiss dome exhumation. While the exhumation of the Kongur Shan commenced during the late Miocene, extensional processes in the Muztagh Ata massif began earlier and have slowed down since the late Miocene. We present a new model of synorogenic extension suggesting that thermal and density effects associated with a lithospheric tear fault along the eastern margin of the subducting Alai slab localize extensional upper plate deformation along the KES and decouple crustal motion between the central/western Pamir and eastern Pamir/Tarim basin.
Erosional exhumation and topography in mountain belts are temporally and spatially variable over million year timescales because of changes in both the location of deformation and climate. We investigate spatiotemporal variations in exhumation across a 150 x 250 km compartment of the NW Himalaya, India. Twenty-four new and 241 previously published apatite and zircon fission track and white mica Ar-40/Ar-39 ages are integrated with a 1-D numerical model to quantify rates and timing of exhumation alongstrike of several major structures in the Lesser, High, and Tethyan Himalaya. Analysis of thermochronometer data suggests major temporal variations in exhumation occurred in the early middle Miocene and at the Plio-Pleistocene transition. (1) Most notably, exhumation rates for the northern High Himalayan compartments were high (2-3 mm a(-1)) between similar to 23-19 and similar to 3-0 Ma and low (0.5-0.7 mm a(-1)) in between similar to 19-3 Ma. (2) Along the southern High Himalayan slopes, however, high exhumation rates of 1-2 mm a(-1) existed since 11 Ma. (3) Our thermochronology data sets are poorly correlated with present-day rainfall, local relief, and specific stream power which may likely result from (1) a lack of sensitivity of changes in crustal cooling to spatial variations in erosion at high exhumation rates (>similar to 1 mm a(-1)), (2) spatiotemporal variation in erosion not mimicking the present-day topographic or climatic conditions, or (3) the thermochronometer samples in this region having cooled under topography that only weakly resembled the modern-day topography.
Along the Southern Himalayan Front (SHF), areas with concentrated precipitation coincide with rapid exhumation, as indicated by young mineral cooling ages. Twenty new, young ( < 1-5 Ma) apatite fission track (AFT) ages have been obtained from the Himalayan Crystalline Core along the Sutlej Valley, NW India. The AFT ages correlate with elevation, but show no spatial relationship to tectonic structures, such as the Main Central Thrust or the Southern Tibetan Fault System. Monsoonal precipitation in this region exerts a strong influence on erosional surface processes. Fluvial erosional unloading along the SHF is focused on high mountainous areas, where the orographic barrier forces out > 80% of the annual precipitation. AFT cooling ages reveal a coincidence between rapid erosion and exhumation that is focused in a similar to 50-70-km-wide sector of the Himalaya, rather than encompassing the entire orogen. Assuming simplified constant exhumation rates, the rocks of two age vs. elevation transects were exhumed at similar to 1.4 +/- 0.2 and similar to 1.1 +/- 0.4 mm/a with an average cooling rate of similar to 40-50degreesC/Ma during Pliocene-Quarternary time. Following other recently published hypotheses regarding the relation between tectonics and climate in the Himalaya, we suggest that this concentrated loss of material was accommodated by motion along a back-stepping thrust to the south and a normal fault zone to the north as part of an extruding wedge. Climatically controlled erosional processes focus on this wedge and suggest that climatically controlled surface processes determine tectonic deformation in the internal part of the Himalaya. (C) 2004 Elsevier B.V. All rights reserved
Metamorphic dome complexes occur within the internal structures of the northern Himalaya and southern Tibet. Their origin, deformation, and fault displacement patterns are poorly constrained. We report new field mapping, structural data, and cooling ages from the western flank of the Leo Pargil dome in the northwestern Himalaya in an attempt to characterize its post-middle Miocene structural development. The western flank of the dome is characterized by shallow, west-dipping pervasive foliation and WNW-ESE mineral lineation. Shear-sense indicators demonstrate that it is affected by east-west normal faulting that facilitated exhumation of high-grade metamorphic rocks in a contractional setting. Sustained top-to-northwest normal faulting during exhumation is observed in a progressive transition from ductile to brittle deformation. Garnet and kyanite indicate that the Leo Pargil dome was exhumed from the mid-crust. Ar- 40/Ar-39 mica and apatite fission track (AFT) ages constrain cooling and exhumation pathways front 350 to 60 degrees C and suggest that the dome cooled in three stages since the middle Miocene. Ar-40/Ar-39 white mica ages of 16-14 Ma suggest a first phase of rapid cooling and provide minimum estimates for the onset of dome exhumation. AFT ages between 10 and 8 Ma suggest that ductile fault displacement had ceased by then, and AFT track-length data from high-elevation samples indicate that the rate of cooling had decreased significantly. We interpret this to indicate decreased fault displacement along the Leo Pargil shear zone and possibly a transition to the Kaurik-Chango normal fault system between 10 and 6 Ma. AFT ages from lower elevations indicate accelerated cooling since the Pliocene that cannot be related to pure fault displacement, and therefore may reflect more pronounced regionally distributed and erosion-driven exhumation
Whether variations in the spatial distribution of erosion influence the location, style, and magnitude of deformation within the Himalayan orogen is a matter of debate. We report new Ar-40/Ar-39 white mica and apatite fission- track (AFT) ages that measure the vertical component of exhumation rates along an similar to 120-km-wide NE-SW transect spanning the greater Sutlej region of northwest India. The Ar-40/Ar-39 data indicate that first the High Himalayan Crystalline units cooled below their closing temperature during the early to middle Miocene. Subsequently, Lesser Himalayan Crystalline nappes cooled rapidly, indicating southward propagation of the orogen during late Miocene to Pliocene time. The AFT data, in contrast, imply synchronous exhumation of a NE-SW-oriented similar to 80 x 40 km region spanning both crystalline nappes during the Pliocene-Quaternary. The locus of pronounced exhumation defined by the AFT data correlates with a region of high precipitation, discharge, and sediment flux rates during the Holocene. This correlation suggests that although tectonic processes exerted the dominant control on the denudation pattern before and until the middle Miocene; erosion may have been the most important factor since the Pliocene
We investigated deep-seated gravitational slope deformation (DSGSD) and slow mass movements in the southern Tien Shan Mountains front using synthetic aperture radar (SAR) time-series data obtained by the ALOS/PALSAR satellite. DSGSD evolves with a variety of geomorphological changes (e.g. valley erosion, incision of slope drainage networks) over time that affect earth surfaces and, therefore, often remain unexplored. We analysed 118 interferograms generated from 20 SAR images that covered about 900 km(2). To understand the spatial pattern of the slope movements and to identify triggering parameters, we correlated surface dynamics with the tectono-geomorphic processes and lithologic conditions of the active front of the Alai Range. We observed spatially continuous, constant hillslope movements with a downslope speed of approximately 71 mm year(-1) velocity. Our findings suggest that the lithological and structural framework defined by protracted deformation was the main controlling factor for sustained relief and, consequently, downslope mass movements. The analysed structures revealed integration of a geological/structural setting with the superposition of Cretaceous-Paleogene alternating carbonatic and clastic sedimentary structures as the substratum for younger, less consolidated sediments. This type of structural setting causes the development of large-scale, gravity-driven DSGSD and slow mass movement. Surface deformations with clear scarps and multiple crest lines triggered planes for large-scale deep mass creeps, and these were related directly to active faults and folds in the geologic structures. Our study offers a new combination of InSAR techniques and structural field observations, along with morphometric and seismologic correlations, to identify and quantify slope instability phenomena along a tectonically active mountain front. These results contribute to an improved natural risk assessment in these structures.
Orogenic plateaus are extensive, high-elevation areas with low internal relief that have been attributed to deep-seated and/or climate-driven surface processes. In the latter case, models predict that lateral plateau growth results from increasing aridity along the margins as range uplift shields the orogen interior from precipitation. We analyze the spatiotemporal progression of basin isolation and filling at the eastern margin of the Puna Plateau of the Argentine Andes to determine if the topography predicted by such models is observed. We find that the timing of basin filling and reexcavation is variable, suggesting nonsystematic plateau growth. Instead, the Airy isostatically compensated component of topography constitutes the majority of the mean elevation gain between the foreland and the plateau. This indicates that deep-seated phenomena, such as changes in crustal thickness and/or lateral density, are required to produce high plateau elevations. In contrast, the frequency of the uncompensated topography within the plateau and in the adjacent foreland that is interrupted by ranges appears similar, although the amplitude of this topographic component increases east of the plateau. Combined with sedimentologic observations, we infer that the low internal relief of the plateau likely results from increased aridity and sediment storage within the plateau and along its eastern margin.
The timing of the late Cenozoic collision between the Pamir salient and the Tien Shan as well as changes in the relative motion between the Pamir and Tarim are poorly constrained. The northern margin of the Pamir salient indented northward by similar to 300 km during the late Cenozoic, accommodated by south-dipping intracontinental subduction along the Main Pamir Thrust (MPT) coupled to strike-slip faults on the eastern flank of the orogen and both strike-slip and thrust faults on the western margin. The Kashgar-Yecheng transfer system (KYTS) is the main dextral slip shear zone separating Tarim from the Eastern Pamir, with an estimated cumulative offset of similar to 280 km at an average late Cenozoic dextral slip rate of 11-15 mm/a (Cowgill, 2010). In order to better constrain the slip history of the KYTS, we collected thermochronologic samples along the eastward-flowing, deeply incised, antecedent Tashkorgan-Yarkand River, which crosses the fault system on the eastern flank of the orogen. We present 29 new biotite (40)Ar/(39)Ar ages, apatite and zircon (U-Th-Sm)/He ages, and apatite fission track (AFT) analysis, combined with published muscovite and biotite (40)Ar/(39)Ar and AFT data, to create a unique thermochronologic dataset in this poorly studied and remote region. We constrain the timing of four major N-trending faults: the latter three are strands of the KYTS. The westernmost, the Kuke fault, experienced significant dip-slip, west-side-up displacement between > 12 and 6 Ma. To the east, within the KYTS, our new thermochronologic data and geomorphic observations suggest that the Kumtag and Kusilaf dextral slip faults have been inactive since at least 3-5 Ma. Long-term incision rates across the Aertashi dextral slip fault, the easternmost strand of the KYTS, are compatible with slow horizontal slip rates of 1.7-5.3 mm/a over the past 3 to 5 Ma. In summary, these data show that the slip rate of the KYTS decreased substantially during the late Miocene or Pliocene. Furthermore, Miocene-present regional kinematic reconstructions suggest that this deceleration reflects the substantial increase of northward motion of Tarim rather than a significant decrease of the northward velocity of the Pamir. (C) 2011 Elsevier B.V. All rights reserved.
The northern part of the Pamir orogen is the preeminent example of an active intracontinental subduction zone in the early stages of continent-continent collision. Such zones are the least understood type of plate boundaries because modern examples are few and of limited access, and ancient analogs have been extensively overprinted by subsequent tectonic and erosion processes. In the Pamir, it has been assumed that most of the plate convergence was accommodated by overthrusting along the plate-bounding Main Pamir Thrust (MPT), which forms the principal northern mountain and deformation front of the Pamir. However, the synopsis of our new and previously published thermochronologic data from this region shows that the hanging wall of the MPT experienced relatively minor amounts of late Cenozoic exhumation. The Pamir orogen as a whole is an integral part of the overriding plate in a subduction system, while the remnant basin to the north constitutes the downgoing plate, with the bulk of the convergence accommodated by underthrusting. Herein, we demonstrate that the observed deformation of the upper and lower plates within the Pamir-Alai convergence zone resembles highly arcuate oceanic subduction systems characterized by slab rollback, subduction erosion, subduction accretion, and marginal slab-tear faults. We suggest that the curvature of the North Pamir is genetically linked to the short width and rollback of the south-dipping Alai slab; northward motion (indentation) of the Pamir is accommodated by crustal processes related to this rollback. The onset of south-dipping subduction is tentatively linked to intense Pamir contraction following break-off of the north-dipping Indian slab beneath the Karakoram.
The Kenya rift revisited
(2017)
We present three-dimensional (3-D) models that describe the present-day thermal and rheological state of the lithosphere of the greater Kenya rift region aiming at a better understanding of the rift evolution, with a particular focus on plume-lithosphere interactions. The key methodology applied is the 3-D integration of diverse geological and geophysical observations using gravity modelling. Accordingly, the resulting lithospheric-scale 3-D density model is consistent with (i) reviewed descriptions of lithological variations in the sedimentary and volcanic cover, (ii) known trends in crust and mantle seismic velocities as revealed by seismic and seismological data and (iii) the observed gravity field. This data-based model is the first to image a 3-D density configuration of the crystalline crust for the entire region of Kenya and northern Tanzania. An upper and a basal crustal layer are differentiated, each composed of several domains of different average densities. We interpret these domains to trace back to the Precambrian terrane amalgamation associated with the East African Orogeny and to magmatic processes during Mesozoic and Cenozoic rifting phases. In combination with seismic velocities, the densities of these crustal domains indicate compositional differences. The derived lithological trends have been used to parameterise steady-state thermal and rheological models. These models indicate that crustal and mantle temperatures decrease from the Kenya rift in the west to eastern Kenya, while the integrated strength of the lithosphere increases. Thereby, the detailed strength configuration appears strongly controlled by the complex inherited crustal structure, which may have been decisive for the onset, localisation and propagation of rifting.
Though orogen-parallel shortening and vertical extension have dominated the tectonic evolution of the central Andes, a significant kinematic shift from horizontal contraction to extension appears to have occurred within the high Puna-Altiplano Plateau, with the establishment of extension oblique to the orogen since late Miocene time. We present data from the southern margin of the Puna Plateau, NW Argentina, where new normal faults have been documented in the Fiambala, Punta Negra, and La Quebrada areas. The unifying characteristics of these areas are that young normal faults reactivate or crosscut older thrust and reverse faults. The relationship between the faults and the late Miocene- Pliocene Punaschotter conglomerate suggests that the extensional faulting must be younger than 3.5 to 7 Ma. Existing data are incomplete but indicate that similar horizontal extension has occurred in many regions throughout the Puna- Altiplano Plateau, while shortening continues along the plateau margins. Given the spatial and temporal distribution of this late Miocene to Pliocene kinematic shift, both lithospheric loss in the Puna Plateau and plateau-wide gravitational extensional spreading enhanced by slowing of plate convergence rate could be responsible. The young, disorganized, horizontal extension in the Andes today may be the precursor to more pronounced extension such as observed on the Tibetan Plateau since mid-Miocene time.
Modification of the landscape by glacial erosion reflects the dynamic interplay of climate through temperature, precipitation, and prevailing wind direction, and tectonics through rock uplift and exhumation rate, lithology, and range and fault geometry. We investigate these relationships in the northeast Pamir Mountains using mapping and dating of moraines and terraces to determine the glacial history. We analyze modem glacial morphology to determine glacier area, spacing, headwall relief, debris cover, and equilibrium line altitude (ELA) using the area x altitude balance ratio (AABR), toe-to-headwall altitude ratio (THAR) and toe-to-summit altitude method (TSAM) for 156 glaciers and compare this to lithologic, tectonic, and climatic data We observe a pronounced asymmetry in glacial ELA, area, debris cover, and headwall relief that we interpret to reflect both structural and climatic control: glaciers on the downwind (eastern) side of the range are larger, more debris covered, have steeper headwalls, and tend to erode headward, truncating the smaller glaciers of the upwind, fault-controlled side of the range. We explain this by the transfer of moisture deep into the range as wind-blown or avalanched snow and by limitations imposed on glacial area on the upwind side of the range by the geometry of the Kongur extensional system (KES). The correspondence between rapid exhumation along the KES and maxima in glacier debris cover and headwall relief and minimums in all measures of ELA suggest that taller glacier headwalls develop in a response to more rapid exhumation rates. However, we find that glaciers in the Muji valley did not extend beyond the range front until at least 43 ka, in contrast to extensive glaciation since 300 ka in the south around the high peaks, a pattern which does not clearly reflect uplift rate. Instead, the difference in glacial history and the presence of large peaks (Muztagh Ata and Kongur Shan) with flanking glaciers likely reflects lithologic control (i.e., the location of crustal gneiss domes) and the formation of peaks that rise above the ELA and escape the glacial buzzsaw. (C) 2014 Elsevier B.V. All rights reserved.
Fluvial fill terraces preserve sedimentary archives of landscape responses to climate change, typically over millennial timescales. In the Humahuaca Basin of NW Argentina (Eastern Cordillera, southern Central Andes), our 29 new optically stimulated luminescence ages of late Pleistocene fill terrace sediments demonstrate that the timing of past river aggradation occurred over different intervals on the western and eastern sides of the valley, despite their similar bedrock lithology, mean slopes, and precipitation. In the west, aggradation coincided with periods of increasing precipitation, while in the east, aggradation coincided with decreasing precipitation or more variable conditions. Erosion rates and grain size dependencies in our cosmogenic Be-10 analyses of modern and fill terrace sediments reveal an increased importance of landsliding compared to today on the west side during aggradation, but of similar importance during aggradation on the east side. Differences in the timing of aggradation and the Be-10 data likely result from differences in valley geometry, which causes sediment to be temporarily stored in perched basins on the east side. It appears as if periods of increasing precipitation triggered landslides throughout the region, which induced aggradation in the west, but blockage of the narrow bedrock gorges downstream from the perched basins in the east. As such, basin geometry and fluvial connectivity appear to strongly influence the timing of sediment movement through the system. For larger basins that integrate subbasins with differing geometries or degrees of connectivity (like Humahuaca), sedimentary responses to climate forcing are likely attenuated.
Sr isotope records from marginal marine basins track the mixing between seawater and local continental runoff, potentially recording the effects of sea level, tectonic, and climate forcing in marine fossils and sediments. Our 110 new Sr-87/Sr-86 analyses on oyster and foraminifera samples from six late Miocene stratigraphic sections in southern Turkey, Crete, and Sicily show that Sr-87/Sr-86 fell below global seawater values in the basins several million years before the Messinian Salinity Crisis, coinciding with tectonic uplift and basin shallowing. 87Sr/86Sr from more centrally located basins (away from the Mediterranean coast) drop below global seawater values only during the Messinian Salinity Crisis. In addition to this general trend, 55 new Sr-87/Sr-86 analyses from the astronomically tuned Lower Evaporites in the central Apennines (Italy) allow us to explore the effect of glacio-eustatic sea level and precipitation changes on Sr-87/Sr-86. Most variation in our data can be explained by changes in sea level, with greatest negative excursions from global seawater values occurring during relative sea level lowstands, which generally coincided with arid conditions in the Mediterranean realm. We suggest that this greater sensitivity to lowered sea level compared with higher runoff could relate to the inverse relationship between Sr concentration and river discharge. Variations in the residence time of groundwater within the karst terrain of the circum-Mediterranean region during arid and wet phases may help to explain the single (robust) occurrence of a negative excursion during a sea level highstand, but this explanation remains speculative without more detailed paleoclimatic data for the region.
The southwest margin of the Central Anatolian Plateau has experienced multiple phases of topographic growth, including the formation of localized highs prior to the Late Miocene that were later affected by wholesale uplift of the plateau margin. Our new biostratigraphic data limit the age of uplifted marine sediments at the southwest plateau margin at 1.5 km elevation to <7.17 Ma, and regional lithostratigraphic correlations imply that the age is <6.7 Ma. Single-grain CA-TIMS U-Pb zircon analyses from a reworked ash within the marine sediments yield dates as young as 10.6 Ma, indicating a maximum age that is consistent with the biostratigraphy. Our structural measurements within the uplifted region and fault inversion modeling agree with previous findings in surrounding regions, with early contraction followed by strike-slip and extensional deformation during uplift. Focal mechanisms from shallow earthquakes show that the extensional phase has continued to the present. Broad similarities in the change in the tectonic stress regime (after 8 Ma) and the onset of surface uplift (after 7 Ma) imply that deep-seated process(es) caused post-7 Ma uplift. The geometry of lithospheric slabs beneath the plateau margin, Pliocene to recent alkaline volcanism, and the uplift pattern with accompanying normal faulting point toward slab tearing and localized heating at the base of the lithosphere as a probable mechanism for post-7 Ma uplift of the southwest margin. Considering previous work in the region, there appears to be an important link between slab dynamics and surface uplift throughout the Anatolian Plateau’s southern margin.
Uplifted Neogene marine sediments and Quaternary fluvial terraces in the Mut Basin, southern Turkey, reveal a detailed history of surface uplift along the southern margin of the Central Anatolian plateau from the Late Miocene to the present. New surface exposure ages (Be-10, Al-26, and Ne-21) of gravels capping fluvial strath terraces located between 28 and 135 m above the Goksu River in the Mut Basin yield ages ranging from ca. 25 to 130 ka, corresponding to an average incision rate of 0.52 to 0.67 mm/yr. Published biostratigraphic data combined with new interpretations of the fossil assemblages from uplifted marine sediments reveal average uplift rates of 0.25 to 0.37 mm/yr since Late Miocene time (starting between 8 and 5.45 Ma), and 0.72 to 0.74 mm/yr after 1.66 to 1.62 Ma. Together with the terrace abandonment ages, the data imply 0.6 to 0.7 mm/yr uplift rates from 1.6 Ma to the present. The different post-Late Miocene and post-1.6 Ma uplift rates can imply increasing uplift rates through time, or multi-phased uplift with slow uplift or subsidence in between. Longitudinal profiles of rivers in the upper catchment of the Mut and Ermenek basins show no apparent lithologic or fault control on some knickpoints that occur at 1.2 to 1.5 km elevation, implying a transient response to a change in uplift rates. Projections of graded upper relict channel segments to the modern outlet, together with constraints from uplifted marine sediments, show that a slower incision/uplift rate of 0.1 to 0.2 mm/yr preceded the 0.7 mm/yr uplift rate. The river morphology and profile projections therefore reflect multi-phased uplift of the plateau margin, rather than steadily increasing uplift rates. Multi-phased uplift can be explained by lithospheric slab break-off and possibly also the arrival of the Eratosthenes Seamount at the collision zone south of Cyprus.
Despite global warming and unlike their Himalayan neighbours, glaciers in the Karakoram mountains do not show signs of significant retreat. Here we report high velocity variations of Biafo Glacier, central Karakoram, which occurred between 2001 and 2009 and which indicate considerable dynamics in its flow behaviour. We have generated a dense time series of glacier surface velocities, based on cross-correlation of optical satellite images, which clearly shows seasonal and interannual velocity variations, reaching 50% in some places. The interannual velocity variations resemble the passing of a broad wave of high velocities, with peak velocities during 2005 and some diffusion down-glacier over a period of at least 4 years. High interannual velocity variations are also observed at other glaciers in the vicinity, suggesting a common cause, although these appear to partly comprise longer acceleration phases. Analysis of weather station data provides some indications of meteorological conditions that could have promoted sustained sliding events during this period, but this does not explain the wave-like nature of the acceleration at Biafo Glacier, and the regular, protracted velocity changes.
Some of the largest and most erosive floods on Earth result from the failure of glacial dams. While potentially cataclysmic ice dams are recognized to have repeatedly formed along ice-sheet margins, much less is known about the frequency and longevity of ice dams caused by mountain glaciers, and their impact on landscape evolution. Here we present field observations and results from cosmogenic nuclide dating that allow reconstructing a > 100-k.y.-long history of glacial damming in the Shyok Valley, eastern Karakoram (South Asia). Our field observations provide evidence that Asia's second-longest glacier, the Siachen, once extended for over 180 km and blocked the Shyok River during the penultimate glacial period, leading to upstream deposition of a more than 400-m-thick fluvio-lacustrine valley fill. Be-10-depth profile modeling indicates that glacial damming ended with the onset of the Eemian interglacial and that the Shyok River subsequently incised the valley fill at an average rate of similar to 4-7 m k.y.(-1). Comparison with contemporary ice-dammed lakes in the Karakoram and elsewhere suggests recurring outburst floods during the aggradation period, while over 25 cycles of fining-upward lake deposits within the valley fill indicate impounding of floods from farther upstream. Despite prolonged damming, the net effect of this and probably earlier damming episodes by the Siachen Glacier is dominated by glacial erosion in excess of fluvial incision, as evidenced by a pronounced overdeepening that follows the glaciated valley reach. Strikingly similar overdeepened valleys at all major confluences of the Shyok and Indus Rivers with Karakoram tributaries indicate that glacial dams and subsequent outburst floods have been widespread and frequent in this region during the Quaternary. Our study suggests that the interaction of Karakoram glaciers with the Shyok and Indus Rivers promoted valley incision and headward erosion into the western margin of the Tibetan Plateau.
The response of surface processes to climatic forcing is fundamental for understanding the impacts of climate change on landscape evolution. In the Himalaya, most large rivers feature prominent fill terraces that record an imbalance between sediment supply and transport capacity, presumably due to past fluctuations in monsoon precipitation and/or effects of glaciation at high elevation. Here, we present volume estimates, chronological constraints, and Be-10-derived paleo-erosion rates from a prominent valley fill in the Yamuna catchment, Garhwal Himalaya, to elucidate the coupled response of rivers and hillslopes to Pleistocene climate change. Although precise age control is complicated due to methodological problems, the new data support formation of the valley fill during the late Pleistocene and its incision during the Holocene. We interpret this timing to indicate that changes in discharge and river-transport capacity were major controls. Compared to the present day, late Pleistocene hillslope erosion rates were higher by a factor of similar to 2-4, but appear to have decreased during valley aggradation. The higher late Pleistocene erosion rates are largely unrelated to glacial erosion and could be explained by enhanced sediment production on steep hillslopes due to increased periglacial activity that declined as temperatures increased. Alternatively, erosion rates that decrease during valley aggradation are also consistent with reduced landsliding from threshold hillslopes as a result of rising base levels. In that case, the similarity of paleo-erosion rates near the end of the aggradation period with modern erosion rates might imply that channels and hillslopes are not yet fully coupled everywhere and that present-day hillslope erosion rates may underrepresent long-term incision rates. (C) 2015 Elsevier B.V. All rights reserved.
High Asian glacial landscapes have large variations in topographic relief and the size and steepness of snow accumulation areas. Associated differences in glacial cover and dynamics allow a first-order determination of the dominant processes shaping these landscapes. Here we provide a regional synthesis of the topography and flow characteristics of 287 glaciers across High Asia using digital elevation analysis and remotely sensed glacier surface velocities. Glaciers situated in low-relief areas on the Tibetan Plateau are mainly nourished by direct snowfall, have little or no debris cover, and have a relatively symmetrical distribution of velocities along their length. In contrast, avalanche-fed glaciers with steep accumulation areas, which occur at the deeply incised edges of the Tibetan Plateau, are heavily covered with supraglacial debris, and flow velocities are highest along short segments near their headwalls but greatly reduced along their debris-mantled lower parts. The downstream distribution of flow velocities suggests that the glacial erosion potential is progressively shifted upstream as accumulation areas get steeper and hillslope debris fluxes increase. Our data suggest that the coupling of hillslopes and glacial dynamics increases with topographic steepness and debris cover. The melt-lowering effect of thick debris cover allows the existence of glaciers even when they are located entirely below the snow line. However, slow velocities limit the erosion potential of such glaciers, and their main landscape-shaping contribution may simply be the evacuation of debris from the base of glacial headwalls, which inhibits the formation of scree slopes and thereby allows ongoing headwall retreat by periglacial hillslope processes. We propose a conceptual model in which glacially influenced plateau margins evolve from low-relief to high-relief landscapes with distinctive contributions of hillslope processes and glaciers to relief production and decay.
Controversy about the current state and future evolution of Himalayan glaciers has been stirred up by erroneous statements in the fourth report by the Intergovernmental Panel on Climate Change(1,2). Variable retreat rates(3-6) and a paucity of glacial mass-balance data(7,8) make it difficult to develop a coherent picture of regional climate-change impacts in the region. Here, we report remotely-sensed frontal changes and surface velocities from glaciers in the greater Himalaya between 2000 and 2008 that provide evidence for strong spatial variations in glacier behaviour which are linked to topography and climate. More than 65% of the monsoon-influenced glaciers that we observed are retreating, but heavily debris-covered glaciers with stagnant low-gradient terminus regions typically have stable fronts. Debris-covered glaciers are common in the rugged central Himalaya, but they are almost absent in subdued landscapes on the Tibetan Plateau, where retreat rates are higher. In contrast, more than 50% of observed glaciers in the westerlies-influenced Karakoram region in the northwestern Himalaya are advancing or stable. Our study shows that there is no uniform response of Himalayan glaciers to climate change and highlights the importance of debris cover for understanding glacier retreat, an effect that has so far been neglected in predictions of future water availability(9,10) or global sea level(11).
Erosion in the Himalaya is responsible for one of the greatest mass redistributions on Earth and has fueled models of feedback loops between climate and tectonics. Although the general trends of erosion across the Himalaya are reasonably well known, the relative importance of factors controlling erosion is less well constrained. Here we present 25 Be-10-derived catchment-averaged erosion rates from the Yamuna catchment in the Garhwal Himalaya, northern India. Tributary erosion rates range between similar to 0.1 and 0.5mmyr(-1) in the Lesser Himalaya and similar to 1 and 2mmyr(-1) in the High Himalaya, despite uniform hillslope angles. The erosion-rate data correlate with catchment-averaged values of 5 km radius relief, channel steepness indices, and specific stream power but to varying degrees of nonlinearity. Similar nonlinear relationships and coefficients of determination suggest that topographic steepness is the major control on the spatial variability of erosion and that twofold to threefold differences in annual runoff are of minor importance in this area. Instead, the spatial distribution of erosion in the study area is consistent with a tectonic model in which the rock uplift pattern is largely controlled by the shortening rate and the geometry of the Main Himalayan Thrust fault (MHT). Our data support a shallow dip of the MHT underneath the Lesser Himalaya, followed by a midcrustal ramp underneath the High Himalaya, as indicated by geophysical data. Finally, analysis of sample results from larger main stem rivers indicates significant variability of Be-10-derived erosion rates, possibly related to nonproportional sediment supply from different tributaries and incomplete mixing in main stem channels.
Climatic controls on debris-flow activity and sediment aggradation: The Del Medio fan, NW Argentina
(2016)
In the Central Andes, several studies on alluvial terraces and valley fills have linked sediment aggradation to periods of enhanced sediment supply. However, debate continues over whether tectonic or climatic factors are most important in triggering the enhanced supply. The Del Medio catchment in the Humahuaca Basin (Eastern Cordillera, NW Argentina) is located within a transition zone between subhumid and arid climates and hosts the only active debris-flow fan within this intermontane valley. By combining Be-10 analyses of boulder and sediment samples within the Del Medio catchment, with regional morphometric measurements of nearby catchments, we identify the surface processes responsible for aggradation in the Del Medio fan and their likely triggers. We find that the fan surface has been shaped by debris flows and channel avulsions during the last 400 years. Among potential tectonic, climatic, and autogenic factors that might influence deposition, our analyses point to a combination of several favorable factors that drive aggradation. These are in particular the impact of occasional abundant rainfall on steep slopes in rock types prone to failure, located in a region characterized by relatively low rainfall amounts and limited transport capacity. These characteristics are primarily associated with the climatic transition zone between the humid foreland and the arid orogen interior, which creates an imbalance between sediment supply and sediment transfer. The conditions and processes that drive aggradation in the Del Medio catchment today may provide a modern analog for the conditions and processes that drove aggradation in other nearby tributaries in the past.
In recent decades, slope instability in high-mountain regions has often been linked to increase in temperature and the associated permafrost degradation and/or the increase in frequency/intensity of rainstorm events. In this context we analyzed the spatiotemporal evolution and potential controlling mechanisms of small- to medium-sized mass movements in a high-elevation catchment of the Italian Alps (Sulden/Solda basin). We found that slope-failure events (mostly in the form of rockfalls) have increased since the 2000s, whereas the occurrence of debris flows has increased only since 2010. The current climate-warming trend registered in the study area apparently increases the elevation of rockfall-detachment areas by approximately 300 m, mostly controlled by the combined effects of frost-cracking and permafrost thawing. In contrast, the occurrence of debris flows does not exhibit such an altitudinal shift, as it is primarily driven by extreme precipitation events exceeding the 75th percentile of the intensity-duration rainfall distribution. Potential debris-flow events in this environment may additionally be influenced by the accumulation of unconsolidated debris over time, which is then released during extreme rainfall events. Overall, there is evidence that the upper Sulden/Solda basin (above ca. 2500 m above sea level [a.s.l.]), and especially the areas in the proximity of glaciers, have experienced a significant decrease in slope stability since the 2000s, and that an increase in rockfalls and debris flows during spring and summer can be inferred. Our study thus confirms that "forward-looking" hazard mapping should be undertaken in these increasingly frequented, high-elevation areas of the Alps, as environmental change has elevated the overall hazard level in these regions.
Rainfall in the central Andes associated with the South American Monsoon and the South American Low-Level Jet results from orographic effects on atmospheric circulation exerted by the Andean Plateau and the Eastern Cordillera. However, despite its importance for South American climate, no reliable records exist that allow decoding the evolution of thresholds and interactions between Andean topography and atmospheric circulation, especially regarding the onset of humid conditions in the inherently dry southern central Andes. Here, we employ multi-proxy isotope data of lipid biomarkers, pedogenic carbonates and volcanic glass from the Eastern Cordillera of NW Argentina and present the first long-term evapotranspiration record. We find that regional eco-hydrology and vegetation changes are associated with initiation of moisture transport via the South American Low-Level Jet at 7.6 Ma, and subsequent lateral growth of the orogen at 6.5 Ma. Our results highlight that topographically induced changes in atmospheric circulation patterns, not global climate change, were responsible for late Miocene environmental change in this part of the southern hemisphere. This suggests that mountain building over time fundamentally controlled habitat evolution along the central Andes.
We examined the relationship between the mechanical strength of the lithosphere and the distribution of seismicity within the overriding continental plate of the southern Central Andes (SCA, 29 degrees-39 degrees S), where the oceanic Nazca Plate changes its subduction angle between 33 degrees S and 35 degrees S, from subhorizontal in the north (<5 degrees) to steep in the south (similar to 30 degrees). We computed the long-term lithospheric strength based on an existing 3D model describing variations in thickness, density, and temperature of the main geological units forming the lithosphere of the SCA and adjacent forearc and foreland regions. The comparison between our results and seismicity within the overriding plate (upper-plate seismicity) shows that most of the events occur within the modeled brittle domain of the lithosphere. The depth where the deformation mode switches from brittle frictional to thermally activated ductile creep provides a conservative lower bound to the seismogenic zone in the overriding plate of the study area. We also found that the majority of upper-plate earthquakes occurs within the realm of first-order contrasts in integrated strength (12.7-13.3 log Pam in the Andean orogen vs. 13.5-13.9 log Pam in the forearc and the foreland). Specific conditions characterize the mechanically strong northern foreland of the Andes, where seismicity is likely explained by the effects of slab steepening.
In an ocean-continent subduction zone, the assessment of the lithospheric thermal state is essential to determine the controls of the deformation within the upper plate and the dip angle of the subducting lithosphere. In this study, we evaluate the degree of influence of both the configuration of the upper plate (i.e., thickness and composition of the rock units) and variations of the subduction angle on the lithospheric thermal field of the southern Central Andes (29 degrees-39 degrees S). Here, the subduction angle increases from subhorizontal (5 degrees) north of 33 degrees S to steep (similar to 30 degrees) in the south. We derived the 3D temperature and heat flow distribution of the lithosphere in the southern Central Andes considering conversion of S wave tomography to temperatures together with steady-state conductive thermal modeling. We found that the orogen is overall warmer than the forearc and the foreland and that the lithosphere of the northern part of the foreland appears colder than its southern counterpart. Sedimentary blanketing and the thickness of the radiogenic crust exert the main control on the shallow thermal field (<50km depth). Specific conditions are present where the oceanic slab is relatively shallow (<85 km depth) and the radiogenic crust is thin. This configuration results in relatively colder temperatures compared to regions where the radiogenic crust is thick and the slab is steep. At depths >50km, the temperatures of the overriding plate are mainly controlled by the mantle heat input and the subduction angle. The thermal field of the upper plate likely preserves the flat subduction angle and influences the spatial distribution of shortening.
The southern Central Andes (SCA) (between 27 degrees S and 40 degrees S) is bordered to the west by the convergent margin between the continental South American Plate and the oceanic Nazca Plate. The subduction angle along this margin is variable, as is the deformation of the upper plate. Between 33 degrees S and 35 degrees S, the subduction angle of the Nazca plate increases from sub-horizontal (< 5 degrees) in the north to relatively steep (similar to 30 degrees) in the south. The SCA contain inherited lithological and structural heterogeneities within the crust that have been reactivated and overprinted since the onset of subduction and associated Cenozoic deformation within the Andean orogen. The distribution of the deformation within the SCA has often been attributed to the variations in the subduction angle and the reactivation of these inherited heterogeneities. However, the possible influence that the thickness and composition of the continental crust have had on both short-term and long-term deformation of the SCA is yet to be thoroughly investigated. For our investigations, we have derived density distributions and thicknesses for various layers that make up the lithosphere and evaluated their relationships with tectonic events that occurred over the history of the Andean orogeny and, in particular, investigated the short- and long-term nature of the present-day deformation processes. We established a 3D model of lithosphere beneath the orogen and its foreland (29 degrees S-39 degrees S) that is consistent with currently available geological and geophysical data, including the gravity data. The modelled crustal configuration and density distribution reveal spatial relationships with different tectonic domains: the crystalline crust in the orogen (the magmatic arc and the main orogenic wedge) is thicker (similar to 55 km) and less dense (similar to 2900 kg/m(3)) than in the forearc (similar to 35 km, similar to 2975 kg/m(3)) and foreland (similar to 30 km, similar to 3000 kg/m(3)). Crustal thickening in the orogen probably occurred as a result of stacking of low-density domains, while density and thickness variations beneath the forearc and foreland most likely reflect differences in the tectonic evolution of each area following crustal accretion. No clear spatial relationship exists between the density distribution within the lithosphere and previously proposed boundaries of crustal terranes accreted during the early Paleozoic. Areas with ongoing deformation show a spatial correlation with those areas that have the highest topographic gradients and where there are abrupt changes in the average crustal-density contrast. This suggests that the short-term deformation within the interior of the Andean orogen and its foreland is fundamentally influenced by the crustal composition and the relative thickness of different crustal layers. A thicker, denser, and potentially stronger lithosphere beneath the northern part of the SCA foreland is interpreted to have favoured a strong coupling between the Nazca and South American plates, facilitating the development of a sub-horizontal slab.
Magmatic continental rifts often constitute nascent plate boundaries, yet long-term extension rates and transient rate changes associated with these early stages of continental breakup remain difficult to determine. Here, we derive a time-averaged minimum extension rate for the inner graben of the Northern Kenya Rift (NKR) of the East African Rift System for the last 0.5 m.y. We use the TanDEM-X science digital elevation model to evaluate fault-scarp geometries and determine fault throws across the volcano-tectonic axis of the inner graben of the NKR. Along rift-perpendicular profiles, amounts of cumulative extension are determined, and by integrating four new Ar-40/Ar-39 radiometric dates for the Silali volcano into the existing geochronology of the faulted volcanic units, time-averaged extension rates are calculated. This study reveals that in the inner graben of the NKR, the long-term extension rate based on mid-Pleistocene to recent brittle deformation has minimum values of 1.0-1.6 mm yr(-1), locally with values up to 2.0 mm yr(-1). A comparison with the decadal, geodetically determined extension rate reveals that at least 65% of the extension must be accommodated within a narrow, 20-km-wide zone of the inner rift. In light of virtually inactive border faults of the NKR, we show that extension is focused in the region of the active volcano-tectonic axis in the inner graben, thus highlighting the maturing of continental rifting in the NKR.
The structural evolution of calderas in rifts helps to characterize the spatiotemporal relationships between magmatism, long wavelength crustal deformation and the formation of tectonic deformation zones along the rift axis. We document the structural characteristics of the c. 36 ka old Menengai Caldera located within a young zone of extension in the central Kenya Rift. Field mapping and high-resolution digital surface models show that NNE-striking Holocene normal faults perpendicular to the regional ESE-WNWextension direction dominate the interior sectors of the rift. Inside the caldera, these structures are overprinted by post-collapse doming and faulting of the magmatic centre, resulting in obliquely slipping normal faults bounding a resurgence horst. Radiocarbon dating of faulted units as young as 5 ka cal BP and the palaeo-shorelines of a lake formed during the African Humid Period in the Nakuru Basin indicate that volcanism and fault activity inside and in the vicinity of Menengai must have been sustained during the Holocene. Our analysis confirms that the caldera is located at the centre of an extending rift segment and suggests that other magmatic centres and young zones of faulting along the volcano-tectonic axis of the Kenya Rift may constitute nucleation points of faulting that ultimately foster future continental break-up.
Complex, time-dependent, and asymmetric rift geometries are observed throughout the East African Rift System (EARS) and are well documented, for instance, in the Kenya Rift. To unravel asymmetric rifting processes in this region, we conduct 2D geodynamic models. We use the finite element software ASPECT employing visco-plastic rheologies, mesh-refinement, distributed random noise seeding, and a free surface. In contrast to many previous numerical modeling studies that aimed at understanding final rifted margin symmetry, we explicitly focus on initial rifting stages to assess geodynamic controls on strain localization and fault evolution. We thereby link to geological and geophysical observations from the Southern and Central Kenya Rift. Our models suggest a three-stage early rift evolution that dynamically bridges previously inferred fault-configuration phases of the eastern EARS branch: (1) accommodation of initial strain localization by a single border fault and flexure of the hanging-wall crust, (2) faulting in the hanging-wall and increasing upper-crustal faulting in the rift-basin center, and (3) loss of pronounced early stage asymmetry prior to basinward localization of deformation. This evolution may provide a template for understanding early extensional faulting in other branches of the East African Rift and in asymmetric rifts worldwide. By modifying the initial random noise distribution that approximates small-scale tectonic inheritance, we show that a spectrum of first-order fault configurations with variable symmetry can be produced in models with an otherwise identical setup. This approach sheds new light on along-strike rift variability controls in active asymmetric rifts and proximal rifted margins.
This paper is mainly based on field work carried out on the Messinian deposits of the Adana Basin ( southern Turkey), as well as on the interpretation of seismic reflection profiles to understand 3D geometries of the basin fill. Chronostratigraphic constraints for the Messinian deposits are from micropaleontological studies on foraminifera, ostracods, and calcareous nannofossils, recently carried out on the Messinian deposits of the Adana Basin. Our results indicate that this basin developed in a marginal area strictly related to the Mediterranean realm. The Messinian deposits of the Adana Basin record all the main steps of the Messinian Salinity Crisis ( MSC) that affected the Mediterranean area at the end of the Miocene. The new stratigraphic model for the Messinian deposits of the Adana Basin provided in this work gives new insights into both the MSC and the Taurus petroleum system. Despite their complete correspondence with the MSC, the Messinian deposits of the Adana Basin show some differences with respect to the current conceptual model for the MSC. For example, in the current conceptual model for the MSC, only one regional erosional surface ( MES) characterizes the MSC deposits. In the Adana Basin, two regional erosional surfaces, named MES1 and MES2, separate the Messinian deposits related to the MSC in Lower Evaporites, Resedimented Lower Evaporites ( RLE), and upper Messinian continental deposits containing a late Lago-Mare ostracod assemblage ( mainly fluvial coarse-grained and fine-grained sediments). In some places, Brecciated Limestones lie just above the MES1 and beneath the RLE. In addition, the RLE are thought to be related to the same step that brought to the Messinian halite deposition throughout the Mediterranean, pointing to a hyperhaline environment. In contrast, the fine-grained deposits of the RLE of the Adana Basin show the occurrence of Parathetyan brackish ostracod fauna ( early Lago-Mare ostracod assemblages), which defines an oligohaline depositional environment for the RLE. In terms of hydrocarbon prospecting, the Messinian evaporites of the Adana Basin have been considered as a perfect seal for the active Taurus petroleum system. Our results show that due to the complex stratigraphic architecture of the basin fill and the occurrence of two regional erosional surfaces ( MES1 and MES2), the Messinian evaporites are discontinuously present both in surface and in the subsurface of the Adana Basin. However, seal properties in the Adana Basin could be found in the Lower Pliocene deep marine clays of the Avadan Formation. This work leads to suggest a new stratigraphical model for the Messinian deposits of the Adana Basin, allowing us to amend the classical scheme with respect to the Messinian, and to officially define some new formations within the stratigraphy of the Adana Basin.
The intermontane Quebrada de Humahuaca Basin (Humahuaca Basin) in the Eastern Cordillera of the southern Central Andes of NW Argentina (23 degrees-24 degrees S) records the evolution of a formerly contiguous foreland-basin setting to an intermontane depositional environment during the late stages of Cenozoic Andean mountain building. This basin has been and continues to be subject to shortening and surface uplift, which has resulted in the establishment of an orographic barrier for easterly sourced moisture-bearing winds along its eastern margin, followed by leeward aridification. We present new U-Pb zircon ages and palaeocurrent reconstructions suggesting that from at least 6Ma until 4.2Ma, the Humahuaca Basin was an integral part of a largely contiguous depositional system that became progressively decoupled from the foreland as deformation migrated eastward. The Humahuaca Basin experienced multiple cycles of severed hydrological conditions and subsequent re-captured drainage, fluvial connectivity with the foreland and sediment evacuation. Depositional and structural relationships among faults, regional unconformities and deformed landforms reveal a general pattern of intrabasin deformation that appears to be associated with different cycles of alluviation and basin excavation in which deformation is focused on basin-internal structures during or subsequent to phases of large-scale sediment removal.
The intermontane Humahuaca Basin in the Eastern Cordillera of the northwest Argentine Andes lies leeward of an orographic barrier to easterly derived moisture. An average of >2000 mm/yr of rainfall along the eastern flanks of the barrier contrasts with <200 mm/yr in the orogen interior. Paleoenvironmental reconstructions suggest that the basin became disconnected from the foreland during the Miocene-Pliocene by the growth of fault-bounded mountain ranges. Fossil records, sedimentology, and stable isotope data imply that rerouting of the fluvial network by 4.2 Ma and reduced rainfall by ca. 3 Ma were consequences of that range uplift. Here, we present cosmogenic nuclide-derived (Be-10) paleodenudation rates from 6 to 2 Ma fluvial deposits collected from the Humahuaca Basin. Despite increased tectonic activity, our Be-10 data show a tenfold decrease in denudation rates at ca. 3 Ma, documenting a link between uplift-induced semiarid conditions and decreasing hillslope denudation rates. This new data set thus demonstrates the influence of hydrological change on spatiotemporal denudation patterns in tectonically active mountain areas.
Stable-isotopic and sedimentary records from the orogenic Puna Plateau of NW Argentina and adjacent intermontane basins to the east furnish a unique late Cenozoic record of range uplift and ensuing paleoenvironmental change in the south-central Andes. Today, focused precipitation in this region occurs along the eastern, windward flanks of the Eastern Cordillera and Sierras Pampeanas ranges, while the orogen interior constitutes high-elevation regions with increasingly arid conditions in a westward direction. As in many mountain belts, such hydrologic and topographic gradients are commonly mirrored by a systematic relationship between the oxygen and hydrogen stable isotope ratios of meteoric water and elevation. The glass fraction of isotopically datable volcanic ash intercalated in sedimentary sequences constitutes an environmental proxy that retains a signal of the hydrogen-isotopic composition of ancient precipitation. This isotopic composition thus helps to elucidate the combined climatic and tectonic processes associated with topographic growth, which ultimately controls the spatial patterns of precipitation in mountain belts. However, between 25.5 and 27 degrees S present-day river-based hydrogen isotope lapse rates are very low, possibly due to deep-convective seasonal storms that dominate runoff. If not accounted for, the effects of such conditions on moisture availability in the past may lead to misinterpretations of proxy-records of rainfall. Here, we present hydrogen-isotope data of volcanic glass (delta Dg), extracted from 34 volcanic ash layers in different sedimentary basins of the Eastern Cordillera and the Sierras Pampeanas. Combined with previously published delta Dg records and our refined U-Pb and (U-Th)/He zircon geochronology on 17 tuff samples, we demonstrate hydrogen-isotope variations associated with paleoenvironmental change in the Angastaco Basin, which evolved from a contiguous foreland to a fault-bounded intermontane basin during the late Mio-Pliocene. We unravel the environmental impact of Mio-Pliocene topographic growth and associated orographic effects on long-term hydrogen-isotope records of rainfall in the south-central Andes, and potentially identify temporal variations in regional isotopic lapse rates that may also apply to other regions with similar topographic boundary conditions. (C) 2016 Elsevier B.V. All rights reserved.
Sedimentary basin fills along the windward flanks of orogenic plateaus are valuable archives of paleoenvironmental change with the potential to resolve the history of surface uplift and orographic barrier formation. The intermontane basins of the southern Central Andes contain thick successions of sedimentary material that are commonly interbedded with datable volcanic ashes. We relate variations in the hydrogen isotopic composition of hydrated volcanic glass (delta D-g) of Neogene to Quaternary fills in the semiarid intermontane Humahuaca Basin (Eastern Cordillera, northwest Argentina) to spatiotemporal changes in topography and associated orographic effects. delta D values from volcanic glass in the basin strata (-117 parts per thousand to -98 parts per thousand) show two main trends that accompany observed tectonosedimentary events in the study area. Between 6.0 and 3.5 Ma, delta D-g values decrease by similar to 17 parts per thousand; this is associated with surface uplift in the catchment area. After 3.5 Ma, delta D-g values show abrupt deuterium enrichment, which we associate with (1) the attainment of threshold elevations for blocking moisture transport in the basin-bounding ranges to the east, and (2) the onset of semiarid conditions in the basin. Such orographic barriers throughout the eastern flanks of the Central Andes have impeded moisture transport into the orogen interior; this has likely helped maintain aridity and internal drainage conditions on the adjacent Andean Plateau.
The Andean Plateau of NW Argentina is a prominent example of a high-elevation orogenic plateau characterized by internal drainage, arid to hyper-arid climatic conditions and a compressional basin-and-range morphology comprising thick sedimentary basins. However, the development of the plateau as a geomorphic entity is not well understood. Enhanced orographic rainout along the eastern, windward plateau flank causes reduced fluvial run-off and thus subdued surface-process rates in the arid hinterland. Despite this, many Puna basins document a complex history of fluvial processes that have transformed the landscape from aggrading basins with coalescing alluvial fans to the formation of multiple fluvial terraces that are now abandoned. Here, we present data from the San Antonio de los Cobres (SAC) area, a sub-catchment of the Salinas Grandes Basin located on the eastern Puna Plateau bordering the externally drained Eastern Cordillera. Our data include: (a) new radiometric U-Pb zircon data from intercalated volcanic ash layers and detrital zircons from sedimentary key horizons; (b) sedimentary and geochemical provenance indicators; (c) river profile analysis; and (d) palaeo-landscape reconstruction to assess aggradation, incision and basin connectivity. Our results suggest that the eastern Puna margin evolved from a structurally controlled intermontane basin during the Middle Miocene, similar to intermontane basins in the Mio-Pliocene Eastern Cordillera and the broken Andean foreland. Our refined basin stratigraphy implies that sedimentation continued during the Late Mio-Pliocene and the Quaternary, after which the SAC area was subjected to basin incision and excavation of the sedimentary fill. Because this incision is unrelated to baselevel changes and tectonic processes, and is similar in timing to the onset of basin fill and excavation cycles of intermontane basins in the adjacent Eastern Cordillera, we suspect a regional climatic driver, triggered by the Mid-Pleistocene Climate Transition, caused the present-day morphology. Our observations suggest that lateral orogenic growth, aridification of orogenic interiors, and protracted plateau sedimentation are all part of a complex process chain necessary to establish and maintain geomorphic characteristics of orogenic plateaus in tectonically active mountain belts.
The northern edge of the western central Tien Shan range is bounded by the Issyk-Ata fault situated south of Bishkek, the capital of Kyrgyzstan. Contraction in this thick-skinned orogen occurs with low-strain accumulation and long earthquake recurrence intervals. In the nineteenth to twentieth centuries, a sequence of large earthquakes with magnitudes between 6.9 and 8 affected the northern Tien Shan but left nearly the entire extent of the Issyk-Ata fault unruptured. Here, the only known historic earthquake ruptured in A.D. 1885 (M6.9) along the western end of the Issyk-Ata fault. Because earthquakes in low-strain regions often tend to cluster in time and may promote failure along nearby structures, the earthquake history of the northern Tien Shan represents an exceptional structural setting for studying fault behavior affected by an intraplate earthquake sequence. We present a paleoseismological study from one site (Belek) along the Issyk-Ata fault located east of the A.D. 1885 epicentral area. Our analysis combines a range of tools, including photogrammetry, differential Global Positioning System, 3D visualization, and age modeling with different dating methods (infrared stimulated luminescence, radiocarbon, U-series) to improve the reliability of an event chronology for the trench stratigraphy and fault geometry. We were able to distinguish three different surfacerupturing paleoearthquakes; these affected the area before 10.5 +/- 1.1 cal ka B.P., at similar to 5.6 +/- 1.0 cal ka B.P., and at similar to 630 +/- 100 cal B.P., respectively. Associated paleomagnitudes for the last two earthquakes range between M6.7 and 7.4, with a cumulative slip rate of 0.7 +/- 0.32 mm/yr. We did not find evidence for the A.D. 1885 event at Belek. Our study yielded two main overall results: first, it extends the regional historic and paleoseismic record; second, the documented rupture events along the Issyk-Ata fault suggest that this fault was not affected in its entirety; instead, these events indicate segmented rupture behavior.
New thermochronometric data from the Eastern Cordillera of the Colombian Andes reveal diachronous exhumation associated with Cenozoic contractional deformation in this sector of the northern Andes. We present a comprehensive account of exhumation patterns along a 150-km-long, across-strike transect between similar to 4 degrees and 6 degrees N by integrating 29 new apatite fission track (AFT) ages and 17 new zircon fission track (ZFT) ages with sparse published thermochronological data from this area. Our data reveal episodic eastward migration of the orogenic front at an average rate of 2.5-2.7 mm/a during the Late Cretaceous-Cenozoic. We identify three major stages of orogen propagation: (1) slow propagation (0.5-3.1 mm/a) until early Eocene; (2) rapid orogenic advance (4.0-18.0 mm/a) during middle-late Eocene, which accounts for similar to 86% of the orogen's present width; and (3) slow orogen propagation (1.2-2.1 mm/a) from Oligocene to Holocene times. Our data demonstrate that in the course of changes in plate kinematics, the presence of inherited crustal anisotropies, such as the former rift-bounding faults of the Eastern Cordillera, favor a nonsystematic progression of foreland basin deformation through time by preferentially concentrating accommodation of slip and thrust loading along these zones of weakness.
Foreland basin development in the Andes of central Colombia has been suggested to have started in the Late Cretaceous through tectonic loading of the Central Cordillera. Eastward migration of the Cenozoic orogenic front has also been inferred from the foreland basin record west of the Eastern Cordillera. However, farther east, limited data provided by foreland basin strata and the adjacent Eastern Cordillera complicate any correlation among mountain building, exhumation, and foreland basin sedimentation. In this study, we present new data from the Medina Basin in the eastern foothills of the Eastern Cordillera of Colombia. We report sedimentological data and palynological ages that link an eastward-thinning early Oligocene to early Miocene syntectonic wedge containing rapid facies changes with an episode of fast tectonic subsidence starting at ca. 31 Ma. This record may represent the first evidence of topographic loading generated by slip along the principal basement-bounding thrusts in the Eastern Cordillera to the southwest of the basin. Zircon fission-track ages and paleo-current analysis reveal the location of these thrust loads and illustrate a time lag between the sedimentary signal of topographic loading and the timing of exhumation (ca. 18 Ma). This lag may reflect the period between the onset of range uplift and significant removal of overburden. Vitrinite reflectance data document northward along-strike propagation of the deformation front and folding of the Oligocene syntectonic wedge. This deformation was coupled with a nonuniform incorporation of the basin into the wedge-top depozone. Thus, our data set constitutes unique evidence for the early growth and propagation of the deformation front in the Eastern Cordillera, which may also improve our understanding of spatiotemporal patterns of foreland evolution in other mountain belts.
The Tuz Golu Basin is the largest sedimentary depression located at the center of the Central Anatolian Plateau, an extensive, low-relief region with elevations of ca. 1 km located between the Pontide and Tauride mountains. Presently, the basin morphology and sedimentation processes are mainly controlled by the extensional Tuz Golu Fault Zone in the east and the transtensional Inonu-Eskisehir Fault System in the west. The purpose of this study is to contribute to the understanding of the Plio-Quaternary deformation history and to refine the timing of the latest extensional phase of the Tuz Golu Basin. Field observations, kinematic analyses, interpretations of seismic reflection lines, and Ar-40/Ar-39 dating of a key ignimbrite layer suggest that a regional phase of NNW-SSE to NE-SW contraction ended by 6.81 +/- 0.24 Ma and was followed by N-S to NE-SW extension during the Pliocene-Quaternary periods. Based on sedimentological and chronostratigraphic markers, the average vertical displacement rates over the past 5 or 3 Ma with respect to the central part of Tuz Golu Lake are 0.03 to 0.05 mm/year for the fault system at the western flank of the basin and 0.08 to 0.13 mm/year at the eastern flank. Paleo-shorelines of the Tuz Golu Lake, vestiges of higher lake levels related to Quaternary climate change, are important strain markers and were formed during Last Glacial Maximum conditions as indicated by a radiocarbon age of 21.8 +/- 0.4 ka BP obtained from a stromatolitic crust. Geomorphic observations and deformed lacustrine shorelines suggest that the main strand of the Tuz Golu Fault Zone straddling the foothills of the Sereflikochisar-Aksaray range has not been active during the Holocene. Instead, deformation appears to have migrated towards the interior of the basin along an offshore fault that runs immediately west of Sereflikochisar Peninsula. This basinward migration of deformation is probably associated with various processes acting at the lithospheric scale, such as plateau uplift and/or microplate extrusion.
Vegetation has long been hypothesized to influence the nature and rates of surface processes. We test the possible impact of vegetation and climate on denudation rates at orogen scale by taking advantage of a pronounced along-strike gradient in rainfall and vegetation density in the Himalaya. We combine 12 new 10Be denudation rates from the Sutlej Valley and 123 published denudation rates from fluvially- dominated catchments in the Himalaya with remotely-sensed measures of vegetation density and rainfall metrics, and with tectonic and lithologic constraints. In addition, we perform topographic analyses to assess the contribution of vegetation and climate in modulating denudation rates along strike. We observe variations in denudation rates and the relationship between denudation and topography along strike that are most strongly controlled by local rainfall amount and vegetation density, and cannot be explained by along-strike differences in tectonics or lithology. A W–E along-strike decrease in denudation rate variability positively correlates with the seasonality of vegetation density (R = 0.95, p < 0.05), and negatively correlates with mean vegetation density (R = −0.84, p < 0.05). Vegetation density modulates the topographic response to changing denudation rates, such that the functional relationship between denudation rate and topographic steepness becomes increasingly linear as vegetation density increases. We suggest that while tectonic processes locally control the pattern of denudation rates across strike of the Himalaya (i.e., S–N), along strike of the orogen (i.e., E–W) climate exerts a measurable influence on how denudation rates scatter around long-term, tectonically-controlled erosion, and on the functional relationship between topography and denudation
Vegetation has long been hypothesized to influence the nature and rates of surface processes. We test the possible impact of vegetation and climate on denudation rates at orogen scale by taking advantage of a pronounced along-strike gradient in rainfall and vegetation density in the Himalaya. We combine 12 new Be-10 denudation rates from the Sutlej Valley and 123 published denudation rates from fluvially-dominated catchments in the Himalaya with remotely-sensed measures of vegetation density and rainfall metrics, and with tectonic and lithologic constraints. In addition, we perform topographic analyses to assess the contribution of vegetation and climate in modulating denudation rates along strike. We observe variations in denudation rates and the relationship between denudation and topography along strike that are most strongly controlled by local rainfall amount and vegetation density, and cannot be explained by along-strike differences in tectonics or lithology. A W-E along-strike decrease in denudation rate variability positively correlates with the seasonality of vegetation density (R = 0.95, p < 0.05), and negatively correlates with mean vegetation density (R = -0.84, p < 0.05). Vegetation density modulates the topographic response to changing denudation rates, such that the functional relationship between denudation rate and topographic steepness becomes increasingly linear as vegetation density increases. We suggest that while tectonic processes locally control the pattern of denudation rates across strike of the Himalaya (i.e., S-N), along strike of the orogen (i.e., E-W) climate exerts a measurable influence on how denudation rates scatter around long-term, tectonically-controlled erosion, and on the functional relationship between topography and denudation. (C) 2016 Elsevier B.V. All rights reserved.
Understanding the rates and pattern of erosion is a key aspect of deciphering the impacts of climate and tectonics on landscape evolution. Denudation rates derived from terrestrial cosmogenic nuclides (TCNs) are commonly used to quantify erosion and bridge tectonic (Myr) and climatic (up to several kiloyears) time scales. However, how the processes of erosion in active orogens are ultimately reflected in Be-10 TCN samples remains a topic of discussion. We investigate this problem in the Arun Valley of eastern Nepal with 34 new Be-10-derived catchment-mean denudation rates. The Arun Valley is characterized by steep north-south gradients in topography and climate. Locally, denudation rates increase northward, from <0.2mmyr(-1) to similar to 1.5mmyr(-1) in tributary samples, while main stem samples appear to increase downstream from similar to 0.2mmyr(-1) at the border with Tibet to 0.91mmyr(-1) in the foreland. Denudation rates most strongly correlate with normalized channel steepness (R-2=0.67), which has been commonly interpreted to indicate tectonic activity. Significant downstream decrease of Be-10 concentration in the main stem Arun suggests that upstream sediment grains are fining to the point that they are operationally excluded from the processed sample. This results in Be-10 concentrations and denudation rates that do not uniformly represent the upstream catchment area. We observe strong impacts on Be-10 concentrations from local, nonfluvial geomorphic processes, such as glaciation and landsliding coinciding with areas of peak rainfall rates, pointing toward climatic modulation of predominantly tectonically driven denudation rates.
Elevation-dependent changes in n-alkane delta D and soil GDGTs across the South Central Andes
(2016)
Surface uplift of large plateaus may significantly influence regional climate and more specifically precipitation patterns and temperature, sometimes complicating paleoaltimetry interpretations. Thus, understanding the topographic evolution of tectonically active mountain belts benefits from continued development of reliable proxies to reduce uncertainties in paleoaltimetry reconstructions. Lipid biomarker-based proxies provide a novel approach to stable isotope paleoaltimetry and complement authigenic or pedogenic mineral proxy materials, in particular outside semi-arid climate zones where soil carbonates are not abundant but (soil) organic matter has a high preservation potential. Here we present delta D values of soil-derived n-alkanes and mean annual air temperature (MAT) estimates based on branched glycerol dialkyl glycerol tetraether (brGDGT) distributions to assess their potential for paleoelevation reconstructions in the southern central Andes. We analyzed soil samples across two environmental and hydrological gradients that include a hillslope (26-28 degrees S) and a valley (22-24 degrees S) transect on the windward flanks of Central Andean Eastern Cordillera in NW Argentina. Our results show that present-day n-alkane delta D values and brGDGT-based MAT estimates are both linearly related with elevation and in good agreement with present-day climate conditions. Soil n-alkanes show a delta D lapse rate (A(delta D)) of -1.64 parts per thousand/100 m (R-2 = 0.91, p < 0.01) at the hillslope transect, within the range of delta D lapse rates from precipitation and surface waters in other tropical regions in the Andes like the Eastern Cordillera in Colombia and Bolivia and the Equatorial and Peruvian Andes. BrGDGT-derived soil temperatures are similar to monitored winter temperatures in the region and show a lapse rate of Delta T = -0.51 degrees C/100 m (R-2 = 0.91, p < 0.01), comparable with lapse rates from in situ soil temperature measurements, satellite derived land-surface temperatures at this transect, and weather stations from the Eastern Cordillera at similar latitude. As a result of an increasing leeward sampling position along the valley transect lapse rates are biased towards lower values and display higher scatter (Delta(delta D) = -0.9 parts per thousand/100 m, R-2 = 0.76, p < 0.01 and Delta T = -0.19 degrees C/100 m, R-2 = 0.48, p < 0.05). Despite this higher complexity, they are in line with lapse rates from stream-water samples and in situ soil temperature measurements along the same transect. Our results demonstrate that both soil n-alkane delta D values and MAT reconstructions based on brGDGTs distributions from the hillslope transect (Delta(delta D) = -1.64 parts per thousand/100 m, R-2 = 0.91, p < 0.01 and Delta T = -0.51 degrees C/100 m, R-2 = 0.91, p < 0.01) track the direct effects of orography on precipitation and temperature and hence the combined effects of local and regional hydrology as well as elevation. (C) 2016 Elsevier B.V. All rights reserved.
Resolving Earth's surface at the meter scale is essential for an improved understanding of the dynamics of mass-movement processes. In this study, we explore the applicability and potential of digital elevation models (DEMs) derived from stereophotogrammetry to detect debris-flow channels in the Quebrada del Toro in the northwestern Argentine Andes. Our analysis relies on a high-resolution (3 m) DEM created from SPOT-7 tri-stereo satellite data. We carefully validated DEM quality with ∼6,000 differential GPS points and identified optimal parameters for DEM generation in high-relief terrain. After multiple processing steps, we achieved an accuracy of 0.051 ± 1.915 m (1σ) using n = 3,139 control points with cm precision. Previous studies have used the drainage area and slope framework to identify topographic signatures of debris flows within a catchment. We built upon this and investigated individual river-channel segments using connected-component (CC) analysis on meter-scale topographic data. We define CC as segments of similar slope along the channel profile. Based on seven manually identified debris-flow catchments, we developed a debris-flow similarity index using component length and mean channel-segment slope and identified channel segments that have likely been shaped by debris flows. The presented approach has the potential to resolve intra-catchment variability of transport processes, allows to constrain the extent of debris-flow channels more precisely than slope-area analysis, and highlights the versatility of combined space- and field-based observations for natural-hazard assessments.
The evolution of through-going normal-fault arrays from initial nucleation to growth and subsequent interaction and mechanical linkage is well documented in many extensional provinces. Over time, these processes. lead to predictable spatial and temporal variations in the amount and rate of displacement accumulated along strike of individual fault segments, which should be manifested in the patterns of footwall exhumation. Here, we investigate the along-strike and vertical distribution of low-temperature apatite (U-Th)/He (AHe) cooling ages along the bounding fault system, the Livingstone fault, of the Karonga Basin of the northern Malawi Rift. The fault evolution and linkage from rift initiation to the present day has been previously constrained through investigations of the hanging wall basin fill. The new cooling ages from the footwall of the Livingstone fault can be related to the adjacent depocentre evolution and across a relay zone between two palaeo-fault segments. Our data are complimented by published apatite fission track (AFT) data and reveal significant variation in rock cooling history along-strike: the centre of the footwall yields younger cooling ages than the former tips of earlier fault segments that are now linked. This suggests that low-temperature thermochronology can detect fault interactions along strike. That these former segment boundaries are preserved within exhumed footwall rocks is a function of the relatively recent linkage of the system. Our study highlights that changes in AHe (and potentially AFT) ages associated with the along-strike displacement profile can occur over relatively short horizontal distances (of a few kilometres). This is fundamentally important in the assessment of the vertical cooling history of footwalls in extensional systems: temporal differences in the rate of tectonically driven exhumation at a given location along fault strike may be of greater importance in controlling changes in rates of vertical exhumation than commonly invoked climatic fluctuations. (C) 2016 Elsevier B.V. All rights reserved.
The Cenozoic East African Rift System (EARS) is an exceptional example of active continental extension, providing opportunities for furthering our understanding of hydrocarbon plays within rifts. It is divided into structurally distinct western and eastern branches. The western branch comprises deep rift basins separated by transfer zones, commonly localised onto pre-existing structures, offering good regional scale hydrocarbon traps. At a basin-scale, local discrete inherited structures might also play an important role on fault localisation and hydrocarbon distribution. Here, we consider the evolution of the Central basin of the Malawi Rift, in particular the influence of pre-existing structural fabrics. Integrating basin-scale multichannel 2D, and high resolution seismic datasets we constrain the border, Mlowe-Nkhata, fault system (MNF) to the west of the basin and smaller Mbamba fault (MF) to the east and document their evolution. Intra basin structures define a series of horsts, which initiated as convergent transfers, along the basin axis. The horsts are offset along a NE SW striking transfer fault parallel to and along strike of the onshore Karoo (Permo-Triassic) Ruhuhu graben. Discrete pre-existing structures probably determined its location and, oriented obliquely to the extension orientation it accommodated predominantly strike-slip deformation, with more slowly accrued dip-slip. To the north of this transfer fault, the overall basin architecture is asymmetric, thickening to the west throughout; while to the south, an initially symmetric graben architecture became increasingly asymmetric in sediment distribution as strain localised onto the western MNF. The presence of the axial horst increasingly focussed sediment supply to the west. As the transfer fault increased its displacement, so this axial supply was interrupted, effectively starving the south-east while ponding sediments between the western horst margin and the transfer fault. This asymmetric bathymetry and partitioned sedimentation continues to the present-day, overprinting the early basin symmetry and configuration. Sediments deposited earlier become increasingly dissected and fault juxtapositions changed at a small (10-100 m) scale. The observed influence of basin-scale transfer faults on sediment dispersal and fault compartmentalization due to pre-existing structures oblique to the extension orientation is relevant to analogous exploration settings. (C) 2016 Elsevier Ltd. All rights reserved.
[ 1] The Eastern Cordillera of Colombia is key to understanding the role of inherited basement anisotropies in the evolution of active noncollisional mountain belts. In particular, the Rio Blanco-Guatiquia region of the Eastern Cordillera is exemplary in displaying a variety of phenomena that document the importance of the orientation, geometry, and segmentation of preorogenic anisotropies. We document the first unambiguous evidence that extensional basement structures played an important role in determining the locus of deformation during contractional reactivation in the Eastern Cordillera. Detailed structural field mapping and analysis of industry seismic reflection profiles have helped to identify the inherited San Juanito, Naranjal, and Servita normal faults and associated transfer faults as important structures that were inverted during the Cenozoic Andean orogeny. Apparently, the more internal faults in the former rift basin were not properly oriented for an efficient reactivation in contraction. However, these faults have a fundamental role as strain risers, as folding is concentrated west of them. In contrast, reactivated normal faults such as the more external Servita fault are responsible for uplifting the eastern flank of the Eastern Cordillera. In addition, these structures are adjacent and intimately linked to the development of thin-skinned faults farther east. In part, the superimposed compression in this prestrained extensional region is compensated by lateral escape. The dominant presence of basement involved buckling and thrusting, and the restricted development of thin-skinned thrusting in this inversion orogen makes the Eastern Cordillera a close analog to the intraplate Atlas Mountains of Morocco and other inverted sectors of the Andean orogen farther south
Lower Cretaceous early syn-rift facies along the eastern flank of the Eastern Cordillera of Colombia, their provenance, and structural context, reveal the complex interactions between Cretaceous extension, spatio-temporal trends in associated sedimentation, and subsequent inversion of the Cretaceous Guatiquia paleo-rift. South of 4 degrees 30'N lat, early syn-rift alluvial sequences in former extensional footwall areas were contemporaneous with fan- delta deposits in shallow marine environments in adjacent hanging-wall areas. In general, footwall erosion was more pronounced in the southern part of the paleorift. In contrast, early syn-rift sequences in former footwall areas in the northern rift sectors mainly comprise shallow marine supratidal sabkha to intertidal strata, whereas hanging-wall units display rapid transitions to open-sea shales. In comparison with the southern paleo-rift sector, fan-delta deposits in the north are scarce, and provenance suggests negligible footwall erosion. The southern graben segment had longer, and less numerous normal faults, whereas the northern graben segment was characterized by shorter, rectilinear faults. To the east, the graben system was bounded by major basin-margin faults with protracted activity and greater throw as compared with intrabasinal faults to the west. Intrabasinal structures grew through segment linkage and probably interacted kinematically with basin-margin faults. Basin-margin faults constitute a coherent fault system that was conditioned by pre-existing basement fabrics. Structural mapping, analysis of present-day topography, and balanced cross sections indicate that positive inversion of extensional structures was focused along basin-bounding faults, whereas intrabasinal faults remained unaffected and were passively transported by motion along the basin-bounding faults. Thus, zones of maximum subsidence in extension accommodated maximum elevation in contraction, and former topographic highs remained as elevated areas. This documents the role of basin-bounding faults as multiphased, long-lived features conditioned by basement discontinuities. Inversion of basin-bounding faults was more efficient in the southern than in the northern graben segment, possibly documenting the inheritance and pivotal role of fault-displacement gradients. Our observations highlight similarities between inversion features in orogenic belts and intra-plate basins, emphasizing the importance of the observed phenomena as predictive tools in the spatiotemporal analysis of inversion histories in orogens, as well as in hydrocarbon and mineral deposits exploration.
A valley-filling ignimbrite re-exposed through subsequent river incision at the southern margin of the Andean (Puna) plateau preserves pristine geological evidence of pre-late Miocene palaeotopography in the north western Argentine Andes. Our new Ar-40/(39) Ar dating of the Las Papas Ignimbrites yields a plateau age of 9.24 +/- 0.03 Ma, indicating valley-relief and orographic-barrier conditions comparable to the present-day. A later infill of Plio-Pleistocene coarse conglomerates has been linked to wetter conditions, but resulted in no additional net incision of the Las Papas valley, considering that the base of the ignimbrite remains unexposed in the valley bottom. Our observations indicate that at least 550 m of local plateau margin relief (and likely > 2 km) existed by 9 Ma at the southern Puna margin, which likely aided the efficiency of the orographic barrier to rainfall along the eastern and south eastern flanks of the Puna and causes aridity in the plateau interior.
The structural and topographic evolution of orogenic plateaus is an important research topic because of its impact on atmospheric circulation patterns, the amount and distribution of rainfall, and resulting changes in surface processes. The Puna region in the north-western Argentina (between 13 degrees S and 27 degrees S) is part of the Andean Plateau, which is the world's second largest orogenic plateau. In order to investigate the deformational events responsible for the initial growth of this part of the Andean plateau, we carried out structural and stratigraphic investigations within the present-day transition zone between the northern Puna and the adjacent Eastern Cordillera to the east. This transition zone is characterized by ubiquitous exposures of continental middle Eocene redbeds of the Casa Grande Formation. Our structural mapping, together with a sedimentological analysis of these units and their relationships with the adjacent mountain ranges, has revealed growth structures and unconformities that are indicative of syntectonic deposition. These findings support the notion that tectonic shortening in this part of the Central Andes was already active during the middle Paleogene, and that early Cenozoic deformation in the region that now constitutes the Puna occurred in a spatially irregular manner. The patterns of Paleogene deformation and uplift along the eastern margin of the present-day plateau correspond to an approximately north-south oriented swath of reactivated basement heterogeneities (i.e. zones of mechanical weakness) stemming from regional Paleozoic mountain building that may have led to local concentration of deformation belts.
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.
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.
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.
A comparison of deformation rates in active rifts over different temporal scales may help to decipher variations in their structural evolution, controlling mechanisms, and evolution of sedimentary environments through time. Here we use deformed lake shorelines in the Suguta and Turkana basins in northern Kenya as strain markers to estimate deformation rates at the 10(3)-10(4) yr time scale and compare them with rates spanning 10(1)-10(7) yr. Both basins are internally drained today, but until 7 to 5 kyr lake levels were 300 and 100 m higher, respectively, maintained by the elevation of overflow sills connecting them with the Nile drainage. Protracted high lake levels resulted in formation of a maximum highstand shoreline - a distinct geomorphic feature virtually continuous for several tens of kilometers. We surveyed the elevation of this geomorphic marker at 45 sites along >100 km of the rift, and use the overflow sills as vertical datum. Thin-shell elastic and thermomechanical models for this region predict up to similar to 10 m of rapid isostatic rebound associated with lake-level falls lasting until similar to 2 kyr ago. Holocene cumulative throw rates along four rift-normal profiles are 6.8-8.5 mm/yr, or 7.5-9.6 mm/yr if isostatic rebound is considered. Assuming fault dips of 55-65, inferred from seismic reflection profiles, we obtained extension rates of 3.2-6 mm/yr (including uncertainties in field measurements, fault dips, and ages), or 3.5-6.7 mm/yr considering rebound. Our estimates are consistent, within uncertainties, with extension rates of 4-5.1 mm/yr predicted by a modern plate-kinematic model and plate reconstructions since 3.2 Myr. The Holocene strain rate of 10(-15) s(-1) is similar to estimates on the similar to 10(6) yr scale, but over an order of magnitude higher than on the similar to 10(7) yr scale. This is coherent with continuous localization and narrowing of the plate boundary, implying that the lithospheric blocks limiting the Kenya Rift are relatively rigid. Increasing strain rate under steady extension rate suggests that, as the magnitude of extension and crustal thinning increases, the role of regional processes such as weakening by volcanism becomes dominant over far-field plate tectonics controlling the breakup process and the transition from continental rifting to oceanic spreading.
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.
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.
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.
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.
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.
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 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.
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 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.
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.
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.
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.
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.
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.