Refine
Year of publication
Document Type
- Doctoral Thesis (12)
- Article (8)
- Postprint (1)
Is part of the Bibliography
- yes (21) (remove)
Keywords
- Andes (21) (remove)
Institute
The removal, redistribution, and transient storage of sediments in tectonically active mountain belts is thought to exert a first-order control on shallow crustal stresses, fault activity, and hence on the spatiotemporal pattern of regional deformation processes. Accordingly, sediment loading and unloading cycles in intermontane sedimentary basins may inhibit or promote intrabasinal faulting, respectively, but unambiguous evidence for this potential link has been elusive so far. Here we combine 2D numerical experiments that simulate contractional deformation in a broken-foreland setting (i.e., a foreland where shortening is diachronously absorbed by spatially disparate, reverse faults uplifting basement blocks) with field data from intermontane basins in the NW Argentine Andes. Our modeling results suggest that thicker sedimentary fills (>0.7-1.0 km) may suppress basinal faulting processes, while thinner fills (<0.7 km) tend to delay faulting. Conversely, the removal of sedimentary loads via fluvial incision and basin excavation promotes renewed intrabasinal faulting. These results help to better understand the tectono-sedimentary history of intermontane basins that straddle the eastern border of the Andean Plateau in northwestern Argentina. For example, the Santa Maria and the Humahuaca basins record intrabasinal deformation during or after sediment unloading, while the Quebrada del Toro Basin reflects the suppression of intrabasinal faulting due to loading by coarse conglomerates. We conclude that sedimentary loading and unloading cycles may exert a fundamental control on spatiotemporal deformation patterns in intermontane basins of tectonically active broken forelands. (C) 2018 Elsevier B.V. All rights reserved.
Most hydrological studies rely on a model calibrated using discharge alone. However, judging the model reliability based on such calibration is problematic, as it does not guarantee the correct representation of internal hydrological processes. This study aims (a) to develop a comprehensive multi-objective calibration framework using remote sensing vegetation data and hydrological signatures (flow duration curve - FDC, and baseflow index) in addition to discharge, and (b) to apply this framework for calibration of the Soil and Water Assessment Tool (SWAT) in a typical Andean catchment. Overall, our calibration approach outperformed traditional discharge-based and FDC signature-based calibration strategies in terms of vegetation, streamflow, and flow partitioning simulation. New hydrological insights for the region are the following: baseflow is the main component of the streamflow sustaining the long dry-season flow, and pasture areas offer higher water yield and baseflow than other land-cover types. The proposed approach could be used in other data-scarce regions with complex topography.
Most hydrological studies rely on a model calibrated using discharge alone. However, judging the model reliability based on such calibration is problematic, as it does not guarantee the correct representation of internal hydrological processes. This study aims (a) to develop a comprehensive multi-objective calibration framework using remote sensing vegetation data and hydrological signatures (flow duration curve - FDC, and baseflow index) in addition to discharge, and (b) to apply this framework for calibration of the Soil and Water Assessment Tool (SWAT) in a typical Andean catchment. Overall, our calibration approach outperformed traditional discharge-based and FDC signature-based calibration strategies in terms of vegetation, streamflow, and flow partitioning simulation. New hydrological insights for the region are the following: baseflow is the main component of the streamflow sustaining the long dry-season flow, and pasture areas offer higher water yield and baseflow than other land-cover types. The proposed approach could be used in other data-scarce regions with complex topography.
Combined structural and magnetotelluric investigation across the West Fault Zone in northern Chile
(2002)
The characterisation of the internal architecture of large-scale fault zones is usually restricted to the outcrop-based investigation of fault-related structural damage on the Earth's surface. A method to obtain information on the downward continuation of a fault is to image the subsurface electrical conductivity structure. This work deals with such a combined investigation of a segment of the West Fault, which itself is a part of the more than 2000 km long trench-linked Precordilleran Fault System in the northern Chilean Andes. Activity on the fault system lasted from Eocene to Quaternary times. In the working area (22°04'S, 68°53'W), the West Fault exhibits a clearly defined surface trace with a constant strike over many tens of kilometers. Outcrop condition and morphology of the study area allow ideally for a combination of structural geology investigation and magnetotelluric (MT) / geomagnetic depth sounding (GDS) experiments. The aim was to achieve an understanding of the correlation of the two methods and to obtain a comprehensive view of the West Fault's internal architecture. Fault-related brittle damage elements (minor faults and slip-surfaces with or without striation) record prevalent strike-slip deformation on subvertically oriented shear planes. Dextral and sinistral slip events occurred within the fault zone and indicate reactivation of the fault system. Youngest deformation increments mapped in the working area are extensional and the findings suggest a different orientation of the extension axes on either side of the fault. Damage element density increases with approach to the fault trace and marks an approximately 1000 m wide damage zone around the fault. A region of profound alteration and comminution of rocks, about 400 m wide, is centered in the damage zone. Damage elements in this central part are predominantly dipping steeply towards the east (70-80°). Within the same study area, the electrical conductivity image of the subsurface was measured along a 4 km long MT/GDS profile. This main profile trends perpendicular to the West Fault trace. The MT stations of the central 2 km were 100 m apart from each other. A second profile with 300 m site spacing and 9 recording sites crosses the fault a few kilometers away from the main study area. Data were recorded in the frequency range from 1000 Hz to 0.001 Hz with four real time instruments S.P.A.M. MkIII. The GDS data reveal the fault zone for both profiles at frequencies above 1 Hz. Induction arrows indicate a zone of enhanced conductivity several hundred meters wide, that aligns along the WF strike and lies mainly on the eastern side of the surface trace. A dimensionality analysis of the MT data justifies a two dimensional model approximation of the data for the frequency range from 1000 Hz to 0.1 Hz. For this frequency range a regional geoelectric strike parallel to the West Fault trace could be recovered. The data subset allows for a resolution of the conductivity structure of the uppermost crust down to at least 5 km. Modelling of the MT data is based on an inversion algorithm developed by Mackie et al. (1997). The features of the resulting resistivity models are tested for their robustness using empirical sensitivity studies. This involves variation of the properties (geometry, conductivity) of the anomalies, the subsequent calculation of forward or constrained inversion models and check for consistency of the obtained model results with the data. A fault zone conductor is resolved on both MT profiles. The zones of enhanced conductivity are located to the east of the West Fault surface trace. On the dense MT profile, the conductive zone is confined to a width of about 300 m and the anomaly exhibits a steep dip towards the east (about 70°). Modelling implies that the conductivity increase reaches to a depth of at least 1100 m and indicates a depth extent of less than 2000 m. Further conductive features are imaged but their geometry is less well constrained. The fault zone conductors of both MT profiles coincide in position with the alteration zone. For the dense profile, the dip of the conductive anomaly and the dip of the damage elements of the central part of the fault zone correlate. This suggests that the electrical conductivity enhancement is causally related to a mesh of minor faults and fractures, which is a likely pathway for fluids. The interconnected rock-porosity that is necessary to explain the observed conductivity enhancement by means of fluids is estimated on the basis of the salinity of several ground water samples (Archie's Law). The deeper the source of the water sample, the more saline it is due to longer exposure to fluid-rock interaction and the lower is the fluid's resistivity. A rock porosity in the range of 0.8% - 4% would be required at a depth of 200 m. That indicates that fluids penetrating the damaged fault zone from close to the surface are sufficient to explain the conductivity anomalies. This is as well supported by the preserved geochemical signature of rock samples in the alteration zone. Late stage alteration processes were active in a low temperature regime (<95°C) and the involvement of ascending brines from greater depth is not indicated. The limited depth extent of the fault zone conductors is a likely result of sealing and cementation of the fault fracture mesh due to dissolution and precipitation of minerals at greater depth and increased temperature. Comparison of the results of the apparently inactive West Fault with published studies on the electrical conductivity structure of the currently active San Andreas Fault, suggests that the depth extent and conductivity of the fault zone conductor may be correlated to fault activity. Ongoing deformation will keep the fault/fracture mesh permeable for fluids and impede cementation and sealing of fluid pathways.
The Central Andes region in South America is characterized by a complex and heterogeneous deformation system. Recorded seismic activity and mapped neotectonic structures indicate that most of the intraplate deformation is located along the margins of the orogen, in the transitions to the foreland and the forearc. Furthermore, the actively deforming provinces of the foreland exhibit distinct deformation styles that vary along strike, as well as characteristic distributions of seismicity with depth. The style of deformation transitions from thin-skinned in the north to thick-skinned in the south, and the thickness of the seismogenic layer increases to the south. Based on geological/geophysical observations and numerical modelling, the most commonly invoked causes for the observed heterogeneity are the variations in sediment thickness and composition, the presence of inherited structures, and changes in the dip of the subducting Nazca plate. However, there are still no comprehensive investigations on the relationship between the lithospheric composition of the Central Andes, its rheological state and the observed deformation processes. The central aim of this dissertation is therefore to explore the link between the nature of the lithosphere in the region and the location of active deformation. The study of the lithospheric composition by means of independent-data integration establishes a strong base to assess the thermal and rheological state of the Central Andes and its adjacent lowlands, which alternatively provide new foundations to understand the complex deformation of the region. In this line, the general workflow of the dissertation consists in the construction of a 3D data-derived and gravity-constrained density model of the Central Andean lithosphere, followed by the simulation of the steady-state conductive thermal field and the calculation of strength distribution. Additionally, the dynamic response of the orogen-foreland system to intraplate compression is evaluated by means of 3D geodynamic modelling.
The results of the modelling approach suggest that the inherited heterogeneous composition of the lithosphere controls the present-day thermal and rheological state of the Central Andes, which in turn influence the location and depth of active deformation processes. Most of the seismic activity and neo--tectonic structures are spatially correlated to regions of modelled high strength gradients, in the transition from the felsic, hot and weak orogenic lithosphere to the more mafic, cooler and stronger lithosphere beneath the forearc and the foreland. Moreover, the results of the dynamic simulation show a strong localization of deviatoric strain rate second invariants in the same region suggesting that shortening is accommodated at the transition zones between weak and strong domains. The vertical distribution of seismic activity appears to be influenced by the rheological state of the lithosphere as well. The depth at which the frequency distribution of hypocenters starts to decrease in the different morphotectonic units correlates with the position of the modelled brittle-ductile transitions; accordingly, a fraction of the seismic activity is located within the ductile part of the crust. An exhaustive analysis shows that practically all the seismicity in the region is restricted above the 600°C isotherm, in coincidence with the upper temperature limit for brittle behavior of olivine. Therefore, the occurrence of earthquakes below the modelled brittle-ductile could be explained by the presence of strong residual mafic rocks from past tectonic events. Another potential cause of deep earthquakes is the existence of inherited shear zones in which brittle behavior is favored through a decrease in the friction coefficient. This hypothesis is particularly suitable for the broken foreland provinces of the Santa Barbara System and the Pampean Ranges, where geological studies indicate successive reactivation of structures through time. Particularly in the Santa Barbara System, the results indicate that both mafic rocks and a reduction in friction are required to account for the observed deep seismic events.
Sudden glacier advances in the Cachapoal Valley, Southern Central Andes of Chile (34 degrees S)
(2021)
Throughout the Andes Mountains of South America, a general trend of glacier shrinkage has taken place in modern times. However, a few glaciers have undergone considerable temporally advances or even surged during the mid-19th to 20th century CE. These valley glaciers are mainly located in the Central Andes of Chile and Argentina. The research presented here focuses on the changes of the Cachapoal Glacier in the Southern Central Andes of Chile. Spectacular glacier advances occurred at least three times in historical times, which lead to river blockages and successive lake outburst floods. The glacier advances were reconstructed with a multi-method approach including geomorphological mapping, Be-10 cosmogenic exposure dating of moraines, multi-temporal comparison of historical and recent photographs and paintings as well as the interpretation of aerial photographs and satellite images and the analysis of early travel reports. The article highlights the diversity of environmental conditions for the formation of glaciers in terms of the topographical and climatic setting and the resulting distinct glacier behavior along the Andes Mountains. It is argued for the Cachapoal Glacier that the glacier advances are intrinsic to the glacier type and may not be necessarily climate-dependent. This is characteristic for avalanche-fed glaciers of which the glacier dynamic is strongly controlled by the topographic setting and sudden inputs of ice and rock avalanches as well as by the specific debris transfer system and hydrological drainage pattern. At the regional level, the fluctuations of the Cachapoal Glacier are compared with glaciers of neighboring mountain ranges in the Southern Central Andes and at the global scale with those of the Karakoram Mountains in High Asia with a similar dynamic glacier behavior.
Mountain ranges can fundamentally influence the physical and and chemical processes that shape Earths’ surface. With elevations of up to several kilometers they create climatic enclaves by interacting with atmospheric circulation and hydrologic systems, thus leading to a specific distribution of flora and fauna. As a result, the interiors of many Cenozoic mountain ranges are characterized by an arid climate, internally drained and sediment-filled basins, as well as unique ecosystems that are isolated from the adjacent humid, low-elevation regions along their flanks and forelands. These high-altitude interiors of orogens are often characterized by low relief and coalesced sedimentary basins, commonly referred to as plateaus, tectono-geomorphic entities that result from the complex interactions between mantle-driven geological and tectonic conditions and superposed atmospheric and hydrological processes. The efficiency of these processes and the fate of orogenic plateaus is therefore closely tied to the balance of constructive and destructive processes – tectonic uplift and erosion, respectively. In numerous geological studies it has been shown that mountain ranges are delicate systems that can be obliterated by an imbalance of these underlying forces. As such, Cenozoic mountain ranges might not persist on long geological timescales and will be destroyed by erosion or tectonic collapse. Advancing headward erosion of river systems that drain the flanks of the orogen may ultimately sever the internal drainage conditions and the maintenance of storage of sediments within the plateau, leading to destruction of plateau morphology and connectivity with the foreland. Orogenic collapse may be associated with the changeover from a compressional stress field with regional shortening and topographic growth, to a tensional stress field with regional extensional deformation and ensuing incision of the plateau. While the latter case is well-expressed by active extensional faults in the interior parts of the Tibetan Plateau and the Himalaya, for example, the former has been attributed to have breached the internally drained areas of the high-elevation sectors of the Iranian Plateau.
In the case of the Andes of South America and their internally drained Altiplano-Puna Plateau, signs of both processes have been previously described. However, in the orogenic collapse scenario the nature of the extensional structures had been primarily investigated in the northern and southern terminations of the plateau; in some cases, the extensional faults were even regarded to be inactive. After a shallow earthquake in 2020 within the Eastern Cordillera of Argentina that was associated with extensional deformation, the state of active deformation and the character of the stress field in the central parts of the plateau received renewed interest to explain a series of extensional structures in the northernmost sectors of the plateau in north-western Argentina. This study addresses (1) the issue of tectonic orogenic collapse of the Andes and the destruction of plateau morphology by studying the fill and erosion history of the central eastern Andean Plateau using sedimentological and geochronological data and (2) the kinematics, timing and magnitude of extensional structures that form well-expressed fault scarps in sediments of the regional San Juan del Oro surface, which is an integral part of the Andean Plateau and adjacent morphotectonic provinces to the east.
Importantly, sediment properties and depositional ages document that the San Juan del Oro Surface was not part of the internally-drained Andean Plateau, but rather associated with a foreland-directed drainage system, which was modified by the Andean orogeny and that became successively incorporated into the orogen by the eastward-migration of the Andean deformation front during late Miocene – Pliocene time. Structural and geomorphic observations within the plateau indicate that extensional processes must have been repeatedly active between the late Miocene and Holocene supporting the notion of plateau-wide extensional processes, potentially associated with Mw ~ 7 earthquakes. The close relationship between extensional joints and fault orientations underscores that 3 was oriented horizontally in NW-SE direction and 1 was vertical. This unambiguously documents that the observed deformation is related to gravitational forces that drive the orogenic collapse of the plateau. Applied geochronological analyses suggest that normal faulting in the northern Puna was active at about 3 Ma, based on paired cosmogenic nuclide dating of sediment fill units. Possibly due to regional normal faulting the drainage system within the plateau was modified, promoting fluvial incision.
Subduction zones are regions of intense earthquake activity up to great depth. Sources are located inside the subducting lithosphere and, as a consequence, seismic radiation from subduction zone earthquakes is strongly affected by the interior slab structure. The wave field of these intraslab events observed in the forearc region is profoundly influenced by a seismically slow layer atop the slab surface. This several kilometer thick low-velocity channel (wave guide) causes the entrapment of seismic energy producing strong guided wave phases that appear in P onsets in certain regions of the forearc. Observations at the Chile-Peru subduction zone presented here, as well as observations at several other circum-pacific subduction zones show such signals. Guided wave analysis contributes details of immense value regarding the processes near the slab surface, such as layering of subducted lithosphere, source locations of intraslab seismicity and most of all, range and manner of mineralogical phase transitions. Seismological data stem from intermediate depth events (depth range 70 km - 300 km) recorded in northern Chile near 21 Grad S during the collaborative research initiative " Deformation Processes in the Andes" (SFB 267). A subset of stations - all located within a slab-parallel transect close to 69 Grad W - show low-frequency first arrivals (2 Hz), sometimes followed by a second high-frequency phase. We employ 2-dimensional finite-difference simulations of complete P-SV wave propagation to explore the parameter space of subduction zone wave guides and explain the observations. Key processes underlying the guided wave propagation are studied: Two distinct mechanisms of decoupling of trapped energy from the wave guide are analyzed - a prerequisite to observe the phases at stations located at large distances from the wave guide (up to 100 km). Variations of guided wave effects perpendicular to the strike of the subduction zone are investigated, such as the influence of phases traveling in the fast slab. Further, the merits and limits of guided wave analysis are assessed. Frequency spectra of the guided wave onsets prove to be a robust quantity that captures guided wave characteristics at subduction zones including higher mode excitation. They facilitate the inference of wave guide structure and source positioning: The peak frequency of the guided wave fundamental mode is associated with a certain combination of layer width and velocity contrast. The excitation strength of the guided wave fundamental mode and higher modes is associated with source position and orientation relative to the low-velocity layer. The guided wave signals at the Chile-Peru subduction zone are caused by energy that leaks from the subduction zone wave guide. On the one hand, the bend shape of the slab allows for leakage at a depth of 100 km. On the other, equalization of velocities between the wave guide and the host rocks causes further energy leakage at the contact zone between continental and oceanic crust (70 km depth). Guided waves bearing information on deep slab structure can therefore be recorded at specific regions in the forearc. These regions are determined based on slab geometry, and their locations coincide with the observations. A number of strong constraints on the structure of the Chile-Peru slab are inferred: The deep wave guide for intraslab events is formed by a layer of 2 km average width that remains seismically slow (7 percent velocity reduction compared to surrounding mantle). This low-velocity layer at the top of the Chile-Peru slab is imaged from a depth of 100 km down to at least 160 km. Intermediate depth events causing the observed phases are located inside the layer or directly beneath it in the slab mantle. The layer is interpreted as partially eclogized lower oceanic crust persisting to depth beyond the volcanic arc.
Understanding the influence of climate change and population pressure on human conflict remains a critically important topic in the social sciences. Long-term records that evaluate these dynamics across multiple centuries and outside the range of modern climatic variation are especially capable of elucidating the relative effect of-and the interaction between-climate and demography. This is crucial given that climate change may structure population growth and carrying capacity, while both climate and population influence per capita resource availability. This study couples paleoclimatic and demographic data with osteological evaluations of lethal trauma from 149 directly accelerator mass spectrometry C-14-dated individuals from the Nasca highland region of Peru. Multiple local and supraregional precipitation proxies are combined with a summed probability distribution of 149 C-14 dates to estimate population dynamics during a 700-y study window. Counter to previous findings, our analysis reveals a precipitous increase in violent deaths associated with a period of productive and stable climate, but volatile population dynamics. We conclude that favorable local climate conditions fostered population growth that put pressure on the marginal and highly circumscribed resource base, resulting in violent resource competition that manifested in over 450 y of internecine warfare. These findings help support a general theory of intergroup violence, indicating that relative resource scarcity-whether driven by reduced resource abundance or increased competition-can lead to violence in subsistence societies when the outcome is lower per capita resource availability.