1 / 14

RE(CALIBRATION) OF MAGNITUDE SCALES FOR THE CROATIAN SEISMographic NETWORK

M. Herak, S. Markušić, D. Herak Department of Geophysics, Faculty of Science University of Zagreb, Zagreb. RE(CALIBRATION) OF MAGNITUDE SCALES FOR THE CROATIAN SEISMographic NETWORK. History and motivation.

guy
Download Presentation

RE(CALIBRATION) OF MAGNITUDE SCALES FOR THE CROATIAN SEISMographic NETWORK

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. M. Herak, S. Markušić, D. Herak Department of Geophysics, Faculty of Science University of Zagreb, Zagreb RE(CALIBRATION) OF MAGNITUDE SCALES FOR THE CROATIAN SEISMographic NETWORK

  2. History and motivation... The first magnitude calibration in Croatia was done by D. Skoko in his Ph. D. dissertation in 1969. Two calibrating functions for recordings on 1000 kg Wiechert horizontal seismograph were derived – for distances of 2.2˚–20˚, and 20˚–170˚. In the ‘Balkan project’ (1974), a new formula was derived. It was retained even after introduction of electromagnetic and digital seismographs: After a significant enlargement of the network in the few first years of the 21st century, Futač (2007) checked its appropriatness, and concluded that data do not indicate the need for change. As digital high-quality data accumulated, it became clear that the formula does not work well for short distances, significantly underestimating the magnitude. This prompted recalibration of the ML-scale. The need to use moment magnitudes in hazard analyses, motivated efforts to routinely observe also Mw thus initiating this work.

  3. Goals and earthquakes • Goals: • Calibration of MLand MWA using recent digital data from the Croatian Seismographic Network, and • Introduction of an automatic routine procedure for determination of MW using spectral analyses of local and regional earthquakes recorded by the Croatian Seismographic Network • Selection of earthquakes: • Events from the Croatian Earthquake Catalogue satisfying any of the following: • Year > 2002 • M > 2.5 & D < 200 km (for any station) • M > 3.0 & D < 300 km (for any station) • M > 3.5 & D < 400 km (for any station) • M > 4.0 & D < 500 km (for any station) • M > 5.0 & D < 700 km (for any station) • 11900 3-component BB seimograms found!

  4. Goals and earthquakes 1. Determination of seismic moment (M0) and moment magnitude (Mw) The source displacement spectrum looks like: v – velocity of seismic wavesf – frequency f0 – corner frequency  – density low frequency spectral level amplitude corner frequency, f0 frequency

  5. M0 and Mw 1. Determination of seismic moment (M0) and moment magnitude (Mw) Displacement spectrum at receiver at the epicentral distance :  – near surface attenuationQ – quality factor (for P or S-waves) RP – average radiation pattern coefficient ( 0.6) FS – free surface amplification ( 2 for SH) G = 1/gd – geometrical spreading function, gd – geodistance If we know (or invert for!) attenuation (, Q), and assume reasonable values representing the effects caused by the radiation pattern, free surface amplification and geometrical spreading, we can fit the formula to the observed spectrum and get the scalar seismic moment (M0), and then compute the moment magnitude: Mw= 2/3 log(M0) – 6.07

  6. M0and Mw 1. Determination of seismic moment (M0) and moment magnitude (Mw) Determine theoretical onset times of Pg, Pn, Sg, Sn. LP-filter, read max. S-wave amplitudes. Wood-Anderson filter, read max. S-wave amplitudes Determine the windows for: noise, P-waves, S-waves. Compute spectra, smooth, correct for the instrument and noise. Fitting theoretical curve to spectra for P- and S-waves yields estimates for: Q(f),, f0, M0, both for P- and S-waves Compute MW as average of MW(P) and MW(S) Save everything, next earthquake S-waves P-waves

  7. M0 and Mw; ML and MWA 2. Calibration of ML and the Wood-Anderson magnitude MWA Compute MLand MWA for every recording (all stations) using the best calibrating function so far, in the form: Compute representative earthquake magnitude for all events for which at least three magnitude estimates exist as the median of the station magnitudes. Check if there is any distance dependence – if yes adjust coefficients observing all possible constraints (e.g. anchoring of MWA at 100 km). Compute station corrections (SC) and start again... The best estimates are:Amax is the trace amplitude in nm on the simulated Wood-Anderson seismograph (magnification 2080). This is the same MWAformulaas obtained for Central Europe by Stange (2006).

  8. M = (Station magnitude) – (Earthquake magnitude) vs. time MWMLMWA MWMLMWA

  9. M = (Station magnitude) – (Earthquake mgnitude) vs. time MWMLMWA MWMLMWA

  10. M0and Mw – results Individual (station) magnitudes Corner frequency vs MW Earthquake magnitudes S-waves P-waves MWS vs MWP

  11. M0 and Mw; ML and MWA (Event magnitude) – (Station magnitude) (with added station corrections)

  12. Relationships... Individual (station) magnitudes Earthquake magnitudes MW vs. MWA ML vs. MWA MW vs. ML

  13. Station corrections Sta MLMWAMwSMwP ---------------------------- KIJV 0.02 0.06 -0.03 -0.03 CACV 0.04 0.06 -0.12 -0.08 PTJ 0.09 0.06 -0.05 -0.00 DBR 0.06 0.15 0.01 0.06 STON 0.05 0.08 0.01 0.06 STA 0.10 0.14 -0.01 -0.04 NVLJ 0.05 0.05 -0.05 -0.08 RIY -0.06 -0.12 -0.13 0.02 ZAG -0.27 -0.32 -0.21 -0.16 SLUN 0.06 0.07 -0.05 -0.04 UDBI -0.06 -0.11 -0.12 -0.05 BRJN 0.04 0.04 -0.03 -0.10 SISC -0.40 -0.41 -0.20 -0.24 HVAR 0.04 0.00 -0.02 0.12 MORI -0.04 -0.01 -0.07 -0.08 ZIRJ -0.03 -0.01 -0.03 -0.04 OZLJ 0.15 0.16 -0.13 -0.01 KALN 0.07 0.01 -0.05 -0.03 RIC -0.04 -0.11 -0.18 -0.01 KSY -0.10 -0.11 -0.17 -0.10 SLNJ 0.14 0.04 -0.01 0.14 ---------------------------- The largest (negative) station corrections are found on thick alluvium, in river valleys close to the Pannonian basin. ML station corrections Corrections for the moment magnitudes are smaller than for local magnitudes (!)

  14. Conclusions Regression of Mw vs. ML using both magnitudes determined independently on the same set of stations and earthquakes is far better than regressions of locally determined ML vs. Mwdetermined (mostly for large and distant earthquakes!) by other agencies. Proposed inversion of observed spectral shapes of local and regional earthquakes results in consistent Mw estimation using data collected by the Croatian Seismographic Network. New calibrating functions for local magnitudes (ML and MWA) yield distance-independent estimates. The correspondence between the three magnitudes is close to 1:1 relationship. This is encouraging, promising an easy conversion from MLto MW for older events. Results obtained in studies done so far using MLas proxy for MWwill most probably not significantly change after proper conversion.

More Related