Routine Earthquake Locations Research Articles

Reference seismic crustal model of the Dinarides

Abstract. Continental collision zones are structurally one of the most heterogeneous areas intermixing various different units within a relatively small space. A good example of this is the Dinarides, a mountain chain situated in the central Mediterranean, where thick carbonates cover older crystalline basement units and remnants of subducted oceanic crust. This is further complicated by the highly variable crustal thickness ranging from 20 to almost 50 km. In terms of spatial extension, this area is relatively small but covers tectonically differentiated domains making, any seismic or geological analysis complex, with significant challenges in areas that lack seismic information on crustal structure. Presently there is no comprehensive 3D crustal model of the Dinarides (and surrounding areas). Using the compilations of previous studies and employing kriging interpolation, we created a vertically and laterally varying crustal model defined on a regular grid for the wider area of the Dinarides, also covering parts of Adriatic Sea and the SW part of the Pannonian Basin. The model is divided by three interfaces, Neogene deposit bottom, carbonate rock complex bottom and Moho discontinuity, with seismic velocities (P and S waves) and density defined at each grid point. To validate the newly derived model, we calculated travel times for an earthquake recorded on several seismic stations in the Dinarides area. The calculated travel times show significant improvement when compared to the simple 1D model used for routine earthquake location in Croatia. The model derived in this work represents the first step towards improving our knowledge of the crustal structure in the complex area of the Dinarides. We hope that the newly assembled model will be useful for all forthcoming studies (e.g., as a starting model for seismic tomography, as a model for earthquake simulations) which require knowledge of the crustal structure.

Minimum 1D VP and VP/VS Models and Hypocentral Determinations in the Central Mediterranean Area

Abstract 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.

One-dimensional velocity structure modeling of the Earth's crust in the northwestern Dinarides

Abstract. The studied area of the northwestern (NW) Dinarides is located in the northeastern (NE) corner of the Adriatic microplate and is bordered by the Adriatic foreland, the Southern Alps, and the Pannonian basin. Its complex crustal structure is the result of interactions among different tectonic units, the most important of which are the Eurasian plate and the Adriatic microplate. Despite numerous seismic studies in this tectonically complex area, there is still a need for a detailed, small-scale study focusing mainly on the upper, brittle part of the crust. In this work, we investigated the velocity structure of the crust with one-dimensional (1-D) simultaneous hypocenter–velocity inversion using routinely picked P- and S-wave arrival times. Most of the models computed in the combined P and S inversion converged to a stable solution in the depth range between 0 and 26 km. We further evaluated the inversion results with hypocenter shift tests, high- and low-velocity tests, and relocations. This helped us to select the best performing velocity model for the entire study area. Based on these results and the seismicity distribution, we divided the study area into three subregions, reselected earthquakes and stations, and performed the combined P and S inversion for each subregion separately to gain better insight into the crustal structure. In the eastern subregion, the P velocities in the upper 8 km of the crust are lower compared to the regional velocities and the velocities of the other two subregions. The P velocities between 8 and 23 km depth are otherwise very similar for all three models. Conversely, the S velocities between 2 and 23 km depth are highest in the eastern subregion. The NW and southwestern (SW) subregions are very similar in terms of the crustal structure between 0 and 23 km depth, with slightly higher P velocities and lower S velocities in the SW subregion. High vP/vS values were obtained for the layers between 0 and 4 km depth. Below that, no major deviations of vP/vS in the regional model from the value of 1.73 are observed, but in each subregion we can clearly distinguish two zones separated by a decrease in vP/vS at 16 km depth. Compared to the model currently used by the Slovenian Environment Agency to locate earthquakes, the obtained velocity models show higher velocities and agree very well with some of the previous studies. In addition to the general structural implications and the potential to improve the results of seismic tomography, the new 1-D P and S velocity models can also be used for reliable routine earthquake location and for detecting systematic travel time errors in seismological bulletins.

ARRU Phase Picker: Attention Recurrent-Residual U-Net for Picking SeismicP- andS-Phase Arrivals

AbstractSeismograms 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.

A new 3-D P-wave velocity model for the Gulf of Cadiz and adjacent areas derived from controlled-source seismic data: application to nonlinear probabilistic relocation of moderate earthquakes

SUMMARYThe SW Iberian margin is well known for its complex tectonic setting and crustal structure and by the occurrence of moderate magnitude earthquakes and some great tsunamigenic earthquakes. Fortunately, many seismic reflection and refraction profiles have been carried out, providing detailed information about the crustal structure of the main geologic domains in this region. These studies show a first-order variation due to the transition from oceanic to continental domain, large-scale heterogeneities within the crust and an irregular Moho topography. Routine earthquake locations in this area have been usually computed using a general 1-D velocity model which is clear that cannot account for such a heterogeneous structure. In addition, regional seismic stations used to locate the Gulf of Cadiz seismicity are on land and far away to the east, implying large azimuthal gaps and distances. In this context, a 3-D approach seems necessary to properly solve the crustal velocity field and improve earthquake location in this area. With this purpose, we present a new digital 3-D P-wave velocity distribution for the crust and uppermost mantle derived from previously published controlled-source seismic experiments carried out in SW Iberia and the Gulf of Cadiz over the last 40 yr. We have reviewed more than 50 wide-angle and multichannel seismic reflection and refraction profiles and digitized the most significant published 2-D seismic velocity models, performing an updated compilation of crustal parameters (P-wave velocities and geometry and depth of the main crustal interfaces). These velocities as a function of position and depth have been interpolated using ordinary kriging algorithm to obtain, in the form of a regular georeferenced 20 × 20 × 1 km grid spacing, a high-resolution 3-D P-wave velocity distribution for the crust and uppermost mantle and a continuous Moho depth map of the whole area of this study (33°N–41°N latitude and 15°W–5°W longitude). Since current seismic location tools allow the implementation of 3-D grid structures, we have applied our 3-D model to relocate a selection of moderate earthquakes occurred in the studied region using a probabilistic nonlinear method. In the Gulf of Cadiz area the probabilistic approximation provides maximum likelihood hypocentres located within the uppermost mantle with the majority of depths ranging between 20 and 45 km. This model would subsequently be implemented at the Spanish Seismic Network for the routine relocation of the seismicity of the area.

Tsunami Potential of Moderate Earthquakes: The July 1, 2009 Earthquake (Mw 6.45) and its Associated Local Tsunami in the Hellenic Arc

On July 1, 2009, a Mw 6.45 earthquake ruptured south of Crete Island (Greece) along the Hellenic Arc triggering a local tsunami. Eyewitness reported the tsunami from Myrtos and Arvi (south-eastern Crete) and from the north of Chrisi Islet, located to the southeast of Crete. The earthquake occurred offshore, about 80 km south of Crete, where routine earthquake locations are poor. The hypocentre is relocated using a 2-D velocity model and several local 1-D velocity models. Epicentral locations obtained by using the different velocity models show very minor variations. Instead, relocated hypocentres can be grouped into two sets of solutions: (1) those with a shallower depth (depth 25 km) either in the upper-plate or along the plate interface. Models accounting for either homogeneous or heterogeneous slip on the causative fault are tested, with the heterogeneous one better reproducing the observations in terms of number of tsunami waves reaching the shoreline and duration of the sea disturbance. The short travel time, about 10 min, of the first tsunami arrival at the southern coast of Crete represents a big challenge for tsunami early warning systems operating in the area.

Regional minimum 1-D P-wave velocity model for a new seismicity catalogue with precise and consistent earthquake locations in southern Iran

For faster and more robust ray tracing in 1-D velocity models and also due to the lack of reliable 3-D models, most seismological centers use 1-D models for routine earthquake locations. In this study, as solution to the coupled hypocenter-velocity problem, we compute a regional P-wave velocity model for southern Iran that can be used for routine earthquake location and also a reference initial model for 3-D seismic tomography. The inversion process was based on travel time data from local earthquakes paired reports obtained by merging the catalogues of Iranian Seismic Center (IRSC, 6422 events) and by the Broadband Iranian National Seismic Network (BIN, 4333 events) for southern Iran in the period 2006 through July 2017. After cleaning the data set from large individual reading errors and by identifying event reports from both networks belonging to same earthquake (a process called event pairing), we obtained a data set of 1115 well-locatable events with a total number of 24,606 P-wave observations. This data set was used to calculate a regional minimum 1-D model for southern Iran as result of an extensive model search by trial-and-error process including several dozens of inversions. Significantly different from previous models, we find a smoothly increasing P-velocity by depth with velocities of 5.8 km/s at shallow and velocities of 6.4 km/s at deepest crustal levels. For well-locatable events, location uncertainties are estimated in the order of ± 3 km for epicenter and double this uncertainty for hypocentral depth. The use of the minimum 1-D model with appropriate station delays in routine hypocenter location processing will yield a high-quality seismic catalogue with consistent uncertainty estimates across the region and it will also allow detection of outlier observations. Based on the two catalogues by IRSC and BIN and using the minimum 1-D model and station delays for all stations in the region, we established a new combined earthquake catalogue for southern Iran. While the general distribution of the seismicity corresponds well with that of the two individual catalogues by IRSC and BIN, the new catalogue significantly enhances the correlation of seismicity with the regional fault systems within and between the major crustal blocks that as an assembly build this continental region. Furthermore, the unified seismic catalogue and the minimum 1-D model resulting from this study provide important ingredients for seismic hazard studies.

Improving Seismic Surveillance at Mt. Etna Volcano by Probabilistic Earthquake Location in a 3D Model

The increasing accuracy of 3D velocity models developed recently for Mt. Etna has enabled their use today in routine earthquake locations. In this work, we tested the potential and performance of a global-search probabilistic earthquake lo- cation method (NonLinLoc) in a 3D velocity model, to improve earthquake locations for seismic surveillance. In addition, NonLinLoc hypocenter locations and those ob- tained by standard iterative-linear 3D locations, SimulPS-14, have also been com- pared. To this end, a dataset of 328 selected earthquakes, occurring during the 2002-2003 Etna flank eruption, and the recent highest resolution 3D velocity model, have been used. The results revealed that the differences in hypocentral coordinates between the two methods are typically of the same order or smaller than the spatial location uncertainty. To evaluate the consistency of results between the two 3D location algorithms, synthetic datasets with real source-receiver configuration are also considered. Furthermore, by using NonLinLoc we estimated the influence of the source-receiver geometry on the quality of hypocenter locations. If we vary the net- work geometry in a dense and well-distributed network like at Etna, reducing the number of stations (by 20% and 50%), it is significant that no large systematic hypo- central shifts of the relocated earthquakes are observed if they occur within the net- work. NonLinLoc is a fast and promising approach for automatic earthquake locations and surveillance purposes at Mt. Etna, because (1) it works well with a reduced num- ber of seismic pickings, which are usually available in the automatic locations; (2) it is not particularly sensitive to tolerable levels of random noise in arrival times; and (3) it produces full location uncertainty and resolution information with respect to standard iterative-linear 3D locations.

Constraints on the crustal structure beneath the Sinai subplate, SE Mediterranean, from analysis of local and regional travel times

The Sinai Peninsula has been recognized as a subplate of the African Plate located at the triple junction of the Gulf of Suez rift, the Dead Sea Transform fault, and the Red Sea rift. The upper and lower crustal structures of this tectonically active, rapidly developing region are yet poorly understood because of many limitations. For this reason, a set of P- and S-wave travel times recorded at 14 seismic stations belonging to the Egyptian National Seismographic Network (ENSN) from 111 local and regional events are analyzed to investigate the crustal structures and the locations of the seismogenic zones beneath central and southern Sinai. Because the velocity model used for routine earthquake location by ENSN is one-dimensional, the travel-time residuals will show lateral heterogeneity of the velocity structures and unmodeled vertical structures. Seismic activity is strong along the eastern and southern borders of the study area but low to moderate along the northern boundary and the Gulf of Suez to the west. The crustal Vp/Vs ratio is 1.74 from shallow (depth ≤ 10 km) earthquakes and 1.76 from deeper (depth > 10 km) crustal events. The majority of the regional and local travel-time residuals are positive relative to the Preliminary Reference Earth Model (PREM), implying that the seismic stations are located above widely distributed, tectonically-induced low-velocity zones. These low-velocity zones are mostly related to the local crustal faults affecting the sedimentary section and the basement complex as well as the rifting processes prevailing in the northern Red Sea region and the ascending of hot mantle materials along crustal fractures. The delineation of these low-velocity zones and the locations of big crustal earthquakes enable the identification of areas prone to intense seismotectonic activities, which should be excluded from major future development projects and large constructions in central and southern Sinai.

A Catalog of Regional Moment Tensors in Utah from 1998 to 2011

Online Material: The electronic supplement includes details for all 48 moment tensor solutions in this study, including waveform fits and a table of event parameters. The University of Utah Seismograph Stations (UUSS) operates regional seismic networks in and around Utah and Yellowstone National Park and provides routine earthquake locations and earthquake catalogs for these two areas (Pechmann et al. , 2007). In this study, we expand the reporting to include moment tensors for seismic events occurring in the Utah region (Fig. 1). This region is a complex seismic environment with diverse source types: tectonic earthquakes, mining‐induced seismicity (MIS) including mine collapses, and man‐made sources such as explosions. Most tectonic earthquakes in the region occur within the Intermountain Seismic Belt (ISB; Smith and Arabasz, 1991). In northern Utah, the ISB roughly delineates the boundary between the Basin and Range and the Middle Rocky Mountains, whereas in central and southern Utah, the ISB separates the actively extending Basin and Range from the more stable Colorado Plateau. The MIS is mostly associated with long‐wall coal mining in central Utah (Fig. 1; Arabasz et al. , 2007). The man‐made blasts occur throughout the state and are related to both surface mining operations and military munitions detonations (Arrowsmith et al. , 2008). Figure 1. UUSS Utah monitoring area (dark gray box) and the surrounding region. Solid black line separates the UUmod velocity model regions between Basin and Range to the northwest and Colorado Plateau to the southeast. Dashed line separates physiographic provinces: Basin and Range (BR), Colorado Plateau (CP), and Middle Rocky Mountains (MR). The transition zone is roughly the overlap between the physiographic and velocity model boundaries between the Basin and Range and Colorado Plateau in Central Utah. Crescent‐shaped and rectangular polygons in central‐eastern Utah are coal mining areas. ECSZ is the Eastern California …

Characteristics of the seismicity of Vesuvius and Campi Flegrei during the year 2000

This paper describes the characteristics of the seismicity in the volcanic Neapolitan area during the year 2000 recorded by the monitoring seismic network of the Osservatorio Vesuviano. In particular, a detailed analysis of the seismicity of Vesuvius is presented. We compared the seismic velocity models available for the Vesuvius area locating the earthquakes recorded in the year 2000 and on the basis of the results, we introduce for routine earthquake location the new velocity model obtained by the seismic tomography experiments (TomoVes) performed in the area. We also determined the focal mechanisms and analysed the seismicity rate, comparing the results with those obtained for the past years. After the introduction of the new acquisition system at the Osservatorio Vesuviano, a re-calibration of the duration magnitude scale was necessary to avoid biases related to the different instrumental response. Consequently, we re-calibrated the magnitude relation used for the Vesuvius earthquakes, obtaining a new formula to be used for the earthquakes recorded by the new acquisition system. Finally, we give a description of the seismic activity in the Campi Flegrei area during the summer of 2000.

Accurate hypocentre locations in the Middle-Durance Fault Zone, South-Eastern France

Abstract A one-dimensional velocity model and station corrections for the Middle-Durance fault zone (south-eastern France) were computed by inverting P-wave arrival times recorded on a local seismic network of 8 stations. A total of 93 local events with a minimum of 6 P-phases, RMS 0.4 s and a maximum gap of 220° were selected. Comparison with previous earthquake locations shows an improvement for the relocated earthquakes. Tests were carried out to verify the robustness of inversion results in order to corroborate the conclusions drawn from our findings. The obtained minimum 1-D velocity model can be used to improve routine earthquake locations and represents a further step toward more detailed seismotectonic studies in this area of south-eastern France.

Repeating earthquakes recorded by Liaoning Regional Seismograph Network

In the list of ‘repeating pairs’ or ‘doublets’ of earthquakes in China identified by Schaff and Richards using teleseismic waveform cross-correlation, there were 23 repeating pairs located in Liaoning Province. In this study the waveforms of these events were cross-correlated using records from Liaoning Regional Seismograph Network (LRSN), and the ‘repeating events’ in the sense of regional waveform cross-correlation were obtained. The result was compared with that of Schaff and Richards and was used for the assessment of the seismic phase picking and event location practice of LRSN. The result shows that ‘repeating events’ in the sense of teleseismic waveform cross-correlation and those in the sense of regional waveform cross-correlation have significant difference, although with some overlap. However, the overall assessment of the location accuracy and the phase pick errors of LRSN by using these two sets of ‘repeating events’, respectively, provides similar results, while ‘repeating events’ in the sense of regional waveform cross-correlation seem to be better performing in such an assessment. With the assumption that the separation between the ‘repeaters’ be less than 1 km, the uncertainty in routine earthquake location of LRSN is estimated to be below 5 km, with the average of 2 km. In the observational bulletins of LRSN the time error in phase picking is estimated to be within ±1s for 94% Pg readings and for 88% Sg readings.

Crustal seismic velocity in the Marche region (Central Italy): computation of a minimum 1-D model with seismic station corrections

A 1-D velocity model for the Marche region (Central Italy) was computed by inverting P- and S-wave arrival times of local earthquakes. A total of 160 seismic events with a minimum of ten observations, a travel time residual ≤0.8s and an azimuthal gap lower than 180° have been selected. This “minimum 1-D velocity model” is complemented by station corrections, which can be used to take into account possible near-surface velocity heterogeneities beneath each station. Using this new P-wave velocity model and the program HYPOELLIPSE (Lahr 1999), the selected local events were relocated. Earthquake locations in this study are of higher quality with respect to the original ones. The obtained minimum 1-D velocity model can be used to improve the routine earthquake locations and represents a further step towards more detailed seismotectonic studies of the area.

A Simple Algorithm for Local Earthquake Location Using 3D VP and VS Models: Test Examples in the Central United States and in Central Eastern Taiwan

Traditional local-earthquake location using a horizontally layered homogeneous velocity model is limited in its resolution and reliability due to the existence of frequently overlooked 3D complexity of the real Earth. During traditional 3D seismic tomography, simultaneous earthquake relocation using the resultant 3D velocity model has produced reliable earthquake locations; however, only a small subset of events are typically used and thus relocated in the inversion. The rest of the events in a catalog must then be relocated using the 3D models. The repeated calculation of travel times across 3D V P and V S models is also not efficient and not practical for a routine network earthquake location when the very time-consuming exact 3D raytracing is used. Because high-resolution earthquake data are now available from many modern seismic networks, representative high-resolution 3D V P and V S models for a region can be better determined. By taking advantage of recently available high-speed computer technology and large disk space, we implemented a simple algorithm to efficiently locate every local earthquake using the best available regional 3D V P and V S models. Once the V P and V S information for all cubic cells in a 3D grid model are determined, P and S travel times from each grid point to all seismic stations can be calculated and stored on disk files for later usage. During the iteration process for earthquake location, travel times from a trial hypocenter to all recording stations can be determined simply by a linear interpolation from those of the adjacent eight grid points available in the previously stored disk files without the need for raytracing. The iterations continue until the hypocenter adjustments at the end of the last iteration are below the given criteria and the travel-time residual, or the difference between the observed and the calculated travel times, is a minimum. Therefore, any local earthquake can be efficiently and reliably located using the available 3D velocity models. This simple location program has been applied to relocate earthquakes in the New Madrid Seismic Zone (nmsz) of the central United States and in the central eastern Taiwan region. Preliminary results in both regions reveal that earthquake hypocenters can be efficiently relocated in spite of the very significant lateral structural variations. Tests with data from Taiwan further demonstrate that the resolution of seismic tomography and the relocated seismicity is sensitive to relative distribution of seismic-network stations and background seismicity. Thus, this single-event location program can be applied to relocate all earthquakes in a seismic-network catalog and, more importantly, to allow routine earthquake location for any seismic network using the available 3D velocity models.

Relocation of the 1998 Zhangbei-Shangyi earthquake sequence using the double difference earthquake location algorithm

On January 10, 1998, at 11h50min Beijing Time (03h50min UTC), an earthquake ofML=6.2 occurred in the border region between the Zhangbei County and Shangyi County of Hebei Province. This earthquake is the most significant event to have occurred in northern China in the recent years. The earthquake-generating structure of this event was not clear due to no active fault capable of generating a moderate earthquake was found in the epicentral area, nor surface ruptures with any predominate orientation were observed, no distinct orientation of its aftershock distribution given by routine earthquake location was shown. To study the seismogenic structure of the Zhangbei-Shangyi earthquake, the main shock and its aftershocks withML≥3.0 of the Zhangbei-Shangyi earthquake sequence were relocated by the authors of this paper in 2002 using the master event relative relocation technique. The relocated epicenter of the main shock was located at 41.145°N, 114.462°E, which was located 4 km to the NE of the macro-epicenter of this event. The relocated focal depth of the main shock was 15 km. Hypocenters of the aftershocks distributed in a nearly vertical plane striking 180°200° and its vicinity. The relocated results of the Zhangbei-Shangyi earthquake sequence clearly indicated that the seismogenic structure of this event was a NNE-SSW-striking fault with right-lateral and reverse slip. In this paper, a relocation of the Zhangbei-Shangyi earthquake sequence has been done using the double difference earthquake location algorithm (DD algorithm), and consistent results with that obtained by the master event technique were obtained. The relocated hypocenters of the main shock are located at 41.131°N, 114.456°E, which was located 2.5 km to the NE of the macro-epicenter of the main shock. The relocated focal depth of the main shock was 12.8 km. Hypocenters of the aftershocks also distributed in a nearly vertical N10°E-striking plane and its vicinity. The relocated results using DD algorithm clearly indicated that the seismogenic structure of this event was a NNE-striking fault again.

One-dimensional P-wave velocity structure in the northern Red Sea area, deduced from travel time data

Travel time inversion of local earthquake data was used to improve a one-dimensional (1-D) P-wave velocity structure model for the northern Red Sea area, including the entrance of the Gulf of Suez and the Gulf of Aqaba. The 1-D velocity model and station corrections were obtained by using 1538 P-arrivals from 216 well-located earthquakes recorded at ten local seismic stations from August 1994 to January 1999. Earthquakes used in the present study were located in depths less than 28 km so that the upper 20 km of the 1-D velocity model was well constrained. The resulting P-velocity model indicated velocities from 5.0 to 6.0 km/sec down to a depth of 10 km, and velocities from 6.0 to 6.8 km/sec with a nearly constant gradient in a depth range from 10 to 20 km. This 1-D velocity model reduced the RMS residual by 47% from 0.21 to 0.11 sec and improved routine earthquake locations, which clearly indicated clusters of hypocenters in the southern tip of the Sinai Peninsula and the entrance of the Gulf of Suez. The resulting station corrections suggested irregular surface geology in the area. Negative and positive station corrections were obtained; the minimum value (−0.35) was observed at the station MAZR, while the maximum value (0.24) at SHRM.

Improved seismic velocity reference model from local earthquake data in Northwestern Italy

ABSTRACTResolution and reliability estimates of results obtained by seismic tomography strongly depend on the reference model. Inadequate initial reference models may severely distort tomographic images or introduce artefacts that lead to misinterpretations of the results. Reference models are usually obtained by means of a priori near‐surface geological information or by geophysical information derived by controlled‐source seismology.Starting from the idea that a reference model must approximate the weighted average of data selected for the three‐dimensional (3D) inversion, one‐dimensional (1D) model for Northwestern Italy is derived that is able to minimize mean of RMS of a set of well‐locatable earthquakes, by computing a solution of the coupled hypocentre 1D velocity problem.Such a model, termed the Minimum 1D model, can be used both as an initial reference model for 3D inversion and as a reference velocity model for high‐quality routine earthquake location.

Geotomography with local earthquake data

The inversion of local earthquake data (LED) for three‐dimensional velocity structure requires the simultaneous solution of the coupled hypocenter‐model problem. The Aki‐Christoffersson‐Husebye method (ACH) involves the inversion of large matrices, a task that is often performed by approximative solutions when the matrices become too big, as is the case for most LED, considering the coupled inverse problem. Such an approximate method (herein referred to as approximate geotomographic method) is used to perform tests with LED to obtain the best suited inversion parameters, such as velocity damping and number of iteration steps. The ACH method has been proposed for use of teleseismic data. Several adjustments to the original ACH method, which are necessary for use of LED, have been developed and are discussed. Such adjustments are the separation of the unknown hypocentral from the velocity model parameters for the inversion, the use of geometric weighting and step length weighting, the calculation of a minimum one‐dimensional (1D) model as the starting three‐dimensional (3D) model for the model inversion, and the display of an approximate resolution matrix (ray density tensors) before the inversion is performed. The ray density tensors allow the block cutting, e.g., the definition of the 3D velocity grid, to better correspond with the resolution capability of the specific data set. The adjustments to the method are tested by inversion of realistic LED of known variance. Synthetic LED are also used to demonstrate the effects of systematic errors, such as mislocations of seismic stations, on the resulting velocity field. Using the data sets from Long Valley, California, Yellowstone National Park, Wyoming, and Borah Peak, Idaho, the effects of improvements to the ACH method and of the data filtering process are shown. The use of the minimum 1D models for routine earthquake location improves this location procedure, as shown with the relocation of shots for the Long Valley and Yellowstone areas. The three‐dimensional velocity fields obtained by the ACH method for the Long Valley and Yellowstone areas show local anomalies in the p velocity that can be correlated with tectonic and volcanic features. A pronounced anomaly of low p velocity below the Yellowstone caldera can be interpreted as a large magma chamber. However, the bulk of the paper addresses problems of the inversion method. The LED from the areas mentioned above are used to numerically and theoretically tune the inversion method for the defects that all real data contain. It is shown that one of the most important steps for any inversion of LED is the selection of the data for quality and for geometrical distribution.