Wave propagation in a seven-story reinforced concrete building: III. Damage detection via changes in wavenumbers

This paper presents low frequency wavenumbers in a seven-storey reinforced concrete building estimated from its recorded response to eleven earthquakes, one of which (1994 Northridge) caused visible structural damage, and two of which are its aftershocks. The wavenumbers, K_i,j(f), are estimated from pairs (i,j) of records at neighboring recording sites in the building, distributed vertically or horizontally. Changes in K_i,j(f) from one event to another are compared in the undamaged (lower) and in the damaged (upper) part of the building, with the aim to find whether trends in K_i,j(f) can indicate damage. The results suggest significant and permanent increase of the wavenumbers in the damaged parts for the 1994 Northridge earthquake and its aftershocks, which is not the case for the other events in the damaged parts, and for all eleven events in the undamaged parts of the building. This increase in wavenumbers in the damaged parts can be explained by reduced wave velocities through the damaged structural members, and by scattering of waves from the discontinuities created by the damage. It is concluded from this qualitative analysis that wavenumbers estimated from strong motion recordings in a building can indicate location of damage, and that it would be useful to refine further this method (extend it to higher frequencies, and add the capability to quantify the damage). However, this would require more dense strong motion instrumentation in buildings than currently available. Deployment of dense arrays in selected buildings would provide data for further work on this subject.     

Diatoms from the genus Achnanthidium in flowing waters of the Appalachian Mountains (North America): Ecology, distribution and taxonomic notes

Diatoms from the genus Achnanthidium are abundant in rivers, streams, and springs of the Appalachian Mountains. They inhabit clean and polluted waters, including those affected by acid mine drainage. The identification of Achnanthidium taxa is difficult due to their small cell size and insufficient information in the diatom floras. We studied the taxonomy and ecology of Achnanthidium in Appalachian rivers by analyzing a data set of benthic diatom samples and corresponding water chemistry data collected during several water-quality surveys from 181 sampling sites. Ten species were identified using scanning electron and light microscopy: A. alpestre (Lowe & Kociolek) Lowe & Kociolek, A. atomus (Hustedt) Monnier, Lange-Bertalot, & Ector, A. deflexum (Reimer) Kingston, A. duthii (Sreenivasa) Edlund, A. eutrophilum (Lange-Bertalot) Lange-Bertalot, A. cf. gracillimum (Meister) Lange-Bertalot, A. cf. latecephalum Kobayasi, A. minutissimum (Kützing) Czarnecki (sensu lato), A. reimeri (Camburn) comb. nov., and A. rivulare Potapova & Ponader. The distribution of common taxa in relation to water chemistry was studied by fitting non-parametric regression models (generalized additive models, GAM, and non-parametric multiplicative regression models, NPMR) to species relative abundances. Studied Achnanthidium species differed considerably in their responses to water chemistry. These results suggest that species-level identifications will lead to more accurate bioassessments.     

Interception losses in three non‐native Hawaiian forest stands

Interception losses in stands of non‐native trees in Hawaiian forests and their potential negative impacts on fresh water availability are poorly understood. In this study, a canopy water balance analysis was conducted to estimate interception losses using measurements of rainfall (RF), throughfall (TF), and stemflow (SF) at three locations, each dominated by one or more of the following non‐native tree species: Psidium cattleianum Sabine (Strawberry guava), Schinus terebinthifolius Raddi (Christmas berry), Syzygium cumini (L.) Skeels (Java plum), and Coffea arabica L. (Coffee). Mean TF expressed as percentage of total RF was the lowest (43.3%) under a monotypic stand of P. cattleianum and the highest (56.5%) under mixture of S. terebinthifolius, P. cattleianum, and S. cumini. Observed SF was highest (33.9%) under P. cattleianum and lowest (3.6%) under a mixture of S. terebinthifolius, P. cattleianum, and S. cumini. The relatively high SF under P. cattleianum can be attributed to its smooth bark, stem density, and steep branching. The mean observed canopy interception varied between 23% under P. cattleianum and 45% at the site dominated by C. arabica. Mean direct TF coefficients from individual events at each location ranged from a low of 0.36 under the canopy dominated by C. arabica to a high of 0.51 under the canopy dominated by S. terebinthifolius, P. cattleianum, and S. cumini. In contrast, the mean SF partitioning coefficients from individual storm events at each location ranged from a low of 0.05 under the canopy dominated by S. terebinthifolius, P. cattleianum, and S. cumini to a high of 0.37 under P. cattleianum. Mean canopy storage capacity was highest (1.90) at the site dominated by S. terebinthifolius, P. cattleianum, and S. cumini whereas trunk storage capacity was highest (0.54) under the P. cattleianum. Copyright © 2012 John Wiley & Sons, Ltd.     

Effect of soil stiffness in the attenuation of Rayleigh-wave motions from field measurements

The amplitude of vertical ground surface vibrations generated by impact tests on the ground surface was measured at various radial distances from the point of impact at locations of Greece. The results of measurements were analyzed in the frequency domain (in the range from 0–100 Hz) and the attenuation characteristics of soil materials were studied in terms of a frequency-independent attenuation coefficient, a₀, of the empirical Bornitz equation. The aim of the study was to investigate the effect of soil stiffness (expressed by the value of low-amplitude shear wave velocity of soil, V_SO) on the value of attenuation coefficient, a₀. Values of V_SO for the tested soils were estimated by applying the methodology of Spectral Analysis of Surface Waves (SASW) technique and utilizing the surface vibration data. An empirical relationship between a₀ and V_SO1 (V_SO1 is the representative value of V_SO for the soil profile up to a depth of one wavelength) was established for values of V_SO1 ranging from 140 to 1000 m/s. A similar relationship in terms of the low-amplitude shear modulus of soil, G_O1, was also established by converting the V_SO1 values to G_O1 values. The experimental results were compared to values reported in the literature for comparable soil types and frequencies of vibration and a reasonable agreement was found to exist. The proposed empirical relationship can be utilized in many practical applications of soil dynamics requiring the knowledge of the attenuation rate of Rayleigh waves with distance in various types of soils.     

Seasonal climate patterns and their influence on calibration of the Hargreaves-Samani equation

Abstract

Annual patterns in climate parameters were studied to evaluate how these influence the quality of reference evapotranspiration (ETo) estimates obtained from the Hargreaves-Samani (HS) equation, since the method only uses the measured temperature directly. The work evaluates how these patterns can be used to improve the HS ETo estimates. Ten-year moving averages from a set of California Irrigation Management Information System (CIMIS) stations were used to evaluate the relationships between solar radiation (R_s), temperature (T) and ETo. The results indicate that T treads behind solar radiation and its value peaks some 25 days later. Thus, the main irrigation season in the Mediterranean climate (1 May–30 September) can be divided into three phases: increasing R_s and T; decreasing R_s with increasing T; and decreasing R_s and T. Non-univocal annual cycles were observed between R_s and T, ETo and R_s, and ETo and T. These annual patterns result in important seasonal changes in the ratio between the HS and Penman-Monteith (FAO PM) ETo estimates. The changes are particularly important during the irrigation season, where the FAO PM initially calculates greater ETo values than the HS methodology, and from the end of May to early September, where the HS equation overestimates the ETo values (by 17 mm, or 3%). These patterns obtained from 2000–2009 data were used to calibrate and improve HS ETo estimates at new sites for the 2010–2011 period. Calibration based on the proposed seasonal region-wide FAO PM/HS ETo ratios improved both the bias (decreased from 0.40 to 0.36 mm d^-1) and r² (increased from 0.67 to 0.87) of the ETo estimates for the irrigation season. The proposed methodology can be easily applied to other regions, even when the existing weather stations are sparse.

Editor Z.W. Kundzewicz     

A simple two‐parameter model for scaling hillslope surface runoff

The objective of this research was to develop and parameterise a physically justified yet low‐parameter model to quantify observed changes in surface runoff ratios with hillslope length. The approach starts with the assumption that a unit of rainfall‐excess runoff generated at a point is a fraction β of precipitation P (m) which travels some variable distance down a slope before reinfiltrating, depending on the local rainfall, climate, soils, etc. If this random distance travelled Y is represented by a distribution, then a survival function will describe the probability of this unit of runoff travelling further than some distance x (m). The total amount of per unit width runoff Q (m²) flowing across the lower boundary of a slope of length λ (m) may be considered the sum of all the proportions of the units of rainfall excess runoff integrated from the lower boundary x = 0 to the upper boundary x = λ of the slope. Using these assumptions we derive a model Q(λ) = βPμλ/(μ + λ), (μ > 0, 0 ≤ β ≤ 1, λ ≥ 0) that describes the change in surface runoff with slope length, where μ (m) is the mean of the random variable Y. Dividing both sides of this equation by Pλ yields a simple two‐parameter equation for the dimensionless hillslope runoff ratio Q_h(λ) = βμ/(μ + λ). The model was parameterised with new rainfall and runoff data collected from three replicates of bounded 2 m wide plots of four different lengths (0.5, 1.0, 2.0 and 4.0 m) for 2 years from a forested SE Australian site, and with 32 slope length–runoff data sets from 12 other published studies undertaken between 1934 and 2010. Using the parameterised model resulted in a Nash and Sutcliffe statistic between observed and predicted runoff ratio (for all data sets combined) of 0.93, compared with –2.1 when the runoff ratio was fixed at the value measured from the shortest plot. Copyright © 2014 John Wiley & Sons, Ltd.     

Converted-wave seismology in anisotropic media Revisited,Part II: Application to parameter estimation

In transversely isotropic media with a vertical symmetry axis (VTI), the converted-wave (C-wave) moveout over intermediate-to-far offsets is determined by four parameters. These are the C-wave stacking velocity V _C2, the vertical and effective velocity ratios γ ₀and γ _eff, and the anisotropic parameter X _eff. We refer to the four parameters as the C-wave stacking velocity model. The purpose of C-wave velocity analysis is to determine this stacking velocity model. The C-wave stacking velocity model V _C2, γ ₀, γ _geff, and X _eff can be determined from P- and C-wave reflection moveout data. However, error propagation is a severe problem in C-wave reflection-moveout inversion. The current short-spread stacking velocity as deduced from hyperbolic moveout does not provide sufficient accuracy to yield meaningful inverted values for the anisotropic parameters. The non-hyperbolic moveout over intermediate-offsets (x/z from 1.0 to 1.5) is no longer negligible and can be quantified using a background γ. Non-hyperbolic analysis with a γ correction over the intermediate offsets can yield V _C2 with errors less than 1% for noise free data. The procedure is very robust, allowing initial guesses of γ with up to 20% errors. It is also applicable for vertically inhomogeneous anisotropic media. This improved accuracy makes it possible to estimate anisotropic parameters using 4C seismic data. Two practical work flows are presented for this purpose: the double-scanning flow and the single-scanning flow. Applications to synthetic and real data show that the two flows yield results with similar accuracy but the single-scanning flow is more efficient than the double-scanning flow. This work is funded by the Edinburgh Anisotropy Project of the British Geological Survey. First Author Li Xiangyang, he is currently a professorial research seismologist (Grade 6) and technical director of the Edinburgh Anisotropy Project in the British Geological Survey. He also holds a honorary professorship multicomponent seismology at the School of Geosciences, University of Edinburgh. He received his BSc(1982) in Geophysics from Changchun Geological Institute, China, an MSc (1984) in applied geophysics from East China Petroleum Institute (now known as the China University of Petroleum), and a PhD (1992) in seismology from the University of Edinburgh. During 1984–1987, he worked as a lecturer with the East China Petroleum Institute. Since 1991, he has been employed by the British Geological Survey. His research interests include seismic anisotropy and multicomponent seismology.     

Hydrogen isotope composition of natural gases from the Tarim Basin and its indication of depositional environments of the source rocks

By measuring carbon and hydrogen isotope compositions for C1, C2 and C3 of 74 gas samples, natural gases from the Tarim Basin can be divided into six groups on the basis of their origins: (1) coal-type gas derived from coal measures; (2) coal-type gas generated from the T-J lacustrine mudstones; (3) oil-type gas derived from the Cambrian and low Ordovician marine source rocks; (4) oil-type gas from the source rocks deposited in the marine-transitional facies; (5) mixing gas between gas derived from the Carboniferous transitional source rocks and the Mesozoic humic gas, and (6) mixing gases of thermal genetic gas and little deep gas in the Southwest depression of the Tarim Basin. The δ D values of methane in natural gases originating from different type kerogens are affected by both palaeo-environments of the source rock formation (kerogen types) and thermal maturity, with sedimentary environment (kerogen type) as the main controlling factor. Under the similar thermal maturity, the hydrogen isotope composition of methane is more enriched in deuterium in marine environments than lacustrine one. With the increase of thermal maturity and the increase of carbon atomic numbers of gaseous alkanes, the hydrogen isotopes become enriched in deuterium. The δ D values of ethane and propane (δ D2, δ D3) are controlled mainly by thermal maturity and to a lesser degree by sedimentary environment of the source rock formation. The partial reversal of hydrogen isotopes for gaseous alkanes would be related to the microbial oxidation, mixing of sapropelic and humic gases and / or mixing of gases from similar kerogen sources with various thermal maturities. In the oil-type gas, the sulfate reduction reaction would result in the reversed order of δ D1 and δ D2 (e.g. δ D1>δ D2).     

Lg coda Q0 value and its relation with the tectonics in Chinese mainland and adjacent regions

We have collected 432 vertical component records from 45 stations of new CENC （China Earthquake Network Center） in Chinese mainland and adjacent regions. These records were used to calculate Q0 （Q at 1Hz） and η values of Lg coda from each station by the stack spectral ratio （SSR） method. Then the tomography method was applied to obtaining lateral variation of Q0 and η values in Chinese mainland and adjacent regions. The result indicates that Q0 value varies between 150 and 600 in the studied areas. Yunnan, southwest Sichuan, and northwest Myanmar show the lowest Q0 value （Q0〈240） and the crust of these regions is characterized by complicated crack and strong hydrothermal activity. The highest Q0 value （Qo〉510） exists in the border of southern Mongolia, Alxa and Ordos block. The η value varies between 0.45 and 0.75 in Chinese mainland and its adjacent regions.     

Modeling gross primary production of a temperate grassland ecosystem in Inner Mongolia,China, using MODIS imagery and climate data

Carbon fluxes in temperate grassland ecosystems are characterized by large inter-annual variations due to fluctuations in precipitation and land water availability. Since an eddy flux tower has been in operation in the Xilin Gol grassland, which belongs to typical temperate grassland in North China, in this study, observed eddy covariance flux data were used to critically evaluate the biophysical per- formance of different remote sensing vegetation indices in relation to carbon fluxes. Furthermore, vegetation photosynthesis model (VPM) was introduced to estimate gross primary production (GPP) of the grassland ecosystem for assessing its dependability. As defined by the input variables of VPM, Moderate Resolution Imaging Spectroradimeter (MODIS) and standard data product MOD09A1 were downloaded for calculating enhanced vegetation index (EVI) and land surface water index (LSWI). Measured air temperature (Ta) and photosynthetically active radiation (PAR) data were also included for model simulating. Field CO2 flux data, during the period from May, 2003 to September, 2005, were used to estimate the "observed" GPP (GPPobs) for validation. The seasonal dynamics of GPP predicted from VPM (GPPVPM) was compared quite well (R2=0.903, N=111, p<0.0001) with the observed GPP. The ag- gregate GPPVPM for the study period was 641.5 g C·m?2, representing a ~6% over-estimation, compared with GPPobs. Additionally, GPP predicted from other two typical production efficiency model (PEM) represents either higher overestimation or lower underestimation to GPPobs. Results of this study demonstrate that VPM has potential for estimating site-level or regional grassland GPP, and might be an effective tool for scaling-up carbon fluxes.     

Determining Q of near-surface materials from Rayleigh waves

High-frequency (≥2 Hz) Rayleigh wave phase velocities can be inverted to shear (S)-wave velocities for a layered earth model up to 30 m below the ground surface in many settings. Given S-wave velocity (V_S), compressional (P)-wave velocity (V_P), and Rayleigh wave phase velocities, it is feasible to solve for P-wave quality factor Q_P and S-wave quality factor Q_S in a layered earth model by inverting Rayleigh wave attenuation coefficients. Model results demonstrate the plausibility of inverting Q_S from Rayleigh wave attenuation coefficients. Contributions to the Rayleigh wave attenuation coefficients from Q_P cannot be ignored when Vs/V_P reaches 0.45, which is not uncommon in near-surface settings. It is possible to invert Q_P from Rayleigh wave attenuation coefficients in some geological setting, a concept that differs from the common perception that Rayleigh wave attenuation coefficients are always far less sensitive to Q_P than to Q_S. Sixty-channel surface wave data were acquired in an Arizona desert. For a 10-layer model with a thickness of over 20 m, the data were first inverted to obtain S-wave velocities by the multichannel analysis of surface waves (MASW) method and then quality factors were determined by inverting attenuation coefficients.     

Diagnosing vertical motion in the Equatorial Atlantic

Estimating the vertical velocity (w) in the oceanic upper-layers is a key issue for understanding the cold tongue development in the Eastern Equatorial Atlantic. In this methodological paper, we develop an expanded and general formulation of the vertical velocity equation based on the primitive equation (PE) system, in order to gain new insight into the physical processes responsible for the Equatorial and Angola upwellings. This approach is more accurate for describing the real ocean than simpler considerations based on just the wind-driven patterns of surface layer divergence. The w-sources/forcings are derived from the PE w-equation and diagnosed from a realistic ocean simulation of the Equatorial Atlantic. Sources of w are numerous and express the high complexity of terms related to the turbulent momentum flux, to the circulation and to the mass fields, some of them depending explicitly on w and others not. The equatorial upwelling is found to be mainly induced by the (i) the zonal turbulent momentum flux, (ii) the curl of turbulent momentum flux and (iii) the imbalance between the circulation and the pressure fields. The Angola upwelling in the eastern part of the basin is controlled by strong curl of turbulent momentum flux. A strong cross-regulation is evidenced between the w-forcings independent of w and dependent on w, which suggests an equatorial balanced-dynamics. The w-forcing depending on w represents the negative feedback of the ocean to the w-forcing independent of w: in the equatorial band, this adjustment is led by non-linear processes and by vortex stretching outside.     

Changes in the yield of chlorophyll a from dissolved available inorganic nitrogen after an enrichment event—applications for predicting eutrophication in coastal waters

An understanding of the dynamic relationship between nitrogen supply and the formation of phytoplankton biomass is important in predicting and avoiding marine eutrophication. This relationship can be expressed as the short-term yield q of chlorophyll from dissolved available inorganic nitrogen (DAIN), the sum of nitrate, nitrite and ammonium. This paper communicates the results of a continuous culture nitrate enrichment experiment undertaken to investigate the cumulative yield of chlorophyll from DAIN (q). The purposes of the study were: to acquire a better understanding of the relationship between chlorophyll formation and DAIN; to obtain values that could be used in models for predicting eutrophication. The results of a time series experiment carried out using microplankton (all organisms <200 μm in size) indicate that the parameter q does not have a single value but is affected by the ecophysiological response of phytoplankton to changing nutrient status after an enrichment event. It is also dependent on changes in the allocation of nitrogen between autotrophs and heterotrophs. The value of yield obtained at the height of the bloom can be represented by q (max) (2.35 μg chl (μmol N)⁻¹). The post-bloom, steady state value of q can be represented by q_eq (0.95 μg chl (μmol N)⁻¹). The microcosm steady state yield was not significantly different from the median value obtained from synoptic studies of Scottish west coast waters. It is proposed that q_eq is the most appropriate value for assessing the general potential for eutrophication resulting from continuous nutrient enrichment into coastal waters. It is further proposed that q (max) be used for cases of sporadic enrichment and where a short burst of unrestricted growth may be detrimental.     

The evolution from Miocene potassic to Quaternary sodic magmatism in western Turkey: implications for enrichment processes in the lithospheric mantle

New K-Ar dating and major- and trace-element analyses from the U ak-Selendi-Emet (USE) area constrain the timing of changes in the nature of volcanism in the Miocene in western Turkey. The data reveal a change from dominantly calc-alkaline and silicic in the Early Miocene to largely alkaline and more mafic volcanism in the Middle Miocene. This probably reflects a decreasing amount of crustal contamination with time, a result of extensional tectonics. High levels of various incompatible elements (including K) in the more mafic members, suggest an enriched subcontinental lithospheric source region for the Middle Miocene USE lavas. Highly variable Nb/Y, Ti/Y and Th/Nb ratios suggest a lithospheric mantle heterogeneously enriched by two processes: (1) enrichment by subduction-related processes producing high Th/Nb but low Nb/Y and Ti/Y; and (2) enrichment by small degree melts of depleted upper mantle producing low Th/Nb but high Nb/Y and Ti/Y. Both of these enrichment processes have variably contributed to Middle Miocene K-rich lavas in the USE area. The mechanism which initiated the melting of the enriched lithosphere is considered to be extension which produced decompression melting. Comparisons with the nearby Kula lavas reveals that by the Pliocene to Quaternary, volcanism, although still enriched in incompatible elements, had become sodic. It seems likely that continued extension up to this time thinned the lithosphere to such an extent that asthenospheric melts were produced which ascended and mixed with previously enriched lithosphere.     

Comparison of techniques that use the single scattering model to compute the quality factorQ from coda waves

Using simulated coda waves, the resolution of the single-scattering model to extract codaQ (Q _c ) and its power law frequency dependence was tested. The back-scattering model ofAki andChouet (1975) and the single isotropic-scattering model ofSato (1977) were examined. The results indicate that: (1) The inputQ _c models are reasonably well approximated by the two methods; (2) almost equalQ _c values are recovered when the techniques sample the same coda windows; (3) lowQ _c models are well estimated in the frequency domain from the early and late part of the coda; and (4) models with highQ _c values are more accurately extracted from late code measurements.     

再论各向异性介质转换波地震学,第二部分:参数计算研究

在具有垂直对称轴横向各向同性介质中，利用四种参数来确定中间至远偏移距转换波（C-波）动校正。它们是C-波叠加速度Vc2，垂直速度比和有效速度比γ0和γeff以及各向异性参数χeff。我们将这四种参数作为C波叠加速度模型。C-波速度分析的目的就是确定这种叠加速度模型。C-波叠加速度模型Vc2，γ0，γeff，和χeff可以由P-波和C-波反射动校正资料获得。然而错误的传播是C-波反射动校正反演中的严重问题。当前短排列叠加速度由于是从双曲线动校正推算而得，因而其精度不足以为各向异性参数提供有意义的反演值。中间偏移非双曲线动校正不再被人们所勿略，而是可以用一个背景γ加以量化。非双曲线分析通过中间偏移距的γ校正量可以产生Vc2，若数据不含燥音，其误差小于1％。方法稳健，允许γ启始假定值的误差达20％。该方法也适用垂直非均匀各向异性介质。精度的提高使能够用4分量地震资料计算各向异性参数。为此提出了两种工作流程：双扫描和单扫描流程。理论数据和实际数据的应用表明这两种流程得出的结果其精度相似，但是单扫描流程比双扫描更有效。     

Eddy correlations for water flow in a single fracture with abruptly changing aperture

Fluid flow in single fractures with non‐uniform apertures is an important research subject in many disciplines. The abruptly changing aperture is a special case of such non‐uniformity. This paper simulates water flow in a single fracture with abruptly changing aperture (SF‐ACA) using the Lattice Boltzmann Method (LBM) and the Finite Volume Method (FVM). The flow occurs with the Reynolds number (Re) ranging from 5 to 900 and a ratio of aperture change (E) of 3 (E = D/d, where D and d are the larger and smaller apertures, respectively). For Re values between 5 and 100, both LBM and FVM can successfully simulate the eddy development in the expansion regime of an SF‐ACA. Flow with high Re values (up to 900) is simulated by FVM, which appears to be numerically more stable than LBM for high‐Re flow problems studied here. The flow symmetry in the expansion regime breaks at the Re value between 400 and 500. Our simulation result shows a linear relationship between l₁/d and Re at low Re (5–100) or higher Re (110–900) values, where defined as the length from the location of abrupt expansion to the right edge of the first eddy along the flow direction. If considering the simulation results for the entire simulated range of Re (5–900), the l₁/d–Re relationship is better described by a non‐linear logarithmical function. The l₁/d approaches an asymptotic constant at large Re. Copyright © 2011 John Wiley & Sons, Ltd.     

Reliability of rayleigh wave dispersion curve obtained from f–k spectral analysis of microtremor array measurement

The paper utilizes previously developed microtremor simulation technique to evaluate the reliability of Rayleigh wave dispersion curve estimated by f–k spectral analysis of microtremor array measurement. The simulated microtremors are obtained for a fictitious (TEST) site. Attempt is also made to obtain the dispersion curve for two real sites (OHDATE and SKC) by inverse analysis of the microtremor array measurement using f–k spectra method. The estimated dispersion curve from simulated microtremors (TEST site) compares well with the theoretical dispersion curve, demonstrating the reliability of f–k spectra method and indicating that the estimated dispersion curve from microtremor measurements could be adequately used as the target for inverse analysis purposes. It is also demonstrated that the dispersion curve from microtremor measurements can be utilized to estimate the soil profile at OHDATE and SKC sites by inverse analysis. Results show that the theoretical dispersion curve of the fundamental mode of Rayleigh wave after the end of inverse analysis are in good agreement with the dispersion curve obtained by f–k spectral analysis of microtremor array measurement.     

On the ecology of Azolla filiculoides Lam. in Damietta district, Egypt

Azolla filiculoides Lam. is a tiny aquatic fern growing as free floating mats covering the stagnant water of the irrigation and drainage canals of Damietta district, Egypt. Ecological characteristics and distribution of the fern were investigated and results were discussed.Vegetation analysis of Azolla indicate that the most common associated species with Azolla are: Phragmites australis [presence (P) = 93.3%], Echinochloa stagnina (P = 70%), Eichhornia crassipes (P = 57%), Pluchea dioscoridis (P = 56.7%), Persicaria salicifolia (P = 26.7%), Phylla nodiflora (P = 13.3%) and Ludwigia stolonifera (P = 10%). Submerged species, e.g. Ceratophyllum spp. and Potamogeton spp. are completely absent in the stands dominated with solid mats of Azolla. This could be explained by that these mats prevent light penetration and oxygen which are required to submerged plants which in turn are decayed.The pH of the water where Azolla was recorded ranged from 7.1 to 9.0 (mean = 7.8), total dissolved salts ranged from 280 to 2600 mg/l (mean = 880 mg/l). The total phosphorus showed the range from 0.01 to 0.82 mg/l (mean = 0.22 mg/l). The ammonia-N is relatively high, varying from 2 to 18 mg/l (mean = 7.0 mg/l), however, the nitrate-N varied from 2.1 to 54.1 mg/l with a mean value of 20.5 mg/l.The fresh biomass of Azolla ranged from 114 to 4280 g/m² in irrigation canals and from 506 to 2760 g/m² in drainage canals. The dry biomass varied from 16 to 23 g/m² in irrigation canals. The range from 26 to 320 g/m² was obtained for drainage canals.     

Nd isotopic compositions of the Tethyan Himalayan Sequence in southeastern Tibet

The Himalayan orogen consists of three major lithologic units that are separated by two major north-dipping faults: the Lesser Himalayan Sequence (LHS) below the Main Central Thrust (MCT), the Greater Himalayan Crystalline Complex (GHC) above the MCT, and the Tethyan Himalayan Sequence (THS) juxtaposed by the South Tibet Detachment fault (STD) over the GHC. Due to widespread meta-morphism and intense deformation, differentiating the above three lithologic units is often difficult. This problem has been overcome by the use of Sm-Nd isotopic analysis. The previous studies suggested that the LHS can be clearly distinguished from the GHC and THS by their Nd isotope compositions. However, the lack of detailed and systematic Sm-Nd isotopic studies of the THS across the Himalaya in general has made differentiation of this unit from the nearby GHC impossible, as the two appear to share overlapping Nd compositions and model ages. To address this problem, we systematically sam-pled and analyzed Nd isotopes of the THS in southeastern Tibet directly north of Bhutan. Our study identifies two distinctive fields in a εNd -TDM plot. The first is defined by the εNd(210 Ma) values of -3.45 to -7.34 and TDM values of 1.15 to 1.29 Ga from a Late Triassic turbidite sequence, which are broadly similar to those obtained from the Lhasa block. The second field is derived from the Early Cretaceous meta-sedimentary rocks with εNd(130 Ma) values from -15.24 to -16.61 and TDM values from 1.63 to 2.00 Ga; these values are similar to those obtained from the Greater Himalayan Crystalline Complex in Bhutan directly south of our sampling traverse, which has εNd(130 Ma) values of -10.89 to -16.32 and Nd model ages (TDM) of 1.73 to 2.20 Ga. From the above observations, we suggest that the Late Triassic strata of the southeast Tibetan THS were derived from the Lhasa block in the north, while the Early Cretaceous strata of the THS were derived from a source similar to the High Himalayan Crystalline Complex or Indian craton in the south. Our interpretation is consistent with the existing palaeocurrent data and provenance analysis of the Late Triassic strata in southeastern Tibet, which indicate the sediments derived from a northern source. Thus, we suggest that the Lhasa terrane and the Indian craton were close to one another in the Late Triassic and were separated by a rift valley across which a large submarine fan was transported southward and deposited on the future northern margin of the Indian continent.