Refine
Year of publication
Language
- English (52)
Is part of the Bibliography
- yes (52) (remove)
Keywords
- Chile (5)
- Marine terraces (3)
- Maule earthquake (3)
- TerraceM (3)
- African Humid Period (2)
- Central Chile (2)
- Coastal uplift (2)
- GPS (2)
- Geomorphology (2)
- LiDAR (2)
Institute
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.
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.
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.
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.
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.
Lago Laja is a late Quaternary volcanic‐dammed lake located near the drainage divide of the south central Andes. Field observations, lake reflection seismic profiles, bathymetry, and remote sensing data reveal an active fault system that runs parallel to the volcanic arc along the axis of the Main Cordillera, the Lago Laja fault system (LLFS). Normal faults of this extensional system cut late Pleistocene volcanics, <7.1 ka still water lacustrine sediments, 6.3 ka pyroclastic deposits, and Holocene alluvial fans. We divide the LLFS in three segments on the basis of fault geometry, width, and slip magnitude. The underwater faults of the central segment in the lake's deepest part have the maximum Holocene vertical slip rate of >2.7 mm/yr. Since 7.1 ka, the LLFS accounts for ∼0.7% of arc‐normal extension at an average minimum rate of 1.2 mm/yr and strain rate of ∼10−14 s−1. Seismites and surface ruptures evidence M>6 paleoearthquakes. The Main Cordillera at ∼37°S is a large‐scale pop‐up structure uplifted by thrusting along its foothills. In this light, we interpret extension in the axial and highest part of the Andes as incipient synorogenic gravitational collapse in response to uplift and crustal thickening. Thermal weakening due to elevated heat flow and postglacial lithospheric rebound and unbending have probably contributed to the arc‐limited collapse and Holocene acceleration of deformation rates. The lack of significant strike‐slip offsets along the LLFS as well as along both foothills‐thrust systems at 37°S contrasts with the intra‐arc dextral fault zone south of 38°S. Regional structural data indicates that north of 38°S, diffusely distributed strain reflects low partitioning of oblique subduction, while to the south deformation is localized in a discrete strike‐slip fault zone along the volcanic arc, reflecting a higher degree of partitioning. We relate this strain partitioning gradient to favorable fault orientations in the fore arc north of the Arauco Peninsula, a major seismotectonic boundary.
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.
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.
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.