A Least-squares Approach to Depth Determination from Numerical Horizontal Self-potential Gradients

We have developed a least-squares minimization approach to depth determination of a buried ore deposit from numerical horizontal gradients obtained from self-potential (SP) data using filters of successive window lengths (graticule spacings). The problem of depth determination from SP gradients has been transformed into the problem of finding a solution to a nonlinear equation of the form f(z)=0. Formulas have been derived for vertical and horizontal cylinders and spheres. Procedures are also formulated to estimate the electrical dipole moment and the polarization angle. The method is applied to synthetic data with and without random noise. Finally, the validity of the method is tested on two field examples. In both cases, the depth obtained is found to be in a very good agreement with that obtained from drilling information.     

Determination of S-wave velocity structure using microtremors and spac method applied in Thessaloniki (Greece)

Array measurements of microtremors at 16 sites in the city of Thessaloniki were performed to estimate the V_s velocity of soil formations for site effect analysis. The spatial autocorrelation method was used to determine phase velocity dispersion curves in the frequency range from 0.8–1.5 to 6–7 Hz. A Rayleigh wave inversion technique (stochastic method) was subsequently applied to determine the Vs profiles at all the examined sites. The determination of V_s profiles reached a depth of 320 m. Comparisons with V_s values from cross-hole tests at the same sites proved the reliability of the SPAC method. The accuracy of the V_s profiles, the ability to reach large penetration depths in densely populated urban areas and its low cost compared to conventional geophysical prospecting, make Mictrotremor Exploration Method very attractive and useful for microzonation and site effects studies. An example of its application for the site characterization in Thessaloniki is presented herein.     

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.     

New palaeointensity results from Cretaceous basalt of Inner Mongolia, China

We present new Thellier-Thellier palaeointensity results from three cooling units (32 samples) of Inner Mongolia lava flows (91.7 Ma) emplaced during the Cretaceous Normal Superchron (CNS). Based on rock-magnetic and microscopy observations the magneto-mineralogy of all samples is determined to be primary and unaltered high-Ti titanomagnetite. Accepted palaeointensity determinations, obtained in the 80-200 °C temperature interval, are of good technical quality with positive standard partial thermoremanent magnetisation (pTRM) checks and pTRM-tail checks. Obtained palaeointensity estimates range from 14.7 to 28.0 μT, with virtual axial dipole moments (VADM) of 2.4 to 4.6 (× 10²² Am²). The data agree well with recently published results from the same region and, combining the two datasets, we obtain independent estimates from six different cooling units yielding a time-averaged VADM of 3.2 ± 1.6 (× 10²² Am²). These data suggest a relatively low dipole moment towards the end of the Cretaceous Normal Superchron.     

Toxic Cyanobacteria in Greek Freshwaters, 1987—2000: Occurrence,Toxicity, and Impacts in the Mediterranean Region

In a survey in Greece from 1987 to 2000 hepatotoxic cyanobacterial blooms were observed in 9 out of 33 freshwaters. Microcystins (MCYSTs) were detected by HPLC in 7 of these lakes, and the total MCYST concentration per scum dry weight ranged from 50.3 to 1638 ± 464 μg g^—1. Cyanobacterial genera (Microcystis, Anabaena, Anabaenopsis, Aphanizomenon, Cylindrospermopsis) with known toxin producing taxa were present in 31 freshwaters. From our data and a review of the literature, it would appear that Mediterranean countries are more likely 1) to have toxic cyanobacterial blooms consisting of Microcystis spp. and 2) to have higher intracellular MCYST concentrations. A case study in Lake Kastoria is used to highlight seasonal patterns of cyanobacterial and MCYST‐LR occurrence and to assess cyanotoxin risk. Cyanobacterial biovolume was high (> 11 μL L^—1) throughout the year and was in excess of Guidance Level 2 (10 μL L^—1) proposed by WHO for recreational waters and Alert Level 2 for drinking water. Further, surface water samples from April to November exceeded Guidance Level 3, with the potential for acute cyanobacterial poisoning. Intracellular MCYST‐LR concentrations (max 3186 μg L^—1) exceeded the WHO guideline for drinking water (1 μg L^—1) from September to November with a high risk of adverse health effects. Preliminary evidence indicates that in 3 lakes microcystins are accumulated in some aquatic organisms. Generally, a high risk level can be deduced from the data for the Mediterranean region.     

Bedload transport resistance in rough open‐channel flows

During bedload movement by saltation, streamwise momentum is transferred from the ?ow to the saltating grains. When the grains collide with other grains on the bed or in the ?ow, streamwise momentum is reduced, and there is a decrease in streamwise ?ow velocity and an increase in ?ow resistance, herein termed bedload transport resistance f_bt. Based on experiments in two ?umes with ?xed and mobile plane beds and previously published data, an equation is developed that may be used to predict f_bt for both capacity and non‐capacity ?ows. The variables in this equation are identi?ed by dimensional analysis and the coef?cients are determined by non‐linear regression. This equation applies to rough turbulent open‐channel ?ows, where the relative submergence is between 1 and 20 and the entire sediment load moves by saltation. An investigation of the relative magnitudes of f_bt and grain resistance f_c suggests that where dimensionless shear stress θ is less than 1 and saltation is the dominant mode of bedload transport, f_bt/f_c increases with θ but never exceeds 1. Copyright © 2004 John Wiley & Sons, Ltd.     

Hydrological connectivity of soil pipes determined by ground‐penetrating radar tracer detection

Soil pipes are common and important features of many catchments, particularly in semi‐arid and humid areas, and can contribute a large proportion of runoff to river systems. They may also signi?cantly in?uence catchment sediment and solute yield. However, there are often problems in ?nding and de?ning soil pipe networks which are located deep below the surface. Ground‐penetrating radar (GPR) has been used for non‐destructive identi?cation and mapping of soil pipes in blanket peat catchments. While GPR can identify subsurface cavities, it cannot alone determine hydrological connectivity between one cavity and another. This paper presents results from an experiment to test the ability of GPR to establish hydrological connectivity between pipes through use of a tracer solution. Sodium chloride was injected into pipe cavities previously detected by the radar. The GPR was placed downslope of the injection points and positioned on the ground directly above detected soil pipes. The resultant radargrams showed signi?cant changes in re?ectance from some cavities and no change from others. Pipe waters were sampled in order to check the radar results. Changes in electrical conductivity of the pipe water could be detected by the GPR, without data post‐processing, when background levels were increased by more than approximately twofold. It was thus possible to rapidly determine hydrological connectivity of soil pipes within dense pipe networks across hillslopes without ground disturbance. It was also possible to remotely measure travel times through pipe systems; the passing of the salt wave below the GPR produced an easily detectable signal on the radargram which required no post‐processing. The technique should allow remote sensing of water sources and sinks for soil pipes below the surface. The improved understanding of ?owpath connectivity will be important for understanding water delivery, solutional and particulate denudation, and hydrological and geomorphological model development. Copyright © 2004 John Wiley & Sons, Ltd.     

A hybrid method for size characterization of coarse subsurface fluvial sediments

This paper introduces and evaluates a novel method for ascertaining the grain‐size distribution of subsurface sediments that involves profoundly less sampling effort than standard methods. It is based on hybrid sampling principles previously applied to the construction of synthetic surface grain‐size distributions. The method is developed from an empirical demonstration of the approximate similarity of surface and subsurface grain‐size distributions when compared over a common range of sizes. Subsurface hybrid models are found to provide good facsimiles of grain‐size distributions de?ned using standard criteria and to yield distribution percentiles with millimetre accuracy. The technique is presented as an expedient alternative to standard methods for large, perennial gravel‐bed rivers. As this is a new technique, prudent application is advised in lieu of further investigation. Copyright © 2004 John Wiley & Sons, Ltd.     

Comparison of hydraulic geometry between sand‐ and gravel‐bed rivers in relation to channel pattern discrimination

A comparison has been made between the hydraulic geometry of sand‐ and gravel‐bed rivers, based on data from alluvial rivers around the world. The results indicate a signi?cant difference in hydraulic geometry among sand‐ and gravel‐bed rivers with different channel patterns. On this basis, some diagrams for discrimination of meandering and braided channel patterns have been established. The relationships between channel width and water discharge, between channel depth and water discharge, between width–depth ratio and water discharge and between channel slope and water discharge can all be used for channel pattern discrimination. The relationship between channel width and channel depth can also be used for channel pattern discrimination. However, the accuracy of these relationships for channel pattern discrimination varies, and the depth–discharge relationship is a better discriminator of pattern type than the classic slope–discharge function. The cause for this difference has been explained qualitatively. To predict the development of channel patterns under different natural conditions, the pattern discriminator should be searched on the basis of independent or at least semi‐independent variables. The relationship between stream power and bed material grain size can be used to discriminate channel patterns, which shows a better result than the discriminator using the slope–discharge relationship. Copyright © 2004 John Wiley & Sons, Ltd.