首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Recently Nagata et al. (J Geophys Res 117:B02314, 2012) have proposed a new version of rate- and state-dependent friction law (RSF) that seems to have eventually resolved all the previously known discrepancies in the existing RSFs from laboratory observations. The values of a and b, empirical RSF parameters determined by fitting the same laboratory experiments, have been revised to be five times greater and a newly noticed weakening effect by shear stress with a coefficient c has been introduced. By using this revised RSF, we reinvestigated a problem of 2D quasi-static nucleation on faults. A crack-like nucleation-zone expansion known for the ‘aging’ version of RSF is not sustainable with the ‘Nagata’ law, which is understandable as the Nagata law does not produce a slip-weakening distance proportional to the involved strength reduction, an aging law’s feature that contradicts laboratory observations. The later stage of Nagata-law nucleation shows localization of quasi-static slip within a limited spatial extent, but the localization is much milder than that predicted by the ‘slip’ version of RSF. With an appropriate c parameter of the Nagata law, the nucleation size seems to be reduced only by a factor from that of the aging law.  相似文献   

2.
We estimate seismological fracture energies from two subsets of events selected from the seismic sequences of L’Aquila (2009), and Northridge (1994): 57 and 16 selected events, respectively, including the main shocks. Following Abercrombie and Rice (Geophys J Int 162: 406–424, 2005), we postulate that fracture energy (G) represents the post-failure integral of the dynamic weakening curve, which is described by the evolution of shear traction as a function of slip. Following a direct-wave approach, we compute mainshock-/aftershock-source spectral ratios, and analyze them using the approach proposed by Malagnini et al. (Pure Appl. Geophys., this issue, 2014) to infer corner frequencies and seismic moment. Our estimates of source parameters (including fracture energies) are based on best-fit grid-searches performed over empirical source spectral ratios. We quantify the source scaling of spectra from small and large earthquakes by using the MDAC formulation of Walter and Taylor (A revised Magnitude and Distance Amplitude Correction (MDAC2) procedure for regional seismic discriminants, 2001). The source parameters presented in this paper must be considered as point-source estimates representing averages calculated over specific ruptured portions of the fault area. In order to constrain the scaling of fracture energy with coseismic slip, we investigate two different slip-weakening functions to model the shear traction as a function of slip: (i) a power law, as suggested by Abercrombie and Rice (Geophys J Int 162: 406–424, 2005), and (ii) an exponential decay. Our results show that the exponential decay of stress on the fault allows a good fit between measured and predicted fracture energies, both for the main events and for their aftershocks, regardless of the significant differences in the energy budgets between the large (main) and small earthquakes (aftershocks). Using the power-law slip-weakening function would lead us to a very different situation: in our two investigated sequences, if the aftershock scaling is extrapolated to events with large slips, a power law (a la Abercrombie and Rice) would predict unrealistically large stress drops for large, main earthquakes. We conclude that the exponential stress evolution law has the advantage of avoiding unrealistic stress drops and unbounded fracture energies at large slip values, while still describing the abrupt shear-stress degradation observed in high-velocity laboratory experiments (e.g., Di Toro et al., Fault lubrication during earthquakes, Nature 2011).  相似文献   

3.
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.  相似文献   

4.
The earthquake occurred in Lushan County on 20 April, 2013 caused heavy casualty and economic loss. In order to understand how the seismic energy propagates during this earthquake and how it causes the seismic hazard, we simulated the strong ground motions from a representative kinematic source model by Zhang et al. (Chin J Geophys 56(4):1408–1411, 2013) for this earthquake. To include the topographic effects, we used the curved grids finite difference method by Zhang and Chen (Geophys J Int 167(1):337–353, 2006), Zhang et al. (Geophys J Int 190(1):358–378, 2012) to implement the simulations. Our results indicated that the majority of seismic energy concentrated in the epicentral area and the vicinal Sichuan Basin, causing the XI and VII degree intensity. Due to the strong topographic effects of the mountain, the seismic intensity in the border area across the northeastern of Boxing County to the Lushan County also reached IX degree. Moreover, the strong influence of topography caused the amplifications of ground shaking at the mountain ridge, which is easy to cause landslides. These results are quite similar to those observed in the Wenchuan earthquake of 2008 occurred also in a strong topographic mountain area.  相似文献   

5.
It is well known that quantitative estimation of slip distributions on fault plane is one of the most important issues for earthquake source inversion related to the fault rupture process. The characteristics of slip distribution on the main fault play a fundamental role to control strong ground motion pattern. A large amount of works have also suggested that variable slip models inverted from longer period ground motion recordings are relevant for the prediction of higher frequency ground motions. Zhang et al. (Chin J Geophys 56:1412–1417, 2013) and Wang et al. (Chin J Geophys 56:1408–1411, 2013) published their source inversions for the fault rupturing process soon after the April 20, 2013 Lushan earthquake in Sichuan, China. In this study, first, we synthesize two forward source slip models: the value of maximum slip, fault dimension, size, and dimension of major asperities, and corner wave number obtained from Wang’s model is adopted to constrain the generation of k ?2 model and crack model. Next, both inverted and synthetic slip models are employed to simulate the ground motions for the Lushan earthquake based on the stochastic finite-fault method. In addition, for a comparison purpose, a stochastic slip model and another k ?2 model (k ?2 model II) with 2 times value of corner wave number of the original k ?2 model (k ?2 model I) are also employed for simulation for Lushan event. The simulated results characterized by Modified Mercalli Intensity (MMI) show that the source slip models based on the inverted and synthetic slip distributions could capture many basic features associated with the ground motion patterns. Moreover, the simulated MMI distributions reflect the rupture directivity effect and the influence of the shallow velocity structure well. On the other hand, the simulated MMI by stochastic slip model and k ?2 model II is apparently higher than observed intensity. By contrast, our simulation results show that the higher frequency ground motion is sensitive to the degree of slip roughness; therefore, we suggest that, for realistic ground‐motion simulations due to future earthquake, it is imperative to properly estimate the slip roughness distribution.  相似文献   

6.
We performed three-dimensional (3-D) finite difference simulations of long-period ground motions (2–10 s) in the Kanto basin using the Japan Seismic Hazard Information Station (J-SHIS 2009), Yamada and Yamanaka (Exploration Geophysics 65(3):139–150, 2012) (YY), and Head Quarter for Earthquake Research Promotion (HERP 2012) velocity models for two intermediate depth (68–80 km) moderate earthquakes (Mw 5.8–5.9), which occurred beneath the Kanto basin. The models primarily differ in the basic data set used in the construction of the velocity models. The J-SHIS and HERP models are the results of integration of mainly geological, geophysical, and earthquake data. On the other hand, the YY model is oriented towards the microtremor-array-observation data. We obtained a goodness of fit between the observed and synthetic data based on three parameters, peak ground velocities (PGVs), smoothed Fourier spectra (FFT), and cross-correlations, using an algorithm proposed by Olsen and Mayhew (Seism Res Lett 81:715–723, 2010). We found that the three models reproduced the PGVs and FFT satisfactorily at most sites. However, the models performed poorly in terms of cross-correlations especially at the basin edges. We found that the synthetics using the YY model overestimate the observed waveforms at several sites located in the areas having V s 0.3 km/s in the top layer; on the other hand, the J-SHIS and HERP models explain the waveforms better at the sites and perform similarly at most sites. We also found that the J-SHIS and HERP models consist of thick sediments beneath some sites, where the YY model is preferable. Thus, we have concluded that the models require revisions for the reliable prediction of long-period ground motions from future large earthquakes.  相似文献   

7.
The 23 October 2011 Van (Mw 7.1) earthquake that occurred in Eastern Turkey resulted in heavy damage particularly in the city of Van and town of Ercis. This paper presents ground motion simulations of Van earthquake by using stochastic finite fault method (EXSIM, Motazedian and Atkinson in Bull Seismol Soc Am 95:995–1010, 2005; Boore in Bull Seismol Soc Am 99:3202–3216, 2009) that provides a simple and effective tool to generate high frequency strong motion. The input parameters related to source, path, and site effects are calibrated on the basis of minimizing the error functions between simulations and observations both in time and frequency domain. Validated model parameters are used to produce synthetics in regional extent with the aim of understanding the level and distribution of the ground shaking particularly in the near fault region where no recordings are available within the 40 km of the epicenter. This paper evaluates the effect of two different slip models on ground motion intensity measures over the area of interest and addresses the variability in the near fault region associated with the source effect. The synthetics are compared with the corresponding estimations of ground motion prediction equations by Boore and Atkinson (Earthq Spectra 24:99–138, 2008), Akkar and Bommer (Seismol Res Lett 81:195–206, 2010) and Akkar and Cagnan (Bull Seismol Soc Am 100:2978–2995, 2010). Our results indicate that despite the limitation of the method for incorporating the directivity effect and inadequate representation of the soil conditions at the individual stations, a satisfactory match between synthetics and observations are obtained both in time and frequency domain. Spatial distributions of the synthetics in regional level also show reasonable correlation with ground motion prediction equations and damage observations.  相似文献   

8.
It has been two decades since the last comprehensive standard model of ambient earth noise was published Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The PETERSON model was updated by analyzing the absolute quietest conditions for stations within the GSN (Berger et al. in J Geophys Res 109, 2005; Mcnamara and Buland in Bull Seism Soc Am 94:1517–1527, 2004; Ringler et al. in Seismol Res Lett 81(4) doi:10.1785/gssrl.81.4.605, 2010). Unfortunately, both the original model and the updated models did not include any deployed station in North Africa and Middle East, which reflects the noise levels within the desert environment of those regions. In this study, a survey was conducted to create a new seismic noise model from very broadband stations which recently deployed in North Africa. For this purpose, 1 year of continuous recording of seismic noise data of the Egyptian National Seismic Network (ENSN) was analyzed in order to create a new noise model. Seasonal and diurnal variations in noise spectra were recorded at each station. Moreover, we constructed a new noise model for each individual station. Finally, we obtained a new cumulative noise model for all the stations. We compared the new high-noise model (EHNM) and new low-noise model (ELNM) with both the high-noise model (NHNM) and low-noise model (NLNM) of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The obtained noise levels are considerably lower than low-noise model of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993) at ultra long period band (ULP band), but they are still below the high-noise model of Peterson (Observations and modelling of seismic background noise, US Geological Survey, open-file report 93–322, 1993). The results of this study could be considered as a first step to create permanent seismic noise models for North Africa and Middle East regions.  相似文献   

9.
Historically, surface subsidence as a result of subsurface groundwater fluctuations have produced important and, at times, catastrophic effects, whether natural or anthropogenic. Over the past 30?years, numerical and analytical techniques for the modeling of this surface deformation, based upon elastic and poroelastic theory, have been remarkably successful in predicting the magnitude of that deformation (Le Mouélic and Adragna in Geophys Res Lett 29:1853, 2002). In this work we have extended the formula for a circular-shaped aquifer (Geertsma in J Petroleum Tech 25:734–744, 1973) to a more realistic elliptical shape. We have improved the accuracy of the approximation by making use of the cross terms of the expansion for the elliptic coordinates in terms of the eccentricity, e, and the mean anomaly angle, M, widely used in astronomy. Results of a number of simulations, in terms of e and M developed from the transcendental Kepler equation, are encouraging, giving realistic values for the elliptical approximation of the vertical deformation due to groundwater change. Finally, we have applied the algorithm to modeling of groundwater in southern California.  相似文献   

10.
We analyze the impact of a linear trend in the mean log-conductivity on the transport of a conservative tracer in a bounded domain. The effects of such a linear trend on solute transport were analyzed in depth for unbounded domains (Rubin and Seong, Water Resour Res 30(11):2901–2911, 1994; Indelman and Rubin, Water Resour Res 31(5):1257–1265, 1995; Water Resour Res 32(5):1257–1265, 1996), whereas studies concerning this special case of medium nonstationarity in finite domains usually focus on head or flow statistics (Guadagnini et al., Stoch Environ Res Risk Assess, 17:394–407, 2003). In this study both ensemble and effective plume moments are provided for an instantaneous release of a solute through a linear source normal to the mean flow direction, by taking into account different sizes of the source. The analysis involving a steady velocity field spatially nonstationary is developed by using the stochastic finite element method. Results show that ensemble moments are affected by increasing trends both parallel and normal to the mean flow direction, but the impact on effective plume moments is very different. A parallel trend does not seem to influence the effective second moments; while a normal trend, although modifies the transverse effective moment only weakly, strongly increases the longitudinal one, especially for large initial sizes of the source. Furthermore, the increase of the particle displacement variance produced by a parallel trend in the finite domain disagrees with the results obtained in an unbounded domain, due to the boundary conditions here considered making both head and velocity moments nonstationary and nonsymmetric.  相似文献   

11.
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).  相似文献   

12.
“Repeating earthquakes” identified by waveform cross-correlation, with inter-event separation of no more than 1 km, can be used for assessment of location precision. Assuming that the network-measured apparent inter-epicenter distance X of the “repeating doublets” indicates the location precision, we estimated the regionalized location quality of the China National Seismograph Network by comparing the “repeating events” in and around China by Schaff and Richards (Science 303: 1176–1178, 2004; J Geophys Res 116: B03309, 2011) and the monthly catalogue of the China Earthquake Networks Center. The comparison shows that the average X value of the China National Seismograph Network is approximately 10 km. The mis-location is larger for the Tibetan Plateau, west and north of Xinjiang, and east of Inner Mongolia, as indicated by larger X values. Mis-location is correlated with the completeness magnitude of the earthquake catalogue. Using the data from the Beijing Capital Circle Region, the dependence of the mis-location on the distribution of seismic stations can be further confirmed.  相似文献   

13.
On the selection of GMPEs for Vrancea subcrustal seismic source   总被引：2，自引：0，他引：2  
The Vrancea subcrustal seismic source is characterized by large magnitude ( $M_{W} \ge 7$ ) intermediate-depth earthquakes that occur two or three times during a century on average. In this study several procedures are used to grade four candidate ground motion prediction equations proposed for Vrancea source in the SHARE project. In the work of Delavaud et al. (J Seismol 16(3):451–473, 2012) four ground motion prediction models developed for subduction zones (Zhao et al. in Bull Seism Soc Am 96(3):898–913, 2006; Atkinson and Boore in Bull Seism Soc Am 93(4):1703–1729, 2003; Youngs et al. in Seism Res Lett 68(1):58–73, 1997; Lin and Lee in Bull Seism Soc Am 98(1):220–240, 2008) are suggested as suitable for Vrancea subcrustal seismic source. The paper presents the appropriateness analysis of the four suggested ground motion prediction equations done using a dataset of 109 triaxial accelerograms recorded during seven Vrancea seismic events with moment magnitude $M_{W}$ between 5.4 and 7.4, occurred in the past 35 years. The strong ground motions were recorded in Romania, as well as in Bulgaria, Republic of Moldova and Serbia. Based on the ground motion dataset several goodness-of-fit measures are used in order to quantify how well the selected models match with the recorded data. The compatibility of the four ground motion prediction models with respect to magnitude scaling and distance scaling implied by strong ground motion dataset is investigated as well. The analyses show that the Youngs et al. (Seism Res Lett 68(1):58–73, 1997) and Zhao et al. (Bull Seism Soc Am 96(3):898–913, 2006) ground motion prediction models have a better fit with the data and can be candidate models for Probabilistic Seismic Hazard Assessment.  相似文献   

14.
Following the 2007, August 15th, M w 8.0, Pisco earthquake in central Peru, Sladen et al. (J Geophys Res 115: B02405, 2010) have derived several slip models of this event. They inverted teleseismic data together with geodetic (InSAR) measurements to look for the co-seismic slip distribution on the fault plane, considering those data sets separately or jointly. But how close to the real slip distribution are those inverted slip models? To answer this crucial question, the authors generated some tsunami records based on their slip models and compared them to DART buoys, tsunami records, and available runup data. Such an approach requires a robust and accurate tsunami model (non-linear, dispersive, accurate bathymetry and topography, etc.) otherwise the differences between the data and the model may be attributed to the slip models themselves, though they arise from an incomplete tsunami simulation. The accuracy of a numerical tsunami simulation strongly depends, among others, on two important constraints: (i) A fine computational grid (and thus the bathymetry and topography data sets used) which is not always available, unfortunately, and (ii) a realistic tsunami propagation model including dispersion. Here, we extend Sladen’s work using newly available data, namely a tide gauge record at Callao (Lima harbor) and the Chilean DART buoy record, while considering a complete set of runup data along with a more realistic tsunami numerical that accounts for dispersion, and also considering a fine-resolution computational grid, which is essential. Through these accurate numerical simulations we infer that the InSAR-based model is in better agreement with the tsunami data, studying the case of the Pisco earthquake indicating that geodetic data seems essential to recover the final co-seismic slip distribution on the rupture plane. Slip models based on teleseismic data are unable to describe the observed tsunami, suggesting that a significant amount of co-seismic slip may have been aseismic. Finally, we compute the runup distribution along the central part of the Peruvian coast to better understand the wave amplification/attenuation processes of the tsunami generated by the Pisco earthquake.  相似文献   

15.
The 2008 Wenchuan earthquake, a major intraplate earthquake with M w 7.9, occurred on the slowly deforming Longmenshan fault. To better understand the causes of this devastating earthquake, we need knowledge of the regional stress field and the underlying geodynamic processes. Here, we determine focal mechanism solutions (FMSs) of the 2008 Wenchuan earthquake sequence (WES) using both P-wave first-motion polarity data and SH/P amplitude ratio (AR) data. As P-wave polarities are more reliable information, they are given priority over SH/P AR, the latter of which are used only when the former has loose constraint on the FMSs. We collect data from three categories: (1) permanent stations deployed by the China Earthquake Administration (CEA); (2) the Western Sichuan Passive Seismic Array (WSPSA) deployed by Institute of Geology, CEA; (3) global stations from Incorporated Research Institutions for Seismology. Finally, 129 events with magnitude over M s 4.0 in the 2008 WES are identified to have well-constrained FMSs. Among them, 83 are well constrained by P-wave polarities only as shown by Cai et al. (Earthq Sci 24(1):115–125, 2011), and the rest of which are newly constrained by incorporating SH/P AR. Based on the spatial distribution and FMSs of the WES, we draw following conclusions: (1) the principle compressional directions of most FMSs of the WES are subhorizontal, generally in agreement with the conclusion given by Cai et al. (2011) but with a few modifications that the compressional directions are WNW–ESE around Wenchuan and ENE–WSW around Qingchuan, respectively. The subhorizontal compressional direction along the Longmenshan fault from SW to NE seems to have a left-lateral rotation, which agrees well with regional stress field inverted by former researchers (e.g., Xu et al., Acta Seismol Sin 30(5), 1987; Acta Geophys Sin 32(6), 1989; Cui et al., Seismol Geol 27(2):234–242, 2005); (2) the FMSs of the events not only reflected the regional stress state of the Longmenshan region, but also were obviously controlled by the faults to some extent, which was pointed out by Cai et al. (2011) and Yi et al. (Chin J Geophys 55(4):1213–1227, 2012); (3) while the 2008 Wenchuan earthquake and some of its strong aftershocks released most of the elastic energy accumulated on the Longmenshan fault, some other aftershocks seem to occur just for releasing the elastic energy promptly created by the 2008 Wenchuan earthquake and some of its strong aftershocks. (4) Our results further suggest that the Longmenshan fault from Wenchuan to Beichuan was nearly fully destroyed by the 2008 Wenchuan earthquake and accordingly propose that there is less probability for great earthquakes in the middle part of the Longmenshan fault in the near future, although there might be a barrier to the southwest of Wenchuan and it is needed to pay some attention on it in the near future.  相似文献   

16.
We apply the Bakun and Wentworth (Bull Seism Soc Am 87:1502–1521, 1997) method to determine the location and magnitude of earthquakes occurred in Central Asia using MSK-64 intensity assignments. The attenuation model previously derived and validated by Bindi et al. (Geophys J Int, 2013) is used to analyse 21 earthquakes that occurred over the period 1885–1964, and the estimated locations and magnitudes are compared to values available in literature. Bootstrap analyses are performed to estimate the confidence intervals of the intensity magnitudes, as well as to quantify the location uncertainty. The analyses of seven significant earthquakes for the hazard assessment are presented in detail, including three large historical earthquakes that struck the northern Tien-Shan between the end of the nineteenth and the beginning of the twentieth centuries: the 1887, M 7.3 Verny, the 1889, M 8.3 Chilik and the 1911, M 8.2 Kemin earthquakes. Regarding the 1911, Kemin earthquake the magnitude values estimated from intensity data are lower (i.e. MILH?=?7.8 and MIW?=?7.6 considering surface wave and moment magnitude, respectively) than the value M?=?8.2 listed in the considered catalog. These values are more in agreement with the value M S?=?7.8 revised by Abe and Noguchi (Phys Earth Planet In, 33:1–11, 1983b) for the surface wave magnitude. For the Kemin earthquake, the distribution of the bootstrap solutions for the intensity centre reveal two minima, indicating that the distribution of intensity assignments do not constrain a unique solution. This is in agreement with the complex source rupture history of the Kemin earthquake, which involved several fault segments with different strike orientations, dipping angles and focal mechanisms (e.g. Delvaux et al. in Russ Geol Geophys 42:1167–1177, 2001; Arrowsmith et al. in Eos Trans Am Geophys Union 86(52), 2005). Two possible locations for the intensity centre are obtained. The first is located on the easternmost sub-faults (i.e. the Aksu and Chon-Aksu segments), where most of the seismic moment was released (Arrowsmith et al. in Eos Trans Am Geophys Union 86(52), 2005). The second location is located on the westernmost sub-faults (i.e. the Dzhil'-Aryk segment), close to the intensity centre location obtained for the 1938, M 6.9 Chu-Kemin earthquake (MILH?=?6.9 and MIW?=?6.8).  相似文献   

17.
A previous study (Lyman et al., Nature 465:334–337, 2010) showed a robust warming signal of the global upper ocean (0–700 m). They examined several sources of uncertainty that contribute to differences among heat content estimations. However, their focus was limited to globally averaged estimation. This study presents the spatial pattern of the global heat content change based on observed gridded datasets (Levitus et al., Geophys Res Lett 36:L07608, 2009). The western Pacific, Atlantic, and Indian Oceans showed significant warming trends, whereas eastern Pacific and some areas of the Gulf Stream experienced negative trends during 1993–2009. Steady warming trend was obtained from the first EOF mode when El Nino and Southern Oscillation (ENSO)-related signals were removed. This result implies that the rapid increase in heat content of the upper ocean around 2000–2005 is not related to a sampling transition from XBT to Argo observations but is associated with a natural variability dominated by strong ENSO-related signals.  相似文献   

18.
Earthquake early warning systems (EEWS) are considered to be an effective, pragmatic, and viable tool for seismic risk reduction in cities. While standard EEWS approaches focus on the real-time estimation of an earthquake’s location and magnitude, innovative developments in EEWS include the capacity for the rapid assessment of damage. Clearly, for all public authorities that are engaged in coordinating emergency activities during and soon after earthquakes, real-time information about the potential damage distribution within a city is invaluable. In this work, we present a first attempt to design an early warning and rapid response procedure for real-time risk assessment. In particular, the procedure uses typical real-time information (i.e., P-wave arrival times and early waveforms) derived from a regional seismic network for locating and evaluating the size of an earthquake, information which in turn is exploited for extracting a risk map representing the potential distribution of damage from a dataset of predicted scenarios compiled for the target city. A feasibility study of the procedure is presented for the city of Bishkek, the capital of Kyrgyzstan, which is surrounded by the Kyrgyz seismic network by mimicking the ground motion associated with two historical events that occurred close to Bishkek, namely the 1911 Kemin (M?=?8.2; ±0.2) and the 1885 Belovodsk (M?=?6.9; ±0.5) earthquakes. Various methodologies from previous studies were considered when planning the implementation of the early warning and rapid response procedure for real-time risk assessment: the Satriano et al. (Bull Seismol Soc Am 98(3):1482–1494, 2008) approach to real-time earthquake location; the Caprio et al. (Geophys Res Lett 38:L02301, 2011) approach for estimating moment magnitude in real time; the EXSIM method for ground motion simulation (Motazedian and Atkinson, Bull Seismol Soc Am 95:995–1010, 2005); the Sokolov (Earthquake Spectra 161: 679–694, 2002) approach for estimating intensity from Fourier amplitude spectra; and the Tyagunov et al. (Nat Hazard Earth Syst Sci 6:573–586, 2006) approach for risk computation. Innovatively, all these methods are jointly applied to assess in real time the seismic risk of a particular target site, namely the city of Bishkek. Finally, the site amplification and vulnerability datasets considered in the proposed methodology are taken from previous studies, i.e., Parolai et al. (Bull Seismol Soc Am, 2010) and Bindi et al. (Soil Dyn Earthq Eng, 2011), respectively.  相似文献   

19.
The Princeton Ocean Model (POM) with generalized coordinate system (POMgcs) is used to study the summer surface-layer thermal response to surface gravity waves in the Yellow Sea (YS). The parameterization schemes of wave breaking developed by Mellor and Blumberg (J Phys Oceanogr 34:693–698, 2004) and Kantha and Clayson (Ocean Model 6:101–124, 2004), respectively, and Stokes production developed by Kantha and Clayson (Ocean Model 6:101–124, 2004) are both included in the Mellor–Yamada turbulence closure model Mellor and Yamada (Rev Geophys 20:851–875, 1982) of POMgcs. Numerical results show that surface gravity waves impact the depth of surface mixed layer of temperature in the YS in summer. The surface mixed layer in the YS cannot be reproduced well and has a visible difference from the observation if the parameterization schemes are not included. A diagnostic analysis of turbulent kinetic energy suggests that both Stokes production and wave breaking play key roles in enhancing the turbulent mixing near the sea surface in the YS. Stokes production seems to have a greater impact throughout the upper mixed layer in the YS in summer than that of wave breaking. In addition, a diagnostic analysis of the momentum balance shows that Coriolis–Stokes forcing has a significant effect on the momentum budget in the upper layer in the YS, and surface gravity waves are able to reduce the velocity of mean flow near the surface and make the mean flow near the surface more homogeneous vertically in the YS.  相似文献   

20.
To investigate the nature of Pn propagation, we have implemented the spectral-element method (SEM) for vertically and laterally varying media with and without attenuation. As a practical measure, essential features of the Pn waves are distilled into seismic attributes including arrival times, amplitudes and pulse frequencies. To validate the SEM simulations, we first compare the SEM results with reflectivity calculations of Braile and Smith (Geophys. J.R. Astr. Soc. 40, 145–176, 1975) and then to the asymptotic results of ?erveny and Ravindra (Theory of Seismic Headwaves, University of Toronto Press, pp. 235–250, 1971). Models with random, laterally varying Moho structures are then simulated, where the amplitude and pulse frequency characteristics are found to be stable to small Moho interface perturbations. SEM calculations for models with different upper-mantle velocity gradients are next performed where it is found that interference effects can strongly influence the Pn amplitudes and pulses frequencies. For larger-scale, laterally varying structures, SEM models similar to that found along the Hi-CLIMB array in Tibet are then performed. It is observed that large-scale structures, along with small-scale structures, upper-mantle velocity gradients and attenuation, can all significantly affect the Pn attributes. Ambiguities between upper-mantle velocity gradients and attenuation are also found when using Pn amplitudes and pulse frequency attributes. These ambiguities may be resolved, to some degree, by using the curvature of the travel times at longer regional distance, however, this would also be complicated by lateral variability.  相似文献   

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

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