Seismic source and structure estimation in the western Mediterranean using a sparse broadband network

Abstract

Abstract. We present a study of regional earthquakes in the western Mediterranean geared toward the development of methodologies and path calibrations for source characterization using regional broadband stations. The results of this study are useful for the monitoring and discrimination of seismic events under a comprehensive test ban treaty, as well as the routine analysis of seismicity and seismic hazard using a sparse array of stations. The area consists of several contrasting geological provinces with distinct seismic properties, which complicates the modeling of seismic wave propagation. We started by analyzing surface wave group velocities throughout the region and developed a preliminary model for each of the major geological provinces. We found variations of crustal thickness ranging from 45 km under the Atlas and Betic mountains and 37 km under the Saharan shield, to 20 km for the oceanic crust of the western Mediterranean Sea, which is consistent with earlier works. Throughout most of the region, the upper mantle velocities are low which is typical for tectonically active regions. The most complex areas in terms of wave propagation are the Betic Cordillera in southern Spain and its north African counterparts, the Rif and Tell Atlas mountains, as well as the Alboran Sea, between Spain and Morocco. The complexity of the wave propagation in these regions is probably due to the sharp velocity contrasts between the oceanic and continental regions as well as the the existence of deep sedimentary basins that have a very strong influence on the surface wave dispersion. We used this preliminary regionalized velocity model to correct the surface wave source spectra for propagation effects which we then inverted for source mechanism. We found that this method, which is in use in many parts of the word, works very well, provided that data from several stations are available. In order to study the events in the region using very few broadband stations or even a single station, we developed a hybrid inversion method which combines Phi waveforms synthesized with the traditional body wave methods, with surface waves that are computed using normal modes. This procedure facilitates the inclusion of laterally varying structure in the Green's functions for the surface waves and allows us to determine source mechanisms for many of the larger earthquakes (M > 4) throughout the region with just one station. We compared our results with those available from other methods and found that they agree quite well. The epicentral depths that we have obtained from regional waveforms are consistent with observed teleseismic depth phases, as far as they are available. We also show that the particular upper mantle structure under the region causes the various Pn and Sn phases to be impulsive, which makes them a useful tool for depth determination as well. Thus we conclude that with proper calibration of the seismic structure in the region and high-quality broadband data, it is now possible to characterize and study events in this regio...

Calibration of the Regional Crustal Waveguide and the Retrieval of Source Parameters Using Waveform Modeling

“…This characteristic of P nl waves is due to their primarily sampling the lower crust and the Moho interface, and consequently being relatively insensitive to lateral crustal variations, which occur mainly in the upper crust. Even for mixed paths traversing continental and oceanic crust in the western Mediterranean, THIO et al (1999) found P nl waveforms to remain quite stable. Hence this portion of the wavetrain can be adequately modeled over long distances with a single model, although its relative timing with respect to the rest of the wavetrain may vary.…”

“…This is especially true for Rayleigh waves (WOODS and HARKRIDER, 1995;SONG et al, 1996;THIO et al, 1999). Therefore, it becomes dicult to model the entire wavetrain waveform with a single 1-D model.…”

“…Therefore, it becomes dicult to model the entire wavetrain waveform with a single 1-D model. An example of such behavior is provided in Figure 3, taken from THIO et al (1999), which compares the vertical and tangential displacement records for an earthquake in Gafsa in Tunisia (June 12, 1992, origin time: 19 h 16 m 45.7 s) recorded by the stations MEB (Medea, Algeria) and VSL (Villasalto, on the Island of Sardina). The path to MEB is a pure continental one whereas the one to VSL consists of both continental and oceanic sections almost equal in length.…”

“…Recent studies have demonstrated that regional path calibration is essential for obtaining accurate locations of seismic events at all magnitudes (RUSSIAN FEDER-ATION/UNITED STATES CALIBRATION WORKING GROUP, 1998;SAIKIA et al, 1996a, b;THIO et al, 1999;ZHU et al, 1997;RYABOY, 1999). Accurate locations are integral to Comprehensive Test-Ban Treaty (CTBT) seismic monitoring eorts, particularly for smaller events, and are also crucial to many types of calibration studies.…”

“…M 0 is also an important parameter for accurate singlestation path calibrations of regional seismic-phase amplitudes used for discrimination (PRIESTLEY and PATTON, 1997) and reliable magnitude measurements. Furthermore, regional analysis of shallow crustal events can improve source depth estimates as teleseismic P waves can be contaminated by source-end uppercrustal reverberations, making it dicult to identify emergent ®rst arrivals (STEIN and WIENS, 1986;THIO et al, 1999). Well calibrated regional crustal waveguides are also important to modeling or inversion of teleseismic P waves (P, pP and sp) to determine source depth (NELSON et al, 1987;WOODS et al, 1998), as teleseismic Pwave based hypocenters can have signi®cant error in depth, even for larger events (M w > 5).…”