Attenuation of P,S, and coda waves in Koyna region,India

The attenuation properties of the crust in the Koyna region of the Indian shield have been investigated using 164 seismograms from 37 local earthquakes that occurred in the region. The extended coda normalization method has been used to estimate the quality factors for P waves and S waves , and the single back-scattering model has been used to determine the quality factor for coda waves (Q _c). The earthquakes used in the present study have the focal depth in the range of 1–9 km, and the epicentral distance vary from 11 to 55 km. The values of and Q _c show a dependence on frequency in the Koyna region. The average frequency dependent relationships (Q = Q ₀ f ⁿ) estimated for the region are , and . The ratio is found to be greater than one for the frequency range considered here (1.5–18 Hz). This ratio, along with the frequency dependence of quality factors, indicates that scattering is an important factor contributing to the attenuation of body waves in the region. A comparison of Q _c and in the present study shows that for frequencies below 4 Hz and for the frequencies greater than 4 Hz. This may be due to the multiple scattering effect of the medium. The outcome of this study is expected to be useful for the estimation of source parameters and near-source simulation of earthquake ground motion, which in turn are required in the seismic hazard assessment of a region.     

Moment inequality and complete convergence of moving average processes under asymptotically linear negative quadrant dependence assumptions

Let {Y, Y _i , −∞ < i < ∞} be a doubly infinite sequence of identically distributed and asymptotically linear negative quadrant dependence random variables, {a _i , −∞ < i < ∞} an absolutely summable sequence of real numbers. We are inspired by Wang et al. (Econometric Theory 18:119–139, 2002) and Salvadori (Stoch Environ Res Risk Assess 17:116–140, 2003). And Salvadori (Stoch Environ Res Risk Assess 17:116–140, 2003) have obtained Linear combinations of order statistics to estimate the quantiles of generalized pareto and extreme values distributions. In this paper, we prove the complete convergence of under some suitable conditions. The results obtained improve and generalize the results of Li et al. (1992) and Zhang (1996). The results obtained extend those for negative associated sequences and ρ^*-mixing sequences. CIC Number O211, AMS (2000) Subject Classification 60F15, 60G50 Research supported by National Natural Science Foundation of China     

Estimation of coda wave attenuation in East Central Iran

The attenuation of coda waves, Q _c , has been estimated in Zarand, Jiroft, and Bam regions of east central Iran using a single back-scattering model of S-coda envelopes. For this purpose, the recordings of 97 earthquakes by three seismic networks and a local strong ground motion network have been used. In this research, the frequency-dependent Q _c values are estimated at central frequencies of 1.5, 3, 6, 8, 12, 16, and 24 Hz using different lapse time windows from 20 to 60 s. The frequency-dependent relationships obtained are for Zarand, for Jiroft, and for Bam region. From the strong ground motion data, we obtain the relation . The Q _c frequency-dependent relationship for the entire region of east central Iran from all data (both seismograms and accelerograms) is . The average Q _c values estimated and their frequency dependent relationships correlate well with a highly heterogeneous and highly tectonically active region. Results also show that the attenuation is higher in Bam region compared to Zarand and Jiroft regions.     

Airborne Gradiometry Error Analysis

Gravity gradiometry is one of the older methods of determining the Earth’s local gravitational field, but lies in the shadow of more conventional static and moving-base gravimeter-based systems. While the static torsion balance appears to have been relegated to the museum, support for the airborne and space-borne differential accelerometer (gradiometer) continues so as to overcome limitations in spatial resolution and accuracy inherent in ordinary moving-base gravimetry. One airborne system exists, building on 30 year old technology concepts, and new technologies (e.g., cold-atom interferometry) promise significant improvements. Concomitant advances are required to measure accurately the angular velocity and angular acceleration of the platform, which inseparably combine (in an absolute sense) with the Earth’s gravitational gradients. A numerical analysis of instrument errors, with simulated aircraft dynamics, shows that navigation-grade gyros are just sufficient to account for these effects in gradiometers with 1E/ sensitivity. More accurate instruments, with 0.1 E/ sensitivity, require commensurate sensitivity in the gyros, of the order of 0.01°/h/ = 1.5\times10^{−4 °} \ for typical survey aircraft dynamics. On the other hand, typical orientation errors in the platform, which are problematic for vector gravimetry, are much less of a concern in gradiometry. They couple to the gradient signals and affect only the very low frequencies of the total gradient error.     

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.     

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.     

Generous statistical tests

A common statistical problem is deciding which of two possible sources, A and B, of a contaminant is most likely the actual source. The situation considered here, based on an actual problem of polychlorinated biphenyl contamination discussed below, is one in which the data strongly supports the hypothesis that source A is responsible. The problem approach here is twofold: One, accurately estimating this extreme probability. Two, since the statistics involved will be used in a legal setting, estimating the extreme probability in such a way as to be as generous as is possible toward the defendant’s claim that the other site B could be responsible; thereby leaving little room for argument when this assertion is shown to be highly unlikely. The statistical testing for this problem is modeled by random variables {X _i } and the corresponding sample mean the problem considered is providing a bound ɛ for which for a given number a ₀. Under the hypothesis that the random variables {X _i } satisfy E(X _i ) ≤ μ, for some 0  < μ < 1, statistical tests are given, described as “generous”, because ɛ is maximized. The intent is to be able to reject the hypothesis that a ₀ is a value of the sample mean while eliminating any possible objections to the model distributions chosen for the {X _i } by choosing those distributions which maximize the value of ɛ for the test used.     

Ground-motion Attenuation Relation from Strong-motion Records of the 2001 Mw 7.7 Bhuj Earthquake Sequence (2001–2006), Gujarat, India

Predictive relations are developed for peak ground acceleration (PGA) from the engineering seismoscope (SRR) records of the 2001 M_w 7.7 Bhuj earthquake and 239 strong-motion records of 32 significant aftershocks of 3.1 ≤ M_w ≤ 5.6 at epicentral distances of 1 ≤ R ≤ 288 km. We have taken advantage of the recent increase in strong-motion data at close distances to derive new attenuation relation for peak horizontal acceleration in the Kachchh seismic zone, Gujarat. This new analysis uses the Joyner-Boore’s method for a magnitude-independent shape, based on geometrical spreading and anelastic attenuation, for the attenuation curve. The resulting attenuation equation is,

Calibration of an M L scale for South Africa using tectonic earthquake data recorded by the South African National Seismograph Network: 2006 to 2009

A relation to determine local magnitude (M _L) based on the original Richter definition is empirically derived from synthetic Wood–Anderson seismograms recorded by the South African National Seismograph Network. In total, 263 earthquakes in the distance range 10 to 1,000 km, representing 1,681 trace amplitudes measured in nanometers from synthesized Wood–Anderson records on the vertical channel were considered to derive an attenuation relation appropriate for South Africa through multiple regression analysis. Additionally, station corrections were determined for 26 stations during the regression analysis resulting in values ranging between ?0.31 and 0.50. The most appropriate M _L scale for South Africa from this study satisfies the equation: $$ {M_{{{\bf L}}}} = {\text{lo}}{{\text{g}}_{{10}}}(A) + 1.149\;{\text{lo}}{{\text{g}}_{{10}}}(R) + 0.00063R + 2.04 - S $$ The anelastic attenuation term derived from this study indicates that ground motion attenuation is significantly different from Southern California but comparable with stable continental regions.     

Lyapunov exponent and dimension of the strange attractor of elastic frictional system

ＬｙａｐｕｎｏｖｅｘｐｏｎｅｎｔａｎｄｄｉｍｅｎｓｉｏｎｏｆｔｈｅｓｔｒａｎｅａｔｔｒａｃｔｏｒｏｆｅｌａｓｔｉｃｆｒｉｃｔｉｏｎａｌｓｙｓｔｅｍＺｈｉ－ＲｅｎＮＩＵ（牛志仁）ａｎｄＤａｎｇ－ＭｉｎＣＨＥＮ（陈党民）（ＳｅｉｓｍｏｌｏｇｉｃａｌＢｕｒ...     

The surface roughness and planetary boundary layer

Applications of the entrainment process to layers at the boundary, which meet the self similarity requirements of the logarithmic profile, have been studied. By accepting that turbulence has dominating scales related in scale length to the height above the surface, a layer structure is postulated wherein exchange is rapid enough to keep the layers internally uniform. The diffusion rate is then controlled by entrainment between layers. It has been shown that theoretical relationships derived on the basis of using a single layer of this type give quantitatively correct factors relating the turbulence, wind and shear stress for very rough surface conditions. For less rough surfaces, the surface boundary layer can be divided into several layers interacting by entrainment across each interface. This analysis leads to the following quantitatively correct formula compared to published measurements. 1 $$\begin{gathered} \frac{{\sigma _w }}{{u^* }} = \left( {\frac{2}{{9Aa}}} \right)^{{1 \mathord{\left/ {\vphantom {1 4}} \right. \kern-\nulldelimiterspace} 4}} \left( {1 - 3^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} \frac{a}{k}\frac{{d_n }}{z}\frac{{\sigma _w }}{{u^* }}\frac{z}{L}} \right)^{{1 \mathord{\left/ {\vphantom {1 4}} \right. \kern-\nulldelimiterspace} 4}} \hfill \\ = 1.28(1 - 0.945({{\sigma _w } \mathord{\left/ {\vphantom {{\sigma _w } {u^* }}} \right. \kern-\nulldelimiterspace} {u^* }})({z \mathord{\left/ {\vphantom {z L}} \right. \kern-\nulldelimiterspace} L})^{{1 \mathord{\left/ {\vphantom {1 4}} \right. \kern-\nulldelimiterspace} 4}} \hfill \\ \end{gathered} $$ where \(u^* = \left( {{\tau \mathord{\left/ {\vphantom {\tau \rho }} \right. \kern-0em} \rho }} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}} \) , σ _w is the standard deviation of the vertical velocity,z is the height andL is the Obukhov scale lenght. The constantsa, A, k andd _n are the entrainment constant, the turbulence decay constant, Von Karman's constant, and the layer depth derived from the theory. Of these,a andA, are universal constants and not empirically determined for the boundary layer. Thus the turbulence needed for the plume model of convection, which resides above these layers and reaches to the inversion, is determined by the shear stress and the heat flux in the surface layers. This model applies to convection in cool air over a warm sea. The whole field is now determined except for the temperature of the air relative to the water, and the wind, which need a further parameter describing sea surface roughness. As a first stop to describing a surface where roughness elements of widely varying sizes are combined this paper shows how the surface roughness parameter,z ₀, can be calculated for an ideal case of a random distribution of vertical cylinders of the same height. To treat a water surface, with various sized waves, such an approach modified to treat the surface by the superposition of various sized roughness elements, is likely to be helpful. Such a theory is particularly desirable when such a surface is changing, as the ocean does when the wind varies. The formula, 2 $$\frac{{0.118}}{{a_s C_D }}< z_0< \frac{{0.463}}{{a_s C_D (u^* )}}$$ is the result derived here. It applies to cylinders of radius,r, and number,m, per unit boundary area, wherea _s =2rm, is the area of the roughness elements, per unit area perpendicular to the wind, per unit distance downwind. The drag coefficient of the cylinders isC _D . The smaller value ofz _o is for large Reynolds numbers where the larger scale turbulence at the surface dominates, and the drag coefficient is about constant. Here the flow between the cylinders is intermittent. When the Reynolds number is small enough then the intermittent nature of the turbulence is reduced and this results in the average velocity at each level determining the drag. In this second case the larger limit forz ₀ is more appropriate.     

Steady flows driven by sources of random strength in heterogeneous aquifers with application to partially penetrating wells

Average steady source flow in heterogeneous porous formations is modelled by regarding the hydraulic conductivity K(x) as a stationary random space function (RSF). As a consequence, the flow variables become RSFs as well, and we are interested into calculating their moments. This problem has been intensively studied in the case of a Neumann type boundary condition at the source. However, there are many applications (such as well-type flows) for which the required boundary condition is that of Dirichlet. In order to fulfill such a requirement the strength of the source must be proportional to K(x), and therefore the source itself results a RSF. To solve flows driven by sources whose strength is spatially variable, we have used a perturbation procedure similar to that developed by Indelman and Abramovich (Water Resour Res 30:3385–3393, 1994) to analyze flows generated by sources of deterministic strength. Due to the linearity of the mathematical problem, we have focused on the explicit derivation of the mean head distribution G _d (x) generated by a unit pulse. Such a distribution represents the fundamental solution to the average flow equations, and it is termed as mean Green function. The function G _d (x) is derived here at the second order of approximation in the variance σ² of the fluctuation (where K _A is the mean value of K(x)), for arbitrary correlation function ρ(x), and any dimensionality d of the flow domain. We represent G _d (x) as product between the homogeneous Green function G _d ⁽⁰⁾(x) valid in a domain with constant K _A , and a distortion term Ψ _d (x) = 1 + σ²ψ _d (x) which modifies G _d ⁽⁰⁾(x) to account for the medium heterogeneity. In the case of isotropic formations ψ _d (x) is expressed via one quadrature. This quadrature can be analytically calculated after adopting specific (e.g.. exponential and Gaussian) shape for ρ(x). These general results are subsequently used to investigate flow toward a partially-penetrating well in a semi-infinite domain. Indeed, we construct a σ²-order approximation to the mean as well as variance of the head by replacing the well with a singular segment. It is shown how the well-length combined with the medium heterogeneity affects the head distribution. We have introduced the concept of equivalent conductivity K ^eq(r,z). The main result is the relationship where the characteristic function ψ^(w)(r,z) adjusts the homogeneous conductivity K _A to account for the impact of the heterogeneity. In this way, a procedure can be developed to identify the aquifer hydraulic properties by means of field-scale head measurements. Finally, in the case of a fully penetrating well we have expressed the equivalent conductivity in analytical form, and we have shown that (being the effective conductivity for mean uniform flow), in agreement with the numerical simulations of Firmani et al. (Water Resour Res 42:W03422, 2006).     

The viscosity of hydrous dacitic liquids: implications for the rheology of evolving silicic magmas

The viscosity of a series of six synthetic dacitic liquids, containing up to 5.04 wt% dissolved water, was measured above the glass transition range by parallel-plate viscometry. The temperature of the 10¹¹ Pa s isokom decreases from 1065 K for the anhydrous liquid, to 864 K and 680 K for water contents of 0.97 and 5.04 wt% H₂O. Including additional measurements at high temperatures by concentric-cylinder and falling-sphere viscometry, the viscosity (η) can be expressed as a function of temperature and water content w according to: where η is in Pa s, T is temperature in K, and w is in weight percent. Within the conditions of measurement, this parameterization reproduces the 76 viscosity data with a root-mean square deviation (RMSD) of 0.16 log units in viscosity, or 7.8 K in temperature. The measurements show that water decreases the viscosity of the dacitic liquids more than for andesitic liquids, but less than for rhyolites. At low temperatures and high water contents, andesitic liquids are more viscous than the dacitic liquids, which are in turn more viscous than rhyolitic liquids, reversing the trend seen for high temperatures and low water contents. This suggests that the relative viscosity of different melts depends on temperature and water content as much as on bulk melt composition and structure. At magmatic temperatures, rhyolites are orders of magnitude more viscous than dacites, which are slightly more viscous than andesites. During degassing, all three liquids undergo a rapid viscosity increase at low water contents, and both dacitic and andesitic liquids will degas more efficiently than rhyolitic liquids. During cooling and differentiation, changing melt chemistry, decreasing temperature and increasing crystal content all lead to increases in the viscosity of magma (melt plus crystals). Under closed system conditions, where melt water content can increase during crystallization, viscosity increases may be small. Conversely, viscosity increases are very abrupt during ascent and degassing-induced crystallization.     

Abiogenic hydrocarbons in commercial gases from the Songliao Basin, China

This paper discusses the kinetic fractionation, composition and distribution characteristics of carbon and hydrogen isotopes for various alkane gases formed in different environments, by different mecha- nisms and from different sources in nature. It is demonstrated that the biodegradation or thermode- gradation of complex high-molecule sedimentary organic material can form microbial gas or thermogenic gas. The δ 13C1 value ranges from -110‰ to -50‰ for microbial gases but from -50‰ to -35‰ (even heavier) f...     

Empirische Untersuchung zur Frage der Beziehung zwischen durchschnittlicher und kennzeichnender Wellenperiode im Seegang

Zusammenfassung Als Zusammenhang zwischen der kennzeichnenden Wellenperiode und der durchschnittlichen Periode im Seegang wird die Formel angesetzt. Mit Hilfe empirischer Unterlagen wird nachgewiesen, daßc eine Funktion des von D. E. Cartwright und M. S. Longuet-Higgins [1956] eingeführten Spektralparameters ist. Es wird eine vorläufige quantitative Beziehung zwischenc und abgeleitet.

---

Empirical investigations of the relation between the mean and the significant wave period in the sea

Summary It is supposed that the formula represents the relation between the significant wave period and the mean period in the sea. With the aid of empirical data it is demonstrated thatc is a function of the spectral parameter introduced by D. E. Cartwright and M. S. Longuet-Higgins [1956]. A preliminary quantitative relation betweenc and is derived.

---

Etudes empiriques de la relation entre la période moyenne et la période significative des vagues dans la houle

Résumé On suppose que la formule représente la relation entre la période significative des vagues et la période moyenne dans la houle. A l'aide des données empiriques on montre quec est une fonction du paramètre spectral , introduit par D. E. Cartwright et M. S. Longuet-Higgins [1956]. Une relation quantitative préliminaire entrec et est dérivée.

---

     

Lg magnitude scale in the mainland of China

The data of short-period seismograms had been collected widely in the mainland area of China not including Xinjiang and Tibet. The physical quantities of Lg wave are determined respectively in the five subregions. The group velocities of priminary arrival and maximum amplitude of Lg wave are equal to 3.54±0.02 km/s and 3.30±0.05 km/s, respectively. The periods of Lg waves are between 0.2s to 1.2s, averaging 0.7s. The γ-values of Lg waves in the five subregions are equal to 0.0034±0.0001 km⁻¹ for East, 0.0031±0.0004 km⁻¹ for Southwest, 0.0027±0.0004 km⁻¹ for Northeast, 0.0022±0.0001 km⁻¹ for South, and 0.0021±0.0002 km⁻¹ for Northwest subreqion, respectively. The average γ-value for the five subregions, γ=0.0027±0.0006 km⁻¹. The relations among the amplitude ratioH/Z, the station correctionD _z andD _h of amplitudes, and among them and station site condition are discussed. The subregional magnitude calibration functions ofm _Lg had been established according to each regional γ-value. From these, the unified magnitude calibration function of Chinese mainland not including Xinjiang and Tibet was given by

Variation of Seismic Coda Wave Attenuation in the Garhwal Region, Northwestern Himalaya

Seismic coda wave attenuation ( $ Q_{\text{c}}^{ - 1} $ ) characteristics in the Garhwal region, northwestern Himalaya is studied using 113 short-period, vertical component seismic observations from local events with hypocentral distance less than 250?km and magnitude range between 1.0 to 4.0. They are located mainly in the vicinity of the Main Boundary Thrust (MBT) and the Main Central Thrust (MCT), which are well-defined tectonic discontinuities in the Himalayas. Coda wave attenuation ( $ Q_{\text{c}}^{ - 1} $ ) is estimated using the single isotropic scattering method at central frequencies 1.5, 3, 5, 7, 9, 12, 16, 20, 24 and 28?Hz using several starting lapse times and coda window lengths for the analysis. Results show that the ( $ Q_{\text{c}}^{ - 1} $ ) values are frequency dependent in the considered frequency range, and they fit the frequency power law ( $ Q_{\text{c}}^{ - 1} \left( f \right) = Q_{0}^{ - 1} f^{ - n} $ ). The Q ₀ (Q _c at 1?Hz) estimates vary from about 50 for a 10?s lapse time and 10?s window length, to about 350 for a 60?s lapse time and 60?s window length combination. The exponent of the frequency dependence law, n ranges from 1.2 to 0.7; however, it is greater than 0.8, in general, which correlates well with the values obtained in other seismically and tectonically active and highly heterogeneous regions. The attenuation in the Garhwal region is found to be lower than the Q _c ^?1 values obtained for other seismically active regions of the world; however, it is comparable to other regions of India. The spatial variation of coda attenuation indicates that the level of heterogeneity decreases with increasing depth. The variation of coda attenuation has been estimated for different lapse time and window length combinations to observe the effect with depth and it indicates that the upper lithosphere is more active seismically as compared to the lower lithosphere and the heterogeneity decreases with increasing depth.     

Comparative soil CO2 flux measurements and geostatistical estimation methods on Masaya volcano, Nicaragua

We present a comparative study of soil CO₂ flux () measured by five groups (Groups 1–5) at the IAVCEI-CCVG Eighth Workshop on Volcanic Gases on Masaya volcano, Nicaragua. Groups 1–5 measured using the accumulation chamber method at 5-m spacing within a 900 m² grid during a morning (AM) period. These measurements were repeated by Groups 1–3 during an afternoon (PM) period. Measured ranged from 218 to 14,719 g m⁻² day⁻¹. The variability of the five measurements made at each grid point ranged from ±5 to 167%. However, the arithmetic means of fluxes measured over the entire grid and associated total CO₂ emission rate estimates varied between groups by only ±22%. All three groups that made PM measurements reported an 8–19% increase in total emissions over the AM results. Based on a comparison of measurements made during AM and PM times, we argue that this change is due in large part to natural temporal variability of gas flow, rather than to measurement error. In order to estimate the mean and associated CO₂ emission rate of one data set and to map the spatial distribution, we compared six geostatistical methods: arithmetic and minimum variance unbiased estimator means of uninterpolated data, and arithmetic means of data interpolated by the multiquadric radial basis function, ordinary kriging, multi-Gaussian kriging, and sequential Gaussian simulation methods. While the total CO₂ emission rates estimated using the different techniques only varied by ±4.4%, the maps showed important differences. We suggest that the sequential Gaussian simulation method yields the most realistic representation of the spatial distribution of , but a variety of geostatistical methods are appropriate to estimate the total CO₂ emission rate from a study area, which is a primary goal in volcano monitoring research.Editorial responsibility: H Shinohara