首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
We use GPS displacements collected in the 15 months after the 1999 Chi-Chi, Taiwan earthquake  ( M w 7.6)  to evaluate whether post-seismic deformation is better explained by afterslip or viscoelastic relaxation of the lower crust and upper mantle. We find that all viscoelastic models tested fail to fit the general features in the post-seismic GPS displacements, in contrast to the satisfactory fit obtained with afterslip models. We conclude that afterslip is the dominant mechanism in the 15-month period, and invert for the space–time distribution of afterslip, using the Extended Network Inversion Filter. Our results show high slip rates surrounding the region of greatest coseismic slip. The slip-rate distribution remains roughly stationary over the 15-month period. In contrast to the limited coseismic slip on the décollement, afterslip is prominent there. Maximum afterslip of 0.57 m occurs downdip and to the east of the hypocentral region. Afterslip at hypocentral depths is limited to the southern part of the main shock rupture, with little or no slip on the northern section where coseismic slip was greatest. Whether this results from along strike variations in frictional properties or dynamic conditions that locally favour stable sliding is not clear. In general, afterslip surrounds the area of greatest coseismic slip, consistent with post-seismic slip driven by the main shock stress change. The total accumulated geodetic afterslip moment is  3.8 × 1019 N m  , significantly more than the seismic moment released by aftershocks,  6.6 × 1018 N m  . Afterslip and aftershocks appear to have different temporal evolutions and some spatial correlations, suggesting that aftershock rates may not be completely controlled by the rate of afterslip.  相似文献   

2.
We present an approximate method to estimate the resolution, covariance and correlation matrix for linear tomographic systems Ax = b that are too large to be solved by singular value decomposition. An explicit expression for the approximate inverse matrix A − is found using one-step backprojections on the Penrose condition AA −≈ I , from which we calculate the statistical properties of the solution. The computation of A − can easily be parallelized, each column being constructed independently.
The method is validated on small systems for which the exact covariance can still be computed with singular value decomposition. Though A − is not accurate enough to actually compute the solution x , the qualitative agreement obtained for resolution and covariance is sufficient for many purposes, such as rough assessment of model precision or the reparametrization of the model by the grouping of correlating parameters. We present an example for the computation of the complete covariance matrix of a very large (69 043 × 9610) system with 5.9 × 106 non-zero elements in A . Computation time is proportional to the number of non-zero elements in A . If the correlation matrix is computed for the purpose of reparametrization by combining highly correlating unknowns x i , a further gain in efficiency can be obtained by neglecting the small elements in A , but a more accurate estimation of the correlation requires a full treatment of even the smaller A ij . We finally develop a formalism to compute a damped version of A −.  相似文献   

3.
The deformation at the core–mantle boundary produced by the 2004 Sumatra earthquake is investigated by means of a semi-analytic theoretical model of global coseismic and postseismic deformation, predicting a millimetric coseismic perturbation over a large portion of the core–mantle boundary. Spectral features of such deformations are analysed and discussed. The time-dependent postseismic evolution of the elliptical part of the gravity field ( J 2) is also computed for different asthenosphere viscosity models. Our results show that, for asthenospheric viscosities smaller than 1018 Pa s, the postseismic J 2 variation in the next years is expected to leave a detectable signal in geodetic observations.  相似文献   

4.
The Pisco earthquake ( M w 8.0; 2007 August 15) occurred offshore of Peru's southern coast at the subduction interface between the Nazca and South American plates. It ruptured a previously identified seismic gap along the Peruvian margin. We use Wide Swath InSAR observations acquired by the Envisat satellite in descending and ascending orbits to constrain coseismic slip distribution of this subduction earthquake. The data show movement of the coastal regions by as much as 85 cm in the line-of-sight of the satellite. Distributed-slip model indicates that the coseismic slip reaches values of about 5.5 m at a depth of ∼18–20 km. The slip is confined to less than 40 km depth, with most of the moment release located on the shallow parts of the interface above 30 km depth. The region with maximum coseismic slip in the InSAR model is located offshore, close to the seismic moment centroid location. The geodetic estimate of seismic moment is 1.23 × 1021 Nm ( M w 8.06), consistent with seismic estimates. The slip model inferred from the InSAR observations suggests that the Pisco earthquake ruptured only a portion of the seismic gap zone in Peru between 13.5° S and 14.5° S, hence there is still a significant seismic gap to the south of the 2007 event that has not experienced a large earthquake since at least 1687.  相似文献   

5.
Oxygen consumption by rotifers Macrotrachela musculosa and Trichotria truncata from Spitsbergen tundra (77°N) was measured using the method of Cartesian divers. The metabolic rate of M. musculosa was: 0.205 10−3mm3 02 per g 10−6 per hour at 2°C, 0.201 10−6mm3 at 6°C and 0.616 10−3mm3 02 per g 10−6 per hour at 10°C. The metabolic rate of Trichotria truncata at 6° was 0.103 10−3mm3 per g 10−6 per hour. The relation between body weight and oxygen consumption by M. musculosa at 2°C is expressed with the equation R = 0.18W0.67, with R – oxygen consumption in mm310−3 per individual per hour and W – wet weight of an animal in g 10−6.  相似文献   

6.
Summary. Asymptotic expressions for components of the electromagnetic field of a grounded electric dipole are considered for the model consisting of a thin surface-layer overlapping a stratified medium with a highly resistive screen on the roof. It is shown that the method of spatial derivatives makes it possible to obtain proper estimates of the impedance at distances of r ≥|λ0| from the nearest edge of the surface anomaly (|λ0| being the effective depth of the field penetration in the underlying section). the magnetotelluric methods allow one to obtain the true values of impedance, provided r ≥ max {|λ0|, |/( S −1+ Z −0|1/2} where S is the integrated conductivity of the surface layer, is the transverse resistance of the screen, and Z −0 is the Tikhonov—Cagniard impedance for the medium underlying the surface layer.  相似文献   

7.
Earthquake populations have recently been shown to have many similarities with critical-point phenomena, with fractal scaling of source sizes (energy or seismic moment) corresponding to the observed Gutenberg–Richter (G–R) frequency–magnitude law holding at low magnitudes. At high magnitudes, the form of the distribution depends on the seismic moment release rate M˙ and the maximum magnitude m max . The G–R law requires a sharp truncation at an absolute maximum magnitude for finite M˙ . In contrast, the gamma distribution has an exponential tail which allows a soft or 'credible' maximum to be determined by negligible contribution to the total seismic moment release. Here we apply both distributions to seismic hazard in the mainland UK and its immediate continental shelf, constrained by a mixture of instrumental, historical and neotectonic data. Tectonic moment release rates for the seismogenic part of the lithosphere are calculated from a flexural-plate model for glacio-isostatic recovery, constrained by vertical deformation rates from tide-gauge and geomorphological data. Earthquake focal mechanisms in the UK show near-vertical strike-slip faulting, with implied directions of maximum compressive stress approximately in the NNW–SSE direction, consistent with the tectonic model. Maximum magnitudes are found to be in the range 6.3–7.5 for the G–R law, or 7.0–8.2 m L for the gamma distribution, which compare with a maximum observed in the time period of interest of 6.1 m L . The upper bounds are conservative estimates, based on 100 per cent seismic release of the observed vertical neotectonic deformation. Glacio-isostatic recovery is predominantly an elastic rather than a seismic process, so the true value of m max is likely to be nearer the lower end of the quoted range.  相似文献   

8.
Summary. Fault zones in wet Westerly granite deformed at temperatures of 300° and 400°C require markedly lower shear stresses for sliding than when dry, even when the effective confining pressure is held constant between the wet and dry tests, provided that the strain rate is lower than 10−7s−1. The rate of strength reduction is enhanced by increasing the pore water pressure. The deformation rate is a power function of the applied stress where the stress exponent is approximately 7 for pore water pressure of 100 MPa and 21 for pore water pressure of 20 MPa.
The experimental results are extrapolated to conditions believed to occur at depths of 10 km along the San Andreas Fault Zone. It is suggested that for slow tectonic deformation at strain rates of 10−11 and 10−14s−1 the shear stress for sliding on faults in granite would be approximately 60 and 20 MPa, respectively, at pore water pressures equal to the hydrostatic head. Fluid overpressures of c. 0.8 lithostatic pressure are required to lower the shear stress for sliding into the 10 MPa range at the slower strain rate.  相似文献   

9.
Summary. The statistical capability of the m b: M s discriminant for the discrimination of earthquake and explosion populations is examined by application of discriminant functions to a group of 83 explosions and 72 earthquakes in Eurasia. Equations are derived for the probability that an event is an earthquake or an explosion. The positive sign of DIS in the decision index equation, DIS i = 34.3383 – 11.9569 mb t + 7.1161 M si , indicates that the event i is an earthquake. Its negative sign indicates that event i is an explosion. The probability of correct classification for an event, P i , is related to its DIS i value, by P i = [1-exp (DIS i )]−1, where a large, positive DIS indicates a high probability that an event is an earthquake and a large, negative DIS indicates a high probability that an event is an explosion. The discrimination line M s = 1.680 m b– 4.825, or m b= 0.595 M s+ 2.872 very successfully separates the explosion population from the earthquake population. The points on this line have an equal chance of being an earthquake or an explosion; moreover, for any event, the distance parallel to the M s-axis from the point representing that event in the m b: M s plane to this line is a measure of the probability for the correct classification of that event.  相似文献   

10.
The attenuation mechanism of seismic waves in northwestern Himalayas   总被引：3，自引：0，他引：3  
We analysed local earthquake waveforms recorded on a broad-band seismic network in northwestern Himalayas to compute the intrinsic and scattered attenuation parameters from coda waves. Similar to other tectonically active and heterogeneous regions, attenuation-frequency relation for western Himalaya is   Q −1 c = (113 ± 7)  f (1.01±0.05)  where   Qc   is the coda Q parameter. Intrinsic  ( Q −1 i )  and scattering  ( Q −1 s )  attenuations was separated using   Qc   and direct S -wave Q data  ( Qd )  . It is observed that estimated   Q −1 c   is close to   Q −1 i   and both of them are much larger than   Q −1 s   suggesting that coda decay is predominantly caused by intrinsic attenuation. At higher frequencies, both the attenuation parameters   Qc   and,   Qd   are similar indicating that coda is predominantly composed of back-scattered S waves at these frequencies.  相似文献   

11.
By inversion analysis of the baseline changes and horizontal displacements observed with GPS (Global Positioning System) during 1990–1994, a high-angle reverse fault was detected in the Shikoku-Kinki region, southwest Japan. The active blind fault is characterized by reverse dip-slip (0.7±0.2  m yr−1 within a layer 17–26  km deep) with a length of 208±5  km, a (down-dip) width of 9±2  km, a dip-angle of 51°±2° and a strike direction of 40°±2° (NE). Evidence from the geological investigation of subfaults close to the southwestern portion of the fault, two historical earthquakes ( M L=7.0, 1789 and 6.4, 1955) near the centre of the fault, and an additional inversion analysis of the baseline changes recorded by the nationwide permanent GPS array from 18 January to 31 December 1995 partially demonstrates the existence of the fault, and suggests that it might be a reactivation of a pre-existing fault in this region. The fact that hardly any earthquakes ( M L>2.0) occurred at depth on the inferred fault plane suggests that the fault activity was largely aseismic. Based on the parameters of the blind fault estimated in this study, we evaluated stress changes in this region. It is found that shear stress concentrated and increased by up to 2.1 bar yr−1 at a depth of about 20  km around the epicentral area of the 1995 January 17  Kobe earthquake ( M L=7.2, Japan), and that the earthquake hypocentre received a Coulomb failure stress of about 5.6 bar yr−1 during 1990–1994. The results suggest that the 1995  Kobe earthquake could have been induced or triggered by aseismic fault movement.  相似文献   

12.
We invert surface-wave and geodetic data for the spatio-temporal complexity of slip during the M w =8.1 Chile 1995 event by simulated annealing. This quasi-global inversion method allows for a wide exploration of model space, and retains the non-linearity of the source tomography problem. Complex source spectra are obtained from 5 to 45 mHz from first- and second-orbit fundamental-mode Rayleigh waves using an empirical Green's function cross-correlation technique. Coseismic displacement vectors were measured at 10 GPS sites near Antofagasta. They are part of a French-Chilean experiment which monitors the Northern Chile seismic gap. The spectra, together with the geodetic data, are inverted for the moment distribution on a 2-D dipping fault, under the physical constraints of slip positivity and causality. Marginal a posteriori distributions of the model parameters are obtained from several independently inverted solutions. In general, features of the slip model are well resolved. Data are well fitted by a purely unilateral southward rupture with a nearly uniform velocity around 2.5–3.0 km s−1, and a total duration of 65 s. Several regions of moment release were imaged, one near the hypocentre, a major one 80 km south of it and a minor one 160 km south of it. The major patch of moment release seemed to have propagated to relatively shallow depths near the trench, 100 km SSW of the epicentre. The region of major slip is located updip of the 1987, M w =7.5 earthquake, suggesting a causal relationship. Most of the slip occurred updip of the hypocentre (36 km), but the entire coupled plate interface (20–40 km) ruptured during the Chile 1995 event.  相似文献   

13.
Summary. During 1977 March and April, three Sacks-Evertson borehole dilatometers were installed at the ends of boreholes drilled into the sidewall of an experimental tunnel at a depth of 3.1 km in the ERPM gold mine near Johannesburg. In the following year coseismic strain changes ranging from 5 ± 10−10 to values exceeding 5 ± 10−6 were recorded for hundreds of mine tremors in the magnitude range -1 to 3.7 and at hypocentral distances of 50 m to about 2 km. Hypocentral coordinates and magnitudes were determined from seismograms recorded from an underground array of geophones. Amplitudes and polarities of the coseismic strain steps are generally in excellent agreement with theoretical expectations based on point-source dislocation theory; specifically, the strain steps are proportional to the seismic moment divided by the cube of hypocentral distance. At a strain level of 5 ± 10−9 or greater the tremors do not appear to be preceded by any short-term indications of instability even for tremors producing coseismic steps greater than 5 ± 10−6 and for which the strainmeters were within a source radius of the hypocentre. Continuous strain changes observed at the times when the mine excavation, at a distance of about 100 m, is extended are in good agreement with calculated changes based on the theory of elasticity. A similar calculation is consistent with post-seismic strain changes observed to follow some of the closer tremors. These post-seismic strains show a logarithmic dependence on time following the tremor and appear to be due to the interaction of a tremor with the adjacent mine excavation rather than to deformation within the actual seismic source region.  相似文献   

14.
Seismic sources with observable glut moments of spatial degree two   总被引：1，自引：0，他引：1  
Let ζΛ and r Λ. be the hypocentral position and time of an extended indigenous seismic source. Backus showed that the force moment tensors of the source, Γ( m +1, n )(ζΛ, r Λ), determine and are determined by the motion which the source produces. For small m + n , only the long-period motion is relevant. The glut moment tensor Λ( m,n ) (ζΛ, r Λ.) can be calculated uniquely from γ( m +1, n )(ζΛ r Λ) only if m = 0 or m = 1. The tensor G =Λ(2,0)(ζΛ) gives the spatial variance tensor WΛ of the source, and WΛ. roughly describes the size, shape and orientation of the source region. Therefore the failure of the observed F =Γ(3,0)(ζΛ) to determine G uniquely is of seismological interest. In the present paper we show that F determines G uniquely if we assume the source to be a simple straight line source (SSLS) or an ideal fault in an isotropic medium with isotropic prestress (IFIMIP). We give tests on F which determine whether it can come from a SSLS, from an IFIMIP or from a simple plane surface source (SPSS). If we assume the source to be a SPSS then knowing F and the fault plane determines G to within an unknown scalar multiple of a certain tensor tangent to the fault plane. Moreover F determines the fault plane uniquely unless F can come from a SSLS. If it can, then F determines this virtual source line uniquely, and F permits the fault plane to be any plane containing the virtual source line.  相似文献   

15.
In the austral summer seasons 2001/02 and 2002/03, Global Positioning System (GPS) data were collected in the vicinity of Vostok Station to determine ice flow velocities over Lake Vostok. Ten GPS sites are located within a radius of 30 km around Vostok Station on floating ice as well as on grounded ice to the east and to the west of the lake. Additionally, a local deformation network around the ice core drilling site 5G-1 was installed.
The derived ice flow velocity for Vostok Station is  2.00 m a−1± 0.01 m a−1  . Along the flowline of Vostok Station an extension rate of about 10−5 a−1 (equivalent to 1 cm km−1 a−1) was determined. This significant velocity gradient results in a new estimate of 28 700 years for the transit time of an ice particle along the Vostok flowline from the bedrock ridge in the southwest of the lake to the eastern shoreline. With these lower velocities compared to earlier studies and, hence, larger transit times the basal accretion rate is estimated to be 4 mm a−1 along a portion of the Vostok flowline. An assessment of the local accretion rate at Vostok Station using the observed geodetic quantities yields an accretion rate in the same order of magnitude. Furthermore, the comparison of our geodetic observations with results inferred from ice-penetrating radar data indicates that the ice flow may not have changed significantly for several thousand years.  相似文献   

16.
Summary. This note presents an exact analytical formula for determining the magnitude of coseismic surface volume change (δ V ) of earthquake faults in a half-space. For a Poisson solid, the formula is remarkably simple; δ V = M zz |8μ, where M zz is one of the moment tensor elements of the source. Maximum δ V values derive from dip slip on faults plunging 45°. For these events, surface volume changes of 0.0001 and 4.3 km3 can be expected for magnitude 5 and 8 earthquakes respectively. All of the coseismic surface volume change is recovered in the interseismic period through relaxation of the Earth and rebound of the surface. A useful rule of thumb for estimating the magnitude of vertical rebound in 45° dip slip events is δ h p=Δ s /24, where Δ s is the coseismic slip on the fault.  相似文献   

17.
Summary. Using nine IDA records for the Indonesian earthquake of 1977 August 19, we have formed an optimal linear combination of the records and have measured the frequency and Q of 0 S 0 and 1 S 0. The frequency was measured using the moment ratio method. The attenuation was measured by the minimum width method and by the time-lapse method. The frequency and attenuation were measured simultaneously by varying them to obtain a best fit to the data. A 2000-hr stack, the sum of nine individual records, for 0 S 0 gave a frequency of 0.814664 mHz±4 ppm. The values for the Q of 0 S 0 for the three different methods of measurement were 5600,5833 and 5700, respectively. The error in the estimates of Q -1 is about 5 per cent for the minimum power method. For 1 S 0 a 300-hr stack yielded a frequency of 1.63151 mHz±30 ppm. The values of Q for this mode were 1960, 1800 and 1850, respectively, with an error in Q -1 of about 12 per cent for the minimum power method.  相似文献   

18.
We evaluate the seismic moment–frequency relation for the Harvard catalogue in the period 1977–1994. This catalogue is composed of about 12 000 earthquakes. After selection of events in terms of depth and energy, we retain about 8000 data points. We estimate two parameters of the seismic moment distribution: the power exponent β and the cut-off value M m . The method used is a least-squares linear fit on a log–log scale performed over a range selected on the basis of the standard deviation from the histogram. The analysis is carried out for different subdivisions of the Earth in square grids of different sizes. Neither parameter exhibits a dependence on cell size, suggesting the universality of their values and the interpretation of the existence of a cut-off as a finite size effect linked to a finite catalogue length. The variations of the parameters are investigated as a function of time (duration of the catalogue) and versus the number of events used for building up the distribution. Again, β and M m do not depend on time, but M m depends on the number of events, reaching a stable value for N ≈ 1000. The only significant change in the parameters is observed for different values of M 0upper in the catalogue, revealing the existence of universality classes.  相似文献   

19.
Summary. We have implemented an algorithm which is based on Bailey's solution of the inverse problem of electromagnetic induction in the Earth. The study was motivated by recent determinations of very long period data and also benefited from recent redeterminations of high frequency data. The algorithm has been successfully tested to provide reliable estimates of conductivity down to a depth of 2000 km, using synthetic data in the period range from 4 days to 11 years. Smooth data sets, which are required for the inversion, were constructed from various sources. At a given depth, the range of inverted models is less than one order of magnitude. Due to the lack of high frequency data, the conductivity of the upper 600 km of the mantle, which is found to be of the order of 10−1Ω−1 m−1, may be overestimated. The algorithm performs well in the middle mantle, where conductivity rises steadily from 1 to 50 Ω−1 m−1. The lack of very low frequency data and limitations of the algorithm prevent one from obtaining meaningful estimates in the lower mantle. However, the study of the propagation of the late 1960s secular variation acceleration provides an estimate of the mean conductivity of the whole mantle. Thus, a complete mantle profile can be constructed. It is found that deep mantle conductivity probability does not exceed a few hundred Ω−1 m−1.  相似文献   

20.
Summary. The method proposed by Mendiguren to determine the source parameters from free oscillation data is applied to the 1970 July 31 deep Colombian earthquake. The results indicate a source propagating horizontally for about 150 km along the lithosphere and cutting across its width. The slab behaves as a guide for source propagation. The horizontal propagation velocity is determined as 3.8 km/s. The intensity of the source grew proportionately to the second power of the propagation distance. This rate of source intensity growth may be interpreted either by a fan-shaped fault model or by a cone-shaped volume source. The average slip and stress drop are estimated as 360 cm and 300 bar for the fault model. For the volume source model the transformational shear strain and stress are estimated as 11 × 10−5 and 160 bar. There is no evidence of a double couple radiation preceding the P origin time. It is shown that the isotropic and deviatoric components of the moment tensor cannot be uniquely resolved when only observations of a single mode are available. It is observed that, from a statistical basis, the available 0 Sn data for Colombian shock can be equally well explained by a pure deviatoric source model or by a source model including an isotropic component. Numerical experiments indicate that the inclusion of higher mode data does not change this situation. But, on the other hand, numerical experiments show that the available data and the scheme used for the inversion would not result in a solution including an artificial implosive component if the actual source were pure deviatoric. If the departure from a pure deviatoric source is produced by noise, it has to be non-random, as it could be produced by lateral heterogeneities not included in the inversion scheme.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号