Visualization of Volcanic Rock Geochemical Data and Classification with Artificial Neural Networks

An unsupervised neural network technique, Growing Cell Structures (GCS) was used to visualize geochemical differences between four different island arc volcanic rock types: basalts, andesites, dacites and rhyolites. The output of the method shows the cluster structure of the dataset clearly, and the relevant geochemical patterns and relationships between its variables. The data can be separated into four clusters, each associated with a specific volcanic rock type (basalt, andesite, dacite and rhyolite), according to a unique combination of major element concentrations. Following clustering, performance of the trained GCS network as a classifier of volcanic rock type was tested using two test datasets with major element concentration data for 312 and 496 island arc volcanic rock samples of known volcanic type. Preliminary classification results are promising. In the first test dataset 94% of basalts, 76% of andesites, 83% of dacites and 100% of the rhyolites were classified correctly. Successful classification rates in the second dataset were 100%, 80%, 77%, and 98% respectively. The success of the analysis suggests that neural networks analysis constitutes a useful analytical tool for identification of natural clusters and examination of the relationships between numeric variables in large datasets, and that can be used for automatic classification of new data.     

Hydraulic evaluation of aquifer storage and recovery (ASR) with urban stormwater in a brackish limestone aquifer

A 5-year aquifer storage and recovery trial at Andrews Farm in South Australia involving the injection of more than 250 ML (250,000 m³) of fresh but turbid stormwater into a brackish limestone aquifer over 4 years and recovery of 150 ML in the fifth provided the opportunity to evaluate rates of clogging and unclogging and the potential to recover water suitable for irrigation supplies. Results reveal there is some clogging by injected sediment, but only to a relatively small degree considering the high suspended solid concentrations and moderate aquifer transmissivity. This clogging was offset by increased matrix porosity through calcite dissolution and by routine well redevelopments after each 40 ML of injection. No significant microbial clogging occurred. Breakthrough responses at three observation wells and the proportion of injectant in the recovered water were determined from chloride data. Temperature and caliper profiles clearly indicate the heterogeneous nature of the aquifer that is attributed, in part, to sand removal during the initial well development. The recovery efficiency was greater than 60%. The trial demonstrates that urban stormwater containing high and variable particulate levels, which receives only passive pre-treatment and is not disinfected, can be used to freshen a heterogeneous brackish aquifer to create a useful water resource.     

Effect of Spectral Resolution on the Mg II Index as a Measure of Solar Variability

The solar Mgii core-to-wing ratio is a useful index of UV variability throughout the solar cycle because it has been measured since 1978 in a series of successive satellite missions: Nimbus 7, Solar Mesosphere Explorer (SME), the NOAA 9–14 series, Upper Atmosphere Research Satellite (UARS), and ERS-2. Eventual construction of a single time series from 1978 to the present by combining these measurements will give a long record of almost daily UV variability to serve as a surrogate for estimating both UV and EUV solar radiation. Here we address the effect of spectral resolution on determination of both long-term and short-term solar variability from this index. We use UARS/SOLSTICE measurements of the Mgii line from October 1991 to December 1996 to study the effect of two spectral resolution regimes characteristic of existing measurements, 0.20 to 0.25 nm and 1.10 to 1.15 nm, on determination of the amplitude of 27-day rotational modulation and the more gradual change in chromospheric radiation in the declining phase of solar cycle 22. The two Mgii indices give solar variations that differ by a scaling factor of 2× for both the solar cycle change from 1992 to 1997 and the amplitude of 27-day modulation over the same period. Both types of measurements appear to yield solar signal equally well except at solar minimum when the solar changes become quite small.     

A generalized genetic framework for the development of sinkholes and Karst in Florida,U.S.A.

Karst topography in Florida is developed on the Tertiary limestones of the Floridan aquifer Post-depositional diagenesis and solution have made these limestones highly permeable, T=ca. 50,000 m²/d. Zones of megaporosity have formed at unconformities, and dissolution has enlarged joints and fractures Erosion of the overlying clastic Miocene Hawthorn group strata on one flank of a structural arch has exposed the limestone The elevated edge of the Hawthorn cover forms the Cody scarp Ubiquitous solution pipes have previously formed at joint intersections and are now filled Downwashing of the fill deeper into solution cavities in the limestone and subsidence of the overlying unconsolidated sediments causes surface collapse a subsidence doline or sinkhole This process may penetrate up to 60 m of the semi-consolidated Hawthorn cover, as occurred when the Winter Park sinkhole developed Dense clusters of solution pipes may have formed cenotes which are now found on the exposed limestone terrain Groundwater moves laterally as diffuse flow except where input or outflow is concentrated. At sinking streams, vertical shafts, and springs, karst caves have formed, but only the major sinking streams form through-flowing conduit systems Shaft recharge dissipates diffusely. Spring discharge is concentrated from diffuse flow In both cases, conduits taper and merge into a zone of megaporosity     

Effects of a variety of impacts on seagrass restoration in Jamaica

For the first time seagrass rehabilitation was attempted in a tropical nation, Jamaica. Repairs to a variety of developmental impact types were attempted with test plots: dredge and fill, bauxite and oil spills, channelization, urban run-off, urban sewage, river bulkheading, jetty construction and artificial beach, industrial lagoons, saline lagoons, thermal effluents, and cement tailings. Three species were routinely planted (Thalassia testudinum, Halodule wrightii, Syringodium filiforme) plus Ruppia maritima in saline lagoons, each by two planting techniques. Types of damage which were successfully attempted to be restored for the first time were from heated wastes, bauxite spills, jettied river mouths, small scale oil spills, and diluted urban wastes. Types of pollution damage previously repaired in subtropics and temperate zones, which could be successfully repaired in the tropics also were dredge and fill and urban run-off. Erosion seagrass restoration repair was partially successful at some sites. Attempts to repair high salinity lagoons (up to 160‰) and cement tailings were unsuccessful. Functional species were found which could tolerate ongoing pollutant loads in some cases.     

The effect of a shallow low viscosity zone on the apparent compensation of mid-plate swells

A simple model for mid-plate swells is that of convection in a fluid which has a low viscosity layer lying between a rigid bed and a constant viscosity region. Finite element calculations have been used to determine the effects of the viscosity contrast, the layer thickness and the Rayleigh number on the flow and on the perceived compensation mechanism for the resulting topographic swell. As the viscosity decreases in the low viscosity zone, the effective local Rayleigh number for the top boundary layer of the convecting cell increases. Also, because the lower viscosity facilitates greater velocities in the low viscosity zone, the low viscosity layer produces proportionally greater horizontal flow near the conducting lid, causing the base of the conducting lid to appear like a free boundary. The change in the local Rayleigh number and in the effective boundary condition both cause the top boundary layer to thin. Through a Green's function analysis, we have found that the low viscosity zone damps the response of the surface topography to the temperature anomalies at depth, whereas it causes the gravity and geoid response functions to change sign at depth counteracting the positive contributions from the shallower temperature variations. By increasing the viscosity contrast, the conbined effects of the thinning of the boundary layer and the behaviour of the response functions allow the apparent depth of compensation to become arbitrarily small. Therefore, shallow depths of compensation cannot be used to argue against dynamic support of mid-plate swells. Furthermore, we compared the distribution of the effective compensating densities, which is used to obtain the geoid, to that of Pratt compensation, which is often used to calculate the depth of compensation from geoid and topography data for mid-plate swells. For all of our calculations including those with no low viscosity layer, the effective gravitational mass distribution is more complex than assumed in simple Pratt models, so that the Pratt models are not an appropriate gauge of the compensation mechanism.     

Close interval sampling of interstitial silicate and porosity in marine sediments

A method was developed to measure porosity and dissolved interstitial silicate at millimeter intervals or less in a sediment core. In cores from Emerald Basin (Scotian Shelf), interstitial concentrations near the sediment surface did not drop rapidly to bottom-water concentrations as measured in bottle casts (28 μM) but remained as high as 166 μM in the upper 0.5 mm of sediment High rates of benthic silicate release were measured which could not be accounted for by interstitial concentration gradients or by ventilation of macro-invertebrate burrows. The silicate discontinuity observed between the sediments and water column suggests that a diffusive sublayer exists in a zone of viscous flow above the sediment surface. This is possible only if a surface reaction is primarily responsible for silicate release. By assuming a linear concentration gradient across this diffusive sublayer, the silicate release rates were used to estimate the thickness of the sublayer to be about 2 mm.     

The Petrology and Geochemistry of Cone-sheets from the Cuillin Igneous Complex, Isle of Skye: Evidence for Combined Assimilation and Fractional Crystallization during Lithospheric Extension

The origin and evolution of the magma(s) involved in the formationof the olivine tholeiite cone-sheets which invade the 60-MaCuillin Igneous Complex, Isle of Skye, have been investigatedusing mineral chemistry data together with whole-rock major-and trace-element and isotope (Sr and Nd) geochemistry data.The most primitive compositions identified are almost identicalto those of the basalts being erupted at the present day alongthe spreading centre in Iceland. Rare examples of lavas fromthe slightly older lava pile of west-central and north Skye,together with a large number of dykes from the axial zone ofthe Skye Main Swarm, are of similar composition and are referredto as being of the Preshal More Basalt type (spelling is accordingto current Ordnance Survey of Great Britain maps). The intra-suitecompositional variation exhibited by the cone-sheets can beexplained in terms of relatively low-pressure fractionationof the three-phase assemblage olivine+clinopyroxene+plagioclasein their cotectic proportions of 10: 35: 55. Modelling of major-and compatible trace-element data indicates that the most evolvedcomposition may be derived by 60% crystallization of the leastevolved composition. Incompatible trace-element modelling impliesgreater degrees of crystal-liquid fractionation: Y and Zr indicate64 and 68% crystallization, respectively, whereas the rare earthelements (Eu, Yb, Gd, and Ce) indicate between 71 and 75% crystallization.This decoupled behaviour of compatible and incompatible elementsis attributed to the cone-sheet magma(s) evolving in a subjacentmagma chamber, before final emplacement in the overlying crust.Sr- and Nd-isotope data indicate that the cone-sheet magma(s)assimilated upper-crustal lithologies during fractional crystallization,most likely involving amphibolite facies gneiss of the LewisianGneiss Complex. This style of simultaneous assimilation andfractional crystallization (AFC) process in the upper crustfor the Skye cone-sheet magma(s) is in direct contrast to theprocesses identified for the magmas which produced the bulkof the lava field of west-central and north Skye, where assimilationoccurred after fractionation. The cone-sheet magmas were generatedby relatively large degrees of partial melting of a depletedmantle source associated with significant lithospheric stretching.The cause of this depletion, together with the temporal relationshipsbetween the cone-sheet magmas and the dominant transitionalmagmas of the Skye lava field and sill complex, are also discussed. ^* Reprint request to B. R. Bell