首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
In this paper, I introduce a novel approach to modelling the individual random component (also called the intra-event uncertainty) of a ground-motion relation (GMR), as well as a novel approach to estimating the corresponding parameters. In essence, I contend that the individual random component is reproduced adequately by a simple stochastic mechanism of random impulses acting in the horizontal plane, with random directions. The random number of impulses was Poisson distributed. The parameters of the model were estimated according to a proposal by Raschke J Seismol 17(4):1157–1182, (2013a), with the sample of random difference ξ?=?ln(Y 1 )-ln(Y 2 ), in which Y 1 and Y 2 are the horizontal components of local ground-motion intensity. Any GMR element was eliminated by subtraction, except the individual random components. In the estimation procedure, the distribution of difference ξ was approximated by combining a large Monte Carlo simulated sample and Kernel smoothing. The estimated model satisfactorily fitted the difference ξ of the sample of peak ground accelerations, and the variance of the individual random components was considerably smaller than that of conventional GMRs. In addition, the dependence of variance on the epicentre distance was considered; however, a dependence of variance on the magnitude was not detected. Finally, the influence of the novel model and the corresponding approximations on PSHA was researched. The applied approximations of distribution of the individual random component were satisfactory for the researched example of PSHA.  相似文献   

2.
The problem of estimating the time derivatives of the horizontal components of the geomagnetic field and forecasting the probability of the occurrence of perturbations that exceed a given threshold level (the over-threshold perturbations) arises in the applications concerned with the geomagnetically induced currents (GICs). In this work, we consider the temporal and spatial structure of the Pi3 pulsations with quasi-periods of 102 to 103 s during which the auroral and subauroral stations of the IMAGE network record over-threshold values in the derivatives of the meridional (along the longitudinal circle) BX component and latitudinal (along the latitudinal circle) BY component. The extreme |dBX/dt| values mainly develop against the background of the Pi3 pulsations with a complex frequency content, whereas the extreme |dBY/dt| values appear when the buildup (decay) phases of the bay-like disturbance associated with the evolution of a substorm coincide with the respective phases of the field of pulsations. The conditions under which the derivatives |dBX/dt| and |dBY/dt| reach their over-threshold values are studied for subauroral latitudes by the technique of superposed epoch analysis. The extreme values of the derivatives most frequently occur during the main phase of moderate magnetic storms or beyond the storm—during high substorm activity under the conditions of a negative vertical component of the interplanetary magnetic field. The probability of the occurrence of over-threshold values increases at high amplitudes of the Pi3 pulsations and depends on their spectral content. The problem of analyzing and forecasting the over-threshold |dBY/dt| perturbations is complicated by the fact that the scale of the perturbations is small along the lines of latitude and large along the meridians. This can result in GIC excitation in the North–South oriented electric power lines by the geomagnetic perturbations localized within a narrow band in longitude which can be missed during the measurements.  相似文献   

3.
Accurate and precise estimation of return levels is often a key goal of any extreme value analysis. For example, in the UK the British Standards Institution (BSI) incorporate estimates of ‘once-in-50-year wind gust speeds’—or 50-year return levels—into their design codes for new structures; similarly, the Dutch Delta Commission use estimates of the 10,000-year return level for sea-surge to aid the construction of flood defence systems. In this paper, we briefly highlight the shortcomings of standard methods for estimating return levels, including the commonly-adopted block maxima and peaks over thresholds approach, before presenting an estimation framework which we show can substantially increase the precision of return level estimates. Our work allows explicit quantification of seasonal effects, as well as exploiting recent developments in the estimation of the extremal index for handling extremal clustering. From frequentist ideas, we turn to the Bayesian paradigm as a natural approach for building complex hierarchical or spatial models for extremes. Through simulations we show that the return level posterior mean does not have an exceedance probability in line with the intended encounter risk; we also argue that the Bayesian posterior predictive value gives the most satisfactory representation of a return level for use in practice, accounting for uncertainty in parameter estimation and future observations. Thus, where feasible, we propose a Bayesian estimation strategy for optimal return level inference.  相似文献   

4.
Bayesian probability theory is an appropriate and useful method for estimating parameters in seismic hazard analysis. The analysis in Bayesian approaches is based on a posterior belief, also their special ability is to take into account the uncertainty of parameters in probabilistic relations and a priori knowledge. In this study, we benefited the Bayesian approach in order to estimate maximum values of peak ground acceleration (Amax) also quantiles of the relevant probabilistic distributions are figured out in a desired future interval time in Iran. The main assumptions are Poissonian character of the seismic events flow and properties of the Gutenberg-Richter distribution law. The map of maximum possible values of Amax and also map of 90% quantile of distribution of maximum values of Amax on a future interval time 100 years is presented. According to the results, the maximum value of the Amax is estimated for Bandar Abbas as 0.3g and the minimum one is attributed to Esfahan as 0.03g. Finally, the estimated values in Bayesian approach are compared with what was presented applying probabilistic seismic hazard (PSH) methods based on the conventional Cornel (1968) method. The distribution function of Amax for future time intervals of 100 and 475 years are calculated for confidence limit of probability level of 90%.  相似文献   

5.
We used CHAMP satellite vector data and the latest IGRF12 model to investigate the regional magnetic anomalies over mainland China. We assumed satellite points on the same surface (307.69 km) and constructed a spherical cap harmonic model of the satellite magnetic anomalies for elements X, Y, Z, and F over Chinese mainland for 2010.0 (SCH2010) based on selected 498 points. We removed the external field by using the CM4 model. The pole of the spherical cap is 36N° and 104°E, and its half-angle is 30°. After checking and comparing the root mean square (RMS) error of ΔX, ΔY, and ΔZ and X, Y, and Z, we established the truncation level at K max = 9. The results suggest that the created China Geomagnetic Referenced Field at the satellite level (CGRF2010) is consistent with the CM4 model. We compared the SCH2010 with other models and found that the intensities and distributions are consistent. In view of the variation of F at different altitudes, the SCH2010 model results obey the basics of the geomagnetic field. Moreover, the change rate of X, Y, and Z for SCH2010 and CM4 are consistent. The proposed model can successfully reproduce the geomagnetic data, as other data-fitting models, but the inherent sources of error have to be considered as well.  相似文献   

6.
In this paper we present a stochastic model reduction method for efficiently solving nonlinear unconfined flow problems in heterogeneous random porous media. The input random fields of flow model are parameterized in a stochastic space for simulation. This often results in high stochastic dimensionality due to small correlation length of the covariance functions of the input fields. To efficiently treat the high-dimensional stochastic problem, we extend a recently proposed hybrid high-dimensional model representation (HDMR) technique to high-dimensional problems with multiple random input fields and integrate it with a sparse grid stochastic collocation method (SGSCM). Hybrid HDMR can decompose the high-dimensional model into a moderate M-dimensional model and a few one-dimensional models. The moderate dimensional model only depends on the most M important random dimensions, which are identified from the full stochastic space by sensitivity analysis. To extend the hybrid HDMR, we consider two different criteria for sensitivity test. Each of the derived low-dimensional stochastic models is solved by the SGSCM. This leads to a set of uncoupled deterministic problems at the collocation points, which can be solved by a deterministic solver. To demonstrate the efficiency and accuracy of the proposed method, a few numerical experiments are carried out for the unconfined flow problems in heterogeneous porous media with different correlation lengths. The results show that a good trade-off between computational complexity and approximation accuracy can be achieved for stochastic unconfined flow problems by selecting a suitable number of the most important dimensions in the M-dimensional model of hybrid HDMR.  相似文献   

7.
The Aki-Utsu method of Gutenberg-Richter (G-R) b value estimation is often misapplied so that estimations not using the G-R histogram are often meaningless because they are not based on adequate samples. We propose a method to estimate the likelihood Pr(b?b m , N, M 1, M 2) that an observed b m estimate, based on a sample of N magnitudes within an [M 1????≤?ΔM/2,?M 2?+?ΔM/2) range, where ΔM?=?0.1 is the usual rounding applied to magnitudes, is due to a “true” source b value, b, and use these likelihoods to estimate source b ranges corresponding to various confidence levels. As an example of application of the method, we estimate the b values before and after the occurrence of a 7.4-magnitude earthquake in the Mexican subduction zone, and find a difference of 0.82 between them with 100% confidence that the b values are different.  相似文献   

8.
The 2017 Guptkashi earthquake occurred in a segment of the Himalayan arc with high potential for a strong earthquake in the near future. In this context, a careful analysis of the earthquake is important as it may shed light on source and ground motion characteristics during future earthquakes. Using the earthquake recording on a single broadband strong-motion seismograph installed at the epicenter, we estimate the earthquake’s location (30.546° N, 79.063° E), depth (H?=?19 km), the seismic moment (M0?=?1.12×1017 Nm, M w 5.3), the focal mechanism (φ?=?280°, δ?=?14°, λ?=?84°), the source radius (a?=?1.3 km), and the static stress drop (Δσ s ~22 MPa). The event occurred just above the Main Himalayan Thrust. S-wave spectra of the earthquake at hard sites in the arc are well approximated (assuming ω?2 source model) by attenuation parameters Q(f)?=?500f0.9, κ?=?0.04 s, and fmax?=?infinite, and a stress drop of Δσ?=?70 MPa. Observed and computed peak ground motions, using stochastic method along with parameters inferred from spectral analysis, agree well with each other. These attenuation parameters are also reasonable for the observed spectra and/or peak ground motion parameters in the arc at distances ≤?200 km during five other earthquakes in the region (4.6?≤?M w ?≤?6.9). The estimated stress drop of the six events ranges from 20 to 120 MPa. Our analysis suggests that attenuation parameters given above may be used for ground motion estimation at hard sites in the Himalayan arc via the stochastic method.  相似文献   

9.
Intense quasimonchromatic geomagnetic pulsations with a period of ～15 min, observed on the Earth’s surface in the near-noon sector at the beginning of the recovery phase of a very strong (Dst min = ?260 nT) magnetic storm of May 15, 2005, are analyzed. The variations were registered at auroral latitudes only in the X field component, and wave activity shifted into the postnoon sector of the polar cap an hour later; in this case pulsations were observed in the X and Y field components. Within the magnetosphere the source of magnetic pulsations could be the surface waves on the magnetopause caused by the pulse of the solar wind magnetic pressure. Geomagnetic pulsations in the polar cap, observed in phase at different latitudes, could apparently reflect quasiperiodic variations in the NBZ system of field-aligned currents. Such variations can originate due to the series of pulsed reconnections in the postnoon outer cusp at large (～20 nT) positive B z values and large (about ?40 nT) negative values of IMF B x .  相似文献   

10.
Actual evapotranspiration(ET_a) over the Tibetan Plateau(TP) is an important component of the water cycle,and greatly influences the water budgets of the TP lake basins.Quantitative estimation of ET_a within lake basins is fundamental to physically understanding ET_a changes,and thus will improve the understanding of the hydro logical processes and energy balance throughout the lake basins.In this study,the spatiotemporal dynamic changes of ET_α within the Lake Selin Co(the TP's largest lake) and its surrounding small lakes and land area during 2003-2012 are examined at the basin scale.This was carried out using the well-established Water and Energy Budget-based Distributed Hydrological Model(WEB-DHM) for the land area,the Penman method for the water area when unfrozen,and a simple sublimation estimation approach for the water area when frozen.The relationships between ET_a changes and controlling factors are also discussed.Results indicate that the simulated land ET_a from the WEB-DHM reasonably agrees with the estimated ET_a values from the nonlinear complementary relationship model using appropriately calibrated parameter values at a point scale.Land ET_a displayed a non-significant increase of 7.03 mm year~(-1),and largely depends on precipitation.For the water area,the combined effects of reduced wind speed and net radiation offset the effect of rising temperature and vapor pressure deficit,and contributed to a non-significant decrease in evaporation of 4.17 mm year~(-1).Sensitivity analysis shows that vapor pressure deficit and wind speed are the most sensitive variables to the changes of evaporation from the water area.  相似文献   

11.
The purpose of this work is to define a seismic regionalization of Mexico for seismic hazard and risk analyses. This seismic regionalization is based on seismic, geologic, and tectonic characteristics. To this end, a seismic catalog was compiled using the more reliable sources available. The catalog was made homogeneous in magnitude in order to avoid the differences in the way this parameter is reported by various agencies. Instead of using a linear regression to converts from m b and M d to M s or M w , using only events for which estimates of both magnitudes are available (i.e., paired data), we used the frequency-magnitude relations relying on the a and b values of the Gutenberg-Richter relation. The seismic regions are divided into three main categories: seismicity associated with the subduction process along the Pacific coast of Mexico, in-slab events within the down-going COC and RIV plates, and crustal seismicity associated to various geologic and tectonic regions. In total, 18 seismic regions were identified and delimited. For each, the a and b values of the Gutenberg-Richter relation were determined using a maximum likelihood estimation. The a and b parameters were repeatedly estimated as a function of time for each region, in order to confirm their reliability and stability. The recurrence times predicted by the resulting Gutenberg-Richter relations obtained are compared with the observed recurrence times of the larger events in each region of both historical and instrumental earthquakes.  相似文献   

12.
Aftershock sequences of some strong earthquakes of Kamchatka, the Kurile Islands, and Japan are examined. Such source parameters as the length L, along-dip width W, motion on fault D, and stress drop Δσ are determined from the aftershock sequences considered. The values of these parameters were obtained by the formal estimation of linear source parameters (lower bound estimates) and visually (upper bound estimates). The correlation dependences of the obtained parameters on the surface wave (M S ) and seismic moment (M W ) magnitudes are calculated.  相似文献   

13.
In this study a new method is presented to determine model parameters from magnetic anomalies caused by dipping dikes. The proposed method is applied by employing only the even component of the anomaly. First, the maximum of the even component is divided to its value at any distance x in order to obtain S1. Then, theoretical even component values are computed for the minimal depth (h) and half-width (b) values. S2 is obtained by dividing their maximum to the value computed for the same distance x. A set of S2 values is calculated by slowly increasing the half-width, and h and b for the S2 closest to S1 are determined. The same procedure is repeated by increasing the depth. The determined b values are plotted against the corresponding values of h. After repeating the process and plotting curves for different distances, it is possible to determine the actual depth and half-width values.  相似文献   

14.
Recent publications on the regression between earthquake magnitudes assume that both magnitudes are affected by error and that only the ratio of error variances is known. If X and Y represent observed magnitudes, and x and y represent the corresponding theoretical values, the problem is to find the a and b of the best-fit line \(y = a x + b\). This problem has a closed solution only for homoscedastic errors (their variances are all equal for each of the two variables). The published solution was derived using a method that cannot provide a sum of squares of residuals. Therefore, it is not possible to compare the goodness of fit for different pairs of magnitudes. Furthermore, the method does not provide expressions for the x and y. The least-squares method introduced here does not have these drawbacks. The two methods of solution result in the same equations for a and b. General properties of a discussed in the literature but not proved, or proved for particular cases, are derived here. A comparison of different expressions for the variances of a and b is provided. The paper also considers the statistical aspects of the ongoing debate regarding the prediction of y given X. Analysis of actual data from the literature shows that a new approach produces an average improvement of less than 0.1 magnitude units over the standard approach when applied to \(M_{w}\) vs. \(m_{b}\) and \(M_{w}\) vs. \(M_{S}\) regressions. This improvement is minor, within the typical error of \(M_{w}\). Moreover, a test subset of 100 predicted magnitudes shows that the new approach results in magnitudes closer to the theoretically true magnitudes for only 65 % of them. For the remaining 35 %, the standard approach produces closer values. Therefore, the new approach does not always give the most accurate magnitude estimates.  相似文献   

15.
In this study, continuous wavelet transform is applied to estimate the frequency-dependent quality factor of shear waves, Q S , in northwestern Iran. The dataset used in this study includes velocigrams of more than 50 events with magnitudes between 4.0 and 6.5, which have occurred in the study area. The CWT-based method shows a high-resolution technique for the estimation of S-wave frequency-dependent attenuation. The quality factor values are determined in the form of a power law as Q S (f)?=?(147?±?16)f 0.71?±?0.02 and (126?±?12)f 0.73?±?0.02 for vertical and horizontal components, respectively, where f is between 0.9 and 12 Hz. Furthermore, in order to verify the reliability of the suggested Q S estimator method, an additional test is performed by using accelerograms of Ahar-Varzaghan dual earthquakes on August 11, 2012, of moment magnitudes 6.4 and 6.3 and their aftershocks. Results indicate that the estimated Q S values from CWT-based method are not very sensitive to the numbers and types of waveforms used (velocity or acceleration).  相似文献   

16.
A key aim of most extreme value analyses is the estimation of the r-year return level; the wind speed, or sea-surge, or rainfall level (for example), we might expect to see once (on average) every r years. There are compelling arguments for working within the Bayesian setting here, not least the natural extension to prediction via the posterior predictive distribution. Indeed, for practitioners the posterior predictive return level has been cited as perhaps the most useful point summary from a Bayesian analysis of extremes, and yet little is known of the properties of this statistic. In this paper, we attempt to assess the performance of predictive return levels relative to their estimative counterparts obtained directly from the return level posterior distribution; in particular, we make comparisons with the return level posterior mean, mode and 95% credible upper bound. Differences between the predictive return level and standard summaries from the return level posterior distribution, for wind speed extremes observed in the UK, motivates this work. A large scale simulation study then reveals the superiority of the predictive return level over the other posterior summaries in many cases of practical interest.  相似文献   

17.
To study the prospective areas of upcoming strong-to-major earthquakes, i.e., M w  ≥ 6.0, a catalog of seismicity in the vicinity of the Thailand-Laos-Myanmar border region was generated and then investigated statistically. Based on the successful investigations of previous works, the seismicity rate change (Z value) technique was applied in this study. According to the completeness earthquake dataset, eight available case studies of strong-to-major earthquakes were investigated retrospectively. After iterative tests of the characteristic parameters concerning the number of earthquakes (N) and time window (T w ), the values of 50 and 1.2 years, respectively, were found to reveal an anomalous high Z-value peak (seismic quiescence) prior to the occurrence of six out of the eight major earthquake events studied. In addition, the location of the Z-value anomalies conformed fairly well to the epicenters of those earthquakes. Based on the investigation of correlation coefficient and the stochastic test of the Z values, the parameters used here (N = 50 events and T w  = 1.2 years) were suitable to determine the precursory Z value and not random phenomena. The Z values of this study and the frequency-magnitude distribution b values of a previous work both highlighted the same prospective areas that might generate an upcoming major earthquake: (i) some areas in the northern part of Laos and (ii) the eastern part of Myanmar.  相似文献   

18.
Accelerometric data from the well-studied valley EUROSEISTEST are used to investigate ground motion uncertainty and variability. We define a simple local ground motion prediction equation (GMPE) and investigate changes in standard deviation (σ) and its components, the between-event variability (τ) and within-event variability (φ). Improving seismological metadata significantly reduces τ (30–50%), which in turn reduces the total σ. Improving site information reduces the systematic site-to-site variability, φ S2S (20–30%), in turn reducing φ, and ultimately, σ. Our values of standard deviations are lower than global values from literature, and closer to path-specific than site-specific values. However, our data have insufficient azimuthal coverage for single-path analysis. Certain stations have higher ground-motion variability, possibly due to topography, basin edge or downgoing wave effects. Sensitivity checks show that 3 recordings per event is a sufficient data selection criterion, however, one of the dataset’s advantages is the large number of recordings per station (9–90) that yields good site term estimates. We examine uncertainty components binning our data with magnitude from 0.01 to 2 s; at smaller magnitudes, τ decreases and φ SS increases, possibly due to κ and source-site trade-offs Finally, we investigate the alternative approach of computing φ SS using existing GMPEs instead of creating an ad hoc local GMPE. This is important where data are insufficient to create one, or when site-specific PSHA is performed. We show that global GMPEs may still capture φ SS , provided that: (1) the magnitude scaling errors are accommodated by the event terms; (2) there are no distance scaling errors (use of a regionally applicable model). Site terms (φ S2S ) computed by different global GMPEs (using different site-proxies) vary significantly, especially for hard-rock sites. This indicates that GMPEs may be poorly constrained where they are sometimes most needed, i.e., for hard rock.  相似文献   

19.
An explanation is given to the fact that the cumulative number of damage cases is a decreasing function of the form Y?2/3. This is because the inundated area S depends on precipitation volume V as S ~ V2/3. Such dependence is confirmed by the data on the area of mushroom-shaped plumes at river mouths in the sites of river influx into the ocean.  相似文献   

20.
The hydrogen isotopic composition(δD) of leaf wax long-chain n-alkanes(C27, C29, and C31) from lacustrine sediments has been widely applied to reconstruct terrestrial paleoclimatic and paleohydrological changes. However, few studies have addressed whether the aquatic-derived n-alkanes can affect the δD values of lake sedimentary long-chain n-alkanes, which are usually regarded as a recorder of the terrestrial hydrological signals. Here we systematically investigated δD values of long-chain n-alkanes from modern aquatic plants, both near-shore and off-shore surface sediments, surrounding terrestrial plant litters, as well as river water and lake water in Lake Qinghai and its satellite lakes on the northeastern Qinghai-Tibet Plateau. Our data showed that(i) δD values of long-chain n-alkanes from aquatic plants varied from-184‰ to-132‰ for n-C27, from-183‰ to-138‰ for n-C29, and from-189‰ to-130‰ for n-C31, respectively, with no significant differences among the three n-alkanes homologues;(ii) δD values of long-chain n-alkanes from aquatic plants were generally more positive than those from surrounding terrestrial plants, possibly because that they recorded the D-enrichment of lake water in this semi-arid region;(iii) δD values of long-chain n-alkanes from surface sediments showed significant differences among the three n-alkanes homologues, due to the larger aquatic input of n-C27 to the sedimentary lipid pool than that of n-C31, and(iv) n-C27 δD values of near-shore aquatic plants and near-shore sediments are more negative than those from off-shore as a result of lower δD values of near-shore lake water. Our findings indicate that in this region(i) the offset between sedimentary n-C27 and n-C31 δD values(ΔδDC27-C31) could potentially be used to evaluate if sedimentary long-chain n-alkanes are derived from a single source;(ii) while δD values of n-C27 may be influenced by lake water hydrological changes, sedimentary n-C31 is derived predominantly from terrestrial plants and thus its δD can serve as a relatively reliable indicator for terrestrial paleoclimatic and paleohydrological reconstructions.  相似文献   

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

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