Refine
Year of publication
Language
- English (51)
Is part of the Bibliography
- yes (51) (remove)
Keywords
- Chile (5)
- Marine terraces (3)
- Maule earthquake (3)
- TerraceM (3)
- African Humid Period (2)
- Central Chile (2)
- Coastal uplift (2)
- GPS (2)
- LiDAR (2)
- Maule (2)
- earthquake (2)
- extensional tectonics (2)
- marine terraces (2)
- megathrust earthquake (2)
- postseismic deformation (2)
- 2010 M8.8 maule (1)
- Arauco Bay (1)
- Baragoi (1)
- Be-10 (1)
- Black Sea (1)
- Cenozoic (1)
- Central Anatolia (1)
- Central Anatolian Plateau (1)
- Central Pontides (1)
- Chile convergent margin (1)
- Chile earthquake (1)
- Coastal sedimentation (1)
- Coastal sequence of shorelines (1)
- Continental neotectonics (1)
- Coseismic coastal uplift (1)
- Coseismic slip distribution (1)
- DEM analysis (1)
- Earthquake (1)
- Earthquake hazards (1)
- East Africa (1)
- East African Rift System (1)
- Ethiopia rift (1)
- FEM (1)
- FEM models (1)
- Fluvial Incision (1)
- Geodynamic (1)
- Geomorphology (1)
- Holocene (1)
- Holocene extension (1)
- Intertidal organisms (1)
- Kenya Rift (1)
- Lake Turkana (1)
- M(w)8.8 Maule earthquake (1)
- Middle Holocene (1)
- New-Zealand (1)
- Orogenic Plateaus (1)
- Orogenic wedges (1)
- Palaeo-shorelines (1)
- Permanent uplift (1)
- Perumytilus mussels (1)
- Plateau margins (1)
- Rasa (1)
- Sea-level change (1)
- Seismic cycle (1)
- Seismic gap (1)
- Seismotectonic segmentation (1)
- Seismotektonik (1)
- Seismoturbidites (1)
- Slip deficit (1)
- Slip distribution (1)
- Strandlines (1)
- Subduction earthquakes (1)
- Surface Exposure Age (1)
- TanDEM-X DEM (1)
- Tectonic Geomorphology (1)
- Tectonic deformation (1)
- Tectonic uplift (1)
- Tectonics and climatic interactions (1)
- Terrace reoccupation (1)
- Tsunami effect (1)
- Turbidite paleoseismology (1)
- Tuz Golu Basin (1)
- Uplift Rate (1)
- Uplift rate (1)
- Wave-built terraces (1)
- afterslip (1)
- andean subduction zone (1)
- coastal geomorphology (1)
- conservation (1)
- continental rifting (1)
- continuous GPS measurements (1)
- crustal fault (1)
- cycle (1)
- disturbance (1)
- earthquake cycle (1)
- earthquake deformation (1)
- erosion (1)
- faulting behavior (1)
- finite-element model (1)
- fore arc (1)
- geochronology (1)
- geomorphic markers (1)
- habitats (1)
- isostatic rebound (1)
- kinematic analysis (1)
- knickpoints (1)
- lake shoreline (1)
- lake shorelines (1)
- locking degree (1)
- megathrust (1)
- morphometric indices (1)
- morphometry (1)
- neogene (1)
- neotectonics (1)
- normal faults (1)
- northern Kenya Rift (1)
- orogenic plateau evolution (1)
- paleo-delta (1)
- protection (1)
- renya Rift (1)
- resolution (1)
- rift transfer zone (1)
- salt spray (1)
- sandy-beach (1)
- sea level change (1)
- seismotectonic segmentation (1)
- seismotectonics (1)
- slip (1)
- slip model (1)
- south-central Chile (1)
- splay fault (1)
- splay faults (1)
- subduction zone earthquakes (1)
- tectonic geomorphology (1)
- tsunami hazards (1)
- vegetation (1)
- zonation (1)
Institute
The Gofa Province and the Chew Bahir Basin of southern Ethiopia constitute tectonically active regions, where the Southern Main Ethiopian Rift converges with the Northern Kenya Rift through a wide zone of extensional deformation with several north to northeast-trending, left-stepping en-e & PRIME;chelon basins. This sector of the Southern Main Ethiopian Rift is characterized by a semi-arid climate and a largely uniform lithology, and thus provides ideal conditions for studying the different parameters that define the tectonic and geomorphic features of this complex kinematic transfer zone. In this study, the degree of tectonic activity, spatiotemporal variations in extension, and the nature of kinematic linkage between different fault systems of the transfer zone are constrained by detailed quantitative geomorphic analysis of river catchments and focused field work. We analyzed fluvial and landscape morphometric characteristics in combination with structural, seismicity, and climatic data to better evaluate the tectono-geomorphic history of this transfer zone. Our data reveal significant north-south variations in the degree of extension from the Sawula Basin in the north (mature) to the Chew Bahir Basin in the south (juvenile). First, normalized channel-steepness indices and the spatial arrangement of knickpoints in footwall-draining streams suggest a gradual, southward shift in extensional deformation and recent tectonic activity. Second, based on 1-k(m) radius local relief and mean-hillslope maximum values that are consistent with ksn anomalies, we confirm strain localization within zones of fault interaction. Third, morphometric indices such as hypsometry, basin asymmetry factor, and valley floor width to valley height ratio also indicate a north to south gradient in tectonic activity, highlighting the importance of such a wide transfer zone with diffuse extension linking different rift segments during the break-up of continental crust.
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 first step towards assessing hazards in seismically active regions involves mapping capable faults and estimating their recurrence times. While the mapping of active faults is commonly based on distinct geologic and geomorphic features evident at the surface, mapping blind seismogenic faults is complicated by the absence of on-fault diagnostic features. Here we investigated the Pichilemu Fault in coastal Chile, unknown until it generated a Mw 7.0 earthquake in 2010. The lack of evident surface faulting suggests activity along a partly-hidden blind fault. We used off-fault deformed marine terraces to estimate a fault-slip rate of 0.52 ± 0.04 m/ka, which, when integrated with satellite geodesy suggests a 2.12 ± 0.2 ka recurrence time for Mw~7.0 normal-faulting earthquakes. We propose that extension in the Pichilemu region is associated with stress changes during megathrust earthquakes and accommodated by sporadic slip during upper-plate earthquakes, which has implications for assessing the seismic potential of cryptic faults along convergent margins and elsewhere.
The first step towards assessing hazards in seismically active regions involves mapping capable faults and estimating their recurrence times. While the mapping of active faults is commonly based on distinct geologic and geomorphic features evident at the surface, mapping blind seismogenic faults is complicated by the absence of on-fault diagnostic features. Here we investigated the Pichilemu Fault in coastal Chile, unknown until it generated a Mw 7.0 earthquake in 2010. The lack of evident surface faulting suggests activity along a partly-hidden blind fault. We used off-fault deformed marine terraces to estimate a fault-slip rate of 0.52 ± 0.04 m/ka, which, when integrated with satellite geodesy suggests a 2.12 ± 0.2 ka recurrence time for Mw~7.0 normal-faulting earthquakes. We propose that extension in the Pichilemu region is associated with stress changes during megathrust earthquakes and accommodated by sporadic slip during upper-plate earthquakes, which has implications for assessing the seismic potential of cryptic faults along convergent margins and elsewhere.
Tectonically active coasts are dynamic environments characterized by the presence of multiple marine terraces formed by the combined effects of wave erosion, tectonic uplift, and sea-level oscillations at glacialcycle timescales. Well-preserved erosional terraces from the last interglacial sea-level highstand are ideal marker horizons for reconstructing past sea-level positions and calculating vertical displacement rates. We carried out an almost continuous mapping of the last interglacial marine terrace along similar to 5000 km of the western coast of South America between 1 degrees N and 40 degrees S. We used quantitatively replicable approaches constrained by published terrace-age estimates to ultimately compare elevations and patterns of uplifted terraces with tectonic and climatic parameters in order to evaluate the controlling mechanisms for the formation and preservation of marine terraces and crustal deformation. Uncertainties were estimated on the basis of measurement errors and the distance from referencing points. Overall, our results indicate a median elevation of 30.1 m, which would imply a median uplift rate of 0.22 m kyr(-1) averaged over the past similar to 125 kyr. The patterns of terrace elevation and uplift rate display high-amplitude (similar to 100-200 m) and long-wavelength (similar to 10(2) km) structures at the Manta Peninsula (Ecuador), the San Juan de Marcona area (central Peru), and the Arauco Peninsula (south-central Chile). Medium-wavelength structures occur at the Mejillones Peninsula and Topocalma in Chile, while short-wavelength (< 10 km) features are for instance located near Los Vilos, Valparaiso, and Carranza, Chile. We interpret the long-wavelength deformation to be controlled by deep-seated processes at the plate interface such as the subduction of major bathymetric anomalies like the Nazca and Carnegie ridges. In contrast, short-wavelength deformation may be primarily controlled by sources in the upper plate such as crustal faulting, which, however, may also be associated with the subduction of topographically less pronounced bathymetric anomalies. Latitudinal differences in climate additionally control the formation and preservation of marine terraces. Based on our synopsis we propose that increasing wave height and tidal range result in enhanced erosion and morphologically well-defined marine terraces in south-central Chile. Our study emphasizes the importance of using systematic measurements and uniform, quantitative methodologies to characterize and correctly interpret marine terraces at regional scales, especially if they are used to unravel the tectonic and climatic forcing mechanisms of their formation. This database is an integral part of the World Atlas of Last Interglacial Shorelines (WALIS), published online at https://doi.org/10.5281/zenodo.4309748 (Freisleben et al., 2020).
The southern margin of the Central Anatolian Plateau (CAP) records a strong uplift phase after the early Middle Pleistocene, which has been related to the slab break-off of the subducting Arabian plate beneath the Anatolian microplate. During the last 450 kyr the area underwent an uplift phase at a mean rate of similar to 3.2 m/kyr, as suggested by Middle Pleistocene marine sediments exposed at similar to 1,500 m above sea level. These values are significantly higher than the 1.0-1.5 m/kyr estimated since the Late Pleistocene, suggesting temporal variations in uplift rate. To estimate changes in uplift rate during the Pleistocene we studied the marine terraces along the CAP southern margin, mapping the remnants of the platforms and their associated deposits in the field, and used the TerraceM software to identify the position and elevation of associated shoreline angles. We used shoreline angles and the timing of Quaternary marine sedimentation as constrains for a Landscape Evolution Model that simulates wave erosion of an uplifting coast. We applied random optimization algorithms and minimization statistics to find the input parameters that better reproduce the morphology of CAP marine terraces. The best-fitting uplift rate history suggests a significative increase from 1.9 to 3.5 m/kyr between 500 and 200 kyr, followed by an abrupt decrease to 1.4 m/kyr until the present. Our results agree with slab break-off models, which suggest a strong uplift pulse during slab rupture followed by a smoother decrease.
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.
Geomorphic strain markers accumulating the effects of many earthquake cycles help to constrain the mechanical behaviour of continental rift systems as well as the related seismic hazards. In the Corinth Rift (Greece), the unique record of onshore and offshore markers of Pleistocene similar to 100-ka climate cycles provides an outstanding possibility to constrain rift mechanics over a range of timescales. Here we use high-resolution topography to analyse the 3D geometry of a sequence of Pleistocene emerged marine terraces associated with flexural rift-flank uplift. We integrate this onshore dataset with offshore seismic data to provide a synoptic view of the flexural deformation across the rift. This allows us to derive an average slip rate of 4.5-9.0 mm.yr(-1) on the master fault over the past similar to 610 ka and an uplift/ subsidence ratio of 1:1.1-2.4. We reproduce the observed flexure patterns, using 3 and 5-layered lithospheric scale finite element models. Modelling results imply that the observed elastic flexure is produced by coseismic slip along 40-60 degrees planar normal faults in the elastic upper crust, followed by postseismic viscous relaxation occurring within the basal lower crust or upper mantle. We suggest that such a mechanism may typify rapid localised extension of continental lithosphere.
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.
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.
Great megathrust earthquakes arise from the sudden release of energy accumulated during centuries of interseismic plate convergence. The moment deficit (energy available for future earthquakes) is commonly inferred by integrating the rate of interseismic plate locking over the time since the previous great earthquake. But accurate integration requires knowledge of how interseismic plate locking changes decades after earthquakes, measurements not available for most great earthquakes. Here we reconstruct the post-earthquake history of plate locking at Guafo Island, above the seismogenic zone of the giant 1960 (M-w = 9.5) Chile earthquake, through forward modeling of land-level changes inferred from aerial imagery (since 1974) and measured by GPS (since 1994). We find that interseismic locking increased to similar to 70% in the decade following the 1960 earthquake and then gradually to 100% by 2005. Our findings illustrate the transient evolution of plate locking in Chile, and suggest a similarly complex evolution elsewhere, with implications for the time- and magnitude-dependent probability of future events.
Great megathrust earthquakes arise from the sudden release of energy accumulated during centuries of interseismic plate convergence. The moment deficit (energy available for future earthquakes) is commonly inferred by integrating the rate of interseismic plate locking over the time since the previous great earthquake. But accurate integration requires knowledge of how interseismic plate locking changes decades after earthquakes, measurements not available for most great earthquakes. Here we reconstruct the post-earthquake history of plate locking at Guafo Island, above the seismogenic zone of the giant 1960 (M-w = 9.5) Chile earthquake, through forward modeling of land-level changes inferred from aerial imagery (since 1974) and measured by GPS (since 1994). We find that interseismic locking increased to similar to 70% in the decade following the 1960 earthquake and then gradually to 100% by 2005. Our findings illustrate the transient evolution of plate locking in Chile, and suggest a similarly complex evolution elsewhere, with implications for the time- and magnitude-dependent probability of future events.
Postseismic surface deformation associated with great subduction earthquakes is controlled by asthenosphere rheology, frictional properties of the fault, and structural complexity. Here by modeling GPS displacements in the 6 years following the 2010 M-w 8.8 Maule earthquake in Chile, we investigate the impact of heterogeneous viscosity distribution in the South American subcontinental asthenosphere on the 3-D postseismic deformation pattern. The observed postseismic deformation is characterized by flexure of the South America plate with peak uplift in the Andean mountain range and subsidence in the hinterland. We find that, at the time scale of observation, over 2 orders of magnitude gradual increase in asthenosphere viscosity from the arc area toward the cratonic hinterland is needed to jointly explain horizontal and vertical displacements. Our findings present an efficient method to estimate spatial variations of viscosity, which clearly improves the fitting to the vertical signal of deformation. Lateral changes in asthenosphere viscosity can be correlated with the thermomechanical transition from weak subvolcanic arc mantle to strong subcratonic mantle, thus suggesting a stationary heterogeneous viscosity structure. However, we cannot rule out a transient viscosity structure (e.g., power law rheology) with the short time span of observation.
Along a subduction zone, great megathrust earthquakes recur either after long seismic gaps lasting several decades to centuries or over much shorter periods lasting hours to a few years when cascading successions of earthquakes rupture nearby segments of the fault. We analyze a decade of continuous Global Positioning System observations along the South American continent to estimate changes in deformation rates between the 2010 Maule (M8.8) and 2015 Illapel (M8.3) Chilean earthquakes. We find that surface velocities increased after the 2010 earthquake, in response to continental-scale viscoelastic mantle relaxation and to regional-scale increased degree of interplate locking. We propose that increased locking occurs transiently during a super-interseismic phase in segments adjacent to a megathrust rupture, responding to bending of both plates caused by coseismic slip and subsequent afterslip. Enhanced strain rates during a super-interseismic phase may therefore bring a megathrust segment closer to failure and possibly triggered the 2015 event.
The African Humid Period (AHP) between similar to 15 and 5.5 cal. kyr BP caused major environmental change in East Africa, including filling of the Suguta Valley in the northern Kenya Rift with an extensive (similar to 2150 km(2)), deep (similar to 300 m) lake. Interfingering fluvio-lacustrine deposits of the Baragoi paleo-delta provide insights into the lake-level history and how erosion rates changed during this time, as revealed by delta-volume estimates and the concentration of cosmogenic Be-10 in fluvial sand. Erosion rates derived from delta-volume estimates range from 0.019 to 0.03 mm yr(-1). Be-10-derived paleo-erosion rates at similar to 11.8 cal. kyr BP ranged from 0.035 to 0.086 mm yr(-1), and were 2.7 to 6.6 times faster than at present. In contrast, at similar to 8.7 cal. kyr BP, erosion rates were only 1.8 times faster than at present. Because Be-10-derived erosion rates integrate over several millennia; we modeled the erosion-rate history that best explains the 10Be data using established non-linear equations that describe in situ cosmogenic isotope production and decay. Two models with different temporal constraints (15-6.7 and 12-6.7 kyr) suggest erosion rates that were 25 to 300 times higher than the initial erosion rate (pre-delta formation). That pulse of high erosion rates was short (similar to 4 kyr or less) and must have been followed by a rapid decrease in rates while climate remained humid to reach the modern Be-10-based erosion rate of,similar to 0.013 mm yr(-1). Our simulations also flag the two highest Be-10-derived erosion rates at 11.8 kyr BP related to nonuniform catchment erosion. These changes in erosion rates and processes during the AHP may reflect a strong increase in precipitation, runoff, and erosivity at the arid-to-humid transition either at 15 or similar to 12 cal. kyr BP, before the landscape stabilized again, possibly due to increased soil production and denser vegetation.
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.
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.
Surface movements during the largest subduction zone earthquakes commonly drown coastlines. Yet, on geological timescales, coastlines above subduction zones uplift. Here I use a morphometric analysis combined with a numerical model of landscape evolution to estimate uplift rates along the central Andean rasa-a low-relief coastal surface bounded by a steep cliff formed by wave erosion. I find that the rasa has experienced steady uplift of 0.13 +/- 0.04 mm per year along a stretch of more than 2,000 km in length, during the Quaternary. These long-term uplift rates do not correlate with Global Positioning System (GPS) measurements of interseismic movements over the decadal scale, which implies that permanent uplift is not predominantly accumulated during the interseismic period. Instead, the rate of rasa uplift correlates with slip during earthquakes straddling the crust-mantle transition, the Moho. Such deeper earthquakes with magnitude 7 to 8 that occurred between 1995 and 2012 resulted in decimetres of coastal uplift. Slip during these earthquakes is located below the locked portion of the plate interface, and therefore may translate into permanent deformation of the overlying plate, where it causes uplift of the coastline. Thus, lower parts of the plate boundary are stably segmented over hundreds to millions of years. I suggest the coastline marks the surface expression of the transition between the shallow, locked seismogenic domain and the deeper, conditionally stable domain where modest earthquakes build up topography.
The subduction of bathymetric anomalies at convergent margins can profoundly affect subduction dynamics, magmatism, and the structural and geomorphic evolution of the overriding plate. The Northern Patagonian Icefield (NPI) is located east of the Chile Triple Junction at similar to 47 degrees S, where the Chile Rise spreading center collides with South America. This region is characterized by an abrupt increase in summit elevations and relief that has been controversially debated in the context of geodynamic versus glacial erosion effects on topography. Here we present geomorphic, thermochronological, and structural data that document neotectonic activity along hitherto unrecognized faults along the flanks of the NPI. New apatite (U-Th)/He bedrock cooling ages suggest faulting since 2-3 Ma. We infer the northward translation of an similar to 140 km long fore-arc sliver-the NPI block-results from enhanced partitioning of oblique plate convergence due to the closely spaced collision of three successive segments of the Chile Rise. In this model, greater uplift occurs in the hanging wall of the Exploradores thrust at the northern leading edge of the NPI block, whereas the Cachet and Liquine-Ofqui dextral faults decouple the NPI block along its eastern and western flanks, respectively. Localized extension possibly occurs at its southern trailing edge along normal faults associated with margin-parallel extension, tectonic subsidence, and lower elevations along the Andean crest line. Our neotectonic model provides a novel explanation for the abrupt topographic variations inland of the Chile Triple Junction and emphasizes the fundamental effects of local tectonics on exhumation and topographic patterns in this glaciated landscape.