Search

We deployed a dense geodetic and seismological network in the Atacama seismic gap in Chile. We derive a microseismicity catalog of >30,000 events, time series from 70 GNSS stations, and utilize a transdimensional Bayesian inversion to estimate interplate locking. We identify two highly locked regions of different sizes whose geometries appear to control seismicity patterns. Interface seismicity concentrates beneath the coastline, just downdip of the highest locking. A region with lower locking (27.5°S–27.7°S) coincides with higher seismicity levels, a high number of repeating earthquakes and events extending toward the trench. This area is situated where the Copiapó Ridge is subducted and has shown previous indications of both seismic and aseismic slip, including an earthquake sequence in 2020. While these findings suggest that the structure of the downgoing oceanic plate prescribes patterns of interplate locking and seismicity, we note that the Taltal Ridge further north lacks a similar signature. Plain Language Summary: Deformation along plate boundaries can occur seismically (i.e. through earthquakes) as well as aseismically (i.e. slipping slowly), and it is important to understand where each of these modes is dominant. Along the Chilean subduction contact, North‐Central Chile is the only place where aseismic deformation episodes have been observed so far. In order to study these processes in detail, we deployed and operated dense geodetic and seismological networks in this region. Analyzing the data collected by these networks, we find notable relationships between seismic and aseismic processes. Thousands of small earthquakes are found at the boundaries of locked regions, whereas no small earthquakes are found at their interior. Thus, implying such regions are mechanically coupled, that is, currently accumulating elastic deformation energy that will 1 day be released during a large earthquake. Along the North‐Central Chilean plate boundary, there is one region (around 27.5°S) that shows many signs of aseismic deformation. It is located where a chain of seamounts is being subducted, which is likely responsible for the different behavior of this segment. Key Points: Microseismicity catalog and map of interplate locking derived for the Atacama 1922 seismic gap in North‐Central ChileSeismicity in vicinity of plate interface coincides with downdip edge of high couplingSeismo‐geodetic signals due to the subduction of the Copiapó ridge are prominent but negligible for the subducting Taltal Ridge [ABSTRACT FROM AUTHOR]

At the Blanco transform fault system (BTFS) off Oregon, 138 local earthquakes and 84 double‐couple focal mechanisms from ocean‐bottom‐seismometer recordings jointly discussed with bathymetric features reveal a highly segmented transform system without any prominent fracture zone traces longer than 100 km. In the west, seismicity is focused at deep troughs (i.e., the West and East Blanco, and Surveyor Depressions). In the east, the BTFS lacks a characteristic transform valley and instead developed the Blanco Ridge, which is the most seismically active feature, showing strike‐slip and dip‐slip faulting. Sandwiched between the two main segments of the BTFS is the Cascadia Depression, representing a short intra‐transform spreading segment. Seismic slip vectors reveal that stresses at the eastern BTFS are roughly in line with plate motion. In contrast, stresses to the west are clockwise skewed, indicating ongoing reorganization of the OTF system. As we observed no prominent fracture zones at the BTFS, plate tectonic reconstructions suggest that the BTFS developed from non‐transform offsets rather than pre‐existing transform faults during a series of ridge propagation events. Our observations suggest that the BTFS can be divided into two oceanic transform systems. The eastern BTFS is suggested to be a mature transform plate boundary since ∼0.6 Ma. In contrast, the western BTFS is an immature transform system, which is still evolving to accommodate far‐field stress change. The BTFS acts as a natural laboratory to yield processes governing the development of oceanic transform faults. Plain Language Summary: The Blanco transform fault system (BTFS) northwest off the coast of Oregon is seismically very active. We used 1 year of ocean bottom seismometer data collected between September 2012 and October 2013 to locate 138 local earthquakes. The events align perfectly with the morphologic features of the BTFS, dividing the BTFS into five transform segments and two short intra‐transform spreading centers. Furthermore, we observe different seismotectonic behaviors of the western and eastern BTFS based on the along‐strike variation in morphology, magnetization, focal depth distribution, and strain partitioning. Although many segmented oceanic transform systems were formed from a single transform fault in response to rotations in plate motion, the BTFS turns out to be originated from non‐transform offsets between ridge segments, as we observed no prominent fracture zone traces neither in morphology nor gravity field data. A clockwise shift in the Juan de Fuca/Pacific pole of rotation at ∼5 Ma followed by a series of ridge propagation events initiated the formation of the BTFS, integrated each segment of the BTFS by shortening the ridge segments in between. Our observations suggest that the Blanco Ridge and the Gorda transform segment in the eastern BTFS were formed at ∼1.6 and 0.6 Ma, respectively, and ever since, the eastern BTFS became a mature transform boundary. In contrast, seismic slip vectors comparing to plate motion directions reveal that stresses in the western BTFS are systematically skewed, suggesting the immature transform plate boundary is still adjusting to the new stress regime. Key Points: Local seismicity of the Blanco transform fault system (BTFS) reveals along‐strike variations dominated by strike‐slip and oblique dip‐slipThe BTFS developed from non‐transform offsets rather than discrete transform faults in response to plate rotation and ridge propagationThe BTFS consists of a mature plate boundary in the east and an immature system in the west, separated by a central spreading center [ABSTRACT FROM AUTHOR]

SUMMARY: Volcano seismology is an essential tool for monitoring volcanic processes in the advent and during eruptions. A variety of seismic signals can be recorded at volcanoes, of which some are thought to be related to the migration of fluids which is of primary importance for the anticipation of imminent eruptions. We investigate the volcanic crises at Villarrica volcano in 2015 and report on a newly discovered very-long-period (VLP) signal that accompanies phases of periodic long period (LP) signal burst. Despite their low amplitude emergent character, we can locate the source region of the 1 Hz LP signals to the close vicinity of the volcano using a network-based correlation method. The source of the VLP signal with a period of about 30–100 s appears to locate in the vicinity of two stations a few kilometres from the summit. Both stations record very similar VLP waveforms that are correlated with the envelope of the LP bursts. A shallow magma reservoir was inferred by Contreras from surface deformation as the source of inflation following the eruption in 2015. Cyclic volume changes of 6 m3 in this reservoir at 3 km depth can explain the observed amplitudes of the vertical VLP signal. We propose that the LP signal is generated by the migration of gas or gas-rich magma that is periodically released from the inflating reservoir through a non-linear valve structure which modulates the flux, and thereby causes bursts of flow-related LP signals and pressure changes observed as VLP deformation. Our model predicts that the correlated occurrence of LP bursts and VLP surface motion depends on the intensity of the fluid flux. A weaker flux of fluids may not exceed the opening pressure of valve structure, and higher rates might maintain pressure above the closing pressure. In both cases, the VLP signal vanishes. Our observation provides constrains for models of fluid transport inside volcanoes. At Villarrica the VLP signal, and its relation to the LP activity, reveal additional information about fluxes in the magmatic reservoir that might aide forecasting of volcanic activity. [ABSTRACT FROM AUTHOR]

Seismic event detection and phase picking are the base of many seismological workflows. In recent years, several publications demonstrated that deep learning approaches significantly outperform classical approaches, achieving human‐like performance under certain circumstances. However, as studies differ in the datasets and evaluation tasks, it is unclear how the different approaches compare to each other. Furthermore, there are no systematic studies about model performance in cross‐domain scenarios, that is, when applied to data with different characteristics. Here, we address these questions by conducting a large‐scale benchmark. We compare six previously published deep learning models on eight data sets covering local to teleseismic distances and on three tasks: event detection, phase identification and onset time picking. Furthermore, we compare the results to a classical Baer‐Kradolfer picker. Overall, we observe the best performance for EQTransformer, GPD and PhaseNet, with a small advantage for EQTransformer on teleseismic data. Furthermore, we conduct a cross‐domain study, analyzing model performance on data sets they were not trained on. We show that trained models can be transferred between regions with only mild performance degradation, but models trained on regional data do not transfer well to teleseismic data. As deep learning for detection and picking is a rapidly evolving field, we ensured extensibility of our benchmark by building our code on standardized frameworks and making it openly accessible. This allows model developers to easily evaluate new models or performance on new data sets. Furthermore, we make all trained models available through the SeisBench framework, giving end‐users an easy way to apply these models. Plain Language Summary: The first step in many seismological workflows is identifying if a signal contains an earthquake, and at which time which type of seismic wave arrived. These steps are known as event detection, phase identification and phase picking. In recent years, machine learning methods, in particular deep learning methods have been developed, showing promising performance on these tasks. However, so far these models have not been compared systematically in a quantitative way. Here we evaluate the performance of six deep learning models on eight datasets. Additionally, we compare them to a traditional picking algorithm not using machine learning. From our results we identify that the models EQTransformer, GPD and PhaseNet perform best. As in many use cases no picker trained on the target region will be available, we further evaluated how well models are transferable across regions. We identified that transfer across regions works well as long as the distance ranges stay similar. To foster application of the results, we make all our trained models available through the SeisBench framework. Key Points: We conducted a large scale benchmark of machine learning pickers using six models and eight datasetsBest overall performance is observed for EQTransformer, GPD and PhaseNet, with advantages for EQTransformer on teleseismic distancesModels transfer well between different regions with similar distances, but not between regional and teleseismic distances [ABSTRACT FROM AUTHOR]

The Liguro-Provençal basin was formed as a back-arc basin of the retreating Calabrian–Apennines subduction zone during the Oligocene and Miocene. The resulting rotation of the Corsica–Sardinia block is associated with rifting, shaping the Ligurian Basin. It is still debated whether oceanic or atypical oceanic crust was formed or if the crust is continental and experienced extreme thinning during the opening of the basin. We perform ambient noise tomography, also taking into account teleseismic events, using an amphibious network of seismic stations, including 22 broadband ocean bottom seismometers (OBSs), to investigate the lithospheric structure of the Ligurian Basin. The instruments were installed in the Ligurian Basin for 8 months between June 2017 and February 2018 as part of the AlpArray seismic network. Because of additional noise sources in the ocean, OBS data are rarely used for ambient noise studies. However, we carefully pre-process the data, including corrections for instrument tilt and seafloor compliance and excluding higher modes of the ambient-noise Rayleigh waves. We calculate daily cross-correlation functions for the AlpArray OBS array and surrounding land stations. We also correlate short time windows that include teleseismic earthquakes, allowing us to derive surface wave group velocities for longer periods than using ambient noise only. We obtain group velocity maps by inverting Green's functions derived from the cross-correlation of ambient noise and teleseismic events, respectively. We then used the resulting 3D group velocity information to calculate 1D depth inversions for S-wave velocities. The group velocity and shear-wave velocity results compare well to existing large-scale studies that partly include the study area. In onshore France, we observe a high-velocity area beneath the Argentera Massif, roughly 10 km below sea level. We interpret this as the root of the Argentera Massif. Our results add spatial resolution to known seismic velocities in the Ligurian Basin, thereby augmenting existing seismic profiles. In agreement with existing seismic studies, our shear-wave velocity maps indicate a deepening of the Moho from 12 km at the south-western basin centre to 20–25 km at the Ligurian coast in the north-east and over 30 km at the Provençal coast. The maps also indicate that the south-western and north-eastern Ligurian Basin are structurally separate. The lack of high crustal vP/vS ratios beneath the south-western part of the Ligurian Basin preclude mantle serpentinisation there. [ABSTRACT FROM AUTHOR]

The northern margin of the Ligurian Basin shows notable seismicity at the Alpine front, including frequent magnitude 4 events. Seismicity decreases offshore towards the Basin centre and Corsica, revealing a diffuse distribution of low-magnitude earthquakes. We analyse data of the amphibious AlpArray seismic network with focus on the offshore component, the AlpArray ocean bottom seismometer (OBS) network, consisting of 24 broadband OBSs deployed for 8 months, to reveal the seismicity and depth distribution of micro-earthquakes beneath the Ligurian Sea. Two clusters occurred between ∼ 10 km to ∼ 16 km depth below the sea surface, within the lower crust and uppermost mantle. Thrust faulting focal mechanisms indicate compression and an inversion of the Ligurian Basin, which is an abandoned Oligocene–Miocene rift basin. The basin inversion is suggested to be related to the Africa–Europe plate convergence. The locations and focal mechanisms of seismicity suggest reactivation of pre-existing rift-related structures. Slightly different striking directions of presumed rift-related faults in the basin centre compared to faults further east and hence away from the rift basin may reflect the counter-clockwise rotation of the Corsica–Sardinia block. High mantle S-wave velocities and a low Vp/Vs ratio support the hypothesis of strengthening of crust and uppermost mantle during the Oligocene–Miocene rifting-related extension and thinning of continental crust. [ABSTRACT FROM AUTHOR]

We resolve a previously unrecognized shallow subducting seamount from a re-processed multichannel seismic profile crossing the 1994 Mw 7.8 Java tsunami earthquake rupture area. Seamount subduction occurs where the overriding plate experiences uplift by lateral shortening and vertical thickening. Pronounced back-thrusting at the landward slope of the forearc high and the formation of splay faults branching off the landward flank of the subducting seamount are observed. The location of the seamount in relation to the 1994 earthquake hypocentre and its co-seismic slip model suggests that the seamount acted as a seismic barrier to the up-dip co-seismic rupture propagation of this moderate-size earthquake. [ABSTRACT FROM AUTHOR]

The aftershock distribution of the 2014 Mw 8.1 Iquique earthquake offshore northern Chile, identified from a long‐term deployment of ocean bottom seismometers installed eight months after the mainshock, in conjunction with seismic reflection imaging, provides insights into the processes regulating the updip limit of coseismic rupture propagation. Aftershocks updip of the mainshock hypocenter frequently occur in the upper plate and are associated with normal faults identified from seismic reflection data. We propose that aftershock seismicity near the plate boundary documents subduction erosion that removes mass from the base of the wedge and results in normal faulting in the upper plate. The combination of very little or no sediment accretion and subduction erosion over millions of years has resulted in a very weak and aseismic frontal wedge. Our observations thus link the shallow subduction zone seismicity to subduction erosion processes that control the evolution of the overriding plate. Plain Language Summary: To better understand the controls on shallow seismicity and subduction erosion following large subduction earthquakes, we use marine recordings of the Mw 8.1 2014 Iquique earthquake aftershocks and long‐offset multi‐channel seismic data. By comparing the aftershock locations and seismic imaging, we observe that most aftershocks occurred in the upper continental plate and abruptly stopped in the frontal forearc. The amplitude characteristics of upper‐crust reflections indicate a fractured and fluid‐filled outer forearc, which is associated with the absence of aftershocks. Large‐scale faulting, as evidenced by disrupted reflections in the seismic image, can be correlated to upper plate seismicity. We propose that the aftershocks updip of the main earthquake area reflect active subduction erosion processes. Key Points: We investigate structure and seismicity at the updip end of the 2014 Iquique earthquake rupture using amphibious seismic dataSeismicity updip of the 2014 Iquique earthquake occurs over a broad range likely interpreted to be related to the basal erosion processesCoseismic stress changes and aftershocks activate extensional faulting of the upper plate and subduction erosion [ABSTRACT FROM AUTHOR]

The Liguro-Provençal basin was formed as a back-arc basin of the retreating Calabrian-Apennines subduction zone during the Oligocene and Miocene. The resulting rotation of the Corsica-Sardinia block is associated with rifting, shaping the Ligurian Sea. It is still debated whether oceanic or atypical oceanic crust was formed or if the crust is continental and experienced extreme thinning during the opening of the basin. We invert velocity models using an amphibious network of seismic stations, including 22 broadband Ocean Bottom Seismometers (OBS) to investigate the lithospheric structure of the Ligurian sea. The instruments were installed in the Ligurian Sea for eight months between June 2017 and February 2018 as part of the AlpArray seismic network. Because of additional noise sources in the ocean, OBS data are rarely used for ambient noise studies. However, we attentively pre-process the data, including corrections for instrument tilt and seafloor compliance. We took extra care to exclude higher modes of the ambient-noise Rayleigh waves. We calculate daily cross-correlation functions for the LOBSTER array and surrounding land stations. Additionally, we correlate short time windows that include teleseismic earthquakes that allow us to derive surface wave group velocities for longer periods than using ambient noise only. Group velocity maps are obtained by inverting Green's functions derived from the cross-correlation of ambient noise and teleseismic events, respectively. We then used the resulting 3D group velocity information to calculate 1D depth inversions for S-wave velocities. The shear-wave velocity results show a deepening of the Moho from 12 km at the southwestern basin centre to 20-25 km at the Ligurian coast in the northeast and over 30 km at the Provençal coast. We find no hint on mantle serpentinisation and no evidence for an Alpine slab, at least down to depths of 25 km. However, we see a separation of the southwestern and northeastern Ligurian Basin that coincides with the promoted prolongation of the Alpine front. [ABSTRACT FROM AUTHOR]

What controls subduction megathrust seismogenesis downdip of the mantle wedge corner (MWC)? We propose that, in the region of the 2010 Mw = 8.8 Maule, Chile, earthquake, serpentine minerals derived from the base of the hydrated mantle wedge exert a dominant control. Based on modeling, we predict that the megathrust fault zone near the MWC contains abundant lizardite/chrysotile‐rich serpentinite that transforms to antigorite‐rich serpentinite at greater depths. From the MWC at 32–40 km depth to at least 55 km, the predominantly velocity‐strengthening megathrust accommodated dynamic propagation of the 2010 rupture but with small slip and negative stress drop. The downdip distribution of interplate aftershocks exhibits a gap around the MWC that can be explained by the velocity‐strengthening behavior of lizardite/chrysotile. Interspersed velocity‐weakening and dynamic weakening antigorite‐rich patches farther downdip may be responsible for increased abundance of aftershocks and possibly for some of the high‐frequency energy radiation during the 2010 rupture. Plain Language Summary: A subduction megathrust rupture may extend from the trench to beneath populated coastal area. To understand what controls the megathrust seismogenic behavior of its deeper part, we conduct a case study of the 2010 magnitude 8.8 Maule, Chile, earthquake and its aftershocks. With numerical thermal modeling, we find that the deep seisomogenic behavior may be controlled by serpentine materials derived from the base of the overlying hydrated mantle wedge. The deep part of the megathrust ruptured "passively" against increasing resistance (negative stress drop). Lower‐temperature serpentines lizardite and chrysotile that are known to facilitate aseismic slip may be responsible for an observed gap in aftershock distribution in the downdip direction. Higher‐temperature serpentine antigorite that is known to facilitate seismic slip may be responsible for increased aftershocks farther downdip and possibly radiation of high‐frequency energy during the 2010 Maule earthquake. Key Points: Seismogenic behavior of subduction megathrust may be affected by materials derived from the base of the serpentinized forearc mantle wedgeThermopetrologic modeling suggests abundance of lizardite and chrysotile around the mantle wedge corner but antigorite farther downdipSerpentine frictional behavior may explain distributions of stress drop and seismic energy radiation in 2010 rupture and aftershocks [ABSTRACT FROM AUTHOR]

The Ligurian Basin is located in the Mediterranean Sea to the north-west of Corsica at the transition from the Western Alpine orogen to the Apennine system and was generated by the south-eastward trench retreat of the Apennines–Calabrian subduction zone. Late-Oligocene-to-Miocene rifting caused continental extension and subsidence, leading to the opening of the basin. Yet it remains unclear if rifting caused continental break-up and seafloor spreading. To reveal its lithospheric architecture, we acquired a 130 km long seismic refraction and wide-angle reflection profile in the Ligurian Basin. The seismic line was recorded in the framework of SPP2017 4D-MB, a Priority Programme of the German Research Foundation (DFG) and the German component of the European AlpArray initiative, and trends in a NE–SW direction at the centre of the Ligurian Basin, roughly parallel to the French coastline. The seismic data were recorded on the newly developed GEOLOG recorder, designed at GEOMAR, and are dominated by sedimentary refractions and show mantle Pn arrivals at offsets of up to 70 km and a very prominent wide-angle Mohorovičić discontinuity (Moho) reflection. The main features share several characteristics (e.g. offset range, continuity) generally associated with continental settings rather than documenting oceanic crust emplaced by seafloor spreading. Seismic tomography results are complemented by gravity data and yield a ∼ 6–8 km thick sedimentary cover and the seismic Moho at 11–13 km depth below the sea surface. Our study reveals that the oceanic domain does not extend as far north as previously assumed. Whether Oligocene–Miocene extension led to extremely thinned continental crust or exhumed subcontinental mantle remains unclear. A low grade of mantle serpentinisation indicates a high rate of syn-rift sedimentation. However, rifting failed before oceanic spreading was initiated, and continental crust thickens towards the NE within the northern Ligurian Basin. [ABSTRACT FROM AUTHOR]

The development of German proficiency by a three-year-old Australian boy living in Germany was monitored for a five-month period. His command of German negation is reported. The study is seen as bearing on issues in first and second language acquisition, such as competence and interference. (JB)

The Ligurian Basin is located in the Mediterranean Sea to the north-west of Corsica at the transition from the western Alpine orogen to the Apennine system and was generated by the south-eastward trench retreat of the Apennines-Calabrian subduction zone. Late Oligocene to Miocene rifting caused continental extension and subsidence, leading to the opening of the basin. Yet, it still remains enigmatic if rifting caused continental break-up and seafloor spreading. To reveal its lithospheric architecture, we acquired a state of the art seismic refraction and wide-angle reflection profile in the Ligurian Basin. The seismic line was recorded in the framework of SPP2017 4D-MB, the German component of the European AlpArray initiative, and trends in a NE-SW direction at the centre of the Ligurian Basin, roughly parallel to the French coastline. The seismic data recorded on the newly developed GEOLOG recorder, designed at GEOMAR, are dominated by sedimentary refractions and show mantle Pn arrivals at offsets of up to 70 km and a very prominent wide-angle Moho reflection. The main features share several characteristics (i.e. offset range, continuity) generally associated with continental settings rather than documenting oceanic crust emplaced by seafloor spreading. Seismic tomography results are augmented by gravity data and yield a 7.5-8 km thick sedimentary cover which is directly underlain by serpentinised mantle material at the south-western end of the profile. The acoustic basement at the north-eastern termination is interpreted to be continental crust, thickening towards the NE. Our study reveals that the oceanic domain does not extend as far north as previously assumed and that extension led to extreme continental thinning and exhumation of sub-continental mantle which eventually became serpentinised. [ABSTRACT FROM AUTHOR]

We examine the intra‐arc crustal seismicity of the Andean Southern Volcanic Zone. Our aim is to resolve interseismic deformation in an active magmatic arc dominated by both margin‐parallel (Liquiñe‐Ofqui fault system, LOFS) and Andean transverse faults. Crustal seismicity provides information about the schizosphere tectonic state, delineating the geometry and kinematics of high strain domains driven by oblique‐subduction. Here, we present local seismicity based on 16‐month data collected from 34 seismometers monitoring a ~200‐km‐long section of the Southern Volcanic Zone, including the Lonquimay and Villarrica volcanoes. We located 356 crustal events with magnitudes between Mw 0.6 and Mw 3.6. Local seismicity occurs at depths down to 40 km in the forearc and consistently shallower than 12 km beneath the volcanic chain, suggesting a convex shape of the crustal seismogenic layer bottom. Focal mechanisms indicate strike‐slip faulting consistent with ENE‐WSW shortening in line with the long‐term deformation history revealed by structural geology studies. However, we find regional to local‐scale variations in the shortening axes orientation as revealed by the nature and spatial distribution of microseismicity, within three distinctive latitudinal domains. In the northernmost domain, seismicity is consistent with splay faulting at the northern termination of the LOFS; in the central domain, seismicity distributes along ENE‐ and WNW‐striking discrete faults, spatially associated with, hitherto seismic Andean transverse faults. The southernmost domain, in turn, is characterized by activity focused along a N15°E striking master branch of the LOFS. These observations indicate a complex strain compartmentalization pattern within the intra‐arc crust, where variable strike‐slip faulting dominates over dip‐slip movements. Plain Language Summary: In active volcanic chains, there is a strong interplay between deformation and volcanism. In this research, we take the "pulse" of active tectonics in a volcanic arc setting by measuring natural seismicity. We installed a network of 34 seismometers over 200 km along the volcanic arc in Southcentral Chile. Our results indicate active faulting in coherence with long‐lived faults and the overall regional‐scale stress regime. In fact, crustal regions in which faults are more active are spatially associated with inherited faults and with higher temperature gradient domains, as inferred from volcanic and geothermal activity. The maximum depth of seismicity is shallow (<12 km) in the central part of the volcanic arc, whereas is much deeper (~40 km) toward the forearc and back‐arc regions. This observation suggests that elevated geothermal gradients lead to a thinner brittle portion of crust, which is then mechanically easier to (re‐)break up. Our results can be used for understanding the way by which faulting interacts with crustal fluids, which in turn may help improving geothermal exploration strategies and seismic hazard assessment. Furthermore, seismicity reveals detailed information on how ongoing deformation accommodates within complex fault systems characterized by different orientations and kinematics. Key Points: Crustal seismicity in the SVZ of the Andes reveals the nature of active upper crustal faulting consistent with long‐term arc tectonicsStrain partitioning is compartmentalized into arc‐parallel (LOFS) and Andean transverse faults (ATF)Crustal seismicity occurs at all depths down to 40 km depth in the forearc, but shallower than 12 km along the volcanic chain [ABSTRACT FROM AUTHOR]

The Sumatran subduction zone exhibits strong seismic and tsunamogenic potential with the prominent examples of the 2004, 2005 and 2007 earthquakes. Here, we invert travel-time data of local earthquakes for vp and vp=vs velocity models of the central Sumatran forearc. Data were acquired by an amphibious seismometer network consisting of 52 land stations and 10 ocean-bottom seismometers located on a segment of the Sumatran subduction zone that had not ruptured in a great earthquake since 1797 but witnessed recent ruptures to the north in 2005 (Nias earthquake, Mw D 8:7) and to the south in 2007 (Bengkulu earthquake, Mw D 8:5). The 2-D and 3-D vp velocity anomalies reveal the downgoing slab and the sedimentary basins. Although the seismicity pattern in the study area appears to be strongly influenced by the obliquely subducting Investigator Fracture Zone to at least 200 km depth, the 3-D velocity model shows prevailing trench-parallel structures at depths of the plate interface. The tomographic model suggests a thinned crust below the basin east of the forearc islands (Nias, Pulau Batu, Siberut) at ~180 km distance to the trench. vp velocities beneath the magmatic arc and the Sumatran fault zone (SFZ) are around 5 km s-1 at 10 km depth and the vp=vs ratios in the uppermost 10 km are low, indicating the presence of felsic lithologies typical for continental crust. We find moderately elevated vp=vs values of 1.85 at ~150 km distance to the trench in the region of the Mentawai Fault. vp=vs ratios suggest an absence of large-scale alteration of the mantle wedge and might explain why the seismogenic plate interface (observed as a locked zone from geodetic data) extends below the continental forearc Moho in Sumatra. Reduced vp velocities beneath the forearc basin covering the region between the Mentawai Islands and the Sumatra mainland possibly reflect a reduced thickness of the overriding crust. [ABSTRACT FROM AUTHOR]

On 2016 December 25, an Mw 7.6 earthquake broke a portion of the Southern Chilean subduction zone south of Chiloé Island, located in the central part of the Mw 9.5 1960 Valdivia earthquake. This region is characterized by repeated earthquakes in 1960 and historical times with very sparse interseismic activity due to the subduction of a young (~15 Ma), and therefore hot, oceanic plate.We estimate the coseismic slip distribution based on a kinematic finite-fault source model, and through joint inversion of teleseismic body waves and strong motion data. The coseismic slip model yields a total seismic moment of 3.94 x 1020 N·m that occurred over ~30 s, with the rupture propagating mainly downdip, reaching a peak slip of ~4.2 m. Regional moment tensor inversion of stronger aftershocks reveals thrust type faulting at depths of the plate interface. The fore- and aftershock seismicity is mostly related to the subduction interface with sparse seismicity in the overriding crust. The 2016 Chilo´e event broke a region with increased locking and most likely broke an asperity of the 1960 earthquake. The updip limit of the main event, aftershocks, foreshocks and interseismic activity are spatially similar, located ~15 km offshore and parallel to Chiloé Islands west coast. The coseismic slip model of the 2016 Chiloé earthquake suggests a peak slip of 4.2 m that locally exceeds the 3.38 m slip deficit that has accumulated since 1960. Therefore, the 2016 Chiloé earthquake possibly released strain that has built up prior to the 1960 Valdivia earthquake. [ABSTRACT FROM AUTHOR]

The Sumatran subduction zone exhibits strong seismic and tsunamogenic potential with the prominent examples of the 2004, 2005 and 2007 earthquakes. Here, we invert travel time data of local earthquakes for vp and vp/vs velocity models of the central Sumatran forearc. Data were acquired by an amphibious seismometer network consisting of 52 land stations and 10 ocean bottom seismometers located on a segment of the Sumatran subduction zone that had not ruptured in a great earthquake since 1797 but witnessed recent ruptures to the north in 2005 (Nias earthquake, Mw=8.7) and to the south in 2007 (Bengkulu earthquake, Mw=8.5). 2D and 3D vp velocity anomalies reveal the downgoing slab and the sedimentary basins. Although the seismicity pattern in the study area appears to be strongly influenced by the obliquely subducting Investigator Fracture Zone to at least 200 km depth, the 3D velocity model shows prevailing trench parallel structures at depths of the plate interface. The tomographic model suggests a thinned crust below the basin east of the forearc islands (Nias, Pulau Batu, Siberut) at ~ 180 km distance to the trench. Vp velocities beneath the magmatic arc and the Sumatran fault zone SFZ are around 5 km/s at 10 km depth and the vp/vs ratios in the uppermost 10 km are low, indicating the presence of felsic lithologies typical for continental crust. We find moderately elevated vp/vs values of 1.85 at ~ 150 km distance to the trench in the region of the Mentawai fault. Vp/vs ratios suggest absence of large scale alteration of the mantle wedge and might explain why the seismogenic plate interface (observed as a locked zone from geodetic data) extends below the continental forearc Moho in Sumatra. Reduced vp velocities beneath the forearc basin covering the region between Mentawai Islands and the Sumatra mainland possibly reflect a reduced thickness of the overriding crust. [ABSTRACT FROM AUTHOR]

Receiver functions (RF) have been used for several decades to study structures beneath seismic stations. Although most available stations are deployed on shore, the number of ocean bottom station (OBS) experiments has increased in recent years. Almost all OBSs have to deal with higher noise levels and a limited deployment time (∼1 year), resulting in a small number of usable records of teleseismic earthquakes. Here we use OBSs deployed as midaperture array in the deep ocean (4.5-5.5 km water depth) of the eastern mid-Atlantic. We use evaluation criteria for OBS data and beamforming to enhance the quality of the RFs. Although some stations show reverberations caused by sedimentary cover, we are able to identify the Moho signal, indicating a normal thickness (5-8 km) of oceanic crust. Observations at single stations with thin sediments (300-400 m) indicate that a probable sharp lithosphere-asthenosphere boundary (LAB) might exist at a depth of ∼70-80 km which is in line with LAB depth estimates for similar lithospheric ages in the Pacific. The mantle discontinuities at ∼410 km and ∼660 km are clearly identifiable. Their delay times are in agreement with PREM. Overall the usage of beam-formed earthquake recordings for OBS RF analysis is an excellent way to increase the signal quality and the number of usable events. [ABSTRACT FROM AUTHOR]

Our knowledge of the absolute S-wave velocities of the oceanic lithosphere is mainly based on global surface wave tomography, local active seismic or compliance measurements using oceanic infragravity waves. The results of tomography give a rather smooth picture of the actual S-wave velocity structure and local measurements have limitations regarding the range of elastic parameters or the geometry of the measurement. Here, we use the P-wave polarization (apparent P-wave incidence angle) of teleseismic events to investigate the S-wave velocity structure of the oceanic crust and the upper tens of kilometres of the mantle beneath single stations. In this study, we present an up to our knowledge new relation of the apparent Pwave incidence angle at the ocean bottom dependent on the half-space S-wave velocity. We analyse the angle in different period ranges at ocean bottom stations (OBSs) to derive apparent S-wave velocity profiles. These profiles are dependent on the S-wave velocity as well as on the thickness of the layers in the subsurface. Consequently, their interpretation results in a set of equally valid models. We analyse the apparent P-wave incidence angles of an OBS data set which was collected in the Eastern Mid Atlantic. We are able to determine reasonable S-wave-velocity-depth models by a three-step quantitative modelling after a manual data quality control, although layer resonance sometimes influences the estimated apparent S-wave velocities. The apparent S-wave velocity profiles are well explained by an oceanic PREM model in which the upper part is replaced by four layers consisting of a water column, a sediment, a crust and a layer representing the uppermost mantle. The obtained sediment has a thickness between 0.3 and 0.9 km with S-wave velocities between 0.7 and 1.4 km s-1. The estimated total crustal thickness varies between 4 and 10 km with S-wave velocities between 3.5 and 4.3 km s-1. We find a slight increase of the total crustal thickness from ~5 to ~8 km towards the South in the direction of a major plate boundary, the Gloria Fault. The observed crustal thickening can be related with the known dominant compression in the vicinity of the fault. Furthermore, the resulting mantle S-wave velocities decrease from values around 5.5 to 4.5 km s-1 towards the fault. This decrease is probably caused by serpentinization and indicates that the oceanic transform fault affects a broad region in the uppermost mantle. Conclusively, the presented method is useful for the estimation of the local S-wave velocity structure beneath ocean bottom seismic stations. It is easy to implement and consists of two main steps: (1) measurement of apparent P-wave incidence angles in different period ranges for real and synthetic data, and (2) comparison of the determined apparent S-wave velocities for real and synthetic data to estimate S-wave velocity-depth models. [ABSTRACT FROM AUTHOR]

Powerful subduction zone earthquakes rupture thousands of square kilometres along continental margins but at certain locations earthquake rupture terminates. To date, detailed knowledge of the parameters that govern seismic rupture and aftershocks is still incomplete. On 2015 September 16, the Mw 8.3 Illapel earthquake ruptured a 200 km long stretch of the Central Chilean subduction zone, triggering a tsunami and causing significant damage. Here, we analyse the temporal and spatial pattern of the coseismic rupture and aftershocks in relation to the tectonic setting in the earthquake area. Aftershocks cluster around the area of maximum coseismic slip, in particular in lateral and downdip direction. During the first 24 hr after the main shock, aftershocks migrated in both lateral directions with velocities of approximately 2.5 and 5 kmhr-1. At the southern rupture boundary, aftershocks cluster around individual subducted seamounts that are related to the downthrusting Juan Fernández Ridge. In the northern part of the rupture area, aftershocks separate into an upper cluster (above 25 km depth) and a lower cluster (below 35 km depth). This dual seismic-aseismic transition in downdip direction is also observed in the interseismic period suggesting that it may represent a persistent feature for the Central Chilean subduction zone. [ABSTRACT FROM AUTHOR]

The very shallow part of subduction megathrusts occasionally hosts tsunami earthquakes, with unusually slow rupture propagation. The aftershock sequence of the 2010 Mw 8.8 Maule earthquake, offshore Chile, provides us with the opportunity to study systematic changes in source properties for smaller earthquakes within a single segment of a subduction zone. We invert amplitude spectra for double-couple moment tensors and centroid depths of 71 aftershocks of the Maule earthquake down to magnitudes Mw 4.0. In addition, we also derive average source durations. We find that shallower earthquakes tend to have longer normalized source durations on average, similar to the pattern observed previously for larger magnitude events. This depth dependence is observable for thrust and normal earthquakes. The normalized source durations of normal- faulting earthquakes are at the lower end of those for thrust earthquakes, probably because of the higher stress drops of intraplate earthquakes compared to interplate earthquakes. We suggest from the similarity of the depth dependence of normal and thrust events and between smaller and larger magnitude earthquakes that the depth-dependent variation of rigidity, rather than frictional conditional stability at the plate interface, is primarily responsible for the observed pattern. Tsunami earthquakes probably require both low rigidity and conditionally stable frictional conditions; the presence of longduration moderate-magnitude events is therefore a helpful but not sufficient indicator for areas at risk of tsunami earthquakes. [ABSTRACT FROM AUTHOR]

Abstract: On 27 February 2010 the M w 8.8 Maule earthquake in Central Chile ruptured a seismic gap where significant strain had accumulated since 1835. Shortly after the mainshock a dense network of temporary seismic stations was installed along the whole rupture zone in order to capture the aftershock activity. Here, we present the aftershock distribution and first motion polarity focal mechanisms based on automatic detection algorithms and picking engines. By processing the seismic data between 15 March and 30 September 2010 from stations from IRIS, IPGP, GFZ and University of Liverpool we determined 20,205 hypocentres with magnitudes M w between 1 and 5.5. Seismic activity occurs in six groups: 1.) Normal faulting outer rise events 2.) A shallow group of plate interface seismicity apparent at 25–35km depth and 50–120km distance to the trench with some variations between profiles. Along strike, the aftershocks occur largely within the zone of coseismic slip but extend ~50km further north, and with predominantly shallowly dipping thrust mechanisms. Along dip, the events are either within the zone of coseismic slip, or downdip from it, depending on the coseismic slip model used. 3.) A third band of seismicity is observed further downdip at 40–50km depth and further inland at 150–160km trench perpendicular distance, with mostly shallow dipping (~28°) thrust focal mechanisms indicating rupture of the plate interface significantly downdip of the coseismic rupture, and presumably above the intersection of the continental Moho with the plate interface. 4.) A deep group of intermediate depth events between 80 and 120km depth is present north of 36°S. Within the Maule segment, a large portion of events during the inter-seismic phase originated from this depth range. 5.) The magmatic arc exhibits a small amount of crustal seismicity but does not appear to show significantly enhanced activity after the M w 8.8 Maule 2010 earthquake. 6.) Pronounced crustal aftershock activity with mainly normal faulting mechanisms is found in the region of Pichilemu (~34.5°S). These crustal events occur in a ~30km wide region with sharp inclined boundaries and oriented oblique to the trench. The best-located events describe a plane dipping to the southwest, consistent with one of the focal planes of the large normal-faulting aftershock (M w =6.9) on 11 March 2010. [Copyright &y& Elsevier]

Abstract: The Sumatran margin suffered three great earthquakes in recent years (Aceh-Andaman 26 December 2004 Mw=9.1, Nias 28 March 2005 Mw=8.7, Bengkulu 12 September 2007 Mw=8.5). Here we present local earthquake data from a dense, amphibious local seismic network covering a segment of the Sumatran margin that last ruptured in 1797. The occurrence of forearc islands along this part of the Sumatran margin allows the deployment of seismic land-stations above the shallow part of the thrust fault. In combination with ocean bottom seismometers this station geometry provides high quality hypocentre location for the updip end of the seismogenic zone in an area where geodetic data are also available. In this region, the Investigator Fracture Zone (IFZ), which consists of 4 sub-ridges, is subducted below the Sunda plate. This topography appears to influence seismicity at all depth intervals. A well-defined linear streak of seismicity extending from 80 to 200 km depth lies along the prolongation of closely spaced IFZ sub-ridges. More intermediate depth seismicity is located to the southeast of this string of seismicity and is related to subducted rough oceanic seafloor. The plate interface beneath Siberut Island which ruptured last in 1797 is characterised by almost complete absence of seismicity. [ABSTRACT FROM AUTHOR]

The article present the discussion on the depth of seismic faulting at slow- and ultraslow-spreading ridges. Topic include (Grevemeyer et al., 2019) investigated micro-seismicity of the MidCayman Spreading Center and compiled maximum depths of brittle deformation of ten additional slow- and ultraslow-spreading ridges; and Schlindwein (2020) suggests that the aseismic region at the top of the lithosphere might be controlled by hydration of ultramafic rocks.

On April 2014 the Iquique Mw 8.1 earthquake ruptured the interpolate contact between the oceanic Nazca and continental South American plates offshore northern Chile between 19.5°S to 21°S in April 2014. This earthquake did not fully release the strain accumulated since the last great megathrust (Mw 8.8) in 1877 and had left an unbroken segment in the south. From December 2014 to November 2016, we deployed an offshore network of 15 Ocean Bottom Seismometers (OBS) that covered the rupture area and the unbroken southern segment using the Chilean Navy ship OPV Toro and R/V Sonne. That data set is supplemented by five weeks of data from 67 OBS installed for a controlled source seismic experiment during cruise MGL1610 of the R/V Marcus Langseth in late 2016. Data acquired onshore by stations from of the IPOC (Integrated Plate Boundary Observatory Chile) and CSN (Chilean Seismological Service) networks are also included.We present first results of this ongoing project, which include double-difference hypocenter relocations based on waveform cross-correlation. Most of the seismicity occurs between 19.5 and 21°S up-dip of the patch of maximum coseismic slip during the 2014 earthquake, while the seismicity in seismogenic depths is highly concentrated forming well-defined clusters. The observed seismicity provides constraints on the structure of the marine forearc and enables us to relate the seismicity distribution to the background seismicity, seafloor morphology, and regional tectonics. [ABSTRACT FROM AUTHOR]

Large intraplate earthquakes in oceanic lithosphere are rare and usually related to regions of diffuse deformation within the oceanic plate. The 23 January 2018 MW 7.9 strike-slip Gulf of Alaska earthquake ruptured an oceanic fracture zone system offshore Kodiak Island. Bathymetric compilations show a muted topographic expression of the fracture zone due to the thick sediment that covers oceanic basement but the fracture zone system can be identified by offset N-S magnetic anomalies and E-W linear zones in the vertical gravity gradient. Back-projection from global seismic stations reveals that the initial rupture at first propagated from the epicenter to the north, likely rupturing along a weak zone parallel to the ocean crustal fabric. The rupture then changed direction to eastward directed with most energy emitted on Aka fracture zone resulting in an unusual multi-fault earthquake. Similarly, the aftershocks show complex behavior and are related to two different tectonic structures: (1) events along N-S trending oceanic fabric, which ruptured mainly strike-slip and additionally, in normal and oblique slip mechanisms and (2) strike-slip events along E-W oriented fracture zones. To explain the complex faulting behavior we adopt the classical stress and strain partitioning concept and propose a generalized model for large intra-oceanic strike-slip earthquakes of trench-oblique oriented fracture zones/ocean plate fabric near subduction zones. Taking the Kodiak asperity position of 1964 maximum afterslip and outer-rise Coulomb stress distribution into account, we propose that the unusual 2018 Gulf of Alaska moment release was stress transferred to the incoming oceanic plate from co- and post-processes of the nearby great 1964 MW 9.2 megathrust earthquake. [ABSTRACT FROM AUTHOR]

The depth of earthquakes along mid-ocean ridges is restricted by the relatively thin brittle lithosphere that overlies a hot, upwelling mantle. With decreasing spreading rate, earthquakes may occur deeper in the lithosphere, accommodating strain within a thicker brittle layer. New data from the ultraslow-spreading Mid-Cayman Spreading Center (MCSC) in the Caribbean Sea illustrate that earthquakes occur to 10 km depth below seafloor and, hence, occur deeper than along most other slow-spreading ridges. The MCSC spreads at 15 mm/yr full rate, while a similarly well-studied obliquely opening portion of the Southwest Indian Ridge (SWIR) spreads at an even slower rate of ~8 mm/yr if the obliquity of spreading is considered. The SWIR has previously been proposed to have earthquakes occurring as deep as 32 km, but no shallower than 5 km. These characteristics have been attributed to the combined effect of stable deformation of serpentinized mantle and an extremely deep thermal boundary layer. In the context of our MCSC results, we reanalyze the SWIR data and find a maximum depth of seismicity of 17 km, consistent with compilations of spreading-rate dependence derived from slow- and ultraslow-spreading ridges. Together, the new MCSC data and SWIR reanalysis presented here support the hypothesis that depth-seismicity relationships at mid-ocean ridges are a function of their thermal-mechanical structure as reflected in their spreading rate. [ABSTRACT FROM AUTHOR]

Abstract The Earth's ocean floor deforms continuously under the influence of plate tectonic processes. In recent years, the development of deep-sea instruments using acoustic direct-path ranging allows observations of ocean floor deformation with unprecedented spatial and temporal resolution. Due to rapid technological progress, acoustic ranging emerged as a central research field to monitor seafloor deformation. Here we review recent developments and the progress of direct-path ranging applications. We discuss the methodology and examine the effects of the oceanographic environment on the measurement precision. Comparing the resolution of previous deployments, we find that the baseline uncertainty increases linearly with baseline length, at least for distances up to 3 km, but with different linear relations for each deployment. Measurements of displacement at millimeter-level precision across normal, thrust or strike-slip faults are discussed to evaluate the influence of dedicated network designs appropriate for the discrete fault geometries. Furthermore, tectonically quiet areas, such as flanks of coastal or ocean island volcanoes and passive continental margins pose substantial hazards that often lack in-situ monitoring and are therefore a significant target for the application of seafloor geodetic techniques. [ABSTRACT FROM AUTHOR]

Plate tectonics describes oceanic transform faults as conservative strike-slip boundaries, where lithosphere is neither created nor destroyed. Therefore, seafloor accreted at ridge-transform intersections should follow a similar subsidence trend with age as lithosphere that forms away from ridge-transform intersections. Yet, recent compilations of high-resolution bathymetry show that the seafloor is significantly deeper along transform faults than at the adjacent fracture zones. We present residual mantle Bouguer anomalies, a proxy for crustal thickness, for 11 transform fault systems across the full range of spreading rates. Our results indicate that the crust is thinner in the transform deformation zone than in either the adjacent fracture zones or the inside corner regions. Consequently, oceanic transform faulting appears not only to thin the transform valley crust but also leads to a secondary phase of magmatic addition at the transition to the passive fracture zones. These observations challenge the concept of transform faults being conservative plate boundaries. [ABSTRACT FROM AUTHOR]

Splay faults, large thrust faults emerging from the plate boundary to the seafloor in subduction zones, are considered to enhance tsunami generation by transferring slip from the very shallow dip of the megathrust onto steeper faults, thus increasing vertical displacement of the seafloor. These structures are predominantly found offshore, and are therefore difficult to detect in seismicity studies, as most seismometer stations are located onshore. The Mw (moment magnitude) 8.8 Maule earthquake on 27 February 2010 affected ~500 km of the central Chilean margin. In response to this event, a network of 30 ocean-bottom seismometers was deployed for a 3 month period north of the main shock where the highest coseismic slip rates were detected, and combined with land station data providing onshore as well as offshore coverage of the northern part of the rupture area. The aftershock seismicity in the northern part of the survey area reveals, for the first time, a well-resolved seismically active splay fault in the submarine forearc. Application of critical taper theory analysis suggests that in the northernmost part of the rupture zone, coseismic slip likely propagated along the splay fault and not the subduction thrust fault, while in the southern part it propagated along the subduction thrust fault and not the splay fault. The possibility of splay faults being activated in some segments of the rupture zone but not others should be considered when modeling slip distributions. [ABSTRACT FROM AUTHOR]

Subduction megathrust earthquakes show complex rupture behaviour and large lateral variations of slip. However, the factors controlling seismic slip are still under debate. Here, we present 2-D velocity-depth tomographic models across four trench-perpendicular wide angle seismic profiles complemented with high resolution bathymetric data in the area of maximum coseismic slip of the M w 8.8 Maule 2010 megathrust earthquake (central Chile, 34°–36°S). Results show an abrupt lateral velocity gradient in the trench-perpendicular direction (from 5.0 to 6.0 km/s) interpreted as the contact between the accretionary prism and continental framework rock whose superficial expression spatially correlates with the slope-shelf break. The accretionary prism is composed of two bodies: (1) an outer accretionary wedge (5–10 km wide) characterized by low seismic velocities of 1.8–3.0 km/s interpreted as an outer frontal prism of poorly compacted and hydrated sediment, and (2) the middle wedge (∼50 km wide) with velocities of 3.0–5.0 km/s interpreted as a middle prism composed by compacted and lithified sediment. In addition, the maximum average coseismic slip of the 2010 megathrust event is fairly coincident with the region where the accretionary prism and continental slope are widest (50–60 km wide), and the continental slope angle is low (<5°). We observe a similar relation along the rupture area of the largest instrumentally recorded Valdivia 1960 M w 9.5 megathrust earthquake. For the case of the Maule event, published differential multibeam bathymetric data confirms that coseismic slip must have propagated up to ∼6 km landwards of the deformation front and hence practically the entire base of the middle prism. Sediment dewatering and compaction processes might explain the competent rheology of the middle prism allowing shallow earthquake rupture. In contrast, the outer frontal prism made of poorly consolidated sediment has impeded the rupture up to the deformation front as high resolution seismic reflection and multibeam bathymetric data have not showed evidence for new deformation in the trench region. [ABSTRACT FROM AUTHOR]

Epithermal gold (Au) deposits result from the combination of a sustained flux of metal-rich fluids and an efficient precipitation mechanism. Earthquakes may trigger gold precipitation by rapid loss of fluid pressure but their efficiency and time-integrated contribution to gold endowment are poorly constrained. In order to quantify the feedbacks between earthquake-driven fracturing and gold precipitation in the shallow crust, we studied the gold-rich fluids in the active Tolhuaca geothermal system, located in the highly seismic Southern Andes of Chile. We combined temperature measurements in the deep wells with fluid inclusion data, geochemical analyses of borehole fluids and numerical simulations of coupled heat and fluid flow to reconstruct the physical and chemical evolution of the hydrothermal reservoir. The effect of seismic perturbations on fluid parameters was constrained using a thermo-mechanical piston model that simulates the suction pump mechanism occurring in dilational jogs. Furthermore, we evaluated the impact of fluid parameters on gold precipitation by calculating the solubility of gold in pressure ( P )–enthalpy ( H ) space. The reconstructed fluid conditions at Tolhuaca indicate that single-phase convective fluids feeding the hydrothermal reservoir reach the two-phase boundary with a high gold budget (~ 1–5 ppb) at saturated liquid pressures between 20 and 100 bar (210 °C < T sat < 310 °C). We show that if hydrothermal fluids reach this optimal threshold for gold precipitation at a temperature near 250 °C, small adiabatic pressure drops (~ 10 bar) triggered by transient fault-rupture can produce precipitation of 95% of the dissolved gold. Our results at the active Tolhuaca geothermal system indicate that subtle, externally-forced perturbations – equivalent to low magnitude earthquakes (Mw < 2) of a hydrothermal reservoir under optimal conditions – may significantly enhance gold precipitation rates in the shallow crust and lead to overall increases in metal endowment over time. [ABSTRACT FROM AUTHOR]