Walter, William R.

Robustness of source parameter estimates is a fundamental issue in understanding the relationships between small and large events; however, it is difficult to assess how much of the variability of the source parameters can be attributed to the physical source characteristics or to the uncertainties of the methods and data used to estimate the values. In this study, we apply the coda method by Mayeda et al. using the coda calibration tool (CCT), a freely available Java-based code (https://github.com/LLNL/coda-calibration-tool) to obtain a regional calibration for Central Italy for estimating stable source parameters. We demonstrate the power of the coda technique in this region and show that it provides the same robustness in source parameter estimation as a data-driven methodology [generalized inversion technique (GIT)], but with much fewer calibration events and stations. The Central Italy region is ideal for both GIT and coda approaches as it is characterized by high-quality data, including recent well-recorded seismic sequences such as L'Aquila (2009) and Amatrice–Norcia–Visso (2016–2017). This allows us to apply data-driven methods such as GIT and coda-based methods that require few, but high-quality data. The data set for GIT analysis includes ∼5000 earthquakes and more than 600 stations, while for coda analysis we used a small subset of 39 events spanning 3.5 < Mw < 6.33 and 14 well-distributed broad-band stations. For the common calibration events, as well as an additional 247 events (∼1.7 < Mw < ∼5.0) not used in either calibration, we find excellent agreement between GIT-derived and CCT-derived source spectra. This confirms the ability of the coda approach to obtain stable source parameters even with few calibration events and stations. Even reducing the coda calibration data set by 75 per cent, we found no appreciable degradation in performance. This validation of the coda calibration approach over a broad range of event size demonstrates that this procedure, once extended to other regions, represents a powerful tool for future routine applications to homogeneously evaluate robust source parameters on a national scale. Furthermore, the coda calibration procedure can homogenize the Mw estimates for small and large events without the necessity of introducing any conversion scale between narrow-band measures such as local magnitude (ML) and Mw, which has been shown to introduce significant bias.

The Amatrice–Norcia–Visso sequence is characterized by complex behaviour that is somewhat atypical of main-shock–aftershock sequences, as there were multiple large main shocks that continued for months. In this study we focus on the Amatrice sequence (main shock 2016 August 24, Mw = 5.97) to evaluate the apparent stress values and magnitude-dependent scaling in order to improve our knowledge of processes that control small and large earthquakes within this active region of Italy. Apparent stress is proportional to the ratio of radiated seismic energy and seismic moment, and as such, these stress parameters play an important role in hazard prediction as they have a strong effect on the observed and predicted ground shaking. We analyse 83 events of the sequence from 2016 August 24 to October 16, within a radius of 20 km from the main shock and with an Mw ranging between 5.97 and 2.72. Taking advantage of the averaging nature of coda waves, we analyse coda-envelope-based spectral ratios between neighbouring event pairs.We use equations proposed byWalter et al. to consider stable, low-frequency and high-frequency spectral ratio levelswhich provide measures of the corner frequency and apparent stress ratios of the events within the sequence. The results demonstrate non-self-similar behaviour within the sequence, suggesting a change in dynamics between the largest events and the smaller aftershocks. The apparent stress and corner frequency estimates are compared to those obtained by Malagnini and Munaf`o who utilized hundreds of direct S-wave spectral ratio measurements to obtain their results. Although our analysis is based only on 83 events, our results are in very good agreement, demonstrating once more that the use of coda waves is very stable and provides lower variance measures than those using direct waves. A comparison with recent Central Apennines source scaling models derived from various seismic sequences (1997–1998 Colfiorito, 2002 San Giuliano di Puglia, 2009 L’Aquila) shows that the Amatrice sequence source scaling in this study is well represented by the models proposed by Pacor et al. and Malagnini and Mayeda.

A stable estimate of the earthquake source spectra in the western Alps is obtained using an empirical method based on coda envelope amplitude measurements described by Mayeda for events ranging between MW∼ 1.0 and ∼5.0. Path corrections for consecutive narrow frequency bands ranging between 0.3 and 25.0 Hz were included using a simple 1-D model for five three-component stations of the Regional Seismic network of Northwestern Italy (RSNI). The 1-D assumption performs well, even though the region is characterized by a complex structural setting involving strong lateral variations in the Moho depth. For frequencies less than 1.0 Hz, we tied our dimensionless, distance-corrected coda amplitudes to an absolute scale in units of dyne cm by using independent moment magnitudes from long-period waveform modelling for three moderate magnitude events in the region. For the higher frequencies, we used small events as empirical Green's functions, with corner frequencies above 25.0 Hz. For each station, the procedure yields frequency-dependent corrections that account for site effects, including those related to fmax, as well as to S-to-coda transfer function effects. After the calibration was completed, the corrections were applied to the entire data set composed of 957 events. Our findings using the coda-derived source spectra are summarized as follows: (i) we derived stable estimates of seismic moment, M0, (and hence MW) as well as radiated S-wave energy, (ES), from waveforms recorded by as few as one station, for events that were too small to be waveform modelled (i.e. events less than MW∼ 3.5); (ii) the source spectra were used to derive an equivalent local magnitude, ML(coda), that is in excellent agreement with the network averaged values using direct S waves; (iii) scaled energy, graphic, where ER, the radiated seismic energy, is comparable to results from other tectonically active regions (e.g. western USA, Japan) and supports the idea that there is a fundamental difference in rupture dynamics between small and large crustal earthquakes in tectonically active regions.