A hybrid indirect boundary element—discrete wave number method applied to simulate the seismic response of stratified alluvial valleys

A hybrid indirect boundary element – discrete wavenumber method is presented and applied to model the ground motion on stratified alluvial valleys under incident plane SH waves from an elastic half-space. The method is based on the single-layer integral representation for diffracted waves. Refracted waves in the horizontally stratified region can be expressed as a linear superposition of solutions for a set of discrete wavenumbers. These solutions are obtained in terms of the Thomson–Haskell propagators formalism. Boundary conditions of continuity of displacements and tractions along the common boundary between the half-space and the stratified region lead to a system of equations for the sources strengths and the coefficients of the plane wave expansion. Although the regions share the boundary, the discretization schemes are different for both sides: for the exterior region, it is based on the numerical and analytical integration of exact Green's functions for displacements and tractions whereas for the layered part, a collocation approach is used. In order to validate this approach results are compared for well-known cases studied in the literature. A homogeneous trapezoidal valley and a parabolic stratified valley were studied and excellent agreement with previous computations was found. An example is given for a stratified inclusion model of an alluvial deposit with an irregular interface with the half-space. Results are displayed in both frequency and time domains. These results show the significant influence of lateral heterogeneity and the emergence of locally generated surface waves in the seismic response of alluvial valleys.     

Food-web manipulation of drinking water reservoirs with salmonids: vertical distribution of prey and predator

In a novel biomanipulation experiment salmonids were used as a tool to improve water quality. The manipulation was initiated in spring 2000 as a response to non-point sources of phosphorus in a drinking water reservoir in Saxony, Germany. Salmonids (brown trout, Salmo trutta forma lacustris) were chosen as predators as the reservoir has a large hypolimnic water body and surface temperatures rarely exceed 20 °C. The vertical distributions of prey fish and brown trout were analysed with a fleet of vertical gill nets set in the pelagic zone of the reservoir. Consumption of brown trout was estimated by means of a bioenergetic model and the diet analyses of the trout. While the dominant planktivore (roach, Rutilus rutilus) was caught almost exclusively in the epilimnion during the stratification period trout were caught mainly below a depth of 10 m. Diet analysis revealed that the trout performed vertical migrations to consume food in the epilimnic layer, as an important food component were adult terrestrial and aquatic insects. The amount of fish in the food increased strongly with the size of the brown trout. The consumption estimate suggested that the trout had consumed 2-3% of the total roach stock during the study period (May-November 2000) of the first year of biomanipulation. We conclude that in general salmonids are suitable for food-web manipulation in deep reservoirs, but the stocked fish should be as large as possible (> 300 mm) and the proportion of large trout (> 500 mm) should be as high as possible.     

System identification of the Vincent Thomas suspension bridge using earthquake records

The Vincent Thomas Bridge in the Los Angeles metropolitan area, is a critical artery for commercial traffic flow in and out of the Los Angeles Harbor, and is at risk in the seismically active Southern California region, particularly because it straddles the Palos Verdes fault zone. A combination of linear and non‐linear system identification techniques is employed to obtain a complete reduced‐order, multi‐input–multi‐output (MIMO) dynamic model of the Vincent Thomas Bridge based on the dynamic response of the structure to the 1987 Whittier and 1994 Northridge earthquakes. Starting with the available acceleration measurements (which consists of 15 accelerometers on the bridge structure and 10 accelerometers at various locations on its base), an efficient least‐squares‐based time‐domain identification procedure is applied to the data set to develop a reduced‐order, equivalent linear, multi‐degree‐of‐freedom model. Although not the main focus of this study, the linear system identification method is also combined with a non‐parametric identification technique, to generate a reduced‐order non‐linear mathematical model suitable for use in subsequent studies to predict, with good fidelity, the total response of the bridge under arbitrary dynamic environments. Results of this study yield measurements of the equivalent linear modal properties (frequencies, mode shapes and non‐proportional damping) as well as quantitative measures of the extent and nature of non‐linear interaction forces arising from strong ground shaking. It is shown that, for the particular subset of observations used in the identification procedure, the apparent non‐linearities in the system restoring forces are quite significant, and they contribute substantially to the improved fidelity of the model. Also shown is the potential of the identification technique under discussion to detect slight changes in the structure's influence coefficients, which may be indicators of damage and degradation in the structure being monitored. Difficulties associated with accurately estimating damping for lightly damped long‐span structures from their earthquake response are discussed. The technical issues raised in this paper indicate the need for added spatial resolution in sensor instrumentation to obtain identified mathematical models of structural systems with the broadest range of validity. Copyright © 2003 John Wiley & Sons, Ltd.     

Observing artificially induced strain: implications for subglacial deformation

By analysing a series of four successive thin‐sections from a ceramic clay that was subjected to uniaxial compression, we were able to monitor the development of microstructures in a fine‐grained sediment. The artificially induced microstructures, such as unidirectional clay reorientations and linear and circular grain arrangements, are identical to features that have been observed in thin‐sections of subglacially deformed tills, and therefore may be used as representative analogues. We argue that the structures, reflecting slip, planar shear displacements as well as rotational movements, can be explained by assuming a Coulomb‐plastic response to imposed shear. We conclude that sediments subjected to subglacial deformation behave as Coulomb materials, at least during the final stages of the deformation. The present study bridges the gap between field studies, experimental studies and theoretical modelling. The microscopic observations assist in visualising inferred subglacial processes and facilitate up‐ and downscaling between diverse methodological approaches. Copyright © 2003 John Wiley & Sons, Ltd.     

The use of inclined hemisphere projections for analyzing failure mechanisms in discontinuous rocks

This paper demonstrates the advantages of using inclined stereographic projections in kinematic analysis of rock blocks in discontinuous rock masses. Some examples of limiting cases are presented. The application of inclined projections is illustrated by its use in a mine slope in Brazil. It is clear from the discussion of these examples that inclined hemisphere projections provide better results than horizontal projections. It is also demonstrated that horizontal projections can lead to incorrect results in limiting cases.     

Engineering geology of loess and loess-like deposits: a commentary on the Russian literature

Loess and loess-like deposits were much studied in the Soviet Union, and are currently under investigation in Russia and surrounding countries. There is a vast literature in Russian, which touches on all aspects of loess science and technology. In particular, the studies of the origin of collapsibility are almost totally in Russian, and of course studies on the various regions of Russia and the countries of the Former Soviet Union FSU appear in Russian. This review looks at the literature in Russian and attempts to pick out key contributors, major topics and works and to identify the critical regions and zones of investigation. Because so many regions of the FSU had people living on loess ground, there is a vast literature on engineering geology and ground engineering topics, and this tends to dominate all the literature on loess in Russian. Following Russian practice, the fine-grained deposits under consideration are divided into loess and loess-like deposits. Three main topics are recognised across the whole spectrum of loess research: formation and distribution of loess deposits; stratigraphy, cyclicity and palaeoclimatology; and engineering topics, in particular hydrocollapse and subsidence, and we concentrate on the engineering geology topics. An attempt is made, based on the map of Abelev and Abelev [Abelev, Yu.M., Abelev, M.Yu., 1968. Fundamentals of design and construction on collapsible macroporous soils, 2nd ed. Stroiizdat, Moscow, 431 pp. (in Russian)] of collapsing loess deposits, to define seven loess regions within the geographical limits of the old USSR. The seven regions are those where geotechnical problems might be expected.     

Digital model of corrected accumulated flow for peak discharge data acquisition and drainage system design

The design of a drainage system for a roofing slate quarry was implemented by the enhancement of discharge peak estimation, and the uncertainty inevitably associated with the engineering model was reduced.

The development of a topographical, geological, and vegetation cover database developed from a Geographical Information System (GIS) allowed for the definition of the drainage network for a hydraulic system, along with the calculation of the runoff coefficient. This is applied to the digital model of accumulated flow (DMF) as a weight correction coefficient, using a matrix-based model at 5×5 m resolution. The new digital model of corrected accumulated flow (DMCF) is the result of combining the thematic maps with the map of slope <3%, which was previously created from the slope model. It is demonstrated that this new model allows to apply the “Rational Method” on cartographic units defined by the GIS.

The DMCF is compared with other traditional applications of the Rational Method based on the calculation of the discharge peak considering: (1) the drainage basin as a single watershed or (2) defining an average runoff coefficient in each sub-watershed. Both approaches have bigger discharge peaks than those obtained by the DMCF since the slope, lithology, and vegetation cover have average values, and the runoff coefficient is poorly defined, increasing the uncertainty in the discharge peak.     

Groundwater residence times and flow paths in fractured rock determined using environmental tracers in the Mission Tunnel; Santa Barbara County,California, USA

Multiple geochemical tracers [ion chemistry, stable isotopes of water, chlorofluorocarbons (CFC), tritium] and a 25-year-long record of discharge were used to understand residence times and flow paths of groundwater seeps in the fractured rock aquifer surrounding the Mission Tunnel, Santa Barbara, California. Tritium data from individual seeps indicate that seep waters are a mixture of >45-year-old (recharged prior to the nuclear bomb tests) and young groundwater. CFC data support this interpretation, however, a two-end member mixing model cannot completely explain the age tracer data. Microbial degradation and partial re-equilibration complicate the CFC signal. Spectral analysis of precipitation and groundwater seepage records shows that seepage lags precipitation by 3 months. This delay is related to the advancement of the wetting front and increasing the number of active flow paths. Additionally, the amount of seepage produced by precipitation is less during extended periods of drought than during normal or wet periods, suggesting antecedent conditions strongly affect flow through this fractured rock aquifer.     

Influence of accretion-disk models on the structure of the Iron Κ<Subscript>α</Subscript> line

Broad iron Κ_α emission with a characteristic two-peaked profile is observed in most Seyfert galaxies in the X-ray. We have calculated the profiles of such lines emitted by an accretion disk in a Schwartzschild metric. The dependence of the temperature distribution in the disk on the line shape is demonstrated. All the calculations include general relativistic effects. The disk material is assumed to move in circular geodesics in the equatorial plane. The line profile is extremely complex, even in a traditional model for the radial temperature distribution, complicating interpretation of the observational data.