Q β structure of the crust and upper mantle in the eastern Sino-Korean paraplatform

Based on the long period surface wave data recorded by the China Digital Seismograph Network (CDSN), theQ _R of fundamental mode Rayleigh wave with periods from 10 s to 146 s is determined for the eastern Sino-Korean paraplatform in this paper. TheQ _β models of the crust and upper mantle are respectively obtained for the 4 paths, with the aid of stochastic inverse method. It shows that in the eastern Sino-Korean paraplatform, the average crustalQ _β is about 200, and that there exists a weak attenuation layer in the middle crust (about 10–20 km deep) which is possibly related to earthquake-prone layer. A strong attenuation layer (lowQ) of 70 km thick extensively exists in the uppermost mantle, with the buried depth about 80 km. The averageQ _R of fundamental mode Rayleigh wave is between the value of stable tectonic region and that of active tectonic region, and much close to the latter.     

Great-circle Rayleigh wave attenuation and group velocity,Part III: Inversion of global average group velocities and attenuation coefficients

Average shear-velocity models for the upper mantle have been derived by controlled Monte Carlo inversion of global average Rayleigh wave group velocity (GAGV) data for periods between 50 and 300 seconds. GAGV data have been corrected for attenuative dispersion using a method based on the theory of Liu, Anderson and Kanamori. Two types of model bounds have been used with one- or two-layer low-velocity zones beginning at depths of 70 and 100 km. All models fitting GAGV data within one standard deviation have low-velocity zones in the 100–200 km depth range. Models with low-velocity zones beginning at 70 km, as well as 100 km, fit GAGV data within one standard deviation, so the average thickness of the lithosphere (taken as the depth to the top of the low-velocity zone) cannot be determined with precision.Global average models for shear-wave attenuation (Q^?1_β) have been derived from global average Rayleigh wave attenuation coefficients for periods between 50 and 300 s and average shear-velocity models. Zones of high Q^?1_β coincide with the low-velocity zones of all shear-velocity models, however, models with low-velocity zones beginning at a depth of 70 km have the highest-attenuation layer in the lower half of the low-velocity zone. Resolution kernels for these attenuation models show that parameters for layers shallower than the lower part of the low-velocity-high-attenuation zone are strongly coupled but are distinct from the lower part of this zone. This suggests that the deeper part of the low-velocity-high-attenuation zone is the most mobile part of the zone or that on the average, the top of the zone is deeper than 70 km.The average Q_β of the lithosphere, low-velocity zone, and sub-low-velocity layer (asthenosphere) are approximately 200, 85–110 and 170–200, respectively.     

Tomography and velocity structure of the crust and uppermost mantle in southeastern Europe obtained from surface wave analysis

A set of two hundred shear-wave velocity models of the crust and uppermost mantle in southeast Europe is determined by application of a sequence of methods for surface-waves analysis. Group velocities for about 350 paths have been obtained after analysis of more than 600 broadband waveform records. Two-dimensional surface-wave tomography is applied to the group-velocity measurements at selected periods and after regionalisation, two sets of local dispersion curves (for Rayleigh and Love waves) are constructed in the period range 8–40 s. The shear-wave velocity models are derived by applying non-linear iterative inversion of local dispersion curves for grid cells predetermined by the resolving power of data. The period range of observations limits the velocity models to depths of 70 km in accordance to the penetration of the surface waves with a maximum period of 40 s. Maps of the Moho boundary depth, velocity distribution above and below Moho boundary, as well as velocity distribution at different depths are constructed. Well-known geomorphologic units (e.g. the Pannonian basin, southeastern Carpathians, Dinarides, Hellenides, Rodophean massif, Aegean Sea, western Turkey) are delineated in the obtained models. Specific patterns in the velocity models characterise the southeast Carpathians and adjacent areas, coast of Albania, Adriatic coast of southern Italy and the southern coast of the Black Sea. The models obtained in this study for the western Black Sea basin shows the presence of layers with shear-wave velocities of 3.5 km/s–3.7 km/s in the crust and thus do not support the hypothesis of existence of oceanic structure in this region.     

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.     

Q-structure beneath the Tibetan Plateau from the inversion of Love- and Rayleigh-wave attenuation data

The fundamental mode Love and Rayleigh waves generated by ten earthquakes and recorded across the Tibet Plateau, at QUE, LAH, NDI, NIL, KBL, SHL, CHG, SNG and HKG are analysed. Love- and Rayleigh-wave attenuation coefficients are obtained at time periods of 5–120 s using the spectral amplitudes of these waves for 23 different paths. Love wave attenuation coefficient varies from 0.0021 km^?1, at a period of 10 s, to 0.0002 km^?1 at a period of 90 s, attaining two maxima at time periods of 10 and 115 s, and two minima at time periods of 25 and 90 s. The Rayleigh-wave attenuation coefficient also shows a similar trend. The very low value for the dissipation factor, Q_β, obtained in this study suggests high dissipation across the Tibetan paths. Backus-Gilbert inversion theory is applied to these surface wave attenuation data to obtain average Q_β^?1 models for the crust and uppermost mantle beneath the Tibetan Plateau. Independent inversion of Love- and Rayleigh-wave attenuation data shows very high attenuation at a depth of ～50–120 km ($\text{Q}{}_{\mathrm{\beta }}\text{? 10}$). The simultaneous inversion of the Love and Rayleigh wave data yields a model which includes alternating regions of high and low Q_β^?1 values. This model also shows a zone of high attenuating material at a depth of ～40–120 km. The very high inferred attenuation at a depth of ～40–120 km supports the hypothesis that the Tibetan Plateau was formed by horizontal compression, and that thickening occurred after the collision of the Indian and Eurasian plates.     

Crustal and upper mantle structures of the brazilian coast

The Rayleigh wave phase and group velocities in the period range of 24–39 sec, obtained from two earthquakes which occurred in northeastern brazil and which were recorded by the Brazilian seismological station RDJ (Rio de Janeiro), have been used to study crustal and upper mantle structures of the Brazilian coastal region. Three crustal and upper mantle models have been tried out to explain crustal and upper mantle structures of the region. The upper crust has not been resolved, due basically to the narrow period range of the phase and group velocities data. The phase velocity inversions have exhibited good resolutions for both lower crust and upper mantle, with shear wave velocities characteristic of these regions. The group velocity data inversions for these models have showed good results only for the lower crust. The shear wave velocities of the lower crust (3.86 and 3.89 km/sec), obtained with phase velocity inversions, are similar to that (=3.89 km/sec) found byHwang (1985) to the eastern South American region, while group velocity inversions have presented shear velocity (=3.75 km/sec) similar to that (=3.78 km/sec) found byLazcano (1972) to the Brazilian shield. It was not possible to define sharply the crust-mantle transition, but an analysis of the phase and group velocity inversions results has indicated that the total thickness of the crust should be between 30 and 39 km. The crustal and upper mantle model, obtained with phase velocity inversion, can be used as a preliminary model for the Brazilian coast.     

Shear wave velocity models retrieved using Rg wave dispersion data in shallow crust in some regions of southern Ontario, Canada

Many crucial tasks in seismology, such as locating seismic events and estimating focal mechanisms, need crustal velocity models. The velocity models of shallow structures are particularly important in the simulation of ground motions. In southern Ontario, Canada, many small shallow earthquakes occur, generating high-frequency Rayleigh (Rg) waves that are sensitive to shallow structures. In this research, the dispersion of Rg waves was used to obtain shear-wave velocities in the top few kilometers of the crust in the Georgian Bay, Sudbury, and Thunder Bay areas of southern Ontario. Several shallow velocity models were obtained based on the dispersion of recorded Rg waves. The Rg waves generated by an m _N 3.0 natural earthquake on the northern shore of Georgian Bay were used to obtain velocity models for the area of an earthquake swarm in 2007. The Rg waves generated by a mining induced event in the Sudbury area in 2005 were used to retrieve velocity models between Georgian Bay and the Ottawa River. The Rg waves generated by the largest event in a natural earthquake swarm near Thunder Bay in 2008 were used to obtain a velocity model in that swarm area. The basic feature of all the investigated models is that there is a top low-velocity layer with a thickness of about 0.5 km. The seismic velocities changed mainly within the top 2 km, where small earthquakes often occur.     

利用ML振幅研究地壳横波Q值Ⅱ:Q横向变化特征

利用中国地震年报(1984～2004年)测定M_L所用的54000余条区域地震最大振幅数据,采用二维层析成像方法反演中国大陆地壳介质品质因子Q₀结构.反演结果表明,中国大陆地壳介质品质因子Q₀平均为370,Q₀横向变化量可达±200.在整个研究区解的分辨率能够达到3°×3°,大部分地区能够达到2°×2°.Q₀的变化显示出与地质构造分区相关的特征：低衰减 (高Q₀值) 异常区主要分布在塔里木盆地、四川盆地、鄂尔多斯高原及华南地区,鄂尔多斯高原北部Q₀值最高;高衰减异常区(低Q₀值)主要分布在青藏高原以东的川滇地区、华北平原、渤海湾,最低Q₀值在渤海湾地区,天山、祁连山和汾渭地堑Q₀也呈较低异常变化.总体看来,构造上稳定的克拉通地台的介质品质因子Q₀较高,而构造上活动的地区往往Q₀较低,显现出区域的衰减变化特征与地壳结构和区域构造之间的密切关系.     

Anisotropy in the Earth's crust and uppermost mantle in southeastern Europe obtained from Rayleigh and Love surface waves

Digital seismograms from 25 earthquakes located in the southeastern part of Europe, recorded by three-component very broadband seismometers at the stations Vitosha (Bulgaria) and Muntele Rosu (Romania), were processed to obtain the dispersion properties of Rayleigh and Love surface waves. Rayleigh and Love group-velocity dispersion curves were obtained by frequency–time analysis (FTAN). The path-averaged shear-wave velocity models were computed from the obtained dispersion curves. The inversion of the dispersion curves was performed using an approach based on the Backus–Gilbert inversion method. Finally, 70 path-averaged velocity models (35 R-models computed from Rayleigh dispersion curves and 35 L-models computed from Love dispersion curves) were obtained for southeastern Europe. For most of the paths, the comparison between each pair of models (R-model and L-models for the same path) shows that for almost all layers the shear-wave velocities in the L-models are higher than in the R-models. The upper sedimentary layers are the only exception. The analysis of both models shows that the depth of the Moho boundary in the L-models is shallower than its depth in the R-models. The existence of an anisotropic layer associated with the Moho boundary at depths of 30–45 km may explain this phenomenon. The anisotropy coefficient was calculated as the relative velocity difference between both R- and L-models at the same depths. The value of this coefficient varies between 0% and 20%. Generally, the anisotropy of the medium caused by the polarization anisotropy is up to 10–12%, so the maximum observed discrepancies between both types of models are also due to the lateral heterogeneity of the shear-wave velocity structure of the crust and the upper mantle in the region.     

Crustal seismic structure and depth distribution of earthquakes in the Archean Kuusamo region,Fennoscandian Shield

Two-dimensional crustal velocity models are derived from passive seismic observations for the Archean Karelian bedrock of north-eastern Finland. In addition, an updated Moho depth map is constructed by integrating the results of this study with previous data sets. The structural models image a typical three-layer Archean crust, with thickness varying between 40 and 52 km. P wave velocities within the 12–20 km thick upper crust range from 6.1 to 6.4 km/s. The relatively high velocities are related to layered mafic intrusive and volcanic rocks. The middle crust is a fairly homogeneous layer associated with velocities of 6.5–6.8 km/s. The boundary between middle and lower crust is located at depths between 28 and 38 km. The thickness of the lower crust increases from 5–15 km in the Archean part to 15–22 km in the Archean–Proterozoic transition zone. In the lower crust and uppermost mantle, P wave velocities vary between 6.9–7.3 km/s and 7.9–8.2 km/s. The average V_p/V_s ratio increases from 1.71 in the upper crust to 1.76 in the lower crust.The crust attains its maximum thickness in the south-east, where the Archean crust is both over- and underthrust by the Proterozoic crust. A crustal depression bulging out from that zone to the N–NE towards Kuusamo is linked to a collision between major Archean blocks. Further north, crustal thickening under the Salla and Kittilä greenstone belts is tentatively associated with a NW–SE-oriented collision zone or major shear zone. Elevated Moho beneath the Pudasjärvi block is primarily explained with rift-related extension and crustal thinning at ∼2.4–2.1 Ga.The new crustal velocity models and synthetic waveform modelling are used to outline the thickness of the seismogenic layer beneath the temporary Kuusamo seismic network. Lack of seismic activity within the mafic high-velocity body in the uppermost 8 km of crust and relative abundance of mid-crustal, i.e., 14–30 km deep earthquakes are characteristic features of the Kuusamo seismicity. The upper limit of seismicity is attributed to the excess of strong mafic material in the uppermost crust. Comparison with the rheological profiles of the lithosphere, calculated at nearby locations, indicates that the base of the seismogenic layer correlates best with the onset of brittle to ductile transition at about 30 km depth.We found no evidence on microearthquake activity in the lower crust beneath the Archean Karelian craton. However, a data set of relatively well-constrained events extracted from the regional earthquake catalogue implies a deeper cut-off depth for earthquakes in the Norrbotten tectonic province of northern Sweden.     

Structure of the crust and upper mantle beneath Australia from Rayleigh- and Love-wave observations

Fundamental and first higher modes of the Rayleigh- and Love-wave group velocities along seven paths in Australia were jointly inverted by a controlled Monte Carlo procedure to obtain regional shear-wave velocity structures of the crust and upper mantle. Our data support the results of Gonez and Cleary which show an S-wave low velocity zone centred near 110 km depth in eastern Australia. However, the thickness-velocity contrast of the low velocity zone is significantly smaller. The crustal models for eastern Australia are characterized by upper crusts which are both thicker and have lower velocities than those in western Australia and have a less sharp crust-upper mantle boundary. The S-wave velocities for the upper mantle appear to be similar (～ 4.55 km s^?1) throughout the continent, with no obvious dependence on the age of cratonization or crustal thickness.     

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.     

Realistic Models of Anisotropic Laminated Lower Crust

—The genesis of the laminated lower crust has been attributed to extensional processes leading to structural and textural ordering. This implies that the lower crust might be anisotropic. Laboratory measurements of lower crustal rock samples and xenolithes show evidence of anisotropy in these rocks due to oriented structure.¶In this paper we investigate the seismic shear-wave response of realistic anisotropic lower crustal models using the anisotropic reflectivity method. Our models are based on representative petrophysical data obtained from exposed lower crustal sections in Calabria (South Italy), Val Strona and Val Sesia (Ivrea Zone, Northern Italy). The models consist of stacks of anisotropic layers characterized by quantified elastic tensors derived from representative rock samples which provide alternating high and low velocity layers.¶The seismic signature of the data is comparable to seismic observations of in situ lower crust. For the models based on the Calabria and Val Strona sequences shear-wave splitting occurs for the Moho reflection at offsets beyond 70 km with travel-time delays up to 300 and 500 ms, respectively. The leading shear wave is predominantly horizontally polarized and followed by a predominantly vertically polarized shear wave. Contrastingly, the Val Sesia model shows no clear evidence of birefringence. Isotropic versus anisotropic modelling demonstrates that the shear-wave splitting is clearly related to the intrinsic anisotropy of the lower crustal rocks for the Val Strona sequence. No evidence of birefringence caused by thin layering is found.     

Great-circle Rayleigh wave attenuation and group velocity,Part I: Observations for periods between 150 and 600 seconds for seven great-circle paths

Rayleigh wave attenuation coefficients and group velocities have been estimated for seven great-circle paths. The attenuation coefficient measurements cover the period range from 100 to 500 s, and group velocities the range from 100 to 600 s. Global average group velocities and attenuation coefficients have also been estimated for these period ranges. The spread of the individual path group velocities for 20-s averaging windows centred at 290, 250, 210, 180 and 150 s is less than 0.034, 0.028, 0.024, 0.048 and 0.071 km/s, respectively. Global average attenuation coefficients, when combined with global average group velocities, show that Q for Rayleigh waves has an approximately constant value of about 145 for periods between 150 and 220 s and slowly increases to a value of about 200 at a period of 400 s.     

A Summary of Attenuation Measurements from Borehole Recordings of Earthquakes: The 10 Hz Transition Problem

—Earthquake seismograms recorded by instruments in deep boreholes have low levels of background noise and wide signal bandwidth. They have been used to extend our knowledge of crustal attenuation both in the near-surface and at seismogenic depths. Site effects are of major importance to seismic hazard estimation, and the comparison of surface, shallow and deep recordings allows direct determination of the attenuation in the near-surface. All studies to date have found that Q is very low in the near-surface (～ 10 in the upper 100 m), and increases rapidly with depth. Unlike site amplification, attenuation at shallow depths exhibits little dependence on rock-type. These observations are consistent with the opening of fractures under decreasing lithostatic pressure being the principal cause of the severe near-surface attenuation. Seismograms recorded in deep boreholes are relatively unaffected by near-surface effects, and thus can be used to measure crustal attenuation to higher frequencies (≥ 100 Hz) than surface recordings. Studies using both direct and coda waves recorded at over 2 km depth find Q to be high (～ 1000) at seismogenic depths in California, increasing only weakly with frequency between 10 and 100 Hz. Intrinsic attenuation appears to be the dominant mechanism. These observations contrast with those of the rapidly increasing Q with frequency determined from surface studies in the frequency range 1 to 10 Hz. Further work is necessary to constrain the factors responsible for this apparent change in the frequency dependence of Q, but it is clearly unwise to extrapolate Q estimates made below about 10 Hz to higher frequencies.     

Frequency dependence of crustalQ β in stable and tectonically active regions

Fundamental-mode Rayleigh wave attenuation data for stable and tectonically active regions of North America, South America, and India are inverted to obtain several frequency-independent and frequency-dependentQ models. Because of trade-offs between the effect of depth distribution and frequency-dependence ofQ on surface wave attenuation there are many diverse models which will satisfy the fundamental-mode data. Higher-mode data, such as 1-Hz Lg can, however, constrain the range of possible models, at least in the upper crust. By using synthetic Lg seismograms to compute expected Lg attenuation coefficients for various models we obtained frequency-dependentQ models for three stable and three tectonically active regions, after making assumptions concerning the nature of the variation ofQ with frequency.In stable regions, ifQ varies as , where is a constant, models in which =0.5, 0.5, and 0.75 satisfy fundamental-mode Rayleigh and 1-Hz Lg data for eastern North America, eastern South America, and the Indian Shield, respectively. IfQ is assumed to be independent of frequency (=0.0) for periods of 3 s and greater, and is allowed to increase from 0.0 at 3 s to a maximum value at 1 s, then that maximum value for is about 0.7, 0.6, and 0.9, respectively, for eastern North America, eastern South America, and the Indian Shield. TheQ models obtained under each of the above-mentioned two assumptions differ substantially from one another for each region, a result which indicates the importance of obtaining high-quality higher-mode attenuation data over a broad range of periods.Tectonically active regions require a much lower degree of frequency dependence to explain both observed fundamental-mode and observed Lg data. Optimum values of for western North America and western South America are 0.0 if is constant (Q is independent of frequency), but uncertainty in the Lg attenuation data allows to be as high as about 0.3 for western North America and 0.2 for western South America. In the Himalaya, the optimum value of is about 0.2, but it could range between 0.0 and 0.5. Frequency-independent models (=0.0) for these regions yield minimumQ values in the upper mantle of about 40, 70, and 40 for western North America, western South America, and the Himalaya, respectively.In order to be compatible with the frequency dependence ofQ observed in body-wave studies,Q in stable regions must be frequency-dependent to much greater depths than those which can be studied using the surface wave data available for this study, andQ in tectonically active regions must become frequency-dependent at upper mantle or lower crustal depths.On leave from the Department of Geophysics, Yunnan University, Kunming Yunnan, People's Republic of China     

Structure of a Volcanic Continental Margin Derived from Ocean-bottom Seismographic Data: The Northern Vøring Margin, off Norway

—This paper presents a crustal model derived from an Ocean Bottom Seismograph (OBS) study along the northern Vøring margin off Norway. The profile was acquired to map the crustal structure in the northernmost part of the Vøring Basin, and to link crustal models of the Lofoten and central Vøring Basin obtained by previous OBS studies. The Vøring margin, as well as the Lofoten margin to the north, was created by continental breakup between Norway and Greenland in late Paleocene-early Eocene. The rifting and continental breakup process were accompanied by intense extrusive and intrusive magmatic activities. The OBS data provide the whole crustal structure along the northern Vøring margin, in the area where the deep crustal structure cannot be resolved by conventional multichannel reflection data due to sill intrusions in the sedimentary sequence. The shallow part of the crustal model is characterized by up to 10 km thick sediments, a sequence of flood basalts and sill intrusions. The P-wave velocities in the flood basalts and sill intrusions are estimated to 5.0 km/s and 4.7–5.8 km/s, respectively. The model indicates an abrupt thickening of the upper crystalline crust from approx.3 km in the NE, to about 10 km towards the SE, with velocities of 6.0–6.2 km/s. The lower crustal velocities are not well resolved due to lack of clear refraction arrivals from the lower crust. However, the observed amplitude versus offsets are best explained by a model with a change in lower crustal velocities from 6.8 to 7.2 km/s beneath the Bivrost lineament. The modelling infers the presence of a lower crustal reflector beneath the lineament, which represents the landward continuation of the Bivrost lineament. Reflection arrivals from the Moho reveal a Moho depth of 23 km in the middle of the profile and 18– 20 km in the northeastern part of the profile. A 370 km long crustal section from the central part of the Vøring Basin to the Lofoten margin, obtained by the results of this study and previous OBS studies, shows a simple thinned continental crust on the Lofoten margin, and a high velocity lower crust underlying an upper crust of varying thickness in the Vøring Basin. The transition between these structures is situated beneath the Bivrost lineament in the lower crust, and beneath the basement high about 40 km south of the lineament in the upper crust.     

Mapping Crustal Heterogeneity Using Lg Propagation Efficiency Throughout the Middle East, Mediterranean, Southern Europe and Northern Africa

—?In this paper we describe a technique for mapping the lateral variation of Lg characteristics such as Lg blockage, efficient Lg propagation, and regions of very high attenuation in the Middle East, North Africa, Europe and the Mediterranean regions. Lg is used in a variety of seismological applications from magnitude estimation to identification of nuclear explosions for monitoring compliance with the Comprehensive Nuclear-Test-Ban Treaty (CTBT). These applications can give significantly biased results if the Lg phase is reduced or blocked by discontinuous structure or thin crust. Mapping these structures using quantitative techniques for determining Lg amplitude attenuation can break down when the phase is below background noise. In such cases Lg blockage and inefficient propagation zones are often mapped out by hand. With our approach, we attempt to visually simplify this information by imaging crustal structure anomalies that significantly diminish the amplitude of Lg. The visualization of such anomalies is achieved by defining a grid of cells that covers the entire region of interest. We trace Lg rays for each event/station pair, which is simply the great circle path, and attribute to each cell a value equal to the maximum value of the Lg/P-coda amplitude ratio for all paths traversing that particular cell. The resulting map, from this empirical approach, is easily interpreted in terms of crustal structure and can successfully image small blockage features often missed by analysis of raypaths alone. This map can then be used to screen out events with blocked Lg prior to performing Q tomography, and to avoid using Lg-based methods of event identification for the CTBT in regions where they cannot work.¶For this study we applied our technique to one of the most tectonically complex regions on the earth. Nearly 9000 earthquake/station raypaths, traversing the vast region comprised of the Middle East, Mediterranean, Southern Europe and Northern Africa, have been analyzed. We measured the amplitude of Lg relative to the P-coda and mapped the lateral variation of Lg propagation efficiency. With the relatively dense coverage provided by the numerous crossing paths we are able to map out the pattern of crustal heterogeneity that gives rise to the observed character of Lg propagation. We observe that the propagation characteristics of Lg within the region of interest are very complicated but are readily correlated with the different tectonic environments within the region. For example, clear strong Lg arrivals are observed for paths crossing the stable continental interiors of Northern Africa and the Arabian Shield. In contrast, weakened to absent Lg is observed for paths crossing much of the Middle East, and Lg is absent for paths traversing the Mediterranean. Regions that block Lg transmission within the Middle East are very localized and include the Caspian Sea, the Iranian Plateau and the Red Sea. Resolution is variable throughout the region and strongly depends on the distribution of seismicity and recording stations. Lg propagation is best resolved within the Middle East where regions of crustal heterogeneity on the order of 100?km are imaged (e.g., South Caspian Sea and Red Sea). Crustal heterogeneity is resolvable but is poorest in seismically quiescent Northern Africa.