Boundary integral equation method for the diffraction of elastic waves using simplified green's functions

An alternate formulation of the ‘substructure deletion method’ suggested by Dasgupta in 1979¹ has been successfully implemented. The idea is to utilize simple Green's functions developed for a surface problem to replace the more complicated Green's functions required for embedded problems while still being able to generate an accurate solution. Since the exterior medium is usually represented by a continuum model, the interior medium in the present approach will also be represented by a continuum model rather than a finite element model as suggested originally, thereby eliminating the incompatibility between the solutions of the interior and exterior media. Detailed studies of the method's accuracy and limitations were performed using two-dimensional examples in wave scattering of canyons and alluvial valleys, problems which are more suitable for this method than the embedded foundation problem. The results obtained indicate that the alternate formulation gives accurate results only when the vertical dimension of the scattering object is not too large; if the aspect ratio (vertical over lateral) exceeds a certain limit, the results will not approach the known results given by boundary integral equation solutions or indirect boundary integral equations no matter what the refinement of the model may be. The greatest advantage of the present method is that the task of calculating Green's functions is reduced significantly; computational time using this new formulation is approximately five times less than for conventional boundary integral equation methods.     

Consistency of the spatial autocorrelation method with seismic interferometry and its consequence

We have cross‐checked the conventional theory of the spatial autocorrelation method and the consequence of seismic interferometry: the retrieval of the elastodynamic Green's function. Their mutual consistency is almost complete. The basic formulas of the conventional spatial autocorrelation theory can be derived by an alternative approach based on the retrieval of the elastodynamic Green's function. The only discrepancy is found with the average of the complex coherence function over azimuth in a wavefield dependent on azimuth. It is hypothesized, in discussion, that this discrepancy is due to the way of representing the wavefield in the background theory of seismic interferometry that can produce only wavefields moderately dependent on azimuth and that the mentioned consequence of seismic interferometry can also only make sense in a wavefield moderately dependent on azimuth. Our field experiment with a wavefield dependent on azimuth showed that the consequence of seismic interferometry in the logical framework of the conventional spatial autocorrelation theory is appropriate under such degrees of approximation as the measure proposed in this study, i.e., the deviation of the total dispersion curves is between about 10 and 16 per cent at the maximum from those averaged over azimuth. The acceptance of the retrieval of Green's function gives a proper physical meaning to the complex coherence function: the real part of the elastodynamic Green's function normalized by its zero‐offset version. This makes it possible to take a deterministic approach rather than the statistical one on which the conventional spatial autocorrelation method is based and gives fruitful new aspects and perspectives. For example, the formula for the multi‐mode case is given and the possibility of exploration of two or three dimensional velocity structures is suggested.     

基于波场分离的弹性波逆时偏移

尽管叠前逆时偏移成像精度高,但仅针对单一纵波的成像也可能形成地下介质成像盲区,由于基于弹性波方程的逆时偏移成像可形成多波模式的成像数据,因此弹性波逆时偏移成像可提供更为丰富的地下构造信息.本文依据各向同性介质的一阶速度-应力方程组构建震源和检波点矢量波场,再利用Helmholtz分解提取纯纵波和纯横波波场,使用震源归一化的互相关成像条件获得纯波成像,避免了直接使用坐标分量成像而引起的纵横波串扰问题.针对转换波成像的极性反转问题,文中提出一种共炮域极性校正方法.为有效节约存储成本,也提出一种适用于弹性波逆时偏移的震源波场逆时重建方法,在震源波场正传过程中,仅保存PML边界内若干层的速度分量波场,进而逆时重建出所有分量的震源波场.本文分别对地堑模型和Marmousi2模型进行了弹性波逆时偏移成像测试,结果表明:所提出的共炮域极性校正方法正确有效,基于波场分离的弹性波逆时偏移成像的纯波数据能够对复杂地下构造准确成像.     

Regional Seismic Wavefield Computation on a 3-D Heterogeneous Earth Model by Means of Coupled Traveling Wave Synthesis

-- I present a new algorithm for calculating seismic wave propagation through a three-dimensional heterogeneous medium using the framework of mode coupling theory originally developed to perform very low frequency (f < ~0.01т.05 Hz) seismic wavefield computation. It is a Greens function approach for multiple scattering within a defined volume and employs a truncated traveling wave basis set using the locked mode approximation. Interactions between incident and scattered wavefields are prescribed by mode coupling theory and account for the coupling among surface waves, body waves, and evanescent waves. The described algorithm is, in principle, applicable to global and regional wave propagation problems, but I focus on higher frequency (typically f S ~0.25 Hz) applications at regional and local distances where the locked mode approximation is best utilized and which involve wavefields strongly shaped by propagation through a highly heterogeneous crust. Synthetic examples are shown for P-SV-wave propagation through a semi-ellipsoidal basin and SH-wave propagation through a fault zone.     

A numerical model for the dynamic analysis of inclined pile groups

The paper presents a numerical model for the dynamic analysis of pile groups with inclined piles in horizontally layered soil deposits. Piles are modelled with Euler–Bernoulli beams, while the soil is supposed to be constituted by independent infinite viscoelastic horizontal layers. The pile–soil–pile interaction as well as the hysteretic and geometric damping is taken into account by means of two‐dimensional elastodynamic Green's functions. Piles cap is considered by introducing a rigid constraint; the condensation of the problem permits a consistent derivation of both the dynamic impedance matrix of the soil–foundation system and the foundation input motion. These quantities are those used to perform inertial soil–structure interaction analyses in the framework of the substructure approach. Furthermore, the model allows evaluating the kinematic stress resultants in piles resulting from waves propagating in the soil deposit, taking into account the pile–soil–pile interactions. The model validation is carried out by performing accuracy analyses and comparing results in terms of dynamic impedance functions, kinematic response parameters and pile stress resultants, with those furnished by 3D refined finite element models. To this purpose, classical elastodynamic solutions are adopted to define the soil–pile interaction problem. The model results in low computational demands without significant loss of precision, compared with more rigorous approaches or refined finite element models. Copyright © 2015 John Wiley & Sons, Ltd.     

A spatial windowing technique to account for finite dimensions in 2.5D dynamic soil–structure interaction problems

The dynamic interaction between a layered halfspace and quasi translationally invariant structures such as roads, railway tracks, tunnels, dams, and lifelines can be modelled using a computationally efficient 2.5D approach, assuming invariance of the geometry in the longitudinal direction. This assumption is not always fulfilled in practice, however. Even for elongated structures, full 3D computations may be required for an accurate solution of the dynamic soil–structure interaction problem. This paper presents a spatial windowing technique for elastodynamic transmission and radiation problems that allows accounting for the finite length of a structure, still maintaining the computational efficiency of a 2.5D formulation. The proposed technique accounts for the diffraction occurring at the structure's edges, but not for its modal behaviour resulting from reflections of waves at its boundaries. Numerical examples of a barrier for vibration transmission and a surface foundation are discussed to demonstrate the accuracy and applicability of the proposed methodology. Full 3D calculations are performed to provide a rigorous validation for each of these examples. It is demonstrated that the proposed technique is appropriate as long as the response is not dominated by the resonant behaviour of individual modes of the structure.     

Elastodynamic boundary element formulation employing discontinuous curved elements

A boundary element formulation having discontinuous curved quadratic elements is presented for 2D elastodynamics. The first fundamental solution for static case is subtracted from and added to the first fundamental solution for dynamic case. As both kernels have the same order of singularity, the integral involving the regular expression arising from the subtraction can be calculated. matrix is calculated by employing the well-known rigid-body motion technique. The formulation is performed in Fourier transform space. Based on the formulation presented in this study, a general purpose computer program is developed for elastic or visco-elastic 2D elastodynamic problems. The program performs the analysis in Fourier transform space and can also be used for static analysis by assigning a very small value close to zero for the frequency. The results of some elastodynamic and dynamic soil–structure interaction problems obtained using the present study are compared with those in the literature.     

Influences of grass and moss on runoff and sediment yield on sloped loess surfaces under simulated rainfall

It is important to evaluate the impacts of grasses on soil erosion process so as to use them effectively to control soil and water losses on the Loess Plateau. Laboratory-simulated rainfall experiments were conducted to investigate the runoff and sediment processes on sloped loess surfaces with and without the aboveground parts of grasses and moss (GAM: grass and moss; NGAM: no grass and moss) under slope gradients of 5°, 10°, 15°, 20°, 25° and 30°. The results show that runoff from GAM and NGAM plots increased up to a slope gradient of 10° and decreased thereafter, whereas the runoff coefficients increased with gradient. The average runoff rates and runoff coefficients of NGAM plots were less than those of GAM plots except for the 5° slope. This behaviour may be due to the reduction in water infiltration under moss. The difference between GAM and NGAM plots in average runoff rates varied from 1·4 to 8%. At the same gradients, NGAM plots yielded significantly (α = 0·05) more sediment than GAM plots. Average sediment deliveries for different slopes varied from 0·119 to 3·794 g m⁻² min⁻¹ from GAM plots, and from 0·765 to 16·128 g m⁻² min⁻¹ from NGAM plots. Sediment yields from GAM plots were reduced by 45 to 85%, compared with those from the NGAM plots. Plots at 30° yielded significantly higher sediments than at the other gradients. Total sediments S increased with slope gradients G in a linear form, i.e. S = 9·25G − 39·6 with R² = 0·77*, for the GAM plots, and in an exponential model, i.e. S = 40·4 exp(0·1042G) with R² = 0·93**, for the NGAM plots. In all cases, sediment deliveries decreased with time, and reached a relative steady state at a rainfall duration of 14 min. Compared with NGAM plots, the final percentage reductions in sediment delivery from GAM plots were higher than those at the initial time of rainfall at all slopes. Copyright © 2006 John Wiley & Sons, Ltd.     

Sediment concentration in interrill flow: interactions between soil surface conditions,vegetation and rainfall

A database composed of 673 natural rainfall events with sediment concentration measurements at the field or plot scale was analysed. Measurements were conducted on similar soil type (loess soils prone to sealing phenomenon) to apprehend the variability and complexity involved in interrill erosion processes attributable to soil surface conditions. The effects of the dominant controlling factors are not described by means of equations; rather, we established a classification of potential sediment concentration domain according to combination of the dominant parameters. Thereby, significant differences and evolution trends of mean sediment concentration between the different parameter categories are identified. Further, when parameter influences interact, it allows us to discern the relative effects of factors according to their respective degree of expression. It was shown that crop cover had a major influence on mean sediment concentration, particularly when soil surface roughness is low and when maximum 6‐min intensity of rainfall events exceeds 10 mm h^?1: mean sediment concentration decreases from 8·93 g l^?1 for 0–20 per cent of coverage to 0·97 g l^?1 for 21–60 per cent of coverage. The established classification also indicates that the increase of the maximum 6‐min intensity of the rainfall factor leads to a linear increase of mean sediment concentration for crop cover over 21 per cent (e.g. from 2·96 g l^?1 to 14·44 g l^?1 for the 1–5 cm roughness class) and to an exponential increase for low crop cover (e.g. from 3·92 g l^?1 to 58·76 g l^?1 for the 1–5 cm roughness class). The implication of this work may bring perspective for erosion prediction modelling and give references for the development of interrill erosion equation. Copyright © 2002 John Wiley & Sons, Ltd.     

Spatial variation of isotope composition in precipitation in a tropical environment: a case study from the Deduru Oya river basin,Sri Lanka

The stable isotope composition (¹⁸O and ²H) in the tropical precipitation collected from 18 locations throughout the Deduru Oya river basin in Sri Lanka, has been studied during August and September 2001, in order to characterize the isotopic composition of precipitation in the dry and intermediate climatic zones of Sri Lanka. The isotope compositions are described with respect to the distance from the coast and the altitude. The analyses show that δ¹⁸O vary from ? 5·11 to 1·39‰ and δD vary from ? 35·71 to 12·55‰. The d‐excess values range from ? 0·65 to 13·17 with an average value of ～7. Regression for the δ¹⁸O ? δD is y = 6·8x + 4·9 (R² = 0·9) which is compatible with the precipitation in other tropical regions. The lower slope in the regression line and the lower d‐excess value indicate high temperature events which were possibly aided by concentration through successive evaporation within the atmosphere. The spatial variation of isotope composition indicates two different cloud contributions for the rain events, of which one may be linked to the Indian Ocean contribution and the other to the high altitude condensation. Copyright © 2008 John Wiley & Sons, Ltd.     

Zur ökonomischen Bedeutung der biologischen Selbstreinigung

In twelve flowing water (Q 1.04… 30.4 m³/s, 0.5… 1.2 m/s) the self-purification efficiency in 110 river sections (L 1.3… 26.4 km) is determined as the load difference in kg · d^?1 COD-Mn as well as in g · m^?3 referred to the daily passage. In the economical comparison, the self-purification efficiency is valued as substitution for wastewater treatment plants of the same capacity. At a mean specific capacity of 8.6 g · m^?3, the self-purification efficiency of the investigated waters is equivalent to economical values of 124 · 10³… 534 · 10³ M km^?1 investment costs, 11 · 10³… 80 · 10³ M · a^?1km^?1 operating costs and 5… 81 MWh · a^?1 km^?1 expenditure of energy. The specific capacity in g · m^?3 COD-Mn shows an exponential regression to the degree of saprobity (L = 0.015 · exp (1.7358 · S)). From this empirical model the limits of the self-purification capacity of aerobic waters by oxygen input can be detected: in respective examples more than 50% of the required oxygen input are due to weirs.     

Applying local Green's functions to study the influence of the crustal structure on hydrological loading displacements

The influence of the elastic Earth properties on seasonal or shorter periodic surface deformations due to atmospheric surface pressure and terrestrial water storage variations is usually modeled by applying a local half-space model or an one dimensional spherical Earth model like PREM from which a unique set of elastic load Love numbers, or alternatively, elastic Green's functions are derived. The first model is valid only if load and observer almost coincide, the second model considers only the response of an average Earth structure. However, for surface loads with horizontal scales less than 2500 km², as for instance, for strong localized hydrological signals associated with heavy precipitation events and river floods, the Earth elastic response becomes very sensitive to inhomogeneities in the Earth crustal structure.We derive a set of local Green's functions defined globally on a 1° × 1° grid for the 3-layer crustal structure TEA12. Local Green's functions show standard deviations of ±12% in the vertical and ±21% in the horizontal directions for distances in the range from 0.1° to 0.5°. By means of Green's function scatter plots, we analyze the dependence of the load response to various crustal rocks and layer thicknesses. The application of local Green's functions instead of a mean global Green's function introduces a variability of 0.5–1.0 mm into the hydrological loading displacements, both in vertical and in horizontal directions. Maximum changes due to the local crustal structures are from −25% to +26% in the vertical and −91% to +55% in the horizontal displacements. In addition, the horizontal displacement can change its direction significantly. The lateral deviations in surface deformation due to local crustal elastic properties are found to be much larger than the differences between various commonly used one-dimensional Earth models.     

Long-range Automaton Models of Earthquakes: Power-law Accelerations,Correlation Evolution,and Mode-switching

—?We introduce a conceptual model for the in-plane physics of an earthquake fault. The model employs cellular automaton techniques to simulate tectonic loading, earthquake rupture, and strain redistribution. The impact of a hypothetical crustal elastodynamic Green's function is approximated by a long-range strain redistribution law with a r ^?p dependance. We investigate the influence of the effective elastodynamic interaction range upon the dynamical behaviour of the model by conducting experiments with different values of the exponent (p). The results indicate that this model has two distinct, stable modes of behaviour. The first mode produces a characteristic earthquake distribution with moderate to large events preceeded by an interval of time in which the rate of energy release accelerates. A correlation function analysis reveals that accelerating sequences are associated with a systematic, global evolution of strain energy correlations within the system. The second stable mode produces Gutenberg-Richter statistics, with near-linear energy release and no significant global correlation evolution. A model with effectively short-range interactions preferentially displays Gutenberg-Richter behaviour. However, models with long-range interactions appear to switch between the characteristic and GR modes. As the range of elastodynamic interactions is increased, characteristic behaviour begins to dominate GR behaviour. These models demonstrate that evolution of strain energy correlations may occur within systems with a fixed elastodynamic interaction range. Supposing that similar mode-switching dynamical behaviour occurs within earthquake faults then intermediate-term forecasting of large earthquakes may be feasible for some earthquakes but not for others, in alignment with certain empirical seismological observations. Further numerical investigation of dynamical models of this type may lead to advances in earthquake forecasting research and theoretical seismology.     

Dissolved organic carbon export from a cutover and restored peatland

High demand for horticultural peat has increased peatland drainage and peat extraction in Canada. The hydrology and carbon cycling of these cutover peatlands is greatly altered, necessitating active restoration efforts to permit the regeneration of Sphagnum mosses and the re‐establishment of natural peatland function. The effect of peatland extraction and restoration on the export of dissolved organic carbon (DOC) was examined for three successive seasons (May to October, 1999 to 2001) at two different sites (cutover and restored) in eastern Québec. A shift towards higher DOC concentrations was observed following peatland extraction (maximum: 182·6 mg L^?1) and concentrations remained high post‐restoration (maximum: 191·0 mg L^?1). The cutover site exported more DOC than the restored site in all three study seasons. The highest exports occurred during the wettest year (1999), with cutover and restored site export of 10·3 and 4·8 g m^?2, respectively. In 2000, 8·5 g C m^?2 was released from the cutover site, while the restored site released less than half that amount (3·4 g C m^?2). In 2001, the restored site released about the same amount of DOC as in the previous year (3·5 g C m^?2), while the cutover site load dropped to 6·2 g C m^?2. Both sites were net exporters of DOC in all years. Copyright © 2007 John Wiley & Sons, Ltd.     

Atmospheric chloride deposition in continental Spain

The atmospheric bulk deposition rate of chloride in continental Spain was studied to get basic information in order to help in the evaluation of diffuse recharge to aquifers through an environmental chemical balance. Both new, recent data and bibliographic data have been used. Most sampling records are less than 5 years long and often only 1 year long. This means that the calculated mean yearly bulk deposition rate of chloride is quite uncertain by 30% on average, and larger than the values derived form records up to 15 years long. A map of atmospheric bulk deposition of chloride has been drawn using ordinary kriging. The mean bulk deposition rate of chloride varies from 1 to 30 g m^?2 year^?1 in coastal areas, with strong negative landward gradients between 0·1 and 1 g m^?2 year^?1 km^?1. In the centre of the Iberian Peninsula, chloride deposition rates vary from 0·2 to 0·5 g m^?2 year^?1, with gradients around or less than 5 × 10^?3 g m^?2 year^?1 km^?1. The coefficient of variation of the mean bulk atmospheric deposition rate of chloride, for any place, ranges from 0·1 to 1. Values larger than ～0·5 are not a good indicator of natural uncertainty for this series of data that has a skewed distribution. The map of bulk deposition rate and its error is one of the terms needed for aquifer recharge estimation by means of the chloride ion balance. Copyright © 2008 John Wiley & Sons, Ltd.     

Mapping of erosion rates in marly badlands based on a coupling of anatomical changes in exposed roots with slope maps derived from LiDAR data

Black marls form very extensive outcrops in the Alps and constitute some of the most eroded terrains, thus causing major problems of sedimentation in artificial storage systems (e.g. reservoirs) and river systems. In the experimental catchments near Draix (France), soil erosion rates have been measured in the past at the plot scale through a detailed monitoring of surface elevation changes and at the catchment scale through continuous monitoring of sediment yield in traps at basin outlets. More recently, erosion rates have been determined by means of dendrogeomorphic techniques in three monitored catchments of the Draix basin. A total of 48 exposed roots of Scots pine have been sampled and anatomical variations in annual growth rings resulting from denudation analysed. At the plot scale, average medium‐term soil erosion rates derived from exposed roots vary between 1·8 and 13·8 mm yr^?1 (average: 5·9 mm yr^?1) and values are significantly correlated with slope angle. The dendrogeomorphic record of point‐scale soil erosion rates matches very well with soil erosion rates measured in the Draix basins. Based on the point‐scale measurements and dendrogeomorphic results obtained at the point scale, a linear regression model involving slope angle was derived and coupled to high‐resolution slope maps obtained from a LiDAR‐generated digital elevation model so as to generate high‐resolution soil erosion maps. The resulting regression model is statistically significant and average soil erosion rates obtained from the areal erosion map (5·8, 5·2 and 6·2 mm yr^?1 for the Roubine, Moulin and Laval catchments, respectively) prove to be well in concert with average annual erosion rates measured in traps at the outlet of these catchments since 1985 (6·3, 4·1 and 6·4 mm yr^?1). This contribution demonstrates that dendrogeomorphic analyses of roots clearly have significant potential and that they are a powerful tool for the quantification and mapping of soil erosion rates in areas where measurements of past erosion is lacking. Copyright © 2011 John Wiley & Sons, Ltd.     

Factors controlling stable oxygen,hydrogen and carbon isotope ratios in regional groundwater of the eastern United Arab Emirates (UAE)

Increasing groundwater salinity and depletion of the aquifers are major concerns in the UAE. Isotopes of oxygen, hydrogen, and carbon concentrations in groundwater were used to estimate evaporation loss using the isotopes of oxygen and hydrogen, and using a carbon isotope to trace inorganic carbon cycling in two main aquifers in the eastern part of the United Arab Emirates. The δD‐δ¹⁸O of groundwater samples plotted on a line given by: δD = 4 δ¹⁸O + 4 ·4 (r² = 0·4). In comparison, the local meteoric water line (LMWL) has been defined by the line: δD = 8 δ¹⁸O + 15. In order to better understand the system investigated, samples were separated into two groups based on the δD‐δ¹⁸O relationship. These are (1) samples that plot above the LMWL (δD = 6·1 δ¹⁸O + 12·4, r² = 0·8) and which are located predominantly in the north of the study area, and (2) samples that plot below the LMWL (δD = 5·6 δ¹⁸O + 6·2, r² = 0·8) and which are mostly distributed in the south. Slopes for both the groups are similar and lower than that for LMWL indicating potential evaporation of recharging water. However, the y‐intercept, which differs between the two groups, suggests evaporation of return flow and evapotranspiration in the unsaturated zone to be more significant in the south. This is attributed to intense agricultural activities in the region. Samples from the eastern Gravel Plain aquifer have δ¹³C and dissolved inorganic carbon (DIC) values in the range from ? 10 to 17‰, and 12–100 mg C/l, respectively, while the range for those from the Ophiolite aquifer is from ? 11 to ? 16.4‰, and 16–114 mg C/l respectively. This suggests the control of C‐3 and C‐4 plants on DIC formation, an observation supported by the range δ¹³C of soil organic matter (from ? 18·5 to ? 22·1‰.) Copyright © 2007 John Wiley & Sons, Ltd.     

Isotopic composition (δ18O and δD) of precipitation and groundwater in a semi‐arid,mountainous area (Guadiana Menor basin,Southeast Spain)

We characterize the precipitation and groundwater in a mountainous (peaks slightly above 3000 m a.s.l.), semi‐arid river basin in SE Spain in terms of the isotopes ¹⁸O and ²H. This basin, with an extension of about 7000 km², is an ideal site for such a study because fronts from the Atlantic and the Mediterranean converge here. Much of the land is farmed and irrigated both by groundwater and runoff water collected in reservoirs. A total of approximately 100 water samples from precipitation and 300 from groundwater have been analysed. To sample precipitation we set up a network of 39 stations at different altitudes (800–1700 m a.s.l.), with which we were able to collect the rain and snowfall from 29 separate events between July 2005 and April 2007 and take monthly samples during the periods of maximum recharge of the aquifers. To characterize the groundwater we set up a control network of 43 points (23 springs and 20 wells) to sample every 3 months the main aquifers and both the thermal and non‐thermal groundwater. We also sampled two shallow‐water sites (a reservoir and a river). The isotope composition of the precipitation forms a local meteoric water line (LMWL) characterized by the equation δD = 7·72δ¹⁸O + 9·90, with mean values for δ¹⁸O and δD of − 10·28‰ and − 69·33‰, respectively, and 12·9‰ for the d‐excess value. To correlate the isotope composition of the rainfall water with groundwater we calculated the weighted local meteoric water line (WLMWL), characterized by the equation δD = 7·40δ¹⁸O + 7·24, which takes into account the quantity of water precipitated during each event. These values of (dδD/dδ¹⁸O)< 8 and d‐excess (δD–8δ¹⁸O)< 10 in each curve bear witness to the ‘amount effect’, an effect which is more manifest between May and September, when the ground temperature is higher. Other effects noted in the basin were those of altitude and the continental influence. The isotopic compositions of the groundwater are represented by the equation δD = 4·79δ¹⁸O − 18·64. The groundwater is richer in heavy isotopes than the rainfall, with mean values of − 8·48‰ for δ¹⁸O and − 59·27‰ for δD. The isotope enrichment processes detected include a higher rate of evaporation from detrital aquifers than from carbonate ones, the effects of recharging aquifers from irrigation return flow and/or from reservoirs' leakage and enrichment in δ¹⁸O from thermal water. Copyright © 2010 John Wiley & Sons, Ltd.     

The role of cryptogams in runoff and erosion control on bariland in the Nepal Middle Hills of the Southern Himalaya

Cryptogams are communities of non-vascular plants that live on the soil surface. Numerous functions have been attributed to these crusts, including changes in soil fertility and nutrient status, soil hydrology and soil erosion. Most significant for this paper is the reported benefit of cryptogams in reducing soil erosion by water in semi-arid areas. However, to date there have been few attempts to understand the soil conservation value of cryptogams in subsistence agricultural systems or in humid mountain environments. This paper investigates the potential of cryptogams in soil erosion by water on agricultural hillslope terraces (bariland) in the Nepal Middle Hills of the southern monsoonal Himalaya. The research is significant because the loss of fertile topsoil is considered by some to be the biggest threat to the livelihoods of subsistence farmers in the area in the medium and long term. The current study was conducted in the field between two of the weeding events that take place under maize cover, grown in the traditional manner. Three groundcover types which represented (i) maize only (types A), (ii) maize and weed cover (types B), and (iii) maize and cryptogam cover (types C) were monitored utilizing multiple microerosion plots. Measurements of runoff and soil loss data were collected sequentially on a storm-by-storm basis throughout the monitored period from 24 July 1997 to 29 August 1997. Measurements of infiltration rates were also taken on each of the groundcover types at selected times. Results collected from the erosion plots demonstrate that runoff and soil losses over distances of <2 m can be significantly reduced by up to 50 per cent with cryptogam cover, compared to maize-only canopies. Mean runoff for all storm events sampled from plot types A, B and C were 3·4 l m⁻², 1·6 l m⁻² and 1·5 l m⁻² respectively. For soil loss, the results were 21·7 g m⁻², 11·3 g m⁻² and 10·2 g m⁻² respectively. Therefore, cryptogams would appear to offer a similar degree of protection to the soil surface from runoff and raindrop erosion, to that afforded by weed cover. Weed and cryptogam covers protect the soil surface from rainfall kinetic energies and work to preserve surface microtopographies, depressional storage and surface water detention. Terminal infiltration rates taken at the end of the monitored period showed that well developed maize- and cryptogam-covered soil surfaces (types C) have a mean terminal infiltration rate of 35·0 mm h⁻¹ compared to 44·5 mm h⁻¹ for comparable maize- and weed-covered soil surfaces (types B), and 15·5 mm h⁻¹ for maize-only soil surfaces (types A). These results show that cryptogams and weeds also have relatively higher infiltration rates than comparable maize-only covered plots, devoid of groundcover. The findings in this study may have implications for traditional weed management practices used by local hill farmers, which often destroy cryptogam soil coatings two to three times during the maize growing period. However, further work needs to be done to ascertain farmers' understandings of cryptogams. It is hoped that conservationists will benefit from incorporating cryptogams into the design of future soil erosion studies relating to development programmes. Copyright © 2001 John Wiley & Sons, Ltd.     

The direct boundary element method: 2D site effects assessment on laterally varying layered media (methodology)

The Direct Boundary Element Method (DBEM) is presented to solve the elastodynamic field equations in 2D, and a complete comprehensive implementation is given. The DBEM is a useful approach to obtain reliable numerical estimates of site effects on seismic ground motion due to irregular geological configurations, both of layering and topography.The method is based on the discretization of the classical Somigliana's elastodynamic representation equation which stems from the reciprocity theorem. This equation is given in terms of the Green's function which is the full-space harmonic steady-state fundamental solution. The formulation permits the treatment of viscoelastic media, therefore site models with intrinsic attenuation can be examined. By means of this approach, the calculation of 2D scattering of seismic waves, due to the incidence of P and SV waves on irregular topographical profiles is performed. Sites such as, canyons, mountains and valleys in irregular multilayered media are computed to test the technique. The obtained transfer functions show excellent agreement with already published results.