Moment magnitudes (MW) are computed for small and moderate earthquakes using a spectral fitting method. 40 of the resulting values are compared with those from broadband moment tensor solutions and found to match with negligible offset and scatter for available MW values of between 2.8 and 5.0. Using the presented method, MW are computed for 679 earthquakes in Switzerland with a minimum ML= 1.3. A combined bootstrap and orthogonal L1 minimization is then used to produce a scaling relation between ML and MW. The scaling relation has a polynomial form and is shown to reduce the dependence of the predicted MW residual on magnitude relative to an existing linear scaling relation. The computation of MW using the presented spectral technique is fully automated at the Swiss Seismological Service, providing real-time solutions within 10 minutes of an event through a web-based XML database. The scaling between ML and MW is explored using synthetic data computed with a stochastic simulation method. It is shown that the scaling relation can be explained by the interaction of attenuation, the stress-drop and the Wood–Anderson filter. For instance, it is shown that the stress-drop controls the saturation of the ML scale, with low-stress drops (e.g. 0.1–1.0 MPa) leading to saturation at magnitudes as low as ML= 4.
We estimate moment magnitudes MW for earthquakes in Switzerland recorded between 1998 and 2009. Compared to previous studies the MW determination in Switzerland could be extended to lower magnitudes and scaling relations can be investigated from ML = 0 to 5.4. Three different spectral methods are applied to estimate MW and the scaling relations resulting from the different methods are compared with one another. Above ML = 4, the obtained estimates are consistent with the previously used scaling relation of MW = ML−0.3. BelowML = 4, all three methods indicate that a 1:1-type relationship is inappropriate. Therefore, w propose a new empirical piecewise ML to MW scaling relation for earthquakes in Switzerland. The scali ng is linear belowML = 2 and aboveML = 4. To obtain a smooth transition between the two linear scales we fit a quadratic relation in between ( 2 ≤ ML < 4). This scaling relation is also consistent with MW estimates from moment tensor solutions based on broadband waveform fit ting of local earthquakes with ML > 3.0.
We investigate the global variation of earthquake stress drops using spectra of about 2000 events of mb ≥ 5.5 between 1990 and 2007. We use an iterative least squares method to isolate source displacement spectra from travel path and receiver contributions, based on a convolutional model. The observed P wave source spectra are corrected with a globally averaged empirical correction spectrum and estimates of near‐source attenuation. Assuming a Brune‐type source model, we estimate corner frequencies and compute stress drops. Stress drop estimates for individual earthquakes range from about 0.3 to 50 MPa, but the median stress drop of about 4 MPa does not vary with moment, implying earthquake self‐similarity over the Mw = 5.2 to 8.3 range of our data. A comparison of our results with previous studies confirms this observation over most of the instrumentally observable magnitude range. While the absolute values of our estimated stress drops depend upon the assumed source model, we identify relative regional variations of stress drop that are robust with respect to the processing parameters and modeling assumptions, which includes an inherent assumption of constant rupture velocity. We find a dependence of median stress drop on focal mechanism, with a factor of 3–5 times higher stress drops for strike‐slip earthquakes and also find a factor of 2 times higher stress drops for intraplate earthquakes compared to interplate earthquakes.
By using seismic records of the 2004 magnitude 6.0 Parkfield earthquake, we identified a burst of high-frequency seismic radiation that occurred about 13 kilometers northwest of the hypocenter and 5 seconds after rupture initiation. We imaged this event in three dimensions by using a waveform back-projection method, as well as by timing distinct arrivals visible on many of the seismograms. The high-frequency event is located near the south edge of a large slip patch seen in most seismic and geodetic inversions, indicating that slip may have grown abruptly at this point. The time history obtained from full-waveform back projection suggests a rupture velocity of 2.5 kilometers per second. Energy estimates for the subevent, together with long-period slip inversions, indicate a lower average stress drop for the northern part of the Parkfield earthquake compared with that for the region near its hypocenter, which is in agreement with stress-drop estimates obtained from small-magnitude aftershocks.
We estimate source parameters from spectra of 42367 earthquakes between 1984 and 2005 that occurred in the Parkfield segment of the San Andreas Fault in central California. We use a method that isolates the source term of the displacement spectra based on a convolutional model and correct the observed P wave source spectra with a spatially varying empirical Green's function (EGF). Our Brune‐type stress drop estimates vary from 0.1 to over 100 MPa with a median value of 6.75 MPa, which is nearly constant with moment, implying self‐similarity over the ML = 0.5 to 3.0 range of our data. The corner frequency decreases for earthquakes at shallower depths, consistent with slower rupture velocities and reduced shear wave velocities in local velocity models. The estimated median stress drops show significant lateral variations: we find lower stress drops in the Middle Mountain asperity and along the creeping fault section, and higher stress drops in the hypocentral region of the 2004 M6.0 Parkfield earthquake. The main shock did not alter the overall pattern of high and low stress drop regions. However, a statistical test reveals areas with significant changes in computed stress drops after the main shock, which we compare to estimated absolute shear stress changes from a main shock slip model. By calculating Δt* from the spectral EGF ratio, we also identify areas with increased attenuation after the main shock, and we are able to distinguish source effects and near‐source attenuation effects in the spectral analysis. These results are confirmed independently from spectral ratios of repeating microearthquake clusters.