Refine
Year of publication
Language
- English (52)
Is part of the Bibliography
- yes (52)
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
Individual great earthquakes are posited to release the elastic strain energy that has accumulated over centuries by the gradual movement of tectonic plates(1,2). However, knowledge of plate deformation during a complete seismic cycle-two successive great earthquakes and the intervening interseismic period-remains incomplete(3). A complete seismic cycle began in south-central Chile in 1835 with an earthquake of about magnitude 8.5 (refs 4,5) and ended in 2010 with a magnitude 8.8 earthquake(6). During the first earthquake, an uplift of Isla Santa Maria by 2.4 to 3m was documented(4,5). In the second earthquake, the island was uplifted(7) by 1.8 m. Here we use nautical surveys made in 1804, after the earthquake in 1835 and in 1886, together with modern echo sounder surveys and GPS measurements made immediately before and after the 2010 earthquake, to quantify vertical deformation through the complete seismic cycle. We find that in the period between the two earthquakes, Isla Santa Maria subsided by about 1.4 m. We simulate the patterns of vertical deformation with a finite-element model and find that they agree broadly with predictions from elastic rebound theory(2). However, comparison with geomorphic and geologic records of millennial coastline emergence(8,9) reveal that 10-20% of the vertical uplift could be permanent.
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.
Major earthquakes ( M > 8) have repeatedly ruptured the Nazca-South America plate interface of south-central Chile involving meter scale land-level changes. Earthquake recurrence intervals, however, extending beyond limited historical records are virtually unknown, but would provide crucial data on the tectonic behavior of forearcs. We analyzed the spatiotemporal pattern of Holocene earthquakes on Santa Maria Island (SMI; 37 degrees S), located 20 km off the Chilean coast and approximately 70 km east of the trench. SMI hosts a minimum of 21 uplifted beach berms, of which a subset were dated to calculate a mean uplift rate of 2.3 +/- 0.2 m/ky and a tilting rate of 0.022 +/- 0.002 degrees/ky. The inferred recurrence interval of strandline-forming earthquakes is similar to 180 years. Combining coseismic uplift and aseismic subsidence during an earthquake cycle, the net gain in strandline elevation in this environment is similar to 0.4 m per event
The architecture of coastal sequences in tectonically-active regions results mostly from a combination of sea-level and land-level changes. The objective of this study is to unravel these signals by combining sequence stratigraphy and sedimentology of near-shore sedimentary sequences in wave-built terraces. We focus on Santa Maria Island at the south-central Chile margin, which hosts excellent exposures of coastal sediments from Marine Isotope Stage 3. A novel method based on statistical analysis of grain-size distributions coupled with fades descriptions provided a detailed account of transgressive-regressive cycles. Radiocarbon ages from paleosols constrain the chronology between >53 and similar to 31 cal ka BP. Because the influence of glaciations can be neglected, we calculated relative sea-level curves by tying the onset of deposition on a bedrock abrasion platform to a global sea-level curve. The observed depositional cycles match those predicted for uplift rates between 1.2 and 1.8 m/ka. The studied sedimentary units represent depositional cycles that resulted in reoccupation events of an existing marine terrace. Our study demonstrates wave-built marine terrace deposits along clastic shorelines in temperate regions can be used to distinguish between tectonic uplift and climate-induced sea-level changes.
Much progress has been made in estimating recurrence intervals of great and giant subduction earthquakes using terrestrial, lacustrine, and marine paleoseismic archives. Recent detailed records suggest these earthquakes may have variable recurrence periods and magnitudes forming supercycles. Understanding seismic supercycles requires long paleoseismic archives that record timing and magnitude of such events. Turbidite paleoseismic archives may potentially extend past earthquake records to the Pleistocene and can thus complement commonly shorter-term terrestrial archives. However, in order to unambiguously establish recurring seismicity as a trigger mechanism for turbidity currents, synchronous deposition of turbidites in widely spaced, isolated depocenters has to be ascertained. Furthermore, characteristics that predispose a seismically active continental margin to turbidite paleoseismology and the correct sample site selection have to be taken into account.
Here we analyze 8 marine sediment cores along 950 km of the Chile margin to test for the feasibility of compiling detailed and continuous paleoseismic records based on turbidites. Our results suggest that the deposition of areally widespread, synchronous turbidites triggered by seismicity is largely controlled by sediment supply and, hence, the climatic and geomorphic conditions of the adjacent subaerial setting. The feasibility of compiling a turbidite paleoseismic record depends on the delicate balance between sufficient sediment supply providing material to fail frequently during seismic shaking and sufficiently low sedimentation rates to allow for coeval accumulation of planktonic foraminifera for high-resolution radiocarbon dating.
We conclude that offshore northern central Chile (29-32.5 degrees S) Holocene turbidite paleoseismology is not feasible, because sediment supply from the semi-arid mainland is low and almost no Holocene turbidity-current deposits are found in the cores. In contrast, in the humid region between 36 and 38 degrees S frequent Holocene turbidite deposition may generally correspond to paleoseismic events. However, high terrigenous sedimentation rates prevent high-resolution radiocarbon dating. The climatic transition region between 32.5 and 36 degrees S appears to be best suited for turbidite paleoseismology. (C) 2015 Elsevier Ltd. All rights reserved.
The Maule earthquake of 27th February 2010 (M-w = 8.8) affected similar to 500 km of the Nazca-South America plate boundary in south-central Chile producing spectacular crustal deformation. Here, we present a detailed estimate of static coseismic surface offsets as measured by survey and continuous GPS, both in near- and far-field regions. Earthquake slip along the megathrust has been inferred from a Joint inversion of our new data together with published GPS, InSAR, and land-level changes data using Green's functions generated by a spherical finite-element model with realistic subduction zone geometry. The combination of the data sets provided a good resolution, indicating that most of the slip was well resolved. Coseismic slip was concentrated north of the epicenter with up to 16 m of slip, whereas to the south it reached over 10 m within two minor patches. A comparison of coseismic slip with the slip deficit accumulated since the last great earthquake in 1835 suggests that the 2010 event closed a mature seismic gap. Slip deficit distribution shows an apparent local overshoot that highlight cycle-to-cycle variability, which has to be taken into account when anticipating future events from interseismic observations. Rupture propagation was obviously not affected by bathymetric features of the incoming plate. Instead, splay faults in the upper plate seem to have limited rupture propagation in the updip and along-strike directions. Additionally, we found that along-strike gradients in slip are spatially correlated with geometrical inflections of the megathrust. Our study suggests that persistent tectonic features may control strain accumulation and release along subduction megathrusts.
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 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.
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.