共查询到20条相似文献，搜索用时 31 毫秒
1.
A blind thrust fault with a unique strike, which is orthogonal to the strike of most tectonic structures in Taiwan, triggered the Jiashian earthquake on March 4, 2010 (M = 6.4; 22.96°N, 120.70°E). This study utilizes 100 global positioning system stations to examine changes of surface displacements during the Jiashian earthquake. We mitigate effects of short-term noise and long-term plate movements from surface displacement data using a frequency dependent filter via the Hilbert–Huang transform and compute the horizontal azimuth (i.e. GPS-azimuth) using residual data at the NS component relative to residual data at the EW component. Analytical results show that orientations of horizontal azimuths were aligned and orthogonal to the strike of the blind thrust fault. Meanwhile, inverse orientations are observed before and after the earthquake that agrees well with the seismic rebound theory. As stress disturbed on strata a few days before the earthquake, an impeded region can be clearly identified by disordered orientations of horizontal azimuths for anticipating the mainshock. These results provide an additional view to explore stress disturbance associated with earthquakes and offer more information to examine diverse models of tectonic evolution in this region. 相似文献
2.
The Mw 9.0 Tohoku-Oki earthquake that occurred off the Pacific coast of Japan on March 11, 2011, was followed by thousands of aftershocks, both near the plate interface and in the crust of inland eastern Japan. In this paper, we report on two large, shallow crustal earthquakes that occurred near the Ibaraki-Fukushima prefecture border, where the background seismicity was low prior to the 2011 Tohoku-Oki earthquake. Using densely spaced geodetic observations (GPS and InSAR datasets), we found that two large aftershocks in the Iwaki and Kita-Ibarake regions (hereafter referred to as the Iwaki earthquake and the Kita-Ibarake earthquake) produced 2.1 m and 0.44 m of motion in the line-of-sight (LOS), respectively. The azimuth-offset method was used to obtain the preliminary location of the fault traces. The InSAR-based maximum offset and trace of the faults that produced the Iwaki earthquake are consistent with field observations. The fault location and geometry of these two earthquakes are constrained by a rectangular dislocation model in a multilayered elastic half-space, which indicates that the maximum slips for the two earthquakes are 3.28 m and 0.98 m, respectively. The Coulomb stress changes were calculated for the faults following the 2011 Mw 9.0 Tohoku-Oki earthquake based on the modeled slip along the fault planes. The resulting Coulomb stress changes indicate that the stresses on the faults increased by up to 1.1 MPa and 0.7 MPa in the Iwaki and Kita-Ibarake regions, respectively, suggesting that the Tohoku-Oki earthquake triggered the two aftershocks, supporting the results of seismic tomography. 相似文献
3.
《Journal of Asian Earth Sciences》2011,40(6):620-626
The Mw 9.3 Sumatra earthquake of December 26, 2004 caused extensive coseismic displacements globally, measurements of which were made essentially using modern geodetic techniques. This earthquake induced considerable perturbation in stress distribution as far as ∼8000 km away from the epicenteral region, which is tending to relax to its normal rates as seen from postseismic transient deformation. The monitoring of crustal displacements from strategically located sites using GPS provides coseismic as well as postseismic deformation that facilitates the understanding of the fault geometry, elastic thickness, postseismic relaxation mechanisms, rheology and earthquake recurrence time interval.We investigated coseismic and postseismic GPS derived displacements in Indian region together with the GPS data collected from Andaman and Sumatra region. It is found that while EW displacements are significantly large in peninsular India, those in the region to the north of Central India Tectonic Zone (CITZ) are relatively small. We could delineate the postseismic transients from position time series and interpreted them in terms of viscoelastic relaxation. It is inferred that the postseismic deformation is characterized by a power-law viscoelastic flow in the mantle. In Indian peninsula region, the timescale parameter of the exponential decay (τ = 250 days) would require an extremely low viscosity for the upper mantle. Relying on the prevailing coseismic and postseismic displacement fields, the present study also reflects upon the contemporary litho-tectonics of the Indian sub-continent. 相似文献
4.
The Philippine Fault results from the oblique convergence between the Philippine Sea Plate and the Sunda Block/Eurasian Plate. The fault exhibits left-lateral slip and transects the Philippine archipelago from the northwest corner of Luzon to the southeast end of Mindanao for about 1200 km. To better understand fault slip behavior along the Philippine Fault, eight GPS surveys were conducted from 1996 to 2008 in the Luzon region. We combine the 12-yr survey-mode GPS data in the Luzon region and continuous GPS data in Taiwan, along with additional 15 International GNSS Service sites in the Asia-Pacific region, and use the GAMIT/GLOBK software to calculate site coordinates. We then estimate the site velocity from position time series by linear regression. Our results show that the horizontal velocities with respect to the Sunda Block gradually decrease from north to south along the western Luzon at rates of 85–49 mm/yr in the west–northwest direction. This feature also implies a southward decrease of convergence rate along the Manila Trench. Significant internal deformation is observed near the Philippine Fault. Using a two dimensional elastic dislocation model and GPS velocities, we invert for fault geometries and back-slip rates of the Philippine Fault. The results indicate that the back-slip rates on the Philippine Fault increase from north to south, with the rates of 22, 37 and 40 mm/yr, respectively, on the northern, central, and southern segments. The inferred long-term fault slip rates of 24–40 mm/yr are very close to back-slip rates on locked fault segments, suggesting the Philippine Fault is fully locked. The stress tensor inversions from earthquake focal mechanisms indicate a transpressional regime in the Luzon area. Directions of σ1 axes and maximum horizontal compressive axes are between 90° and 110°, consistent with major tectonic features in the Philippines. The high angle between σ1 axes and the Philippine Fault in central Luzon suggests a weak fault zone possibly associated with fluid pressure. 相似文献
5.
The Aegean region including western Turkey, mainland Greece, and the Hellenic Arc is the most seismological and geodynamical active domain in the Alpine Himalayan Belt. In this study, we processed 3 years of survey-mode GPS data and present the analysis of a combination of geodetic and seismological data around Izmir, which is the third most populated city in Turkey. The velocities obtained from 15 sites vary between 25 mm/yr and 28 mm/yr relative to the Eurasian plate. The power law exponent of earthquake size distribution (b-value) ranges from 0.8 to 2.8 in the Izmir region between 26.2°E and 27.2°E. The lowest b-value zones are found along Karaburun Fault (b = 0.8) and, between Seferihisar and Tuzla Faults (b = 0.8). A localized stress concentration is expected from numerical models of seismicity along geometrical locked fault patches. Therefore, areas with lowest b-values are considered to be the most likely location for a strong earthquake, a prediction that is confirmed by the 2005 Mw = 5.9 Seferihisar earthquake sequences, with epicentres located to the south of the Karaburun Fault. The north–south extension of the Izmir area is corroborated by extension rates up to 140 nanostrain/yr as obtained from our GPS data. We combined the 3-year GPS velocity field with the published velocity field to determine the strain rate pattern in the area. The spatial distribution of b-value reflects the normal background due to the tectonic framework and is corroborated by the geodetic data. b-Values correlate with strain pattern. This relationship suggests that decrease of b-values signifies accumulating strain. 相似文献
6.
Evidence of precursor phenomena in the Kobe earthquake obtained from atmospheric radon concentration
《Applied Geochemistry》2006,21(6):1064-1072
Atmospheric 222Rn concentrations were determined over a 10a period, which included the date of the Kobe, Japan earthquake, on January 17th 1995. It was found that the seismically related 222Rn anomaly was higher than the 99% confidence limits for the residual value of atmospheric 222Rn which had been observed 2 months before. The residual 222Rn concentration, in which residual values of the daily minimum are the difference between each normal 222Rn concentration (calculated from January 1984 to December 1993) and the daily minimum 222Rn concentration (January 1994 to January 1995), was calculated by applying the exponential smoothing method to the residual values for each day. It was found that the fluctuations of the residual values can be fitted very well to a log-periodic oscillation model. The real residual values stopped increasing at 1994.999 (December 31st 1994), which corresponds with the critical point (tc) of best fit model. This anomalous 222Rn variation can be seen as the result of local stresses, not primary stresses which directly lead to the Kobe earthquake. On the other hand, when the critical exponent (z) and the radial frequency (ω) of the model were simultaneously fixed 0.2 ⩽ z ⩽ 0.6 and 6 ⩽ ω ⩽ 12, tc (critical point) was between January 13th 1995 and January 27th 1995. The Kobe earthquake occurrence date (January 17th 1995) is within this range. Therefore this anomalous 222Rn variation can also be seen as the result of primary stresses which possibly led to the Kobe earthquake. There is a distinct possibility that similar statistical oscillations will be detected in other measurements such as microseismicity, tectonic strain, fluctuation in the ground level, or changes in groundwater elevations and composition. 相似文献
7.
To investigate subsurface structure and seismogenic layers, 3D velocity inversion was carried out in the source zone of 1905 Kangra earthquake (M8.0) in the northwestern Himalaya. P-wave and S-wave phase data of 159 earthquakes recorded by a network of 21 stations were used for this purpose. Inverted velocity tomograms up to a depth range of 18 km show significant variations of 14% in Vp and Vs and 6% in the Vp/Vs across the major tectonic zones in the region. Synthesis of seismicity pattern, velocity structure, distinctive focal mechanisms coupled with nature of stress distribution allows mapping of three different source regions that control regional seismotectonics. Accumulating strains are partly consumed by sliding of Chamba Nappe to the southwest through reverse-fault movements along Chamba/Panjal/Main Boundary Thrusts. This coupled with normal-fault type displacements along Chenab Normal Fault in the north account for low magnitude widespread seismicity in upper 8–10 km of the crust. At intermediate depths from 8 to 15 km, adjusting to residual compressive stresses, the detachment or lower end of the MBT slips to produce thrust dominated seismicity. Nucleation of secondary stresses in local NE–SW oriented structure interacts in complex manner with regional stresses to generate normal type earthquakes below the plane of detachment and therefore three seismic regimes at different depths produce intense seismicity in a block of 30 × 30 km2 centered NE to the epicenter of Kangra earthquake. 相似文献
8.
Static stress changes caused by megathrust slip of the 2011 Mw 9.0 Tohoku-Oki earthquake considerably affected the seismicity patterns in inland areas, resulting in the occurrence of numerous earthquakes along several active faults in Japan. On June 30, 2011, the Mj 5.4 central Nagano earthquake occurred at a shallow depth of 5 km, indicating the reactivation of the Gofukuji fault in Central Japan. This study was undertaken to elucidate spatial and temporal changes of 3He/4He ratios around a source region before and after an inland earthquake using both existing and new and helium isotope data from hot spring and drinking water wells. Gas samples near the Gofukuji fault and its surrounding active faults are characterized by an increase in postseismic 3He/4He ratios. In contrast, the postseismic ratios decreased by up to about 30% away from the mainshock epicenter. Episodic faulting could either release stored crustal (radiogenic) helium from host rocks, or enhance the transfer of mantle volatiles through permeable fault zones, such that subsequent fluid flow near to the source region could then explain the spatio-temporal variations in 3He/4He ratios. 相似文献
9.
Esmaeil Shabanian Olivier Bellier Mohammad R. Abbassi Lionel Siame Yassaman Farbod 《Tectonophysics》2010,480(1-4):280-304
NE Iran, including the Kopeh Dagh and Allah Dagh-Binalud deformation domains, comprises the northeastern boundary of the Arabia–Eurasia collision zone. This study focuses on the evolution of the Plio-Quaternary tectonic regimes of northeast Iran. We present evidence for drastic temporal changes in the stress state by inversion of both geologically and seismically determined fault slip vectors. The inversions of fault kinematics data reveal distinct temporal changes in states of stress during the Plio-Quaternary (since ～ 5 Ma). The paleostress state is characterized by a regional transpressional tectonic regime with a mean N140 ± 10°E trending horizontal maximum stress axis (σ1). The youngest (modern) state of stress shows two distinct strike-slip and compressional tectonic regimes with a regional mean of N030 ± 15°E trending horizontal σ1. The change from the paleostress to modern stress states has occurred through an intermediate stress field characterized by a mean regional N trending σ1. The inversion analysis of earthquake focal mechanisms reveals a homogeneous, transpressional tectonic regime with a regional N023 ± 5°E trending σ1. The modern stress state, deduced from the youngest fault kinematics data, is in close agreement with the present-day stress state given by the inversions of earthquake focal mechanisms. According to our data and the deduced results, in northeast Iran, the Arabia–Eurasia convergence is taken up by strike-slip faulting along NE trending left-lateral and NNW trending right-lateral faults, as well as reverse to oblique-slip reverse faulting along NW trending faults. Such a structural assemblage is involved in a mechanically compatible and homogeneous modern stress field. This implies that no strain and/or stress partitioning or systematic block rotations have occurred in the Kopeh Dagh and Allah Dagh-Binalud deformation domains. The Plio-Quaternary stress changes documented in this paper call into question the extrapolation of the present-day seismic and GPS-derived deformation rates over geological time intervals encompassing tens of millions of years. 相似文献
10.
The Jiashian earthquake (ML 6.4) occurred on 4 March 2010. It was the largest inland event in southern Taiwan of 2010. The mainshock location was unexpected since it occurred in an area with relatively low background seismicity. In addition, reports of earthquake focal mechanisms do not fit with any known active fault geometry. In order to understand the origin of this earthquake, especially its rupture process, we perform a joint source inversion by using teleseismic body wave, GPS coseismic displacements and near field ground motion data. In this study, we considered a northwest–southeast trending fault with a northeast dip retrieved from GPS coseismic data and aftershocks distribution. To analyze the detailed slip distribution in space and time, we used near field 3D Green’s functions provided by spectral-element method and a full time–space inversion technique. We find a complex rupture process with several slip patches distributed inside two main asperities. The slip map reveals a mean slip of 12.9 cm for a maximum slip of 27.3 cm leading to a Mw 6.47 for this event. The rupture initiates in the deepest portion of the fault at 20 km depth, and propagated upward up to 2 km depth to form the two asperities. The source time function of this event revealed two pulses corresponding to the two asperities, for a total duration time of about 16 s. Most aftershocks occurred near the upper boundary of the deepest asperity while no aftershocks are located close to the shallowest one. We infer that the locations of these slip patches are related to the surrounding fault systems that may have restricted the rupture propagation during the earthquake. 相似文献
11.
18th June, 2010 5.9 Mw earthquake at North Andaman triggered along NW–SE pre-existing fault with reverse fault mechanism. Macroseismic survey and GPS geodesy reveal maximum damages following NE–SW trend due to normal fault mechanism. Coulomb stress modeling for post- and inter-seismic earthquakes after the 2004 mega-earthquake show different stages of fault segment linkage at North Andaman. The present earthquake has been explained as co-shock due to asiesmic soft linkage of fault propagation. 相似文献
12.
We estimate interseismic coupling on the subducting plate interface in the Tokai area, central Japan, by inverting two geodetic data sets. The data record surface motion between March 1996 to May 2000; one represents vertical motion deduced from the leveling observations and the other is the horizontal velocity field deduced from GPS observations. In the inversion, we employed the analytical solutions of surface displacement due to a triangular dislocation element embedded in a homogeneous elastic half space in order to represent the curved plate interface. The vertical data show that the most strongly coupled portion of the subduction interface is concentrated beneath Omaezaki Cape, while the horizontal data show strongest coupling in the shallower region of the subducting plate interface. The estimated maximum value of coupling from the horizontal data is 40 mm/year, while that from vertical data is 25 mm/year. 相似文献
13.
This study analyzed the rupture directivity of the 2011 Tohoku earthquake by using 100-s Rayleigh-wave travel-times, influenced by the finite source, to derive the fault parameters of the earthquake. The results demonstrated that the earthquake exhibited a slow rupture propagation with a rupture velocity of approximately 1.5–2.0 km/s and asymmetric bilateral faulting. The two rupture directions were N60°E and N127°E, with rupture lengths of approximately 276 km and 231 km, respectively. The rupture toward N60°E had a source duration of approximately 183 s, longer than that toward N127°E (approximately 156 s). Overall, the entire source duration of the earthquake faulting lasted approximately 183 s. Regarding historical seismicity in eastern Japan, the 2011 Tohoku earthquake not only ruptured a locked area in which large earthquakes have rarely occurred, but also ruptured the source regions of several historical earthquakes. With the exception of its slow rupture velocity and generation of a tsunami, the rupture features of the 2011 Tohoku earthquake were inconsistent with those of typical tsunami earthquakes. 相似文献
14.
《Gondwana Research》2014,25(1):204-213
Bounded by the western and eastern syntaxes, the Himalayan region has experienced at least five M ~ 8 earthquakes during a seismically very active phase from 1897 through 1952. However, there has been a paucity of M ~ 8 earthquakes since 1952. Examining of various catalogues and seismograms from the Gottingen Observatory, it is established that this quiescence of M ~ 8 earthquakes is real. While it has not been possible to forecast earthquakes, there has been a success in making a medium term forecast of an M 7.3 earthquake in the adjoining Indo-Burmese arc. Similarly we find that in the central Himalayan region, earthquakes of M > 6.5 have been preceded by seismic swarms and quiescences. In the recent past, based on GPS data, estimates have been made of the accumulated strains and it is postulated that a number of M ~ 8 earthquakes are imminent in the Himalayan region. We examine these estimates and find that while earthquakes of M ~ 8 may occur in the region, however, the available GPS data and their interpretation do not necessarily suggest their size and time of occurrence and whether an earthquake in a particular segment will occur sooner in comparison to that in the neighboring segment. We also comment on the inference of occurrence of M ~ 8 earthquakes based on M8 algorithm for the region. We conclude that while an M ~ 8 earthquake could occur any time anywhere in the Himalayan region, there is no indication as of now as to where and when it would occur. We impress on the need for preparedness to mitigate the pending earthquake disaster in the region. 相似文献
15.
《Quaternary Science Reviews》2005,24(12-13):1479-1498
Multiple peat-silt couplets preserved in tidal marsh sediment sequences suggest that numerous great plate boundary earthquakes caused the coast around Cook Inlet, Alaska, to subside over the past 3500 years. Field and laboratory analyses of the two youngest couplets record the well-documented earthquake of AD 1964 and the penultimate one, approximately 850 cal yr BP. Diatom assemblages from a range of modern day estuarine environments from tidal flat through salt marsh to acidic bog produce quantitative diatom transfer function models for elevation reconstructions based on fossil samples. Only nine out of 124 fossil assemblages analysed, including previously published data for the AD 1964 earthquake, have a poor modern analogue. Calibration of fossil samples indicate co-seismic subsidence of 1.50±0.32 m for AD 1964, similar to measurements taken after the earthquake, and 1.45±0.34 m for the ∼850 cal yr BP earthquake. Elevation standard errors for individual fossil samples range from ∼0.08 m in peat layers to ∼0.35 m in silt units. Lack of a chronology within fossil silt units prevents identification of changes in the rate of recovery and land uplift between the post-seismic and inter-seismic periods. However, preservation of multiple peat-silt couplets indicates no net emergence over multiple earthquake cycles. Glacio-isostatic movements from Little Ice Age glacier advance and retreat explains a ∼0.15 m relative sea-level oscillation recorded within the peat layer subsequently submerged as a result of the AD 1964 earthquake. Before both this and the ∼850 cal yr BP earthquake, diatom assemblages suggest pre-seismic relative sea-level rise of ∼0.12±0.13 m, representing possible precursors to great earthquakes. 相似文献
16.
Friction data used in modelling studies of subduction zone seismogenesis are often poorly representative of in situ conditions. We investigated the influence of in situ effective stresses and temperatures on the frictional properties of (simulated) fault gouges, prepared either from Nankai ODP material or illite shale, at sliding velocities approaching those relevant for earthquake nucleation and SSEs. Biaxial (double direct shear) experiments were performed at room temperature, normal stresses of 5–30 MPa, and sliding velocities of 0.16–18 μm/s. All materials exhibited velocity strengthening under these conditions, along with an increase in the friction coefficient and slip hardening rate with increasing normal stress. Illite gouge showed increased velocity strengthening towards higher normal stresses. The effect of temperature was investigated by means of ring shear experiments on illite gouge at 200–300 °C, an effective normal stress of 170 MPa, a pore-fluid pressure of 100 MPa and sliding velocities of 1–100 μm/s. These experiments showed a transition from velocity strengthening to velocity weakening at ∼250 °C. Our results provide a possible explanation for the updip seismogenic limit within subduction zone megathrusts and imply an enhanced tendency for earthquake nucleation and SSEs at low effective normal stresses. 相似文献
17.
Coseismic displacements play a key role in understanding earthquake dynamics. To derive displacement fields from optical and microwave remote sensing datasets, various methods are available. This study evaluated in detail the offset tracking technique on optical ASTER data for 2005 Kashmir earthquake. This technique required input parameters like resampling methods, correlator types, window sizes and step sizes. For accurate displacement field calculation, careful selection of these parameters is imperative which depends on the study area and dataset characteristics. In the study, we made relative comparisons of coseismic displacement fields calculated by using different combinations of input parameters. The results were validated by field based displacement data of vertical separation. Validation was based on the hypothesis that horizontal displacement component may also have vertical component contribution depending upon local characteristics of the fault. Validation results showed that general trend of the measured displacements was in agreement with the field data. Field measurements were bounded within the uncertainty limits of the technique however at some locations significant deviations were also observed. All the coseismic displacement results obtained by using different input parameter were within the uncertainty limit ±1/10 of the pixel size, except for window size 4 × 4 and 8 × 8. The measured component of the fault rupture for northwest of Muzaffarabad is irregular. It may be due to rugged topography as compared to southeast part. The measured fault rupture also coincided to surface rupture mapped in the field. Analysis of the results showed that in comparison to standard parameter set, defined in the literature (Sinc resampling method, Frequential correlator with window size 32 × 32 and step size 8), selection of resampling method and correlator type had no significant effect on the calculated displacement field. However, window size and step size were found to be the most important parameters of the technique. This study showed that changing the window size does not affect the detection limit of the technique but induces bias in the calculations. Step size is used in relation to window size, however relatively larger step size is found to provide better sensitivity and accuracy of the technique. 相似文献
18.
Zhi Wang Jian Wang Zhiliang Chen Yuping Liu Runqiu Huang Shunping Pei Qingzhi Zhang Wenqing Tang 《Gondwana Research》2011,19(1):202-212
Seismic imaging together with global positioning system (GPS) and crustal stress data analyses show that the Mw7.9 2008 Wenchuan earthquake occurred within a distinct area of high crustal stress (～ 17.5 MPa) and high Poisson's ratio (7–10%) anomalies centered on the Longmen-Shan (Shan means Mountain in Chinese) tectonic fault belt. Low P-wave and S-wave velocities in the southwest (SW) segment contrast with high-velocity anomalies in the central portion (CP) and northeast (NE) segment within the uppermost ～ 15 km depths along the tectonic fault belt, though a presumably ductile zone with low-velocity anomalies separates the CP and NE segment. The rupture initiated near the southwestern end of the CP at a zone of high Poisson's ratio (σ) which extends down into the lower crust. These low-velocity and high-σ anomalies immediately below the source hypocenter, together with the high crustal stress, indicate the presence of high-pressure fluids from the lower crust, which might have reduced the mechanical strength of the fractured rock matrix and triggered the earthquake. Our study suggests that the structural heterogeneity and high crustal stress played an important role in the nucleation of the Wenchuan earthquake and its rupture process. 相似文献
19.
The 10 June 2012 Mw 6.0 aftershock sequence in southwestern Anatolia is examined. Centroid moment tensors for 23 earthquakes with moment magnitudes (Mw) between 3.7 and 6.0 are determined by applying a waveform inversion method. The mainshock is a shallow focus strike-slip with reverse component event at a depth of 30 km. The seismic moment (Mo) of the mainshock is estimated as 1.28 × 1018 Nm and rupture duration of the Fethiye mainshock is 38 s. The focal mechanisms of the aftershocks are mainly strike-slip faulting with a reverse component. The geometry of the focal mechanisms reveals a strike-slip faulting regime with NE–SW trending direction of T-axis in the entire activated region. A stress tensor inversion of focal mechanism data is performed to obtain a more accurate picture of the Fethiye earthquake stress field. The stress tensor inversion results indicate a predominant strike-slip stress regime with a NW–SE oriented maximum horizontal compressive stress (SH). According to variance of the stress tensor inversion, to first order, the Fethiye earthquake area is characterized by a homogeneous interplate stress field. The Coulomb stress change associated with the mainshock and the largest aftershock are also investigated to evaluate any significant enhancement of stresses along the Gulf of Fethiye and surrounding region. Positive lobes with stress more than 0.4 bars are obtained, indicating that these values are large enough to increase the Coulomb stress failure towards NNW–SSE and E–W directions. 相似文献
20.
Bor-Shouh Huang Yi-Ling Huang Peih-Lin Leu Shiann-Jong Lee 《Journal of Asian Earth Sciences》2011,40(3):762-769
The rupture process of the disastrous Sumatra–Andaman earthquake of 26 December 2004 was analyzed by array processes for teleseismic P-waves recorded by a dense broadband seismic array in Taiwan with epicentral distances of close to 31°. The azimuthal variation from the BATS array center to both ends of the rupture fault is approximately 21°, which is larger than that reported previously for seismic arrays used to image the rupture process of this earthquake, thereby providing a high spatial resolution in studying the source rupture behavior. Two array-processing methods were used to analyze teleseismic P-wave trains. Both analyses were based on data recorded by a broadband network, covering a region of 200 × 400 km, with the aim of evaluating the rupture behavior of the earthquake. Consistent results from both analyses indicate that the earthquake had a rupture duration exceeding 500 s, with major asperities encountered at 80, 260, and 330 s after the initiation of rupturing. We traced the ruptured fault for more than 1200 km from the point of initial rupture. The average rupture velocity was approximately 3.0 km/s and the major northward rupture propagation began at 80 s after the initiation of rupturing. 相似文献