Modeling the strong ground motion and rupture characteristics of the March 31, 2006, Darb-e-Astane earthquake,Iran, using a hybrid of near-field SH-wave and empirical Green’s function method

We analyze the strong motion accelerograms of the moderate (M _w = 6.1), March 31, 2006, Darb-e-Astane earthquake of western Iran and also those of one of its prominently recorded, large (M _w = 5.1) foreshock and (M _w = 4.9) aftershock. (1) Using derived SH-wave spectral data, we first objectively estimate the parameters W _o\mathit{\Omega} _{\rm o} (long period spectral level), f _c (corner frequency) and Q(f) (frequency dependent, average shear wave quality factor), appropriate for the best-fit Brune ω ^− 2 spectrum of each of these three events. We then perform a non-linear least square analysis of the SH-wave spectral data to provide approximate near-field estimates of the strike, dip, and rake of the causative faults and also the seismic moment, moment magnitude, source size, and average stress drop of these three events. (2) In the next step, we use these approximate values and an empirical Green’s function approach, in an iterative manner, to optimally model the strong ground motion and rupture characteristics of the main event in terms of peak ground acceleration/velocity/displacement and duration of ground shaking and thereby provide improved, more reliable estimates of the causative fault parameters of the main event and its asperities. Our near-field estimates for both the main moderate event and the two smaller events are in good conformity with the corresponding far-field estimates reported by other studies.     

Lapse time and frequency-dependent attenuation characteristics of coda waves in the Northwestern Himalayas

We analyzed the local earthquakes waveform recorded on a broadband seismic network in the northwestern Himalayan Region to compute lapse time and frequency dependence of coda Q (Q _c). The observed Q _c values increase with increasing lapse time at all frequency bands. The increase in Q _c values with lapse time is attributed to an increase in Q _c with depth. This implies that attenuation decreases with increasing depth. The approximate radius of medium contributing to coda generation varies from 55 to 130 km. By comparing the Q _c values with those from other regions of the world, we find that they are similar to those obtained from tectonically active regions. The estimated Q _c values show a frequency-dependent relationship, Q _c = Q ₀ f ⁿ , where Q ₀ is Q _c at 1 Hz and n represents degree of frequency dependence. They represent the level of heterogeneity and tectonic activity in an area. Our results show that northwest Himalayas are highly heterogeneous and tectonically very active. Q ₀ increases from 113 ± 7 to 243 ± 10 and n decreases from 1.01 ± 0.05 to 0.85 ± 0.03 when lapse time increases from 30 to 70 s. As larger time window sees the effect of deeper part of the Earth, it is concluded that Q ₀ increases and n decreases with increasing depth; i.e., heterogeneity decreases with depth in the study area.     

A steady state of the disc dynamo

Abstract

We discuss the steady states of the αω-dynamo in a thin disc which arise due to α-quenching. Two asymptotic regimes are considered, one for the dynamo numberD near the generation thresholdD ₀, and the other for |D| ? 1. Asymptotic solutions for |D—D ₀| ? |D ₀| have a rather universal character provided only that the bifurcation is supercritical. For |D| ? 1 the asymptotic solution crucially depends on whether or not the mean helicity α, as a function ofB, has a positive root (hereB is the mean magnetic field). When such a root exists, the field value in the major portion of the disc is O(l), while near the disc surface thin boundary layers appear where the field rapidly decreases to zero (if the disc is surrounded by vacuum). Otherwise, when α = O(|B|^?s) for |B| → ∞, we demonstrate that |B| = O(|D|^1/s ) and the solution is free of boundary layers. The results obtained here admit direct comparison with observations of magnetic fields in spiral galaxies, so that an appropriate model of nonlinear galactic dynamos hopefully could be specified.     

Coda wave attenuation in the Parecis Basin,Amazon Craton,Brazil: sensitivity to basement depth

Small local earthquakes from two aftershock sequences in Porto dos Gaúchos, Amazon craton—Brazil, were used to estimate the coda wave attenuation in the frequency band of 1 to 24 Hz. The time-domain coda-decay method of a single backscattering model is employed to estimate frequency dependence of the quality factor (Q _c) of coda waves modeled using Q_c = Q₀ f^hQ_{\rm c} =Q_{\rm 0} f^\eta , where Q ₀ is the coda quality factor at frequency of 1 Hz and η is the frequency parameter. We also used the independent frequency model approach (Morozov, Geophys J Int, 175:239–252, 2008), based in the temporal attenuation coefficient, χ(f) instead of Q(f), given by the equation c(f)=g+\fracpfQ_e \chi (f)\!=\!\gamma \!+\!\frac{\pi f}{Q_{\rm e} }, for the calculation of the geometrical attenuation (γ) and effective attenuation (Q_e^-1 )(Q_{\rm e}^{-1} ). Q _c values have been computed at central frequencies (and band) of 1.5 (1–2), 3.0 (2–4), 6.0 (4–8), 9.0 (6–12), 12 (8–16), and 18 (12–24) Hz for five different datasets selected according to the geotectonic environment as well as the ability to sample shallow or deeper structures, particularly the sediments of the Parecis basin and the crystalline basement of the Amazon craton. For the Parecis basin Q_c = (98±12)f^(1.14±0.08)Q_{\rm c} =(98\pm 12)f^{(1.14\pm 0.08)}, for the surrounding shield Q_c = (167±46)f^(1.03±0.04)Q_{\rm c} =(167\pm 46)f^{(1.03\pm 0.04)}, and for the whole region of Porto dos Gaúchos Q_c = (99±19)f^(1.17±0.02)Q_{\rm c} =(99\pm 19)f^{(1.17\pm 0.02)}. Using the independent frequency model, we found: for the cratonic zone, γ = 0.014 s^− 1, Q_e^-1 = 0.0001Q_{\rm e}^{-1} =0.0001, ν ≈ 1.12; for the basin zone with sediments of ~500 m, γ = 0.031 s^− 1, Q_e^-1 = 0.0003Q_{\rm e}^{-1} =0.0003, ν ≈ 1.27; and for the Parecis basin with sediments of ~1,000 m, γ = 0.047 s^− 1, Q_e^-1 = 0.0005Q_{\rm e}^{-1} =0.0005, ν ≈ 1.42. Analysis of the attenuation factor (Q _c) for different values of the geometrical spreading parameter (ν) indicated that an increase of ν generally causes an increase in Q _c, both in the basin as well as in the craton. But the differences in the attenuation between different geological environments are maintained for different models of geometrical spreading. It was shown that the energy of coda waves is attenuated more strongly in the sediments, Q_c = (78±23)f^(1.17±0.14)Q_{\rm c} =(78\pm 23)f^{(1.17\pm 0.14)} (in the deepest part of the basin), than in the basement, Q_c = (167±46)f^(1.03±0.04)Q_{\rm c} =(167\pm 46)f^{(1.03\pm 0.04)} (in the craton). Thus, the coda wave analysis can contribute to studies of geological structures in the upper crust, as the average coda quality factor is dependent on the thickness of sedimentary layer.     

Monitoring the West Bohemian earthquake swarm in 2008/2009 by a temporary small-aperture seismic array

The most recent intense earthquake swarm in West Bohemia lasted from 6 October 2008 to January 2009. Starting 12 days after the onset, the University of Potsdam monitored the swarm by a temporary small-aperture seismic array at 10 km epicentral distance. The purpose of the installation was a complete monitoring of the swarm including micro-earthquakes (M _L < 0). We identify earthquakes using a conventional short-term average/long-term average trigger combined with sliding-window frequency-wavenumber and polarisation analyses. The resulting earthquake catalogue consists of 14,530 earthquakes between 19 October 2008 and 18 March 2009 with magnitudes in the range of − 1.2 ≤ M _L ≤ 2.7. The small-aperture seismic array substantially lowers the detection threshold to about M _c = − 0.4, when compared to the regional networks operating in West Bohemia (M _c > 0.0). In the course of this work, the main temporal features (frequency–magnitude distribution, propagation of back azimuth and horizontal slowness, occurrence rate of aftershock sequences and interevent-time distribution) of the recent 2008/2009 earthquake swarm are presented and discussed. Temporal changes of the coefficient of variation (based on interevent times) suggest that the swarm earthquake activity of the 2008/2009 swarm terminates by 12 January 2009. During the main phase in our studied swarm period after 19 October, the b value of the Gutenberg–Richter relation decreases from 1.2 to 0.8. This trend is also reflected in the power-law behavior of the seismic moment release. The corresponding total seismic moment release of 1.02×10¹⁷ Nm is equivalent to M _L,max = 5.4.     

Characteristics of coda wave attenuation in Yunnan area

The characteristic of seismic coda wave attenuation in Yunnan area in 7 frequency-bands range from 1 Hz to 20 Hz was estimated by using the local earthquake's waveform data recorded from 22 Yunnan digital seismic stations.Coda attenuation Q-c1 of each station was firstly calculated by single scattering method. Then, mean free path Le and seismic albedo Bo of each station were calculated, and scattering attenuation Q-1s and intrinsic attenuation Q-1i were separated from total attenuation Q-1t by multiple lapse time window analysis based on the multiple scattering model in uniform random isotropic scattering medium. The attenuating characteristics in Yunnan show that most value of Le are in 10～30 km, with maximal within 2～6 Hz;Bo are about 0.5 at 1～2 Hz, but less than 0.5at other frequency-bands, which means Q-1i is comparable with Q-1s at 1～2 Hz, and after 1～2 Hz, Q-1i is greater than Q-1s and dominates the attenuation process. Q-1c is close to Q-1i at other frequency bands except 1～2 Hz.Results show that Q-1 especially Qs-1 varies spatially, Q-1 in eastern Yunnan zone is a bit higher than in northwestern Yunnan zone;northwestern Yunnan zone higher than southwestern Yunnan zone. Comparing with other results in global, Qs-1 in Yunnan is lower than the global average value among these results, Q-1i is higher than the global average value, and Q-1t lies the middle among these results.     

Numerical simulations of transport of non-ergodic solute plumes in heterogeneous aquifers

Transport of non-ergodic solute plumes by steady-state groundwater flow with a uniform mean velocity, μ, were simulated with Monte Carlo approach in a two-dimensional heterogeneous and statistically isotropic aquifer whose transmissivity, T, is log-normally distributed with an exponential covariance. The ensemble averages of the second spatial moments of the plume about its center of mass, <S _{i i} (t)>, and the plume centroid covariance, R _{i i} (t) (i=1,2), were simulated for the variance of Y=log T, σ _Y ²=0.1, 0.5 and 1.0 and line sources normal or parallel to μ of three dimensionless lengths, 1, 5, and 10. For σ _Y ²=0.1, all simulated <S _{i i} (t)>−S _{i i} (0) and R _{i i} (t) agree well with the first-order theoretical values, where S _{i i} (0) are the initial values of S _{i i} (t). For σ _Y ²=0.5 and 1.0 and the line sources normal to μ, the simulated longitudinal moments, <S ₁₁(t)>−S ₁₁(0) and R ₁₁(t), agree well with the first-order theoretical results but the simulated transverse moments <S ₂₂(t)>−S ₂₂(0) and R ₂₂(t) are significantly larger than the first-order values. For the same two larger values of σ _Y ² but the line sources parallel to μ, the simulated <S ₁₁(t)>−S ₁₁(0) are larger than but the simulated R ₁₁ are smaller than the first-order values, and both simulated <S ₂₂(t)>−S ₂₂(0) and R ₂₂(t) stay larger than the first-order values. For a fixed value of σ _Y ², the summations of <S _{i i} (t)>−S _{i i} (0) and R _{i i} , i.e., X _{i i} (i=1,2), remain almost the same no matter what kind of source simulated. The simulated X ₁₁ are in good agreement with the first-order theory but the simulated X ₂₂ are significantly larger than the first-order values. The simulated X ₂₂, however, are in excellent agreement with a previous modeling result and both of them are very close to the values derived using Corrsin's conjecture. It is found that the transverse moments may be significantly underestimated if less accurate hydraulic head solutions are used and that the decreasing of <S ₂₂(t)>−S ₂₂(0) with time or a negative effective dispersivity, defined as , may happen in the case of a line source parallel to μ where σ _Y ² is small.     

Numerical simulations of transport of non-ergodic solute plumes in heterogeneous aquifers

Transport of non-ergodic solute plumes by steady-state groundwater flow with a uniform mean velocity, μ, were simulated with Monte Carlo approach in a two-dimensional heterogeneous and statistically isotropic aquifer whose transmissivity, T, is log-normally distributed with an exponential covariance. The ensemble averages of the second spatial moments of the plume about its center of mass, <S _{i i} (t)>, and the plume centroid covariance, R _{i i} (t) (i=1,2), were simulated for the variance of Y=log T, σ _Y ²=0.1, 0.5 and 1.0 and line sources normal or parallel to μ of three dimensionless lengths, 1, 5, and 10. For σ _Y ²=0.1, all simulated <S _{i i} (t)>−S _{i i} (0) and R _{i i} (t) agree well with the first-order theoretical values, where S _{i i} (0) are the initial values of S _{i i} (t). For σ _Y ²=0.5 and 1.0 and the line sources normal to μ, the simulated longitudinal moments, <S ₁₁(t)>−S ₁₁(0) and R ₁₁(t), agree well with the first-order theoretical results but the simulated transverse moments <S ₂₂(t)>−S ₂₂(0) and R ₂₂(t) are significantly larger than the first-order values. For the same two larger values of σ _Y ² but the line sources parallel to μ, the simulated <S ₁₁(t)>−S ₁₁(0) are larger than but the simulated R ₁₁ are smaller than the first-order values, and both simulated <S ₂₂(t)>−S ₂₂(0) and R ₂₂(t) stay larger than the first-order values. For a fixed value of σ _Y ², the summations of <S _{i i} (t)>−S _{i i} (0) and R _{i i} , i.e., X _{i i} (i=1,2), remain almost the same no matter what kind of source simulated. The simulated X ₁₁ are in good agreement with the first-order theory but the simulated X ₂₂ are significantly larger than the first-order values. The simulated X ₂₂, however, are in excellent agreement with a previous modeling result and both of them are very close to the values derived using Corrsin's conjecture. It is found that the transverse moments may be significantly underestimated if less accurate hydraulic head solutions are used and that the decreasing of <S ₂₂(t)>−S ₂₂(0) with time or a negative effective dispersivity, defined as , may happen in the case of a line source parallel to μ where σ _Y ² is small.     

Shallow water tidal currents in close proximity to the seafloor and boundary-induced turbulence

Velocity measurements with vertical resolution 0.02 m were conducted in the lowest 0.5 m of the water column using acoustic Doppler current profiler (ADCP) at a test site in the western part of the East China Sea. The friction velocity u _* and the turbulent kinetic energy dissipation rate ε _wl(ζ) profiles were calculated using log-layer fits; ζ is the height above the bottom. During a semidiurnal tidal cycle, u _* was found to vary in the range (1–7) × 10⁻³ m/s. The law-of-the-wall dissipation profiles ε _wl(ζ) were consistent with the dissipation profiles ε _mc(ζ) evaluated using independent microstructure measurements of small-scale shear, except in the presence of westward currents. It was hypothesized that an isolated bathymetric rise (25 m height at a 50-m seafloor) located to the east of the measurement site is responsible for the latter. Calculation of the depth integrated internal tide generating body force in the region showed that the flanks of the rise are hotspots of internal wave energy that may locally produce a significant turbulent zone while emitting tidal and shorter nonlinear internal waves. This distant topographic source of turbulence may enhance the microstructure-based dissipation levels ε _mc(ζ) in the bottom boundary layer (BBL) beyond the dissipation ε _wl(ζ) associated with purely locally generated turbulence by skin currents. Semi-empirical estimates for dissipation at a distance from the bathymetric rise agree well with the BBL values of ε _mc measured 15 km upslope.     

Model of vibrational level populations of Herzberg states of oxygen molecules at heights of the lower thermosphere and mesosphere

The paper presents a model of the kinetics of electronically excited O₂(c¹Σ_u⁻,v), O₂(A′³Δ_u,v), O₂(A³Σ_u⁺,v) molecules at heights of the lower thermosphere and mesosphere with allowance for electronic excitation transfer processes during molecular collisions. The model is used to calculate the relative O₂(A³Σ_u⁺,v) and O₂(A′³Δ_u,v) populations at heights of 80–110 km. The calculated populations are compared with the available literature results on experimental estimates, and good agreement is obtained. It is shown how the increase in the quenching rates of the considered states by oxygen atoms affects the calculation results.     

Structure refinement of astrophyllite

The crystal structure of astrophyllite K₂Na(Fe, Mn, Mg,□)₇[Ti₂(Si₄O₁₂)₂|O₃](OH, F)₄ has been refined. The dimensions of the triclinic unit cell are: a = 0.5359(2) nm,b = 1.1614(4) nm, c = 1.1861(4) nm, α= 113.16(2)°, β= 103.04(2)°,γ= 94.56(2)°,V = 0.6495(5) nm³, Z= 1, space group P1, R=0.057 for 5308 reflections |F_o|>3σ|F_o|. According to structural and compositional differences the monoclinic astrophyllite K₂NaNa(Fe, Mn)₄Mg₂Ti2[Si₄O₁₂]₂(OH)₄(OH, F)₂ and astrophyllite should be considered as two different mineral species. Astrophyllite, monoclinic astrophyllite, bafertisite and lamprophyllite contain heteropolyhedral sheets which topologically are related with Si, O sheets of mica where one or several SiO₄ tetrahedra are replaced by TiO _n polyhedra. Therefore this heterophyllotitanosilicate series represents a kind of functional substitution in inorganic crystals.     

Lateral variations of coda Q and attenuation of seismic waves in Southwest Anatolia

The seismic quality factor (Q _c) and the attenuation coefficient (δ) in the earth’s crust in southwest (SW) Anatolia are estimated by using the coda wave method based on the decrease of coda wave amplitude by time on the seismogram. The quality factor Q _o, the value of Q _c at 1 Hz, and its frequency dependency η are determined from this method depending on the attenuation properties of scattered coda waves. δ is determined from the observations of amplitude variations of seismic waves. In applying the coda wave method, firstly, a type curve representing the average pattern of the individual coda decay curves for 0.75, 1.5, 3.0, 6.0, 12.0, and 24.0 Hz values was estimated. Secondly, lateral variation of coda Q and the attenuation coefficients for three main tectonic patterns are estimated. The shape of the type curve is controlled by the scattering and attenuation in the crustal volume sampled by the coda waves. The Q _o and η values vary from 30 to 180 and from 0.55 to 1.25, respectively for SW Anatolia. In SW Anatolia, coda Q–f relation is described by and δ = 0.008 km⁻¹. These results are expected to help in understanding the degree of tectonic complexity of the crust in SW Anatolia.     

Dependence of geomagnetic activity during magnetic storms on the solar wind parameters for different types of streams

The dependence of the maximal values of the |Dst| and AE geomagnetic indices observed during magnetic storms on the value of the interplanetary electric field (E _y ) was studied based on the catalog of the large-scale solar wind types created using the OMNI database for 1976–2000 [Yermolaev et al., 2009]. An analysis was performed for eight categories of magnetic storms caused by different types of solar wind streams: corotating interaction regions (CIR, 86 storms); magnetic clouds (MC, 43); Sheath before MCs (Sh_MC, 8); Ejecta (95); Sheath (Sh_E, 56); all ICME events (MC + Ejecta, 138); all compression regions Sheaths before MCs and Ejecta (Sh_MC + Sh_E, 64); and an indeterminate type of storm (IND, 75). It was shown that the |Dst| index value increases with increasing electric field E _y for all eight types of streams. When electric fields are strong (E _y > 11 mV m⁻¹), the |Dst| index value becomes saturated within magnetic clouds MCs and possibly within all ICMEs (MC + Ejecta). The AE index value during magnetic storms is independent of the electric field value E _y for almost all streams except magnetic clouds MCs and possibly the compressed (Sheath) region before them (Sh_MC). The AE index linearly increases within MC at small values of the electric field (E _y < 11 mV m⁻¹) and decrease when these fields are strong (E _y > 11 mV m⁻¹). Since the dynamic pressure (Pd) and IMF fluctuations (σB) correlate with the E _y value in all solar wind types, both geomagnetic indices (|Dst| and AE) do not show an additional dependence on Pd and IMF δB. The nonlinear relationship between the intensities of the |Dst| and AE indices and the electric field E _y component, observed within MCs and possibly all ICMEs during strong electric fields E _y , agrees with modeling the magnetospheric-ionospheric current system of zone 1 under the conditions of the polar cap potential saturation.     

Turbulent transport of magnetic fields. V. Distribution of magnetic energy in a simple α2-dynamo

Abstract

In this paper we analyse the stationary mean energy density tensor T_ij = B_iB_j for the x ²-sphere. This model is one of the simplest possible turbulent dynamos, originally due to Krause and Steenbeck (1967): a conducting sphere of radius R with homogeneous, isotropic and stationary turbulent convection, no differential rotation and negligible resistivity. The stationary solution of the (linear) equation for T_ij is found analytically. Only T_rr , T _θθ and T _φφ are unequal to zero, and we present their dependence on the radial distance r.

The stationary solution depends on two coefficients describing the turbulent state: the diffusion coefficient β≈?u²?τ_c/3 and the vorticity coefficient γ ≈ ?|?×u|²?τ_c/3 where u(r, t) is the turbulent velocity and _c its correlation time. But the solution is independent of the dynamo coefficient α≈??u·?×u?τ_c/3 although α does occur in the equation for T_ij . This result confirms earlier conclusions that helicity is not required for magnetic field generation. In the stationary state, magnetic energy is generated by the vorticity and transported to the boundary, where it escapes at the same rate. The solution presented contains one free parameter that is connected with the distribution of B over spatial scales at the boundary, about which T_ij gives no information. We regard this investigation as a first step towards the analysis of more complicated, solar-type dynamos.     

Vesiculation path of ascending magma in the 1983 and the 2000 eruptions of Miyakejima volcano, Japan

The vesiculation of magma during the 1983 eruption of Miyakejima Volcano, Japan, is discussed based on systematic investigations of water content, vesicularity, and bubble size distribution for the products. The eruption is characterized by simultaneous lava effusion and explosive sub-plinian (‘dry’) eruptions with phreatomagmatic (‘wet’) explosions. The magmas are homogeneous in composition (basaltic andesite) and in initial water content (H₂O = 3.9±0.9 wt%), and residual groundmass water contents for all eruption styles are low (H₂O <0.4 wt%) suggestive of extensive dehydration of magma. For the scoria erupted during simultaneous ‘dry’ and ‘wet’ explosive eruptions, inverse correlation was observed between vesicularity and residual water content. This relation can be explained by equilibrium exsolution and expansion of ca. 0.3 wt% H₂O at shallow level with different times of quenching, and suggests that each scoria with different vesicularity, which was quenched at a different time, provides a snapshot of the vesiculation process near the point of fragmentation. The bubble size distribution (BSD) varies systematically with vesicularity, and total bubble number density reaches a maximum value at vesicularity Φ ∼ 0.5. At Φ  ∼ 0.5, a large number of bubbles are connected with each other, and the average thickness of bubble walls reaches the minimum value below which they would rupture. These facts suggest that vesiculation advanced by nucleation and growth of bubbles when Φ < 0.5, and then by expansion of large bubbles with coalescence of small ones for Φ > 0.5, when bubble connection becomes effective. Low vesicularity and low residual water content of lava and spatter (Φ  < 0.1, H₂O  < 0.1 wt%), and systematic decrease in bubble number density from scoria through spatter to lava with decrease in vesicularity suggest that effusive eruption is a consequence of complete degassing by bubble coalescence and separation from magma at shallow levels when magma ascent rate is slow.

Physical properties of the 1980 Mount St. Helens cryptodome magma

 Physical properties of cryptodome and remelted samples of the Mount St. Helens grey dacite have been measured in the laboratory. The viscosity of cryptodome dacite measured by parallel–plate viscometry ranges from 10.82 to 9.94 log₁₀ η (Pa s) (T=900–982  °C), and shrinkage effects were dilatometrically observed at T>900  °C. The viscosity of remelted dacite samples measured by the micropenetration method is 10.60–9.25 log₁₀ η (Pa s) (T=736–802  °C) and viscosities measured by rotational viscometry are 3.22–1.66 log₁₀ η (Pa s) (T=1298–1594  °C). Comparison of the measured viscosity of cryptodome dacitic samples with the calculated viscosity of corresponding water-bearing melt demonstrates significant deviations between measured and calculated values. This difference reflects a combination of the effect of crystals and vesicles on the viscosity of dacite as well as the insufficient experimental basis for the calculation of crystal-bearing vesicular melt viscosities at low temperature. Assuming that the cryptodome magma of the 18 May 1980 Mount St. Helens eruption was residing at 900  °C with a phenocryst content of 30 vol.%, a vesicularity of 36 vol.% and a bulk water content of 0.6 wt.%, we estimate the magma viscosity to be 10^10.8 Pa s. Received: 25 August 1996 / Accepted: 19 July 1997     

Coda Q in the Kachchh Basin, Western India Using Aftershocks of the Bhuj Earthquake of January 26, 2001

Q_C-estimates of Kachchh Basin in western India have been obtained in a high frequency range from 1.5 to 24.0 Hz using the aftershock data of Bhuj earthquake of January 26, 2001 recorded within an epicentral distance of 80 km. The decay of coda waves of 30 sec window from 186 seismograms has been analysed in four lapse time windows, adopting the single backscattering model. The study shows that Q_c is a function of frequency and increases as frequency increases. The frequency dependent Q_c relations obtained for four lapse-time windows are: Q_c=82 f^1.17 (20–50 sec), Q_c=106 f^1.11 (30–60 sec), Q_c=126f^1.03 (40–70 sec) and Q_c=122f^1.02 (50–80 sec). These empirical relations represent the average attenuation properties of a zone covering the surface area of about 11,000, 20,000, 28,000 and 38,000 square km and a depth extent of about 60, 80, 95, 110 km, respectively. With increasing window length, the degree of frequency dependence, n, decreases marginally from 1.17 to 1.02, whereas Q₀ increases significantly from 82 to 122. At lower frequencies up to 6 Hz, Q_c⁻¹ of Kachchh Basin is in agreement with other regions of the world, whereas at higher frequencies from 12 to 24 Hz it is found to be low.     

CO<Subscript>2</Subscript> and He degassing at El Chichón volcano,Chiapas, Mexico: gas flux,origin and relationship with local and regional tectonics

During 2007–2008, three CO₂ flux surveys were performed on El Chichón volcanic lake, Chiapas, Mexico, with an additional survey in April 2008 covering the entire crater floor (including the lake). The mean CO₂ flux calculated by sequential Gaussian simulation from the lake was 1,190 (March 2007), 730 (December 2007) and 1,134 g m⁻² day⁻¹ (April 2008) with total emission rates of 164 ± 9.5 (March 2007), 59 ± 2.5 (December 2007) and 109 ± 6.6 t day⁻¹ (April 2008). The mean CO₂ flux estimated from the entire crater floor area was 1,102 g m⁻² day⁻¹ for April 2008 with a total emission rate of 144 ± 5.9 t day⁻¹. Significant change in CO₂ flux was not detected during the period of survey, and the mapping of the CO₂ flux highlighted lineaments reflecting the main local and regional tectonic patterns. The ³He/⁴He ratio (as high as 8.1 R _A) for gases in the El Chichón crater is generally higher than those observed at the neighbouring Transmexican Volcanic Belt and the Central American Volcanic Arc. The CO₂/³He ratios for the high ³He/⁴He gases tend to have the MORB-like values (1.41 × 10⁹), and the CO₂/³He ratios for the lower ³He/⁴He gases fall within the range for the arc-type gases. The high ³He/⁴He ratios, the MORB-like CO₂/³He ratios for the high ³He/⁴He gases and high proportion of MORB-CO₂ (M = 25 ±15%) at El Chichón indicate a greater depth for the generation of magma when compared to typical arc volcanoes.     

Regional characteristics of coda attenuation after the lancang-Gengma earthquake

In the light of the single scattering model of coda originating from local earthquakes, and based on the aftershock coda registered respectively at the 4 short period stations installed near the foci shortly after theM7.6 Lancang andM7.2 Gengma earthquakes, this paper has tentatively calculated the rate of amplitude attenuation and theQ _c-value of the coda in the Lancang and Gengma areas using a newly-founded synthetic determination method. Result of the study shows the rate of coda amplitude attenuation demonstrates remarkable regional differences respectively in the southern and northern areas. The southern area presents a faster attenuation (Q _c=114), whereas the northern area shows a slower attenuation (Q _c=231). The paper also discusses the reasons causing such differences. Result of the study also suggests a fairly good linear relation between the coda source factorA _o(f) and the seismic moment and the magnitude. Using the earthquake scaling law, the following formulas can be derived: lgM ₀=lgA ₀(f)+17.6,M _D=0.67lgA ₀(f)+1.21 and logM ₀=1.5M _D+15.79. In addition, the rates of amplitude attenuationβ _s andβ _m are respectively calculated using the single scattering and multiple scattering models, and the ratioβ _s/β_m=1.20−1.50 is found for the results respectively from the two models. Finally, the mean free pathL of the S-wave scattering in the southern and northern areas are determined to be 54 km and 122 km respectively by the relations which can distinguish between the inherentQ _i and scatteringQ _s, testify to this areas having lowQ-values correspond to stronger scatterings. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, 71–82, 1992. This study is partly supported by the Seismological Science Foundation of the State Seismological Bureau of China, and the present English version of the paper is translated from its Chinese original by Wenyi Xia, Seismological Bureau of Yunnan Province.     

Retrospective forecasting test of a statistical physics model for earthquakes in Sichuan-Yunnan region

Pattern informatics (PI) model is one of the recently developed predictive models of earthquake phys- ics based on the statistical mechanics of complex systems. In this paper, retrospective forecast test of the PI model was conducted for the earthquakes in Sichuan-Yunnan region since 1988, exploring the possibility to apply this model to the estimation of time-dependent seismic hazard in continental China. Regional earthquake catalogue down to ML3.0 from 1970 to 2007 was used. The ‘target magnitude’ for the forecast test was MS5.5. Fifteen-year long ‘sliding time window’ was used in the PI calculation, with ‘anomaly training time window’ being 5 years and ‘forecast time window’ being 5 years, respectively. Receiver operating characteristic (ROC) test was conducted for the evaluation of the forecast result, showing that the PI forecast outperforms not only random guess but also the simple number counting approach based on the clustering hypothesis of earthquakes (the RI forecast). If the ‘forecast time window’ was shortened to 3 years and 1 year, respectively, the forecast capability of the PI model de- creased significantly, albeit outperformed random forecast. For the one year ‘forecast time window’, the PI result was almost comparable to the RI result, indicating that clustering properties play a more important role at this time scale.