Stress Rotation in the Kachchh Rift Zone, Gujarat, India

Double difference relocations of the 1402 Kachchh events (2001–2006) clearly delineate two fault zones viz. south-dipping North Wagad fault (NWF) and almost vertical Gedi fault (GF). The relocated focal depths delineate a marked variation of 4 and 7 km in the brittle-ductile transition depths beneath GF and NWF, respectively. The focal mechanism solutions of 464 aftershocks (using 8–12 first motions) show that the focal mechanisms ranged between pure reverse and pure strike-slip except for a few pure dip-slip solutions. The stress inversions performed for five rectangular zones across the Kachchh rift reveal both clockwise and anticlockwise rotation (7–32°) in the σ₁ orientation within the rupture zone, favoring a heterogeneous stress regime with an average N-S fault normal compression. This rotation may be attributed to the presence of crustal mafic intrusives (5–35 km depth) in the rupture zone of the 2001 Bhuj main shock. Results suggest a relatively homogeneous stress regime in the GF zone favoring strike-slip motion, with a fault normal N-S compression.     

Strong earthquake activity all over the world and strong-moderate earthquake activity within and near China(June 2008)

Illustration All the data in this catalog are chosen from the "Preliminary Seismological Report of Chinese Seismic Stations" (Its abbreviation is "Monthly Report").The catalog includes the events of     

Receiver Functions for Three-layer Media

We have extended the H-k stacking method of receiver functions applicable to a three-layer model, which is useful in studying detailed crustal structures. We have demonstrated its application with two sample sites in Taiwan, making use of travel times of converted phases from the direct P waves of teleseismic events as the P reflects and refracts at different discontinuities in the crust. This three-layer extension allows a closer examination on the crust, as well as the relationship between crustal layers and the associated V_p/V_s ratios. Data were processed using the multiple-taper correlation technique to obtain the radial receiver functions (RRFs). The relative time delays of the converted phases measured from the RRFs were used to estimate the depths of the crustal discontinuities. Results not only yield the depths to principal crustal layers, including the Moho, but also give the corresponding ranges of V_p/V_s ratio which are comparable with findings from other tomographic studies.     

Strong earthquake activity all over the world and strong-moderate earthquake activity within and near China (December 2007~January 2008)

Illustration All the data in this catalog are chosen from the ″Preliminary Seismological Report of Chi-nese Seismic Stations″ (Its abbreviation is ″Monthly Report″). The catalog includes the events of M≥4.7 in and near China and M≥6 all over the world. The ″Monthly Report″ is monthly compiled     

Body-wave Attenuation in the Region of Garda, Italy

We analyzed the spectral amplitude decay with hypocentral distance of P and S waves generated by 76 small magnitude earthquakes (M_L 0.9–3.8) located in the Garda region, Central-Eastern Alps, Italy. These events were recorded by 18 stations with velocity sensors, in a distance range between 8 and 120 km. We calculated nonparametric attenuation functions (NAF) and estimated the quality factor Q of both body waves at 17 different frequencies between 2 and 25 Hz. Assuming a homogeneous model we found that the Q frequency dependence of P and S can be approximated with the functions Q _P = 65 f ^0.9 and Q _S = 160 f ^0.6 , respectively. At 2 Hz the Q _S /Q _P ratio reaches the highest value of 2.8. At higher frequencies Q _S /Q _P varies between 0.7 and 1.7, suggesting that for this frequency band scattering may be an important attenuation mechanism in the region of Garda. To explore the variation of Q in depth, we estimated Q at short (r ≤ 30 km) and intermediate (35–90 km) distance paths. We found that in the shallow crust P waves attenuate more than S (1.3 < Q _S /Q _P < 2.5). Moreover, P waves traveling along paths in the lower crust (depths approximately greater than 30 km) attenuate more than S waves. To quantify the observed variability of Q in depth we considered a three-layer model and inverted the NAF to estimate Q in each layer. We found that in the crust Q increases with depth. However, in the upper mantle (~40–50 km depth) Q decreases and in particular the high frequency Q _S (f > 9 Hz) has values similar to those estimated for the shallow layer of the crust.     

Strong earthquake activity all over the world and strong-moderate earthquake activity within and near China (April 2008～May 2008)

Illustration All the data in this catalog are chosen from the "Preliminary Seismological Report of Chinese Seismic Stations" (Its abbreviation is "Monthly Report").The catalog includes the events of     

Strong earthquake activity all over the world and strong-moderate earthquake activity within and near China (October 2007～November 2007)

Illustration All the data in this catalog are chosen from the ″Preliminary Seismological Report of Chi-nese Seismic Stations″ (Its abbreviation is ″Monthly Report″). The catalog includes the events of M≥4.7 in and near China and M≥6 all over the world. The ″Monthly Report″ is monthly compiled     

Strong earthquake activity all over the world and strong-moderate earthquake activity within and near China (February 2008～March 2008)

Illustration All the data in this catalog are chosen from the ″Preliminary Seismological Report of Chi-nese Seismic Stations″ (Its abbreviation is ″Monthly Report″). The catalog includes the events of     

Epicentral distribution in 2007

For showing the epicentral distribution in and near China as well as all over the world,two epicentral maps for the earthquakes occurred last year are published annually in the 6-th issue     

Data-resolution Matrix and Model-resolution Matrix for Rayleigh-wave Inversion Using a Damped Least-squares Method

Inversion of multimode surface-wave data is of increasing interest in the near-surface geophysics community. For a given near-surface geophysical problem, it is essential to understand how well the data, calculated according to a layered-earth model, might match the observed data. A data-resolution matrix is a function of the data kernel (determined by a geophysical model and a priori information applied to the problem), not the data. A data-resolution matrix of high-frequency (≥2 Hz) Rayleigh-wave phase velocities, therefore, offers a quantitative tool for designing field surveys and predicting the match between calculated and observed data. We employed a data-resolution matrix to select data that would be well predicted and we find that there are advantages of incorporating higher modes in inversion. The resulting discussion using the data-resolution matrix provides insight into the process of inverting Rayleigh-wave phase velocities with higher-mode data to estimate S-wave velocity structure. Discussion also suggested that each near-surface geophysical target can only be resolved using Rayleigh-wave phase velocities within specific frequency ranges, and higher-mode data are normally more accurately predicted than fundamental-mode data because of restrictions on the data kernel for the inversion system. We used synthetic and real-world examples to demonstrate that selected data with the data-resolution matrix can provide better inversion results and to explain with the data-resolution matrix why incorporating higher-mode data in inversion can provide better results. We also calculated model-resolution matrices in these examples to show the potential of increasing model resolution with selected surface-wave data.     

Compaction of a Rock Fracture Moderated by Competing Roles of Stress Corrosion and Pressure Solution

Unusually rapid closure of stressed fractures, observed in the initial stages of loading and at low temperatures, is examined using models for subcritical crack growth and pressure solution. The model for stress corrosion examines tensile stress concentrations induced at the Hertzian contact of propping fracture asperities, and mediates fracture growth according to a kinetic rate law. Conversely, pressure solution is described by the rate-limiting process of dissolution, resulting from the elevated stresses realized at the propping asperity contact. Both models are capable of following the observed compaction of fractures in novaculite. However, closure rates predicted for stress corrosion cracking are orders of magnitudes faster than those predicted for pressure dissolution. For consistent kinetic parameters, predictions from stress corrosion better replicate experimental observations, especially in the short-term and at low temperature when mechanical effects are anticipated to dominate. Rates and magnitudes of both stress corrosion and pressure solution are dependent on stresses exerted over propping asperities. Rates of closure due to stress corrosion cracking are shown to be always higher than for pressure solution, except where stress corrosion ceases as contact areas grow, and local stresses drop below an activation threshold. A simple rate law is apparent for the progress of fracture closure, defined in terms of a constant and an exponent applied to the test duration. For current experimental observations, this rate law is shown to replicate early progress data, and shows promise to define the evolution of transport properties of fractures over extended durations.     

First Order Seismic Microzonation of Haldia, Bengal Basin (India) Using a GIS Platform

The seismic microzonation of the Bengal Basin, Haldia region, India is carried out using the Analytical Hierarchy Process (AHP) on the Geographic Information System (GIS). Three themes are used for the seismic microzonation, namely Peak Ground Acceleration (PGA), predominant frequency and elevation map. An analysis of the maximum magnitude (m _max) and the b value is carried out after preparing the earthquake catalogue from various sources. On the basis of the tectonic set up and seismicity of the region, five seismic zones are delineated which can be a threat to Haldia. They are broadly classified as Zone 1: Arakan-Yoma Zone (AYZ), Zone 2: Himalayan Zone (HZ), Zone 3: Shillong Plateau Zone (SPZ), Zone 4: Bay of Bengal Zone (BBZ) and Zone 5: Shield Zone (SZ). The m _max for Zones 1, 2, 3, 4 and 5 are 8.30 ± 0.51, 9.09 ± 0.58, 9.20 ± 0.51, 6.62 ± 0.43 and 6.61 ± 0.43, respectively. The PGA value is computed for Haldia following the attenuation relationship taking the m _max of each source zone. The expected PGA at Haldia varies from 0.09–0.19 g. The predominant frequency of Haldia is also calculated using the H/V ratio with a frequency ranging from 0.1–3.0 Hz. The elevation map of Haldia is also generated using the Shuttle Radar Topography Mission (STRM) data. A first-order seismic microzonation map of Haldia is prepared in which four zones of hazard have been broadly classified for Haldia as very high seismic hazard zone, high seismic hazard zone, moderate seismic hazard zone and less seismic hazard zone. The very high seismic hazard zone is observed along the southern part of Haldia where there are major industrial and port facilities. The PGA for the four hazard zones are: 0.09–0.13 g for low hazard zone, > 0.13–0.15 g for moderate hazard zone, > 0.15–0.16 g for high hazard zone and > 0.16–0.19 g for very high hazard zone.