A Surface Wave Dispersion Study of the Middle East and North Africa for Monitoring the Comprehensive Nuclear-Test-Ban Treaty

—?We present results from a large-scale study of surface-wave group velocity dispersion across the Middle East, North Africa, southern Eurasia and the Mediterranean. Our database for the region is populated with seismic data from regional events recorded at permanent and portable broadband, three-component digital stations. We have measured the group velocity using a multiple narrow-band filter on deconvolved displacement data. Overall, we have examined more than 13,500 seismograms and made good quality dispersion measurements for 6817 Rayleigh- and 3806 Love-wave paths. We use a conjugate gradient method to perform a group-velocity tomography. Our current results include both Love- and Rayleigh-wave inversions across the region for periods from 10 to 60 seconds. Our findings indicate that short-period structure is sensitive to slow velocities associated with large sedimentary features such as the Mediterranean Sea and Persian Gulf. We find our long-period Rayleigh-wave inversion is sensitive to crustal thickness, such as fast velocities under the oceans and slow along the relatively thick Zagros Mts. and Turkish-Iranian Plateau. We also find slow upper mantle velocities along known rift systems. Accurate group velocity maps can be used to construct phase-matched filters along any given path. The filters can improve weak surface wave signals by compressing the dispersed signal. The signals can then be used to calculate regionally determined M _S measurements, which we hope can be used to extend the threshold of m _b :M _S discriminants down to lower magnitude levels. Other applications include using the group velocities in the creation of a suitable background model for forming station calibration maps, and using the group velocities to model the velocity structure of the crust and upper mantle.     

Rayleigh Wave Group Velocity Tomography of Siberia, China and the Vicinity

—Rayleigh waves are used in a tomographic inversion to obtain group velocity maps of East Asia (40° E–160° E and 20° N–70° N). The period range studied is 30 to 70 seconds. Seismograms used for this study were recorded at CDSN stations, at a temporary broadband seismic array in Tibet, at several SRO stations, and Kirnos-equipped stations established in Asia by the former Soviet Union, in Siberia, in the Sakhalin and in Mongolia. Altogether more than 1200 paths were available in the tomographic inversion. The study area includes the Angara craton, the geologically ancient core of Asia, and the subsequently accreted units, the Altaids (a Paleozoic collision complex), the Sino-Korean platform (a chain of Archaen terranes separated by belts of active structures), the south China platform (a collage of Precambrian, Paleozoic and Mesozoic metamorphic and igneous terranes), as well as the Tibetan plateau (an active tectonic feature created in late Cenozoic through collision of the Indian subcontinent and the Asian continent). Many of these main units are recognizable in the tomographic images as distinctive units; Tibet appears as a prominent low velocity (about ?15% from the average) structure, with western and central Tibet often appearing as the areas with the lowest velocities, the Central Asian fold-belt, and the Angara craton are consistently high group velocity areas. Some lesser tectonic features are also recognizable. For example, Lake Baikal is seen as a high velocity feature at periods greater than 40 seconds. However, the high group velocity feature does not stop near the southern end of Lake Baikal; it extends south-southwestward across Mongolia. The North China Plain, a part of the platform where extensional tectonics dominate, is an area of high velocities as a result of relatively thin crust. The south China block, the least tectonically active region of China, is generally an area of high velocity. For periods longer than 40 seconds, a NNE trending high group velocity gradient clearly exists in eastern China; the velocities are noticeably higher in the east. From the group velocity maps, average dispersion curves at twelve locations were determined and inverted to obtain velocity structures. Main results of group velocity inversion include: (1) a Tibetan crust of around 60?km thick, with low crustal and upper mantle shear velocities, at 3.3?km/s and 4.2?km/s, respectively; (2) with the Moho constrained at 40–43?km, the Angara craton and the Central Asian foldbelt have a V _S in excess of 4.6?km/s; (3) relatively low shear velocities are obtained for tectonically active areas. In many parts of the study area, where Precambrian basement is exposed, the process in the crust and upper mantle due to recent tectonic activities have modified the crust and upper mantle velocity structures under the Precambrian terranes, they are no longer underlain by high velocity crust and mantle.     

Deep structure and tomographic imaging of strong earthquake source zones

This work generalizes the results of tomographic imaging performed by the authors for epicentral zones. Seismic events in North Africa (the M _w = 5.8 earthquake of 1985 near the town of Constantine), eastern Anatolia (the Erzincan M _w = 6.7 earthquake of 1992), the Lesser and Greater Caucasus (the 1988 Spitak M _w = 6.8 and the 1991 Racha M _w = 7.0 earthquakes), and northern Sakhalin (the 1995 Neftegorsk M _w = 7.1 earthquake) are examined. It is shown how various morphokinematic types of active faults differ in the resulting tomographic images at various depths. A classification of tomographic images of strong earthquake source zones is proposed in accordance with the rank of their generating faults. The sources of the Spitak, Racha, and Erzincan earthquakes are confined to large boundary faults separating tectonic zones. Lower velocity bands are revealed in the tomographic images, and low velocity “pockets” 1–2 km or somewhat more in width penetrating to a depth of up to 15 km are observed near the fault zones. The Constantine and Neftegorsk earthquakes were generated by faults of a lower rank. The source zones of these events are imaged tomographically as narrow gradient zones.     

Optimization of Surface Wave Identification and Measurement

—?Accurate and reliable measurement of surface waves is important to Comprehensive Nuclear-Test-Ban Treaty (CTBT) monitoring because the M _s :m _b discriminant and its regional variants can in many cases unambiguously identify events as earthquakes or explosions. Surface wave processing at the International Data Center (IDC) is designed to be completely automated and is performed using the program Maxsurf. Maxsurf searches for surface wave characteristics in the expected surface wave arrival time window for all continuous long-period and broadband data in the IDC processing stream. The Prototype IDC GSETT3 Reviewed Event Bulletin (REB) now contains a very large and growing data set of surface wave measurements. Users of this data set need to be aware of processing changes and calibration errors in the GSETT3 experimental bulletin. The prototype International Monitoring System (IMS) surface wave detection threshold is approximately one magnitude unit lower than the detection threshold of other global networks that use visual identification of surface waves. Surface wave identification and measurement can be improved through development of regionalized earth models, phase-matched filtering and the use of path corrected spectral magnitudes in place of M _s . Regionalized earth models are developed through tomographic inversion of a very large data set of phase and group velocity dispersion measurements. Discrimination capability can be improved through the use of maximum likelihood magnitudes and maximum likelihood upper bounds.     

Tomographic Imaging of the Seismic Structure Beneath the East Anatolian Plateau, Eastern Turkey

The high level of seismic activity in eastern Turkey is thought to be mainly associated with the continuing collision of the Arabian and Eurasian tectonic plates. The determination of a detailed three-dimensional (3D) structure is crucial for a better understanding of this on-going collision or subduction process; therefore, a body wave tomographic inversion technique was performed on the region. The tomographic inversion used high quality arrival times from earthquakes occurring in the region from 1999 to 2001 recorded by a temporary 29 station broadband IRIS-PASSCAL array operated by research groups from the Universities of Bo?azi?i (Turkey) and Cornell (USA). The data was inverted and consisted of 3,114 P- and 2,298 S-wave arrival times from 252 local events with magnitudes (M _D) ranging from 2.5 to 4.8. The stability and resolution of the results were qualitatively assessed by two synthetic tests: a spike test and checkerboard resolution test and it was found that the models were well resolved for most parts of the imaged domain. The tomographic inversion results reveal significant lateral heterogeneities in the study area to a depth of ~20?km. The P- and S-wave velocity models are consistent with each other and provide evidence for marked heterogeneities in the upper crustal structure beneath eastern Turkey. One of the most important features in the acquired tomographic images is the high velocity anomalies, which are generally parallel to the main tectonic units in the region, existing at shallow depths. This may relate to the existence of ophiolitic units at shallow depths. The other feature is that low velocities are widely dispersed through the 3D structure beneath the region at deeper crustal depths. This feature can be an indicator of the mantle upwelling or support the hypothesis that the Anatolian Plateau is underlain by a partially molten uppermost mantle.     

Short scale heterogeneity in the lowermost mantle: insights from PcP-P and ScS-S data

We compare lateral variations at the base of the mantle as inferred from a global dataset of PcP-P travel time residuals, measured on broadband records, and existing P and S tomographic velocity models, as well as ScS-S travel time data in some selected regions. In many regions, the PcP-P dataset implies short scale lateral variations that are not resolved by global tomographic models, except under eastern Eurasia, where data and models describe a broad region of fast velocity anomalies across which variations appear to be of thermal origin. In other regions, such as central America and southeastern Africa, correlated short scale lateral variations (several hundred kilometers) are observed in PcP and ScS, implying large but not excessive values for the ratio R=∂ ln V_s/∂ ln V_p (∼2.5). On the other hand, in at least two instances, in the heart of the African Plume and on the edge of the Pacific Plume, variations in P and S velocities appear to be incompatible, implying strong lateral gradients across compositionally different domains, possibly also involving topography on the core-mantle boundary. One should be cautious in estimating R at the base of the mantle from global datasets, as different smoothing and sampling of P and S datasets may result in strong biases and meaningless results.     

Shallow Velocity Structure at the Shagan River Test Site in Kazakhstan

—?During 1997 and 1998, twelve chemical explosions were detonated in boreholes at the former Soviet nuclear test site near the Shagan River (STS) in Kazakhstan. The depths of these explosions ranged from 2.5 to 550 m, while the explosive yield varied from 2 to 25 tons. The purpose of these explosions was for closure of the unused boreholes at STS, and each explosion was recorded at local distances by a network of seismometers operated by Los Alamos National Laboratory and the Institute of Geophysics for the National Nuclear Center (NNC). Short-period, fundamental-mode Rayleigh waves (Rg) were generated by these explosions and recorded at the local stations, resultingly the waves exhibited normal dispersion between 0.2 and 3 seconds. Dispersion curves were generated for each propagation path using the Multiple Filter Analysis and Phase Match Filtering techniques. Tomographic maps of Rg group velocity were constructed and show a zone of relatively high velocities for the southwestern (SW) region of the test site and slow propagation for the northeastern (NE) region. For 0.5?sec Rg, the regions are separated by the 2.1?km/sec contour, as propagation in the SW is greater than 2.1?km/sec and less in the NE region. At 1.0 sec period, the 2.3?km/sec contour separates the two regions. Finally, for 1.5 and 2.0 sec, the separation between the two regions is less distinct as velocities in the NE section begin to approach the SW except for a low velocity region (<2.1?km/sec) near the center of the test site. Local geologic structure may explain the different regions as the SW region is composed predominantly of crystalline intrusive rocks, while the NE region consists of alluvium, tuff deposits, and Paleozoic sedimentary rocks. Low velocities are also observed along the Shagan River as it passes through the SW region of the test site for shorter period Rg (0.5–1.0?sec). Iterative, least-squares inversions of the Rg group velocity dispersion curves show shear-wave velocities for the southwestern section that are on average 0.4?km/sec higher than the NE region. At depths greater than 1.5?km the statistical difference between the models is no longer significant. The observed group velocities and different velocity structures correlate with P-wave complexity and with spatial patterns of magnitude residuals observed from nuclear explosions at STS, and may help to evaluate the mechanisms behind those observations.     

P and S reflection and P refraction: An integration for characterising shallow subsurface

Characterization of shallow structures was performed by using different approaches analysing both P- and S-wave seismic data with different resolution. The refraction tomography provided P and S velocity models of the first 80 m, while the reflection seismic processing gives a reasonable stacking velocity field until 300 m depth for both P- and S-wave data. So, we estimated the V_p/V_s ratio and an empirical relationship between the two velocities. We characterised the shallow layers using tomographic velocity models and the deeper layers using seismic images with different resolution. The seismic images were obtained by conventional CMP reflection seismic processing and by a novel multi-refractor imaging technique.     

Studies of mantle structure of U.S.S.R. territory on long-range seismic profiles

Long-range seismic sounding carried out during the last few years on the territory of the U.S.S.R. has shown a basic inhomogeneity of the uppermost mantle, as well as evidence of regularities in the distribution of its seismic parameters. The following data were used: times and apparent velocities of P- and S-waves for investigation of mantle velocities, converted waves for seismic discontinuity model studies and wave attenuation for Q-factor estimation. Strong regularities were distinguished in the distribution of average seismic velocities for the uppermost mantle, in their dependence on the age and type of geostructure and on their position relative to the central part of the continent. Old platforms and the inner part of the continent are marked by velocities under the Mohorovi?i? discontinuity of more than 8.2–8.3 km s^?1, young platforms and outer parts of the continent by 8.0–8.2 km s^?1, and orogenic and rift zones by 7.8–8.0 km s^?1. The difference becomes more pronounced at a depth of about 100–200 km: for the old platform mantle velocities of 8.5–8.6 km s^?1 are typical; beneath the orogenic and rift areas, inversion zones with velocities less than 7.8 km s^?1 are observed.The converted waves show fine inhomogeneities of the crust and uppermost mantle, the presence of many discontinuities with positive and negative changes of velocity, and anisotropy of seismic waves in some of the layers. Wave attenuation allowed the determination of the Q-factor in the mantle. It varied from one region to another but a close relation between Q and P-wave velocity is the main cause of its variation.     

Subcrustal Low Velocity Layers in Central India and their Implications

A 2-D subcrustal velocity model for the central Indian continental lithosphere has been derived by travel time and relative amplitude modeling of a digitally normalized analog seismic record section of the Hirapur-Mandla DSS profile, using a ray-tracing technique. Some prominent wave groups with apparent velocities slightly higher than the Moho reflection phase (P_MP) are identified on the normalized record sections assembled with a reduction velocity of 6 km s⁻¹. We interpret these phases as the wide-angle reflections from subcrustal lithospheric boundaries. Comparison of synthetic seismograms with the observed record section shows that the observed phases cannot be explained either by multiples or by the P-to-S converted phase (P_MS) from the Moho. Subcrustal velocity models either with a velocity increase or with a single low velocity layer (LVL) also do not provide a satisfactory fit. We infer that a subcrustal velocity model with two alternate LVLs (velocity 7.2 km s⁻¹), separated by a 6-km thick high velocity layer (velocity 8.1 km s⁻¹), can satisfy both the observed travel times and amplitudes. The prominent reflection phases are modeled at depths of 49, 51, 57 and 60 km. It is inferred that the subcrustal lithosphere in the central Indian region has a lamellar structure with varying structural and mechanical properties. The alternating LVLs, occurring at relatively shallow depths below Moho, may be associated with the zones of weakness and lower viscosity suggesting continued mobility, with a possible thermal source in the upper mantle. This explains the source of observed high heat flow values in the central Indian region.     

S-Wave Velocities of the Lithosphere–Asthenosphere System in the Caribbean Region

An overview of the S-wave velocity (V _s) structural model of the Caribbean with a resolution of 2°?×?2° is presented. New tomographic maps of Rayleigh wave group velocity dispersion at periods ranging from 10 to 40?s were obtained as a result of the frequency time analysis of seismic signals of more than 400 ray-paths in the region. For each cell of 2°?×?2°, group velocity dispersion curves were determined and extended to 150?s by adding data from a larger scale tomographic study (Vdovin et al., Geophys. J. Int 136:324–340, 1999). Using, as independent a priori information, the available geological and geophysical data of the region, each dispersion curve has been inverted by the “hedgehog” non-linear procedure (Valyus, Determining seismic profiles from a set of observations (in Russian), Vychislitielnaya Seismologiya 4, 3–14. English translation: Computational Seismology (V.I. Keylis-Borok, ed.) 4:114–118, 1968), in order to compute a set of V _s versus depth models up to 300?km of depth. Because of the non-uniqueness of the solutions for each cell, a local smoothness optimization has been applied to the whole region in order to choose a three-dimensional model of V _s, satisfying this way the Occam's razor concept. Several known and some new main features of the Caribbean lithosphere and asthenosphere are shown on these models such as: the west directed subduction zone of the eastern Caribbean region with a clear mantle wedge between the Caribbean lithosphere and the subducted slab; the complex and asymmetric behavior of the crustal and lithospheric thickness in the Cayman ridge; the predominant oceanic crust in the region; the presence of continental type crust in Central America, and the South and North America plates; as well as the fact that the bottom of the upper asthenosphere gets shallower going from west to east.     

3-D surface wave group velocity distribution in Central-Southern Africa

Group velocities estimated from fundamental mode Love and Rayleigh waves are used in a tomography process in central-southern Africa. The waves were generated by eighteen earthquakes, which occurred along the East African Rift and recorded at BOSA, LBTB and SLR seismic stations in southern Africa. The group velocities from Love and Rayleigh waves were isolated using the Multiple Filter Technique (MFT) at the period range of 10 to 50 seconds. The tomography method developed by Ditmar and Yanovskaya (1987) and Yanovskaya and Ditmar (1990), was applied to calculate the lateral distribution of surface wave group velocities in central-southern Africa. The results of the tomographic inversion were plotted as distribution maps. In addition to the maps, I also produced two velocity cross-sections across the area of study. The velocity distribution maps show the regional tectonic units, though with poor resolution. The azimuthal bias of the surface wave paths is reflected in the distribution of the group velocities. The Moho depth appears to correlate with velocities at a period of about 30 s. A low velocity feature observed beneath the Zimbabwe craton implies a thickening upper asthenosphere and lithospheric thinning beneath the Zimbabwe craton. Also estimated was a shear wave velocity model beneath the Zimbabwe craton.     

Global Models of Surface Wave Group Velocity

—?Measurements of group velocity are derived from phase-velocity dispersion curves and modeled with global laterally-varying isotropic structure. Maps for both Love and Rayleigh waves are created in the period range 35?s to 175?s. The data set of group-velocity measurements includes over 50,000 minor-arc observations and 5,000 major-arc observations. The errors in the measurements are estimated by an empirical method of comparing pairwise-similar paths, resulting in uncertainties which are 20% to 40% of the size of the typical measurement. The models are determined by least-squares inversion for spherical harmonic maps expanded up to degree 40. This parameterization allows for resolution of structures as small as 500?km. The models explain 70–98% of the variance relative to the Preliminary Reference Earth Model (PREM). For the area of Eurasia, the group-velocity maps from this study are compared with those of Ritzwoller and Levshin (1998). The results of the two studies are in very good agreement, particularly in terms of spatial correlation. The models also agree in amplitude at wavelengths longer than 30?degrees. For shorter wavelengths, the agreement is good only for models at short periods. The global maps are useful for prediction of group arrival times, for revealing tectonic structures, for determination of seismic event locations and source parameters, and as a basis for regional group-velocity studies.     

Surface Wave Velocities Across Arabia

— The group-velocity distribution beneath the Arabian Plate is investigated using Love and Rayleigh waves. We obtained a balanced path coverage using seismograms generated by earthquakes located along the plate boundaries. We measured Love- and Rayleigh-wave group-velocity dispersion using multiple filter analysis and then performed a tomographic inversion using these observations to estimate lateral group velocity variations in the period range of 5–60?s. The Love- and Rayleigh-wave results are consistent and show that the average group velocity across Arabia increases with increasing period. The tomographic results also delineate first-order regional structure heterogeneity as well as the sharp transition between the Arabian shield and the Arabian platform. Systematic differences are observed in the distribution of the short-period group velocities across the two provinces, which are consistent with surface geology. The slower velocities in the platform reveal the imprint of its thick sedimentary section, while faster velocities correlate well with the exposed volcanic flows in the shield. Shear-wave velocity models for the two regions, obtained from the inversion of the group velocities, confirm results from previous studies of higher S-wave velocity in the upper crust beneath the shield. This may be due to the present remnants of the oceanic crust (ophiolite belts) associated with the island arcs evolutionary model of the Arabian shield.¶The mapping of the surface-wave group velocity using a large data can be used in constraining the regional structure at existing and planned broadband stations deployed in this tectonically complex region as part of the seismic monitoring under CTBT.     

Regional Magnitude Scaling, Transportability, and Ms:mb Discrimination at Small Magnitudes

—?Data sets of m _b (Pn) and m _b (Lg) measurements are presented for three continental regions in order to investigate scaling relationships with moment magnitude M _w and event discrimination at small magnitudes. Compilations of published measurements are provided for eastern North American and central Asian earthquakes, and new measurements are reported for earthquakes located in western United States. Statistical tests on M _w :m _b relationships show that the m _b (Lg) scale of Nuttli (1973) is transportable between tectonic regions, and a single, unified M _w :m _b (Lg) relationship satisfies observations for M _w ～4.2–6.5 in all regions. A unified relationship is also developed for nuclear explosions detonated at the Nevada Test Site and test sites of the former Soviet Union. Regional m _b for explosions scale at higher rates than for earthquakes, and of significance is the finding that m _b (Pn) for explosions scales at a higher rate than m _b (Lg). A model is proposed where differences in scaling rates are related to effects of spectral overshoot and near-field Rg scattering on the generation of Pn and Lg waves by explosions. For earthquakes, m _b (Pn) and m _b (Lg) scale similarly, showing rates near 1.0 or 2/3?·?log₁₀ M _o (seismic moment).¶M _w :m _b (Lg) scaling results are converted to unified M _s :m _b (Lg) relationships using scaling laws between log M _o and M _s . For earthquakes with M _s greater than 3.0, the scaling rate is 0.69?·?M _s , which is the same as it is for nuclear explosions if M _s is proportional to 1.12?·?log M _o, as determined by NTS observations. Thus, earthquake and explosion populations are parallel and separated by 0.68 m _b units for large events. For small events (M _s ?M _s :m _b (Lg) plots for stable and tectonic regions, respectively. While the scaling rate for explosions is ～0.69, this value is uncertain due to paucity of M _o observations at small yields. Measurements of [m _b (P)???m _b (Lg)] for earthquakes in the western United States have an average value of ?0.33?±?.03 m _b units, in good agreement with Nuttli's estimate of m _b bias for NTS. This result suggests that Nuttli's method for estimating test site bias can be extended to earthquakes to make estimates of bias on regional scales. In addition, a new approach for quick assessments of regional bias is proposed where M _s :m _b (P) observations are compared with M _s :m _b (Lg) relationships. Catalog M _s :m _b (P) data suggest that m _b bias is significant for tectonic regions of southern Asia, averaging about ?0.4 m _b units.     

Application of 3-D Crustal and Upper Mantle Velocity Model of North America for Location of Regional Seismic Events

—?Seismic event locations based on regional 1-D velocity-depth sections can have bias errors caused by travel-time variations within different tectonic provinces and due to ray-paths crossing boundaries between tectonic provinces with different crustal and upper mantle velocity structures. Seismic event locations based on 3-D velocity models have the potential to overcome these limitations. This paper summarizes preliminary results for calibration of IMS for North America using 3-D velocity model. A 3-D modeling software was used to compute Source-Station Specific Corrections (SSSCs(3-D)) for Pn travel times utilizing 3-D crustal and upper mantle velocity model for the region. This research was performed within the framework of the United States/Russian Federation Joint Program of Seismic Calibration of the International Monitoring System (IMS) in Northern Eurasia and North America.¶An initial 3-D velocity model for North America was derived by combining and interpolating 1-D velocity-depth sections for different tectonic units. In areas where no information on 1-D velocity-depth sections was available, tectonic regionalization was used to extrapolate or interpolate. A Moho depth map was integrated. This approach combines the information obtained from refraction profiles with information derived from local and regional network data. The initial 3-D velocity model was tested against maps of Pn travel-time residuals for eight calibration explosions; corrections to the 3-D model were made to fit the observed residuals. Our goal was to find a 3-D crustal and upper mantle velocity model capable predicting Pn travel times with an accuracy of 1.0–1.5 seconds (r.m.s.).¶The 3-D velocity model for North America that gave the best fit to the observed travel times, was used to produce maps of SSSCs(3-D) for seismic stations. The computed SSSCs(3-D) vary approximately from +5 seconds to ?5 seconds for the western USA and the Pre-Cambrian platform, respectively. These SSSCs(3-D) along with estimated modeling and measurement errors were used to relocate, using regional data, an independent set of large chemical explosions (with known locations and origin times) detonated within various tectonic provinces of North America. Utilization of the 3-D velocity model through application of the computed SSSCs(3-D) resulted in a substantial improvement in seismic event location accuracy and in a significant decrease of error ellipse area for all events analyzed in comparison both with locations based on the IASPEI91 travel times and locations based on 1-D regional velocity models.     

Rayleigh wave group velocity tomography of Gujarat region,Western India and its implications to mantle dynamics

In the present study, fundamental Rayleigh waves with varying period from 10 to 80 s are used to obtain group velocity maps in the northwest Deccan Volcanic Province of India. About 350 paths are obtained using 53 earthquakes (4.8 ≤ M ≥ 7.9) recorded by the SeisNetG (Seismic Network of Gujarat). Individual dispersion curves of group velocity of Rayleigh wave for each source-station path are estimated using multiple filter technique. These curves are used to determine lateral distribution of Rayleigh wave group velocity by tomographic inversion method. Our estimated Rayleigh group velocity at varying depths showed conspicuous corroboration with three tectonic blocks [Kachchh Rift Basin (KRB), Saurashtra Horst (SH), and Mainland Gujarat (MG)] in the region. The seismically active KRB with a thicker crust is characterized as a low velocity zone at a period varying from 10 to 30 s as indicative of mantle downwarping or sagging of the mantle beneath the KRB, while the SH and MG are found to be associated with higher group velocities, indicating the existence of the reduced crustal thickness. The trend of higher group velocity was found prevailed adjacent to the Narmada and Cambay rift basins that also correspond to the reduced crust, suggesting the processes of mantle upwarping or uplifting due to mantle upwelling. The low velocities at periods longer than 40 s beneath the KRB indicate thicker lithosphere. The known Moho depth correlates well with the observed velocities at a period of about 30 s in the Gujarat region. Our estimates of relatively lower group velocities at periods varying from 70 to 80 s may correspond to the asthenospheric flow beneath the region. It is interesting to image higher group velocity for the thinner crust beneath the Arabian Sea adjacent to the west coast of Gujarat at the period of 40 s that may correspond to the upwarped or upwelled mantle beneath the Arabian Sea. Our results have better resolution estimated by a radius of equivalent circular averaging area for each period.     

Great-circle Rayleigh wave attenuation and group velocity,Part IV: Regionalization and pure-path models for shear velocity and attenuation

Pure-path averages for group velocities and specific attenuation have been calculated from individual observations and from path averages for two regionalizations; one original to this study and the other previously devised by Wu. Both are based on four upper-mantle provinces: ocean basin, continent, island arc and mid-ocean ridge. Pure-path group velocities and specific attenuation have also been calculated for combinations of regions and provide well separated regional measurements for such composite regions.Shear-velocity models for pure and combined regions have been derived by a controlled Monte Carlo inversion procedure and indicates that a low-velocity zone is required beneath the oceans, but is not required beneath continents. Models have been produced for pure and combined ocean, ocean-ridge, continent and continent-arc provinces.Q^?1_R determined from pure-path average group velocities and attenuation coefficients has been regionalized successfully for 2- and 3-region combinations. The resulting pure-path Q^?1_R for continents is much lower than that for ocean basins and ocean-ridge provinces. Inversion of Q^?1_R for ocean-ridge provinces shows that the average Q_β for the upper 200 km of these regions is between 85 and 100.     

Velocities,elastic moduli and weathering-age relations for pacific layer 2 basalts

Compressional (V_p) and shear (V_s) wave velocities have been measured to 10 kb in 32 cores of basalt from 14 Pacific sites of the Deep Sea Drilling Project. Both V_pandV_s show wide ranges (3.70to6.38km/sec forV_pand1.77to3.40km/sec forV_sat0.5kb) which are linearly related to density and sea floor age, confirming earlier findings by Christensen and Salisbury of decreasing velocity with progressive submarine weathering based on studies of basalts from five sites in the Atlantic. Combined Pacific and Atlantic data give rates of decreasing velocity of ?1.89and?1.35km/sec per100my forV_pandV_s respectively. New analyses of oceanic seismic refraction data indicate a decrease in layer 2 velocities with age similar to that observed in the laboratory, suggesting that weathering penetrates to several hundred meters in many regions and is largely responsible for the extreme range and variability of layer 2 refraction velocities.     

Theoretical and Observed Depth Correction for M s

—?Modal summation technique is used to generate 5000, three-component theoretical seismograms of Love and Rayleigh waves, assuming modified PREM (PREM-C) and AK135F global earth models. The focal depth h and the geometrical fault parameters are randomly chosen so as to uniformly cover possible source mechanisms and obtain uniform distribution of log h in the interval 1?h?h?M _s of the form:¶ΔM _s (h)=0 forh< 20km, ΔM _s (h)=0.314log(h)-0.409 for 20≠h< 60km, ΔM _s (h)=1.351log(h)-2.253 for 60≠h< 100km, ΔM _s (h)=0.400log(h)-0.350 for 100≠h< 600km .¶After applying the above correction, the relationship between the surface wave magnitude and the scalar seismic moment for the observational data set significantly improves, and becomes independent of the source depth. In relation to CTBT, no depth correction is needed for M _S when the m _b ???M _S discriminant is computed, because the proposed correction is zero for earthquakes with foci above 20?km.