@article{ZaccarelliBindiStrolloetal.2019, author = {Zaccarelli, Riccardo and Bindi, Dino and Strollo, Angelo and Quinteros, Javier and Cotton, Fabrice Pierre}, title = {Stream2segment: An Open-Source Tool for Downloading, Processing, and Visualizing Massive Event-Based Seismic Waveform Datasets}, series = {Seismological research letters}, volume = {90}, journal = {Seismological research letters}, number = {5}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0895-0695}, doi = {10.1785/0220180314}, pages = {2028 -- 2038}, year = {2019}, abstract = {The task of downloading comprehensive datasets of event-based seismic waveforms has been made easier through the development of standardized webservices but is still highly nontrivial because the likelihood of temporary network failures or subtle data errors naturally increases when the amount of requested data is in the order of millions of relatively short segments. This is even more challenging because the typical workflow is not restricted to a single massive download but consists of fetching all possible available input data (e.g., with several repeated download executions) for a processing stage producing any desired user-defined output. Here, we present stream2segment, a highly customizable Python 2+3 package helping the user in the entire workflow of downloading, inspecting, and processing event-based seismic data by means of a relational database management system as archiving storage, which has clear performance and usability advantages, and an integrated processing subroutine requiring a configuration file and a single Python function to produce user-defined output. Stream2segment can also produce diagnostic maps or user-defined plots, which, unlike existing tools, do not require external software dependencies and are not static images but instead are interactive browser-based applications ideally suited for data inspection or annotation tasks and subsequent training of classifiers in foreseen supervised machine-learning applications. Stream2segment has already been used as a data quality tool for datasets within the European Integrated Data Archive and to create a weak-motion database (in the form of a so-called flat file) for the stable continental region of Europe in the context of the European Ground Shaking Intensity Model service, in turn an important building block for seismic hazard studies.}, language = {en} } @article{MelnickMorenoQuinterosetal.2017, author = {Melnick, Daniel and Moreno, Marcos and Quinteros, Javier and Carlos Baez, Juan and Deng, Zhiguo and Li, Shaoyang and Oncken, Onno}, title = {The super-interseismic phase of the megathrust earthquake cycle in Chile}, series = {Geophysical research letters}, volume = {44}, journal = {Geophysical research letters}, number = {2}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1002/2016GL071845}, pages = {784 -- 791}, year = {2017}, abstract = {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.}, language = {en} } @article{MelnickGarcinQuinterosetal.2012, author = {Melnick, Daniel and Garcin, Yannick and Quinteros, Javier and Strecker, Manfred and Olago, Daniel and Tiercelin, Jean-Jacques}, title = {Steady rifting in northern Kenya inferred from deformed Holocene lake shorelines of the Suguta and Turkana basins}, series = {Earth \& planetary science letters}, volume = {331}, journal = {Earth \& planetary science letters}, number = {10}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-821X}, doi = {10.1016/j.epsl.2012.03.007}, pages = {335 -- 346}, year = {2012}, abstract = {A comparison of deformation rates in active rifts over different temporal scales may help to decipher variations in their structural evolution, controlling mechanisms, and evolution of sedimentary environments through time. Here we use deformed lake shorelines in the Suguta and Turkana basins in northern Kenya as strain markers to estimate deformation rates at the 10(3)-10(4) yr time scale and compare them with rates spanning 10(1)-10(7) yr. Both basins are internally drained today, but until 7 to 5 kyr lake levels were 300 and 100 m higher, respectively, maintained by the elevation of overflow sills connecting them with the Nile drainage. Protracted high lake levels resulted in formation of a maximum highstand shoreline - a distinct geomorphic feature virtually continuous for several tens of kilometers. We surveyed the elevation of this geomorphic marker at 45 sites along >100 km of the rift, and use the overflow sills as vertical datum. Thin-shell elastic and thermomechanical models for this region predict up to similar to 10 m of rapid isostatic rebound associated with lake-level falls lasting until similar to 2 kyr ago. Holocene cumulative throw rates along four rift-normal profiles are 6.8-8.5 mm/yr, or 7.5-9.6 mm/yr if isostatic rebound is considered. Assuming fault dips of 55-65, inferred from seismic reflection profiles, we obtained extension rates of 3.2-6 mm/yr (including uncertainties in field measurements, fault dips, and ages), or 3.5-6.7 mm/yr considering rebound. Our estimates are consistent, within uncertainties, with extension rates of 4-5.1 mm/yr predicted by a modern plate-kinematic model and plate reconstructions since 3.2 Myr. The Holocene strain rate of 10(-15) s(-1) is similar to estimates on the similar to 10(6) yr scale, but over an order of magnitude higher than on the similar to 10(7) yr scale. This is coherent with continuous localization and narrowing of the plate boundary, implying that the lithospheric blocks limiting the Kenya Rift are relatively rigid. Increasing strain rate under steady extension rate suggests that, as the magnitude of extension and crustal thinning increases, the role of regional processes such as weakening by volcanism becomes dominant over far-field plate tectonics controlling the breakup process and the transition from continental rifting to oceanic spreading.}, language = {en} } @article{DuesterhoeftQuinterosOberhaenslietal.2014, author = {Duesterhoeft, Erik and Quinteros, Javier and Oberh{\"a}nsli, Roland and Bousquet, Romain and de Capitani, Christian}, title = {Relative impact of mantle densification and eclogitization of slabs on subduction dynamics: A numerical thermodynamic/thermokinematic investigation of metamorphic density evolution}, series = {Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth}, volume = {637}, journal = {Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0040-1951}, doi = {10.1016/j.tecto.2014.09.009}, pages = {20 -- 29}, year = {2014}, abstract = {Understanding the relationships between density and spatio-thermal variations at convergent plate boundaries is important for deciphering the present-day dynamics and evolution of subduction zones. In particular, the interaction between densification due to mineralogical phase transitions and slab pull forces is subject to ongoing investigations. We have developed a two-dimensional subduction zone model that is based on thermodynamic equilibrium assemblage calculations and includes the effects of melting processes on the density distribution in the lithosphere. Our model calculates the "metamorphic density" of rocks as a function of pressure, temperature and chemical composition in a subduction zone down to 250 km. We have used this model to show how the hydration, dehydration, partial melting and fractionation processes of rocks all influence the metamorphic density and greatly depend on the temperature field within the subduction system. These processes are largely neglected by other approaches that reproduce the density distribution within this complex tectonic setting. Our model demonstrates that the initiation of edogitization (i.e., when crustal rocks reach higher densities than the ambient mantle) of the slab is not the only significant process that makes the descending slab denser and generates the slab pull force. Instead, the densification of the lithospheric mantle of the sinking slab starts earlier than eclogitization and contributes significantly to slab pull in the early stages of subduction. Accordingly, the complex metamorphic structure of the slab and the mantle wedge has an important impact on the development of subduction zones. (C) 2014 Elsevier B.V. All rights reserved.}, language = {en} } @article{BaesSobolevQuinteros2018, author = {Baes, Marzieh and Sobolev, Stephan Vladimir and Quinteros, Javier}, title = {Subduction initiation in mid-ocean induced by mantle suction flow}, series = {Geophysical journal international}, volume = {215}, journal = {Geophysical journal international}, number = {3}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggy335}, pages = {1515 -- 1522}, year = {2018}, abstract = {Pre-existing weakness zones in the lithosphere such as transform faults/fracture zones and extinct mid-oceanic ridges have been suggested to facilitate subduction initiation in an intra-oceanic environment. Here, we propose that the additional forcing coming from the mantle suction flow is required to trigger the conversion of a fracture zone/transform fault into a converging plate boundary. This suction flow can be induced either from the slab remnants of former converging plate boundaries or/and from slabs of neighbouring active subduction zones. Using 2-D coupled thermo-mechanical models, we show that a sufficiently strong mantle flow is able to convert a fracture zone/transform fault into a subduction zone. However, this process is feasible only if the fracture zone/transform fault is very close to the mid-oceanic ridge. Our numerical model results indicate that time of subduction initiation depends on the velocity, domain size and location of mantle suction flow and age of the oceanic plate.}, language = {en} }