ALBERT

All Library Books, journals and Electronic Records Telegrafenberg

Your email was sent successfully. Check your inbox.

An error occurred while sending the email. Please try again.

Proceed reservation?

Export
Filter
  • Seismological Society of America
Collection
  • 1
    Publication Date: 2024-04-22
    Description: Geodetic data in plate boundary zones reflect the accrual of tectonic strain and stress, which will ultimately be released in earthquakes, and so they can provide valuable insights into future seismic hazards. To incorporate geodetic measurements of contemporary deformation into the 2022 revision of the New Zealand National Seismic Hazard Model 2022 (NZ NSHM 2022), we derive a range of strain-rate models from published interseismic Global Navigation Satellite Systems velocities for New Zealand. We calculate the uncertainty in strain rate excluding strain from the Taupō rift–Havre trough and Hikurangi subduction zone, which are handled separately, and the corresponding moment rates. A high shear strain rate occurs along the Alpine fault and the North Island dextral fault belt, as well as the eastern coast of the North Island. Dilatation rates are primarily contractional in the South Island and less well constrained in the North Island. Total moment accumulation derived using Kostrov-type summation varies from 0.64 to 2.93×1019 N·m/yr depending on method and parameter choices. To account for both aleatory and epistemic uncertainty in the strain-rate results, we use four different methods for estimating strain rate and calculate various average models and uncertainty metrics. The maximum shear strain rate is similar across all methods, whereas the dilatation rate and overall strain rate style differ more significantly. Each method provides an estimate of its own uncertainty propagated from the data uncertainties, and variability between methods provides an additional estimate of epistemic uncertainty. Epistemic uncertainty in New Zealand tends to be higher than the aleatory uncertainty estimates provided by any single method, and epistemic uncertainty on dilatation rate exceeds the aleatory uncertainty nearly everywhere. These strain-rate models were provided to the NZ NSHM 2022 team and used to develop fault-slip deficit rate models and scaled seismicity rate models.
    Type: Article , PeerReviewed
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 2
    Publication Date: 2024-01-11
    Description: Following the Mw 6.0 Amatrice earthquake on 24 August 2016 in central Italy, the Emersito task force of the Istituto Nazionale di Geofisica e Vulcanologia installed a temporary seismic network focusing on the investigation of amplification effects at municipalities located on topographic reliefs. Fourteen stations were installed at three sites: Amandola, Civitella del Tronto, and Montereale. During the operational period, stations recorded about 150 earthquakes with Mw up to 4.7. Recorded signals were analyzed calculating the horizon- tal-to-vertical spectral ratios at single station, using both ambient noise and earthquake waveforms, as well as standard spectral ratios (SSRs) to a reference site. To robustly estimate site amplification at each station of the site amplification effect at each station, the influence of backazimuth and epicentral distance is investigated. With the aim of reproducing the observed amplification pattern, 2D numerical simulations were performed on a section orthogonal to the topography major axis, constrained through in situ geological investiga- tions and geophysical surveys. Although at Montereale site no clear amplification effects were observed, at Amandola site, all stations on the relief consistently detected significant peaks at about 4 Hz and along N120–150 azimuth. At Civitella del Tronto, a proper reference station is missing, implying a misleading of site response evaluation in terms of SSRs. Moreover, even if all stations show amplification in the frequency band 1–3 Hz, the direction of the maximum amplification varies from northeast to northwest. At the three sites, obser- vations were successfully reproduced by 2D numerical models, the latter suggesting that topography alone cannot reproduce data, and the interplay with subsoil velocity structure is needed to produce a clear amplification effect. We conclude that according to the previous articles, rather than the sole topography convex shape, the geophysical structure has often a predominant role in controlling the observed amplification pattern on topography.
    Description: Published
    Description: 1208–1229
    Description: OST4 Descrizione in tempo reale del terremoto, del maremoto, loro predicibilità e impatto
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 3
    Publication Date: 2024-01-19
    Description: The triggering of large earthquakes by anthropic activities is a challenging issue in seis- mology, invoked also for the M L 5.9 and 5.8 Emilia 2012 destructive earthquakes. The interaction between the two earthquakes that propagated along adjacent thrusts is still an open issue. In this study, we used waveform cross correlation and double-difference (DD) location methods to precisely relocate the aftershock sequence and get insights into fault geometry, structure, and rheology by means of DD seismic tomography. Accurate relocations highlight a complex fault system with small-length fault segments coalescing in the Mirandola and Ferrara thrusts. We observe a broad continuous high-V P = V S anomaly at seismogenic depth (about 6.0 km) that suggests a possible hydraulic connec- tion along the entire fault system. A close look at seismicity indicates a quasi-simultane- ous activation of the entire thrust system, with the two mainshocks and large aftershocks occurring within the high-V P = V S , high fluid pressure, and connecting volume.
    Description: Published
    Description: 671–684
    Description: OST3 Vicino alla faglia
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 4
    Publication Date: 2024-01-16
    Description: The stress released by large coseismic ruptures and related aftershock sequences is a strongly heterogeneous process. We show highly resolved images of the normal fault system ruptured during the 2016-2017 Central Italy earthquake sequence, as obtained by high resolution local earthquake tomography allowing relocation of a massive set of aftershocks. We get evidence that lateral changes of elastic properties on the fault planes account for the complexity in the rupture processes during the two M〉6 earthquakes. We observe an emergent phase in the first part of the Mw6.5 rupture and the coseismic slip becomes large when the rupture breaks through high Poisson ratio portions of the fault. Mainshocks break fault portions that, although limited and segmented by inherited structural complexity, were dynamically interfering during the faulting episodes. The close repetition of slip on the same relatively high Poisson ratio patch suggests a dynamic weakening of the fault and/or an incomplete stress release during the first mainshock.
    Description: Published
    Description: 2642–2654
    Description: OST3 Vicino alla faglia
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 5
    Publication Date: 2024-01-16
    Description: One important aspect of the seismicity is the spatiotemporal clustering; hence, the distinction between independent and triggered events is a critical part of the analysis of seismic catalogs. Stochastic declustering of seismicity allows a probabilistic distinction between these two kinds of events. Such an approach, usually performed with the epidemic‐type aftershock sequence (ETAS) model, avoids the bias in the estimation of the frequency–magnitude distribution parameters if we consider a subset of the catalog, that is, only the independent or the triggered events. In this article, we present a framework to properly include the probabilities of any event to be independent (or triggered) both in the temporal variation of the seismic rates and in the estimation of the b‐value of the Gutenberg–Richter law. This framework is then applied to a high‐definition seismic catalog in the central part of Italy covering the period from April 2010 to December 2015. The results of our analysis show that the seismic activity from the beginning of the catalog to March 2013 is characterized by a low degree of clustering and a relatively high b‐value, whereas the following period exhibits a higher degree of clustering and a smaller b‐value.
    Description: Published
    Description: 1566–1578
    Description: OST2 Deformazione e Hazard sismico e da maremoto
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 6
    Publication Date: 2024-01-16
    Description: The measurement of earthquake source parameters is affected by large uncertainties, and different approaches lead to large variability in results. One crucial aspect is the trade-off between attenuation (Q) and corner frequency (fc ) in spectral fitting: The source corner frequency, inversely proportional to the fault size, can be severely masked by attenuation and site effects. In this article, we describe a method to solve the trade-off based on the fit of displacement spectra to find the source characteristics (corner frequency, f c , and the signal moment, Ω0) and the single-station attenuation operator (t ), in addition to the site response. We follow a parametric approach based on the use of 3D Q seismic tomography and a bootstrap-based method for selecting the best spectra fit. The correction of attenu- ation with synthetic values derived by 3D attenuation tomography efficiently deals with the trade-off between source and path terms, leading to small uncertainties in the deter- mination of source unknowns (f c and signal moment, Ω0 ), thus yielding constrained esti- mates of source parameters for low- to medium-magnitude earthquakes. We show an application to the Emilia 2012 seismic sequence, for which we computed the source param- eters for 1240 aftershocks (from an initial dataset of 1748) with local magnitude ranging from 2.0 to 4.7 using the spectral fit from P and S waves. About 80% of stress-drop esti- mations are characterized by relatively low uncertainties (within 20% of the estimated values), with maximum values of about 40% for the remaining 20%. The attenuation cor- rection is effective to determine source parameters for small-magnitude earthquakes; hence, we obtain reliable estimates of source parameters for the entire aftershock sequence. This approach gives the opportunity to infer the mechanical state of a complete fault system by taking advantage of the larger number of low-magnitude events (with respect to the largest ones) that always follow a major earthquake.
    Description: Published
    Description: 1739–1758
    Description: OST3 Vicino alla faglia
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 7
    Publication Date: 2024-03-28
    Description: An earthquake sequence occurred in the Central Adriatic region during March–June 2021. This sequence started on 27 March with a mainshock of moment magnitude (Mw) 5.2 occurring at 13:47 coordinated universal time (UTC). No foreshock was observed before this mainshock. The sequence lasted approximately three months, until the end of June 2021. Approximately 200 seismic events were recorded by the regional seismic network during this time, including four M ≥ 4.0 earthquakes. The 27 March 2021 earthquake was one of the strongest instrumentally recorded events in the area bounded approximately by the Ancona–Zadar line to the north and the Gargano–Dubrovnik line to the south. The mainshock originated at a focal depth of 9.9 km. The seismicity spread from the mainshock up-dip and down-dip along a northeast-dipping plane. Here, we investigate the geometry of the fault activated by this seismic sequence by using sP depth phases. We aim to significantly reduce the large uncertainties associated with the hypocentral locations of offshore earthquakes beneath the Adriatic Sea—an area that plays a fundamental role in the geodynamics of the Mediterranean. These refined earthquake locations also allow us to make inferences with regards to the seismotectonic context responsible for the analyzed seismicity, thus identifying a structure (here referred to as the MidAdriatic fault) consisting of a northwest–southeast-striking thrust fault with a ∼ 35° northeast-dipping plane. The use of depth-phase arrival times to constrain off-network event locations is of particular interest in Italy due to both the peculiar shape of the peninsula and the extreme scarcity of seafloor stations, the cost and management of which are very expensive and complex. Here, we present the first attempt to apply this off-network locating technique to the Italian offshore seismicity research with the aim of improving hazard estimations in these hard-to-monitor regions.
    Description: Published
    Description: 480–493
    Description: OST3 Vicino alla faglia
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 8
    Publication Date: 2024-04-11
    Description: Volcanic activity produces a broad spectrum of seismic and acoustic signals whose characteristics provide important clues on the underlying magmatic processes. Networks and arrays of seismic and acoustic sensors are the backbone of most modern volcano monitoring programs. Investigation of the signals gathered by these instruments requires efficient workflows and specialist software. The high sampling rates, typically 50 Hz or greater, at which seismic and acoustic waveforms are recorded by multistation networks and dense arrays leads to the rapid accumulation of large volumes of data, making the implementation of efficient data analysis workflows for volcano surveillance a challenging task. Here, we present an open‐source MATLAB graphical user interface, MISARA (Matlab Interface for Seismo‐Acoustic aRray Analysis), designed to provide a user‐friendly workflow for the analysis of seismoacoustic data in volcanic environments. MISARA includes efficient algorithm implementations of well‐established techniques for seismic and acoustic data analysis. It is designed to support visualization, characterization, detection, and location of volcano seismoacoustic signals. Its intuitive, modular, structure facilitates rapid, semiautomated, inspection of data and results, thus reducing user effort. MISARA was tested using seismoacoustic data recorded at Etna Volcano (Italy) in 2010, 2011, and 2019, and is intended for use in education and research, and to support routine data analysis at volcano observatories.
    Description: Published
    Description: 1689–1702
    Description: OSV4: Preparazione alle crisi vulcaniche
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 9
    Publication Date: 2024-03-19
    Description: The 1695 Asolano earthquake (Mw 6.4) is the southernmost of the six largest earthquakes to have occurred in NE Italy or nearby (the others being: 1348, Mw 6.6, Eastern-Alps; 1511, Mw 6.3, Friuli-Slovenia; 1873, Mw 6.2, Alpago-Cansiglio; 1936, Mw 6.1, Alpago-Cansiglio; 1976, Mw 6.4, Friuli). The 1695 earthquake is generally associated with the Montello thrust, most recent studies locating it on the eastern slope of the Montello hill. A full-scale reappraisal of all available historical data leads this study to a more robust macroseismic localization of the 1695 earthquake and to open toward other possible locations of the seismic source that produced it. In particular, it becomes feasible to place its epicenter at the foothills of the Monte Grappa massif, the major morphological expression of the Bassano-Valdobbiadene thrust fault. Here we describe the reasons that make this fault a possible alternative to previous hypotheses.
    Description: Published
    Description: 526–538
    Description: OST2 Deformazione e Hazard sismico e da maremoto
    Description: JCR Journal
    Keywords: Asolo earthquake ; Southern Alps ; Historical Seismology, Geology and deformation
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 10
    Publication Date: 2024-03-12
    Description: The macroseismic source parameters of earthquakes occurring within a sequence are strongly influenced by cumulative damage effects. When we deal with historical seismic sequences, in addition to the cumulative intensities, other intrinsic uncertainties due to the scarcity and indeterminacy of sources come into play. These issues imply that the parameterizations of the single earthquakes within a historical seismic sequence are not univocal and that all the uncertainties that are addressed when assessing macroseismic intensity should be carefully considered in the parameter estimation. In the light of these considerations, we performed some tests on the 2016–2017 and 1703 seismic sequences, which occurred in the same area in central Italy, to compute the macroseismic source parameters by means of two independent methods. Results show that the cumulative effects arising from multiple damaging earthquakes can cause biases in the intensity assessments, which affect the computed magnitude and epicentral locations. To reduce bias in macroseismic intensities due to cumulative damage, we illustrate a simple procedure, called cumulative intensity subtraction (CIS), which consists in discarding the localities strongly damaged by the early earthquakes of a sequence from the intensity distributions used for computing the macroseismic source parameters of the subsequent earthquakes. The outcomes show that, for the 2016 seismic sequence, the CIS approach provides locations in agreement with the instrumental epicenters and with the causative faults. For the 1703 sequence, the CIS approach along with explicit accounting for the indeterminacy in intensity assignments give a range of equally plausible solutions. The CIS represents an exploration of a simple strategy that stems from an attempt to give significance to macroseismic intensity in the presence of cumulative damage.
    Description: Published
    Description: 759–774
    Description: OST4 Descrizione in tempo reale del terremoto, del maremoto, loro predicibilità e impatto
    Description: JCR Journal
    Keywords: macroseismic intesity ; cumulative effects ; microseismic source parameters ; 04.06. Seismology
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 11
    Publication Date: 2024-04-08
    Description: The region of central Italy is well known for its moderate-to-large earthquakes. Events such as 2016 Mw 6.2 Amatrice, generated in the shallow extensional tectonic regime, motivate numerical simulations to gain insights into source-related ground-motion complexities. We utilize a hybrid integral–composite kinematic rupture model by Gallovič and Brokešová (2007) to predict ground motions for other hypothetical Amatrice fault rupture scenarios (scenario events). The synthetic seismograms are computed in 1D crustal velocity models, including region-specific 1D profiles for selected stations up to 10 Hz. We create more than ten thousand rupture scenarios by varying source parameters. The resulting distributions of synthetic spectral accelerations at periods 0.2–2 s agree with the empirical nonergodic ground-motion model ofSgobbaetal.(2021)forcentral Italy in terms of the mean and total variability. However, statistical mixed-effect analysis of the residuals indicates that the between-eventvariability of the scenarios exceeds theempirical one significantly. We quan tify the role ofsourcemodelparametersinthemodelinganddemonstratethepivotalroleof theso-called stress parameterthatcontrols high-frequencyradiation. Weproposerestricting thescenariovariability tokeepthebetween-eventvariabilitywithintheempiricalvalue.The presented validation of the scenario variability can be generally utilized in scenario model ing for more realistic physics-based seismic hazard assessment.
    Description: In press
    Description: OST3 Vicino alla faglia
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 12
    Publication Date: 2024-05-21
    Description: One of the challenges of seismicity monitoring is to achieve multiparametric catalogs complete down to small magnitude using automatic procedures. This can be obtained using seismic networks with high performance and robust, automatic algorithms able to process large data sets, limiting the manual operations of the analysts. The characterization of microseismicity is fundamental to study its spatial and temporal evolution and to define the seismic activity of fault systems. Among the source parameters of microseismic events, focal mechanisms are not generally calculated and, when available in the seismic catalog, their reliability may be dubious. We propose a new tool, named Tool for automatic Earthquake low‐frequency Spectral Level estimAtion (TESLA), to automatically calculate the P‐ and S‐wave low‐frequency spectral levels. Indeed, it has been shown that these levels can be inverted together with P‐phase polarities to better constrain the focal mechanism or to estimate the seismic moment. TESLA is designed to invert the P‐ and S‐displacement spectra searching the optimal signal window to use for the spectral analysis. Using a signal window of fixed duration, although variable according to the earthquake magnitude, is not always the appropriate choice, especially when microseismicity is analyzed. TESLA performs three main tasks for both P and S phases: (1) a systematic exploration of several signal windows to use for the computation of displacement spectra, (2) the spectral analysis for all the selected signal windows, and (3) the evaluation of the best‐displacement spectra through quantitative criteria and the estimation of the low‐frequency spectral levels. The tool is first validated and then applied to the 2013 St. Gallen, Switzerland, induced seismic sequence to calculate the P and S low‐frequency spectral level ratios, which are inverted to estimate focal mechanisms. Our results show the robustness of the tool to process microseismicity and the benefit of using it to automatically analyze large waveform data sets.
    Description: Published
    Description: 2441-2455
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: manuscript
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 13
    Publication Date: 2024-05-21
    Description: No abstract
    Description: Published
    Description: 557-561
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 14
    facet.materialart.
    Unknown
    Seismological Society of America
    In:  Das, R., M. L. Sharma, H. R. Wason, D. Choudhury, and G. Gonzales (2019). A seismic moment magnitude scale, Bull Seismol. Soc. Am. 109, no. 4, 1542–1555, doi: 10.1785/0120180338.
    Publication Date: 2024-05-21
    Description: Moment magnitude Mw was first defined by Hiroo Kanamori in the late 1970s, when the availability of new force balance seismometers made it possible to measure the seismic moment M0 with virtually no limits in the frequency passband. For this reason, Mw does not become saturated even for the largest earthquakes ever recorded. Mw has been chosen in such a way that it coincides best with the previous definitions of magnitude (Ms, ML, mb, etc.) on certain ranges of values but can deviate significantly from them within other ranges. A few years ago, Das and colleagues proposed a new moment magnitude scale Mwg with the aim of better reproducing the values of mb and Ms over their entire range and to better predict the energy ES radiated by earthquakes. We show that there was no need to define such a new scale and that Mwg is not even optimal to achieve the goal of matching ES.
    Description: In press
    Description: OST2 Deformazione e Hazard sismico e da maremoto
    Description: JCR Journal
    Keywords: earthquake magnitude ; moment magnitude scale ; 04.06. Seismology
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 15
    facet.materialart.
    Unknown
    Seismological Society of America
    In:  "Earthquake Magnitude Conversion Problem” by Ranjit Das, H. R. Wason, Gabriel Gonzalez, M. L. Sharma, Deepankar Choudhury, Conrad Lindholm, Narayan Roy, and Pablo Salazar
    Publication Date: 2024-05-21
    Description: Similar to the previous ones, the latest paper by Das and Colleagues (Das et al.,2018) on the application of the general orthogonal regression (GOR) method (Fuller, 1987; Castellaro et al.,2006), for the conversions between different types of earthquake magnitudes, is a collection of incorrect or undemonstrated assertions, most of which have already been pointed out in several contributions that have been published in the last few years (Gasperini and Lolli, 2014a, b; Gasperini et al., 2015, 2018; Pujol, 2018). We recall below only some of them. According to the recent seismological literature, we use here the term “GOR” to indicate the errors-in variable regression method described by Fuller (1987), even if such term is not fully in line with mathematical statistics as orthogonality is only given for equal errors of the dependent and independent variables.
    Description: Published
    Description: 1366-1369
    Description: OST2 Deformazione e Hazard sismico e da maremoto
    Description: JCR Journal
    Keywords: 04.06. Seismology
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 16
    Publication Date: 2024-05-16
    Description: The new generation seismic hazard maps use 3D seismotectonic fault models, which are more consistent with the actual nature of faults, whereas the classical models based on earthquake catalogs only utilize a 2D representation of the seismicity. Although the former provides more reliable information on seismogenic structures, the latter can deliver trustworthy seismicity rates easily. Therefore, it is necessary to combine both the approaches to create a high‐quality seismic hazard assessment model. This study proposes an innovative approach using smoothed seismicity methods that can be advantageous in all contexts with available 3D fault models and high‐quality seismic catalogs. We applied our method on the Adriatic Basal Thrust (ABT) in eastern central Italy—a lithospheric‐scale active contractional structure with a well‐constrained 3D geometric–kinematic reconstruction and a related high‐quality catalog. Our new 3D algorithm was applied to smooth the ABT seismicity on the grid, resulting in a 3D earthquake rate model that also provides rupture parameters such as strike, dip, rake, and seismogenic thickness. Our approach is particularly useful for complex seismotectonic settings, such as in cases of lithospheric shear zones, subduction planes, and overlapping multidepth seismogenic volumes.
    Description: Published
    Description: 10-20
    Description: OST2 Deformazione e Hazard sismico e da maremoto
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 17
    Publication Date: 2024-04-30
    Description: The global premodern descriptive catalogs of von Hoff (1840–1841), Perrey (1845–1850), and Mallet (1853–1855), all covering the period from Antiquity to 1842, have been used to populate the preinstrumental section of parametric catalogs since the late 1960s. The earthquake lists of these three authors have been individually analyzed to compile a comprehensive inventory of the sources on which they relied. Conversely to previous analyses and uses that were focused—primarily if not exclusively—on their seismological content, both in the compilation of the inventory and in this article the hundreds of items supplying earthquake records are in the foreground. After having merged these three sets of sources and having obtained a comprehensive list of about 5000 earthquakes and two times more source entries, similarities and differences are evidenced. This comprehensive analysis is meant to catch and explain how the preference given to one source among many available may have affected the interpretation of the collected records and influenced the accuracy and reliability of the earthquake lists by von Hoff, Perrey, and Mallet, and consequently of their contribution to the construction of the global modern parametric earthquake catalogs.
    Description: Published
    Description: 2456-2468
    Description: OST2 Deformazione e Hazard sismico e da maremoto
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 18
    Publication Date: 2024-05-06
    Description: This article presents the first publicly available version of the NExt STrOng Related Earthquake (NESTORE) software (NESTOREv1.0) designed for the statistical analysis of earthquake clusters. NESTOREv1.0 is a MATLAB (www.mathworks.com/products/ matlab , last accessed August 2022) package capable of forecasting strong aftershocks starting from the first hours after the mainshocks. It is based on the NESTORE algorithm, which has already been successfully applied retrospectively to Italian and California seismicity. The code evaluates a set of features and uses a supervised machine learning approach to provide probability estimates for a subsequent large earthquake during a seismic sequence. By analyzing an earthquake catalog, the software identifies clusters and trains the algorithm on them. It then uses the training results to obtain forecasting for a test set of independent data to estimate training performance. After appropriate testing, the software can be used as an Operational Earthquake Forecasting (OEF) method for the next stronger earthquake. For ongoing clusters, it provides near-real-time forecasting of a strong aftershock through a traffic light classification aimed at assessing the level of concern. This article provides information about the NESTOREv1.0 algorithm and a guide to the software, detailing its structure and main functions and showing the application to recent seismic sequences in California. By making the NESTOREv1.0 software available, we hope to extend the impact of the NESTORE algorithm and further advance research on forecasting the strongest earthquakes during seismicity clusters.
    Description: Published
    Description: 2003–2013
    Description: OST4 Descrizione in tempo reale del terremoto, del maremoto, loro predicibilità e impatto
    Description: OST5 Verso un nuovo Monitoraggio
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 19
    Publication Date: 2024-05-29
    Description: The Båth’s law is an empirical seismological relation between the magnitudes of the mainshock and the largest aftershock of a seismic sequence. This empirical law, differently from other seismological laws, could be valid only when the seismic sequence is ended. Indeed, during the sequence, we don’t know if the strongest event has already happened or not, and then inferring something about the magnitude of the largest aftershock is hazardous. In this opinion paper, we discuss some issues related to the Båth’s law: its validity on a global catalog, the use of the terms “mainshock” and “aftershock” in the seismological community and in the public, and their implications in earthquake forecasting communications. We show the uselessness of Båth’s law in earthquake forecasting, and that the words “mainshock” and “aftershock” have different interpretations for the public and for seismologists. We argue that their use during an ongoing seismic sequence, without a proper explanation of their meaning, could be confounding, in particular for the Italian language. Calling all events just “earthquakes” could be a simple but effective solution.
    Description: Published
    Description: 2565–2568
    Description: OST4 Descrizione in tempo reale del terremoto, del maremoto, loro predicibilità e impatto
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 20
    Publication Date: 2024-05-22
    Description: We present a new approach to estimate the predominant direction of rupture propagation during a seismic sequence. A fast estimation of the rupture propagation direction is essential to knowthe azimuthal distribution of shaking around the seismic source and the associated risks for the earthquake occurrence. The main advantage of the proposed method is that it is conceptually reliable, simple, and fast (near real time). The approach uses the empirical Green’s function technique and can be applied directly to the waveforms without requiring the deconvolution of the instrumental response and without knowing a priori the attenuation model and the orientation of the activated fault system. We apply the method to the 2016–2017 Amatrice-Visso-Norcia high-energy and long-lasting earthquake series in central Italy,which affected a large area up to 80 kmalong strike, withmore than 130,000 events of small-to-moderate magnitude recorded until the end of August 2022. Most of the selected events analyzed in this study have a magnitude greater than 4.4 and only four seismic events have a magnitude in the range of 3.3–3.7. Our results show that the complex activated normal fault system has a rupture direction mainly controlled by the pre-existing normal faults and by the orientation of the reactivated faults. In addition, the preferred direction of rupture propagation is also controlled by the presence of fluid in the pre-existing structural discontinuities. We discuss the possible role of fluids as a cause of bimaterial interface. Another important finding from our analysis is that the spatial evolution of seismicity is controlled by the directivity.
    Description: Published
    Description: 1912–1924
    Description: JCR Journal
    Keywords: directivity ; 04.06. Seismology
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 21
    Publication Date: 2024-05-22
    Description: We investigate the variability of Brune stress drop (Δσ), apparent stress (τa), and Savage– Wood radiation efficiency (ηsw   τa= Δσ), in the 2013–2014 Mw 5.0 earthquake sequence that struck the Matese area in the southern Apennines range of Italy. The sequence is clustered in a relatively small crustal volume in the 13–22 km depth range, which is greater than that of background seismicity and normal-faulting sequences that occurred under the range axis, usually located in the first 15 km of the crust.We find high Savage– Wood radiation efficiency values for most of the analyzed earthquakes located in a narrow crustal volume, with values ranging from well above the self-similarity value to very high values as high as 0.55. In addition, a large variability in radiation efficiency (up to 90%) is observed for two similar magnitude events at different depths. Previous studies reported seismic evidence of fluid involvement in the nucleation process of the Matese earthquakes. By integrating our results with crustal geophysical data published recently, we propose that most of the earthquakes characterized by high values of ηsw are nucleated within high pore pressure zones located in the crystalline midcrust of Adria. We reckon that high pore pressure fluids of deep origin played a role in the rupture process and were responsible for themixed shear-tensile sources inferred from the analysis of the S-wave/P-wave spectral amplitude ratio for most of 2013–2014 earthquakes.
    Description: Published
    Description: 299–319
    Description: JCR Journal
    Keywords: stress drop ; apparent stress ; radiation efficiency
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 22
    Publication Date: 2024-05-22
    Description: Mefite d’Ansanto (Italy) is a nonvolcanic field characterized by persistent strong degassing activity. A seismic field monitoring carried out during the Summer 2021 reveals a persis- tent, extended, and complex source of seismic tremor characterized by a spectrum with a frequency content from about 1 Hz to more than 35 Hz. While at frequency smaller than 3 Hz the signal amplitude is stationary, in the intermediate frequency band (3–20 Hz) sud- den changes of amplitude are often observed, suggesting the existence of an intermittent source (every few minutes to tens of minutes). Furthermore, very short bursts of high-fre- quency energy are recognized in the tremor signal. Results of array analysis and seismo- logical observation indicate that the sources of the analyzed tremor are located in a small area centered on the main vent of the degassing area. The persistent low-frequency tremor and the intermediate frequency signals propagate as surface waves to the seismic stations installed around the source and indicate a very shallow source. On the contrary, impulsive signals at frequencies greater than 20 Hz propagate as body waves, revealing a deeper source likely located between 50 and 100 m depth.
    Description: Published
    Description: 1102–1114
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 23
    Publication Date: 2024-06-13
    Description: We aim to compute macroseismic parameters (location and magnitude) using the BOXER code for the first time on the citizen testimonies, that is, individual intensity data points (IDPs) at the global scale collected and made available by the LastQuake system of the European–Mediterranean Seismological Centre (EMSC). IDPs available for different earthquakes are selected to eliminate those that are geographically inconsistent with most data; then they are clustered spatially based on various methods. For each cluster with at least three IDPs, a macroseismic data point (MDP), corresponding to an intensity value assessed for given localities as in classical macroseismic studies, is computed by various central tendency estimators (average, median, and trimmed averages). Finally, macroseismic parameters are obtained by MDP distribution using two location methods of BOXER code. For each earthquake, we used raw and corrected intensities and 132 different combinations of grouping methods, estimators, and BOXER methods. We assigned a ranking to the combinations that best reproduce instrumental parameters and used such a ranking to select preferred combinations for each earthquake. We analyzed retrospectively the reliability of the parameters as a function of time and space. The results are essentially identical using original and corrected intensities and show higher reliability for BOXER’s method 1 than for method 0; they are dependent on the geographical area, and generally improve over time and with the number of IDPs collected. These findings are useful for the future real-time analyses, and for evaluating the location and magnitude of earthquakes whenever a sufficient number of IDPs are available and with a distribution such that MDPs can be derived and the BOXER method applied.
    Description: Published
    Description: 969–996
    Description: OST4 Descrizione in tempo reale del terremoto, del maremoto, loro predicibilità e impatto
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 24
    Publication Date: 2024-06-13
    Description: Underdetermination is a condition affecting all problems in seismic imaging. It manifests mainly in the nonuniqueness of the models inferred from the data. This condition is exacerbated if simplifying hypotheses like isotropy are discarded in favor of more realistic anisotropic models that, although supported by seismological evidence, require more free parameters. Investigating the connections between underdetermination and anisotropy requires the implementation of solvers which explore the whole family of possibilities behind nonuniqueness and allow for more informed conclusions about the interpretation of the seismic models. Because these aspects cannot be investigated using traditional iterative linearized inversion schemes with regularization constraints that collapse the infinite possible models into a unique solution, we explore the application of transdimensional Bayesian Monte Carlo sampling to address the consequences of underdetermination in anisotropic seismic imaging. We show how teleseismic waves of P and S phases can constrain upper‐mantle anisotropy and the amount of additional information these data provide in terms of uncertainty and trade‐offs among multiple fields.
    Description: Published
    Description: 1214–1226
    Description: OST1 Alla ricerca dei Motori Geodinamici
    Description: JCR Journal
    Keywords: 04.06. Seismology
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 25
    Publication Date: 2024-06-13
    Description: The Algerian offshore earthquake of 18 March 2021, Mw 6.0, was felt by people in various Italian regions, also at large epicentral distance. This unusual human perception far from the source prompted us to analyze the waveforms recorded by land seismic stations installed along the Iberian, French, and Italian coasts. On some seismograms of the selected network, prominent T phases are detected. T waves can travel in the SOund Fixing And Ranging (SOFAR) channel over great distances (thousands of kilometers) with little loss in signal strength and be recorded by near‐coastal seismometers after the P (primary) and S (secondary) phases (hence T or tertiary phases). To explain the subjective perception of ground shaking with quantities that are measured on the seismogram, we estimated the empirical macroseismic intensities for both body and T phases and we calculated the body‐wave seismic attenuation. The P‐wave anelastic attenuation analysis shows two main wave propagation patterns that reflect lithosphere heterogeneity of the Algerian, Liguro‐Provençal, and Tyrrhenian basins. We find that in some cases, in particular along the Italian and French coasts, the largest ground shaking is caused by the T phase. Our observations confirm that the central‐western Mediterranean Sea is a favorable site for T‐wave propagation and suggest that the T phases should be taken into account in ground‐shaking hazard assessment for the central‐western Mediterranean.
    Description: Published
    Description: 859–869
    Description: OST2 Deformazione e Hazard sismico e da maremoto
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 26
    Publication Date: 2024-06-13
    Description: Werecomputethecoefficients of the intensity prediction equation (IPE) in Italy using the data of the DBMI15 version 2.0 (v.2.0) intensity database and the instrumental and combined (instrumental plus macroseismic) magnitudes reported by the CPTI15 v.2.0 catalog. Wefollow the same procedure described in the previous article, consisting of a first step in which the attenuation of intensity I with respect to the distance D from macroseismic hypocenter is referred to the expected intensity at the epicenter IE and a second step in which IE is related to the instrumental magnitude Mi, the combined magnitude Mc,the epicentral intensity I0, and the maximum intensity Imax using error-in-variable (EIV) regression methods. The main methodological difference with respect to the original article concerns the estimation of the uncertainty of IE to be used for EIV regressions, which is empirically derived from the standard deviation of regression between IE and Mi and also used for the regressions of IE with Mc, I0,andImax. In summary, the new IPE determined from DBMI15 v.2.0 is I IE−0:0081D−h−1:072 lnD−lnh , in which D p R2 h2 , h = 4.49 km, and IE can be calculated from the intensity data distribution of the earthquake. If the intensity data distribution is not available, IE can be calculated from the following relationships IE −2:578 1:867Mw,IE I0.
    Description: In press
    Description: OST2 Deformazione e Hazard sismico e da maremoto
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 27
    Publication Date: 2024-06-13
    Description: We analyze the interplay between hydrology, deformation and seismicity in the Matese massif, located in the Italian Southern Apennines. We find that this area is characterized by the concurrent action of two hydrologically-driven processes: the first is the deformation detected by GNSS data in the shallowest part (above the elevation of the major springs) of the Earth crust, in phase with the hydrological forcing; the second is the triggering of seismicity at depth with a delay suggesting a downward diffusive process. We study the first process by applying a Principal Component Analysis to the GNSS displacements time series, aiming to identify a common signal describing the largest data variance. We find that the maximum horizontal displacements associated with the first principal component (PC1) are larger than 1 cm in two GNSS sites and the PC1 temporal evolution is well correlated and in phase with the flow of the largest spring of the region, which we consider as proxy of the water content of the massif. This suggests that the main source of horizontal deformation is the water content fluctuations in the shallow portion of the Matese aquifer, in particular within fractures located in correspondence of the main mapped faults. The deformation rates caused by this process are one order of magnitude larger than the tectonic ones. Finally, we infer the second process by observing the correlation between the background seismicity and the spring discharge with a time lag of 121 days. In our interpretation, downward diffusive processes, driven by aquifer water content variations, propagate pore pressure waves that affect the faults strength favoring the occurrence of micro-earthquakes. This is supported by the values of hydraulic diffusivity (1.5 m^2/s) and rock permeability (3.2-3.8⋅10^−13 m^2), which are compatible with what is observed in karstified limestones.
    Description: Published
    Description: 1899–1912
    Description: OST2 Deformazione e Hazard sismico e da maremoto
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 28
    Publication Date: 2023-10-03
    Description: In part 1, we run multiple GIT decomposition for different choices of model assumptions, namely three different window duration for Fourier calculation, two different parametrization of the attenuation, two different site constraints. We also considered different source models (Brune, Boatwright, Brune with kappa_source) and different approaches to estimate uncertainties of source parameters (i.e., considering the covariance matrix, Monte Carlo sampling of the residual distribution, model selection with threshold based on F-test).
    Description: As part of the community stress-drop validation study initiative, we apply a spectral decomposition approach to isolate the source spectra of 556 events occurred during the 2019 Ridgecrest sequence (Southern California). We perform multiple decompositions by introducing alternative choices for some processing and model assumptions, namely: three different S-wave window durations (i.e., 5 s, 20 s, and variable between 5 and 20 s); two attenuation models that account differently for depth dependencies; and two different site amplification constraints applied to restore uniqueness of the solution. Seismic moment and corner frequency are estimated for the Brune and Boatwright source models, and an extensive archive including source spectra, site amplifications, attenuation models, and tables with source parameters is disseminated as the main product of the present study. We also compare different approaches to measure the precision of the parameters expressed in terms of 95% confidence intervals (CIs). The CIs estimated from the asymptotic standard errors and from Monte Carlo resampling of the residual distribution show an almost one-to-one correspondence; the approach based on model selection by setting a threshold for misfit chosen with an F-ratio test is conservative compared to the approach based on the asymptotic standard errors. The uncertainty analysis is completed in the companion article in which the outcomes from this work are used to compare epistemic uncertainty with precision of the source parameters.
    Description: Published
    Description: 1980–1991
    Description: 3T. Fisica dei terremoti e Sorgente Sismica
    Description: JCR Journal
    Keywords: source parameters ; GIT ; uncertainties ; moment magnitude ; corner frequency ; 04.06. Seismology
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 29
    Publication Date: 2023-10-03
    Description: This paper is the second part of a previous publication (Bindi et al. 2023 DOI: 10.1785/0220230019). All the decompositions of part 1 are organized in a sort of logic tree and mixed-effect regressions are performed to partition the variability into contributions related to duration, attenuation and site-constraint grouping factors. Statistical uncertainties computed in part 1 (i.e., coming from the fit) are compared with epistemic uncertainties associated to the logic tree, and Sammon's maps are used to visualize the impact of the grouping factors on the overall shape of the source spectra.
    Description: As part of the community stress‐drop validation study, we evaluate the uncertainties of seismic moment M0 and corner frequency fc for earthquakes of the 2019 Ridgecrest sequence. Source spectra were obtained in the companion article by applying the spectral decomposition approach with alternative processing and model assumptions. The objective of the present study is twofold: first, to quantify the impact of different assumptions on the source parameters; and second, to use the distribution of values obtained with different assumptions to estimate an epistemic contribution to the uncertainties. Regarding the first objective, we find that the choice of the attenuation model has a strong impact on fc results: by introducing a depth‐dependent attenuation model, fc estimates of events shallower than 6 km increase of about 10%. Also, the duration of the window used to compute the Fourier spectra show an impact on fc ⁠: the average ratio between the estimates for 20 s duration to those for 5 s decreases from 1.1 for Mw〈3 to 0.66 for Mw〉4.5. For the second objective, we use a mixed‐effect regression to partition the intraevent variability into duration, propagation, and site contributions. The standard deviation ϕ of the intraevent residuals for log(fc) is 0.0635, corresponding to a corner frequency ratio 102ϕ=1.33. When the intraevent variability is compared to uncertainties on log(fc), we observe that 2ϕ is generally larger than the 95% confidence interval of log(fc), suggesting that the uncertainty of the source parameters provided by the fitting procedure might underestimate the model‐related (epistemic) uncertainty. Finally, although we observe an increase of log(Δσ) with log(M0) regardless of the model assumptions, the increase of Δσ with depth depends on the assumptions, and no significant trends are detected when depth‐dependent attenuation and velocity values are considered.
    Description: Published
    Description: 1992–2002
    Description: 3T. Fisica dei terremoti e Sorgente Sismica
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 30
    facet.materialart.
    Unknown
    Seismological Society of America
    In:  Taroni, M., J. Zhuang, and W. Marzocchi (2021). High-definition mapping of the Gutenberg–Richter b-value and its relevance: A case study in Italy, Seismol. Res. Lett. 92, 3778–3784, doi: 10.1785/0220210017.
    Publication Date: 2023-02-21
    Description: Taroni et al. (2021) published a statistical framework to reliably estimate the b-value and its uncertainties, with the goal being the interpretation in a seismotectonic context and improving earthquake forecasting capabilities. In this comment, we show that the results presented for the Italian region and the conclusions drawn by the authors, are heavily biased due to quarry-blast events in the Italian earthquake catalog used in the analysis. Without removing this anthropogenic component in the data, a meaningful analysis of the earthquake- size distribution for natural seismicity is, in our opinion, not possible. This comment highlights the need for basic data quality analysis before sophisticated statistical tools are applied to a dataset.
    Description: European Union’s Horizon 2020 research and innovation program under Grant Agreement Number 821115 Pianeta Dinamico-Working Earth INGV-MUR project.
    Description: Published
    Description: 1089-1094
    Description: 5T. Sismologia, geofisica e geologia per l'ingegneria sismica
    Description: JCR Journal
    Keywords: 04.06. Seismology
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 31
    Publication Date: 2023-02-21
    Description: Machine‐learning (ML) methods have seen widespread adoption in seismology in recent years. The ability of these techniques to efficiently infer the statistical properties of large datasets often provides significant improvements over traditional techniques when the number of data are large (millions of examples). With the entire spectrum of seismological tasks, for example, seismic picking and detection, magnitude and source property estimation, ground‐motion prediction, hypocenter determination, among others, now incorporating ML approaches, numerous models are emerging as these techniques are further adopted within seismology. To evaluate these algorithms, quality‐controlled benchmark datasets that contain representative class distributions are vital. In addition to this, models require implementation through a common framework to facilitate comparison. Accessing these various benchmark datasets for training and implementing the standardization of models is currently a time‐consuming process, hindering further advancement of ML techniques within seismology. These development bottlenecks also affect “practitioners” seeking to deploy the latest models on seismic data, without having to necessarily learn entirely new ML frameworks to perform this task. We present SeisBench as a software package to tackle these issues. SeisBench is an open‐source framework for deploying ML in seismology—available via GitHub. SeisBench standardizes access to both models and datasets, while also providing a range of common processing and data augmentation operations through the API. Through SeisBench, users can access several seismological ML models and benchmark datasets available in the literature via a single interface. SeisBench is built to be extensible, with community involvement encouraged to expand the package. Having such frameworks available for accessing leading ML models forms an essential tool for seismologists seeking to iterate and apply the next generation of ML techniques to seismic data.
    Description: Published
    Description: 1695–1709
    Description: 3T. Fisica dei terremoti e Sorgente Sismica
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 32
    Publication Date: 2023-03-13
    Description: Taroni et al. (2021; hereafter TZM21) proposed a method to perform a spatial b-value mapping based on the weighted-likelihood estimation and applied this method to the Italian region as a tutorial example. In the accompanying comment, Gulia et al. (2021; hereafter GGW21) did not challenge the TZM21’s method, but they argued that the catalog used by TZM21 is contaminated by quarry blasts, introducing a bias that may impact any seismotectonic or hazard interpretations. Although in TZM21 the application to the Italian territory was only a tutorial example and we purposely did not make any thorough discussion on the meaning of the results in terms of seismotectonic or seismic hazards (that would have required many more analyses), we acknowledge the potential role of the quarry blasts, and we add some further analysis here. We thank GGW21 for giving us this opportunity. Here, removing the part of the catalog contaminated by quarry blasts and applying the same analysis as in TZM21, we obtain results that are very similar to the ones reported in TZM21; specifically, only one region that is characterized by low natural seismicity rate shows a marked effect of the quarry blasts on the b-value.
    Description: Published
    Description: 195-197
    Description: 6T. Studi di pericolosità sismica e da maremoto
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 33
    Publication Date: 2023-01-23
    Description: The detection level of a seismic network is a measure of its effective ability to record small earthquakes in a given area. It can vary in both space and time and depends on several factors such as meteorological conditions, anthropic noise, local soil conditions—all factors that affect the seismic noise level—as well as the quality and operating condition of the instruments. The ability to estimate the level of detection is of tremendous importance both in the design of a new network and in determining whether a given network can recognize seismicity consistently or needs to be improved in some of its parts. In this article, we determine the detection level of the Cuban seismic network using the empirically estimated seismic noise spectral level at each station site and some theoretical relationships to predict the signal amplitude of a seismic event at individual stations. The minimum local detectable magnitude thus depends on some network parameters such as the signal‐to‐noise ratio and the number of stations used in the calculation. We also demonstrate the effectiveness of our predictions by comparing the estimated detection level with those empirically determined from one year of data (i.e., the year 2020) of the Cuban seismic catalog. Our analysis shows, on the one hand, in which areas the current Cuban network should be improved, also depending on the regional pattern of faults, and, on the other hand, indicates the magnitude threshold that can be assumed homogeneously for the catalog of Cuban earthquakes in 2020. Because the adopted method can use current measurements of the seismic noise level (e.g., daily), the proposed analysis can also be configured for continuous monitoring of network state quality.
    Description: Published
    Description: 2048-2062
    Description: 5T. Sismologia, geofisica e geologia per l'ingegneria sismica
    Description: 1IT. Reti di monitoraggio e sorveglianza
    Description: JCR Journal
    Keywords: seismic monitoring ; detection ; cuba ; seismic network ; Event Detection Level ; 04.06. Seismology
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 34
    Publication Date: 2023-01-27
    Description: We study the spectral decay parameter κ using S-wave recordings from the central Italy dense regional array. The data set used consists of 266 earthquakes, 353 stations, and 13,952 observations of κ with a mean value of 0.0412 ± 0.0177 within the distance range of 7.1–168.8 km. We model the variation of κ with hypocenter distance r as κ r κ0 κs ̃κ r , in which κ0 and κs represent the near-site and the near-source decay parameters, respectively, and ̃κ r the average κ along the S-wave source-station paths. We first determine ̃κ r with a nonparametric inversion approach and then we solved for κ0 and κs with a second inversion. We found that ̃κ r increases with distance within the whole distance range analyzed (9.2–80.6 km). The near-source decay parameter takes values in the range 0:0 〈 κs ≤ 0:026 with a mean value of 0.003 ± 0.006, which represents 7.52% of the mean value of the observed κ. The values of the near-site decay parameter vary in the range 0:0035 ≤ κ0 ≤ 0:0823 with a mean value of 0.0298 ± 0.0133, that is, 72.28% of the mean value of the κ observed. We conclude that most of the high- frequency attenuation takes place near the site, because ̃κ r contributes with only 20.2% of the spectral decay. We also investigate the spatial variability of κ by determining ̃κ r within four quadrants that divide the studied region taking as a reference axis the Apennines chain orientation. We found higher values of ̃κ r in the southern quad- rants, where seismicity and faulting are more active, and less attenuation in the more stable northeast quadrant.
    Description: Published
    Description: 2299–2310
    Description: 4T. Sismicità dell'Italia
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 35
    Publication Date: 2023-01-27
    Description: In this article, we apply the reference-rock identification method (RRIM; Lanzano et al., 2020) to the ITalian ACcelerometric Archive, which includes more than 1600 recording stations in Italy and in the neighboring countries, with different levels of site characterization. The RRIM is based on the identification and the evaluation of site parameters representing the reference site conditions (Steidl et al., 1996) and the construction of the scoring scheme to classify the candidate stations. Given the large number of sites, the preselection of can- didates is performed via residual analysis, selecting those characterized by flat site response and amplitude similar or lower than the one for the generic rock (average shear-wave veloc- ity in the uppermost 30 m, V S30 800 m= s). The main results of this study are: (1) a list of reference rock sites in Italy, with an associated score; (2) a scenario-independent generic-to- reference rock corrective factor for the ground-motion model for shallow active crustal events in Italy (ITA18; Lanzano, Luzi, et al., 2019); (3) a model for the generic-to-reference rockcorrectivefactor,parametrizedintermsofVS30 andκ0,thatis,thehigh-frequencydecay parameter (Anderson and Hough, 1984). A collateral product is a set of coefficients for the prediction of 81 ordinates of the Fourier amplitude spectra (FAS) in the frequency interval 0.1–30 Hz, calibrated with the same dataset and functional form of ITA18 for acceleration response spectra (SA). The application of RRIM allowed us to identify 116 stations with aver- age measured V S30 ∼ 900 m= s. The corrective factor allows to scale both SA and FAS spectra, and has a significant effect at high frequencies, reducing the ground motion by up to a factor 1.7 at f = 10 Hz. The introduction of κ0 in the corrective term modeling is effective from 2 Hz onward and results in a reduction of variability at high frequencies (f 〉 10 Hz).
    Description: Published
    Description: 1583–1606
    Description: 5T. Sismologia, geofisica e geologia per l'ingegneria sismica
    Description: 6T. Studi di pericolosità sismica e da maremoto
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 36
    Publication Date: 2023-03-01
    Description: Minimum 1D velocity models and station corrections have been computed for the central Mediterranean area using two main data sets. The first one consists of accurate first arrival‐time readings from 103 seismic events with magnitude (ML)≥3.5 recorded by the Italian National Seismic Network (RSN) and the AlpArray Seismic Network (AASN) in the period 2014–2021. Earthquakes were selected on the basis of their spatial distribution, epicentral distance to the nearest seismic station, and maximum distance traveled by Pn and Sn phases. This fine selection of high‐quality data combined with the spatial density of the AlpArray seismic stations was decisive in obtaining high resolution for upper mantle velocity, especially in the Alpine belt. To obtain a denser coverage of crustal rays, we extended the first data set with P and S arrivals of local earthquakes from Istituto Nazionale di Geofisica e Vulcanologia (INGV) bulletin data (2016–2018). A total of 75,807 seismic phases (47,183 P phases and 28,264 S phases) have been inverted to calculate best‐fit 1D velocity models, at regional and local scales. We then test the performance of the optimized velocity models by relocating the last four years of seismicity recorded by INGV (period 2017–2020). The computed velocity models are very effective for routine earthquake location, seismic monitoring, source parameter modeling, and future 3D seismic tomography.
    Description: Published
    Description: 2670--2685
    Description: 4T. Sismicità dell'Italia
    Description: JCR Journal
    Keywords: geophysics ; velocity models ; Italian seismicity ; central mediterranean area ; 04. Solid Earth ; 04.06. Seismology
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 37
    Publication Date: 2023-03-31
    Description: We investigate the dependence of the Gutenberg–Richter b parameter on the crustal thickness quantified by the Moho depth, for nine different regional catalogs. We find that, for all the catalogs considered in our study, the b‐value is larger in areas presenting a thicker crust. This result appears in apparent contradiction with previous findings of a b decreasing with the focal depth. However, both the results are consistent with acoustic emission experiments, indicating a b‐value inversely proportion to the applied differential stress. Our results can be indeed interpreted as the signature of a larger stress concentration in areas presenting a thinner crust. This is compatible with the scenario where postseismic deformation plays a central role in stress concentration and in aftershock triggering.
    Description: Published
    Description: 1921–1934
    Description: 6T. Studi di pericolosità sismica e da maremoto
    Description: JCR Journal
    Keywords: b-value ; 04.06. Seismology
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 38
    Publication Date: 2023-05-25
    Description: The strategy for data processing in the Engineering Strong-Motion Database (ESM) is to disseminate only manually revised data to ensure the highest quality. However, manual processing is no longer sustainable, due to the ever-increasing rate of digital earth-quake records, from global, regional, and national seismic networks, and a new frame-work for strong-motion data processing is required, so that records are automatically processed and the human revision is restricted to selected significant records. To this end, we present ESMpro—a modular Python software for a renewed processing frame-work of ESM. The software is available in a stand-alone beta version to facilitate testing and sharing among the scientific community. ESMpro provides automatic settings for waveform trimming and filtering, along with the automatic recognition of poor-quality data and multiple events. ESMpro allows classifying each record in different quality classes to reduce manual revision on a subset of the incoming data. ESMpro also allows handling different processing techniques in a modular and flexible structure to facilitate the implementation of new or alternative algorithms and file formats. The testing performed on the ESM database results in a good correspondence between the automatic and manual data processing, supporting the migration toward fully automatic procedures for massive data processing.
    Description: Published
    Description: 961–974
    Description: 5T. Sismologia, geofisica e geologia per l'ingegneria sismica
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 39
    Publication Date: 2023-05-25
    Description: We present a probabilistic seismic hazard analysis (PSHA) for the entire Po Plain sedimentary basin (Italy)—one of the widest Quaternary alluvial basins of Europe, to evaluate the impact of site-response characterization on hazard estimates. A large-scale application of approach 3 of the U.S. Nuclear Regulatory Commission (NRC) to include seismic amplification in the hazard is presented. Both 1D amplification related to stratigraphic conditions and 3D amplification due to basin effects are considered with the associated uncertainties, and their impact on the hazard is analyzed through a sensitivity analysis. Whereas 3D basin effects are considered through the application of an empirical, spatial invariant correction term, 1D amplification was estimated throughout the study area by means of dynamic (equivalent linear) ground-response analysis. To separate aleatory variabilities and epistemic uncertainties related to site response, a partially nonergodic approach is used. The results provide a finer picture of the actual seismic hazard, highlighting those areas where the ground motion is affected by amplification effects due to local or regional geological features. We found that, for a return period of 475 yr, neglecting basin effects produces a 30% underestimation of the seismic hazard in the long-period ( 〉 1 s) range. Moreover, with reference to the hazard model adopted, such effects are responsible for most of the epistemic uncertainty (up to 80%) in the results. Therefore, such effects deserve special attention in future research related to PSHA in the Po Plain sedimentary basin.
    Description: Published
    Description: 1269–1285
    Description: 5T. Sismologia, geofisica e geologia per l'ingegneria sismica
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 40
    Publication Date: 2023-01-04
    Description: Strong‐motion recordings of the 2014 Mw 8.2 Iquique and 2015 Mw 8.3 Illapel, Chile, earthquakes were analyzed to determine rupture propagation and the location, timing, and strength of subevents that produce most of the high‐frequency (≥1 Hz) ground motions. A moving window,cross‐correlation analysis of recordings from a local dense array, band‐pass filtered at 1 Hz, directly shows that the Iquique earthquake ruptured to the southeast over a distance of about 60 km. Array analysis of lower frequency energy (0.03–0.1 Hz) indicates that it occurred updip of the high‐frequency rupture. A methodology was developed for inverting the envelopes of acceleration records (1–5 Hz) to map high‐frequency source factors on the rupture zone and was applied to the two earthquakes. Waveforms of Mw 6 earthquakes were used as empirical Green’s functions in the inversions. High‐frequency subevents within the two Mw 8 earthquakes were located at depths ranging from 25 to 55 km and mostly occurred downdip of the peak slip of these earthquakes. Fourier spectral ratios of the Iquique mainshock with respect to Mw 5–6 aftershocks were fit to determine their stress drops. The stress drops were roughly constant from Mw 5 to 8 at 10–20 MPa. A compound rupture model is described in which subevents occur in areas of spatially heterogeneous strength and stress on the rupture, and produce the high‐frequency radiated energy of the overall earthquake, but are not located in the areas of peak slip. The stress drop of the overall earthquake is shown to equal the root mean square stress drop of subevents averaged over the rupture area.
    Type: Article , PeerReviewed
    Format: text
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 41
    Publication Date: 2022-12-27
    Description: We investigated the temporal variation of the spectral decay parameter κ before and after two main earthquakes that occurred in the central Italy region, namely the Amatrice (Mw 6.0) of 24 August 2016 and the Norcia (Mw 6.5) of 30 October 2016 earthquakes. For this analysis, we used seismograms from the central Italy dense seismic array stations, and earthquakes located at hypocenter distances r 〈 80 km, having magnitudes Mw 3.4–6.5. The dataset consists of 393 events recorded at 92 stations. We estimated, for both earthquake sequences, average functions κ˜(r) that describe the distance dependence of κ along the S-wave source-station paths using acceleration spectra from foreshocks, mainshock, and aftershocks. We observed that there was a regional attenuation drop within approximately two months after the Amatrice earthquake. Then, κ˜(r) tends to return toward the attenuation values observed before the occurrence of the main event, namely to the values of κ˜(r) obtained from the foreshocks, when the earthquake cycle is probably completed. We also estimated the near-source kappa (κs) using aftershocks from 24 August 2016 to 3 September 2016. The results show that the values of κs are lower than those from aftershocks located to the north near the epicenter of the Amatrice earthquake, suggesting that the tectonic stress was probably high near the rupture zone, and that there may be a likely fluid flow of crustal fluids. κ˜(r) obtained from the foreshocks of the Norcia earthquake is like that calculated with the records of the Amatrice aftershocks. Then, κ˜(r) drops to lower attenuation values during the Norcia main event and tends to increase again during the aftershocks. From the analysis of these two earthquake sequences that occurred in a short-time interval in central Italy, we conclude that the temporal variation of κ˜(r) could be a valuable indicator to monitor the earthquake cycle.
    Description: Published
    Description: 3037–3045
    Description: 5T. Sismologia, geofisica e geologia per l'ingegneria sismica
    Description: JCR Journal
    Repository Name: Istituto Nazionale di Geofisica e Vulcanologia (INGV)
    Type: article
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 42
    Publication Date: 2021-03-24
    Description: In exploration and earthquake seismology, most sources used in subsurface structure imaging and rock property estimation are fixed in certain positions. Continuously moving seismic sources, such as vehicles and the metro, are one kind of important passive sources in ambient noise research. Commonly, seismic data acquisition and processing for moving sources are based on the assumption of simple point passive sources, and the dispersion curve inversion is applied to constrain near-surface velocity. This workflow neglects the Doppler effects. Considering the continuously moving properties of the sources, we first derive the analytical solution for the Rayleigh waves excited by heavy vehicles and then analyze their Doppler effects and dispersion curves. We observe that the moving source data have the Doppler effect when compared with the changes in the frequency of the source intensity, but this effect does not affect the frequency dispersion of Rayleigh waves. The dispersion curves computed for moving source records are consistent with the analytical dispersion solutions, which provide a theoretical foundation for velocity estimation using moving source data.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 43
    Publication Date: 2021-03-24
    Description: Seismograms are convolution results between seismic sources and the media that seismic waves propagate through, and, therefore, the primary observations for studying seismic source parameters and the Earth interior. The routine earthquake location and travel-time tomography rely on accurate seismic phase picks (e.g., P and S arrivals). As data increase, reliable automated seismic phase-picking methods are needed to analyze data and provide timely earthquake information. However, most traditional autopickers suffer from low signal-to-noise ratio and usually require additional efforts to tune hyperparameters for each case. In this study, we proposed a deep-learning approach that adapted soft attention gates (AGs) and recurrent-residual convolution units (RRCUs) into the backbone U-Net for seismic phase picking. The attention mechanism was implemented to suppress responses from waveforms irrelevant to seismic phases, and the cooperating RRCUs further enhanced temporal connections of seismograms at multiple scales. We used numerous earthquake recordings in Taiwan with diverse focal mechanisms, wide depth, and magnitude distributions, to train and test our model. Setting the picking errors within 0.1 s and predicted probability over 0.5, the AG with recurrent-residual convolution unit (ARRU) phase picker achieved the F1 score of 98.62% for P arrivals and 95.16% for S arrivals, and picking rates were 96.72% for P waves and 90.07% for S waves. The ARRU phase picker also shown a great generalization capability, when handling unseen data. When applied the model trained with Taiwan data to the southern California data, the ARRU phase picker shown no cognitive downgrade. Comparing with manual picks, the arrival times determined by the ARRU phase picker shown a higher consistency, which had been evaluated by a set of repeating earthquakes. The arrival picks with less human error could benefit studies, such as earthquake location and seismic tomography.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 44
    Publication Date: 2021-03-24
    Description: Two events of magnitude (mb) 3.6–3.8 occurred in southern North Korea (NK) on 27 June 2019 and 11 May 2020. Although these events were located ∼330–400  km from the known nuclear test site, the fact that they occurred within the territory of NK, a country with a recent history of underground nuclear tests, made them events of interest for the monitoring community. We used P/Lg ratios from regional stations to categorize seismic events that occurred in NK from 2006 to May 2020, including these two recent events, the six declared NK nuclear tests, and the cavity collapse and triggered earthquakes that followed the 3 September 2017 nuclear explosion. We were able to separate the cavity collapse from the population of nuclear explosions. However, based on P/Lg ratios, the distinction between the earthquakes and the cavity collapse is ambiguous. The performed discriminant analyses suggest that combining Pg/Lg and Pn/Lg ratios results in improved discriminant power compared with any of the ratio types alone. We used the two ratio types jointly in a quadratic discriminant function and successfully classified the six declared nuclear tests and the triggered earthquakes that followed the September 2017 explosion. Our analyses also confirm that the recent southern events of June 2019 and May 2020 are both tectonic earthquakes that occurred naturally.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 45
    Publication Date: 2021-04-01
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 46
    Publication Date: 2021-03-23
    Description: We use deep learning to predict surface-to-borehole Fourier amplification functions (AFs) from discretized shear-wave velocity profiles. Specifically, we train a fully connected neural network and a convolutional neural network using mean AFs observed at ∼600 KiK-net vertical array sites. Compared with predictions based on theoretical SH 1D amplifications, the neural network (NN) results in up to 50% reduction of the mean squared log error between predictions and observations at sites not used for training. In the future, NNs may lead to a purely data-driven prediction of site response that is independent of proxies or simplifying assumptions.
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 47
    Publication Date: 2021-03-23
    Description: We used teleseismic P and S waves recorded in the course of the 2016 Incorporated Research Institutions for Seismology (IRIS) community-planned experiment in northern Oklahoma, to estimate amplitude correction factors (ACFs) and orientation correction factors (OCFs) for the gradiometer’s three-component Fairfield nodal sensors and two other gradiometer-styled subarray nodal sensors. These subarrays were embedded in the 13 km aperture nodal array that was also fielded during the 2016 IRIS experiment. The array calibration method we used in this study is based on the premise that a common wavefield should be recorded over a small-aperture array using teleseismic observation. In situ estimates of ACF for the gradiometer vary by 2.3% (standard deviation) for the vertical components and, typically, variability is less than 4.3% for the horizontal components; associated OCFs generally dispersed by 3°. For the two subarrays, the vertical-component ACF usually vary up to 2.4%; their horizontal-component ACFs largely spread up to 3.6%. OCFs for the subarrays generally disperse by 6.5°. ACF and OCF estimates for the gradiometer are seen to be stable across frequency bands having high signal coherence and/or signal-to-noise ratio. Gradiometry analyses of calibrated and uncalibrated gradiometer records from a local event revealed notable improvements in accuracy of attributes obtained from analyzing the calibrated horizontal-component waveforms in the light of catalog epicenter-derived azimuth. The improved waveform relative amplitudes after calibration, coupled with the enhanced wave attribute accuracy, suggests that instrument calibration for amplitude statics and orientation errors should be encouraged prior to doing gradiometry analysis in future studies.
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 48
    Publication Date: 2021-03-23
    Description: Dynamic strains have never played a role in determining local earthquake magnitudes, which are routinely set by displacement waveforms from seismic instrumentation (e.g., ML). We present a magnitude scale for local earthquakes based on broadband dynamic strain waveforms. This scale is derived from the peak root-mean-squared strains (A) in 4589 records of dynamic strain associated with 365 crustal earthquakes and 77 borehole strainmeters along the Pacific-North American plate boundary on the west coast of the United States and Canada. In this data set, catalog moment magnitudes range from 3.5≤Mw≤7.2, and hypocentral distances range from 6≤R≤500  km. The 1D representation of geometrical spreading and attenuation of A common to all strain data is logA0(R)=−0.00072R−1.45log(R). After correcting for instrument gain, site terms, and event terms, the magnitude scale, MDS=logA−logA0(R)−log(3×10−9), scales as ≈0.92Mw with a residual standard deviation of 0.19. This close association with Mw holds for events east of the −124° meridian; west of this boundary, however, a constant correction of 0.41 is needed to adjust for additional along-path attenuation effects. As a check on the accuracy of this magnitude scale, we apply it to dynamic strain records from three strainmeters located in the near field of the 2019 M 6.4 and 7.1 Ridgecrest earthquakes. Results from these six records are in agreement to within 0.5 magnitude units, and five out of six records are in agreement to within 0.34 units.
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 49
    Publication Date: 2021-03-24
    Description: We present a Hellenic database of intensity measures from uniformly processed strong ground motion recordings, together with metadata on earthquake source attributes and recording site conditions. The database consists of information from 471 earthquakes between 1973 and 2015 that produced 2993 usable recordings from 333 sites. A key element of this work is a unified presentation of data from two major data providers that operate in Greece (Institute of Engineering Seismology and Earthquake Engineering and the Institute of Geodynamics, National Observatory of Athens) along with a university-operated local urban array (University of Patras). Consistent procedures were applied to develop source parameters that include hypocenter locations, moment magnitudes (directly estimated or derived using a conversion procedure), fault-plane solutions, and finite-fault parameters (generally, for events with M〉6.0). The time-averaged shear-wave velocity in the upper 30 m parameter is provided for all 333 recording sites based on geophysical measurements where available (102) and proxy-based estimates otherwise. Most events are in the magnitude range of 3.8–7, occurred at shallow hypocentral depths (
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 50
    Publication Date: 2021-03-24
    Description: Earthquake early warning (EEW) is becoming an increasingly attractive real-time strategy for mitigating the threats posed by potentially devastating incoming seismic events. As efforts accelerate to develop practical EEW-based solutions for earthquake-prone countries in Europe, it is important to understand and quantify the level of performance that can be achieved by the underlying seismological algorithms. We conduct a conceptual study on EEW performance in Europe, which explicitly focuses on the accuracy and associated uncertainties of selected methodological approaches. Twenty-three events from four diverse European testbeds are used to compare the quality of EEW predictions produced by the Virtual Seismologist and PRobabilistic and Evolutionary early warning SysTem algorithms. We first examine the location and magnitude estimates of the algorithms, accounting for both bias and uncertainty in the resulting predictions. We then investigate the ground-shaking prediction capabilities of the source-parameter estimates, using an error metric that can explicitly capture the propagation of uncertainties in these estimates. Our work highlights the importance of accounting for EEW parameter uncertainties, which are often neglected in studies of EEW performance. Our findings can be used to inform current and future implementations of EEW systems in Europe. In addition, the evaluation metrics presented in this work can be used to determine EEW accuracy in any worldwide setting.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 51
    Publication Date: 2021-03-24
    Description: For observations of vertical-component acceleration in the normal-mode band (0.3–10 mHz), the detection sensitivity for signals from the Earth’s body can be improved to levels below the Peterson low-noise model (PLNM). This is achieved by deterministic procedures that (at least partly) remove the accelerations originating from atmospheric mass fluctuations. The physical models used in such corrections are still too simple and fail at frequencies above 3 mHz. Anticipating improved atmospheric correction procedures, we explore the prospects of lowering the detection level. From recordings of excellent vertical-component sensors operated under exceptional site conditions at the Black Forest Observatory, we select time windows of very low background signal, for which all of the contributing broadband seismometers showed their best performance. Streckeisen seismometers of type STS-1, STS-2, and STS-6A, a Nanometrics Trillium T360, and the superconducting gravimeter (SG) SG056 manufactured by GWR Instruments take part in this comparison. Because of their low level of self-noise, the STS-1 and the SG056-G1 benefit the most from a correction with the best currently available improved Bouguer plate model for atmospherically induced signals at frequencies below 1 mHz. As far as we know, this is the first case in which the background level of a broadband seismometer could be lowered below the PLNM. At signal periods beyond the normal-mode band (investigated up to 12 hr), the gravimeters show the lowest level of self-noise, directly followed by the STS-6A. In the band from 0.3 to 10 mHz, the STS-1 has the lowest level of self-noise, which is at least 4 dB below the PLNM, directly followed by the T360 and the STS-6A. Sensors of lower self-noise than the currently manufactured STS-6A or T360 are needed before improved atmospheric correction procedures lead to a significantly lower vertical-component detection threshold.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 52
    Publication Date: 2021-03-24
    Description: Developing new ways to observe tsunami contributes to tsunami research. Tidal and deep-ocean gauges are typically used for coastal and offshore observations. Recently, tsunami-induced ground tilts offer a new possibility. The ground tilt signal accompanied by 2010 Mw 8.8 Chilean earthquake were observed at a tiltmeter network in Japan. However, tiltmeter stations are usually not as widely installed as broadband seismometers in other countries. Here, we studied broadband seismic records from Japan’s F-net and found ground tilt signals consistent with previously published tiltmeter dataset for this particular tsunamic event. Similar waveforms can also be found in broadband seismic networks in other countries, such as Taiwan, as well as an ocean-bottom seismometer. We documented a consistent time sequence of evolving back-azimuth directions of the tsunami waves at different stages of tsunami propagation through beamforming-frequency–wavenumber analysis and particle-motion analysis; the outcomes are consistent with the tsunami propagation model provided by the Pacific Tsunami Warning Center. These results shown that dense broadband seismic networks can provide a useful complementary dataset, in addition to tiltmeter arrays and other networks, to study or even monitor tsunami propagation using arrayed methods.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 53
    Publication Date: 2021-03-24
    Description: Thunderstorms are a common atmospheric phenomenon that cause abundant acoustic disturbances, which can interact with the ground surface, creating a link between atmospheric and solid Earth processes. This article reports seismological observations of four thunderstorms through the spring and summer of 2019, as recorded by the distributed acoustic sensing fiber-optic array (4.9 km) on the Penn State campus in State College, Pennsylvania. With a dense sensor array in the local region, we are able to construct the seismic full waveform response of the thunderstorm events (hereafter referred to as thunderquakes) and track the wave propagation across the array. We use a time-domain grid search to obtain the back azimuth and slowness of the waves, and a modified Geiger’s method to pinpoint source locations of the thunderquakes. Correlated with the time of the recorded signal, this data allows reconstruction of thunderstorm movement as well as offering measurements of the seismic velocity.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 54
    Publication Date: 2021-03-24
    Description: We present a novel physics-guided neural network to estimate shear-tensile focal mechanisms for microearthquakes using displacement amplitudes of direct P waves. Compared with conventional data-driven fully connected (FC) neural networks, our physics-guided neural network is implemented in an unsupervised fashion and avoids the use of training data, which may be incomplete or unavailable. We incorporate three FC layers and a scaling and shifting layer to estimate shear-tensile focal mechanisms for multiple events. Then, a forward-modeling layer, which generates synthetic amplitude data based on the source mechanisms emerging from the previous layer, is added. The neural network weights are iteratively updated to minimize the mean squared error between observed and modeled normalized P-wave amplitudes. We apply this machine-learning approach to a set of 530 induced events recorded during hydraulic-fracture simulation of Duvernay Shale west of Fox Creek, Alberta, yielding results that are consistent with previously reported source mechanisms for the same dataset. A distinct cluster characterized by more complex mechanisms exhibits relatively large Kagan angles (5°–25°) compared with the previously reported best double-couple solutions, mainly due to model simplification of the shear-tensile focal mechanism. Uncertainty tests demonstrate the robustness of the inversion results and high tolerance of our neural network to errors in event locations, the velocity model, and P-wave amplitudes. Compared with a single-event grid-search algorithm to estimate shear-tensile focal mechanisms, the proposed neural network approach exhibits significantly higher computational efficiency.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 55
    Publication Date: 2021-03-10
    Description: The rapid availability of reliable damage statistics, after the occurrence of a major earthquake, is an essential mitigation strategy to drive and support emergency intervention operations. Unfortunately, the latency in collecting and organizing actual damage information has a substantial impact on the efficiency of the initial phases of the intervention framework. To speedup preliminary management operations, a quick, although, coarse prediction of the expected damage is highly desirable. For this purpose, we have developed a system for rapid damage estimation. The system, presently implemented for the Friuli Venezia Giulia region, relies on the existing seismological monitoring infrastructure of the National Institute of Oceanography and Applied Geophysics (OGS), which is responsible for delivering earthquake alerts in northeastern Italy. In case of a major earthquake event, the predicted damage is automatically computed using the OpenQuake software engine by means of ad hoc structural exposure and fragility models developed for the region. Damage calculations rely on a combination of actual observed ground motion from the stations of the OGS seismological network and empirical prediction using the ShakeMaps software developed by the U.S. Geological Survey. The resulting damage scenario, aggregated at municipality level, is finally delivered to the control room of the regional civil protection in support of early intervention activities. Although, the system is presently still under active development, a number of experimental trials have confirmed the reliability and the usefulness of the proposed approach. We are confident that the current research will contribute in mitigating the impact of possible future damaging earthquakes by (1) guiding targeted postevent emergency interventions, (2) increasing the preparedness and response capacity of emergency teams and population through preparatory training activities, and (3) supporting the decision-making process during the recovery phase, hence enhancing resilience.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 56
    Publication Date: 2021-03-10
    Description: The Istituto Nazionale di Geofisica e Vulcanologia (INGV) is in charge of earthquake monitoring and surveillance in the Italian territory as a part of the civil protection system. Technological improvements in the last years were taken into account for developing new protocols and software to upgrade all the procedures in the monitoring centers. Real-time earthquake evaluation consists of phase picks, preliminary and automatic hypocenters, local magnitudes and ground-motion parameters. The real-time analysis system presently in use at INGV is the starting point for a new multitier compound system that relies on four main components: an automatic earthquake detection and location system based on Earthworm; a new seismological relational database for parametric data; a full set of new webservices application programming interface specifications to share information and provide data at the application level and; finally, a set of multiplatform interactive revision tools developed to analyze, store, use, and distribute the seismic parameters in real time. These last three components are being completely developed ex-novo at INGV in Rome. Such a system has been engineered to communicate with the International Federation of Digital Seismic Networks standard webservices as well as custom home-made INGV services. Through its custom embedded features Caravel will allow the INGV personnel on duty for seismic surveillance to evaluate and review all automatic estimations before they are communicated to the Italian civil protection agency and then published by external tools through e-mail, SMS, Twitter, and webpages.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 57
    Publication Date: 2021-03-10
    Description: The National Observatory of Athens data center for the European Integrated Data Archive (EIDA@NOA) is the national and regional node that supports International Federation of Digital Seismograph Networks and related webservices for seismic waveform data coming from the southeastern Mediterranean and the Balkans. At present, it serves data from eight permanent broadband and strong-motion networks from Greece and Cyprus, individual stations from the Balkans, temporary networks and aftershock deployments, and earthquake engineering experimental facilities. EIDA@NOA provides open and unlimited access from redundant node end points, intended mainly for research purposes (see Data and Resources). Analysis and quality control of the complete seismic data archive is performed initially by calculating waveform metrics and data availability. Seismic ambient noise metrics are estimated based on power spectral densities, and an assessment of each station’s statistical mode is achieved within each network and across networks. Moreover, the minimum ambient noise level expected for strong-motion installations is defined. Sensor orientation is estimated using surface-wave polarization methods to detect stations with misalignment on particular epochs. A single data center that hosts the complete seismic data archives with their respective metadata from networks covering similar geographical areas allows coordination between network operators and facilitates the adhesion to widely used best practices regarding station installation, data curation, and metadata definition. The overall achievement is harmonization among all contributing networks and a wider usage of all data archives, ultimately strengthening seismological research efforts in the region.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 58
    Publication Date: 2021-03-10
    Description: We present the first estimates of moment tensor solutions and probabilistic nonlinear hypocenter locations for the 1968 Belice earthquake sequence, which is the most relevant seismic activity occurred in western Sicily in historical times. This seismic phase, including six earthquakes with magnitude between 5 and 6.4, produced severe damages and fatalities in a sector of the Nubia–Eurasia plate margin, previously considered aseismic. Poorly constrained and often controversial hypocenter locations and focal mechanism solutions available from the literature have led to a long-lasting ambiguity on the possible causative source of the sequence, also having primary effects on the regional seismotectonic modeling and seismic hazard evaluations. The two main fault models proposed in the literature alternatively assigned a primary role to the north-northwest-trending Nubia–Eurasia plate convergence, causing thrust faulting on about east-trending structures or to the differential foreland retreat driving dextral strike-slip movements on about north–south-oriented faults. By focusing on the starting and the most energetic phase of the 1968 sequence, we computed the moment tensor solutions for three of the strongest earthquakes using digitized waveforms and a time-domain waveform inversion technique. Then, we also analyzed, by means of a Bayesian hypocenter location technique, the spatial distribution of the 1968 earthquakes. All the results indicate that the 1968 Belice sequence was characterized by predominant reverse faulting occurring on about east-to-northeast-trending structures, thus solving the dualism between models previously proposed in the literature. Our findings well agree both with the geodynamic framework governed by the Nubia–Eurasia north-northwest-trending convergence and with the geological reconstructions of the regional thrust front in the western Sicily area. The results of moment tensor estimations and nonlinear hypocenter locations furnishing an improved knowledge of the most relevant seismic activity of western Sicily also concur to better constraint the seismotectonic modeling of the region.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 59
    Publication Date: 2021-03-10
    Description: The 2016 Fukushima earthquake (M 7.4) generated a moderate tsunami, which was recorded by the offshore pressure gauges of the Seafloor Observation Network for Earthquakes and Tsunamis (S-net). We used 28 S-net pressure gauge records for tsunami data assimilation and forecasted the tsunami waveforms at four tide gauges on the Sanriku coast. The S-net raw records were processed using two different methods. In the first method, we removed the tidal components by polynomial fitting and applied a low-pass filter. In the second method, we used a real-time tsunami detection algorithm based on ensemble empirical mode decomposition to extract the tsunami signals, imitating real-time operations for tsunami early warning. The forecast accuracy scores of the two detection methods are 60% and 74%, respectively, for a time window of 35 min, but they improve to 89% and 94% if we neglect the stations with imperfect modeling or insufficient offshore observations. Hence, the tsunami data assimilation approach can be put into practice with the help of the real-time tsunami detection algorithm.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 60
    Publication Date: 2021-03-17
    Description: The NOANET network is a continuously operating high-rate Global Navigation Satellite System (GNSS) network in Greece whose primary role is to enhance and support geophysical research employing GNSS data. This network is operated by the Institute of Geodynamics of the National Observatory of Athens and currently (September 2020) consists of 26 stations, most of which are located close to major seismogenic structures of Greece to optimally measure tectonic and seismically induced motions. All NOANET receivers are configured to record and collect data with a sampling rate of 1 Hz, although some of them also collect data every 5 Hz on their ring buffer. The network is committed to free and open data sharing within the scientific community, and the collected data are made available via the NOANET data repository and distribution point for all interested parties with no limitations. Integrity, validity, and quality checks of the acquired data are performed using a variety of software tools along with in-house developed programs to supervise the network performance and detect ill-formed data and/or awkward station behavior. In addition, the conventional low-rate GNSS observation data of all NOANET stations are routinely processed on a daily basis to supervise their performance through their position time series. Since the beginning of its establishment, the NOANET network has recorded a variety of deformation signals, and a large number of published papers have used GNSS data from stations that are part of NOANET to constrain, model, and interpret the nature of the associated geophysical phenomena.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 61
    Publication Date: 2021-03-17
    Description: The National Seismological Network (NSN) of India has a history of more than 120 yr. During the last two decades, the NSN has gone through a significant modernization process, involving installation of seismic stations equipped with a broadband seismograph (BBS) and a strong-motion accelerograph (SMA). Each station has a very-small-aperture terminal connectivity for streaming data in real time to the central receiving station (CRS) in New Delhi. Seismic data recorded by the network are analyzed continuously on 24×7 basis to monitor the earthquakes in India and its adjoining regions. In this article, we present details of BBS and SMA network configurations; data streaming from the field seismic stations to the CRS for analysis; and the automatic and manual publication of the earthquake parameters including location coordinates, focal depth, time of occurrence, and magnitude, etc. Details of historically significant analog seismic charts and the seismic catalog, which includes more than 34,000 events with magnitude Mw 1.7–9.3 since 1505, are provided. The national network of India has been strengthened over the years and is now capable of estimating the main earthquake source parameters within ∼5–10min with an average of about 8.0 min. The spatial analysis of minimum magnitude of completeness further indicates a significant enhancement in minimum threshold magnitude detection capability of the network in recent decades.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 62
    Publication Date: 2021-03-16
    Description: The time-averaged shear-wave velocity in the upper 30 m depth from the ground surface, or VS30, is often used as a predictor to describe local site effects in ground-motion models. Although VS30 is typically determined from in situ measurements, it is not always feasible to obtain such measurements due to project restrictions or site accessibility. This motivates the development and use of proxy-based VS30 predictions that leverage more readily available secondary information such as surface geology, topographic slope, or geomorphic terrain classes to estimate the mean VS30 and associated uncertainty. Traditionally, empirical distributions of VS30 have been observed to have long right tails, leading to high levels of associated uncertainty. In this study, we present a physical framework that is grounded in fundamental principles of geostatistics and probability to explain the uncertainty and skewness associated with VS30 measurements. Specifically, by invoking Lyapunov’s central limit theorem, we hypothesize that the distribution of VS30 can be theoretically approximated by a reciprocal–normal distribution. We show that a non-normal and skewed distribution of VS30 is to be expected and is not a sign of measurement error or sampling bias, although sampling bias can exaggerate such skewness. Our framework also enables us to propose the mode as a characteristic value of VS30 measurements, as opposed to the mean or median, which can overestimate the most probable value.
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 63
    Publication Date: 2021-03-17
    Description: The intrinsic array nature of distributed acoustic sensing (DAS) makes it suitable for applying beamforming techniques commonly used in traditional seismometer arrays for enhancing weak and coherent seismic phases from distant seismic events. We test the capacity of a dark-fiber DAS array in the Sacramento basin, northern California, to detect small earthquakes at The Geysers geothermal field, at a distance of ∼100  km from the DAS array, using beamforming. We use a slowness range appropriate for ∼0.5–1.0  Hz surface waves that are well recorded by the DAS array. To take advantage of the large aperture, we divide the ∼20  km DAS cable into eight subarrays of aperture ∼1.5–2.0  km each, and apply beamforming independently to each subarray using phase-weighted stacking. The presence of subarrays of different orientations provides some sensitivity to back azimuth. We apply a short-term average/long-term average detector to the beam at each subarray. Simultaneous detections over multiple subarrays, evaluated using a voting scheme, are inferred to be caused by the same earthquake, whereas false detections caused by anthropogenic noise are expected to be localized to one or two subarrays. Analyzing 45 days of continuous DAS data, we were able to detect all earthquakes with M≥2.4, while missing most of the smaller magnitude earthquakes, with no false detections due to seismic noise. In comparison, a single broadband seismometer co-located with the DAS array was unable to detect any earthquake of M
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 64
    Publication Date: 2021-03-17
    Description: In anticipation of future spacecraft missions to icy ocean worlds, the Seismometer to Investigate Ice and Ocean Structure (SIIOS) was funded by National Aeronautics and Space Administration, to prepare for seismologic investigations of these worlds. During the summer of 2018, the SIIOS team deployed a seismic experiment on the Greenland ice sheet situated, approximately, 80 km north of Qaanaaq, Greenland. The seismometers deployed included one Trillium 120 s Posthole (TPH) broadband seismometer, 13 Silicon Audio flight-candidate seismometers, and five Sercel L28 4.5 Hz geophones. Seismometers were buried 1 m deep in the firn in a cross-shaped array centered on a collocated TPH and Silicon Audio instrument. One part of the array consisted of Silicon Audio and Sercel geophones situated 1 m from the center of the array in the ordinal directions. A second set of four Silicon Audio instruments was situated 1 km from the center of the array in the cardinal directions. A mock-lander spacecraft was placed at the array center and instrumented with four Silicon Audio seismometers. We performed an active-source experiment and a passive-listening experiment that lasted for, approximately, 12 days. The active–source experiment consisted of 9–12 sledgehammer strikes to an aluminum plate at 10 separate locations up to 100 m from the array center. The passive experiment recorded the ice-sheet ambient background noise, as well as local and regional events. Both datasets will be used to quantify differences in spacecraft instrumentation deployment strategies, and for evaluating science capabilities for single-station and small-aperture seismic arrays in future geophysical missions. Our initial results indicate that the flight-candidate seismometer performs comparably to the TPH at frequencies above 0.1 Hz and that instruments coupled to the mock-lander perform comparably to ground-based instrumentation in the frequency band of 0.1–10 Hz. For future icy ocean world missions, a deck-coupled seismometer would perform similarly to a ground-based deployment across the most frequency bands.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 65
    Publication Date: 2021-03-17
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 66
    Publication Date: 2021-03-16
    Description: We used our previously published Lg-wave spectral ratio model to develop a model of pseudospectral acceleration (PSA) response ratios at sites in the Atlantic and Gulf Coastal Plain, relative to a reference site condition defined as the mean response for site locations outside the Coastal Plain. The model is strongly dependent on sediment thickness. The results of this study can be used to predict PSA response, for linear behavior, at sites in the Atlantic or Gulf Coastal Plain with a known thickness of Coastal Plain sediment, given a ground-motion model for reference site conditions outside the Coastal Plain region of the central and eastern United States.
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 67
    Publication Date: 2021-03-16
    Description: In this article, we address the question of how observed ground-motion data can most effectively be modeled for engineering seismological purposes. Toward this goal, we use a data-driven method, based on a deep-learning autoencoder with a variable number of nodes in the bottleneck layer, to determine how many parameters are needed to reconstruct synthetic and observed ground-motion data in terms of their median values and scatter. The reconstruction error as a function of the number of nodes in the bottleneck is used as an indicator of the underlying dimensionality of ground-motion data, that is, the minimum number of predictor variables needed in a ground-motion model. Two synthetic and one observed datasets are studied to prove the performance of the proposed method. We find that mapping ground-motion data to a 2D manifold primarily captures magnitude and distance information and is suited for an approximate data reconstruction. The data reconstruction improves with an increasing number of bottleneck nodes of up to three and four, but it saturates if more nodes are added to the bottleneck.
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 68
    Publication Date: 2021-03-17
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 69
    Publication Date: 2021-03-17
    Description: The Swedish National Seismic Network (SNSN) was modernized and rapidly expanded during the period 1998–2012. The network currently operates 68 permanent seismic stations, all with broadband instruments supplying real-time continuous data at 100 samples per second. Continuous data from 10 stations are shared with the international community via Orfeus, and approximately 10 stations of their individual choice are shared with institutes in neighboring countries (Denmark, Finland, Norway, and Germany). The SNSN uses the South Iceland Lowland (SIL) system as the primary system for automatic detection and event definition. In addition, an in-house system based on migration and stacking is used for automatic detection of small events, and implementations of SeisComP (SC) and Earthworm are used primarily for rapid detection of larger regional events. Global monitoring is performed with SC, using approximately 250 global stations, and we operate a continuous rapid risk assessment system serving Swedish crisis management authorities. Since the start of automatic processing in August 2000, the SNSN has recorded and interactively analyzed more than 171,000 seismic events, of which 10,700 were earthquakes with local magnitudes ranging from around −1 to 4.3. The microearthquake activity detected in the last 20 yr has significantly improved the identification and understanding of seismically active structures in Sweden.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 70
    Publication Date: 2021-03-17
    Description: Risks assessment and risks comparison are basic concepts for emergency management. In the fields of earthquake engineering and engineering seismology, the operational earthquake loss forecasting (OELF) is the research frontier for the assessment of short-term seismic risk. It combines seismicity models, continuously updated based on ground-motion monitoring (i.e., operational earthquake forecasting), with large-scale vulnerability models for the built environment and exposure data. With the aim of contributing to the discussion about capabilities and limitations of OELF, the study presented aims at comparing the OELF results and the fatality risk (based on fatality data) related to coronavirus 2019 (COVID-19) that, at the time of writing, is perceived as very relevant and required unprecedented risk reduction measures in several countries, most notably Italy. Results show that, at a national scale in Italy, the COVID-19 risk has been higher than the seismic risk during the two pandemic waves even if, at the end of the so-called lockdown, the evolution of the pandemic suggested the possibility (not realized) of reaching a situation of comparable seismic and COVID-19 risks in a few weeks. Because the two risks vary at a local scale, risks comparison was also carried out on a regional basis, showing that, before the beginning of the second wave, in some cases, the seismic risk, as assessed by means of OELF, was larger than the pandemic one.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 71
    Publication Date: 2021-03-17
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 72
    Publication Date: 2021-03-09
    Description: Induced seismicity predominantly occurs along faults that are optimally oriented to the local principal compressive stress direction, and the characterization of these stress orientations is an important component of understanding seismic hazards. The seismicity rate in southern Kansas rapidly increased in 2013 primarily due to the disposal of large volumes of wastewater into the Arbuckle Group. Previously, local stress orientations in this area were poorly constrained, which limited our understanding of the complex faulting and diverse earthquake mechanisms in this region. We use shear-wave splitting and focal mechanism inversion techniques to create multiple, independent estimates of maximum horizontal stress directions (SHmax) and their spatial variation across the study area. We then create an integrated model of stress orientations for southern Kansas and northern Oklahoma using our local results in conjunction with previous, regional stress orientation estimates. We find that SHmax in both southern Kansas and central Oklahoma exhibits an east-northeast (∼N78° E) orientation, and these regions bound a northeast (∼N59° E) rotation within a ∼20  km area in northern Oklahoma near the Nemaha ridge.
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 73
    Publication Date: 2021-03-09
    Description: Ambient seismic noise (ASN) is becoming of interest for geophysical exploration and engineering seismology, because it is possible to exploit its potential for imaging. Theory asserts that the Green’s function can be retrieved from correlations within a diffuse field. Surface waves are the most conspicuous part of Green’s function in layered media. Thus, the velocities of surface waves can be obtained from ASN if the wavefield is diffuse. There is widespread interest in the conditions of emergence and properties of diffuse fields. In the applications, useful approximations of the Green’s function can be obtained from cross correlations of recorded motions of ASN. An elastic field is diffuse if the background illumination is azimuthally uniform and equipartitioned. It happens with the coda waves in earthquakes and has been verified in carefully planned experiments. For one of these data sets, the 1999 Chilpancingo (Mexico) experiment, there are some records of earthquake pre-events that undoubtedly are composed of ASN, so that the processing for coda can be tested on them. We decompose the ASN energies and study their equilibration. The scheme is inspired by the original experiment and uses the ASN recorded in an L-shaped array that allows the computation of spatial derivatives. It requires care in establishing the appropriate ranges for measuring parameters. In this search for robust indicators of diffusivity, we are led to establish that under certain circumstances, the S and P energy equilibration is a process that anticipates the diffusion regime (not necessarily isotropy), which justifies the use of horizontal-to-vertical spectral ratio in the context of diffuse-field theory.
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 74
    Publication Date: 2021-03-09
    Description: Traditional, empirical ground-motion models (GMMs) are developed by prescribing a functional form between predictive parameters and ground-motion intensity measures. Machine-learning techniques may serve as a fully data-driven alternative to widely used regression techniques, as they do not require explicitly defining these relationships. Although, machine-learning methods offer a nonparametric alternative to regression methods, there are few studies that develop and assess performance of traditional versus machine-learning GMMs side by side. We compare the performance and behavior of these two approaches: a mixed-effects maximum-likelihood (MEML) model and a feed-forward artificial neural network (ANN). We develop and train both models on the same dataset from southern California. We subsequently test both models on a dataset from the 2019 Ridgecrest sequence, in a new region and on magnitudes outside the range of the training dataset, to examine model portability. Our models estimate horizontal peak ground acceleration, and the input parameters include moment magnitude (M) and hypocentral distance (Rhyp), and some include a site parameter, either VS30 or κ0. We find that, with our small set of input parameters, the ANN generally shows more site-specific predictions than the MEML model with more variation between sites, and, performs better than their corresponding MEML model, when applied “blind” to our testing dataset (in which the MEML random effects cannot be considered). Although, previous studies have found that κ0 may be a better predictor of site effects than VS30, we found similar performance, suggesting that including a site parameter may be more important than the physical meaning of the parameter. Finally, when applying our models to our Ridgecrest dataset, we find that both methods perform well; however, the MEML models perform better with the new dataset than the ANN models, suggesting that future applications of ANN models may need to consider how to accommodate model portability.
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 75
    Publication Date: 2021-03-30
    Description: Recent efforts to characterize small (Mw
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 76
    Publication Date: 2021-03-30
    Description: The Global Seismographic Network (GSN) is a multiuse, globally distributed seismic network used by seismologists, to both characterize earthquakes and study the Earth’s interior. Most stations in the network have two collocated broadband seismometers, which enable network operators to identify potential metadata and sensor issues. In this study, we investigate the accuracy with which surface waves can be measured across the GSN, by comparing waveforms of vertical-component Rayleigh waves from Mw 6 and larger events between collocated sensor pairs. We calculate both the amplitude deviation and correlation coefficient between waveforms at sensor pairs. In total, we make measurements on over 670,000 event–station pairs from events that occurred from 1 January 2010 to 1 January 2020. We find that the average sensor-pair amplitude deviation, and, therefore, GSN calibration level, is, approximately, 4% in the 25–250 s period band. Although, we find little difference in sensor-pair amplitude deviations as a function of period across the entire network, the amount of useable data decreases rapidly as a function of increasing period. For instance, we determined that just over 12% of records at 250 s period provided useable recordings (e.g., sensor-pair amplitude deviations of less than 20% and sensor-pair correlation greater than 0.95). We then use these amplitude-estimate deviations to identify how data coverage and quality could be limiting our ability to invert for whole Earth 3D attenuation models. We find an increase in the variance of our attenuation models with increasing period. For example, our degree 12 attenuation inversion at 250 s period shows 32% more variance than our degree 12 attenuation model at 25 s. This indicates that discrepancies of deep-mantle tomography between studies could be the result of these large uncertainties. Because these high uncertainties arise from limited, high-quality observations of long-period (〉100  s) surface waves, improving data quality at remote GSN sites could greatly improve ray-path coverage, and facilitate more accurate and higher resolution models of deep Earth structure.
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 77
    Publication Date: 2021-02-10
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 78
    Publication Date: 2021-02-10
    Description: Modern seismic and Global Navigation Satellite Systems stations are nowadays equipped with Internet of Things devices that acquire, process, and transmit various geophysical parameters in near-real time. This technological advance has introduced a new threat paradigm for common seismological devices. Such threats can be assessed with standard information security methods and practices. This article aims to identify security weaknesses, describe weak security points and potential attacks on such environments, and anticipate the countermeasures needed. Real tests and attacks have been applied to demonstrate the lack of data encryption and user authentication processes, the risks posed by unencrypted communication protocols, unsafe practices regarding settings and passwords, and poor design implementations. All these factors may impact and possibly disrupt the daily operation of seismic observatories because they can lead to falsifying data, altering configurations, or producing malicious false alarms. These in turn may cause unnecessary public concern or distrust, financial losses, or even national security issues. For all these reasons, several countermeasures and solutions are also proposed and evaluated to address each of the identified vulnerabilities.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 79
    Publication Date: 2021-02-16
    Description: Empirical transfer functions (ETFs) between seismic records observed at the surface and depth represent a powerful tool to estimate site effects for earthquake hazard analysis. However, conventional modeling of site amplification, with assumptions of horizontally polarized shear waves propagating vertically through 1D layered homogeneous media, often poorly predicts the ETFs, particularly, in which large lateral variations of velocity are present. Here, we test whether more accurate site effects can be obtained from theoretical transfer functions (TTFs) extracted from physics-based simulations that naturally incorporate the complex material properties. We select two well-documented downhole sites (the KiK-net site TKCH05 in Japan and the Garner Valley site, Garner Valley Downhole Array, in southern California) for our study. The 3D subsurface geometry at the two sites is estimated by means of the surface topography near the sites and information from the shear-wave profiles obtained from borehole logs. By comparing the TTFs to ETFs at the selected sites, we show how simulations using the calibrated 3D models can significantly improve site amplification estimates as compared to 1D model predictions. The primary reason for this improvement in 3D models is redirection of scattering from vertically propagating to more realistic obliquely propagating waves, which alleviates artificial amplification at nodes in the vertical-incidence response of corresponding 1D approximations, resulting in improvement of site effect estimation. The results demonstrate the importance of reliable calibration of subsurface structure and material properties in site response studies.
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 80
    Publication Date: 2021-02-02
    Description: Based on the diffuse field concept for a horizontal-to-vertical spectral ratio of earthquakes (eHVSR), the effectiveness of eHVSRs to invert P- and S-wave velocity structures down to the seismological bedrock (with the S-wave velocity of 3  km/s or higher) has been shown in several published works. An empirical method to correct the difference between eHVSR and a horizontal-to-vertical ratio of microtremors (mHVSR), which is called earthquake-to-microtremor ratio (EMR), has also been proposed for strong-motion sites in Japan. However, the applicability of EMR outside of Japan may not be warranted. We test EMR applicability for the Grenoble basin in France with plentiful microtremor data together with observed weak-motion recordings at five sites. We thereby establish a systematic procedure to estimate the velocity structure from microtremors and delineate the fundamental characteristics of the velocity structures. We first calculate the EMR specific for the Grenoble basin (EMRG) and calculate pseudo eHVSR (pHVSR) from EMRG and mHVSR. We compare the pHVSRs with the eHVSRs at five sites and find sufficient similarity to each other. Then, we invert velocity structures from eHVSRs, pHVSRs, and mHVSRs. The velocity structures from eHVSRs are much closer to those from pHVSRs than those from mHVSRs. We need to introduce a number of layers with gradually increasing S-wave velocities below the geological basin boundary from a previous gravity study because the theoretical eHVSR of the model with a large velocity contrast has larger peak amplitudes than the observed. The depth of the S-wave velocity of 1.3  km/s (Z1.3) shows a strong, linear correlation with the geological boundary depth. Finally, we apply our validated methodology and invert velocity structures using pHVSRs at 14 sites where there are no observed earthquakes. The overall picture of Z1.3 at a cross section in the northeastern part of the basin corresponds to the geological boundary.
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 81
    Publication Date: 2021-02-23
    Description: The Wabash Valley seismic zone (WVSZ) is a region of diffuse, modern intraplate seismicity in the central United States with a history of strong, late Quaternary and Holocene seismicity as determined through paleoliquefaction studies. Yet, there are no specific faults linked to these strong WVSZ paleoearthquakes, some of which were as large as Mw 7.2–7.5. A multidisciplinary investigation of a linear, 5-kilometer-long and ∼3-meter-high scarp on the Ohio River floodplain in the southernmost WVSZ in western Kentucky evaluated whether the scarp is a fluvial landform or a tectonic feature. Geomorphic mapping and optically stimulated luminescence geochronology show that the age and orientation of the scarp are inconsistent with surrounding fluvial landforms. Trenching, core drilling, seismic reflection, electrical resistivity profiling, and cross sections of petroleum well logs all indicate a blind fault directly underlies the scarp. The scarp is interpreted to be the fold axis of a down-to-the-west monocline formed in alluvium by slip on the underlying blind fault, herein named the Uniontown fault. The Uniontown fault connects the Hovey Lake fault, striking N20°E and having   ∼0.5  km of documented strike-slip offset, with an unnamed fault complex to the south that strikes N40°E, suggesting the Uniontown fault is part of a larger, Paleozoic structure that has been reactivated with strike-slip deformation. Geomorphic mapping utilizing luminescence and radiocarbon geochronology indicates that folding and faulting occurred ∼3.5  ka. Paleoliquefaction was suppressed by a thick clay cap in the main Ohio Valley, but paleoliquefaction features are widespread on Ohio River tributaries. Gravel dikes at one site had a maximum age of 3.4±0.4  ka, confirming the region has experienced strong, late Holocene shaking. Estimates using vertical displacement and rupture length indicate that slip on the Uniontown scarp could produce an Mw 6.2–7.7 earthquake, which is comparable to other large paleoearthquakes in the WVSZ paleoseismic record.
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 82
    Publication Date: 2021-02-24
    Description: Infrasound sensors record the ambient acoustic field that contains not only signals of interest but also noise and clutter. Noise is defined as atmospheric turbulence that is incoherent over the distances of meters, whereas, signals of interest and clutter are acoustic pressure waves that are coherent over 10s of meters to 100s of kilometers. There is a growing interest in monitoring sources that extend across the acoustic spectrum from infrasound (below 20 Hz) into the low-end audible acoustic (20–1000 Hz). Monitoring of these extended band signals with a single sensor is made possible with the development of contemporary infrasound sensors, such as Hyperion IFS-3000 with a flat response from 0.01 to 1000 Hz. Combining infrasound sensors with seismometers provides opportunity to better assess noise contributions for both sensor types and improve characterization of sources that occur close to the solid earth–atmosphere boundary. Because sensors are installed to target these broadband acoustic sources, considerations need to be made when selecting a mechanical wind filter to mitigate the noise, while minimizing the impacts to the signals of interest across these frequency ranges. Motivated by these opportunities, this article compares traditional infrasound wind filter designs, that is, porous hoses rosettes and domes, in an urban setting for frequencies 0.01–45 Hz. Data analysis compares the filters, in terms of their response to noise and signals with direct comparisons of wind filters, as a function of frequency. The quantification of performance of these filters in an urban setting provides insight into their effects on detection of sources of interest in this environment.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 83
    Publication Date: 2021-02-10
    Description: The 1763 Komárom (now Komárno, Slovakia) earthquake was the biggest historical and recorded seismic event in the Pannonian basin in the last millennium. Contemporary documents including building damage reports from the towns of Komárom, Győr, and other locations, aftershock reports, liquefaction reports, and local surveys of taxpayers’ loss in villages are presented to build all available datasets to evaluate the effects of the quake. Damage in towns was recorded in detail only at the aforementioned two locations. However, the spatial distribution of shaking intensity can be inferred over a larger area from information obtained from other sources. Distributions of secondary environmental effects, aftershocks, and village damage show interesting spatial correlation among them and indicate that the epicenter may have been west of Komárom, north of the Danube River.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 84
    Publication Date: 2021-04-06
    Description: The use of numerical simulations in probabilistic seismic hazard analysis (PSHA) has achieved a promising level of reliability in recent years. One example is the CyberShake project, which incorporates physics-based 3D ground-motion simulations within seismic hazard calculations. Nonetheless, considerable computational time and resources are required due to the significant processing requirements imposed by source-based models on one hand, and the large number of seismic sources and possible rupture variations on the other. This article proposes to use a less computationally demanding simulation-based PSHA framework for CyberShake. The framework can accurately represent the seismic hazard at a site, by only considering a subset of all the possible earthquake scenarios, based on a Monte-Carlo simulation procedure that generates earthquake catalogs having a specified duration. In this case, ground motions need only be simulated for the scenarios selected in the earthquake catalog, and hazard calculations are limited to this subset of scenarios. To validate the method and evaluate its accuracy in the CyberShake platform, the proposed framework is applied to three sites in southern California, and hazard calculations are performed for earthquake catalogs with different lengths. The resulting hazard curves are then benchmarked against those obtained by considering the entire set of earthquake scenarios and simulations, as done in CyberShake. Both approaches yield similar estimates of the hazard curves for elastic pseudospectral accelerations and inelastic demands, with errors that depend on the length of the Monte-Carlo catalog. With 200,000 yr catalogs, the errors are consistently smaller than 5% at the 2% probability of exceedance in 50 yr hazard level, using only ∼3% of the entire set of simulations. Both approaches also produce similar disaggregation patterns. The results demonstrate the potential of the proposed approach in a simulation-based PSHA platform like CyberShake and as a ground-motion selection tool for seismic demand analyses.
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 85
    Publication Date: 2021-04-13
    Description: This study addresses questions about the productivity of Cascadia mainshock–aftershock sequences using earthquake catalogs produced by the Geological Survey of Canada and the Pacific Northwest Seismic Network. Questions concern the likelihood that future moderate to large intermediate depth intraslab earthquakes in Cascadia would have as few detectable aftershocks as those documented since 1949. More broadly, for Cascadia, we consider if aftershock productivities vary spatially, if they are outliers among global subduction zones, and if they are consistent with a physical model in which aftershocks are clock-advanced versions of tectonically driven background seismicity. A practical motivation for this study is to assess the likely accuracy of aftershock forecasts based on productivities derived from global data that are now being issued routinely by the U.S. Geological Survey. For this reason, we estimated productivity following the identical procedures used in those forecasts and described in Page et al. (2016). Results indicate that in Cascadia we can say that the next intermediate depth intraslab earthquake will likely have just a few detectable aftershocks and that aftershock productivity appears to be an outlier among global subduction zones, with rates that on average are lower by more than half, except for mainshocks in the upper plate. Our results are consistent with a clock-advance model; productivities may be related to the proximity of mainshocks to a population of seismogenic fault patches and correlate with background seismicity rates. The latter and a clear correlation between productivities with mainshock depth indicate that both factors may have predictive value for aftershock forecasting.
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 86
    Publication Date: 2021-04-14
    Description: We present a detailed analysis of the 10 June 2019 ML 4.0 Ohio earthquake sequence, which is the largest earthquake that struck Lake County, northeastern Ohio, since 1986. This sequence is well recorded by local seismic networks, which provides an unprecedented opportunity to understand the local seismotectonics. We utilize a waveform-based cross-correlation method to identify ∼12 times more events than reported by the Advanced National Seismic System (ANSS) Comprehensive Earthquake Catalog: the whole sequence started with several small earthquakes (ML 1–2) beginning 12 March 2019, and the last one occurred ∼1min immediately before the ML 4.0 mainshock; many previously unreported aftershocks (ML 0.3–2.2) are found, which were active for the first week after the mainshock; another major sequence with a 7 December 2019 ML 2.6 mainshock occurred and also started with a few smaller events beginning in mid-November and was followed by its own aftershocks. The relocated seismicity delineates a linear feature, orientation of which is consistent with the resolved focal plane that may correspond to the ruptured fault. Our results highlight that closer monitoring of local seismicity is crucial for understanding the seismotectonics and mitigating future seismic hazard around the southern Great Lakes.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 87
    Publication Date: 2021-04-14
    Description: Seismic bulletins, with trustworthy phase picks, origin times, and source locations are key for regional seismic studies, such as travel-time (TT) tomography, attenuation tomography, and anisotropy studies. To lay the groundwork for such studies in Israel, we revised the seismic bulletin of Israel and the surrounding area and obtained a trustworthy TT data set. From the earthquake and explosion bulletins of the Geophysical Institute of Israel, we compiled a starting data set of about 123,000 earthquakes and explosions that occurred during the past 40 yr. After screening out the poorly recorded events, we were left with a data set of ∼38,000 well-recorded events. We then revised the remaining data set in two consecutive steps. In the first, we reviewed and updated station metadata, including changes in station metadata parameters over time. In the second step, we jointly relocated a list of selected seismic events, using the Bayesian hierarchical location software package (BayesLoc) of Myers et al. (2007) that performs joint relocation of multiple events. We observed striking dissimilarities between the spatial distributions of the newly relocated catalog and the initial locations. Although the depth distribution of the starting catalog is trimodal with peaks at 0, 5, and 10 km, the distribution in this study is unimodal, with a broad peak between 7.5 and 12.5 km. By differencing the observed arrival times and the origin times obtained through relocation with BayesLoc, we obtained a revised TT database that consists of 261,336 Pg, 132,876 Pn, 114,816 Sg, and 60,394 Sn arrivals, from a set of 30,458 jointly relocated seismic sources. We compared prerevision and postrevision TTs as a function of epicentral distance and concluded that the revised data set contains far fewer outliers and inconsistencies than the original data set. The revised TT data set may be used for seismic studies, such as TT tomography, attenuation tomography, and anisotropy studies.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 88
    Publication Date: 2021-04-14
    Description: Coulomb stress change is the change in resultant force of shear stress and friction imposed on a receiver fault plane. The resulting stress change is often computed using the Coulomb 3.4 and the postseismic Green’s functions and postseismic components (PSGRN-PSCMP) programs. Notwithstanding both preferences, both have incomplete optimally oriented failure planes (OOPs) and are inconvenient to resolve Coulomb stress changes on various fault planes placed in varying depths. Here, we present an alternative program termed AutoCoulomb. It leverages the shell command-line tool to automatically batch-process Coulomb stress changes on all sorts of receiver fault planes. We first validate the program. We then apply it to the 2020 Mw 7.8 Simeonof Island, Alaska, earthquake, as a case study. Our results show that Coulomb stress changes resolved on fixed receiver faults, using the three programs, are in line with each other. So are those resolved on 3D OOPs using the PSGRN–PSCMP and the AutoCoulomb programs. Nevertheless, Coulomb stress changes on 2D OOPs, generated by the AutoCoulomb program, always outweigh those done by the Coulomb 3.4 program, indicating that 2D OOPs constrained by the latter are not the most optimal. Some nonoptimal 2D OOPs result in the reversal of the signs of Coulomb stress changes, posing a risk of misleading stress shadows with negative Coulomb stress changes. For the case study, the 28 July 2020 Mw 6.1 aftershock received a positive coseismic Coulomb stress change of ∼3.5 bars. In contrast, the compounded coseismic Coulomb stress changes at the hypocenters of the 1946 Mw 8.2, the 1948 Mw 7.2, and the 2020 Mw 7.8 earthquakes are within a range from −1.1 to 0.1 bar, suggesting that coseismic Coulomb stress changes promoted by preceding mainshocks alone are not responsible for these mainshocks. Other factors, such as postseismic viscoelastic relaxation, afterslip, and slow slip, may contribute to promoting their occurrence.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 89
    Publication Date: 2021-04-27
    Description: The structural complexity of active faults and the stress release history along the fault system may exert control on the locus and extent of individual earthquake ruptures. Fault bends, in particular, are often invoked as a possible mechanism for terminating earthquake ruptures. However, there are few records available to examine how these factors may influence the along-fault recurrence of earthquakes. We present a new paleoearthquake chronology for the southern San Andreas fault at Elizabeth Lake and integrate this record with existing paleoearthquake records to examine how the timing and frequency of earthquakes vary through a major restraining bend. This restraining bend features a mature, throughgoing right-lateral strike-slip fault, two major fault intersections, proposed subsurface fault dip changes, and a 〉200  km long section of fault misaligned with the regional plate motion. The Frazier Mountain, Elizabeth Lake, Pallett Creek, Wrightwood, and Pitman Canyon paleoseismic sites are located on this relatively linear surface trace of the San Andreas fault between fault bends. Our paleoseismic investigations at Elizabeth Lake document 4–5 earthquakes, since ∼1100  C.E., similar to the number of earthquakes recorded at Pallett Creek. In contrast, the Frazier Mountain and Wrightwood sites each record 8–9 earthquakes during this same time period. Differences in earthquake frequency demonstrate that fewer earthquakes rupture the central portion of the restraining bend than occur near the fault bends and intersections. Furthermore, the similarity of earthquake records from the Bidart Fan paleoseismic site northwest of the restraining bend and the Frazier Mountain paleoseismic site suggests that the broad, 30° curve of the Big Bend section of the San Andreas fault exerts less influence on fault rupture behavior than the 3D geometry of the Mojave sections of the fault.
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 90
    Publication Date: 2021-04-27
    Description: The NGA-Sub (subduction zone earthquake) database developed by the Pacific Earthquake Engineering Research center is used to derive new correlation coefficients for a number of ground-motion intensity measure (IM) parameters from ground motions in subduction zone earthquakes, considering both interface and intraslab tectonic settings. The IMs include peak ground acceleration, pseudospectral accelerations with periods from 0.01 to 10 s, Arias intensity, cumulative absolute velocity, peak ground velocity, and significant duration. Comparisons of the estimated correlation coefficients for ground motions from the interface and intraslab tectonic settings generally show a good agreement. Our estimations are also in good agreement with correlation coefficients estimated in previous studies that used ground motions from shallow crustal earthquakes, supporting the concept that any variation in correlation coefficients comes from spectral shape (i.e., the distribution of peaks and troughs) rather than tectonic region. This study also explores the influence of parameters such as magnitude, distance, and site conditions on the estimated correlation coefficients. We did not find apparent trends of the correlation coefficients with respect to these parameters. Finally, we propose analytical models to estimate correlation coefficients between the IMs explored in this study, considering both subduction interface and intraslab tectonic settings.
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 91
    Publication Date: 2021-04-27
    Description: Probabilistic seismic hazard assessment relies on long-term earthquake forecasts and ground-motion models. Our aim is to improve earthquake forecasts by including information derived from geodetic measurements, with an application to the Colombia–Ecuador megathrust. The annual rate of moment deficit accumulation at the interface is quantified from geodetically based interseismic coupling models. We look for Gutenberg–Richter recurrence models that match both past seismicity rates and the geodetic moment deficit rate, by adjusting the maximum magnitude. We explore the uncertainties on the seismic rates (a- and b-values, shape close to Mmax) and on the geodetic moment deficit rate to be released seismically. A distribution for the maximum magnitude Mmax bounding a series of earthquake recurrence models is obtained for the Colombia–Ecuador megathrust. Models associated with Mmax values compatible with the extension of the interface segment are selected. We show that the uncertainties mostly influencing the moment-balanced recurrence model are the fraction of geodetic moment released through aseismic processes and the form of the Gutenberg–Richter model close to Mmax. We combine the computed moment-balanced recurrence models with a ground-motion model, to obtain a series of uniform hazard spectra representative of uncertainties at one site on the coast. Considering the recent availability of a massive quantity of geodetic data, our approach could be used in other well-instrumented regions of the world.
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 92
    Publication Date: 2021-04-27
    Description: Ground-motion models (GMMs) are developed for peak ground displacement (PGD) and for bandlimited PGD based on strong-motion data that has been filtered as part of standard processing and the total PGD that includes the tectonic deformation as well as the vibratory ground motion. For the bandlimited PGD, we develop conditional ground-motion models (CGMMs) using subsets of the Pacific Earthquake Engineering Research Center Next Generation Attenuation-West2 Project (NGA-W2) database and the National Center for Research on Earthquake Engineering Taiwan Senior Seismic Hazard Analysis Committee level 3 project database. The CGMM approach includes the observed pseudospectral acceleration (PSA(T)) as an input parameter in addition to magnitude and distance. The period of the PSA(T) is used as an input parameter; it is magnitude dependent and is based on the period for which there is the highest correlation between the ln(PGD) and ln(PSA(T)). Two CGMMs are developed: a global model based on the NGA-W2 data and a region-specific model for Taiwan. The conditional PGD models are combined with traditional GMMs for PSA(T) values to develop GMMs for both the median and standard deviation of PGD without the dependence on PSA. A second set of PGD GMMs are developed to correct for two factors: the effect of the high-pass filtering from standard record processing and the stronger large magnitude (M〉6.5) scaling due to tectonic deformation. For magnitudes greater than 7, the PGD values from the total PGD GMMs are 2–5 times larger than the bandlimited PGD values based on the strong-motion data sets, but the increase is at very long periods. The appropriate PGD model to use, bandlimited PGD or total PGD, depends on the period range of interest for the specific engineering application.
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 93
    Publication Date: 2021-04-27
    Description: We report tremor or local earthquake signals that occurred during the propagation of Love and Rayleigh waves from the 2012 Mw 8.6 Sumatra earthquake in three intraplate regions: Yellowstone, central Utah, and Raton basin (Colorado). These surface waves likely also dynamically triggered seismic activity along the western boundary of the North American plate, and did not trigger seismic activity in the central and eastern United States. We report additional potential dynamic triggering in the three aforementioned intraplate regions by surface waves from 37 additional large earthquakes, recorded between 2004 and 2017. These surface waves’ transient stresses generally appear to trigger tremor in seismically, volcanically, and hydrothermally active regions, such as Yellowstone, if the waves also arrive from favorable directions. These stresses do not appear to be decisive factors for triggering local earthquakes reported for the Raton basin and central Utah, whereas, surface waves’ incidence angles do appear to be important there.
    Print ISSN: 0037-1106
    Electronic ISSN: 1943-3573
    Topics: Geosciences , Physics
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 94
    Publication Date: 2021-03-03
    Description: The aim of this article is to analyze the background, current status, and outlook of seismic monitoring products and services in Bulgaria, Moldova, Romania, and Ukraine. These countries manage seismic networks that contribute to the European Integrated Data Archive node in the framework of the Observatories and Research Facilities for European Seismology, which represents a collaborative effort in coordinating observational seismology across the European region through the collection, archiving, and dissemination of seismic waveform data, metadata, and related products. All of the aforementioned countries share a common threat: strong earthquakes occurring in the Vrancea area located in central-eastern Romania at intermediate depths (usually in the 60–180 km interval). Events such as the ones on 10 November 1940 and 4 March 1977 generated high damage in Romania, northern Bulgaria, and Moldova. In addition to Vrancea, crustal earthquakes in areas such as Shabla or Dulovo can lead to cross-border damage. Therefore, understanding the way national seismic networks are distributed, how they cooperate, and the products and services that they provide in (near) real time and their terms is of significant interest in the context of necessary hazard harmonization and joint emergency intervention and risk mitigation actions.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 95
    Publication Date: 2021-03-03
    Description: In recent years, the French seismological, geodetic, and gravimetric community has been structured within Réseau Sismologique et géodésique Français (RESIF) (French seismological and geodetic network). In addition to instrumental developments, RESIF has structured the work on French seismicity (metropolitan and overseas) within the RESIF transverse seismicity action (ATS). The purpose of this article is to present the ATS and the way it is structured to propose to the community different products: seismicity bulletin and catalog, historical and instrumental macroseismicity data, and ShakeMaps. The places where these products can be found are indicated, as well as the way they are realized and the improvements in progress for a better realization and availability. The link with European plate observing system is also underlined.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 96
    Publication Date: 2021-03-03
    Description: The Istituto Nazionale di Geofisica e Vulcanologia (INGV) is an Italian research institution with focus on earth sciences. Moreover, the INGV is the operational center for seismic surveillance and earthquake monitoring in Italy and is a part of the civil protection system as a center of expertise on seismic, volcanic, and tsunami risks.INGV operates the Italian National Seismic Network and other networks at national scale and is a primary node of the European Integrated Data Archive for archiving and distributing strong-motion and weak-motion seismic recordings. In the control room in Rome, INGV staff performs seismic surveillance and tsunami warning services; in Catania and Naples, the control rooms are devoted to volcanic surveillance. Volcano monitoring includes locating earthquakes in the regions around the Sicilian (Etna, Eolian Islands, and Pantelleria) and the Campanian (Vesuvius, Campi Fregrei, and Ischia) active volcanoes. The tsunami warning is based on earthquake location and magnitude (M) evaluation for moderate to large events in the Mediterranean region and also around the world. The technologists of the institute tuned the data acquisition system to accomplish, in near real time, automatic earthquake detection, hypocenter and magnitude determination, and evaluation of several seismological products (e.g., moment tensors and ShakeMaps). Database archiving of all parametric results is closely linked to the existing procedures of the INGV seismic surveillance environment and surveillance procedures. Earthquake information is routinely revised by the analysts of the Italian seismic bulletin. INGV provides earthquake information to the Department of Civil Protection (Dipartimento di Protezione Civile) to the scientific community and to the public through the web and social media. We aim at illustrating different aspects of earthquake monitoring at INGV: (1) network operations; (2) organizational structure and the hardware and software used; and (3) communication, including recent developments and planned improvements.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 97
    Publication Date: 2021-03-10
    Description: In this article, we describe the infrastructure developed and managed by the Italian National Institute of Oceanography and Applied Geophysics – OGS for the seismological and geodetic monitoring of northeastern Italy. The infrastructure was constituted in response to the ML 6.4 Friuli destructive earthquake of 1976, with the main mandate of supporting civil protection emergency activities. The OGS monitoring infrastructure is presently composed of a seismometric and a strong-motion network, complemented by a number of Global Navigation Satellite Systems stations, each delivering observational data in real time, which are collected and processed by the headquarters of the Center for Seismological Research of OGS in Udine. The OGS networks operate in close cooperation with Italian and international networks from neighboring countries, within the framework of the agreements for real-time data exchange, to obtain improved rapid earthquake location and magnitude estimations. Information regarding seismic events is released to the public through a dedicated web portal and, since 2013, through social media. Aside from the standard monitoring activities (〉30,000 events have been recorded since 1976), the OGS has progressively increased the number of services to the public and to the Civil Protection of the Friuli Venezia Giulia and Veneto regions. The high availability of good quality data has resulted in the enhancement of scientific products, including advanced seismological studies of the area, spanning broadly from seismic source characterization to engineering seismology. In the future, the OGS networks are expected to further contribute to the development of seismological research and monitoring infrastructures of the Central European region.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 98
    Publication Date: 2021-03-10
    Description: The Orfeus European Integrated Data Archive (EIDA) provides a federated approach to the dissemination of seismological waveform data and ensures access to 12 regional seismological data centers—the EIDA nodes. The Istituto Nazionale di Geofisica e Vulcanologia (INGV) is one of the founding partners of this EIDA federation and manages the EIDA data distribution node in Italy. INGV has actively managed the smaller MedNet archive since 1990 and adopted a more comprehensive and systematic approach to seismological data archiving since 2007. The Italian EIDA node data archive currently totals 90 TBytes of waveform data available for download, originating from 25 networks and 974 stations, provided by INGV, MedNet, or contributed by various partner institutions. Geographically, it covers mainly Italy and some stations from the Mediterranean region. The archive is currently growing at a rate of approximately 11 TB/yr. INGV recently strengthened its data management capabilities, resources, and infrastructure to effectively respond to the growing scale of station inventory, archive, and volumes of delivered data, and to acknowledge increasing attention toward open data sharing, appropriate attribution, and FAIR principles (Findability, Accessibility, Interoperability, and Reuse), as well as higher demands on data quality and expectations of the scientific user community. To this end, it established a dedicated internal unit in charge of all relevant activities related to the Italian EIDA node. In this article, we address key aspects of the EIDA node in Italy such as evolution and status of the seismological waveform archive, and we describe the technical, organizational, and operational setup of data and service management. We also outline ongoing activities and future evolutions aiming to further increase the quality of services, data availability, data and metadata quality, resilience, and sustainability.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 99
    Publication Date: 2021-03-10
    Description: In recent years, the Delaware basin of west Texas has seen a sharp rise in earthquake occurrence, driven in large part by increases in unconventional hydrocarbon production. The advent of Texas Seismological Network in 2017 has allowed for the characterization of these events in greater detail. We exploit the recent densification in seismic station coverage to study the spectral properties of earthquakes in this region. We show that the low-frequency moments of S-wave spectra, when corrected for site and distance, can be used to calibrate a consistent moment magnitude scale for small and moderate earthquakes. For a subset of 〉3000 well-recorded events, we compute earthquake stress drop under the assumption of Brune spectral model. Earthquakes in the Delaware basin show coherent spatial patterns in stress drop across the region. Through a reanalysis of independently collated data from the oil and gas industry, we find that earthquakes that are likely associated with hydraulic-fracturing operations have slightly different spectral characteristics than earthquakes that are likely associated with saltwater disposal. In particular, events associated with hydraulic fracturing show greater variability in the statistical distribution of stress drop and have higher median stress-drop values. Although the differences are subtle, they suggest that there may be important distinctions in the underlying physical mechanisms and resulting hazards of distinct classes of induced events, differences that may be unraveled with more detailed joint analyses of industrial and seismic datasets.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 100
    Publication Date: 2021-03-10
    Description: We image the shallow structure across the East Bench segment of the Wasatch fault system in Salt Lake City using ambient noise recorded by a month-long temporary linear seismic array of 32 stations. We first extract Rayleigh-wave signals between 0.4 and 1.1 s period using noise cross correlation. We then apply double beamforming to enhance coherent cross-correlation signals and at the same time measure frequency-dependent phase velocities across the array. For each location, based on available dispersion measurements, we perform an uncertainty-weighted least-squares inversion to obtain a 1D VS model from the surface to 400 m depth. We put all piece-wise continuous 1D models together to construct the final 2D VS model. The model reveals high velocities to the east of the Pleistocene Lake Bonneville shoreline reflecting thinner sediments and low velocities particularly in the top 200 m to the west corresponding to the Salt Lake basin sediments. In addition, there is an ∼400-m-wide low-velocity zone that narrows with depth adjacent to the surface trace of the East Bench fault, which we interpret as a fault-related damage zone. The damage zone is asymmetric, wider on the hanging wall (western) side and with greater velocity reduction. These results provide important constraints on normal-fault earthquake mechanics, Wasatch fault earthquake behavior, and urban seismic hazard in Salt Lake City.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
Close ⊗
This website uses cookies and the analysis tool Matomo. More information can be found here...