Feedback between megathrust earthquake cycle and plate convergence.

Over million years, convergence between the Nazca and South America tectonic plates results in Andean orogeny. Over decades/centuries, it fuels the earthquake cycle of the Andean megathrust. It is well recognised that, over the geologically-long term of million years, Andean orogeny feeds back onto plate convergence rates, generating temporal changes documented throughout the Neogene. In contrast, no feedback mechanism operated over the geologically-short term by the earthquake cycle is currently contemplated. In fact, it is commonly assumed that the rates of contemporary convergence, which are accurately measured via geodesy, remain steady during the megathrust earthquake cycle. Here we investigate whether the contemporary Nazca/South America plate motion varies over year-/decade-long periods in response to megathrust stress variations associated with the earthquake cycle. We focus on the decade preceding the three largest and most recent [Formula: see text] earthquakes (2010 [Formula: see text] Maule, 2014 [Formula: see text] Iquique, 2015 [Formula: see text] Illapel), and find slowdowns of both Nazca and South America whole-plate motions that exceed the impact of data uncertainty or noise. We show that the torque variations required upon Nazca and South America to generate the slowdowns are consistent with that arising from the buildup of interseismic stress preceding the earthquakes.

A comprehensive database of active and potentially-active continental faults in Chile at 1:25,000 scale.

In seismically-active regions, mapping active and potentially-active faults is the first step to assess seismic hazards and site selection for paleoseismic studies that will estimate recurrence rates. Here, we present a comprehensive database of active and potentially-active continental faults in Chile based on existing studies and new mapping at 1:25,000 scale using geologic and geomorphic criteria and digital elevation models derived from TanDEM-X and LiDAR data. The database includes 958 fault strands grouped into 17 fault systems and classified based on activity (81 proved, 589 probable, 288 possible). The database is a contribution to the world compilation of active faults with applications among others in seismic hazard assessments, territorial planning, paleoseismology, geodynamics, landscape evolution processes, geothermal exploration, and in the study of feedbacks between continental deformation and the plate-boundary seismic cycle along subduction zones.

Fast Holocene slip and localized strain along the Liquiñe-Ofqui strike-slip fault system, Chile.

In active tectonic settings dominated by strike-slip kinematics, slip partitioning across subparallel faults is a common feature; therefore, assessing the degree of partitioning and strain localization is paramount for seismic hazard assessments. Here, we estimate a slip rate of 18.8 ± 2.0 mm/year over the past 9.0 ± 0.1 ka for a single strand of the Liquiñe-Ofqui Fault System, which straddles the Main Cordillera in Southern Chile. This Holocene rate accounts for ~ 82% of the trench-parallel component of oblique plate convergence and is similar to million-year estimates integrated over the entire fault system. Our results imply that strain localizes on a single fault at millennial time scale but over longer time scales strain localization is not sustained. The fast millennial slip rate in the absence of historical Mw > 6.5 earthquakes along the Liquiñe-Ofqui Fault System implies either a component of aseismic slip or Mw ~ 7 earthquakes involving multi-trace ruptures and > 150-year repeat times. Our results have implications for the understanding of strike-slip fault system dynamics within volcanic arcs and seismic hazard assessments.

Lithospheric flexure and rheology determined by climate cycle markers in the Corinth Rift.

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 ~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 ~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° 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.

Resilience of an aquatic macrophyte to an anthropogenically induced environmental stressor in a Ramsar wetland of southern Chile.

In mid-2004, anthropogenically induced changes in water quality of the Río Cruces wetland, a Ramsar site located in southern Chile (ca. 40°S), enhanced the resuspension of iron-enriched sediments, which were subsequently deposited over the most abundant aquatic macrophyte of the wetland (Egeria densa Planch. 1849). This event triggered the formation of brownish, necrotic patches and increased iron contents in the leaves and stems of E. densa, which contributed to a significant demise of the plant within the wetland. In this study, we estimate the recovery time as a proxy for resilience of this macrophyte at organismal and population levels. Macro- and micro-optical characteristics, as well as iron contents in tissues of E. densa, were documented in four time windows (2004, 2008, 2012, and 2014). In addition, the size of the macrophyte population and its spatial occurrence were monitored from 2008 to 2016 across 36 study sites within the wetland. Our results suggest necrotic patches and high iron contents in E. densa persisted at least until 2008. After 2013, a significant increase in the spatial occurrence of E. densa was observed within the wetland, reaching full recovery of the population during 2015. The health of plant tissues and iron contents in leaves and stems showed recovery period close to 4 years, while the recovery of the spatial occurrence of E. densa took approximately 9 years. While the monitoring of plant health was not performed on a strict annual basis, the recovery rates estimated here are slower than those described for other macrophytes. This finding might reflect the long-lasting effects of the disturbance from 2004 and the interaction with biotic processes, such as foraging by waterbirds recolonizing the Río Cruces wetland. These results show that full recovery of E. densa was achieved through a cascade of effects starting with abiotic factors (water quality) and passing through physiological and individual levels, to finally reach the population level. A key aspect of this response is the invasive nature of the macrophyte, which likely contributed to its recovery as a consequence of improved water quality. Less successful macrophyte species in other systems may not reach the specific population recovery, and become subdominant species instead, or even be eradicated from the wetland either as the result of herbivory or due to competition with other macrophytes.

Recovery of black-necked swans, macrophytes and water quality in a Ramsar wetland of southern Chile: Assessing resilience following sudden anthropogenic disturbances.

In 2004 migration and mortality for unknown reasons of the herbivorous Black necked swan (Cygnus melancorhyphus (Molina, 1782)) occurred within the Río Cruces wetland (southern Chile), a Ramsar Site and nature sanctuary. Before 2004, this wetland hosted the largest breeding population of this water bird in the Neotropic Realm. The concurrent decrease in the spatial occurrence of the aquatic plant Egeria densa Planch. 1849 - the main food source of swans - was proposed as a cause for swan migration and mortality. Additionally, post-mortem analyses carried out on swans during 2004 showed diminished body weight, high iron loads and histopathological abnormalities in their livers, suggesting iron storage disease. Various hypotheses were postulated to describe those changes; the most plausible related to variations in water quality after a pulp mill located upstream the wetland started to operate in February 2004. Those changes cascaded throughout the stands of E. densa whose remnants had high iron contents in their tissues. Here we present results of a long-term monitoring program of the wetland components, which show that swan population abundance, body weights and histological liver conditions recovered to pre-disturbance levels in 2012. The recovery of E. densa and iron content in plants throughout the wetland, also returned to pre-disturbance levels in the same 8-year time period. These results show the temporal scale over which resilience and natural restoring processes occur in wetland ecosystems of temperate regions such as southern Chile.

Back to full interseismic plate locking decades after the giant 1960 Chile earthquake.

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 (Mw = 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 ~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.

Calibrating coseismic coastal land-level changes during the 2014 Iquique (Mw=8.2) earthquake (northern Chile) with leveling, GPS and intertidal biota.

The April 1st 2014 Iquique earthquake (MW 8.1) occurred along the northern Chile margin where the Nazca plate is subducted below the South American continent. The last great megathrust earthquake here, in 1877 of Mw ~8.8 opened a seismic gap, which was only partly closed by the 2014 earthquake. Prior to the earthquake in 2013, and shortly after it we compared data from leveled benchmarks, deployed campaign GPS instruments, continuous GPS stations and estimated sea levels using the upper vertical level of rocky shore benthic organisms including algae, barnacles, and mussels. Land-level changes estimated from mean elevations of benchmarks indicate subsidence along a ~100-km stretch of coast, ranging from 3 to 9 cm at Corazones (18°30'S) to between 30 and 50 cm at Pisagua (19°30'S). About 15 cm of uplift was measured along the southern part of the rupture at Chanabaya (20°50'S). Land-level changes obtained from benchmarks and campaign GPS were similar at most sites (mean difference 3.7±3.2 cm). Higher differences however, were found between benchmarks and continuous GPS (mean difference 8.5±3.6 cm), possibly because sites were not collocated and separated by several kilometers. Subsidence estimated from the upper limits of intertidal fauna at Pisagua ranged between 40 to 60 cm, in general agreement with benchmarks and GPS. At Chanavaya, the magnitude and sense of displacement of the upper marine limit was variable across species, possibly due to species-dependent differences in ecology. Among the studied species, measurements on lithothamnioid calcareous algae most closely matched those made with benchmarks and GPS. When properly calibrated, rocky shore benthic species may be used to accurately measure land-level changes along coasts affected by subduction earthquakes. Our calibration of those methods will improve their accuracy when applied to coasts lacking pre-earthquake data and in estimating deformation during pre-instrumental earthquakes.

Responses of dune plant communities to continental uplift from a major earthquake: sudden releases from coastal squeeze.

Vegetated dunes are recognized as important natural barriers that shelter inland ecosystems and coastlines suffering daily erosive impacts of the sea and extreme events, such as tsunamis. However, societal responses to erosion and shoreline retreat often result in man-made coastal defence structures that cover part of the intertidal and upper shore zones causing coastal squeeze and habitat loss, especially for upper shore biota, such as dune plants. Coseismic uplift of up to 2.0 m on the Peninsula de Arauco (South central Chile, ca. 37.5º S) caused by the 2010 Maule earthquake drastically modified the coastal landscape, including major increases in the width of uplifted beaches and the immediate conversion of mid to low sandy intertidal habitat to supralittoral sandy habitat above the reach of average tides and waves. To investigate the early stage responses in species richness, cover and across-shore distribution of the hitherto absent dune plants, we surveyed two formerly intertidal armoured sites and a nearby intertidal unarmoured site on a sandy beach located on the uplifted coast of Llico (Peninsula de Arauco) over two years. Almost 2 years after the 2010 earthquake, dune plants began to recruit, then rapidly grew and produced dune hummocks in the new upper beach habitats created by uplift at the three sites. Initial vegetation responses were very similar among sites. However, over the course of the study, the emerging vegetated dunes of the armoured sites suffered a slowdown in the development of the spatial distribution process, and remained impoverished in species richness and cover compared to the unarmoured site. Our results suggest that when released from the effects of coastal squeeze, vegetated dunes can recover without restoration actions. However, subsequent human activities and management of newly created beach and dune habitats can significantly alter the trajectory of vegetated dune development. Management that integrates the effects of natural and human induced disturbances, and promotes the development of dune vegetation as natural barriers can provide societal and conservation benefits in coastal ecosystems.

Ecological implications of extreme events: footprints of the 2010 earthquake along the Chilean coast.

Deciphering ecological effects of major catastrophic events such as earthquakes, tsunamis, volcanic eruptions, storms and fires, requires rapid interdisciplinary efforts often hampered by a lack of pre-event data. Using results of intertidal surveys conducted shortly before and immediately after Chile's 2010 M(w) 8.8 earthquake along the entire rupture zone (ca. 34-38°S), we provide the first quantification of earthquake and tsunami effects on sandy beach ecosystems. Our study incorporated anthropogenic coastal development as a key design factor. Ecological responses of beach ecosystems were strongly affected by the magnitude of land-level change. Subsidence along the northern rupture segment combined with tsunami-associated disturbance and drowned beaches. In contrast, along the co-seismically uplifted southern rupture, beaches widened and flattened increasing habitat availability. Post-event changes in abundance and distribution of mobile intertidal invertebrates were not uniform, varying with land-level change, tsunami height and coastal development. On beaches where subsidence occurred, intertidal zones and their associated species disappeared. On some beaches, uplift of rocky sub-tidal substrate eliminated low intertidal sand beach habitat for ecologically important species. On others, unexpected interactions of uplift with man-made coastal armouring included restoration of upper and mid-intertidal habitat seaward of armouring followed by rapid colonization of mobile crustaceans typical of these zones formerly excluded by constraints imposed by the armouring structures. Responses of coastal ecosystems to major earthquakes appear to vary strongly with land-level change, the mobility of the biota and shore type. Our results show that interactions of extreme events with human-altered shorelines can produce surprising ecological outcomes, and suggest these complex responses to landscape alteration can leave lasting footprints in coastal ecosystems.

Land-level changes produced by the Mw 8.8 2010 Chilean earthquake.

We observed vertically displaced coastal and river markers after the 27 February 2010 Chilean earthquake [moment magnitude (Mw) 8.8]. Land-level changes range between 2.5 and -1 meters, evident along an approximately 500-kilometers-long segment identified here as the maximum length of coseismic rupture. A hinge line located 120 kilometers from the trench separates uplifted areas, to the west, from subsided regions. A simple elastic dislocation model fits these observations well; model parameters give a similar seismic moment to seismological estimates and suggest that most of the plate convergence since the 1835 great earthquake was elastically stored and then released during this event.