首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
It is shown that when the travel-time curve of a refracted wave from a surface source is known and at least one of the following conditions is satisfied, i.e. when
  1. the travel-time curve of a wave reflected from a horizontal interface lying below the deepest low velocity layer is known, or
  2. the travel-time curve of a wave from a deep source situated below the deepest low velocity layer is known, or
  3. the measureH(u)=mes {z∶z≥0,v ?1 (z)≥u} is analytical in some segment [c, d], where \(0< c< d< \infty , c< a_n , H(a_n ) = \bar z_n ,\bar z_n\) is the depth of the lower end of the deepest low velocity layer and in the interval [c, ∞) an analytical functionH(u)) exists which providesH(u)≡H(u)) ifu∈[c, d], then (1) velocityv(z) outside the low velocity layers and (2) the measureH k (u)=mes {z∶z∈L k,v ?1 (z)≥u} for each low velocity layerL k,k=1, 2, ..., n, are defined unambiguously.
  相似文献   

2.
Using seismic data from regional earthquakes in Tibet recorded by the Hi-CLIMB experiment, Pn attributes are used to constrain the velocity gradient and attenuation structure of the Tibetan lithosphere under the Hi-CLIMB array. Numerical modeling is performed using the spectral-element method (SEM) for laterally varying upper-mantle velocity and attenuation, and the seismic attributes considered include the Pn travel-time, envelope amplitude, and pulse frequency. The results from the SEM modeling provide two alternative models for the upper-mantle beneath the Hi-CLIMB array in Tibet. The first model is derived from the 3D velocity model of Griffin et al. (Bull Seism Soc Am 101:1938–1947, 2011) with a constant upper-mantle velocity gradient, and laterally varying upper mantle attenuation. The second model has a laterally varying upper-mantle velocity gradient, and constant upper-mantle attenuation. In both cases, the Qiangtang terrane is distinguished from the Lhasa terrane by a change in Moho depth and upper-mantle velocities. The lower upper-mantle velocities, as well as higher Pn attenuation, suggest hotter temperatures beneath the Qiangtang terrane as compared to the Lhasa terrane. Although the fits to the Pn amplitude and pulse frequency data are comparable between the two models, the first model with the constant upper-mantle velocity gradient fits the travel times somewhat better in relation to the data errors.  相似文献   

3.
Retrospective data of monitoring under conditions of low seismic activity are used to identify free oscillations of the Earth, including the fundamental mode, the oscillation with a central period of 54 min (0 S 2 m multiplet), split into five lines with azimuthal numbers m = ?2, ?1, 0, 1, 2. It is shown that some lines of this oscillation are also recorded in atmospheric pressure variation spectra and group in ensembles of observations around frequencies predicted by the 0 S 2 m splitting theory. This phenomenon is discovered in data recorded both synchronously and in different time intervals. A causal relationship involved in the oscillation under study is determined on the basis of the examination of the direction of the acoustic energy flux. The energy flux in the region of the 0 S 2 m multiplet is shown to be directed from the Earth toward the atmosphere. This suggests that deep processes in the Earth are capable of exciting its upper shells.  相似文献   

4.
We estimate the corner frequencies of 20 crustal seismic events from mainshock–aftershock sequences in different tectonic environments (mainshocks 5.7 < M W < 7.6) using the well-established seismic coda ratio technique (Mayeda et al. in Geophys Res Lett 34:L11303, 2007; Mayeda and Malagnini in Geophys Res Lett, 2010), which provides optimal stability and does not require path or site corrections. For each sequence, we assumed the Brune source model and estimated all the events’ corner frequencies and associated apparent stresses following the MDAC spectral formulation of Walter and Taylor (A revised magnitude and distance amplitude correction (MDAC2) procedure for regional seismic discriminants, 2001), which allows for the possibility of non-self-similar source scaling. Within each sequence, we observe a systematic deviation from the self-similar \( M_{0} \propto \mathop f\nolimits_{\text{c}}^{ - 3} \) line, all data being rather compatible with \( M_{0} \propto \mathop f\nolimits_{\text{c}}^{ - (3 + \varepsilon )} \) , where ε > 0 (Kanamori and Rivera in Bull Seismol Soc Am 94:314–319, 2004). The deviation from a strict self-similar behavior within each earthquake sequence of our collection is indicated by a systematic increase in the estimated average static stress drop and apparent stress with increasing seismic moment (moment magnitude). Our favored physical interpretation for the increased apparent stress with earthquake size is a progressive frictional weakening for increasing seismic slip, in agreement with recent results obtained in laboratory experiments performed on state-of-the-art apparatuses at slip rates of the order of 1 m/s or larger. At smaller magnitudes (M W < 5.5), the overall data set is characterized by a variability in apparent stress of almost three orders of magnitude, mostly from the scatter observed in strike-slip sequences. Larger events (M W > 5.5) show much less variability: about one order of magnitude. It appears that the apparent stress (and static stress drop) does not grow indefinitely at larger magnitudes: for example, in the case of the Chi–Chi sequence (the best sampled sequence between M W 5 and 6.5), some roughly constant stress parameters characterize earthquakes larger than M W ~ 5.5. A representative fault slip for M W 5.5 is a few tens of centimeters (e.g., Ide and Takeo in J Geophys Res 102:27379–27391, 1997), which corresponds to the slip amount at which effective lubrication is observed, according to recent laboratory friction experiments performed at seismic slip velocities (V ~ 1 m/s) and normal stresses representative of crustal depths (Di Toro et al. in Nature in press, 2011, and references therein). If the observed deviation from self-similar scaling is explained in terms of an asymptotic increase in apparent stress (Malagnini et al. in Pure Appl Geophys, 2014, this volume), which is directly related to dynamic stress drop on the fault, one interpretation is that for a seismic slip of a few tens of centimeters (M W ~ 5.5) or larger, a fully lubricated frictional state may be asymptotically approached.  相似文献   

5.
The TKE dissipation rate in the northern South China Sea   总被引:1,自引:0,他引:1  
The microstructure measurements taken during the summer seasons of 2009 and 2010 in the northern South China Sea (between 18°N and 22.5°N, and from the Luzon Strait to the eastern shelf of China) were used to estimate the averaged dissipation rate in the upper pycnocline 〈ε p〉 of the deep basin and on the shelf. Linear correlation between 〈ε p〉 and the estimates of available potential energy of internal waves, which was found for this data set, indicates an impact of energetic internal waves on spatial structure and temporal variability of 〈ε p〉. On the shelf stations, the bottom boundary layer depth-integrated dissipation $ {\widehat{\varepsilon}}_{\mathrm{BBL}} $ reaches 17–19 mW/m2, dominating the dissipation in the water column below the surface layer. In the pycnocline, the integrated dissipation $ {\widehat{\varepsilon}}_{\mathrm{p}} $ was mostly ~10–30 % of $ {\widehat{\varepsilon}}_{\mathrm{BBL}} $ . A weak dependence of bin-averaged dissipation $ \overline{\varepsilon} $ on the Richardson number was noted, according to $ \overline{\varepsilon}={\varepsilon}_0+\frac{\varepsilon_{\mathrm{m}}}{{\left(1+ Ri/R{i}_{\mathrm{cr}}\right)}^{1/2}} $ , where ε 0 + ε m is the background value of $ \overline{\varepsilon} $ for weak stratification and Ri cr?=?0.25, pointing to the combined effects of shear instability of small-scale motions and the influence of larger-scale low frequency internal waves. The latter broadly agrees with the MacKinnon–Gregg scaling for internal-wave-induced turbulence dissipation.  相似文献   

6.
7.
Waves, circulation and vertical dependence   总被引:1,自引:1,他引:0  
Longuet-Higgins and Stewart (J Fluid Mech 13:481–504, 1962; Deep-Sea Res 11:529–562, 1964) and later Phillips (1977) introduced the problem of waves incident on a beach, from deep to shallow water. From the wave energy equation and the vertically integrated continuity equation, they inferred velocities to be Stokes drift plus a return current so that the vertical integral of the combined velocities was nil. As a consequence, it can be shown that velocities of the order of Stokes drift rendered the advective term in the momentum equation negligible resulting in a simple balance between the horizontal gradients of the vertically integrated elevation and wave radiation stress terms; the latter was first derived by Longuet-Higgins and Stewart. Mellor (J Phys Oceanogr 33:1978–1989, 2003a), noting that vertically integrated continuity and momentum equations were not able to deal with three-dimensional numerical or analytical ocean models, derived a vertically dependent theory of wave–circulation interaction. It has since been partially revised and the revisions are reviewed here. The theory is comprised of the conventional, three-dimensional, continuity and momentum equations plus a vertically distributed, wave radiation stress term. When applied to the problem of waves incident on a beach with essentially zero turbulence momentum mixing, velocities are very large and the simple balance between elevation and radiation stress gradients no longer prevails. However, when turbulence mixing is reinstated, the vertically dependent radiation stresses produce vertical velocity gradients which then produce turbulent mixing; as a consequence, velocities are reduced, but are still larger by an order of magnitude compared to Stokes drift. Nevertheless, the velocity reduction is sufficient so that elevation set-down obtained from a balance between elevation gradient and radiation stress gradients is nearly coincident with that obtained by the aforementioned papers. This paper includes four appendices. The first appendix demonstrates the numerical process by which Stokes drift is excluded from the turbulence stress parameterization in the momentum equation. A second appendix determines a bottom slope criterion for the application of linear wave relations to the derivation of the wave radiation stress. The third appendix explores the possibility of generalizing results by non-dimensionalization. The final appendix applies the basic theory to a problem introduced by Bennis and Ardhuin (J Phys Oceanogr 41:2008–2012, 2011).  相似文献   

8.
Rock masses contain ubiquitous multiscale heterogeneities, which (or whose boundaries) serve as the surfaces of discontinuity for some characteristics of the stress state, e.g., for the orientation of principal stress axes. Revealing the regularities that control these discontinuities is a key to understanding the processes taking place at the boundaries of the heterogeneities and for designing the correct procedures for reconstructing and theoretical modeling of tectonic stresses. In the present study, the local laws describing the refraction of the axes of extreme principal stresses T 1 (maximal tension in the deviatoric sense) and T 3 (maximal compression) of the Cauchy stress tensor at the transition over the elementary area n of discontinuity whose orientation is specified by the unit normal n are derived. It is assumed that on the area n of discontinuity, frictional contact takes place. No hypotheses are made on the constitutive equations, and a priori constraints are not posed on the orientation on the stress axes. Two domains, which adjoin area n on the opposite sides and are conventionally marked + and ?, are distinguished. In the case of the two-dimensional (2D) stress state, any principal stress axis on passing from domain ? to domain + remains in the same quadrant of the plane as the continuation of this axis in domain +. The sign and size of the refraction angle depend on the sign and amplitude of the jump of the normal stress, which is tangential to the surface of discontinuity. In the three-dimensional (3D) case, the refraction of axes T 1 and T 3 should be analyzed simultaneously. For each side, + and ?, the projections of the T 1 and T 3 axes on the generally oriented plane n form the shear sectors S + and S ?, which are determined unambiguously and to whose angular domains the possible directions p + and p ? of the shear stress vectors belong. In order for the extreme stress axes T 1 + ,T 3 + and T 1 ? , T 3 ? to be statically compatible on the generally oriented plane n, it is required that sectors S + and S ? had a nonempty intersection. The direction vectors p + and p ? are determined uniquely if, besides axes T 1 ? , T 3 ? and T 1 + , T 3 + , also the ratios of differential stresses R + and R ? (0 ≤ R ± ≤ 1) are known. This is equivalent to specifying the reduced stress tensors T R + and T R ? The necessary condition for tensors T R + and T R ? being statically compatible on plane n is the equality p + = p ?. In this paper, simple methods are suggested for solving the inverse problem of constructing the set of the orientations of the extreme stress axes from the known direction p of the shear stress vector on plane n and from the data on the shear sector. Based on these methods and using the necessary conditions of local equilibrium on plane n formulated above, all the possible orientations of axes T 1 + , T 3 + are determined if the projections of axes T 1 ? , T 3 ? axes on side — are given. The angle between the projections of axes T 1 + , T 1 ? and/or T 3 + , T 3 ? on the plane can attain 90°. Besides the general case, also the particular cases of the contact between the degenerate stress states and the special position of plane n relative to the principal stress axes are thoroughly examined. Generalization of the obtained results makes it possible to plot the local diagram of the orientations of axes T 1 + , T 3 + for a given sector S ?. This diagram is a so-called stress orientation sphere, which is subdivided into three pairs of areas (compression, tension, and compression-extension). The tension and compression zones cannot contain the poles of T 3 + and T 1 + axes, respectively. The compression-extension zones can contain the poles of either T 1 + or T 3 + axis but not both poles simultaneously. In the particular case when the shear stress vector has a unique direction p ? on side ?, the areas of compression-extension disappear and the diagram is reduced to a beach-ball plot, which visualizes the focal mechanism solution of an earthquake. If area n is a generally oriented plane and if the orientation of the pairs of the statically compatible axes T 1 ? , T 3 ? and T 1 + , T 3 + is specified, then, the stress values on side + are uniquely determined from the known stress values on side ?. From the value of differential stress ratio R ?, one can calculate the value of R +, and using the values of the principal stresses on side ?, determine the total stress tensor T + on side +. The obtained results are supported by the laboratory experiments and drilling data. In particular, these results disclose the drawbacks of some established notions and methods in which the possible refraction of the stress axes is unreasonably ignored or taken into account improperly. For example, it is generally misleading to associate the slip on the preexisting fault with the orientation of any particular trihedron of the principal stress axes. The reconstruction should address the potentially statically compatible principal stress axes, which are differently oriented on opposite sides of the fault plane. The fact that, based on the orientation of the intraplate principal stresses at the base of the lithosphere, one cannot make a conclusion on the active or passive influence of the mantle flows on the lithospheric plate motion is another example. The present relationships linking the stress values on the opposite sides of the fault plane on which the orientations of the principal stress axes are known demonstrate the incorrectness of the existing methods, in which the reduced stress tensors within the material domains are reconstructed without allowance for the dynamic interaction of these domains with their neighbors. In addition, using the obtained results, one can generalize the notion of the zone of dynamical control of a fault onto the case of the existence of discontinuities in this region and analyze the stress transfer across the system of the faults.  相似文献   

9.
To provide coastal engineers and scientists with a quantitative evaluation of nearshore numerical wave models in reef environments, we review and compare three commonly used models with detailed laboratory observations. These models are the following: (1) SWASH (Simulating WAves till SHore) (Zijlema et al. 2011), a phase-resolving nonlinear shallow-water wave model with added nonhydrostatic terms; (2) SWAN (Simulating WAve Nearshore) (Booij et al. 1999), a phase-averaged spectral wave model; and (3) XBeach (Roelvink et al. 2009), a coupled phase-averaged spectral wave model (applied to modeling sea-swell waves) and a nonlinear shallow-water model (applied to modeling infragravity waves). A quantitative assessment was made of each model’s ability to predict sea-swell (SS) wave height, infragravity (IG) wave height, wave spectra, and wave setup ( \( \overline{\eta} \) ) at five locations across the laboratory fringing reef profile of Demirbilek et al. (2007). Simulations were performed with the “recommended” empirical coefficients as documented for each model, and then the key wave-breaking parameter for each model (α in SWASH and γ in both SWAN and XBeach) was optimized to most accurately reproduce the observations. SWASH, SWAN, and XBeach were found to be capable of predicting SS wave height variations across the steep fringing reef profile with reasonable accuracy using the default coefficients. Nevertheless, tuning of the key wave-breaking parameter improved the accuracy of each model’s predictions. SWASH and XBeach were also able to predict IG wave height and spectral transformation. Although SWAN was capable of modeling the SS wave height, in its current form, it was not capable of modeling the spectral transformation into lower frequencies, as evident in the underprediction of the low-frequency waves.  相似文献   

10.
A global dataset of more than 3,000 ground motion records from 536 sites from Greece, Italy, Turkey, USA and Japan is compiled and used to propose code-oriented elastic acceleration response spectra and soil amplification factors for a new site classification system, which, besides the classical geotechnical parameters $N_{SPT}, S_{u}$ and PI, uses also the fundamental period of the site, the thickness of soil deposits and the average shear wave velocity to the seismic bedrock, instead of $V_{s,30}$ . We propose a new classification system with the associated amplification factors and normalized response spectra for two seismicity levels, i.e. $M_{s}\le 5.5$ and $M_{s}>5.5$ . Uncertainties in the estimation of soil amplification factors are captured using a logic-tree approach, which allows the efficient use of alternative models and methods. The aim of this work is to improve the present EC8 soil classification. The effectiveness of the proposed classification system is compared to that of EC8 classification system using an error term, which represents the average dispersion of data within all categories of a given classification scheme. Error terms for the new classification system are lower than the error terms for EC8 classification system at all periods.  相似文献   

11.
Unloaded natural rock masses are known to generate seismic signals (Green et al., 2006; Hainzl et al., 2006; Husen et al., 2007; Kraft et al., 2006). Following a 1,000 m3 mass failure into the Mediterranean Sea, centimeter-wide tensile cracks were observed to have developed on top of an unstable segment of the coastal cliff. Nanoseismic monitoring techniques (Wust-Bloch and Joswig, 2006; Joswig, 2008), which function as a seismic microscope for extremely weak seismic events, were applied to verify whether brittle failure is still generated within this unconsolidated sandstone mass and to determine whether it can be detected. Sixteen days after the initial mass failure, three small-aperture sparse arrays (Seismic Navigation Systems-SNS) were deployed on top of this 40-m high shoreline cliff. This paper analyzes dozens of spiky nanoseismic (?2.2 ≥ M L ≥ ?3.4) signals recorded over one night in continuous mode (at 200 Hz) at very short slant distances (3–67 m). Waveform characterization by sonogram analysis (Joswig, 2008) shows that these spiky signals are all short in duration (>0.5 s). Most of their signal energy is concentrated in the 10–75 Hz frequency range and the waveforms display high signal similarity. The detection threshold of the data set reaches M L ?3.4 at 15 m and M L ?2.7 at 67 m. The spatial distribution of source signals shows 3-D clustering within 10 m from the cliff edge. The time distribution of M L magnitude does not display any decay pattern of M L over time. This corroborates an unusual event decay over time (modified Omori’s law), whereby an initial quiet period is followed by regained activity, which then fades again. The polarization of maximal waveform amplitude was used to estimate spatial stress distribution. The orientation of ellipses displaying maximal signal energy is consistent with that of tensile cracks observed in the field and agrees with rock mechanics predictions. The M L– surface rupture length relationship displayed by our data fits a constant-slope extrapolation of empirical data collected by Wells and Coppersmith (1994) for normal fault features at much larger scale. Signal characterization and location as well as the absence of direct anthropogenic noise sources near the monitoring site, all indicate that these nanoseismic signals are generated by brittle failure within the top section of the cliff. The atypical event decay over time that was observed suggests that the cliff material is undergoing post-collapse bulk strain accommodation. This feasibility study demonstrates the potential of nanoseismic monitoring in rapidly detecting, locating and analyzing brittle failure generated within unconsolidated material before total collapse occurs.  相似文献   

12.
Seismic coda wave attenuation ( $ Q_{\text{c}}^{ - 1} $ ) characteristics in the Garhwal region, northwestern Himalaya is studied using 113 short-period, vertical component seismic observations from local events with hypocentral distance less than 250?km and magnitude range between 1.0 to 4.0. They are located mainly in the vicinity of the Main Boundary Thrust (MBT) and the Main Central Thrust (MCT), which are well-defined tectonic discontinuities in the Himalayas. Coda wave attenuation ( $ Q_{\text{c}}^{ - 1} $ ) is estimated using the single isotropic scattering method at central frequencies 1.5, 3, 5, 7, 9, 12, 16, 20, 24 and 28?Hz using several starting lapse times and coda window lengths for the analysis. Results show that the ( $ Q_{\text{c}}^{ - 1} $ ) values are frequency dependent in the considered frequency range, and they fit the frequency power law ( $ Q_{\text{c}}^{ - 1} \left( f \right) = Q_{0}^{ - 1} f^{ - n} $ ). The Q 0 (Q c at 1?Hz) estimates vary from about 50 for a 10?s lapse time and 10?s window length, to about 350 for a 60?s lapse time and 60?s window length combination. The exponent of the frequency dependence law, n ranges from 1.2 to 0.7; however, it is greater than 0.8, in general, which correlates well with the values obtained in other seismically and tectonically active and highly heterogeneous regions. The attenuation in the Garhwal region is found to be lower than the Q c ?1 values obtained for other seismically active regions of the world; however, it is comparable to other regions of India. The spatial variation of coda attenuation indicates that the level of heterogeneity decreases with increasing depth. The variation of coda attenuation has been estimated for different lapse time and window length combinations to observe the effect with depth and it indicates that the upper lithosphere is more active seismically as compared to the lower lithosphere and the heterogeneity decreases with increasing depth.  相似文献   

13.
Recent acoustic Doppler current profiler (ADCP)-measurements in the Scheldt estuary near Antwerp, Belgium, revealed anomalous, i.e. anti-clockwise circulations in a left bend during the major part of the flood period; these circulations were established shortly after the turn of the tide. During ebb, anti-clockwise circulations persisted, as predicted by classical theory. These data were analysed with a 3D and a 1DV-model. The 3D simulations reveal that the anomalous circulations are found when salinity is included in the computations—without salinity “normal” circulations were found. From analytical and 1DV simulations, it is concluded that a longitudinal salinity gradient ${\partial S} \mathord{\left/ {\vphantom {{\partial S} {\partial x}}} \right. \kern-0em} {\partial x}$ may induce a near-bed maximum in flow velocity reversing the direction of the secondary currents. The 1DV-model was then used to assess the contribution of various processes one by one. It was found that because of a reduction in vertical mixing, the vertical velocity profile is not at equilibrium during the first phase of accelerating tide, further enhancing the effects of ${\partial S} \mathord{\left/ {\vphantom {{\partial S} {\partial x}}} \right. \kern-0em} {\partial x}$ . A small vertical salinity gradient ${\partial S} \mathord{\left/ {\vphantom {{\partial S} {\partial z}}} \right. \kern-0em} {\partial z}$ appeared to have a very large effect as the crosscurrents of the secondary circulations induced by ${\partial S} \mathord{\left/ {\vphantom {{\partial S} {\partial x}}} \right. \kern-0em} {\partial x}$ became an order of larger magnitude. However, at the site under consideration, the effects of transverse salinity gradients, generated by differential advection in the river bend, were dominant: adverse directions of the secondary circulations were found even when the vertical velocity profile became more regular with a more or less logarithmic shape, i.e. when the effects of ${\partial S} \mathord{\left/ {\vphantom {{\partial S} {\partial x}}} \right. \kern-0em} {\partial x}$ and ${\partial S} \mathord{\left/ {\vphantom {{\partial S} {\partial z}}} \right. \kern-0em} {\partial z}$ did not play a dominant role anymore. It is argued that data on the secondary velocity structure, which can be measured easily owing to today’s developments in ADCP equipment, may serve as an indicator for the accuracy at which the salinity field is computed with 3D numerical models. Moreover, the large effect of the salinity structure on the velocity field must have a large impact on the morphological development of estuaries, which should therefore be accounted for in morphological modelling studies.  相似文献   

14.
— Surface-wave amplitudes from explosion sources show less variation for a given event han body wave amplitudes, so it is natural to expect that yield estimates derived from surface waves will be more accurate than yield estimates derived from body waves. However, yield estimation from surface waves is complicated by the presence of tectonic strain release, which acts like one or more earthquake sources superimposed on top of the explosion. Moment-tensor inversion can be used to remove the tectonic component of the surface waves, however moment-tensor inversion for shallow sources is inherently non-unique so the explosion isotropic moment cannot be determined with the necessary accuracy by this means. Explosions on an island or near a mountain slope can exhibit anomalous surface waves similar to those caused by tectonic strain release. These complications cause yield estimates derived from surface waves to be less accurate than yield estimates from body waves recorded on a well-calibrated network with good coverage. Surface-wave amplitudes can be expressed as a surface-wave magnitude M s , which is defined as the logarithm of the amplitude plus a distance correction, or as a path corrected spectral magnitude, log $M^{\prime}_0$ , which is derived from the surface-wave spectrum. We derive relations for M s vs. yield and log $M^{\prime}_0$ vs. yield for a large data set and estimate the accuracy of these estimates.  相似文献   

15.
A modified formula of the cumulative frequency-magnitude relation has been formulated and tested in a previous paper by the authors of this study. Based on the modified relationship, the following reoccurrence formulas have been obtained.
  1. For the ‘T-years period’ larger earthquake magnitude,M T $$M_T = \frac{1}{{A_3 }}ln\frac{{A_2 }}{{(1/T) + A_1 }}.$$
  2. For the value of the maximum earthquake magnitude, which is exceeded with probabilityP inT-years period,M PT $$M_{PT} = \frac{{ln(A_2 .T)}}{{A_3 }} - \frac{{ln[A_1 .T - ln(1 - P)]}}{{A_3 }}.$$
  3. For the probability of occurrence of an earthquake of magnitudeM in aT-years period,P MT $$P_{MT} = 1 - \exp [ - T[ - A_1 + A_2 \exp ( - A_3 M)]].$$
The above formulas provide estimates of the probability of reoccurrence of the largest earthquake events which are significantly more realistic than those based on the Gutenberg-Richter relationships; at least for numerous tested earthquake samples from the major area of Greece.  相似文献   

16.
We determined depth variation of the 410- and 660-km discontinuities beneath southeastern China by common-converted-point stacking of \(\rm P\) -wave receiver functions of 121 permanent Chinese seismic stations. We then combined the results with seismic velocity variation to estimate temperature and water content variations in the mantle transition zone of the region. Previous tomographic studies have shown a stagnant slab in the mantle transition zone in eastern Asia that is connected to subduction of the western Pacific. Temperature variations obtained clearly outline the shape of the stagnant slab, with its western edge at 113.5 \(^\circ\) E and the southern edge at 28.5 \(^\circ\) N. The correlation between the location of the stagnant slab and surface tectonics suggests that the Cenozoic extension in eastern China is closely associated with the subduction of the western Pacific and its eastward migration. The water content of the stagnant slab is lower than in surrounding slabs, suggesting that the water has already been released from the subducting slab into the upper mantle.  相似文献   

17.
Three-dimensional frequency dependent S-wave quality factor (Qβ(f)) value for the central Honshu region of Japan has been determined in this paper using an algorithm based on inversion of strong motion data. The method of inversion for determination of three-dimensional attenuation coefficients is proposed by Hashida and Shimazaki (J Phys Earth. 32, 299–316, 1984) and has been used and modified by Joshi (Curr Sci. 90, 581–585, 2006; Nat Hazards. 43, 129–146, 2007) and Joshi et al. (J. Seismol. 14, 247–272, 2010). Twenty-one earthquakes digitally recorded on strong motion stations of Kik-net network have been used in this work. The magnitude of these earthquake ranges from 3.1 to 4.2 and depth ranging from 5 to 20 km, respectively. The borehole data having high signal to noise ratio and minimum site effect is used in the present work. The attenuation structure is determined by dividing the entire area into twenty-five three-dimensional blocks of uniform thickness having different frequency-dependent shear wave quality factor. Shear wave quality factor values have been determined at frequencies of 2.5, 7.0 and 10 Hz from record in a rectangular grid defined by 35.4°N to 36.4°N and 137.2°E to 138.2°E. The obtained attenuation structure is compared with the available geological features in the region and comparison shows that the obtained structure is capable of resolving important tectonic features present in the area. The proposed attenuation structure is compared with the probabilistic seismic hazard map of the region and shows that it bears some remarkable similarity in the patterns seen in seismic hazard map.  相似文献   

18.
Abstract

Temperature profiles (temperature as function of depth) can be used to derive vertical flow velocities or recharge rates. In many cases, solutions to the one-dimensional (1-D) heat transport equation are used, considering steady-state boundary conditions. Factors which can influence the derivation of the mean vertical flow velocity are studied here. Therefore, an explicit finite-difference approximation to the 1-D heat transport equation coupled with an inverse scheme was used to interpret temperature profiles. Measurement error (larger than 0.05°C) can result in important deviation of the derived mean flow velocity. Variation of vertical flow velocity as a function of time leads to asymmetric temperature envelopes. Yearly variation in vertical flow velocities, or temperature variations of the recharge water, also results in asymmetric temperature envelopes. Interpretation of these asymmetric envelopes leads to important differences between derived and actual mean vertical flow velocities.

Citation Vandenbohede, A. &; Lebbe, L. (2010) Vandenbohede, A. and Lebbe, L. 2010. Parameter estimation based on vertical heat transport in the surficial zone. Hydrogeol. J., 18 , 931–943,doi:10.1007/s10040-009–0557–5 [Google Scholar] Recharge assessment by means of vertical temperature profiles: analysis of possible influences. Hydrol. Sci. J. 55(5), 792–804.  相似文献   

19.
Towards fully data driven ground-motion prediction models for Europe   总被引:2,自引:2,他引:0  
We have used the Artificial Neural Network method (ANN) for the derivation of physically sound, easy-to-handle, predictive ground-motion models from a subset of the Reference database for Seismic ground-motion prediction in Europe (RESORCE). Only shallow earthquakes (depth smaller than 25 km) and recordings corresponding to stations with measured $V_{s30}$ properties have been selected. Five input parameters were selected: the moment magnitude $M_{W}$ , the Joyner–Boore distance $R_{JB}$ , the focal mechanism, the hypocentral depth, and the site proxy $V_{S30}$ . A feed-forward ANN type is used, with one 5-neuron hidden layer, and an output layer grouping all the considered ground motion parameters, i.e., peak ground acceleration (PGA), peak ground velocity (PGV) and 5 %-damped pseudo-spectral acceleration (PSA) at 62 periods from 0.01 to 4 s. A procedure similar to the random-effects approach was developed to provide between and within event standard deviations. The total standard deviation ( $\sigma $ ) varies between 0.298 and 0.378 (log $_{10}$ unit) depending on the period, with between-event and within-event variabilities in the range 0.149–0.190 and 0.258–0.327, respectively. Those values prove comparable to those of conventional GMPEs. Despite the absence of any a priori assumption on the functional dependence, our results exhibit a number of physically sound features: magnitude scaling of the distance dependency, near-fault saturation distance increasing with magnitude, amplification on soft soils and even indications for nonlinear effects in softer soils.  相似文献   

20.
The seismic behaviour of caisson foundations supporting typical bridge piers is analysed with 3D finite elements, with due consideration to soil and interface nonlinearities. Single-degree-of freedom oscillators of varying mass and height, simulating heavily and lightly loaded bridge piers, founded on similar caissons are studied. Four different combinations of the static ( $\text{ FS }_\mathrm{V}$ FS V ) and seismic ( $\text{ FS }_\mathrm{E}$ FS E ) factors of safety are examined: (1) a lightly loaded ( $\text{ FS }_\mathrm{V}= 5$ FS V = 5 ) seismically under-designed ( $\text{ FS }_\mathrm{E} < 1$ FS E < 1 ) caisson, (2) a lightly loaded seismically over-designed ( $\text{ FS }_\mathrm{E} >1$ FS E > 1 ) caisson, (3) a heavily loaded ( $\text{ FS }_\mathrm{V} = 2.5$ FS V = 2.5 ) seismically under-designed ( $\text{ FS }_\mathrm{E} < 1$ FS E < 1 ) caisson and (4) a heavily loaded seismically over-designed caisson. The analysis is performed with use of seismic records appropriately modified so that the effective response periods (due to soil-structure-interaction effects) of the studied systems correspond to the same spectral acceleration, thus allowing their inelastic seismic performance to be compared on a fair basis. Key performance measures of the systems are then contrasted, such as: accelerations, displacements, rotations and settlements. It is shown that the performance of the lightly loaded seismically under-designed caisson is advantageous: not only does it reduce significantly the seismic load to the superstructure, but it also produces minimal residual displacements of the foundation. For heavily loaded foundations, however ( $\text{ FS }_{V} = 2.5$ FS V = 2.5 ), the performance of the two systems (over and under designed) is similar.  相似文献   

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

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