TRANSIENT EM ANALOGUE MODELLING FOR KOREAN TREASURE HUNTING WITH THE SIROTEM SYSTEM1

The inductive transient electromagnetic method (TEM) shows good potential for the detection of metallic relics of historical interest, such as buried antique bells, Buddhist idols, or precious metal nuggets. The effectiveness of the method was investigated with analogue models, using transmitter-receiver loop configurations with sizes or receiver spacings slightly different from those usually applied in earth resources exploration with the Sirotem system. The analogue modelling results show that the location and depth of the buried treasure may be obtained from the Sirotem data. A solid metallic object such as an antique bell could be detected to a depth about 12 times greater than its size, and treasure consisting of separate metallic objects, such as gold nuggets, could be detected to a depth 5 times its linear dimensions with a transmitter current of ? 20 A. In scaling down the dimensions of a target, its conductivity should be increased in order to preserve the same TEM conditions found in the field. However, since the buried treasure consists of gold or copper objects, it is not possible to properly scale the conductivity. Hence, in the field, the depth detection limit is expected to be greater than that derived from analogue modelling.     

Lightning fatalities in Australia, 1824–1991

Information on fatalities from lightning strikes has been extracted from a specially compiled database on natural hazards in Australia. Records dating from 1803–1991 indicate that at least 650 persons have been killed by lightning strikes. Maps and charts show the percentages of victims with respect to age, sex, locality, activity, and other circumstances of the strike. The majority of the 650 fatalities recorded have occurred along the more populous southeastern coast. The overall death rate, from 1910–1989, is 0.08 per 100 000 population. The annual number of lightning fatalities has decreased with time, from a death rate of 0.21 in 1910–1919 to 0.01 in 1980–1989. This trend is more pronounced when population figures are taken into account. The group that has been most at risk in Australia is that of males aged 15–19, followed by males aged 20–34. The male:female ratio of victims has decreased with time but is not approaching equality, being 11.6 in 1910–1919 and 5.3 in 1980–1989. The diurnal and monthly occurrences of lightning fatalities peak at 12.00–18.00 hours and November-February respectively. About 86% of fatalities have occurred outdoors and 14% have occurred indoors. Approximately three-fifths of fatalities have been work-related, and the group of workers that has traditionally been most at risk is that of land-workers. Approximately one-fifth of fatalities have been recreation-related, although this proportion has been increasing with time. The recreational activities of water sports, golf, and cricket have had the greatest number of lightning fatalities. Comparisons are made between these data and the results of other similar studies, both in Australia and overseas.     

Ripple formation in combined transdirectional steady and oscillatory flow

Measurements are described of the geometry of ripples formed on beds of sand exposed to a steady current at right angles to an oscillatory flow. Four different sands were studied. The oscillation was produced by an oscillating tray set into the bed of a steady-flow flume. It was observed that straight-crested ripples formed by oscillatory flow would usually develop a ‘serpentine’ form when the superimposed steady current exceeded a certain limit. For amplitudes of the tray velocity U_∞ less than about 0.38 m s^-1 this limit corresponded to U_∞/ū_*c>31, where ū_*c is the shear velocity measured just upstream of the oscillating tray. It is suggested that the serpentine form is caused by the interaction of vortices carried back and forth between adjacent ripples. On this assumption, the wavelength of the serpentine form would be proportional to the product of period of oscillation and near-bed steady current velocity. The present measurements appear to support this hypothesis although there is also evidence that the wavelength is influenced by preferred spacing patterns between vortices. The measurements also show the ratio of the amplitude of the serpentine form to its wavelength to be approximately constant. Empirical relationships are derived relating ripple geometry to flow and sediment properties. It is observed that the influence of Reynolds number and sediment properties on the geometry is very weak. It is suggested that this is typical of ripples formed with relatively low sediment transport rates. It is also found that, under the present experimental conditions, the ripple spacing in the direction of oscillation is almost independent of the magnitude of the steady current and in close agreement with the wavelengths previously measured in an oscillating water tunnel. This suggests that the additional inertia effects associated with oscillating tray rigs were not sufficient to affect bed geometry under the present test conditions.     

Effects of diffusional modification of garnet growth zoning on P-T path calculations

The effect of intragranular diffusion on chemical zoning in garnet and on P-T paths calculated from that zoning was evaluated using a numerical model of multicomponent diffusion in combination with simulations of garnet growth. Syn-and post-growth diffusion of Mg–Fe–Mn–Ca species in garnet was calculated for a model pelitic assemblage over a range of temperatures from 485 to 635°C. Compositions from zoned garnet, as modified by diffusion, hypothetical inclusions of plagioclase within garnet and matrix phases were used to reconstruct pressure-temperature (P-T) paths from isobaric and polybaric model histories. P-T path calculations, based on numerical simulations conducted over an input isobaric heating path that reached peak temperatures between 585 and 635°C, show that relaxation of garnet compositional gradients by diffusion can induce modest to appreciable curvature in the inferred paths. Retrieved paths also indicated somewhat smaller overall temperature changes relative to the actual temperature difference of the input path. The magnitude of these distortions is shown to depend upon the heating and cooling rate and garnet crystal size as well as the actual peak temperature condition. The effect of diffusion on path trajectories in simulations with thermal histories that also included cooling were comparable to heating-only models that reached peak temperatures approximately 15–30°C higher. Compositions of garnets with radii less than 1 mm, that reached actual peak temperatures of 605°C along temperature-time histories characteristic of regional metamorphism, experienced sufficient diffusional relaxation to introduce errors of hundreds of bars to in excess of one kilobar in path trajectories. Path distortions were significant at heating/cooling rates less than 10°C/Ma, but rapidly diminished for rates faster than this. In polybaric simulations diffusion effects were least noticeable when the actual pressure-temperature conditions changed in a clockwise sense (i.e., convex to higher P and higher T), but apprecciable modification was seen in path models that underwent counterclockwise changes in P and T. Reequilibration of garnet rim compositions occurred during cooling on all paths, and temperature maxima obtained from garnet-biotite geothermometry underestimated actual peak conditions by 40 to 70°C. Calculations suggest that P-T path trajectories calculated from garnets of at least 1 mm size, and that experienced actual thermal maxima below 585°C, are not likely to be distorted by diffusional effects during regional metamorphism. However, P-T path reconstructions based on garnet zonation with smaller grains or higher temperatures may lead to misinterpretation of crystallization history. The partitioning record of peak metamorphic temperatures may be destroyed by diffusional reequilibration of garnet rim compositions during cooling, seriously complicating the task of quantitatively estimating diffusion effects on path calculations.     

INTERPRETATION OF VLF-EM IN-PHASE DATA USING CURRENT DENSITY PSEUDOSECTIONS1

The interpretation of VLF-EM surveys in terms of buried conductors can be assisted by the application of a linear filter to the observed in-phase component of the vertical magnetic field. One such scheme is examined critically by using the calculated response from a variety of synthetic models to compute theoretical current density pseudosections. The results confirm that this filter technique provides a useful complementary tool for studying the third (i.e. depth) dimension. For single, steeply dipping plates (> 45°) diagnostic information may be derived concerning the depth, size, lateral location, and direction of dip, since the current density maxima seems always to occur within the conductor or at least within one data interval. However, there are some limitations which do not appear to be widely recognized. For single plates the angle of dip cannot be resolved by the current density transformation. Pattern distortions can occur where targets are in close proximity, where the cross-sectional form of the conductor is complex or where the dip is shallow (< 45°). In these latter cases the current density maxima may not occur within the conductive structure and therefore cannot be used to infer depth of burial and/or conductor shape with the same degree of confidence.     

The agricultural utility of lithic-mulch gardens: Past and present

The mulching of agricultural fields and gardens with stones, pebbles, volcanic ash and similar lithic materials is a variant agricultural strategy uniquely suited to the constraints of a dryland environment. While the agricultural literature contains reports on limited experiments with stone or gravel mulch, data regarding the employment of lithic mulch are lacking. This paper works to bridge this gap by comparing empirical data collected from prehistoric pebble-mulch gardens in the Galisteo Basin (New Mexico) to contemporary agricultural experiments with lithic mulch. Such past-present and experimental-real world comparisons allow us to more fully address the effectiveness of lithic mulching as an agricultural strategy. It is affirmed that lithic mulch is applied to garden plots, especially during periods of drought, in order to reduce soil erosion from wind and water, increase soil temperature to extend the growing season, increase water infiltration, and reduce the evaporative loss of water from wind and sun. The cumulative effect of this soil and water conserving technique is an increase in crop biomass and crop yield.     

The effect of ionic interactions on the oxidation of metals in natural waters

The effect of ionic interactions of the major components of natural waters on the oxidation of Cu(I) and Fe(II) has been examined. The various ion pairs of these metals have been shown to have different rates of oxidation. For Fe(II), the chloride and sulfate ion pairs are not easily oxidized. The measured decrease in the rate constant at a fixed pH in chloride and sulfate solutions agrees very well with the values predicted. The effect of pH (6 to 8) on the oxidation of Fe(II) in water and seawater have been shown to follow the rate equation

$\text{-d in [Fe(II)]/dt = k}{}_1\text{β}{}_1\text{α}{}_{\mathrm{Fe}}\text{/[H}{}^{\mathrm{+}}\text{] + k}{}_2\text{β}{}_2\text{α}{}_{\mathrm{Fe}}\text{/[H}{}^{\mathrm{+}}\text{]}{}^2$

where k₁ and k₂ are the pseudo first order rate constants, β₁ and β₂ are the hydrolysis constants for Fe(OH)⁺ and Fe(OH)⁰. The value of α_FE is the fraction of free Fe²⁺. The value of k₁ (2.0 ±0.5 min^?1) in water and seawater are similar within experimental error. The value of k₂ (1.2 × 10⁵ min^?1) in seawater is 28% of its value in water in reasonable agreement with predictions using an ion pairing model.For the oxidation of Cu(I) a rate equation of the form

$\text{?d ln [Cu(I)]/dt = k}{}_0\text{α}{}_{\mathrm{Cu}}\text{+ k}{}_1\text{β}{}_1\text{α}{}_{\mathrm{Cu}}\text{[Cl]}$

was found where k₀ (14.1 sec^?1) and k₁ (3.9 sec^?1) are the pseudo first order rate constants for the oxidation of Cu⁺ and CuCl⁰, β₁ is the formation constant for CuCl⁰ and α_Cu is the fraction of free Cu⁺. Thus, unlike the results for Fe(II), Cu(I) chloride complexes have measurable rates of oxidation.     

Particulate pollution monitoring in the Elefsis Gulf: The role of mineral magnetic studies

Measurements of ‘quadrature’ susceptibility (X_q) at low and high frequency on seabottom sediments of the Gulf of Elefsis supplement the results derived from the study of magnetic susceptibility (X) and saturation isothermal remanent magnetization (SIRM) (Scoullos et al., 1979) and make possible the distinction between the background soil erosional fraction and that of anthropogenic combustion origin. New devices speed up the on-site appraisal of the impact of discharges on marine sediment quality.     

Thermal histories of the core and mantle

Recognition that the cooling of the core is accomplished by conduction of heat into a thermal boundary layer (D″) at the base of the mantle, partly decouples calculations of the thermal histories of the core and mantle. Both are controlled by the temperature-dependent rheology of the mantle, but in different ways. Thermal parameters of the Earth are more tightly constrained than hitherto by demanding that they satisfy both core and mantle histories. We require evolution from an early state, in which the temperatures of the top of the core and the base of the mantle were both very close to the mantle solidus, to the present state in which a temperature increment, estimated to be ～ 800 K, has developed across D″. The thermal history is not very dependent upon the assumption of Newtonian or non-Newtonian mantle rheology. The thermal boundary layer at the base of the mantle (i.e., D″) developed within the first few hundred million years and the temperature increment across it is still increasing slowly. In our preferred model the present temperature at the top of the core is 3800 K and the mantle temperature, extrapolated to the core boundary without the thermal boundary layer, is 3000 K. The mantle solidus is 3860 K. These temperatures could be varied within quite wide limits without seriously affecting our conclusions. Core gravitational energy release is found to have been remarkably constant at ～ 3 × 10¹¹ W. nearly 20% of the core heat flux, for the past 3 × 10⁹ y, although the total terrestrial heat flux has decreased by a factor of 2 or 3 in that time. This gravitational energy can power the “chemical” dynamo in spite of a core heat flux that is less than that required by conduction down an adiabatic gradient in the outer core; part of the gravitational energy is used to redistribute the excess heat back into the core, leaving 1.8 × 10¹¹ W to drive the dynamo. At no time was the dynamo thermally driven and the present radioactive heating in the core is negligibly small. The dynamo can persist indefinitely into the future; available power 10¹⁰ y from now is estimated to be 0.3 × 10¹¹ W if linear mantle rheology is assumed or more if mantle rheology is non-linear. The assumption that the gravitational constant decreases with time imposes an implausible rate of decrease in dynamo energy. With conventional thermodynamics it also requires radiogenic heating of the mantle considerably in excess of the likely content of radioactive elements.