Intrinsic and scattering attenuation in western India from aftershocks of the 26 January, 2001 Kachchh earthquake

176 vertical-component, short period observations from aftershocks of the M_w 7.7, 26 January, 2001 Kachchh earthquake are used to estimate seismic wave attenuation in western India using uniform and two layer models. The magnitudes (M_w) of the earthquakes are less than 4.5, with depths less than 46 km and hypocentral distances up to 110 km. The studied frequencies are between 1 and 30 Hz. Two seismic wave attenuation factors, intrinsic absorption (Q_i^− 1) and scattering attenuation (Q_s^− 1) are estimated using the Multiple Lapse Time Window method which compares time integrated seismic wave energies with synthetic coda wave envelopes for a multiple isotropic scattering model. We first assume spatial uniformity of Q_i^− 1, Q_s^− 1 and S wave velocity (β). A second approach extends the multiple scattering hypothesis to media consisting of several layers characterized by vertically varying scattering coefficient (g), intrinsic absorption strength (h), density of the media (ρ) and shear wave velocity structure. The predicted coda envelopes are computed using Monte Carlo simulation. Results show that, under the assumption of spatial uniformity, scattering attenuation is greater than intrinsic absorption only for the lowest frequency band (1 to 2 Hz), whereas intrinsic absorption is predominant in the attenuation process at higher frequencies (2 to 30 Hz). The values of Q obtained range from Q_t = 118, Q_i = 246 and Q_s = 227 at 1.5 Hz to Q_t ≈ 4000, Q_i ≈ 4600 and Q_s ≈ 33,300 at 28 Hz center frequencies, being Q_t^− 1 a measure of total attenuation. Results also show that Q_i^− 1, Q_s^− 1 and Q_t^− 1 decrease proportional to f^−ν. Two rates of decay are clearly observed for the low (1 to 6 Hz) and high (6 to 30 Hz) frequency ranges. Values of ν are estimated as 2.07 ± 0.05 and 0.44 ± 0.09 for total attenuation, 1.52 ± 0.21 and 0.48 ± 0.09 for intrinsic absorption and 3.63 ± 0.07 and 0.06 ± 0.08 for scattering attenuation for the low and high frequency ranges, respectively. Despite the lower resolution in deriving the attenuation parameters for a two layered crust, we find that scattering attenuation is comparable to or smaller than the intrinsic absorption in the crust whereas intrinsic absorption dominates in the mantle. Also, for a crustal layer of thickness 42 km, intrinsic absorption and scattering estimates in the crust are lower and greater than those of the mantle, respectively.     

Attenuation of shear-waves in the back-arc region of the Hellenic arc for frequencies from 0.6 to 16 Hz

Q_β for shear-waves is determined for the inner part of the Hellenic arc, the back-arc area, as a function of frequency in the range 0.6–16 Hz. We used 314 digital records from 32 earthquakes with magnitudes (M_w) ranging from 3.9 to 5.1. Epicentral distances ranged from 65 to 515 km. The data were obtained in 1997 during a 6-month operation of a digital portable network in Greece. The Q_β estimates were made for five frequency bands centred at 0.8, 1.5, 3.0, 6.0 and 12.0 Hz and the Q_β values obtained were 47, 79, 143, 271 and 553, respectively. The results show that Q_β for S-waves increases with frequency taking the form Q_β=55f ^0.91 (or Q_β⁻¹0.018f^−0.91). The high attenuation and the strong frequency dependence found, which is close to the frequency dependence of coda Q for Greece, are characteristic of an area with high seismicity, rapid extension, and in agreement with other similar studies in Greece.     

Estimation of the local financial costs of flood damage with different methodologies in Unye (Ordu), Turkey

The aim of this study was to conduct flood analysis with digital modeling systems and estimate the total flood damage in a built-up area. The Ceviz Stream basin, which is located in Unye (Ordu) district in the Eastern Black Sea region of Turkey, was selected as the study area. A 1D/2D coupled flood modeling software MIKE FLOOD was used in flood analysis. According to the modeling results, water depth and velocity maps are produced for maximum Q₅₀, Q₁₀₀, Q₅₀₀ and Q₁₀₀₀ discharge scenarios in the most downstream section of the Ceviz Stream. At the projected Q_50, Q₁₀₀, Q₅₀₀ and Q₁₀₀₀ maximum flow rates, approximately 9.75%, 23.0%, 30.0% and 32.25% of that section of the study area were estimated to be negatively affected by floodwaters. In order to examine the financial magnitude of the effects of flooding, the Huizinga, van Eck and Kok, ICBR and Pistrika and Jonkman damage estimation methodologies were applied, with the highest financial costs estimated with the Pistrika and Jonkman method.

     

Assessment of direct CPT and CPTU methods for predicting the ultimate bearing capacity of single piles

Twelve methods to determine axial pile capacity directly based on cone penetration test (CPT) and piezocone penetration test (CPTU) data are presented, compared and evaluated. Analyses and evaluation were conducted on three types of piles of different size and length. All the tested piles have failed at the end of static load test. Both the CPT methods and the CPTU methods were used to estimate the load bearing capacities of the investigated piles (Q_p). The static load test was performed to determine the measured load bearing capacities (Q_m). The pile capacities determined through different methods were compared with the measured values obtained from the static load tests. Two criteria were selected as bases of evaluation: the best fit line for Q_p versus Q_m and the arithmetic mean and standard deviation for the ratio Q_p/Q_m. Results of the analyses showed that the best methods for determining pile capacity are the two CPTU methods. Furthermore, the CPTU method is simple, easy to apply, and not influenced by the subjective judgements of operating staff. Therefore, it is quite suitable for the application in pile engineering practice.     

Attenuation studies based on local earthquake Coda waves in the southern Indian peninsular shield

The Dharwar craton, Cuddapah basin and the Godavari graben characterise three diverse geological and tectonic settings in the peninsular shield of India. Owing to their contiguous proximity, they offer a unique opportunity to document differences, if any, in their seismic wave attenuation characteristics that might have a bearing on the seismogenic nature of the crust in a craton, basin and a rift-like graben structure. An attempt is made here to bring out these differences using constraints from coda-Q. We considered local earthquakes with epicentral distances ranging from␣14 to 150 km recorded at the digital broadband stations at Dharwar (DHD), Cuddapah (CUD) and Kothagudem (KGD) regions to derive the frequency-dependent coda-Q relations. Using the single scattering method, we obtained the frequency-dependent Q _C relationship (Q _C = Q ₀ f ⁿ )for each of the three geological units separately: DHD: Q _C = (730.62 ± 0.09)f ^{(0.54 ± 0.01)}; CUD: Q _C = (535.06 ± 0.13)f ^{(0.59 ± 0.01)} and KGD: Q _C = (150.56 ± 0.08)f ^{(0.91 ± 0.01)}. The Q _C values obtained for all the three sub-regions show moderate to strong frequency dependence and essentially reflect the level of crustal heterogeneities to varying degrees.     

On the Study of Pile Driving Formula Based on Blow Counts

Accuracy of predicting pile capacities by pile driving formulas have been investigated. Five test piles were driven up to a depth of about 9 m of clay deposit and the penetrations due to final blows were recorded. The pile bearing capacity of each pile was predicted using 6 different pile driving formulas and the predicted pile capacity was compared with measured pile capacity from the pull up tests. Hiley formula, Modified Engineering News Record (ENR) formula, Janbu formula, Dutch formula, Danish formula, and Gates formula were used. The performance and accuracy of each formula was evaluated and the correlation coefficient of each pile driving formula was determined for a more accurate pile capacity prediction. Methods used to evaluate the performance of each formula were; (1) the best fit line for Q _p versus Q _m (2) cumulative probability for Q _p/Q _m and (3) the arithmetic mean and standard deviation for Q _p/Q _m. From the study, it was found that using Dutch formula provided the most accurate pile capacity estimate compared to the other formulas with an average of 7% deviation from value obtained from the field pull up test. It was followed by the Danish formula, Janbu formula, Hiley formula, Modified ENR formula, and Gates formula. The ability to predict the accuracy of estimating pile capacity using an appropriate method is very important and valuable to contractors, developers, geotechnical engineers, and manufacturers.     

The effect of thermal treatment on the 57Fe Mössbauer spectrum of beryl

Mössbauer spectra (MS) of blue, green and yellow beryl (ideally Be₃Al₂Si₆O₁₈) containing approximately 1% of iron were obtained at 295 and 500 K. Room temperature (RT) spectra of both blue and green samples showed the presence of an asymmetric Fe²⁺ doublet (ΔE _Q~2.7 mm/s, δ~1.1 mm/s), with a very broad low-velocity peak. There is no clear evidence for the presence of a ferric component. The MS of the yellow sample at RT consists of an intense central absorption with parameters typical for Fe³⁺ (ΔE _Q~0.4 mm/s, δ~0.29 mm/s), plus an apparently symmetrical Fe²⁺ doublet. This sample acquires a light-blue shade upon heating in air at about 620 K. Thermal treatments at high temperatures caused no significant changes in the MS, but the green and yellow beryl acquire a blue colour. All these results are interpreted in relation to the existence of channel water and the distribution of iron among the available crystallographic sites.     

Temperature dependence of the polarized electronic absorption spectra of olivines. Part I – fayalite

 Polarized electronic absorption spectra of orthorhombic fayalite, Fe₂SiO₄, [E || a(|| Z),E || b(|| X), E || c(|| Y)], space group Pbnm, have been studied in the temperature range 293 ≤T/K ≤ 1273. The spectra were analysed into component bands originating from spin-allowed dd transitions of iron(II) at the different sites, M1 and M2, in the structure. The assignments of bands, made on the basis of the polarization dependence of the spectra and considerations of transition energies, were confirmed by the analysis of the temperature-dependent spectra. The temperature dependencies of integral intensities, half band widths and energy positions of absorptions bands caused by Fe²⁺ on the different octahedral sites, M1 and M2, were evaluated for the individual transitions. Independent of the site symmetry, absorption bands shift to lower energies and half band widths increase on rising temperature. The temperature dependence of band intensities depends on site symmetry. The integral intensities are found to increase with temperature for the transition metal ion on a centrosymmetric site, or remain constant when the site is missing an inversion centre. This is consistent with the general conclusion of Taran et al. (1994). Received: 11 October 2001 / Accepted: 17 January 2002     

A statistical-mechanics model of ordering in aluminosilicate solid solutions

A statistical-thermodynamic approximation for order-disorder phase transitions in aluminosilicate solid solutions is presented. The approximation involves estimating the number of configurations with long range order parameter Q and short range order parameter σ, using an approximation to the probability that a configuration with the correct amount of long range order will also have the correct amount of short range order. This estimate is then used to give a free energy F(Q, σ, x, T), where x is the concentration of Al in the structure, and hence quantities such as T _c (x) are estimated. The predictions of the model for T _c (x) and the critical concentration x _c at which T _c falls to zero are shown to be in good agreement with the results of Monte Carlo simulations. Received: 14 May 1997 / Revised, accepted: 2 June 1997     

CODA Q estimates for Kumaun Himalaya

CodaQ (Q _c) estimates for the Kumaun Himalaya region have been obtained in high frequency range. Local earthquakes, recorded by a digital seismic network in the region, which fall in the epicentral distances range of 10 to 80 km and with a local magnitude range of 1.4 to 2.8, have been used. The coda waves of 30 sec window length, filtered at seven frequency bands centered at 1.5, 3, 6, 9, 12, 18 and 24Hz, have been analysed using the single backscattering model. The values ofQ _c estimates vary from 65 to 283 at 1.5 Hz to 2119 to 3279 at 24.0 Hz which showed thatQ _c is frequency dependent and its value increases as frequency increases. A frequency-dependentQ _c relationship,Q _c = (92 ± 4.73)f^(1.07±.023), is obtained for the region representing the average attenuation characteristics of seismic waves for Kumaun Himalaya region.     

Heat flow and heat generation in the Churchill Province of the Canadian Shield, and their palaeotectonic significance

Six new heat flow determinations are presented for Proterozoic mobile belts of the Churchill Province of the Canadian Shield, an area that was affected by several stages of the Hudsonian orogenic sequence (1.9-1.6 Ga ago). With other, previously published, values the mean of eight determinations considered reliable and representative and corrected for the effects of Pleistocene glaciation is 44 ± 7 mW m⁻². Heat generation measurements have also been made; values range from 0.1–1.04 μW m⁻³.A linear relation between heat flow and heat production is apparent. The heat flow axis intercept is 37 mW m⁻², and the scale depth is 11 km, compared with 28 mW m⁻² and 13.6 km for the Archaean Superior Province. Approximately 20% of the Churchill heat flow appears to be derived from radioactive decay in the upper crust, compared with 30% for the Superior Province and shields as a whole. The observations imply that the heat flow-heat production relation for the Churchill Province should be written as Q = Q_c + Q_e + A₀b where Q_c is equivalent to the reduced heat flow for the Archaean terrain, b is similar for the two, and Q_e is an additional component of heat flow in the Proterozoic mobile belts of the Churchill Province.A speculative tectonic model is presented. It is suggested that rifting along two axes of an original craton, which had lateral variations in near surface radiogenic element concentration, followed by erosion of the radiogenic layer and subsequent reconvergence of the cratonic segments, led to widespread redistribution of radioactive elements into the reactivated inter-rift crustal block. One result would be that crustal temperatures are higher in that part of the Churchill Province than in the Superior.     

河套盆地东部第四纪中晚期介形类特征及其沉积环境意义

根据河套盆地东部DTLC孔介形类组合及其岩性、岩相的分析,着重地探讨了该区第四纪中、晚期介形类的生物地层特征及其沉积古环境意义。研究认为:中更新世早期（Q₂1）,未见介形类,以冲洪积沉积环境为主,兼有河湖交替演化;中更新世晚期（Q₂2）,介形类生物组合为Leucocythere plethora-Eucypris inflate-Limnocyther dubiosa,推测为封闭性的半咸水、咸水湖相沉积环境,气候偏冷;晚更新世早期（Q₃1）,介形类生物组合为Leucocythere plethora-Limnocythere dubios,反映了封闭性的半咸水深湖相沉积环境;晚更新世晚期（Q₃2）,少见淡水介形虫化石,以Ilyocypris manasensis-Candoniella albicans-Cypria luminosa为组合特征,反映了以河流冲洪积环境为主;全新世（Q₄）为河流沉积环境。通过河套盆地区域沉积环境的综合分析认为:中更新世晚期（Q₂2）至晚更新世早期（Q₃1）,河套盆地为封闭性内陆湖盆,黄河为内流河;自晚更新世晚期（Q₃2）以来,黄河逐渐地转为外流河,河套盆地为外流湖盆。      

Parametric studies of physical conditions in the Earth's upper mantle

Seismic body-wave and surface-wave data indicate the existence of a substantial lowvelocity, low-Q zone in the upper mantle beneath western North America. Conditions in this zone are distinctively different from those that are typical of the upper mantle in shield regions. The present study, using Walsh's model for partially melted rock, suggests a common mechanism for low-velocity and low-Q zones. This parametric study also indicates that the pronounced low-velocity, low-Q zone and anomalous travel-time delays of both P and S waves in the Basin and Range Province are consistent with the combined effects of high temperature, chemical composition, phase changes, and partial melting. The observed low P_n velocity in this region is consistent with high temperature, chemical composition, and the presence of a partially molten layer within the upper mantle, but the observed teleseismic delay times result principally from the thickness of the lowvelocity zone. The teleseismic delay-time variations are therefore related to the seismic Q distribution in the asthenosphere. Conditions a few kilometers beneath the Moho boundary influence the P_n velocity; however, the observed correlation among the teleseismic signal amplitudes, travel-time delays, and the upper mantle Q indicates that the P_n velocity is a better indicator of upper mantle Q than suggested by the P_n path alone. This knowledge of the upper mantle seems to account for the anomalous effects of the Basin and Range Province and other regions of similar tectonic nature on observations of teleseismic events. This will provide a practical technique for comparing seismic observations made in unstable tectonic regions with observations made in shield areas.     

Large-scale pile tests in Mercia mudstone: Data analysis and evaluation of current design methods

Full-scale load tests were carried out on six instrumented large diameter bored, cast in-situ piles formed in Mercia mudstone, as part of the design of a new Viaduct in Cardiff, UK. In this paper, the results from six test piles and extensive data from 218 ground investigation boreholes are systematically processed in order to study the load transfer and resistance mechanisms in Mercia mudstone. Data from strain gauges embedded in each pile are first analysed to calibrate the load-deformation relationship of each pile as-built, taking into account (i) the non-linearity of concrete and (ii) the effect of partial steel encasement on pile stiffness at various levels. The shaft and base capacity of the piles are each predicted using 10 calculation methods belonging to the four basic categories: (i) Undrained analysis, (ii) Drained analysis, (iii) Mixed approach and (iv) Empirical correlation. It is found that the shaft capacity prediction methods are moderately consistent. The standard deviations of the ratio Q _sp/Q _sm of predicted to observed shaft capacity lies in the range 0.06–0.24. However, 8 of these methods are over-conservative, giving Q _sp/Q _sm values in the range 0.29–0.67. The remaining two methods yield Q _sp/Q _sm = 1.01 and 1.49. In contrast, the prediction methods for base capacity are found to be much less consistent. The ratio Q _bp/Q _bm of predicted to measured base capacity falls in the interval 0.52–1.93, with corresponding standard deviations of 0.16–0.82. This revised version was published online in July 2006 with corrections to the Cover Date.     

Determination of the degree of Al,Si order Q od in kinetically disordered albite using hard mode infrared spectroscopy

Hard Mode Infrared Spectroscopy (HMIS) is used to correlate the line shifts Δω, the intensity changes ΔA and the variations of spectral line widths Γ of infrared absorption bands with the degree of Al, Si ordering, Q _od, in kinetically disordered Na-feldspar. A simple relationship Δω∝ΔA∝ΔΓ∝Q _od ² was found with the phonon band around 650 cm^?1 being particularly sensitive to small changes of the degree of Al, Si ordering. It is shown that the average degree of Al, Si order can be determined from HMIS with an accuracy of ca. 8 percent using 50 mg of Na-feldspar. The experimental results agree well with recent X-ray determinations using identical samples. The significance of HMIS for the study of kinetic processes in minerals is explained.     

Lateral variations of attenuation coefficients, group and phase velocities of rayleigh waves in Europe

Rayleigh wave group and phase velocities and attenuation coefficients are investigated in the 15–60-s period range for earthquakes occurring in the European area and registered at European stations. Observed group and phase velocities show clear differences in the period range 20–40 s, whereas attenuation coefficients are differentiated from 45 s to larger periods. From the shape of dispersion and attenuation curves, two broad zones have been isolated and designated as Western Europe and Southeastern Europe. Although some scatter is present, specially for attenuation, the data suggest that in general the shear velocity and Q_β values are higher for Southeastern Europe than for Western Europe. For Southeastern Europe there is a channel of low velocity and low Q_β, apparently correlated, located at the lower crust: in both regions there is a channel of low velocity and low Q_β values located at the uppermost mantle. In each case there also exists a channel of high Q_β values located at the lower crust but at slightly different depths.     

松科二井深地震反射数据Q值分析及其对含气层位的确定

地震波在地下介质中传播时,地震波能量会出现一定程度的衰减,品质因子Q作为衡量地下岩石吸收衰减属性的重要参数,对描述岩性特征以及预测油气分布具有重要意义。本文针对深地震反射数据,利用基于S变换谱比法的Q值分析技术,获得了更加准确的Q值。以松辽盆地沙河子组为主要目标层,对穿过松科二井的叠后深地震反射剖面进行Q值计算,生成Q值剖面,总结出沙河子组Q值分布特征,同时结合松科二井测井、分层和气测异常资料进行分析,推断沙河子组为含气储层,为下一步深部储层预测提供思路,为东北地区深地勘查工程提供服务。     

Explosion seismic P and S velocity and attenuation constraints on the lower crust of the North–Central Tibetan Plateau, and comparison with the Tethyan Himalayas: Implications on composition, mineralogy, temperature, and tectonic evolution

P and S velocity and attenuation estimates in the lower crust are obtained from a set of wide angle reflection–refraction profiles in the region of active tectonics at the NE edge of the Tibetan Plateau and discussed together with respect to similar data at its Himalaya–south Tibet edge.The quality factor is estimated in the lower half of the crust by accounting for the differential effect on amplitude–frequency observed between waves of different penetrations, and both in P and S modes. Attenuation values allow to exclude a significant proportion of partial melt and to estimate the homologous temperature, ratio of in situ to solidus absolute temperatures. The latter depend on the physical conditions being of dry, wet or dehydration melting, which are found different among the regions of the northern Bayan Har and northern Qang Tang boundaries between blocks, as well as the Tethyan–Himalayas, south of the Indus–Tsangpo suture. Their in situ temperatures differ also as estimated from their different V_p for a similar felsic composition.Joint measurement of several parameters, V_p, V_s, Q_p and Q_s reveals the composition, mineralogy, temperature and hydration conditions of the lower half of the thickened crust of Tibet that may be discussed in terms of evolution. The material presently in the thickened crust, even its lower part, has a felsic composition, upper to middle crustal lithology, and the temperature conditions estimated suggest that basic material that could have underlain it could be eclogitized and not appear anymore above the seismic Moho.Under northern Qang Tang, the felsic material in the lower half of the crust appears as hot and dry. Its burial may have occurred earlier or may have been moderate in the postcollisional phase. This is consistent with a model of indentation of the Qang Tang crust by an originally thinner Bayan Har crust to bring part of its crust to greater depth, suggested from imaging the crustal architecture. Under northern Bayan Har, the material in the lower half of the crust appears as felsic, at low temperature and not dry conditions. This is evidence that it has been transported from a shallower depth, and this recently enough not to be yet dehydrated and temperature equilibrated in a conductive geotherm. It supports a model of recent overriding of the middle crust of the north Kun Lun block to the north independently suggested from the image of crustal architecture. The Tethyan Himalayas case appears bracketed by these two cases in northern Tibet for V_p and temperature conditions, but shows highest attenuation in the lower crust that is colder but less dry than under northern Qang Tang.     

Attenuation of coda waves in the Nubra-Siachen region,Himalaya, India

The attenuation properties have been estimated in the Nubra-Siachen region situated in the highly mountainous region of Himalayan belt. Coda wave quality factor (Q_c) has been determined for this virgin region by using the single backscattering method. A total of thirty earthquakes recorded in this region, which fall in the epicentral distance range of 3 to 115km have been used for the present work. A 30 sec window length of coda waves at different central frequencies 1.5, 3.0, 6.0, 9.0, 12.0, 18.0 and 24.0 Hz have been studied to determine Q_c at different recording stations. The frequency dependent coda wave quality factor relationships of the form Q_c(f) = Q_ofⁿ, have been computed at each recording stations separately: BASE: Q_c(f) = (137 ± 4.2) f^{(0.99 ± 0.12)}, CHALUNKA: Q_c(f) = (116±3.8)f^(1.0±0.05), PARTA: Q_c(f) = (122±3.0)f^(1.0±0.02), and SASOMA: Q_c(f) = (111±4.1)f^(1.0±0.03). A regional Q_c relation has been developed for the Nubra-Siachen region by using the average value of Q_c at different frequencies obtained at each recording station of the form Q_c(f) = (121±7.2)f^(1.0±0.04). The average Q_c values vary from 183 at 1.5 Hz to 3684 at 24 Hz central frequencies. The present regional relation developed for Nubra-Siachen region indicates heterogeneous and tectonically active region.     

Evidence for lithospheric detachment in the central Andes from local earthquake tomography

Data from three temporary seismic networks were merged for tomographic inversion. Although the deployments did not coincide in time, spatial overlap was achieved by re-occupying existing sites. Travel times and t operators of about 1600 earthquakes were inverted for 3D models of ν_p, ν_p/ν_s and P-wave attenuation (Q_p^− 1). All three attributes provide a consistent image of the entire subduction zone on a lithospheric scale. The tomographic images reveal low velocities and high attenuation in the crust and mantle underlying the Western Cordillera and most of the Puna plateau, indicative of weak rheology and mostly asthenospheric mantle. In contrast, forearc and eastern foreland are characterized by high Q_p values, corresponding to cold temperatures in accordance with thermal models. In the backarc, between 23°S and 24°S, a high velocity, high Q_p structure beneath the Eastern Cordillera and eastern Puna is interpreted as detaching continental lithosphere that has been thickened in the orogenic process. South of this structure, the mantle is characterized by low velocities, high ν_p/ν_s ratios, and low Q_p values. Here it is believed that lithosphere originally underlying Andean crust has already been removed. This is supported by new estimates of crustal thickness and volcanic activity.