The use of GPR and VES in delineating a contamination plume in a landfill site: a case study in SE Brazil

This paper presents the results of the application of the Ground Penetrating Radar (GPR) method, or Georadar, in outlining a zone of contamination due to solid residues at the waste burial site of Rio Claro in the state of São Paulo, SE Brazil. A total of eight GPR profiles with 50- and 100-MHz antennae were surveyed. Six profiles were located within the landfill site and the remaining two were outside. The main objective of the GPR survey was to evaluate the side extension of contamination. A Vertical Electric Sounding (VES) survey was performed at four points within the site in order to investigate the vertical extent of the contamination plume and to define the bottom of the landfill. Two additional VESs were done outside the landfill with the purpose of determining the top of the ground water table and the geoelectric stratigraphy of the background. From the interpretation of the GPR profiles, it was possible to locate the top of the contamination plume and to infer that it was migrating laterally beyond the limits of the waste disposal site. This was observed along the profile situated close to the highway SP-127, which was about 20 m from the limit of the site. The signature of the contaminant appears as a discontinuous reflector that is believed to be a shallow ground water table. The discontinuity is marked by a shadow zone, which is characteristic of conductive contaminant residues. The contamination did not move far enough to reach a sugar cane plantation located at approximately 100 m from the border of the site. In the regions free from contamination, the ground water table was mapped at approximately 10 m of depth, and it was characterized by a strong and continuous reflector. The radar signal penetrated deep enough and enabled the identification of a second reflector at approximately 14 m deep, interpreted as the contact between the Rio Claro and the Corumbataí formations. The contact is marked by the presence of gravel characterized by ferruginous concretes, which cause the strong amplitude reflection in the GPR profile. Within the landfill site, the quantitative interpretation of the VES results showed the contamination zone. The base of the landfill varies between 11 and 15 m deep. Outside the landfill site, the VES results showed no indication of contamination and allowed the determination of the top of the ground water table and the contact between the Rio Claro and the Corumbataí formations. The results of GPR and VES showed a good agreement and the integrated interpretations were supported by local geology and information from several boreholes, about 17 m depth, on average. The bottom of the landfill reaches a maximum of 14.5 m depth.     

SH Wave Propagation in Viscoelastic Media

When comparing solutions for the propagation of SH waves in plane parallel layered elastic and viscoelastic (anelastic) media, one of the first things that becomes apparent is that in the elastic case the location of the saddle points required to obtain a high frequency approximation are located on the real p axis. This is true of the branch points also. In a viscoelastic medium this is not typical. The saddle point corresponding to an arrival lies in the first quadrant of the complex p-plane as do the branch points. Additionally, in the elastic case the saddle point and branch points lie on a straight line drawn through the origin (the positive real axis in the complex p-plane), while in the viscoelastic case this is generally not the case and the saddle point and branch points lie in such a manner as to indicate the degree of their complex values.In this paper simple SH reflected and transmitted particle displacement arrivals due to a point torque source at the surface in a viscoelastic medium composed of a layer over a half space will be considered. The path of steepest descent defining the saddle point in the first quadrant will be parameterized in terms of a real variable and the high frequency solutions and intermediate analytic results obtained will be used to formulate more specific constraints and observations regarding saddle point location relative to branch point locations in the complex p-plane.As saddle point determination for an arrival is, in general, the solution of a non-linear equation in two unknowns (the real and imaginary parts of the complex saddle point p ₀), which must be solved numerically, the use of analytical methods for investigating this problem type is somewhat limited.Numerical experimentation using well documented solution methods, such as Newton's method, was undertaken and some observations were made. Although fairly basic, they did provide for the design of algorithms for the computation of synthetic traces that displayed more efficient convergence and accuracy than those previously employed. This was the primary motivation for this work and the results from the SH problem may be used with minimal modifications to address the more complicated subject of coupled P-SV wave propagation in viscoelastic media.Another reason for revisiting a problem that has received some attention in the literature was to approach it in a fairly comprehensive manner so that a number of specific observations may be made regarding the location of the saddle point in the complex p-plane and to incorporate these into computer software. These have been found to result in more efficient algorithms for the SH wave propagation and a significant enhancement of the comparable software in the P-SV problem.     

Seismic reflection imaging of an oceanic detachment fault: Atlantis megamullion (Mid-Atlantic Ridge, 30°10′N)

We present multichannel seismic reflection data collected over the Atlantis megamullion, at the eastern ridge-transform intersection of Atlantis fracture zone on the northern Mid-Atlantic Ridge, and over its conjugate crust. These data image for the first time the internal structure of a young, well-developed megamullion dome formed by tectonic extension across a long-lived oceanic detachment fault. The exposed, corrugated detachment-fault surface exhibits a sharp, coherent reflection that contrasts with less organized reflectivity of surrounding basaltic seafloor. At the termination of the megamullion the fault is imaged ∼13 km along-strike beneath a volcanic hanging-wall block at a sub-seafloor depth of 0.2-0.5 s two-way travel time, reaching north as far as 30°19′N. The eastward dipping of the fault beneath the hanging-wall block is estimated to be ∼6-14°. The corrugated fault surface is underlain by a continuous, strong, and relatively smooth reflection (D) at 0.2-0.25 s sub-bottom below the central axis of the dome. This reflection deepens up to 0.6 s sub-bottom beneath the western slope and it appears to intersect the seafloor on the eastern slope. We suggest that Atlantis massif formed by sequential slip on two different detachment faults that merged at depth, with breakaways as little as ∼2 km apart. The initial detachment is represented by reflection D, and the second corresponds to the presently exposed fault surface. In this interpretation, much of the sliver between the faults is interpreted to be strongly serpentinized peridotite with reduced seismic velocity; it lies in contact with less altered, higher-velocity mantle below the first detachment, resulting in the strong, smooth character of reflection D. Mantle rocks exposed in the megamullion indicate that the feature formed during a period of extreme tectonic extension and probably limited magmatism. In conjugate crust corresponding to termination of the megamullion, observed sub-bottom reflections are interpreted as base of seismic layer 2A. This layer is as thick as or thicker (∼570-900 m) than layer 2A in normal Atlantic crust, and it suggests that relatively normal crustal accretion occurred by the time the megamullion stopped forming.     

Multicomponent ionic diffusion in porewaters: Coulombic effects revisited

The diffusion of an ion in porewaters cannot occur independently of the other ions in solution as a result of Coulombic coupling, as well as from other effects not considered here. Unfortunately, a longstanding disagreement exists about the correct form and meaning of the equations that describe Coulombic coupling in porewaters, i.e., Ben-Yaakov [Am. J. Sci. 281 (1981) 974] vs. Lasaga [Am. J. Sci. 281 (1981) 981]. This paper re-examines this controversy by reformulating the problem starting from fundamental concepts of mass and charge conservation. We show that these antagonistic formulations are both valid and, in fact, equivalent, when the different interpretations of charge balance are resolved. Most of the disagreements between Ben-Yaakov and Lasaga are then shown to result from differing methods of solution, not fundamental disparities in their models. We note, however, that the explanation for the concept of “stationary” gradients of nonreacting ions as given Ben-Yaakov is inaccurate, and such gradients do lead to diffusive fluxes that are counterbalanced by electrochemical migrational fluxes to produce no net flux (excluding advective flux). We further find that the bicarbonate diffusive flux will not balance the diffusional charge flux of sulfate during its reduction if advection is present. This latter imbalance generates compensating fluxes in the other nonreacting ions. We have applied our theory to a simplified case of sulfate reduction in a marine sediment. The results show that nonreacting ions do diffuse and that with normally expected values of porewater advection, the ratio of the bicarbonate to the sulfate flux can be far different than the ideal value of −2.     

Approximate modal decomposition of inelastic dynamic responses of wall buildings

Two approximate methods for decomposing complicated inelastic dynamic responses of wall buildings into simple modal responses are presented. Both methods are based on the equivalent linear concept, where a non‐linear structure is represented by a set of equivalent linear models. One linear model is used for representing only one vibration mode of the non‐linear structure, and its equivalent linear parameters are identified from the inelastic response time histories by using a numerical optimizer. Several theoretical relations essential for the modal decomposition are derived under the framework of complex modal analysis. Various numerical examinations have been carried out to check the validity of the proposed modal decomposition methods, and the results are quite satisfactory in all cases. Fluctuating bending moment and shear at any location along the wall height contributed by each individual vibration mode can be obtained. Modal contributions to shear and flexural strength demands, as well as the corresponding modal properties, under various seismic loading conditions can also be identified and examined in detail. Furthermore, the effects of higher vibration modes on seismic demands of wall buildings are investigated by using the modal decomposition methods. Several new insights into the complicated inelastic dynamics of multi‐story wall buildings are presented. Copyright © 2004 John Wiley & Sons, Ltd.     

Land cover dynamics in Garhwal Himalayas — a case study of balkhila sub-watershed

The present paper focuses on land cover dynamics pattern in Balkhila sub-watershed situated in Garhwal Himalayas. This study contributes example from human shaped ecosystem in mountainous regions where fragmentation of natural resources is active. The remote sensing and GIS has been used to understand the land cover dynamics along the topography. The results report that the land cover dynamics is dependent on the aspect due to sun illumination. The altitude and slope are no more a barrier for resource extraction and the human activity zone is shifting towards higher altitudes and slopes. The changes are also defined along the road and settlements.     

An automated seamless mosaicing system of multi-charge coupled devices of panchromatic data

Panchromatic data of pixel resolution 5.8 m obtained from IRS-1C and IRS-1D satellites proved to be very useful for mapping purposes. One of the popular data product is the 70 km swath mosaic which is covered by a combination of 3 CCD line sensors, each with 4096 pixels. Each CCD-line sensor with different imaging times causes geometric problems of mosaicing three strips data together. In this paper, we propose the details of the design elements of system that caters to the need for accurate and automatic multi strip image registration without any second resampling of the data. The systematic geometric correction grid mapping is improved to facilitate accurate mosaicing by automatic image registration task that makes use of the overlap data within image strips and image registration is achieved up to sub-pixel level.     

A fast intensity-hue-saturation fusion technique with spectral adjustment for IKONOS imagery

Among various image fusion methods, intensity-hue-saturation (IHS) technique is capable of quickly merging the massive volumes of data. For IKONOS imagery, IHS can yield satisfactory "spatial" enhancement but may introduce "spectral" distortion, appearing as a change in colors between compositions of resampled and fused multispectral bands. To solve this problem, a fast IHS fusion technique with spectral adjustment is presented. The experimental results demonstrate that the proposed approach can provide better performance than the original IHS method, both in processing speed and image quality.