Vp and Vp/Vs structures in the crust and upper mantle of the Taiwan region, China

A tomographic study of the V _p and V _p/V _s structures in the crust and upper mantle beneath the Taiwan region of China is conducted by simultaneous inversion of P and S arrival times. Compared with the previous tomographic results, the spherical finite difference technique is suitable for the strong heterogeneous velocity structure, and may improve the accuracy in the travel time and three-dimensional ray tracing calculations. The V _p and V _p/V _s structures derived from joint inversion and the relocated earthquakes can provide better constraints for analyzing the lateral heterogeneity and deep tectonic characters in the crust and upper mantle. Our tomographic results reveal significant relations between the seismic wavespeed structure and the tectonic characters. In the shallow depth, sedimentary basins and orogen show distinct wavespeed anomalies, with low V _p, high V _p/V _s in basins and high V _p, low V _p/V _s in orogen. As the suture zone of Eurasian Plate and Philippine Sea Plate, Longitudinal Valley is characterized by a significant high V _p/V _s anomaly extending to the middle-lower crust and upper mantle, which reflects the impact of rock cracking, partial melting, and the presence of fluids. In the northeast Taiwan, the V _p, V _p/V _s anomalies and relocated earthquakes depict the subducting Philippine Sea Plate under the Eurasian Plate. The high V _p of oceanic plate and the low V _p, high V _p/V _s atop the subducted oceanic plate extend to 80 km depth. Along the east-west profiles, the thickness of crust reaches 60 km at the east of Central Range with eastward dipping trend, which reveals the eastward subduction of the thickened and deformed crust of the Eurasian continental plate. Supported by Knowledge Innovation Program of the Chinese Academy of Sciences (Grant No. KZCX3-SW-234-2), National Basic Research Program of China (Grant No. 2007CB411701), National High Technology Research and Development Program of China (Grant No. 2006AA09A101-0201) and National Natural Science Foundation of China (Grant Nos. 40804016, 40704013)     

Crust and uppermost mantle structure of the Ailaoshan-Red River fault from receiver function analysis

S-wave velocity structure beneath the Ailaoshan-Red River fault was obtained from receiver functions by using teleseismic body wave records of broadband digital seismic stations. The average crustal thickness, Vp/Vs ratio and Poisson’s ratio were also estimated. The results indicate that the interface of crust and mantle beneath the Ailaoshan-Red River fault is not a sharp velocity discontinuity but a characteristic transition zone. The velocity increases relatively fast at the depth of Moho and then increases slowly in the uppermost mantle. The average crustal thickness across the fault is 36―37 km on the southwest side and 40―42 km on the northeast side, indicating that the fault cuts the crust. The relatively high Poisson’s ratio (0.26―0.28) of the crust implies a high content of mafic materials in the lower crust. Moreover, the lower crust with low velocity could be an ideal position for decoupling between the crust and upper mantle.     

Crustal Structure Along Sunda-Banda Arc Transition Zone from Teleseismic Receiver Functions

We analyzed receiver function of teleseismic events recorded at twelve Indonesian-GEOFON (IA-GE) broadband stations using nonlinear Neighbourhood Algorithm (NA) inversion and H-k stacking methods to estimate crustal thickness, V _p /V _s ratios and S-wave velocity structure along Sunda-Banda arc transition zone. We observed crustal thickness of 34–37 km in Timor Island, which is consistent with the previous works. The thick crust (> 30 km) is also found beneath Sumba and Flores Islands, which might be related to the arc-continent collision causing the thickened crust. In Timor and Sumba Islands, we observed high V _p /V _s ratio (> 1.84) with low velocity zone that might be associated with the presence of mafic and ultramafic materials and fluid filled fracture zone. The high V _p /V _s ratio observed at Sumbawa and Flores volcanic Islands might be an indication of partial melt related to the upwelling of hot asthenosphere material through the subducted slab.     

Crustal velocity structure in the southern edge of the Central Alborz (Iran)

Inversion of local earthquake travel times and joint inversion of receiver functions and Rayleigh wave group velocity measurements were used to derive a simple model for the velocity crustal structure beneath the southern edge of the Central Alborz (Iran), including the seismically active area around the megacity of Tehran. The P and S travel times from 115 well-located earthquakes recorded by a dense local seismic network, operated from June to November 2006, were inverted to determine a 1D velocity model of the upper crust. The limited range of earthquake depths (between 2 km and 26 km) prevents us determining any velocity interfaces deeper than 25 km. The velocity of the lower crust and the depth of the Moho were found by joint inversion of receiver functions and Rayleigh wave group velocity data. The resulting P-wave velocity model comprises an upper crust with 3 km and 4 km thick sedimentary layers with P wave velocities (V_p) of ～5.4 and ～5.8 km s^?1, respectively, above 9 km and 8 km thick layers of upper crystalline crust (V_p ～6.1 and ～6.25 km s^?1 respectively). The lower crystalline crust is ～34 km thick (V_p ～ 6.40 km s^?1). The total crustal thickness beneath this part of the Central Alborz is 58 ± 2 km.     

CrustalP-wave velocity and velocity-ratio study in northeast India by a microearthquake survey

The upper crustal (20 km)P-wave velocity beneath the Shillong Plateau and Nowgong area has been studied by the time-distance plot method. TheP-arrival data of the shallow (20 km) microearthquakes from three temporary networks are used, and the average velocity is found to be 5.55 km/s. The velocity ratio (V _p /V _s ) for the upper crust (0–20 km) as well as for the lower crust (21–40 km) are determined by the Wadati-plot method and station-by-station method. The average value obtained by the two methods is compatible; theV _p /V _s ranges between 1.74 to 1.76. A generalized seismic velocity model of the area is suggested by this study, which has been very useful for microearthquake location.     

Three-dimensional velocity structure and seismicity of Mt. Etna Volcano, Italy

A new set of three-dimensional velocity models beneath Mt. Etna volcano is derived in the present work. We have used P- and S-wave arrivals from local earthquakes recorded at permanent and temporary seismic networks installed since 1980. A set of 1249 earthquakes recorded at more than four seismic stations was selected for traveltime inversion. The velocity models obtained by using different data selection criteria and parametrization display similar basic features, showing a high P-wave velocity at shallow depth in the SE quadrant, in close connection with a high gravimetric Bouguer anomaly. This area shares a low V_p/V_s ratio. High P-wave velocities and high V_p/V_s ratios are obtained along the central conduits, suggesting the presence of dense, intrusive magmatic bodies extending to a depth of about 20 km. The central intrusive core is surrounded by lower P-wave velocities. The relocated earthquake hypocenters also display the presence of an outward dipping brittle region, away from the central conduits, surrounding a ductile zone spatially related to the high P-wave velocity anomalies located in proximity to the central craters.     

用接收函数方法研究云南及其邻区地壳上地幔结构

利用云南和中国地震台网30个台站记录的远震资料,采用接收函数扫描法和线性反演方法对云南及其邻区的壳幔速度结构进行了研究,获得了研究区内地壳厚度、V_p/V_s以及壳幔速度的分布特征.利用接收函数扫描（H-k）法和线性反演方法获取的台站下方地壳厚度结果表明,研究区地壳厚度变化剧烈.速度结构研究结果表明,红河断裂以西的腾冲-保山地块和思茅地块,以及南华板块北部地区的台站（如攀枝花、丽江、东川、永胜等）下方地壳均存在下地壳低速层,且具有高的地壳平均V_p/V_s值.这些不仅暗示研究区的下地壳低速异常可能为高温甚至高温导致的部分熔融所致,同时,也意味着该区下地壳的物质易于发生塑性流动,为地壳的变形和增厚创造了条件.红河断裂带作为云南地区的一个主要边界断裂,其西侧地区地壳厚度变化较东侧剧烈,另一方面,红河断裂西侧的平均地壳V_p/V_s值较其东侧要高.综合前人关于研究区岩石圈速度结构、地热流值、重力场和上地幔各向异性等地球物理场的研究结果.我们推断,现今的思茅块体和保山-腾冲块体在大地构造上应归属于冈瓦纳板块.     

Migration of seismic activity during the 1998 volcanic unrest at Iwate volcano,northeastern Japan,with reference to P and S wave velocity anomaly and crustal deformation

We describe the seismicity at Iwate volcano, northeastern Japan, during the volcanic unrest of 1998 with reference to a three-dimensional P and S wave velocity model from tomographic analysis. The abnormal seismic activity beneath Iwate volcano started under the caldera in February, 1998 and migrated westward in the period February to August, 1998. Previous geodetic modeling [Sato and Hamaguchi, Chikyu Monthly 21 (1999) 312–317] suggested the growth of a dike in the time of the seismic activity. Comparing the seismicity and dike extension with the tomographic images of the P and S wave velocity structure, we find that the trace of the growing dike coincides with the region of the high V_p and high V_p/V_s ratio beneath the volcano. The seismic and geodetic data are consistent with an intrusion of magma or other fluid under the caldera in 1998. Another pressure source causing the predominant crustal deformation at Iwate volcano was detected from geodetic data, which was located in the region with high V_p/V_s ratio under the western end of the volcano through the period from February to August. It is suggested that the activation of the point pressure source probably associated with the inflation of a hot fluid reservoir relate to a geothermal region adjacent to the western edge of the volcano.     

Fine S wave velocity structure beneath Iwate volcano, northeastern Japan, as derived from receiver functions and travel times

A genetic algorithm inversion of receiver functions derived from a dense seismic network around Iwate volcano, northeastern Japan, provides the fine S wave velocity structure of the crust and uppermost mantle. Since receiver functions are insensitive to an absolute velocity, travel times of P and S waves propagating vertically from earthquakes in the subducting slab beneath the volcano are involved in the inversion. The distribution of velocity perturbations in relation to the hypocenters of the low-frequency (LF) earthquakes helps our understanding of deep magmatism beneath Iwate volcano. A high-velocity region (dV_S/V_S=10%) exists around the volcano at depths of 2–15 km, with the bottom depth decreasing to 11 km beneath the volcano’s summit. Just beneath the thinning high-velocity region, a low-velocity region (dV_S/V_S=−10%) exists at depths of 11–20 km. Intermediate-depth LF (ILF) events are distributed vertically in the high-velocity region down to the top of the low-velocity region. This distribution suggests that a magma reservoir situated in the low-velocity region supplies magma to a narrow conduit that is detectable by the hypocenters of LF earthquakes. Another broad low-velocity region (dV_S/V_S=−5 to −10%) occurs at depths of 17–35 km. Additional clusters of deep LF (DLF) events exist at depths of 32–37 km in the broad low-velocity zone. The DLF and ILF events are the manifestations of magma movement near the Moho discontinuity and in the conduit just beneath the volcano, respectively.     

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.     

Flat-Slab Thermal Structure and Evolution Beneath Central Mexico

Recent seismic and magnetotelluric experiments, aimed at better characterizing the shape and state of the subducting slab and continental crust beneath Central Mexico, exposed significant differences with conclusions of previous studies. A new slab geometry is revealed in which the subducting Cocos slab is perfectly flat between 120 to 290?km from the trench, after which it plunges into the asthenosphere at a dip angle of ~65°, in sharp contrast with the previously proposed ~20° dip angle. Seismic tomography studies show negative P-wave velocity anomalies (?2 to ?4%) in the mantle wedge beneath the Mexican Volcanic Belt, and positive anomalies (+2 to +3%) for the subducted Cocos slab. Magnetotelluric experiments exposed a very low-resistivity area (1?C10? ??m) located within the continental crust just below the Mexican Volcanic Arc. Finally, several spots of non-volcanic tremors (NVTs) have been recorded inside the continental crust above the flat-slab segment. While all these experiments provide a better picture of the subduction system beneath Central Mexico, several key processes need further investigation. In this study, we take advantage of these new observations to better constrain the thermal structure beneath Central Mexico. Two different thermal models are computed for a mantle potential temperature (T _p) of 1,350 and 1,450°C, respectively. The new thermal structures are then converted into P-wave velocity anomalies and compared with the observed V _p anomalies. We found that a T _p of 1,450°C produced larger V _p anomalies that do not fit the observations. However, using a T _p of only 1,350°C, our predicted V _p anomalies are positive (+2 to +3%) for the cold slab and negative (?2 to ?4%) in the mantle wedge. These V _p estimates are consistent with the observed seismic tomography from P-wave arrivals, and therefore we conclude that a T _p of 1,350°C is a better estimate for the mantle potential temperature beneath Central Mexico. The new thermal model, in conjunction with phase diagrams for sediments, hydrated basalt and lithospheric mantle, have been used to estimate the amount and location of fluids released from the subducting Cocos slab. Several dehydration pulses have been identified along the slab interface where most of the fluids stored in sediments and oceanic crust are released into the overlying continental crust above the flat-slab. We found a good correlation between the pattern of these dehydration pulses and the location of NVTs, suggesting that slab dehydration is responsible for triggering the tremors. We suggest that NVT bursts localized above the flat slab segment represent the manifestation of ongoing continental crust hydration and weakening, a process that has been going on since 15?Ma ago when the Cocos slab entered into a flat-slab regime. Such continuous weakening would have reduced the suction forces that kept the slab in a flat regime in the last 15?Ma, allowing the slab to easily roll back. The continuous low-resistivity region recorded beneath the volcanic front in Central Mexico might represent the evidence of slab dehydration and crust weakening over time.     

V p /V s anomalies in dilatant rock samples

Summary In a series of triaxial experiments we have measuredV _p ,V _s and volumetric strain simultaneously in dilating dry and saturated rocks. For the first time these data permit quantitative comparison of seismic velocities or their ratio and dilatant volumetric strain. In air-dry samplesV _p /V _s decreases by a few per cent at strains of 10^–3; in saturated materials with high pore pressure,V _p /V _s increases by a comparable amount. Decreases in seismic velocity ratio are difficult to generate in initially saturated rocks even with low pore pressures and at strain rates of 10^–4/sec. A liquid-vapor transition will not produce a significant drop inV _p /V _s . If dilatancy and fluid flow are responsible for seismic travel time anomalies prior to earthquakes, our results suggest that such anomalies will occur only in regions where pore fluid source to sink dimensions are of the order of 10 km or more, or in regions where the rocks are not saturated to begin with.     

Seismic evidence for stepwise thickening of the crust across the NE Tibetan plateau

We performed a receiver function analysis on teleseismic data recorded along two 550 km-long profiles crossing the northeastern Tibetan plateau. Results from time to depth migration, grid-search V_p/V_s determination and simulated annealing inversion of waveforms, reveal that the crust thickens from ∼50 km near the northern edge of the plateau to ∼80 km south of the Jinsha suture in the Qiang Tang block. Crustal thickening occurs in staircase fashion with steps located beneath the main, reactivated sutures. The V_p/V_s ratio, close to the global continental average does not suggest widespread partial melting but rather a more usual separation between an upper felsic and a lower mafic part within the northeastern Tibetan crust.     

Nature of the upper crust beneath central Tibet

A wide-angle reflection/refraction seismic line obtained through Sino-French joint research near the Bangung-Nujiang Suture in Tibet [1] yielded crustal velocity estimates for both P and S waves. AV_p/V_s ratio of 1.65 was derived fo waves propagating in the upper ( < 25 km) crust. According to the high-pressure, high-temperature velocity measurements of Kern [2] and others, such a ratio implies the presence of quartz-rich rocks, e.g., granites and granitic gneisses, in the upper Tibetan crust. Because of the high geothermal gradient in Tibet ( > 50°C/km near the surface), a shallow low-velocity zone is expected to form in a granitic upper crust for both P and S waves. A two-dimensional ray-tracing study indicates that the observed first P and S waves [1] can be interpreted as diving waves above a low-velocity zone where a positive gradient exists. The low-velocity zone commences at a depth of 10 km and extends possibly to 20 km or deeper.     

Modelling shear waves in OBS data from the Vøring basin (Northern Norway) by 2-D ray-tracing

Three component recordings from an array of five ocean bottom seismographs in the northwestern part of the Vøring basin have been used to obtain a 2-D shear-wave (S-wave) velocity-depth model. The shear waves are identified by means of travel-time differences compared to the compressional (P) waves, and by analyzing their particle motions. The model has been obtained by kinematic (travel-time) ray-tracing modelling of the OBS horizontal components.The shear-wave modelling indicates that mode conversions occur at several high velocity interfaces (sills) in the 4–10 km depth range, previously defined by a compressional-wave velocity-depth model using the same data set.An averageV _p /V _s ratio of 2.1 is inferred for the layers above the uppermost sill, indicative of both poorly consolidated sediments and a low sand/shale ratio. A significant decrease in theV _p /V _s ratio (1.7) below the first sill may in part be atributed to well consolidated sediments, and to a change in lithology to more sandy sediments. This layer is interpreted to lie within the lower Cretaceous sequence. At 5–10 km depthV _p /V _s ratios of 1.85 indicate a lower sand/shale ratio consistent with the expected lithologies. The averageV _p /V _s ratio inferred for the crust is 1.75, which is consistent with values obtained north of Vøring, in the Lofoten area. An eastward thinning of the crystalline basement is supported by the shear-wave modelling.     

Roughness effect on the correction factor of surface velocity for rill flows

Flow velocity is one of the most important hydrodynamic variables for both channelized (rill and gullies) and interrill erosive phenomena. The dye tracer technique to measure surface flow velocity V_s is based on the measurement of the travel time of a tracer needed to cover a known distance. The measured V_s must be corrected to obtain the mean flow velocity V using a factor α_v = V/V_s which is generally empirically deduced. The V_s measurement can be influenced by the method applied to time the travel of the dye-tracer and α_v can vary in different flow conditions. Experiments were performed by a fixed bed small flume simulating a rill channel for two roughness conditions (sieved soil, gravel). The comparison between a chronometer-based (CB) and video-based (VB) technique to measure V_s was carried out. For each slope-discharge combination, 20 measurements of V_s, characterized by a sample mean V_m, were carried out. For both techniques, the frequency distributions of V_s/V_m resulted independent of slope and discharge. For a given technique, all measurements resulted normally distributed, with a mean equal to one, and featured by a low variability. Therefore, V_m was considered representative of surface flow velocity. Regardless of roughness, the V_m values obtained by the two techniques were very close and characterized by a good measurement precision. The developed analysis on α_v highlighted that it is not correlated with Reynolds number for turbulent flow regime. Moreover, α_v is correlated neither with the Froude number nor with channel slope. However, the analysis of the empirical frequency distributions of the correction factor demonstrated a slope effect. For each technique (CB, VB)-roughness (soil, gravel) combination, a constant correction factor was statistically representative even if resulted in less accurate V estimations compared to those yielded by the slope-specific correction factor.     

Fracturing and hydrothermal alteration in normal fault zones

Large normal fault zones are characterized by intense fracturing and hydrothermal alteration. Displacement is localized in a slip zone of cataclasite, breccia and phyllonite surrounding corrugated and striated fault surfaces. Slip zone rock grades into fractured, but less comminuted and hydrothermally altered rock in the transition zone, which in turn grades abruptly into the wall rock. Fracturing and fluid flow is episodic, because permeability generated during earthquakes is destroyed by hydrothermal processes during the time between earthquakes.Fracture networks are described by a fracture fabric tensor (F). The permeability tensor (k) is used to estimate fluid transport properties if the trace of F is sufficiently large. Variations in elastic moduli and seismic velocities between fault zone and wall rock are estimated as a function of fracture density (). Fracturing decreases elastic moduli in the transition zone by 50–100% relative to the country rock, and similar or even greater changes presumably occur in the slip zone.P-andS-wave velocity decrease, andV _p /V _s increases in the fault zone relative to the wall rock. Fracture permeability is highly variable, ranging between 10^–13 m² and 10^–19 m² at depths near 10 km. Changes in permeability arise from variations in effective stress and fracture sealing and healing.Hydrothermal alteration of quartzo-feldspathic rock atT>300°C creates mica, chlorite, epidote and alters the quartz content. Alteration changes elastic moduli, but the changes are much less than those caused by fracturing.P-andS-wave velocities also decrease in the hydrothermally altered fault rock relative to the country rock, and there is a slight decrease inV _p /V _s , which partially offsets the increase inV _p /V _s caused by fracturing.Fracturing and hydrothermal alteration affect fault mechanics. Low modulus rock surrounding fault surfaces increases the probability of exceeding the critical slip distance required for the onset of unstable slip during rupture initiation. Boundaries between low modulus fault rock and higher modulus wall rock also act as rupture guides and enhance rupture acceleration to dynamic velocity. Hydrothermal alteration at temperatures in excess of 300°C weakens the deeper parts of the fault zone by producingphyllitic mineral assemblages. Sealing of fracture in time periods between large earthquakes generates pods of abnormally pressured fluid which may play a fundamental role in the initiation of large earthquakes.     

The Himalayan foreland basin crust and upper mantle

Modeling of multimode surface wave group velocity dispersion data sampling the eastern and the western Ganga basins, reveals a three layer crust with an average V_s of 3.7 km s^?1, draped by ～2.5 km foreland sediments. The Moho is at a depth of 43 ± 2 km and 41 ± 2 km beneath the eastern and the western Ganga basins respectively. Crustal V_p/V_s shows a felsic upper and middle crust beneath the eastern Ganga basin (1.70) compared to a more mafic western Ganga basin crust (1.77). Due to higher radiogenic heat production in felsic than mafic rocks, a lateral thermal heterogeneity will be present in the foreland basin crust. This heterogeneity had been previously observed in the north Indian Shield immediately south of the foreland basin and must also continue northward below the Himalaya. The high heat producing felsic crust, underthrust below the Himalayas could be an important cause for melting of midcrustal rocks and emplacement of leucogranites. This is a plausible explanation for abundance of leucogranites in the east-central Himalaya compared to the west. The uppermost mantle V_s is also significantly lower beneath the eastern Ganga basin (4.30 km s^?1) compared to the west (4.44 km s^?1).     

长白山-镜泊湖火山区地壳结构接收函数研究

利用71个远震的波形资料,用接收函数方法提取了布设在长白山—镜泊湖火山区的34个宽频带流动数字地震台站的接收函数,通过对接收函数反演,获得了台站下方的S波速度结构.研究结果表明,沈阳—敦化一线莫霍面深度32～33km,向西地壳厚度加厚,到长春附近地壳厚度约为36km.在天池火山口莫霍面深度为达38km,而镜泊湖火山口森林的莫霍面深度约为39km.总体看研究区的地壳厚度是南浅北深.长白山天池火山口附近地下10km左右有一明显的低速层存在;镜泊湖火山口森林附近30km也可能有低速体存在;研究发现莫霍面上S波速度梯度在火山口附近和远离火山口有明显区别.在火山口附近其莫霍面的S波速度梯度比非火山口地区的S波速度梯度明显小,说明火山口下与一般的地壳莫霍面结构有差别.研究发现沈阳—敦化一线两侧的莫霍面深度有较大变化,其位置与地表的敦化—密山断裂基本一致,说明敦化—密山断裂是研究区的一条非常重要的地质构造带.     

利用区域双差层析成像方法研究川滇南部地壳结构特征

本文联合使用云南、四川和贵州地震台网的85个地震台站在2008年1月—2017年12月期间记录的49130个地震、317366个初至Pg震相绝对到时数据和2674110条高精度的相对到时数据,采用区域双差地震层析成像方法联合反演了川滇南部地壳三维P波速度结构和39621个地震的震源参数,探究了川滇南部中下地壳流和腾冲火山区岩浆囊的分布特征.研究结果表明：（1）川滇南部上地壳的速度异常特征与地表地形密切相关;（2）小江断裂带的中下地壳存在一条绵延近二百多公里的低速异常结构,最南端受到红河断裂带的阻挡而终止于断裂带南段北侧,这可能是川滇南部的一条中下地壳流,低速异常结构在红河断裂带南段转而向南东流动反映了红河断裂带可能为川滇菱形块体的西南边界;（3）红河断裂带各段速度异常存在明显的差异,重定位后的震源分布显示红河断裂带中段和南段虽然不如北段地震活动强烈,但地震震源深度分布较北段深;（4）腾冲火山区西侧和北侧下方10~20 km深度范围内存在的低速异常体推测为通过怒江断裂带形成的岩浆通道从中地壳涌入上地壳的岩浆囊,可能反映了自更新世延续至今的以橄榄玄武岩和安山岩为主要岩性的壳内岩浆活动,持续的岩浆活动为地表热活动提供了主要动力.