Zero velocity surfaces for the general planar three-body problem

The angular momentum and the energy integral of the planar three-body problem are used to establish regions of the physical space where motion is allowed to take place. Although forbidden regions exist for both negative and positive values of the energy of the system, the known integrals of the motion always allow for at least one of the three bodies to escape.     

Multi‐variable relations for soil effects on seismic ground motion

Soil effects on peak ground acceleration, velocity and elastic response spectra (5% damping) are expressed by simple approximate relations in terms of five key parameters: (a) the fundamental vibration period of the non‐linear soil, T_S, (b) the period of a bedrock site of equal thickness, T_b, (c) the predominant excitation period, Te, (d) the peak seismic acceleration at outcropping bedrock, a, and (e) the number of significant excitation cycles, n. Furthermore, another relation is proposed for the estimation of T_S in terms of the soil thickness H, the average shear wave velocity of the soil V?_S,o and a. The aforementioned parameters were first identified through a simplified analytical simulation of the site excitation. The multivariable approximate relations were then formulated via a statistical analysis of relevant data from more than 700 one‐dimensional equivalent‐linear seismic ground response analyses, for actual seismic excitations and natural soil conditions. Use of these relations to back‐calculate the numerical results in the database gives an estimate of their error margin, which is found to be relatively small and unbiased. The proposed relations are also independently verified through a detailed comparison with strong motion recordings from seven well‐documented case studies: (a) two sites in the San Fernando valley during the Northridge earthquake, and (b) five different seismic events recorded at the SMART‐1 accelerometer array in Taiwan. It is deduced that the accuracy of the relations is comparable to that of the equivalent‐linear method. Hence, they can be readily used as a quick alternative for routine applications, as well as for spreadsheet computations (e.g. GIS‐aided seismic microzonation studies) where numerical methods are cumbersome to implement. Copyright © 2003 John Wiley & Sons, Ltd.     

Organochlorine pesticide residues in sediments and biota from the coastal area of Dar es Salaam city,Tanzania

Sediment and biota samples were collected from Msimbazi and Kizinga rivers and from the coastal marine environment of Dar es Salaam during both dry and wet seasons. The samples were analyzed for various organochlorine pesticide residues using GC-ECD and GC-MS. Dieldrin, p,p^′-DDT, p,p^′-DDE, p,p^′-DDD, o,p^′-DDT and γ-HCH were detected at significantly greater concentrations above the method detection limits. Recoveries of pesticide residues ranged 86.5–120% in sediments and 62–102% in biota. The average concentrations of total DDT in sediments for the two seasons were almost the same. Biota samples showed significant difference in levels of residues depending on mode of feeding and age of analyzed biota. p,p^′-DDT to total DDT ratios in all matrices indicated recent use of DDT. The levels of residues in sediments suggest possible adverse effects to humans consuming biota that are directly exposed to the sediments. This effect, which is associated with bio-concentration of residues in the tissues of edible aquatic biota, should not affect wildlife. In spite of concern on the adverse effect to humans posed by pesticide residues in sediments, edible biota from the examined water bodies were found to be safe for human consumption.     

Hydraulic Effects of Shales in Fluvial-Deltaic Deposits: Ground-Penetrating Radar, Outcrop Observations, Geostatistics, and Three-Dimensional Flow Modeling for the Ferron Sandstone, Utah

Ground-penetrating radar (GPR) surveys, outcrop measurements, and cores provide a high-resolution 3D geologic model to investigate the hydraulic effects of shales in marine-influenced lower delta-plain distributary channel deposits within the Cretaceous-age Ferron Sandstone at Corbula Gulch in central Utah, USA. Shale statistics are computed from outcrop observations. Although slight anisotropy was observed in mean length and variogram ranges parallel and perpendicular to pale of low , the anisotropy is not statistically significant and the estimated mean length is 5.4 m. Truncated Gaussian simulation was used to create maps of shales that are placed on variably dipping stratigraphic surfaces interpreted from high-resolution 3D GPR surveys, outcrop interpretations, and boreholes. Sandstone permeability is estimated from radar responses calibrated to permeability measurements from core samples. Experimentally designed flow simulations examine the effects of variogram range, shale coverage fraction, and trends in shale coverage on predicted upscaled permeability, breakthrough time, and sweep efficiency. Approximately 1500 flow simulations examine three different geologic models, flow in the 3 coordinate directions, 16 geostatistical parameter combinations, and 10 realizations for each model. ANOVA and response models computed from the flow simulations demonstrate that shales decrease sweep, recovery, and permeability, especially in the vertical direction. The effect on horizontal flow is smaller. Flow predictions for ideal tracer displacements at Corbula Gulch are sensitive to shale-coverage fraction, but are relatively insensitive to twofold variations in variogram range or to vertical trends in shale coverage. Although the hydraulic effects of shale are statistically significant, the changes in flow responses rarely exceed 20%. As a result, it may be reasonable to use simple models when incorporating analogous shales into models of reservoirs or aquifers.     

On the Dynamics of the Jovian Ionosphere and Thermosphere: III. The Modelling of Auroral Conductivity

Recent work has been concerned with calculating the three-dimensional ion concentrations and Pedersen and Hall conductivities within the auroral region of Jupiter for varying conditions of incident electron precipitation. Using the jovian ionospheric model, we present results that show the auroral ionospheric response to changing the incoming flux of precipitating electrons (for constant initial energy) and also the response to changing the initial energy (for both constant flux and constant energy flux). The results show that, for expected energy fluxes of precipitating particles, the average auroral integrated Pedersen conductivity attains values in excess of 1 mho. In addition, it is shown that electrons with an initial energy of around 60 keV are particularly effective at generating auroral conductivity: Particles of this energy penetrate most effectively to the layer of the jovian ionosphere at which the auroral conductivity is at a maximum.     

Contents to Volume 8

Volume Contents

Contents to Volume 8     

Yielding oscillator subjected to simple pulse waveforms: numerical analysis & closed‐form solutions

Numerical and analytical solutions are presented for the elastic and inelastic response of single‐degree‐of‐freedom yielding oscillators to idealized ground acceleration pulses. These motions are typical of near‐fault earthquake recordings generated by forward rupture directivity and may inflict damage in the absence of substantial structural strength and ductility capacity. Four basic pulse waveforms are examined: (1) triangular; (2) sinusoidal; (3) exponential; and (4) rectangular. In the first part of the article, a numerical study is presented of the effect of oscillator period, strength, damping, post‐yielding stiffness and number of excitation cycles, on inelastic response. Results are presented in the form of dimensionless graphs and regression formulas that elucidate the salient features of the problem. It is shown that conventional R–µ relations may significantly underestimate ductility demand imposed by near‐fault motions. The second part of the article concentrates on elastic‐perfectly plastic oscillators. Closed‐form solutions are derived for post‐yielding response and associated ductility demand. It is shown that all three ground motion histories (i.e. acceleration, velocity, and displacement) control oscillator response—contrary to the widespread view that ground velocity alone is of leading importance. The derived solutions provide insight on the physics of inelastic response, which is often obscured by the complexity of numerical algorithms and actual earthquake motions. The model is evaluated against numerical results from near‐field recordings. A case study is presented. Copyright © 2006 John Wiley & Sons, Ltd.     

3-D shell analysis of cylindrical underground structures under seismic shear (S) wave action

The 3-D shell theory is employed in order to provide a new perspective to earthquake-induced strains in long cylindrical underground structures, when soil-structure interaction can be ignored. In this way, it is possible to derive analytical expressions for the distribution along the cross-section of axial, hoop and shear strains and also proceed to their consistent superposition in order to obtain the corresponding principal and von Mises strains. The resulting analytical solutions are verified against the results of 3-D dynamic FEM analyses. Seismic design strains are consequently established after optimization of the analytical solutions against the random angles which define the direction of wave propagation relative to the longitudinal structure axis, the direction of particle motion and the location on the structure cross-section. The basic approach is demonstrated herein for harmonic shear (S) waves with plane front, propagating in a homogeneous half-space or in a two layer profile, where soft soil overlays the bedrock.     

Global force balance of region 1 current system

At times of strong solar wind forcing such as those that produce major magnetic storms, the region 1 current system dominates over the Chapman–Ferraro current system in mediating the transfer of force between the solar wind and the terrestrial system. The global force balance can be broken into two components, one involving the high-altitude part of the region 1 current system that is in contact with the solar wind (labeled here the HRS) and the other involving the low-altitude part of the region 1 current system that lies in the ionosphere (the LRS). Both communicate their J×B force to the geomagnetic dipole via a gradient in the magnetic field that they generate. In the HRS case the force acts to push the dipole away from the sun. This is the region 1 analog of the Chapman–Ferraro mechanism for transmitting the solar wind's force to the Earth. However, in the LRS case, the force (which is much stronger than in the HRS case) acts to push the dipole toward the sun, seemingly paradoxically. The LRS balances the ‘paradoxical’ sunward force on the dipole with an opposite force on the atmosphere. This paper uses MHD simulations to demonstrate the presence of both the normal force-transmitting gradient generated by the Chapman–Ferraro and the counter-Chapman–Ferraro gradient in the magnetic field generated by the region 1 current system.     

Small-scale seismic heterogeneity and mantle structure

George Helffrich discusses the links between seismology and geochemistry in the Bullerwell Lecture 2005, delivered at the annual meeting of the European Union of Geosciences in Vienna.