Crustal structure of Northwest Zagros (Kermanshah) and Central Iran (Yazd and Isfahan) using teleseismic Ps converted phases

Receiver functions are widely employed to detect P-to-S converted waves and are especially useful to image seismic discontinuities in the crust. In this study we used the P receiver function technique to investigate the velocity structure of the crust beneath the Northwest Zagros and Central Iran and map out the lateral variation of the Moho boundary within this area. Our dataset includes teleseismic data (M _b ≥ 5.5, epicentral distance from 30° to 95°) recorded at 12 three-component short-period stations of Kermanshah, Isfahan and Yazd telemetry seismic networks. Our results obtained from P receiver functions indicate clear Ps conversions at the Moho boundary. The Moho depths were firstly estimated from the delay time of the Moho converted phase relative to the direct P wave beneath each network. Then, we used the P receiver function inversion to find the properties of the Moho discontinuity such as depth and velocity contrast. Our results obtained from PRF are in good agreement with those obtained from the P receiver function modeling. We found an average Moho depth of about 42 km beneath the Northwest Zagros increasing toward the Sanandaj-Sirjan Metamorphic Zone and reaches 51 km, where two crusts (Zagros and Central Iran) are assumed to be superposed. The Moho depth decreases toward the Urmieh-Dokhtar Cenozoic volcanic belt and reaches 43 km beneath this area. We found a relatively flat Moho beneath the Central Iran where, the average crustal thickness is about 42 km. Our P receiver function modeling revealed a shear wave velocity of 3.6 km/s in the crust of Northwest Zagros and Central Iran increasing to 4.5 km/s beneath the Moho boundary. The average shear wave velocity in the crust of UDMA as SSZ is 3.6 km/s, which reaches to 4.0 km/s while in SSZ increases to 4.3 km/s beneath the Moho.     

Erratum to: Lithospheric structure of NW Iran from P and S receiver functions

We computed P and S receiver functions to investigate the lithospheric structure beneath the northwest Iran and compute the Vp/Vs ratio within the crust of this seismologically active area. Our results enabled us to map the lateral variations of the Moho as well as those of the lithosphere–asthenosphere boundary (LAB) beneath this region. We selected data from teleseismic events (M_b > 5.5, epicentral distance between 30° and 95° for P receiver functions and M_b > 5.7, epicentral distance between 60° and 85° for S receiver functions) recorded from 1995 to 2008 at 8 three-component short-period stations of Tabriz Telemetry Seismic Network. Our results obtained from P receiver functions indicate clear conversions at the Moho boundary. The Moho depth was firstly estimated from the delay time of the Moho converted phase relative to the direct P wave. Then we used the H-Vp/Vs stacking algorithm of Zhu and Kanamori to estimate the crustal thickness and Vp/Vs ratio underneath the stations with clear Moho multiples. We found an average Moho depth of 48 km, which varies between 38.5 and 53 km. The Moho boundary showed a significant deepening towards east and north. This may reveal a crustal thickening towards northeast possibly due to the collision between the Central Iran and South Caspian plates. The obtained average Vp/Vs ratio was estimated to be 1.76, which varies between 1.73 and 1.82. The crustal structure was also determined by modeling of P receiver functions. We obtained a three-layered model for the crust beneath this area. The thickness of the layers is estimated to be 6–11, 18–35, and 38–53 km, respectively. The average of the shear wave velocity was calculated to be 3.4 km/s in the crust and reaches 4.3 km/s below the Moho discontinuity. The crustal thickness values obtained from P receiver functions are in good agreement with those derived by S receiver functions. In addition, clear conversions with negative polarity were observed at ~8.7 s in S receiver functions, which could be related to the conversion at the LAB. This may show a relatively thin continental lithosphere of about 85 km implying that the lithosphere was influenced by various geodynamical reworking processes in the past.     

Topography of the Seismic Velocity Discontinuities beneath the Chugoku and Shikoku Districts,Southwest Japan

The first P-arrival-time data from 513 local earthquakes were analyzed to study lateral variation of the depth to the Conrad and Moho discontinuities beneath the Chugoku and Shikoku districts, southwest Japan, as well as to determine earthquake hypocenters and P-wave station corrections. The depth to the discontinuity was estimated by minimizing the travel-time residuals of more than 8700 first P arrivals observed at 55 seismic stations. The Conrad and Moho discontinuities are located within depth ranges of 15–25 km and 30–40 km, respectively. The Moho is deeper under the mountain area than under the Seto Inland Sea area, and especially deep under the Pacific Coast of the Shikoku district and the mountain area in the Chugoku district. The depth variation of the Moho is quite similar to the Bouguer gravity anomaly distribution and the lateral variations of the P-wave velocity. The deep Moho under the southern Shikoku is located at the portion in which the continental Moho under the island arc meets the oceanic Moho that is the boundary interface between the oceanic crust and the Philippine Sea (PHS) plate dipping toward the back arc. Although there are high mountains in the northern and middle Shikoku, the Moho is not so deep because subduction of the PHS plate prevents the Moho from getting deep, while the Moho is deep due to isostatic balance under the mountain area in the Chugoku district. In addition, we indicated the possibility that the upper boundary of the oceanic crust just above the high-velocity PHS plate is in contact with the deep Moho under the western Chugoku. The contact of the Moho with the oceanic crust can explain the markedly negative gravity anomaly observed in the western Chugoku and the later phase that appears just after the first P arrival from local earthquakes.     

Lower-crustal earthquakes caused by magma movement beneath Askja volcano on the north Iceland rift

The lower crust of magmatically active rifts is usually too hot and ductile to allow seismicity. The Icelandic mid-Atlantic rift is characterized by high heat flow, abundant magmatism generating up to 25–30 km thick crust, and seismicity within the upper 8 km of the crust. In a 20-seismometer survey in July-August 2006 within the northern rift zone around the Askja volcano we recorded ~1700 upper-crustal earthquakes cutting off at 7–8 km depth, marking the brittle-ductile boundary. Unexpectedly, we discovered 100 small-magnitude (M_L <1.5) earthquakes, occurring in swarms mostly at 14–26 km depth within the otherwise aseismic lower crust, and beneath the completely aseismic middle crust. A repeat survey during July-August 2007 yielded more than twice as many lower-crustal events. Geodetic and gravimetric data indicate melt drainage from crustal magma chambers beneath Askja. We interpret the microearthquakes to be caused by melt moving through the crust from the magma source feeding Askja. They represent bursts of magma motion opening dykes over distances of a few meters, facilitated by the extensional setting of the active rift zone.     

New constraints on lithospheric thickness of the Iranian plateau using converted waves

The study of mantle lithosphere plays a key role to reveal predominant tectonic setting process of a region. The current geological and tectonic setting of Iran is due to the ongoing continental–continental collision of the Arabian and Eurasian plates. We applied a combined P and S receiver function analysis to the teleseismic data of nine permanent broadband seismic stations of the International Institute of Earthquake Engineering and Seismology located in different tectonic zones of Iranian plateau. More than 4 years of data were used to estimate the thickness of the crust and mantle lithosphere. According to our results, the crust is 50 km thick beneath the Zagros fold and thrust belt (ZFTB). We found the maximum Moho depth of approximately 70 km under the Sanandaj-Sirjan zone (SSZ) indicating the overthrusting of the crust of Central Iran onto the Zagros crust along the main Zagros thrust (MZT). Below the northeasternmost part of the Urumieh–Dokhtar Magmatic Arc (UDMA) and Central Iran, the Moho becomes shallower and lies at 40 km depth. Towards northeast, beneath the Alborz zone, the crust is 55 km thick. Based on S receiver functions, we provided new insights into the thickness of the Arabian and Eurasian lithospheres. The location of the boundary between these plates was estimated to be beneath the SSZ, which is slightly shifted northeastward relative to the surficial expression of the MZT. Furthermore, the Arabian plate is characterized by the relatively thick lithosphere of about 130 km beneath the ZFTB reaching 150 km beneath the SSZ, where the thickest crust was also observed. This may imply that the shortening across the Zagros is accommodated by lithospheric thickening. In contrast, UDMA and Central Iran are recognized by the thin lithosphere of about 80–85 km. This thin lithosphere may be associated with the asthenospheric upwelling caused by either lithospheric delamination or Neo-Tethys slab detachment beneath the Zagros collision zone.     

A detailed receiver function image of the lithosphere beneath the Kopeh-Dagh (Northeast Iran)

The Kopeh-Dagh fold belts are among the most seismically active areas in Northeastern Iran, which build the northern part of the Alpine–Himalayan orogen in western Asia. They act as the abrupt northeastern limit to active deformation in Iran. We perform a combined P and S receiver function analysis to detect the major discontinuities within the lithosphere beneath Northeast Iran. Our results obtained from 12 short period and broadband seismological stations significantly map the lateral variations of the Moho boundary. Based on P receiver functions, we show that the Moho depth varies from ～43 km beneath the southern Kopeh-Dagh foreland basin to ～49 km below the northern part of the basin. S receiver functions reliably reveal an average Moho depth of ～50–55 km beneath the Kopeh-Dagh mountain range showing the regional shortening in response to the collision of Arabia with Eurasia. Furthermore, we observe clear conversions with negative polarity at ～8.5–9.5 s in S receiver functions, which could be related to the conversion at the lithosphere–asthenosphere boundary. This may show a relatively thin continental lithosphere of about 85–95 km beneath the Kopeh-Dagh implying that the lithosphere was influenced by geodynamical reworking processes in the past.     

用转换函数方法研究喜马拉雅地区速度结构

利用流动数字地震台网提供的三分量地震波形记录,应用转换函数及快速模拟退火算法对喜马拉雅山脉地区46个地震站下的地壳横波速度结构进行了反演,为进一步揭示青藏高原喜马拉雅山脉地区的动力学演化过程提供了可靠的地球物理证据.根据本文结果可清晰看到,喜马拉雅山脉地区作为当今地壳活动最活跃的地区,物质交换非常活跃,地下结构远远未达到平衡,地壳速度有很大差异,在板块边界处莫霍界面速度间断不是非常明显,自喜马拉雅南坡向高原腹地,地壳厚度大致从55 km增加到80 km;沿经度方向,莫霍面也有一定的起伏.通过研究得到另外一个证据是,在喜马拉雅的主中央逆冲断裂,由大陆碰撞产生的主要构造,其深度可能要大于80 km.     

Small-scale heterogeneities below the Lanzhou CTBTO seismic array,from seismic wave field fluctuations

Seismic wave field fluctuations below Lanzhou Comprehensive Test Ban Treaty Organization seismic array are determined and interpreted to describe the scattering characteristics of the crust beneath the northeastern margin of Qinghai–Tibet plateau in China. The frequency-dependent intensities of the mean and fluctuation wave fields are analyzed for 21 deep (206–632 km depth) teleseismic events. The observed wave field fluctuations in the frequency range 0.5–2.5 Hz can be explained by the scattering of the teleseismic P wave front at random media-type structures. The stable and narrow range of the fluctuation parameters from these teleseismic events indicates evidence for scattered seismic phases which are generated inside the crust. The reasonable structural models contain heterogeneities with 1–3% velocity fluctuations and 2.4–8.2 km correlation lengths in the crust with 53 km thickness.     

Upper mantle structure in the ocean-continent transition zone: Kamchatka

We consider results from modeling the crustal and upper mantle velocity structure in Kamchatka by seismic tomography and compare these with gravity data and present-day tectonics. We found a well-pronounced (in the physical fields) vertical and lateral variation for the upper mantle and found that it is controlled by fault tectonics. Not only are individual lithosphere blocks moving along faults, but also parts of the Benioff zone. The East Kamchatka volcanic belt (EKVB) is confined to the asthenospheric layer (the asthenosphere lens) at a depth of 70–80 km; this lens is 10–20 km thick and seismic velocity in it is lower by 2–4%. The top of the asthenosphere lens has the shape of a dome uplift beneath the Klyuchevskoi group of volcanoes and its thickness is appreciably greater; overall, the upper mantle in this region is appreciably stratified. A low-velocity heterogeneity (asthenolith) at least 100 km thick has been identified beneath the Central Kamchatka depression; we have determined its extent in the upper mantle and how it is related to the EKVB heterogeneities. Gravity data suggest the development of a rift structure under the Sredinnyi Range volcanic belt. The Benioff zone was found to exhibit velocity inhomogeneity; the anomalous zones that have been identified within it are related to asthenosphere inhomogeneities in the continental and oceanic blocks of the mantle.     

Structure of the crust and uppermost mantle from broadband seismic array in the western Dabie Mountains,east central China

A broadband seismic array of 7 stations was set up in the western Dabie Mountains (31°20′-31°50′N, 114°30′-115°E). Teleseismic events from May 2001 to November 2001 were collected and analyzed by radial receiver function to determine the S-wave velocity structure of the crust and uppermost mantle. The crustal thickness is 32-38 km beneath the array. The crust-mantle boundary appears as a gently north-dipping velocity discontinuity, but turns to be a velocity gradient beneath a station near the Qiliping shea...     

A 3D model of the central asia mantle from dispersion of Rayleigh-wave phase velocities

Frequency-time analysis was applied to records of Rayleigh surface waves due to teleseismic earthquakes (M ≥ 6.0) to obtain dispersion curves of fundamental mode phase velocities for paths between 51 pairs of digital seismic stations in the Asian IRIS networks; the range of periods was 10–200 s. For each of these pairs we derived 1D-shear-velocity sections that fit the corresponding dispersion curves and image the integral earth structure down to depths of ∼650–700 km beneath each of the interstation seismic paths, which traverse the area of study in different directions. These structures were used by Yanovskaya (2001) to derive a 3D-shear-velocity model by 2D tomography, imaging the distribution of the larger horizontal inhomogeneities in the Central Asia mantle for depths of 50 to 600 km, and to construct 2D-velocity sections for five lines passing through the major tectonic features of Central Asia.     

Stress Rotation in the Kachchh Rift Zone, Gujarat, India

Double difference relocations of the 1402 Kachchh events (2001–2006) clearly delineate two fault zones viz. south-dipping North Wagad fault (NWF) and almost vertical Gedi fault (GF). The relocated focal depths delineate a marked variation of 4 and 7 km in the brittle-ductile transition depths beneath GF and NWF, respectively. The focal mechanism solutions of 464 aftershocks (using 8–12 first motions) show that the focal mechanisms ranged between pure reverse and pure strike-slip except for a few pure dip-slip solutions. The stress inversions performed for five rectangular zones across the Kachchh rift reveal both clockwise and anticlockwise rotation (7–32°) in the σ₁ orientation within the rupture zone, favoring a heterogeneous stress regime with an average N-S fault normal compression. This rotation may be attributed to the presence of crustal mafic intrusives (5–35 km depth) in the rupture zone of the 2001 Bhuj main shock. Results suggest a relatively homogeneous stress regime in the GF zone favoring strike-slip motion, with a fault normal N-S compression.     

Marine seismological observations in the Central Kurile segment before catastrophic earthquakes on November, 2006 (<Emphasis Type="Italic">M</Emphasis> = 8.3) and January, 2007 (<Emphasis Type="Italic">M</Emphasis> = 8.1)

The results of detailed seismological observations with bottom seismographs in the Central Kurile segment in August-September, 2006 are discussed. The system of six bottom seismographs was placed on the island slope of the Kurile deep-sea trench southeast of Urup Island and southwest of the Bussol Strait. Over 230 earthquakes with M _LH = 0.5–5.5 were registered in the area with a radius of 150 km around the center of the observation system at depths up to 300 km during 16 days. Records of 80 earthquakes with hypocenters in the earth crust (h = 0–30 km) beneath the island slope of the Kurile deep-sea trench were first obtained by bottom seismographs. These data are inconsistent with previous concepts of aseismicity of this zone. The discovery of the unique morphological structure of the Benioff zone beneath the central Kurile Arc represents the most important result of detailed seismological observations. The zone consists of an inner seismoactive subzone, which is located beneath the island slope of the arc at depths of 15–210 km, being characterized by an angle of incline of 50° under the latter and crosses the ocean bottom approximately 80 km away from the trench axis, and outer low-activity subzone. The latter is traceable beyond the trench almost parallel to the inner zone beginning from a depth of 50 km below the sea bottom up to a depth of approximately 300 km. Due to the slightly lower incline (∼45°) of the outer subzone, both subzones gradually converge downward. The integral thickness of the Benioff zone varies from 150 km in its upper part to 125 km at depths of 210–260 km. The medium sandwiched between these subzones is practically aseismic. The reality of this defined structure is confirmed by the distribution of aftershocks of the earthquake that occurred on November 15, 2006 (M = 8.3). These seismic events served as foreshocks for the subsequent strong earthquake of January 13, 2007 (M = 8.1) with the hypocenter located beyond the trench under the ocean bottom. Such a structure of this zone within the central Kurile Arc segment is unique, having no analogues either in the flanks of the Kurile-Kamchatka Arc or other arcs. The results of detailed seismological observations obtained two months before the first of the catastrophic Central Kurile earthquakes appeared to be typical for the period of foreshocks (the lower seismic activity of the Simushir block, which hosted the hypocenter of the earthquake that occurred on November 15, 2006, particularly at depths of 0–50 km, the gentler incline of the recurrence plot, and other features).     

Path-specific Velocity Structure of Western China from Surface-wave Dispersion

We develop one-dimensional (1-D) path-specific velocity models in western China using new Rayleigh and Love wave group and phase velocity dispersion measurements for 20 events in the region. The earthquakes were grouped into three geographic clusters from which we compute the average phase and group velocity dispersion. We invert the average dispersion curves simultaneously for 1-D shear-velocity models appropriate for the three central Asian paths, using three previous shear-velocity models as initial models. The models are validated by forward modeling waveforms of recent events. The crustal thickness beneath western China in the vicinity of the Lop Nor test site is 50–60 km and our velocity models are consistent with major geologic features (e.g., basins and mountain ranges) and previous structural models for this region.     

Magma storage and ascent at Lipari Island (Aeolian archipelago, Southern Italy) at 223–81 ka: the role of crustal processes and tectonic influence

Fluid inclusion studies together with volcanological and petrochemical data allow reconstruction of the magma feeding system of basaltic-andesitic to andesitic activity during the oldest and intermediate stages of development of Lipari Island (223–81 ka). A major magma storage zone is active during the overall investigated time span at depths of 22 km, close to the crust-mantle Moho transition, at which mantle-derived mafic magmas tend to accumulate due to neutral buoyancy conditions. Beneath central-type volcanoes (M. Mazzacaruso, M. S.Angelo, M. Chirica-Costa d’Agosto), a shallower magma reservoir is located within the upper crust at 5.5–3.5 km, associated with a major lithological discontinuity. For fissural-type volcanoes (Timpone Ospedale, Monterosa, M. Chirica), tectonic structures are suggested to influence further magma ascent and storage at mid-crustal depths (∼14 km), with no ponding at shallower levels. Partial crustal melting processes at the roofs of the deep magma reservoirs (∼17 km) are invoked to explain the origin of cordierite-bearing lavas beneath M. S.Angelo and M. Chirica-Costa d’Agosto volcanoes, which were active during the intermediate stages of development of Lipari (105–81 ka). The generation of felsic anatectic melts in the lower crust could have created density and rheologic barriers to impede the passage of mafic melts and promote their ponding, with influence on the subsequent evolution of Lipari volcano.     

New data on the Indus Kohistan seismic zone and its extension into the Hazara–Kashmir Syntaxis, NW Himalayas of Pakistan

This paper deals with the data obtained from local networks in northern Pakistan for 251 earthquakes of magnitude ≥4.0 for October 8, 2005 to December 31, 2006 period. The study presents focal mechanism solutions (FMS) of 12 pre- (1904–2005) and 17 post- (October 8, 2005–December, 2005) Muzaffarabad Earthquake, their detailed tectonic interpretation, and correlation with surface evidence of co-seismic rupture with published synthetic aperture radar data. Distribution of landslides obtained from National Engineering Services of Pakistan and the earthquake damages are also discussed. Aftershock distribution, which is more prominent in the crystalline zone (northwest of Muzaffarabad), defines a 50-km-wide NW–SE trending zone that extends for 200 km from the main mantle thrust to the center of the Hazara–Kashmir Syntaxis. The FMS of the main shock and 16 aftershocks having magnitude ≥4.0 indicate thrusting to be the dominant mechanism with rupture planes having NW–SE trend and NE dip. In addition, 12 FMS of pre-Muzaffarabad Earthquake (1904–2004) from the same area have been determined and results are compared. This leads to the conclusion that the wedge-shaped NW–SE trending blind zone, referred to by earlier workers as the Indus Kohistan Seismic Zone (IKSZ), has been activated during the Muzaffarabad earthquake. The right-lateral component in all FMS, supported by the surface evidences, suggests the involvement of Balakot–Bagh Fault (BBF). We propose that the IKSZ is the source of the October 8, 2005 Muzaffarabad earthquake that reactivated the BBF. Furthermore, the IKSZ does not end at the nose of the syntaxis but extends further southeast of it. Tectonic complexity seems to be due to a variety of factors. Also, thrust and reverse solutions near the northern collisional boundary (main mantle thrust) have mostly NE/SW-directed P-axis orientations. From the detailed FMS analysis, three conclusions have been drawn: (1) Shallow events (depth ≤10 km) with prominent strike slip solutions (earlier earthquakes) are associated with the surface strike slip faults (e.g., Muzaffarabad Fault) and/or the Besham domal structure; (2) moderate depth events (depth 10–25 km) with thrust/reverse solutions but having minor right-lateral strike slip component (all Muzaffarabad earthquakes and two earlier) are associated with the IKSZ; (3) deeper earthquakes (depth below IKSZ) with pure thrust/reverse solutions may be related to the under-thrusting of the Indian plate beneath the IKSZ, which represents a major thrust zone. Imbricate thrusting and breaking and thickening of the crust are considered to be caused by steep bending of the under-thrusting plate at the collisional boundary. An erratum to this article can be found at     

Crustal thickness variations in the Zagros continental collision zone (Iran) from joint inversion of receiver functions and surface wave dispersion

Variations in crustal thickness in the Zagros determined by joint inversion of P wave receiver functions (RFs) and Rayleigh wave group and phase velocity dispersion. The time domain iterative deconvolution procedure was employed to compute RFs from teleseismic recordings at seven broadband stations of INSN network. Rayleigh wave phase velocity dispersion curves were estimated employing two-station method. Fundamental mode Rayleigh wave group velocities for each station is taken from a regional scale surface wave tomographic imaging. The main variations in crustal thickness that we observe are between stations located in the Zagros fold and thrust belt with those located in the Sanandaj–Sirjan zone (SSZ) and Urumieh–Dokhtar magmatic assemblage (UDMA). Our results indicate that the average crustal thickness beneath the Zagros Mountain Range varies from ～46 km in Western and Central Zagros beneath SHGR and GHIR up to ～50 km beneath BNDS located in easternmost of the Zagros. Toward NE, we observe an increase in Moho depth where it reaches ～58 km beneath SNGE located in the SSZ. Average crustal thickness also varies beneath the UDMA from ～50 km in western parts below ASAO to ～58 in central parts below NASN. The observed variation along the SSZ and UDMA may be associated to ongoing slab steepening or break off in the NW Zagros, comparing under thrusting of the Arabian plate beneath Central Zagros. The results show that in Central Iran, the crustal thickness decrease again to ～47 km below KRBR. There is not a significant crustal thickness difference along the Zagros fold and thrust belt. We found the same crystalline crust of ～34 km thick beneath the different parts of the Zagros fold and thrust belt. The similarity of crustal structure suggests that the crust of the Zagros fold and thrust belt was uniform before subsidence and deposition of the sediments. Our results confirm that the shortening of the western and eastern parts of the Zagros basement is small and has only started recently.     

Lithospheric structure of NW Iran from P and S receiver functions

We computed P and S receiver functions to investigate the lithospheric structure beneath the northwest Iran and compute the Vp/Vs ratio within the crust of this seismologically active area. Our results enabled us to map the lateral variations of the Moho as well as those of the lithosphere–asthenosphere boundary (LAB) beneath this region. We selected data from teleseismic events (M_b?>?5.5, epicentral distance between 30° and 95° for P receiver functions and M_b?>?5.7, epicentral distance between 60° and 85° for S receiver functions) recorded from 1995 to 2008 at 8 three-component short-period stations of Tabriz Telemetry Seismic Network. Our results obtained from P receiver functions indicate clear conversions at the Moho boundary. The Moho depth was firstly estimated from the delay time of the Moho converted phase relative to the direct P wave. Then we used the H-Vp/Vs stacking algorithm of Zhu and Kanamori to estimate the crustal thickness and Vp/Vs ratio underneath the stations with clear Moho multiples. We found an average Moho depth of 48 km, which varies between 38.5 and 53 km. The Moho boundary showed a significant deepening towards east and north. This may reveal a crustal thickening towards northeast possibly due to the collision between the Central Iran and South Caspian plates. The obtained average Vp/Vs ratio was estimated to be 1.76, which varies between 1.73 and 1.82. The crustal structure was also determined by modeling of P receiver functions. We obtained a three-layered model for the crust beneath this area. The thickness of the layers is estimated to be 6–11, 18–35, and 38–53 km, respectively. The average of the shear wave velocity was calculated to be 3.4 km/s in the crust and reaches 4.3 km/s below the Moho discontinuity. The crustal thickness values obtained from P receiver functions are in good agreement with those derived by S receiver functions. In addition, clear conversions with negative polarity were observed at ~8.7 s in S receiver functions, which could be related to the conversion at the LAB. This may show a relatively thin continental lithosphere of about 85 km implying that the lithosphere was influenced by various geodynamical reworking processes in the past.     

Magnetotelluric Studies of Active Continent–Continent Collisions

Continent–continent collisions are an important tectonic process and have played a fundamental role in the evolution of the modern continents. A combination of geological and geophysical data has provided new constraints on the structure and temporal evolution of these orogens. Magnetotelluric (MT) studies have been an important part of these studies since they can constrain the fluid content and thermal structure which are key parameters for defining the rheology of the crust and upper mantle. MT studies of the Himalaya have defined the geometry of active faults associated with continued plateau growth. Orogen scale MT studies have shown that both the India–Asia collision (Tibetan Plateau and Himalaya) and the Arabia–Eurasia collision (Eastern Anatolia) have developed a low resistivity mid-crustal layer with upper surface at 10–20 km that is likely due to a combination of partial melt and associated aqueous fluids. The properties of this layer are consistent with a strength contrast that permits crustal flow over geological timescales. The upper mantle from the Moho to at least 100 km beneath both Northern Tibet and the Anatolian Plateau is characterized by low resistivity values (10–30 Ωm) indicating the presence of shallow asthenosphere. Future integrated seismic and MT studies of collision zones are needed fully to explore the 3D structures associated with deformation and further constrain geodynamic models.     

Topography of the Moho and Conrad discontinuities in the Kyushu district,Southwest Japan

The first P-arrival time data from local earthquakes are inverted for two-dimensional variation of the depths to the Conrad and Moho discontinuities in the Kyushu district, southwest Japan. At the same time, earthquake hypocenters and station corrections are determined from the data. The depths to the discontinuities are estimated by minimizing the travel time residuals of first P-arrival phases for 608 earthquakes observed at 57 seismic stations. In the land area of Kyushu, the Conrad and Moho discontinuities are located within the depth ranges of 16–18 and 34–40 km, respectively. The Conrad discontinuity is not as largely undulated as the Moho discontinuity. The depth to the Moho is deep along the east coast of Kyushu, and the deepest Moho is closely related to markedly low velocity of P wave. We regard the deepest Moho as reflecting the Kyushu–Palau ridge subducting beneath the Kyushu district, together with the Philippine Sea slab. In western Kyushu, the shallow Moho is spreading in the north–northeast–south–southwest direction in the Okinawa trough region. Based on the presence of low-velocity anomaly in three-dimensional velocity structure and seismogenic stress field of shallow crustal earthquakes, the shallow Moho is interpreted as being due to lower crustal erosion associated with a small-scale mantle upwelling in the Okinawa trough region. The velocity discontinuity undulation basically has insignificant effect on hypocenter determination of the local earthquakes, but the Moho topography makes changes in focal depths of some upper mantle earthquakes. The depth variation of the Moho discontinuity has a good correlation with the Bouguer gravity anomaly map; i.e., the shallow Moho of western Kyushu and the deep Moho of eastern Kyushu closely correlate with the positive and negative Bouguer gravity anomalies, respectively.