Pure and Applied Geophysics

Douglas, John; Gehl, Pierre; Bonilla, Luis; Gélis, Céline

doi: 10.1007/s00024-010-0146-5pmid: N/A

An important parameter for the characterization of strong ground motion at high-frequencies (>1 Hz) is kappa, κ, which models a linear decay of the acceleration spectrum, a(f), in log-linear space (i.e. a(f) = A 0 exp(− π κ f) for f > f E where f is frequency, f E is a low frequency limit and A 0 controls the amplitude of the spectrum). κ is a key input parameter in the stochastic method for the simulation of strong ground motion, which is particularly useful for areas with insufficient strong-motion data to enable the derivation of robust empirical ground motion prediction equations, such as mainland France. Numerous studies using strong-motion data from western North America (WNA) (an active tectonic region where surface rock is predominantly soft) and eastern North America (ENA) (a stable continental region where surface rock is predominantly very hard) have demonstrated that κ varies with region and surface geology, with WNA rock sites having a κ of about 0.04 s and ENA rock sites having a κ of about 0.006 s. Lower κs are one reason why high-frequency strong ground motions in stable regions are generally higher than in active regions for the same magnitude and distance. Few, if any, estimates of κs for French sites have been published. Therefore, the purpose of this study is to estimate κ using data recorded by the French national strong-motion network (RAP) for various sites in different regions of mainland France. For each record, a value of κ is estimated by following the procedure developed by Anderson and Hough (Bull Seismol Soc Am 74:1969–1993, 1984): this method is based on the analysis of the S-wave spectrum, which has to be performed manually, thus leading to some uncertainties. For the three French regions where most records are available (the Pyrenees, the Alps and the Côtes-d’Azur), a regional κ model is developed using weighted regression on the local geology (soil or rock) and source-to-site distance. It is found that the studied regions have a mean κ between the values found for WNA and ENA. For example, for the Alps region a κ value of 0.0254 s is found for rock sites, an estimate reasonably consistent with previous studies.

Chopra, Sumer; Kumar, Dinesh; Rastogi, Bal

doi: 10.1007/s00024-010-0132-ypmid: N/A

The strong ground motions for the 2001 Bhuj (M w 7.6) India earthquake have been estimated on hard rock and B/C boundary (NEHRP) levels using a recently modified version of stochastic finite fault modeling based on dynamic corner frequency (Motazedian and Atkinson in Bull Seismol Soc Am 95, 995–1010 2005). Incorporation of dynamic corner frequency removes the limitations of earlier stochastic methods. Simulations were carried out at 13 sites in Gujarat where structural response recorder (SRR) recordings are available. In addition, accelerograms were simulated at the B/C boundary at a large number of points distributed on a grid. The corresponding response spectra have also been estimated. The values of peak ground accelerations and spectral accelerations at three periods (0.4, 0.75 and 1.25 s) are presented in the form of contour maps. The maximum value of peak ground acceleration (PGA) in the center of meizoseismal zone is 550 cm/s2. The response spectral acceleration in same zone is 900 cm/s2 (T = 0.4 s), 600 cm/s2 (T = 0.75 s) and 300 cm/s2 (T = 1.25 s). The innermost PGA contour is on the fault plane. A comparison of the PGA values obtained at 13 sites in this study with those obtained in earlier studies on the same sites, but employing different methods, show that the present PGA values are comparable at most of the sites. The rate of decay of PGA values is fast at short distances as compared to that at longer distances. The PGA values obtained here put some constraints on the expected values from a similar earthquake in the region. A synthetic intensity map has been prepared from the estimated values of PGA using an empirical relation. A comparison with the reported intensity map of the earthquake shows the synthetic MMI values, as expected, are lower by 1 unit compared to reported intensity map. The contour map of PGA along with the contour maps of spectral acceleration at various periods permit the assessment of damage potential to various categories of houses and other structures. Such information will be quite important in planning of mitigation and disaster management programs in the region.

Yadav, Ram; Tripathi, Jayant; Rastogi, Bal; Das, Mridul; Chopra, Sumer

doi: 10.1007/s00024-010-0105-1pmid: N/A

Northeast India and adjoining regions (20°–32° N and 87°–100° E) are highly vulnerable to earthquake hazard in the Indian sub-continent, which fall under seismic zones V, IV and III in the seismic zoning map of India with magnitudes M exceeding 8, 7 and 6, respectively. It has experienced two devastating earthquakes, namely, the Shillong Plateau earthquake of June 12, 1897 (M w 8.1) and the Assam earthquake of August 15, 1950 (M w 8.5) that caused huge loss of lives and property in the Indian sub-continent. In the present study, the probabilities of the occurrences of earthquakes with magnitude M ≥ 7.0 during a specified interval of time has been estimated on the basis of three probabilistic models, namely, Weibull, Gamma and Lognormal, with the help of the earthquake catalogue spanning the period 1846 to 1995. The method of maximum likelihood has been used to estimate the earthquake hazard parameters. The logarithmic probability of likelihood function (ln L) is estimated and used to compare the suitability of models and it was found that the Gamma model fits best with the actual data. The sample mean interval of occurrence of such earthquakes is estimated as 7.82 years in the northeast India region and the expected mean values for Weibull, Gamma and Lognormal distributions are estimated as 7.837, 7.820 and 8.269 years, respectively. The estimated cumulative probability for an earthquake M ≥ 7.0 reaches 0.8 after about 15–16 (2010–2011) years and 0.9 after about 18–20 (2013–2015) years from the occurrence of the last earthquake (1995) in the region. The estimated conditional probability also reaches 0.8 to 0.9 after about 13–17 (2008–2012) years in the considered region for an earthquake M ≥ 7.0 when the elapsed time is zero years. However, the conditional probability reaches 0.8 to 0.9 after about 9–13 (2018–2022) years for earthquake M ≥ 7.0 when the elapsed time is 14 years (i.e. 2009).

Gusev, Alexander

doi: 10.1007/s00024-010-0098-9pmid: N/A

A wavetrain of high-frequency (HF) P waves from a large earthquake, when recorded at a distant station, looks like a segment of modulated noise, with its duration close to the duration of rupture. These wavetrains, with their bursts and fadings, look much more intermittent than a segment of common stationary random noise. We try to describe quantitatively this bursty behavior. To this end, variogram and spectral analyses are applied to time histories of P-wave envelopes (squared-amplitude or instant-power signals) in six HF bands of 1-Hz width. Nine M w = 7.6–9.2 earthquakes were examined, using, in total, 232 records and 992 single-band traces. Variograms of integrated instant power are approximately linear on a log–log scale, indicating that the correlation structure of the instant-power signal is approximately self-similar. Also, estimates of the power spectrum of the instant-power signal look approximately linear on a log–log scale. Log–log slopes of the variograms and spectra deliver estimates of the Hurst exponent H that are mostly in the range 0.6–0.9, markedly above the value H = 0.5 of uncorrelated (white-noise) signals. The preferred estimate over the entire data set is H = 0.83, still, this estimate may include some bias, and must be treated as preliminary. The inter-event scatter of H estimates is about 0.04, reflecting individual event-to-event variations of H. Many of the average log–log spectral plots show slight concavity that perturbs the approximately linear slope; this is a secondary effect that seems to be mostly related to the limited bandwidth of the data. Evidence is given in support of the idea that the observed approximately self-similar correlation structure of the P-wave envelope originates in a similar structure of the body wave instant-power signal radiated by the source, so that the propagation-related distortions can be regarded as limited. The facts presented suggest that the space–time organization of the earthquake rupture process is multiscaled and bears significant fractal features; it deviates from the brittle-crack model with its two well-separated characteristic scales. Phenomenologically, the high-frequency body-wave radiation from an earthquake source can be thought of as a product of stationary noise and the square root of a positive random envelope function with a power-law spectrum. From the viewpoint of applications, the self-similarity of body wave envelopes provides a useful constraint for earthquake source models used to simulate strong ground motions.

Zhou, Hong; Chen, Xiaofei

doi: 10.1007/s00024-010-0106-0pmid: N/A

The classic spectral element method (SEM) is important for seismographical simulation. However, waves subjected to irregular interfaces or surfaces are difficult to simulate accurately using SEM with quadrangular/hexahedral elements. In this paper we propose a new technique to solve this problem. The technique reconstructs some new elements near the surface/interface to substitute for any element crossing the interface, thus making the boundary of some new elements an accurate fit to the interface/surface. Numerical comparisons with the classic SEM show that the technique has improved flexibility when dealing with interface problems without losing accuracy and efficiency. The technique also enables us to vary the size and shape of an element smoothly as the velocity in a medium varies so that this removes the inaccuracy resulting from the high local variation of the grid in the classic SEM. Therefore, the technique widens the application of the classic SEM in seismographic simulation.

Bai, Chao-ying; Zhao, Rui; Greenhalgh, Stewart

doi: 10.1007/s00024-010-0102-4pmid: N/A

A novel hybrid approach to earthquake location is proposed which uses a combined coarse global search and fine local inversion with a minimum search routine. The method exploits the advantages of network ray tracing and robust formulation of the Fréchet derivatives to simultaneously update all sampled initial source parameters in the solution space to determine the best solution. Synthetic examples, involving a three-dimensional (3-D) complex velocity model and a challenging source–receiver layout, are used to demonstrate the advantages over direct grid search algorithms in terms of solution accuracy, computational efficiency, and sensitivity to noise. Therefore, this is a promising scheme for earthquake early warning, tsunami early warning, rapid hazard assessment, and emergency response after strong earthquake occurrence.

Ferahtia, Jalal; Baddari, Kamel; Djarfour, Nouredine; Kassouri, Abdel

doi: 10.1007/s00024-010-0101-5pmid: N/A

Seismic noise is a fundamental part of seismic data which cannot be avoided when conducting any seismic survey. It consists of coherent and random noise. Noise removal or filtering is one of the major concerns in the field of seismic processing. In this paper, we introduce an image filtering technique based on a detection-estimation algorithm for Gaussian and random noise removal in seismic data, namely the trilateral filter, based on a statistic called rank-ordered absolute differences. The non-linear and adaptive behaviour of this filter makes it very robust in the presence of random and coherent noise, in addition to its computational simplicity and its ability to automatically identify noise in data. We have modified the strategy of trilateral filtering by adapting the rank-ordered absolute differences formula in order to extract the signal component. We have successfully used this filter for the removal of surface waves and random spiky noise from synthetic and field data. Results are very encouraging and show the superiority of this filter compared with other filters, particularly when used recursively.

Gélis, C.; Revil, A.; Cushing, M.; Jougnot, D.; Lemeille, F.; Cabrera, J.; Hoyos, A.; Rocher, M.

doi: 10.1007/s00024-010-0097-xpmid: N/A

The Experimental platform of Tournemire (Aveyron, France) developed by IRSN (French Institute for Radiological Protection and Nuclear Safety) is located in a tunnel excavated in a clay–rock formation interbedded between two limestone formations. A well-identified regional fault crosscuts this subhorizontal sedimentary succession, and a subvertical secondary fault zone is intercepted in the clay–rock by drifts and boreholes in the tunnel at a depth of about 250 m. A 2D electrical resistivity survey was carried out along a 2.5 km baseline, and a takeout of 40 m was used to assess the potential of this method to detect faults from the ground surface. In the 300 m-thick zone investigated by the survey, electrical resistivity images reveal several subvertical low-resistivity discontinuities. One of these discontinuities corresponds to the position of the Cernon fault, a major regional fault. One of the subvertical conductive discontinuities crossing the upper limestone formation is consistent with the prolongation towards the ground surface of the secondary fault zone identified in the clay–rock formation from the tunnel. Moreover, this secondary fault zone corresponds to the upward prolongation of a subvertical fault identified in the lower limestone using a 3D high-resolution seismic reflection survey. This type of large-scale electrical resistivity survey is therefore a useful tool for identifying faults in superficial layers from the ground surface and is complementary to 3D seismic reflection surveys.

Dai, Feng; Xia, Kaiwen

doi: 10.1007/s00024-010-0103-3pmid: N/A

Granitic rocks usually exhibit strongly anisotropy due to pre-existing microcracks induced by long-term geological loadings. The understanding of the rock anisotropy in mechanical properties is critical to a variety of rock engineering applications. In this paper, Brazilian tests are conducted statically with a material testing machine and dynamically with a split Hopkinson pressure bar system to measure both static and dynamic tensile strength of Barre granite. To understand the anisotropy in tensile strength, samples are cored and labelled using the three principle directions of Barre granite to form six sample groups. For dynamic tests, a pulse shaping technique is used to achieve dynamic equilibrium in the samples during the dynamic test. The finite element method is then implemented to formulate equations that relate the failure load to the material tensile strength by employing an orthotropic elastic material model. For samples in the same orientation group, the tensile strength shows clear loading rate dependence. The tensile strengths also exhibit clear anisotropy under static loading while the anisotropy diminishes as the loading rate increases, which may be due to the interaction of pre-existing microcracks.

Giannopoulou, K.; Santamouris, M.; Livada, I.; Georgakis, C.; Caouris, Y.

doi: 10.1007/s00024-010-0099-8pmid: N/A

The present paper investigates the impact of canyon geometry on the temperature regime and nocturnal heat island development in the very dense urban area of Athens, Greece. Detailed measurements of air temperature have been carried out within three deep urban canyons of different aspect ratios (H/W = 3, 2.1 and 1.7) during the night period of the summer and autumn of 2007. An analysis was carried out to investigate the relative impact of the canyon geometry, the undisturbed wind velocity, ambient temperature, and cloud cover on the development of a nocturnal heat island. A clear increase of the median, maximum and minimum values of the cooling rates has been observed for decreasing aspect ratios. Under low ambient temperatures, high wind speeds correspond to a substantial rise of the cooling rate in the urban canyons mainly because of the increased convective losses. On the contrary, cooling rates decrease substantially under high-undisturbed wind speeds and ambient temperatures because of the important convective gains. The impact of cloud cover was found to be important as cloudy skies cause a substantial decrease of the cooling rates in the urban canyons. Comparisons were performed between the temperature data collected in the three studied urban canyons and temperatures recorded in an urban as well as a suburban open space station.