首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
GeneratinggroundmotionbytwonewtechniquesofaddingharmonicwaveinthetimedomainandapproximatingtoresponsespectrumasawholeCHANG...  相似文献   

2.
The source parameters of the Bohai Sea earthquake, July 18, 1969 and Yongshan, Yunnan earthquake, May 11, 1974 were determined by full — wave theory synthetic seismograms of teleseismic P waves. P+pP+sP wereform were calculated with WKBJ approximation and real integral paths. One — dimensional unilateral, finite propagation source was also considered. By trail — and — error in comparing the theoretical seismograms with the observational ones of WWSSN stations, the source parameters were obtained as follow: for Bohai earthquake, φ=195°, δ=85°, λ=65°,M o=0.9×1019Nm,L=59.9km.V R=3.5km/s, ∧ R =160°; for Yongshan earthquake, φ=240°, δ=80°, ∧=150°,M o=1.3×1018Nm,L=48.8km,V R=3km/s, ∧ R =−10°, where φ is strike, δ dip angle, λ slip angle,M o seismic moment,L rupture length,V R rupture propagation speed. As III type fractures the faulting propagated along the fault planes, and ∧ R is the angle from the strike to the propagation direction. Yongshan earthquake showed complexity in its focal process, having four sub—ruptures during the first 60 seconds. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 1–8, 1991.  相似文献   

3.
Breakthrough point source model, extended earthquake source model is used to calculate more seismic source parameters in this paper. We express seismic source using higher degree moment tensors, to reduce a large number terms originally presenting in higher degree moment tensor representation, Haskell rupture model is used. We inverted the source parameters of Mani earthquake in Tibet using broad-band body wave of 32 stations of Global Seismograph Network (GSN), the results show that it is a strike-slip fault, rupture direction is 75° , rupture duration is 19 s, the fault plan is f =77° , d =88° , l =0° , the auxiliary plane is f =347° , d =90° , l =178° , and the fault dimension is 47 km′ 28 km. These results will give new quantitative data for earth dynamics and have practical meaning for seismic source tomography research.  相似文献   

4.
We develop a new method for the statistical estimation of the tail of the distribution of earthquake sizes recorded in the Harvard catalog of seismic moments converted to m W -magnitudes (1977–2004 and 1977–2006). For this, we suggest a new parametric model for the distribution of main-shock magnitudes, which is composed of two branches, the pure Gutenberg-Richter distribution up to an upper magnitude threshold m 1, followed by another branch with a maximum upper magnitude bound M max, which we refer to as the two-branch model. We find that the number of main events in the catalog (N = 3975 for 1977–2004 and N = 4193 for 1977–2006) is insufficient for a direct estimation of the parameters of this model, due to the inherent instability of the estimation problem. This problem is likely to be the same for any other two-branch model. This inherent limitation can be explained by the fact that only a small fraction of the empirical data populates the second branch. We then show that using the set of maximum magnitudes (the set of T-maxima) in windows of duration T days provides a significant improvement, in particular (i) by minimizing the negative impact of time-clustering of foreshock/main shock/aftershock sequences in the estimation of the tail of magnitude distribution, and (ii) by providing via a simulation method reliable estimates of the biases in the Moment estimation procedure (which turns out to be more efficient than the Maximum Likelihood estimation). We propose a method for the determination of the optimal choice of the T value minimizing the mean-squares-error of the estimation of the form parameter of the GEV distribution approximating the sample distribution of T-maxima, which yields T optimal = 500 days. We have estimated the following quantiles of the distribution of T-maxima for the whole period 1977–2006: Q 16%(M max) = 9.3, Q 50%(M max) = 9.7 and Q 84%(M max) = 10.3. Finally, we suggest two more stable statistical characteristics of the tail of the distribution of earthquake magnitudes: The quantile Q T (q) of a high probability level q for the T-maxima, and the probability of exceedance of a high threshold magnitude ρ T (m*)  = P{m k  ≥ m*}. We obtained the following sample estimates for the global Harvard catalog and The comparison between our estimates for the two periods 1977–2004 and 1977–2006, where the latter period included the great Sumatra earthquake 24.12.2004, m W  = 9.0 confirms the instability of the estimation of the parameter M max and the stability of Q T (q) and ρ T (m*) = P{m k  ≥ m*}.  相似文献   

5.
Dependence of the energy released during earthquake on ambient shear stress   总被引:1,自引:0,他引:1  
Starting with dislocation model, using the result of the fracture mechanics: the slip displacement at the crack tip is proportional to the length of the crack and the applied ambient shear stressτ 0 2 , we consider the dislocation in the earthquake to be the slip displacement at the crack tip and have obtained the analysis expresses of displacement and velocity pulse for the circular crack and have calculated the seismic wave energy radiated by earthquake. The seismic wave energyEM 0 τ 0 2 f(v) r , i. e.E is proportional to the seismic momentM 0 and the square of the ambient shear stressτ 0 2 and increases with the rupture velocityv r . In frequency domain, integrating the square of source velocity spectrum derived from our the scaling law model, we have also obtained the seismic wave energyE released by earthquake and earthquake radiated effficiencyη.EM 0 τ 0 2 also. If takingτ 0 = 10.0 MPa, E=4.79M 0. This result is consistent with the estimate by Vassiliou and Kanamori (1982). Theη=5.26%. The distribution of the seismic wave energy is that most of the energy contains in the frequency range between the first corner frequencyf c1 and thirdf c3, amount to 92.3% the energy in the rangef<f c1 is about 3.85% and 3.85% whenf>f c3. Thef c3 is about 8Hz forM ⩾ 6, thus most of radiated energy is below 2Hz. This phenomenon had been verified by Vassiliou Kanamori. Previous results show the energy radiated by earthquake to be strongly dependent on ambient shear stress. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 146–152, 1993. This work was supported by the Deutsche Forschungsgemeinschaft, Bonn, F. R. Germany. The support is grateful acknowledged. The authors are also grateful to Professor Klussmann and Mr. Hasthoff for their lots of help.  相似文献   

6.
Transport of non-ergodic solute plumes by steady-state groundwater flow with a uniform mean velocity, μ, were simulated with Monte Carlo approach in a two-dimensional heterogeneous and statistically isotropic aquifer whose transmissivity, T, is log-normally distributed with an exponential covariance. The ensemble averages of the second spatial moments of the plume about its center of mass, <S i i (t)>, and the plume centroid covariance, R i i (t) (i=1,2), were simulated for the variance of Y=log T, σ Y 2=0.1, 0.5 and 1.0 and line sources normal or parallel to μ of three dimensionless lengths, 1, 5, and 10. For σ Y 2=0.1, all simulated <S i i (t)>−S i i (0) and R i i (t) agree well with the first-order theoretical values, where S i i (0) are the initial values of S i i (t). For σ Y 2=0.5 and 1.0 and the line sources normal to μ, the simulated longitudinal moments, <S 11(t)>−S 11(0) and R 11(t), agree well with the first-order theoretical results but the simulated transverse moments <S 22(t)>−S 22(0) and R 22(t) are significantly larger than the first-order values. For the same two larger values of σ Y 2 but the line sources parallel to μ, the simulated <S 11(t)>−S 11(0) are larger than but the simulated R 11 are smaller than the first-order values, and both simulated <S 22(t)>−S 22(0) and R 22(t) stay larger than the first-order values. For a fixed value of σ Y 2, the summations of <S i i (t)>−S i i (0) and R i i , i.e., X i i (i=1,2), remain almost the same no matter what kind of source simulated. The simulated X 11 are in good agreement with the first-order theory but the simulated X 22 are significantly larger than the first-order values. The simulated X 22, however, are in excellent agreement with a previous modeling result and both of them are very close to the values derived using Corrsin's conjecture. It is found that the transverse moments may be significantly underestimated if less accurate hydraulic head solutions are used and that the decreasing of <S 22(t)>−S 22(0) with time or a negative effective dispersivity, defined as , may happen in the case of a line source parallel to μ where σ Y 2 is small.  相似文献   

7.
Transport of non-ergodic solute plumes by steady-state groundwater flow with a uniform mean velocity, μ, were simulated with Monte Carlo approach in a two-dimensional heterogeneous and statistically isotropic aquifer whose transmissivity, T, is log-normally distributed with an exponential covariance. The ensemble averages of the second spatial moments of the plume about its center of mass, <S i i (t)>, and the plume centroid covariance, R i i (t) (i=1,2), were simulated for the variance of Y=log T, σ Y 2=0.1, 0.5 and 1.0 and line sources normal or parallel to μ of three dimensionless lengths, 1, 5, and 10. For σ Y 2=0.1, all simulated <S i i (t)>−S i i (0) and R i i (t) agree well with the first-order theoretical values, where S i i (0) are the initial values of S i i (t). For σ Y 2=0.5 and 1.0 and the line sources normal to μ, the simulated longitudinal moments, <S 11(t)>−S 11(0) and R 11(t), agree well with the first-order theoretical results but the simulated transverse moments <S 22(t)>−S 22(0) and R 22(t) are significantly larger than the first-order values. For the same two larger values of σ Y 2 but the line sources parallel to μ, the simulated <S 11(t)>−S 11(0) are larger than but the simulated R 11 are smaller than the first-order values, and both simulated <S 22(t)>−S 22(0) and R 22(t) stay larger than the first-order values. For a fixed value of σ Y 2, the summations of <S i i (t)>−S i i (0) and R i i , i.e., X i i (i=1,2), remain almost the same no matter what kind of source simulated. The simulated X 11 are in good agreement with the first-order theory but the simulated X 22 are significantly larger than the first-order values. The simulated X 22, however, are in excellent agreement with a previous modeling result and both of them are very close to the values derived using Corrsin's conjecture. It is found that the transverse moments may be significantly underestimated if less accurate hydraulic head solutions are used and that the decreasing of <S 22(t)>−S 22(0) with time or a negative effective dispersivity, defined as , may happen in the case of a line source parallel to μ where σ Y 2 is small.  相似文献   

8.
The intrinsic dissipation and scattering attenuation in southwestern (SW) Anatolia, which is a tectonically active region, is studied using the coda waves. First the coda quality factor (Qc) assuming single scattering is estimated from the slope of the coda-wave amplitude decay. Then the Multiple Lapse Time Window (MLTW) analysis is performed with a uniform earth model. Three non-overlapping temporal data windows are used to calculate the scattered seismic energy densities against the source-receiver distances, which, in turn, are used to calculate separate estimates of the intrinsic and scattering factors. In order to explore the frequency dependency, the observed seismograms are band pass-filtered at the center frequencies of 0.75, 1.5, 3.0, 6.0 and 12.0. The scattering attenuation (Qs−1) is found lower than the intrinsic attenuation (Qi−1) at all frequencies except at 0.75 Hz where the opposite is observed. Overall the intrinsic attenuation dominates over the scattering attenuation in the SW Anatolia region. The integrated energy curves obtained for the first energy window (i.e., 0–15 s) are somewhat irregular with distance while the second (i.e., 15–30 s) and third (i.e., 30–45 s) data windows exhibit more regular change with distance at most frequencies. The seismic albedo B0 is determined as 0.61 at 0.75 Hz and 0.34 at 12.0 Hz while the total attenuation factor denoted by Le−1 changes in the range 0.034–0.017. For the source-station range 20–180 km considered the scattering attenuation is found strongly frequency dependent given by the power law Qs−1 = 0.010*f−1.508. The same relations for Qi−1, Qt−1 (total), Qc−1 and (expected) hold as Qi−1 = 0.0090*f−1.17, Qt−1 = 0.019*f−1.31, Qc−1 = 0.008*f−0.84 and respectively. Compared to the other attenuation factors Qc−1 and are less dependent on the frequency.  相似文献   

9.
An earthquake catalog derived from the detection of seismically-generated T-waves is used to study the time-clustering behavior of moderate-size (?3.0 M) earthquakes between 15 and 35°N along the Mid-Atlantic Ridge (MAR). Within this region, the distribution of inter-event times is consistent with a non-periodic, non-random, clustered process. The highest degrees of clustering are associated temporally with large mainshock-aftershock sequences; however, some swarm-like activity also is evident. Temporal fluctuations characterized by a power spectral density P(f) that decays as 1/fα are present within the time sequence, with α ranging from 0.12 to 0.55 for different regions of the spreading axis. This behavior is negligible at time scales less than ∼5×103 s, and earthquake occurrence becomes less clustered (smaller α) as increasing size thresholds are applied to the catalog. A power-law size-frequency scaling for Mid-Atlantic Ridge earthquakes also can be demonstrated using the distribution of acoustic magnitudes, or source levels. Although fractal seismic behavior has been linked to the structure of the underlying fault population in other environments, power-law fault size distributions have not been observed widely in the mid-ocean ridge setting.  相似文献   

10.
In the light of the single scattering model of coda originating from local earthquakes, and based on the aftershock coda registered respectively at the 4 short period stations installed near the foci shortly after theM7.6 Lancang andM7.2 Gengma earthquakes, this paper has tentatively calculated the rate of amplitude attenuation and theQ c-value of the coda in the Lancang and Gengma areas using a newly-founded synthetic determination method. Result of the study shows the rate of coda amplitude attenuation demonstrates remarkable regional differences respectively in the southern and northern areas. The southern area presents a faster attenuation (Q c=114), whereas the northern area shows a slower attenuation (Q c=231). The paper also discusses the reasons causing such differences. Result of the study also suggests a fairly good linear relation between the coda source factorA o(f) and the seismic moment and the magnitude. Using the earthquake scaling law, the following formulas can be derived: lgM 0=lgA 0(f)+17.6,M D=0.67lgA 0(f)+1.21 and logM 0=1.5M D+15.79. In addition, the rates of amplitude attenuationβ s andβ m are respectively calculated using the single scattering and multiple scattering models, and the ratioβ sm=1.20−1.50 is found for the results respectively from the two models. Finally, the mean free pathL of the S-wave scattering in the southern and northern areas are determined to be 54 km and 122 km respectively by the relations which can distinguish between the inherentQ i and scatteringQ s, testify to this areas having lowQ-values correspond to stronger scatterings. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, 71–82, 1992. This study is partly supported by the Seismological Science Foundation of the State Seismological Bureau of China, and the present English version of the paper is translated from its Chinese original by Wenyi Xia, Seismological Bureau of Yunnan Province.  相似文献   

11.
For the problem of matrix compaction and melt segregation a general mush continuity equation is derived, which explicitly expresses the coupling between the melt percolation and the inelastic matrix deformation and closes the governing equation set. Besides, a general equation is obtained, which describes the change in the volume of pore space due to all the possible reasons (inelastic matrix deformation, the phase transitions, and the advection of porosity by the matrix flow). The features of the isothermal melt segregation inside a partially molten zone are demonstrated using one-dimensional (1D) numerical solutions. It follows from the solutions that the pattern and the characteristic time of the melt segregation inside a partially molten zone of thickness L are controlled by the segregation parameter γ c = (L c )2, where the compaction length δ c = k0)η/(φ0μ) depends on the permeability, k, the value of characteristic porosity, φ0, and the viscosities of the matrix, η, and melt, μ. The solutions demonstrate that at any value of γ c , layers that are highly enriched in melt compared to the maximum initial porosity are formed in the upper part of the zone. At the same time, the evolution of the system and the segregation time differ considerably in the limits of γ c γ* and γ c γ*, where γ* depends on the boundary and initial conditions of the problem, and γ* is about 80 for the problem of melt segregation inside a partially molten zone with the maximum in the initial melt distribution located in the middle of the zone. At γ c γ*, which corresponds to the segregation of low-viscosity ultrabasic melts (kimberlites, carbonatites), all the melt accumulates to the roof of the zone, and the segregation time does not depend on the matrix permeability and melt viscosity and decreases with an increase in the thickness of the zone as L −1. The latter can be the reason for the formation of clusters of the same age and same composition eruptions characteristic of the kimberlite provinces. In the opposite limiting case, γ c γ*, the segregation time does not depend on the matrix viscosity and scales as L with a wave sequence forming in the upper part of the zone, which, probably, elucidates the origin of the rhythmical layering of the large tholeiitic basalt plutons.  相似文献   

12.
Riassunto L'Autore dimostra che, nel sistema di coordinate polari , , , si possono determinare un numeros di funzioni della sola variabile :Q 1,Q 3, ....Q 2s–1 tali che la sommatoria delleQ 2i–1/2i–1 rappresenti il potenzialeV di un geoide di rotazione. La condizione di armonicità determina ciascunaQ (che si riduce a un polinomio nelle potenze di sen ) a meno di una costante arbitraria; si dispone pertanto dis costanti che servono per soddisfare la natura dellaV sulla superficie del geoide. Come esempio l'Autore ha determinato la gravità sul geoide sferico, confermando i risultati delSomigliana, e su uno sferoide generico dove ha ritrovato la relazione diClairaut.
Summary The Author proofs that, in the system of polar coordinates , , , it is possible to determine a numbers of functions only of the variable :Q 1,Q 3 ....Q 2s–1 in such a way as to make the summatory of theQ 2i–1/2i–1 represent the potential function of a rotational geoid. The condition of harmonicity determines, saving an arbitrary constant, each of theQ which is reduced to a polynom developed by the sin powers; therefore one disposes of a number of constants to make use for satisfing theV on the geoid. To illustrate his theory the Author determines the gravity on the spherical geoid, thus confirmingSomigliana's formulas and on a spheroidal on which he pointed outClairaut's relations.
  相似文献   

13.
Numerical algorithms are suggested for reconstructing the potential field in the region Ω = {0 ≤ R 0ρR, 0 ≤ θ ≤ π, 0 ≤ φ ≤ 2π} from its values on the surface Γ1 = {R,0 ≤ α1 ≤ θ ≤ α2 ≤ π, 0 ≤ β1 ≤ φ ≤ β2 ≤ π}.  相似文献   

14.
Based on the daily turbulent heat fluxes and related meteorological variables datasets (1985–2006) from Objectively Analyzed air-sea Fluxes (OAFlux) Project of Woods Hole Oceanographic Institution (WHOI), characteristics of low-frequency oscillation intensity of air-sea turbulent heat fluxes over the northwest Pacific are analyzed by linear perturbation method and correlation analysis. It can be concluded that: 1) the distribution of low-frequency oscillation intensity of latent heat flux (LHF) over the northwest Pacific is mainly affected by that of low-frequency oscillation intensity of anomalous air-sea humidity gradient (Δq′) as well as mean air-sea humidity gradient (), while the distribution of low-frequency oscillation intensity of sensible heat flux (SHF) is mainly affected by that of low-frequency oscillation intensity of anomalous air-sea temperature gradient (ΔT′). 2) The low-frequency oscillation of turbulent heat fluxes over the northwest Pacific is the strongest in winter and the weakest in summer. And the seasonal transition of low-frequency oscillation intensity of LHF is jointly influenced by those of low-frequency oscillation intensity of Δq′, low-frequency oscillation intensity of anomalous wind speed (U′), and mean wind speed (Ū), while the seasonal transition of low-frequency oscillation intensity of SHF is mainly influenced by those of low-frequency oscillation intensity of ΔT′ and Ū. 3) Over the tropical west Pacific and sea areas north of 20°N, the low-frequency oscillation of LHF (SHF) is mainly influenced by atmospheric variables q a ′ (T a ′) and U′, indicating an oceanic response to overlying atmospheric forcing. In contrast, over the tropical eastern and central Pacific south of 20°N, q s ′ (T s ′) also greatly influences the low-frequency oscillation of LHF (SHF). Supported by National Natural Science Foundation of China (Grant No. 40675028) and National Basic Research Program of China (Grant No. 2006CB403600)  相似文献   

15.
Frequency analysis of the Sumatra-Andaman earthquake of 2004, one of the most significant and best-recorded earthquakes, is based on spectral seismograms obtained from their broadband seismograms. The Sumatra-Andaman earthquake is found to have a wide-range frequency content of P-wave radiation during the rupturing process. On the basis of stacking spectral seismograms we distinguished four time events of the rupturing process of a total length of about 540 s. The frequency, f max, is the highest for the first event (0.163 Hz in time interval 0–88 s), lowest for the second — which is the strongest (0.075 Hz in time interval 88–204 s). For third and fourth events frequencies are similar (0.089 and 0.082 Hz in time intervals 204–452 and 452–537 s, respectively). The frequency also shows an azimuthal dependence (±0.02 Hz). Azimuths for which the frequency, f max, has maximum and minimum values are 203–222° and 23–42°, respectively. These observations are discussed in relation to previously published papers on this topic.  相似文献   

16.
This paper presents low frequency wavenumbers in a seven-storey reinforced concrete building estimated from its recorded response to eleven earthquakes, one of which (1994 Northridge) caused visible structural damage, and two of which are its aftershocks. The wavenumbers, Ki,j(f), are estimated from pairs (i,j) of records at neighboring recording sites in the building, distributed vertically or horizontally. Changes in Ki,j(f) from one event to another are compared in the undamaged (lower) and in the damaged (upper) part of the building, with the aim to find whether trends in Ki,j(f) can indicate damage. The results suggest significant and permanent increase of the wavenumbers in the damaged parts for the 1994 Northridge earthquake and its aftershocks, which is not the case for the other events in the damaged parts, and for all eleven events in the undamaged parts of the building. This increase in wavenumbers in the damaged parts can be explained by reduced wave velocities through the damaged structural members, and by scattering of waves from the discontinuities created by the damage. It is concluded from this qualitative analysis that wavenumbers estimated from strong motion recordings in a building can indicate location of damage, and that it would be useful to refine further this method (extend it to higher frequencies, and add the capability to quantify the damage). However, this would require more dense strong motion instrumentation in buildings than currently available. Deployment of dense arrays in selected buildings would provide data for further work on this subject.  相似文献   

17.
ResearchoncharacteristicsofmagnitudestructureofearthquakesequencesPEI-QINGSUN(孙佩卿)QIN-ZULI(李钦祖)YING-HUADAI(戴英华)JUNZHAO(赵军)Se...  相似文献   

18.
The possibility-probability risk calculated using the interior-outer set model is referred to as fuzzy risk. A fuzzy expected value of the possibility-probability distribution is a set with E α(x) and [`(E)]a\bar E_\alpha (x) as its boundaries. The fuzzy expected values E α(x) and [`(E)]a\bar E_\alpha (x) of a possibility-probability distribution represent the fuzzy risk values being calculated. Under such an α level, three risk values can be calculated: conservative risk value, venture risk value and maximum probability risk value. As α adopts all values throughout the set [0, 1], it is possible to obtain a series of risk values. Therefore, the fuzzy risk can be a multi-valued risk or set-valued risk. Calculation of the fuzzy expected value of Yiwu city’s water resource risk has been performed based on the interior-outer set model. We can get a conservative risk value (R C ) of 800 mm for Yiwu city’s water resource risk, a venture risk value (R V ) of 1020 mm, and a maximum probability risk value (R M ) of 988 mm for the α = 0.1 level cut set.  相似文献   

19.
In this paper we suggest that conditional estimator/predictor of rockburst probability (and rockburst hazard, P T (t)) can be approximated with the formula P T (t) = P 1(θ 1)…P N (θ N P dyn T (t), where P dyn T (t) is a time-dependent probability of rockburst given only the predicted seismic energy parameters, while P i (θ i ) are amplifying coefficients due to local geologic and mining conditions, as defined by the Expert Method of (rockburst) Hazard Evaluation (MRG) known in the Polish mining industry. All the elements of the formula are (approximately) calculable (on-line) and the resulting P T value satisfies inequalities 0 ≤ P T (t) ≤ 1. As a result, the hazard space (0–1) can be always divided into smaller subspaces (e.g., 0–10−5, 10−5–10−4, 10−4–10−3, 10−3–1), possibly named with symbols (e.g., A, B, C, D, …) called “hazard states” — which saves the prediction users from worrying of probabilities. The estimator P T can be interpreted as a formal statement of (reformulated) Comprehensive Method of Rockburst State of Hazard Evaluation, well known in Polish mining industry. The estimator P T is natural, logically consistent and physically interpretable. Due to full formalization, it can be easily generalized, incorporating relevant information from other sources/methods.  相似文献   

20.
We analyze the strong motion accelerograms of the moderate (M w = 6.1), March 31, 2006, Darb-e-Astane earthquake of western Iran and also those of one of its prominently recorded, large (M w = 5.1) foreshock and (M w = 4.9) aftershock. (1) Using derived SH-wave spectral data, we first objectively estimate the parameters W o\mathit{\Omega} _{\rm o} (long period spectral level), f c (corner frequency) and Q(f) (frequency dependent, average shear wave quality factor), appropriate for the best-fit Brune ω  − 2 spectrum of each of these three events. We then perform a non-linear least square analysis of the SH-wave spectral data to provide approximate near-field estimates of the strike, dip, and rake of the causative faults and also the seismic moment, moment magnitude, source size, and average stress drop of these three events. (2) In the next step, we use these approximate values and an empirical Green’s function approach, in an iterative manner, to optimally model the strong ground motion and rupture characteristics of the main event in terms of peak ground acceleration/velocity/displacement and duration of ground shaking and thereby provide improved, more reliable estimates of the causative fault parameters of the main event and its asperities. Our near-field estimates for both the main moderate event and the two smaller events are in good conformity with the corresponding far-field estimates reported by other studies.  相似文献   

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

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