Regional Variations of the <Emphasis Type="Italic">ω</Emphasis>-upper Bound Magnitude of GIII Distribution in the Iranian Plateau

The Iranian Plateau does not appear to be a single crustal block, but an assemblage of zones comprising the Alborz—Azerbaijan, Zagros, Kopeh—Dagh, Makran, and Central and East Iran. The Gumbel’s III asymptotic distribution method (GIII) and maximum magnitude expected by Kijko—Sellevoll method is applied in order to check the potentiality of the each seismogenic zone in the Iranian Plateau for the future occurrence of maximum magnitude (M_max). For this purpose, a homogeneous and complete seismicity database of the instrumental period during 1900–2012 is used in 29 seismogenic zones of the examined region. The spatial mapping of hazard parameters (upper bound magnitude (ω), most probable earthquake magnitude in next 100 years (M₁₀₀) and maximum magnitude expected by maximum magnitude estimated by Kijko—Sellevoll method (max M_{K ? S}^max) reveals that Central and East Iran, Alborz and Azerbaijan, Kopeh—Dagh and SE Zagros are a dangerous place for the next occurrence of a large earthquake.     

An evaluation of earthquake hazard parameters in the Iranian Plateau based on the Gumbel III distribution

The Gumbel’s third asymptotic distribution (GIII) of the extreme value method is employed to evaluate the earthquake hazard parameters in the Iranian Plateau. This research quantifies spatial mapping of earthquake hazard parameters like annual and 100-year mode beside their 90 % probability of not being exceeded (NBE) in the Iranian Plateau. Therefore, we used a homogeneous and complete earthquake catalogue during the period 1900–2013 with magnitude M _w ? ≥?4.0, and the Iranian Plateau is separated into equal area mesh of 1° late?×?1° long. The estimated result of annual mode with 90 % probability of NBE is expected to exceed the values of M _w 6.0 in the Eastern part of Makran, most parts of Central and East Iran, Kopeh Dagh, Alborz, Azerbaijan, and SE Zagros. The 100-year mode with 90 % probability of NBE is expected to overpass the value of M _w 7.0 in the Eastern part of Makran, Central and East Iran, Alborz, Kopeh Dagh, and Azerbaijan. The spatial distribution of 100-year mode with 90 % probability of NBE uncovers the high values of earthquake hazard parameters which are frequently connected with the main tectonic regimes of the studied area. It appears that there is a close communication among the seismicity and the tectonics of the region.     

GPS-based crustal deformations in Azerbaijan and their influence on seismicity and mud volcanism

Using Shen’s method (Shen et al., 1996), deformations of the Earth’s crust in Azerbaijan were studied based on GPS measurements. For estimating the rate of deformation, we used the field of velocity vectors for Azerbaijan, Iran, Georgia, and Armenia that were derived from GPS measurements during 1998–2012. It is established that compression is observable along the Greater Caucasus, in Gobustan, the Kura depression, Nakhchyvan Autonomous Republic, and adjacent areas of Iran. The axes of compression/contraction of the crust in the Greater Caucasus region are oriented in the S-NE direction. The maximum strain rate (approximately 200 × 10^?9 per annum) is documented in the zone of mud volcanism at the SHIK site (Shykhlar), which is marked by a sharp change in the direction of the compression axes (SW-NE). It is revealed that the deformation field also includes the zones where strain rates are very low approximating 5 × 10^?9 per annum. These zones include the Caspian-Guba and northern Gobustan areas, characterized by extensive development of mud volcanism. The extension zones are confined to the Lesser Caucasus and are revealed in the Gedabek (GEDA) and Shusha (SHOU) areas, as well as in the zone located between the DAMO and PIRM sites (Iran), where the deformation rate amounts to 100 × 10^?9 per annum. It is concluded that the predominant factor responsible for the eruption of mud volcanoes is the intensity of gas-generation processes in the earth’s interior, while deformation processes play the role of a trigger. The zone of the epicenters of strong earthquakes is correlated to the gradient zone in the crustal strain rates.     

Iranian earthquakes, a uniform catalog with moment magnitudes

A uniform earthquake catalog is an essential tool in any seismic hazard analysis. In this study, an earthquake catalog of Iran and adjacent areas was compiled, using international and national databanks. The following priorities were applied in selecting magnitude and earthquake location: (a) local catalogs were given higher priority for establishing the location of an earthquake and (b) global catalogs were preferred for determining earthquake magnitudes. Earthquakes that have occurred within the bounds between 23–42° N and 42–65° E, with a magnitude range of M _W 3.5–7.9, from the third millennium BC until April 2010 were included. In an effort to avoid the “boundary effect,” since the newly compiled catalog will be mainly used for seismic hazard assessment, the study area includes the areas adjacent to Iran. The standardization of the catalog in terms of magnitude was achieved by the conversion of all types of magnitude into moment magnitude, M _W, by using the orthogonal regression technique. In the newly compiled catalog, all aftershocks were detected, based on the procedure described by Gardner and Knopoff (Bull Seismol Soc Am 64:1363–1367, 1974). The seismicity parameters were calculated for the six main tectonic seismic zones of Iran, i.e., the Zagros Mountain Range, the Alborz Mountain Range, Central Iran, Kope Dagh, Azerbaijan, and Makran.     

Analysis of the stress field and strain rate in Zagros-Makran transition zone

Transition boundary between Zagros continental collision and Makran oceanic-continental subduction can be specified by two wide limits: (a) Oman Line is the seismicity boundary with a sizeable reduction in seismicity rate from Zagros in the west to Makran in the east; and (b) the Zendan-Minab-Palami (ZMP) fault system is believed to be a prominent tectonic boundary. The purpose of this paper is to analyze the stress field in the Zagros-Makran transition zone by the iterative joint inversion method developed by Vavrycuk (Geophysical Journal International 199:69-77, 2014). The results suggest a rather uniform pattern of the stress field around these two boundaries. We compare the results with the strain rates obtained from the Global Positioning System (GPS) network stations. In most cases, the velocity vectors show a relatively good agreement with the stress field except for the Bandar Abbas (BABS) station which displays a relatively large deviation between the stress field and the strain vector. This deviation probably reflects a specific location of the BABS station being in the transition zone between Zagros continental collision and Makran subduction zones.     

The maximum expected earthquake magnitudes in different future time intervals of six seismotectonic zones of Iran and its surroundings

One of the crucial components in seismic hazard analysis is the estimation of the maximum earthquake magnitude and associated uncertainty. In the present study, the uncertainty related to the maximum expected magnitude μ is determined in terms of confidence intervals for an imposed level of confidence. Previous work by Salamat et al. (Pure Appl Geophys 174:763-777, 2017) shows the divergence of the confidence interval of the maximum possible magnitude m_max for high levels of confidence in six seismotectonic zones of Iran. In this work, the maximum expected earthquake magnitude μ is calculated in a predefined finite time interval and imposed level of confidence. For this, we use a conceptual model based on a doubly truncated Gutenberg-Richter law for magnitudes with constant b-value and calculate the posterior distribution of μ for the time interval T_f in future. We assume a stationary Poisson process in time and a Gutenberg-Richter relation for magnitudes. The upper bound of the magnitude confidence interval is calculated for different time intervals of 30, 50, and 100 years and imposed levels of confidence α?=?0.5, 0.1, 0.05, and 0.01. The posterior distribution of waiting times T_f to the next earthquake with a given magnitude equal to 6.5, 7.0, and 7.5 are calculated in each zone. In order to find the influence of declustering, we use the original and declustered version of the catalog. The earthquake catalog of the territory of Iran and surroundings are subdivided into six seismotectonic zones Alborz, Azerbaijan, Central Iran, Zagros, Kopet Dagh, and Makran. We assume the maximum possible magnitude m_max?=?8.5 and calculate the upper bound of the confidence interval of μ in each zone. The results indicate that for short time intervals equal to 30 and 50 years and imposed levels of confidence 1???α?=?0.95 and 0.90, the probability distribution of μ is around μ?=?7.16???8.23 in all seismic zones.     

Activation of seismicity in Central and South Asia after the Makran earthquakes: Possible acceleration of preparation of large seismic events in the Tien Shan region

Two zones of seismicity (ten events with M _w = 7.0–7.7) stretching from Makran and the Eastern Himalaya to the Central and EasternTien Shan, respectively, formed over 11 years after the great Makran earthquake of 1945 (M _w = 8.1). Two large earthquakes (M _w = 7.7) hit theMakran area in 2013. In addition, two zones of seismicity (M ≥ 5.0) occurred 1–2 years after theMakran earthquake in September 24, 2013, stretching in the north-northeastern and north-northwestern directions. Two large Nepal earthquakes struck the southern extremity of the “eastern” zone (April 25, 2015, M _w = 7.8 and May 12, 2015, M _w = 7.3), and the Pamir earthquake (December 7, 2015, M _w = 7.2) occurred near Sarez Lake eastw of the “western” zone. The available data indicate an increase in subhorizontal stresses in the region under study, which should accelerate the possible preparation of a series of large earthquakes, primarily in the area of the Central Tien Shan, between 70° and 79° E, where no large earthquakes (M _w ≥ 7.0) have occurred since 1992.     

Active faulting of the southeastern-most Zagros (Iran): Microearthquake seismicity and crustal structure

The microseismicity of the southeastern-most Zagros is examined by high-resolution data recorded by a temporary dense local seismic network. The seismicity defines a diffuse pattern, mostly located beneath folds in the southern part of the High Zagros Fault (HZF). Seismicity dips gently northward in the depth range 6–25 km, implying slip on a major intracrustal thrust fault extending to the north of the Main Zagros Reverse Fault (MZRF) which seems to connect to the Mountain Frontal Fault (MFF). Furthermore, observed focal mechanisms suggest transpressive motion on the HZF located west of the Zendan-Minab-Palami (ZMP) fault system and striking obliquely to the convergent motion. These observations suggest that the transition zone between the Zagros continental collision zone and the Makran oceanic subduction zone is not confined to the east of the ZMP and some part of the this diffuse transition is accommodated north of the Hormuz Strait in the west by partitioning between strike-slip and shortening components. The Zagros reverse domain is terminated by a transpressive tectonic regime. Moho depth beneath the MZRF, deduced from receiver functions, is almost 45 km thinner than is observed in the central and northern parts of the Zagros. These observations support a model of active underthrusting of the Arabian plate beneath central Iran in the southeastern-most Zagros.     

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.     

Stromboli and its 1975 eruption

On November 4, 1975 in the evening, an eruption took place at Mt. Stromboli. On the following day lava flowed on the Sciara del Fuoco downward to the sea, accompanied by an intense explosive activity at the crater plane. Direct observations on the volcanic activity were carried out since November 6 while a seismic survey was made from Nov. 7 to 12. The total volume of the lava outpoured during this period of activity that lasted 21 days, was estimated to be about 10⁴ m³. This paper reports the results of direct observations, and of the petrological, radioactive disequilibria and seismic activity studies performed for this eruption. The eruption was preceded by an insignificant change of seismic activity, which was monitored by a seismic station located about 2 km East of the crater. A shallow seismicity was strietly related to crater explosions accompanying the eruptive phenomenon. Radioactive disequilibria showed a lack of disequilibrium between²²⁸Ra and²³²Th explainable in terms of a fast rising of magma in the conduit. Chemical analyses of lava samples and deep seismic sounding data indicate a correspondence between the depth (10–15 km) at which crystallization pressure of phenocrysts occurs and a low velocity laver.     

Assessment of seismic hazard in the North Caucasus

The seismicity of the North Caucasus is the highest in the European part of Russia. The detection of potential seismic sources here and long-term prediction of earthquakes are extremely important for the assessment of seismic hazard and seismic risk in this densely populated and industrially developed region of the country. The seismogenic structures of the Iran-Caucasus-Anatolia and Central Asia regions, adjacent to European Russia, are the subjects of this study. These structures are responsible for the specific features of regional seismicity and for the geodynamic interaction with adjacent areas of the Scythian and Turan platforms. The most probable potential sources of earthquakes with magnitudes M = 7.0 ± 0.2 and 7.5 ± 0.2 in the North Caucasus are located. The possible macroseismic effect of one of them is assessed.     

Tectonic implications of the geochemical data from the Makran igneous rocks in Iran

Makran is one of the largest accretionary prisms on Earth, formed by the closure of the Neotethys ocean which is now represented by its remnant, the Gulf of Oman. Tectonic evolution of the Makran island‐arc system is explored within the context of a north dipping subduction zone, with temporal variations in slab dip arrangement. In a Middle Jurassic–Early Paleocene steep slab dip arrangement, the Mesozoic magmatic arc and the Proto‐Jaz Murian depression, which was an intra‐arc extensional basin, were developed. This was associated with development of outer‐arc ophiolitic mélange and oceanward migration of the Bajgan–Durkan continental sliver, which is the continuation of the Sanandaj–Sirjan zone of the Zagros orogenic belt into the Makran region. In a Late Paleocene to Late Pliocene moderate to shallow slab dip arrangement, compression and tectonic inversion of the Proto‐Jaz Murian extensional basin into the Jaz Murian compressive basin was associated with the uplift of the southern part of the Jaz Murian Depression along the South Jaz Murian Fault, and emplacement of the Paleogene–Neogene magmatic arc, behind the Jaz Murian compressive basin. A shallow slab dip arrangement in the Quaternary led to the emplacement of a third magmatic arc inland, over the southern part of the Yazd–Tabas–Lut micro‐continental block. It is envisioned that the Makran island‐arc system will pass through similar tectonic events in the future, as the Zagros island‐arc system did in the past. However, a future remnant and/or residual basin similar to the present Gulf of Oman will continue to survive to the east.     

b-value and fractal dimension variations in Iran

Iran sits on a region with a high intrinsic level of seismic activity due to its tectonic setting. Through statistical examination of the earthquakes listed in the catalogue from International Institute of Earthquake Engineering and Seismology (IIEES), this research attempted to calculate some seismicity factors and find correlation between them. A preliminary analysis indicated changes in the b-value of the Gutenberg-Richter relationship over the study region. Thus, the study area was divided into five zones (Alborz, Zagros, Azerbaijan, Central and East) and b-value was computed for each zone. Considering faulting mechanism styles and the b-values in the region, it was found that the lowest b-values belong to the thrust events and strike-slip faulting earthquakes have intermediate values. These findings support previous studies. Furthermore, results of b-value calculation were used for the estimation of accumulated differential stresses (σ₁–σ₂) over each zone. Overall, the b-value for Iran is averagely low which signifies the high stress tectonic regime in this region. Also, by having calculated fractal dimension (D) in each zone, a correlation obtained showing that in Iran region, the b-value correlates to fractal dimension by D = 4.2b–2 relation which does not support Aki's (1981) speculation of D = 3b/c.     

Investigation of mining-induced earthquakes in Iran within a time window of 2006–2013

Iran is one of the most seismically active areas of the world and frequently suffers from destructive earthquakes. Rare studies on anthropogenic-induced seismicity in Iran may be related to less attention to triggered events and more concern to natural origin of earthquakes. Hence, the present study as a frontier research aims to investigate the mining-induced earthquakes in Iran. For this purpose, distribution of ~?76,000 seismic events was investigated between the years 2006 and 2013. This study considered a correlation test to investigate the possible mining triggering of the seismic events based on a network of 194 geographical pixels (1°?×?1°) in ArcGIS. Results conveniently confirmed a positive meaningful relation between all earthquake events with magnitudes M?>?0.5 and mining activities in Iran (R?=?0.42). Detailed results confirmed that the most of earthquake swarms (at least ~?60%) had mining-induced origin, which were spatially located in same pixels of metallic mineral mining sites. The correlation test between earthquake swarms and mining activities indicated positive and meaningful relationships in four regions of Alborz, Kopet Dag, Kerman, and Zagros, respectively (R?=?0.61, 0.54, 0.51, and 0.50). Hence, aforementioned seismic regions exposed sensitive seismic responses toward mining triggering effects in Iran.     

Earthquake relocation in the Central Alborz region of Iran using a non-linear probabilistic method

In this study, we calculate accurate absolute locations for nearly 3,000 shallow earthquakes (≤20 km depth) that occurred from 1996 to 2010 in the Central Alborz region of northern Iran using a non-linear probabilistic relocation algorithm on a local scale. We aim to produce a consistent dataset with a realistic assessment of location errors using probabilistic hypocenter probability density functions. Our results indicate significant improvement in hypocenter locations and far less scattering than in the routine earthquake catalog. According to our results, 816 earthquakes have horizontal uncertainties in the 0.5–3.0 km range, and 981 earthquakes are relocated with focal-depth errors less than 3.0 km, even with a suboptimal network geometry. Earthquake relocated are tightly clustered in the eastern Tehran region and are mainly associated with active faults in the study area (the Mosha and Garmsar faults). Strong historical earthquakes have occurred along the Mosha and Garmsar faults, and the relocated earthquakes along these faults show clear north-dipping structures and align along east–west lineations, consistent with the predominant trend of faults within the study region. After event relocation, all seismicity lies in the upper 20 km of the crust, and no deep seismicity (>20 km depth) has been observed. In many circumstances, the seismicity at depth does not correlate with surface faulting, suggesting that the faulting at depth does not directly offset overlying sediments.     

Seismicity study in Pournari reservoir area (W. Greece) 1981–2010

The spatial-temporal evolution of seismicity is examined, during the initial impoundment of Pournari reservoir located on Arachthos River (Western Greece), as well as for the next 30 years. The results show that, despite the relatively moderate-to-high seismicity from west to east, there is no remarkable earthquake in the vicinity before the first reservoir impoundment. Immediately after the impoundment (January 1981), and during the first 4 months, a considerable number of low-magnitude seismic events were recorded in the broader area of the dam. Moreover, two independent major events occurred on March 10, 1981 (M _L ?=?5.6) and April 10, 1981 (M _L ?=?4.7) with focal depths 13 and 10 km, respectively. The detailed analysis of the two corresponding aftershock sequences shows that they present different behaviors (e.g., larger b-value and lower magnitude of the main aftershock) than that of other aftershock sequences in Greece. This seismicity is probably due to triggering, via the water loading mechanism and the undrained response due to a flysch appearance on the reservoir basement. The activation of the thrust fault may be attributed to the bulging of evaporites that characterize the disordered structure of W. Greece, via possible water intake. The detailed processing of the recorded seismicity during the period 1982–2010, in comparison with the variations of Pournari Dam water level, shows an increase of shallow seismicity (h?≤?5 km) in the vicinity of the reservoir up to a 10-km distance—in contrast to the initial period, characterized by a number of deeper events due to the background response change from undrained to drained status.     

Liparitic volcanism of some regional volcanic structures in Kazakhstan,Middle Asia and Caucasus

Liparitic volcanism is a typical feature of the orogenic phase giving rise to the Kazakhstan, Middle Asia and Caucasus folded systems. The main characteristics of the liparitic volcanism common to these three regions are the following:

1. Geo-structural zonation of the volcanic structures.
2. Dismembered Moho surface within the volcanic structures.
3. Synchronous, yet independent evolution of liparitic and andesitic volcanisms.
4. Ignimbritic character of the liparitic volcanism.
5. Lateral petrochemical zonation with some features common to the liparitic and andesitic rock series.

Geo-structural and petrochemical zonations are likely governed in the regions studied by a deep-seated plutonic body.     

The spatiotemporal analysis of the minimum magnitude of completeness <Emphasis Type="Italic">Mc</Emphasis> and the Gutenberg–Richter law <Emphasis Type="Italic">b</Emphasis>-value parameter using the earthquake catalog of Greece

Spatiotemporal mapping the minimum magnitude of completeness Mc and b-value of the Gutenberg–Richter law is conducted for the earthquake catalog data of Greece. The data were recorded by the seismic network of the Institute of Geodynamics of the National Observatory of Athens (GINOA) in 1970–2010 and by the Hellenic Unified Seismic Network (HUSN) in 2011–2014. It is shown that with the beginning of the measurements at HUSN, the number of the recorded events more than quintupled. The magnitude of completeness Mc of the earthquake catalog for 1970–2010 varies within 2.7 to 3.5, whereas starting from April 2011 it decreases to 1.5–1.8 in the central part of the region and fluctuates around the average of 2.0 in the study region overall. The magnitude of completeness Mc and b-value for the catalogs of the earthquakes recorded by the old (GINOA) and new (HUSN) seismic networks are compared. It is hypothesized that the magnitude of completeness Mc may affect the b-value estimates. The spatial distribution of the b-value determined from the HUSN catalog data generally agrees with the main geotectonic features of the studied territory. It is shown that the b-value is below 1 in the zones of compression and is larger than or equal to 1 in the zones dominated by extension. The established depth dependence of the b-value is pretty much consistent with the hypothesis of a brittle–ductile transition zone existing in the Earth’s crust. It is assumed that the source depth of a strong earthquake can probably be estimated from the depth distribution of the b-value, which can be used for seismic hazard assessment.     

Seismotectonics of the Bitlis–Zagros Fold and Thrust Belt in Northern Iraq and Surrounding Regions from Moment Tensor Analysis

Northern Iraq represents part of the convergent plate boundary between the Arabian and Eurasian plates. The collision zone between these two plates is manifested by the Bitlis–Zagros Fold and Thrust Belt. This belt is one of the most seismically active regions among the present active belts. This study intends to improve our knowledge on the seismotectonic activities in northern Iraq and the surrounding areas. To reach this goal, we used the waveform moment tensor inversion method to determine the focal depths, moment magnitudes, fault plane solutions, and directions of the principal stress axes of 25 events with magnitudes ≥3.5. The seismic data of these events were collected from 54 broadband stations which belong to the Kandilli Observatory and Earthquake Research Institute, the Incorporated Research Institutions for Seismology, the Observatories and Research Facilities for European Seismology, and the Iraqi Seismological Network. Computer Programs in Seismology, version 3.30 (Herrmann and Ammon 2004), was used for analysis. The results show that the focal depth of these events ranged from 15 to 25 km in general. The fault plane solutions show that the strike-slip mechanism is the most dominant mechanism in the study area, usually with a reverse component. The stress regime shows three major directions; north–south, northeast-southwest, and east–west. These directions are comparable with the tectonic regime in the region.     

A thinned lithosphere beneath coastal area of southeastern China as evidenced by seismic receiver functions

During Mesozoic to Cenozoic time, the large-scale tectono-magmatism had strongly modified the lithosphere beneath the southeastern continent of China, leaving the present-day lithosphere as a new one evolving from the ancient lithosphere that was largely removed and replaced. But this model proposed from geochemical and petrological research is urgently in need of support from seismic observational evidence. In this paper, based on the dataset recorded by the dense stations of two NE oriented broadband seismic profiles deployed in the coastal area of southeastern China(SE China), both P-wave(P-RF) and S-wave(S-RF) receiver functions were isolated. We identified Pls phase converted from the Lithosphere-Asthenosphere Boundary(LAB) in P-RFs of individual stations. Migrated Pls phase indicated a depth of 60–70 km for LAB. Inversions/comparisons of P-RF(Pls phase) and S-RF(Slp phase) waveforms together with Ps and Sp imaging for the crust and upper mantle structure further confirmed this result. P-RF and S-RF migrated images exhibit that a flat LAB is positioned at the depth of 60–70 km spreading along the profile, whereas a distinct structural change of lithospheric base appears at the Min River estuary. Both Ps and PpPs migrated images of P-RFs present an abrupt Moho drop across the Min River fault from south to north, which is consistent with previous result obtained from deep seismic sounding. By taking into consideration other geological and geophysical features such as locally high anomalies of crustal Poisson's ratios and heat flow at the Min River estuary, we infer that the Min River fault penetrates down to the Moho and may, furthermore, interfere in the deeper lithospheric structure.