Pure and Applied Geophysics

Carcione, José; Poletto, Flavio

doi: 10.1007/s00024-013-0643-4pmid: N/A

It is well established that the upper—cooler—part of the crust is brittle, while deeper zones present ductile behaviour. In some cases, this brittle–ductile transition is a single seismic reflector with an associated reflection coefficient. We first develop a stress–strain relation including the effects of crust anisotropy, seismic attenuation and ductility in which deformation takes place by shear plastic flow. Viscoelastic anisotropy is based on the eigenstrain model and the Zener and Burgers mechanical models are used to model the effects of seismic attenuation, velocity dispersion, and steady-state creep flow, respectively. The stiffness components of the brittle and ductile media depend on stress and temperature through the shear viscosity, which is obtained by the Arrhenius equation and the octahedral stress criterion. The P- and S-wave velocities decrease as depth and temperature increase due to the geothermal gradient, an effect which is more pronounced for shear waves. We then obtain the reflection and transmission coefficients of a single brittle–ductile interface and of a ductile thin layer. The PP scattering coefficient has a Brewster angle (a sign change) in both cases, and there is substantial PS conversion at intermediate angles. The PP coefficient is sensitive to the layer thickness, unlike the SS coefficient. Thick layers have a well-defined Brewster angle and show higher reflection amplitudes. Finally, we compute synthetic seismograms in a homogeneous medium as a function of temperature.

Saleh, Salah

doi: 10.1007/s00024-013-0673-ypmid: N/A

Crustal and lithospheric thicknesses of the southeastern Mediterranean Basin region were determined using 3D Bouguer and elevation data analysis. The model is based on the assumption of local isostatic equilibrium. The calculated regional and residual Bouguer anomaly maps were employed for highlighting both deep and shallow structures. Generally, the regional field in the area under study is considered to be mainly influenced by the density contrast between the crust and upper mantle. Use of the gravity and topographic data with earthquake focal depths has improved both the geometry and the density distribution in the 3-D calculated profiles. The oceanic-continental boundary, the basement relief, Moho depth and lithosphere-asthenosphere boundary maps were estimated. The results point to the occurrence of thick continental crust areas with a thickness of approximately 32 km in northern Egypt. Below the coastal regions, the thickness of crust decreases abruptly (transition zone). An inverse correlation between sediment and crustal thicknesses shows up from the study. Furthermore, our density model reveals the existence of a continental crustal zone below the Eratosthenes Seamount block. Nevertheless, the crustal type beneath the Levantine basin is typically oceanic; this is covered by sedimentary sequences more than 14 km thick. The modeled Moho map shows a depth of 28–30 km below Cyprus and a depth of 26–28 km beneath the south Florence Rise in the northern west. However, the Moho lies at a constant shallow depth of 22–24 km below the Levantine Basin, which indicates thinning of the crust beneath this region. The Moho map reveals also a maximum depth of about 33–35 km beneath both the northern Egypt and northern Sinai, both of which are of the continental crust. The resulting mantle density anomalies suggest important variations of the lithosphere-asthenosphere boundary (LAB) topography, indicating prominent lithospheric mantle thinning beneath south Cyprus (LAB ~90 km depth), followed by thickening beneath the Eratosthenes seamount, Florence Rise, Levantine Basin and reaching to maximum thickness below Cyprian Arc (LAB ~115–120 km depth), and further followed by thinning in the north African margin plate and north Sinai subplate (LAB ~90–95 km depth). According to our density model profiles, we find that almost all earthquakes in the study area occurred along the western and central segments of the Cyprian arc while they almost disappear along the eastern segment. The active subduction zone in the Cyprian Arc is associated with large negative anomalies due to its low velocity upper mantle zone, which might be an indication of a serpentinized mantle. This means that collision between Cyprus and the Eratosthenes Seamount block is marked by seismic activity. Additionally, this block is in the process of dynamically subsiding, breaking-up and being underthrusted beneath Cyprus to the north and thrusted onto the Levantine Basin to the south.

Park, Eunjin; Ha, Wansoo; Chung, Wookeen; Shin, Changsoo; Min, Dong-Joo

doi: 10.1007/s00024-013-0651-4pmid: N/A

The wavefield in the Laplace domain has a very small amplitude except only near the source point. In order to deal with this characteristic, the logarithmic objective function has been used in many Laplace domain inversion studies. The Laplace-domain waveform inversion using the logarithmic objective function has fewer local minima than the time- or frequency domain inversion. Recently, the power objective function was suggested as an alternative to the logarithmic objective function in the Laplace domain. Since amplitudes of wavefields are very small generally, a power <1 amplifies the wavefields especially at large offset. Therefore, the power objective function can enhance the Laplace-domain inversion results. In previous studies about synthetic datasets, it is confirmed that the inversion using a power objective function shows a similar result when compared with the inversion using a logarithmic objective function. In this paper, we apply an inversion algorithm using a power objective function to field datasets. We perform the waveform inversion using the power objective function and compare the result obtained by the logarithmic objective function. The Gulf of Mexico dataset is used for the comparison. When we use a power objective function in the inversion algorithm, it is important to choose the appropriate exponent. By testing the various exponents, we can select the range of the exponent from 5 × 10−3 to 5 × 10−8 in the Gulf of Mexico dataset. The results obtained from the power objective function with appropriate exponent are very similar to the results of the logarithmic objective function. Even though we do not get better results than the conventional method, we can confirm the possibility of applying the power objective function for field data. In addition, the power objective function shows good results in spite of little difference in the amplitude of the wavefield. Based on these results, we can expect that the power objective function will produce good results from the data with a small amplitude difference. Also, it can partially be utilized at the sections where the amplitude difference is very small.

Kumar, Parveen; Joshi, A.; Verma, O.

doi: 10.1007/s00024-013-0658-xpmid: N/A

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.

Akyol, Nİhal; Kurtulmuş, Tevfİk; Çamyildiz, Murat; Güngör, Talİp

doi: 10.1007/s00024-013-0661-2pmid: N/A

Recordings of micro- and moderate-size local earthquakes have been used to quantify site effects in the central-west Turkey which contains one of the world’s best examples of a rapid intra-continental extension with its high population and industrial potential. We analyzed 436 earthquakes with local magnitudes ranging between 2.0 and 5.6 using three component digital recordings from 32 stations. Site functions were obtained using two different spectral ratio approaches (horizontal to vertical spectral ratio, HVSR, and standard spectral ratio, SSR). HVSR estimates of transverse and radial S-waves were compared with one another. Epicentral distance, magnitude and back-azimuth dependencies of site functions were also evaluated. In general, HVSR values from transverse and radial S-waves are similar within a factor of 2. The back-azimuth dependencies of transverse S-wave HVSR results are more significant than distance and magnitude dependencies. On the other hand, averaging of transverse and radial S-wave HVSR results eliminates systematic back-azimuth dependencies caused by source radiation effects. Distributions of HVSR estimates along ~N–S linear array, which traversed main grabens in the region with a station spacing of 3–4 km, reflect subsurface geological complexities in the region. The sites located near the basin edges are characterized by broader HVSR curves. Broad HVSR peaks could be attributed to the complexity of wave propagation related to significant 2D/3D velocity variations at the sediment–bedrock interface near the basin edges. The results also show that, even if the site is located on a horst, the presence of weathered zones along the surface could cause moderate frequency dependent site effects. Comparison of HVSR and SSR estimates for the stations on the graben sites showed that SSR estimates give larger values at lower frequencies which could be attributed to lateral variations in regional velocity and attenuation values caused by basin geometry and edge effects.

Sharma, Babita; Chopra, Sumer; Sutar, Anup; Bansal, B.

doi: 10.1007/s00024-013-0647-0pmid: N/A

In the present study ground motions for a Mw 8.5 scenario earthquake are estimated at 13 sites in Kumaun-Garhwal region using the empirical Green’s function technique. The recordings of 1991 Uttarkashi earthquake of Mw 6.8 at these sites are used as an element earthquake. A heterogeneous source model consisting of two asperities is considered for simulating the ground motions. The entire central seismic gap (CSG) can expect acceleration in excess of 100 cm/s2 with NW portion in excess of 400 cm/s2 and SE between 100 and 200 cm/s2. The central portion can expect peak ground acceleration (PGA) between 200 and 400 cm/s2. It has been observed from simulation of strong ground motion that sites located near the rupture initiation point can expect accelerations in excess of 1g. In the present analysis, Bhatwari and Uttarkashi can expect ground accelerations in excess of 1g. The estimates of the PGA are compared with earlier studies in the same region using different methodologies and it was found that the results are comparable. This has put constrains on the expected PGAs in this region. The obtained PGA values can be used in identifying the vulnerable areas in the central Himalaya, thereby facilitating the planning, design and construction of new structures and strengthening of the existing structures in the region.

Mandal, H.; Shukla, A.; Khan, P.; Mishra, O.

doi: 10.1007/s00024-013-0666-xpmid: N/A

The Son-Narmada-Tapti lineament and its surroundings of Central India (CI) is the second most important tectonic regime following the converging margin along Himalayas-Myanmar-Andaman of the Indian sub-continent, which attracted several geoscientists to assess its seismic hazard potential. Our study area, a part of CI, is bounded between latitudes 18°–26°N and longitudes 73°–83°E, representing a stable part of Peninsular India. Past damaging moderate magnitude earthquakes as well as continuing microseismicity in the area provided enough data for seismological study. Our estimates based on regional Gutenberg–Richter relationship showed lower b values (i.e., between 0.68 and 0.76) from the average for the study area. The Probabilistic Seismic Hazard Analysis carried out over the area with a radius of ~300 km encircling Bhopal yielded a conspicuous relationship between earthquake return period (T) and peak ground acceleration (PGA). Analyses of T and PGA shows that PGA value at bedrock varies from 0.08 to 0.15 g for 10 % (T = 475 years) and 2 % (T = 2,475 years) probabilities exceeding 50 years, respectively. We establish the empirical relationships $$ {\text{ZPA}}_{(T = 475)} = 0.1146\;[V_{\text{s}} (30)]^{ - 0.2924}, $$ and $$ {\text{ZPA}}_{(T = 2475)} = 0.2053\;[V_{\text{s}} (30)]^{ - 0.2426} $$ between zero period acceleration (ZPA) and shear wave velocity up to a depth of 30 m [V s (30)] for the two different return periods. These demonstrate that the ZPA values decrease with increasing shear wave velocity, suggesting a diagnostic indicator for designing the structures at a specific site of interest. The predictive designed response spectra generated at a site for periods up to 4.0 s at 10 and 2 % probability of exceedance of ground motion for 50 years can be used for designing duration dependent structures of variable vertical dimension. We infer that this concept of assimilating uniform hazard response spectra and predictive design at 10 and 2 % probability of exceedance in 50 years at 5 % damping at bedrocks of different categories may offer potential inputs for designing earthquake resistant structures of variable dimensions for the CI region under the National Earthquake Hazard Reduction Program for India.

Suárez, Gerardo; Jaramillo, Said; Bandy, W.

doi: 10.1007/s00024-013-0667-9pmid: N/A

So far, the direction and rate of relative motion between the Rivera and the North American plates (RIV-NAM) has been determined by the combination of two Euler poles: Rivera (RIV), with respect to Pacific (PAC), and PAC with respect to North America. Here, we estimate the relative motion of this plate pair (RIV-NAM) assuming that the horizontal projection of the direction of slip of the earthquakes occurring on the RIV-NAM boundaries reflect their relative plate motion. A catalog of earthquakes for which focal mechanisms are reported since 1976 is used in the analysis. Earthquakes were considered in the three segments of the RIV-NAM plate boundary: the subduction zone of the Rivera plate beneath the Jalisco block, the Tres Marias Escarpment and the events associated with the Tamayo Fracture Zone. The best fitting Euler pole is determined using a grid search of 64 potential poles. The slip direction predicted for each grid point is compared to the slip direction of the focal mechanisms of the earthquakes on the plate boundary. The best fitting Euler pole, determined in a root mean square sense (RMS), is located at 21.8°N, 107.6°W. A rate of rotation of 5.3°/year is estimated assuming the seismic earthquake cycle of the 1932 and 1995 great earthquakes represents a lower bound of the rate of plate motion in the subduction zone. The best fitting Euler pole shows that the subduction of the Rivera plate takes place in a direction perpendicular to the trench with a relative velocity of 4.3 cm/year, offshore Manzanillo. The rate of relative motion RIV-NAM decreases from SE to NW. North of approximately 21°N, the subduction of the Rivera plate becomes oblique to the trench and the relative velocity between the two plates decreases to an average of 1.9 cm/year. This slow rate of convergence may explain the rapid decrease of seismicity in the trench and the apparent absence of large earthquakes in this region. In the Tres Marias Escarpment, our best-fitting pole suggests that subduction stops, giving way to high-angle reverse faulting perpendicular to the Tres Marias Escarpment, in agreement with the reverse faulting earthquakes occurring here. To the north of 22.5°N, the slip predicted by the best-fitting pole suggests right-lateral faulting in a direction parallel to the Tamayo Fracture Zone, at a very low velocity (0.5–1.0 cm/year). The best fitting Euler pole determined here lies very close to the RIV-NAM plate boundary in the vicinity of the Tamayo Fracture Zone. This location of our best fitting Euler pole explains the low relative plate velocity, the relatively low level of seismic activity and the presence of a broad zone of deformation that accommodates the RIV-NAM motion.

Kolyukhin, Dmitry; Torabi, Anita

doi: 10.1007/s00024-013-0644-3pmid: N/A

This paper is devoted to statistical analysis of faults’ attributes. The distributions of lengths, widths of damage zones, displacements and thicknesses of fault cores are studied. Truncated power-law (TPL) is considered in comparison with commonly used simple power-law (PL) (or Pareto) distribution. The maximal likelihood and the confidence interval of the exponent for both PL and TPL are estimated by appropriate statistical methods. The Kolmogorov–Smirnov (KS) test and the likelihood ratio test with alternative non-nested hypothesis for exponential distribution are used to verify the statistical approximation. Furthermore, the advantage of TPL is proved by Bayesian information criterion. Our results suggest that a TPL is more suitable for describing fault attributes, and that its condition is satisfied for a wide range of fault scales. We propose that using truncated power laws in general might eliminate or relax the bias related to sampling strategy and the resolution of measurements (such as censoring, truncation, and cut effect) and; therefore, the most reliable range of data can be considered for the statistical approximation of fault attributes.

Mammo, Tilahun

doi: 10.1007/s00024-013-0663-0pmid: N/A

The Oligocene Afar mantle plume resulted in the eruption of a large volume of basaltic magma, including major sequences of rhyolitic ignimbrites, in a short span of time across Ethiopia. In order to assess the impact of these magmatic processes on the crust and to investigate the general crustal configuration beneath the Ethiopian plateau, northern part of the Main Ethiopian Rift and the Afar depression, analysis and modeling of the gravity field have been conducted. The Bouguer gravity map is dominated by long-wavelength anomalies that primarily arise from the isostatic compensation of the topography. Consequently, anomalies within the crust/upper mantle are masked and quantitative interpretation becomes difficult. The long-wavelength anomalies are approximated using admittance technique and subsequently removed from the Bouguer anomalies to obtain the residual isostatic anomalies. The residual map contains both short- and intermediate-wavelength anomalies related to geologic and tectonic features. The long-wavelength regional isostatic field is used to map the crust-mantle interface and the results are in good agreement with those determined by other geophysical methods. Seismic constrained gravity inversion was performed on the isostatic residual field and series of three-dimensional models have been constructed for the structures of the crust and upper mantle beneath the uplifted and rifted flood basalt province of northern Ethiopia. The inversion results have shown that the NW plateau has thick crust that rests on normal lithospheric mantle. Afar, On the other hand, is marked by thin stretched crust resting on a low-density upper mantle indicating a hotter thermal regime and partial melt. No lithospheric mantle is observed beneath Afar. The models further indicate the presence of an extensive sub-crustal thick (~12 km on average) and high-density (~3.06 gm/cc) mafic accreted igneous layer of fractionated cumulate (magmatic underplating) beneath the NW plateau. The study suggests that the underplate was fundamental to the accretion process and may have played a role in compensating most of the plateau uplift and in localizing stresses.