High-speed correction factor to the O+-O resonance charge exchange collision frequency

The high-speed correction factor to the O⁺-O collision frequency, resulting from drift velocities between ions and neutrals, is calculated by solving the integral expression in this factor both numerically and analytically. Although the analytic solution is valid for either small or large drift velocities between ions and neutrals, for temperatures of interest and all drift velocities considered, agreement is found between analytic and detailed numerical integration results within less than 1% error. Let T _r designate the average of the ion and neutral temperatures in K, and u=v _d /<alpha>, where v _d is the relative drift velocity in cm s^?1, and <alpha>=4.56×10³\sqrtT _r cm s^?1 is the thermal velocity of the O⁺-O system. Then, as u ranges from 0 to 2, the correction factor multiplying the collision frequency increases monotonically from 1 to about 1.5. An interesting result emerging from this calculation is that the correction factor for temperatures of aeronomical interest is to a good approximation independent of the temperature, depending only on the scaled velocity u.     

Excess bottom222Rn profiles and their implications in the northwestern Pacific Ocean

Vertical profiles of excess bottom²²²Rn and potential temperature were measured at 23 stations in the northwestern Pacific Ocean. The Rn profiles were classified into the following three types: quasi-exponential (type E), benthic boundary layer (type B), and horizontally disturbed (type H). The ratio among types E, B and H was approximately 2 : 3 : 1.An apparent vertical eddy diffusivity (K) was calculated by applying a one-dimensional diffusion model to the Rn profiles of types E and B. Type E had K values ranging from 15 to 180 cm² s^?1 (average: 70 cm² s^?1). As to type B, K values for the benthic boundary layer (4.5–260 cm² s^?1, average: 120 cm² s^?1) were always more than an order of magnitude larger than those for the upper layer (0.2–35 cm² s^?1, average: 7 cm² s^?1), indicating more active vertical mixing in the benthic boundary layer than in the upper layer.Rn profiles were measured in regions where the bottom topography is known. It was verified that the occurrence of type H related closely with local bottom topographic features accompanied by lateral transient supply of Rn-rich or Rn-poor water.A couple of Rn profiles at the same location, measured at time intervals of several years, were compared with each other for three locations. The general characteristics of Rn profiles were shown to remain unaltered with time, while the fine structure of Rn profiles may have short-term variations caused by local bottom topography and fluctuations of bottom current as indicated in type H.     

Nucleosynthetic osmium isotope anomalies in acid leachates of the Murchison meteorite

We present osmium isotopic results obtained by sequential leaching of the Murchison meteorite, which reveal the existence of very large internal anomalies of nucleosynthetic origin (ε¹⁸⁴Os from ? 108 to 460; ε¹⁸⁶Os from ? 14.1 to 12.6; ε¹⁸⁸Os from ? 2.6 to 1.6; ε¹⁹⁰Os from ? 1.7 to 1.1). Despite these large variations, the isotopic composition of the total leachable osmium (weighted average of the leachates) is close to that of bulk chondrites. This is consistent with efficient large-scale mixing of Os isotopic anomalies in the protosolar nebula. The Os isotopic anomalies are correlated, and can be explained by the variable contributions of components derived from the s, r and p-processes of nucleosynthesis. Surprisingly, much of the s-process rich osmium is released by relatively mild leaching, suggesting the existence of an easily leachable s-process rich presolar phase, or alternatively, of a chemically resistant r-process rich phase. Taken together with previous evidence for a highly insoluble s-process rich carrier, such as SiC, these results argue for the presence of several presolar phases with anomalous nucleosynthetic compositions in the Murchison meteorite. The s-process composition of Os released by mild leaching diverges slightly from that released by aggressive digestion techniques, perhaps suggesting that the presolar phases attacked by these differing procedures condensed in different stellar environments. The correlation between ε¹⁹⁰Os and ε¹⁸⁸Os can be used to constrain the s-process ¹⁹⁰Os/¹⁸⁸Os ratio to be 1.275 ± 0.043. Such a ratio can be reproduced in a nuclear reaction network for a MACS value for ¹⁹⁰Os of ~ 200 ± 22 mbarn at 30 keV. More generally, these results can help refine predictions of the s-process in the Os mass region, which can be used in turn to constrain the amount of cosmoradiogenic ¹⁸⁷Os in the solar system and hence the age of the Galaxy.We also present evidence for extensive internal variation of ¹⁸⁴Os abundances in the Murchison meteorite. A steep anti-correlation is observed between ε¹⁸⁴Os and ε¹⁸⁸Os. Since ¹⁸⁴Os is formed uniquely by the p-process, this anti-correlation cannot be explained by variable addition or subtraction of s-process Os to average solar system material. Instead, this suggests that p-process rich presolar grains (e.g., supernova condensates) may be present in meteorites in sufficient quantities to influence the Os isotopic compositions of the leachates. Nevertheless, ¹⁸⁴Os is a low abundance isotope and we cannot exclude the possibility that the measured anomalies for this isotope reflect unappreciated analytical artifacts.     

The maximum-slope method for the interpretation of electrical soundings over two horizontal layers

Using the Carpenter configuration, in which the distances between the successive electrodes of a four-electrode system are kept at r, s and r (as shown in Fig. 1), where r = ma and s = la, or the Wenner-Carpenter configuration, in which r = s, the corresponding formulae for the apparent resistivities ρ_A^α, ρ_A^β and ρ_A^γ for the α, β and γ configurations are used to compute maximum gradient values for various values of the reflection coefficient K. For any observed curve for these resistivity types, the maximum slope is determined. A set of standard curves is presented which, provided that K is known, then allows the thickness of the top layer to be determined.In the case that K is unknown, another set of standard curves is presented allowing determination of K provided that the maximum slopes of any two of the three resistivity curves are known.The method has been tested for both field and laboratory observations. Discrepancies between the observed and actual values of resistivity are noted and attributed partly to instrumental errors.     

The effect of initial230Th disequilibrium on young UPb ages: the Makalu case,Himalaya

Comparative UPb dating of zircon, xenotime and monazite from two different samples of the Himalayan “Makalu” granite shows the two U decay series to be in disequilibrium, particularly in monazite. This disequilibrium is due to excess or deficit amounts of radiogenic²⁰⁶Pb which originate from an excess or deficit of²³⁰Th, respectively, occurring initially in the mineral. Such an initial disequilibrium is caused by UTh fractionation between the crystallising mineral and the magma. Therefore, the UPb ages of Th-rich minerals such as monazite (and allanite) have to be corrected for excess²⁰⁶Pb due to excess²³⁰Th, whereas Th-poor minerals such as zircon and xenotime require a correction for a deficit of²⁰⁶Pb due to deficiency of²³⁰Th. The extent of this correction depends on the degree of ThU fractionation and on the age of the rock. For the two monazite populations analysed here, these excess amounts of²⁰⁶Pb were, with reference to the amount of radiogenic²⁰⁶Pb, 8–10% and 15–20% respectively, and less than 1% for zircon and xenotime. The varying degrees of Th enrichment relative to U in monazite show that the ThU partition coefficients for this mineral are not constant within a single granite. Furthermore, for monazite there is evidence for excess amounts of radiogenic²⁰⁷Pb originating from the decay of initial excess²³¹Pa, also enriched during crystal growth.The very low Th/U ratios of 0.196 and 0.167, determined for thetwo whole rocks from which the minerals have been extracted, substantiate the view that granite formation is a fundamental mechanism for ThU fractionation in continental crust.The different ages of 21.9 ± 0.2m.y. and24.0 ± 0.4m.y., obtained by averaging the corrected²³⁸U²⁰⁶Pb ages of the monazites, suggest that the apparently homogeneous Makalu granite was generated over a period of at least 2 m.y.     

Radiation damping in the context of one-dimensional wave propagation: A teaching perspective

This paper has two objectives: to connect directly radiation damping and 1D elastic wave propagation, and to create a simple teaching tool to introduce the subject to students and engineers trained in Structural Dynamics. The first objective is achieved by obtaining the equivalent radiation modal damping using time domain solutions for the fundamental mode in shear of an elastic layer (soil) on flexible rock, for the case in which the rock–soil Impedance Ratio in shear, I=(ρ_rV_r)/(ρ_sV_s)>1, where ρ=mass density and V=shear wave velocity. These time domain solutions are developed for the case of steady-state input sinusoidal shear waves propagating vertically in the rock as well as for horizontal free vibration of the layer. Both derivations result in the same approximate expression for the modal radiation damping in the first mode, ζ_r1≈2/(πI), which is in turn identical to the approximate equation obtained by Roesset and Whitman [11] using a frequency domain approach. This expression for ζ_r1 is linked to the fact that, during free vibration, the ratio between two successive positive displacement peaks u_so and u_s1 at the ground surface is u_s0/u_s1=[(1+I)/(1−I)]², associated with two wave reflections at the soil–rock interface. From this ratio, and after applying the standard expression to obtain modal damping from damped free vibration, the same expression for ζ_r1 is reached again, ζ_r1≈[1/(2π)] ln(u_s0/u_s1)≈[1/(2π)] ln[(1+I)/(1−I)]²≈2/(πI). This finding allows development of the simple teaching tool proposed at the end of the paper. While only a crude approximation lacking in rigor, this teaching tool is physically intuitive, links directly wave propagation and modal damping in a simple way and gives the correct result.     

Attenuation of body waves in Southeastern Sicily (Italy)

The attenuation of P- and S-waves in Southeastern Sicily was estimated by applying two different methods in time and frequency domains. We analyzed waveforms from about 290 local events (0.6≤M_L≤4.6) recorded at a three-component digital network.By applying the pulse broadening method to the first P-wave pulse, we found an average Q_p value of ca. 140. The application of the frequency decay method provided a Q_p value of ca. 120, in the low-frequency band (3-9 Hz). Conversely, in the high frequency range (16-27 Hz) the average Q_p is significantly larger (ca. 640). The frequency decay method was also applied to S-waves spectra. In the low frequency range (2-5 Hz) the estimated average Q_s is ca. 190. As for Q_p, also Q_s, in the high frequency range (16-27 Hz), is larger (ca. 700). These results evidenced a frequency dependence of both the quality factors Q_p and Q_s, as commonly observed in tectonically active zones characterized by high degree of heterogeneity.     

Attenuation of High-Frequency Seismic Waves in Eastern Iran

We investigated the frequency-dependent attenuation of the crust in Eastern Iran by analysis data from 132 local earthquakes having focal depths in the range of 5–25 km. We estimated the quality factor of coda waves (Q _c) and body waves (Q _p and Q _s) in the frequency band of 1.5–24 Hz by applying the single backscattering theory of S-coda envelopes and the extended coda-normalization method, respectively. Considering records from recent earthquakes (Rigan M _w 6.5, 2010/12/20, Goharan M _w 6.2, 2013/5/11 and Sirch M _w 5.5, 2013/1/21), the estimated values of Q _c, Q _p and Q _s vary from 151 ± 49, 63 ± 6, and 93 ± 14 at 1.5 Hz to 1,994 ± 124, 945 ± 84 and 1,520 ± 123 at 24 Hz, respectively. The average frequency-dependent relationships (Q = Q _o f ⁿ ) estimated for the region are Q _c = (108 ± 10)f ^(0.96±0.01), Q _p = (50 ± 5)f ^(1.01±0.04), and Q _s = (75 ± 6)f ^(1.03±0.06). These results evidenced a frequency dependence of the quality factors Q _c, Q _p, and Q _s, as commonly observed in tectonically active zones characterized by a high degree of heterogeneity, and the low value of Q indicated an attenuative crust beneath the entire region.     

Study on relationship between features of strain anomaly obtained by borehole strain meter at Wushi station in Xinjiang and strong earthquake: Taking Jiashi <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>=6.8 earthquake as an example

The borehole strain meter at Wushi station recorded obvious anomaly before Jiashi M _s=6.8 earthquake occurred on February 24, 2003. Its features are as follows. (1) Anomaly types are complete. The trend anomaly, short-term anomaly, short-imminent anomaly and exponential anomaly appeared 19 months, 56 days, 4 days, and more than 1 month before the event, respectively; (2) Anomaly is large in magnitude. The maximal magnitude of strain anomaly is 1.7×10^?5, which is rare in the past 20-year observation records at Wushi station; (3) Strain rate fluctuates sharply with obvious alternation of tension and compression. According to the magnitude of strain anomaly, time of exponential anomaly appearance and regional features of recorded anomaly, we could predict the magnitude, occurrence time and potential region to a certain degree.     

Moisture origin and stable isotope characteristics of precipitation in southeast Siberia

The paper presents oxygen and hydrogen isotopes of 284 precipitation event samples systematically collected in Irkutsk, in the Baikal region (southeast Siberia), between June 2011 and April 2017. This is the first high-resolution dataset of stable isotopes of precipitation from this poorly studied region of continental Asia, which has a high potential for isotope-based palaeoclimate research. The dataset revealed distinct seasonal variations: relatively high δ¹⁸O (up to −4‰) and δD (up to −40‰) values characterize summer air masses, and lighter isotope composition (−41‰ for δ¹⁸O and −322‰ for δD) is characteristic of winter precipitation. Our results show that air temperature mainly affects the isotope composition of precipitation, and no significant correlations were obtained for precipitation amount and relative humidity. A new temperature dependence was established for weighted mean monthly precipitation: +0.50‰/°C (r² = 0.83; p <.01; n = 55) for δ¹⁸O and +3.8‰/°C (r² = 0.83, p < 0.01; n = 55) for δD. Secondary fractionation processes (e.g., contribution of recycled moisture) were identified mainly in summer from low d excess. Backward trajectories assessed with the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model indicate that precipitation with the lowest mean δ¹⁸O and δD values reaches Irkutsk in winter related to moisture transport from the Arctic. Precipitation originating from the west/southwest with the heaviest mean isotope composition reaches Irkutsk in summer, thus representing moisture transport across Eurasia. Generally, moisture transport from the west, that is, the Atlantic Ocean predominates throughout the year. A comparison of our new isotope dataset with simulation results using the European Centre/Hamburg version 5 (ECHAM5)-wiso climate model reveals a good agreement of variations in δ¹⁸O (r² = 0.87; p <.01; n = 55) and air temperature (r² = 0.99; p <.01; n = 71). However, the ECHAM5-wiso model fails to capture observed variations in d excess (r² = 0.14; p < 0.01; n = 55). This disagreement can be partly explained by a model deficit of capturing regional hydrological processes associated with secondary moisture supply in summer.     

Zirconium isotope evidence for incomplete admixing of r-process components in the solar nebula

Isotopic anomalies in Mo and Zr have recently been reported for bulk chondrites and iron meteorites and have been interpreted in terms of a primordial nucleosynthetic heterogeneity in the solar nebula. We report precise Zr isotopic measurements of carbonaceous, ordinary and enstatite chondrites, eucrites, mesosiderites and lunar rocks. All bulk rock samples yield isotopic compositions that are identical to the terrestrial standard within the analytical uncertainty. No anomalies in ⁹²Zr are found in any samples including high Nb/Zr eucrites and high and low Nb/Zr calcium-aluminum-rich inclusions (CAIs). These data are consistent with the most recent estimates of <10⁻⁴ for the initial ⁹²Nb/⁹³Nb of the solar system. There exists a trace of isotopic heterogeneity in the form of a small excess of r-process ⁹⁶Zr in some refractory CAIs and some metal-rich phases of Renazzo. A more striking enrichment in ⁹⁶Zr is found in acetic acid leachates of the Allende CV carbonaceous chondrite. These data indicate that the r- and s-process Zr components found in presolar grains were well mixed on a large scale prior to planetary accretion. However, some CAIs formed before mixing was complete, such that they were able to sample a population of r-process-enriched material. The maximum amount of additional r-process component that was added to the otherwise well-mixed Zr in the molecular cloud or disk corresponds to ∼0.01%.     

Cyclic resistance and liquefaction behavior of silt and sandy silt soils

The liquefaction behavior and cyclic resistance ratio (CRR) of reconstituted samples of non-plastic silt and sandy silts with 50% and 75% silt content are examined using constant-volume cyclic and monotonic ring shear tests along with bender element shear wave velocity (V_s) measurements. Liquefaction occurred at excess pore water pressure ratios (r_u) between 0.6 and 0.7 associated with cumulative cyclic shear strains (γ) of 4% to 7%, after which cyclic liquefaction ensued with very large shear strains and excess pore water pressure ratio (r_u>0.8). The cyclic ring shear tests demonstrate that cyclic resistance ratio of silt and sandy silts decreases with increasing void ratio, or with decreasing silt content at a certain void ratio. The results also show good agreement with those from cyclic direct simple shear tests on silts and sandy silts. A unique correlation is developed for estimating CRR of silts and sandy silts (with more than 50% silt content) from stress-normalized shear wave velocity measurements (V_s1) with negligible effect of silt content. The results indicate that the existing CRR–V_s1 correlations would underestimate the liquefaction resistance of silts and sandy silt soils.     

Estimation of Q p and Q s of Kinnaur Himalaya

The attenuation characteristics of the Kinnaur area of the North West Himalayas were studied using local earthquakes that occurred during 2008–2009. Most of the analyzed events are from the vicinity of the Panjal Thrust (PT) and South Tibetan Detachment Thrust, which are well-defined tectonic discontinuities in the Himalayas. The frequency-dependent attenuation of P and S waves was estimated using the extended coda normalization method. Data from 64 local earthquakes recorded at 10 broadband stations were used. The coda normalization of the spectral amplitudes of P and S waves was done at central frequencies of 1.5, 3, 6, 9, and 12 Hz. Q _p increases from about 58 at 1.5 Hz to 706 at 12 Hz, and Q _s increases from 105 at 1.5 Hz to 1,207 at 12 Hz. The results show that the quality factors for both P and S waves (Q _p and Q _s) increase as a function of frequency according to the relation Q?=?Q _o f ⁿ , where Q _o is the corresponding Q value at 1 Hz frequency and “n” is the frequency relation parameter. We obtained Q _p?=?(47?±?2)f ^(1.04±0.04) and Q _s?=?(86?±?4)f ^(0.96±0.03) by fitting power law dependency model for the estimated values of the entire study region. The Q ₀ and n values show that the region is seismically very active and the crust is highly heterogeneous. There was no systematic variation of values of Q _p and Q _s at different frequencies from one tectonic unit to another. As a consequence, average values of these parameters were obtained for each frequency for the entire region, and these were used for interpretation and for comparison with worldwide data. Q _p values lie within the range of values observed for some tectonically active regions of the world, whereas Q _s values were the lowest among the values compared for different parts of the world. Q _s/Q _p values were >1 for the entire range of frequencies studied. All these factors indicate that the crust is highly heterogeneous in the study region. The high Q _s/Q _p values also indicate that the region is partially saturated with fluids.     

Solar almucantar radiance and aerosol scattering: Comparative data from one tropical station and one mid-latitude station

Solar, sky, and aureole radiance data from Kensington, Australia, and Lagos, Nigeria, are analysed. The values of the aerosol optical thickness at Lagos are comparable to measured values elsewhere in the subregion (Abidjan) under light haze. The values at Kensington are close to those of the “background aerosol.” Identical and anomalous values of aerosol optical thickness occurred at 1.67 μm at both stations as less than 15% of the values can be accounted for by aerosols. The high optical depths are therefore explained in terms of absorption, possibly mainly by CO₂ (59–65%) and, to a lesser extent, by H₂O (20–26%). Following correction for molecular multiple scattering, the aerosol scattering phase function increases substantially at 0.5 μm in Kensington and decreases considerably at 0.76 μm in Lagos when compared to the uncorrected function, confirming the relative fineness of the Kensington aerosol. The Lagos aerosol size distribution for moderately dusty air is reasonably bimodally log-normal with modal radiusr _m=0.49 μm (1.36 μm) for the moder≤1 μm (r>1 μm). The Kensington aerosol resembles one withr _m=0.05 μm forr<1 μm. There the Born approximation appears appropriate, and the weighted mean square radiusp ²=0.13 μm². In the near forward direction the aerosol scattering phase funciton for Lagos resembles Deirmendjian’s simulations for a silicate-L haze and, for Kensington, a silicate-H haze.     

Studies of mantle structure of U.S.S.R. territory on long-range seismic profiles

Long-range seismic sounding carried out during the last few years on the territory of the U.S.S.R. has shown a basic inhomogeneity of the uppermost mantle, as well as evidence of regularities in the distribution of its seismic parameters. The following data were used: times and apparent velocities of P- and S-waves for investigation of mantle velocities, converted waves for seismic discontinuity model studies and wave attenuation for Q-factor estimation. Strong regularities were distinguished in the distribution of average seismic velocities for the uppermost mantle, in their dependence on the age and type of geostructure and on their position relative to the central part of the continent. Old platforms and the inner part of the continent are marked by velocities under the Mohorovi?i? discontinuity of more than 8.2–8.3 km s^?1, young platforms and outer parts of the continent by 8.0–8.2 km s^?1, and orogenic and rift zones by 7.8–8.0 km s^?1. The difference becomes more pronounced at a depth of about 100–200 km: for the old platform mantle velocities of 8.5–8.6 km s^?1 are typical; beneath the orogenic and rift areas, inversion zones with velocities less than 7.8 km s^?1 are observed.The converted waves show fine inhomogeneities of the crust and uppermost mantle, the presence of many discontinuities with positive and negative changes of velocity, and anisotropy of seismic waves in some of the layers. Wave attenuation allowed the determination of the Q-factor in the mantle. It varied from one region to another but a close relation between Q and P-wave velocity is the main cause of its variation.     

Pan evaporation modelling and changing attribution analysis on the Tibetan Plateau (1970–2012)

Pan evaporation (E_p) is an important indicator of water and energy and the decline of E_p has been reported in many regions over the last decades. The climate and E_p are dependent on each other. In this study, the temporal trends of E_p and main E_p drivers, namely mean air temperature (T_a), wind speed (u), global solar radiation (R_s), net long‐wave radiation(R_nl) and vapour pressure deficit (D) from 1970 to 2012, were calculated on the basis of 26 meteorological stations on the Tibetan Plateau. The arithmetic average of E_p from 26 stations decreased with the rate of ?11.91 mm a^?2; the trends of R_s, R_nl, T_a, u and D were ?1.434 w m^?2 decade^?1, 0.2511 w m^?2 decade^?1, 0.3590°C decade^?1, ?0.2376 m s^?1 decade^?1 and 9.523 Pa decade^?1, respectively. The diffuse irradiance is an essential parameter to model E_p and quantify the contribution of climatic factors to changing E_p. 60 724 observations of R_s and diffuse solar irradiance (R_d) from seven of the 26 stations were used to develop the correlation between the diffuse fraction (R_d/R_s), and the clearness index (R_s/R_o). On the basis of the estimation of the diffuse component of R_s and climatic data, we modified the PenPan model to estimate Chinese micro‐pan evaporation (E_p) and assess the attribution of E_p dynamics using partial derivatives. The results showed that there was a good agreement between the observed and calculated daily E_p values. The observed decrease in E_p was mostly due to declining wind speed (?13.7 mm a^?2) with some contributions from decreasing solar irradiance (?3.1 mm a^?2); and the increase of temperature had a large positive effect (4.55 mm a^?2) in total whilst the increase of R_nl had insignificant effect (0.35 mm a^?2) on E_p rates. The change of E_p is the net result of all the climatic variables. Copyright © 2014 John Wiley & Sons, Ltd.     

The geophysical significance of oceanic plagiogranite

Rocks ranging in composition from trondhjemite to diorite (plagiogranites) have been recovered from ocean ridges and are common constituents of ophiolites. Velocities and densities of diorite and trondhjemite from the Mid-Atlantic Ridge are shown to differ significantly from similar properties of metadolerite and gabbro. Compressional (V_p) and shear (V_s) velocities of plagiogranites are relatively low (V_p = 4.78–5.91km/s at1kbar,V_s = 2.81–3.37km/s at1kbar), as are densities (2.57–2.64 g/cm³) and Poisson's ratios (0.24–0.27). These data lend strong support to the probable existence of a low-velocity/density zone within layer 3 of the oceanic crust. Based on observations in ophiolites, it is postulated that this zone can be up to 1 km in thickness and is laterally discontinuous.     

Estimation of Strong Ground Motion from a Great Earthquake Mw 8.5 in Central Seismic Gap Region, Himalaya (India) Using Empirical Green’s Function Technique

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

Variation of intrinsic and scattering attenuation of seismic waves with depth in the Bam region,East-Central Iran

The attenuation properties of the lithosphere in the Bam region, East-Central Iran, have been investigated. For this purpose, 42 local earthquakes having focal depths less than 25 km have been used. The quality factor of coda waves (Q_c) has been estimated using the single back-scattering model. The quality factors Q_p, Q_d (P and direct S-waves) have been estimated using the extended-coda normalization method. Q_i and Q_s (the intrinsic and scattering attenuation parameters) have been estimated for the region. The values of Q_p, Q_d, Q_c, Q_i and Q_s show a dependence on frequency in the range of 1.5–24 Hz for the Bam region. The average frequency-dependent relationships estimated for the region are Q_p=(36±6)f^(1.03±0.06), Q_d=(59±8)f^(1.00±0.03), Q_c=(79±5)f^(1.01±0.04), Q_s=(131±4)f^(1.01±0.04) and Q_i=(104±6)f^(1.01±0.05). A comparison between Q_i and Q_s shows that intrinsic absorption is predominant over scattering.The variation of Q has also been estimated at different lapse times to observe heterogeneities variation with depth. The variation of Q with frequency and lapse time shows that the lithosphere becomes more homogeneous with depth.The estimated Q_o values at different stations suggest a low value of Q indicating a heterogeneous and attenuative crust beneath the entire region.     

A semi-analytical model of expected areal-average infiltration under spatial heterogeneity of rainfall and soil saturated hydraulic conductivity

A semi-analytical model for the estimate of expected areal-average infiltration rate at hillslope scale is presented. It accounts for spatial heterogeneity of the saturated hydraulic conductivity, K_s, and rainfall rate, r. The K_s field is characterized by a lognormal probability density function while the rainfall rate r is represented by a uniform distribution between two extreme values. The model formulation relies upon the use of cumulative infiltration as the independent variable which is then expressed as a function of an expected time for use in practical applications. The solution is applicable for those ranges of r and K_s that allow for neglecting the infiltration of surface water running downslope into pervious soils (run-on process). The model was tested by comparisons with Monte Carlo simulations carried out for a variety of coefficients of variation of r and K_s over a clay loam soil and a sandy loam soil. The model was found to be very reliable both with coupled spatial variability of r and K_s and when only one variable is characterized by spatial heterogeneity while the other is uniform.