Identification of buried sinkholes using refraction tomography at Ft. Campbell Army Airfield,Kentucky

Karst aquifers are highly susceptible to contamination, with numerous points of entry for contaminants through recharge features such as sinkholes, swallow holes and solutionally enlarged fractures. These recharge features may be filled or obscured at the surface, requiring the use of geophysical or remote sensing techniques for their identification. This study uses seismic refraction data collected at the Ft. Campbell Army Airfield (CAAF), Kentucky, USA, to test the hypothesis that refraction tomography is a useful tool for imaging bedrock depressions beneath thick overburden (greater than 20 m of unconsolidated sediment). Southeast of the main taxiway of CAAF seismic velocity tomograms imaged a bedrock low, possibly a closed depression, at a depth of 25 m that had been earlier identified through delay-time analysis of the same refraction data. Tomography suggests the bedrock low is about 250-m wide by 10-m deep at its widest point. High rates of contaminant vapor extraction over the western extension of this feature suggest a high concentration of contaminants above, and within, this filled bedrock low, the base of which may contain solutionally enlarged fractures (i.e. karst conduits) that could funnel these contaminants to the upper or lower bedrock aquifers. This study thus demonstrates the viability of seismic refraction tomography as a tool for identification of filled sinkholes and bedrock depressions in karst areas.     

Additive versus multiplicative models in ecologic regression

Much research in environmental epidemiology relies on aggregate-level information on exposure to potentially toxic substances and on relevant covariates. We compare the use of additive (linear) and multiplicative (log-linear) regression models for the analysis of such data. We illustrate how both additive and multiplicative models can be fit to aggregate-level data sets in which disease incidence is the dependent variable, and contrast these results with similar models fitted to individual-level data. We find (1) that for aggregate-level data, multiplicative models are more likely than additive models to introduce bias into the estimation of rates, an effect not found with individual-level data; and (2) that under many circumstances multiplicative models reduce the precision of the estimates, an effect also not found in individual-level models. For both additive and multiplicative models of aggregate-level data, we find that, in the presence of covariates, narrow confidence interval are obtained only when two or more antecedent factors are strongly related to the measured covariate and/or the exposure of primary substantive interest. We conclude that the equivalency of fitting additive versus multiplicative models in studies with individual-level binary data does not carry over to studies that analyze aggregate-level information. For aggregate data, we strongly recommend use of additive models. Supported by Grant #1 U19 EH000102 from the National Center for Environmental Health, Centers for Disease Control and Prevention, Atlanta, GA.     

Rheology and geodynamic modelling: the next step forward

The application of continuum mechanics and microstructural analysis to geological studies over the past 30 years has spurred earth scientists to reassess fundamental tectonic processes such as subduction, collision and rifting in terms of dynamics. Armed with new analytical methods, geologists have returned to the field to look at rock structures with more mechanistic eyes. The advent of sophisticated computers, programs, and laboratory deformation equipment has facilitated the simulation of geodynamic processes that range in scale from the grain to the lithosphere. The result has been specialization, with the concomitant opening of communication gaps between geodynamicists, field geologists and rock mechanicists. Partly, these gaps reflect differences of perception and approach. In order to bridge these gaps, a workshop was organized after the DRM conference to debate how field and laboratory studies of deformed rocks can improve our understanding of lithospheric rheology, and in turn, how this understanding can be used to refine dynamic models of orogenesis. The workshop hosted participants with backgrounds in structural geology, experimental rock mechanics, metamorphic petrology and both numerical and analogue modelling. This paper summarizes the main controversies and conclusions reached during the workshop. For the sake of brevity, referencing in this summary is restricted to literature referred to during the oral presentations and to comments made by speakers themselves (names italicized).     

The Dublin Gulch intrusion-hosted gold deposit, Tombstone plutonic suite, Yukon Territory, Canada

The Dublin Gulch intrusion is a member of the Tombstone plutonic suite, a linear belt of middle Cretaceous intrusions that extend across the Yukon Territory. Like many of the intrusions in this suite, the Dublin Gulch intrusion is associated with several different zones of gold and tungsten mineralization, within and immediately adjacent to the intrusion. The Eagle zone (50.3 Mt @ 0.93 g/t gold), located in the southwestern part of the Dublin Gulch intrusion, hosts the most significant concentration of gold in the area. The gold occurs in a broadly east-west-striking, steeply south-dipping series of sheeted veins. The veins consist of early quartz-scheelite-pyrrhotite-pyrite-arsenopyrite, and are associated with K-feldspar-albite alteration envelopes. These grade out to and are overprinted by sericite-carbonate-chlorite alteration. The same assemblage also occurs in veinlets that refracture sheeted quartz veins and contain the majority of the gold. The gold occurs with molybdenite, lead-bismuth-antimony sulfosalts, galena, and bismuthinite. Gold correlates strongly with bismuth (r²=0.9), a relationship common to several intrusion-related gold deposits, but has a poor correlation with all other elements. Tungsten and molybdenum have a weak inter-element correlation (r²=0.55) and paragenetically pre-date the majority of gold precipitation. Lead, zinc, copper, silver, antimony, and arsenic have moderate to strong inter-element correlations (0.58 to 0.93). The change from tungsten-bearing mineralization through to gold-bismuth-rich ores with elevated syn- to post-ore lead, zinc, copper, silver, antimony, and arsenic can be grossly correlated with a change in hydrothermal fluid composition. Early scheelite-bearing quartz contains primary CO₂-rich fluid inclusions, which are post-dated by secondary inclusions with higher salinities (up to 15 wt% NaCl equiv.) and less CO₂. These latter inclusions are interpreted to coincide with the later gold-bismuth and base metal mineralization. The favored genetic model is one in which early CO₂-rich fluids exsolved from a magma with an initially high CO₂ content, but progressively became more saline and H₂O-rich as the system evolved.     

The beachwatch pollution monitoring programme: Changing priorities to recognize changed circumstances

Since the commissioning of the three major deepwater outfalls in Sydney, there has been a dramatic improvement in general beach water quality at Sydney's ocean beaches. This has required the NSW Environment Protection Authority's Beachwatch Programme to shift from assessing the impacts of gross sewage pollution due to nearshore outfalls to more diffuse pollution from stormwater and sewer overflows.

The visual indicators of gross sewage pollution originally used by Beachwatch for the daily assessment of bacterial water quality are now of limited value. General linear models (with multiple effects) were constructed to identify the secondary indicators which best describe the levels of bacterial contamination due to stormwater pollution. Data from four representative Sydney beaches were modelled separately due to the site-specific effects of pollution sources and hydrodynamics. Factors which were built into the models and were known on the day of prediction were rain amount, presence of stormwater material, the maximum daily stormwater rating, drain or lagoon flows at the time of sampling, and the drain/lagoon bacterial concentration two days before sampling.

Results indicate that there is a general linear relationship between visual indicators and bacterial density; however, the reliability of the prediction on a daily basis is poor and rainfall alone as a predictor is equivalent to or better than visual assessment alone. The primary source of information for the preparation of Beachwatch daily pollution bulletins now is rainfall reports from a telemetered network of rain gauges along the coast. This is a more cost-effective approach and has allowed the reallocation of resources to address more directly the problem of stormwater pollution at higher risk beaches.     

Helium isotope geochemistry of mid-ocean ridge basalts from the South Atlantic

We report new helium isotope results for 49 basalt glass samples from the Mid-Atlantic Ridge between 1°N and 47°S.³He/⁴He in South Atlantic mid-ocean ridge basalts (MORB) varies between 6.5 and 9.0 R_A (R_A is the atmospheric ratio of1.39 × 10⁻⁶), encompassing the range of previously reported values for MORB erupted away from high³He/⁴He hotspots such as Iceland. He, Sr and Pb isotopes show systematic relationships along the ridge axis. The ridge axis is segmented with respect to geochemical variations, and local spike-like anomalies in³He/⁴He, Pb and Sr isotopes, and trace element ratios such as(La/Sm)_N are prevalent at the latitudes of the islands of St. Helena, Tristan da Cunha and Gough to the east of the ridge. The isotope systematics are consistent with injection beneath the ridge of mantle “blobs” enriched in radiogenic He, Pb and Sr, derived from off-axis hotspot sources. The variability in³He/⁴He along the ridge can be used to refine the hotspot source-migrating-ridge sink model.

MORB from the 2–7°S segment are systematically the least radiogenic samples found along the mid-ocean ridge system to date. Here the depleted mantle source is characterized by⁸⁷Sr/⁸⁶Sr of 0.7022, Pb isotopes close to the geochron and with²⁰⁶Pb/²⁰⁴Pb of 17.7, and³He/⁴He of 8.6–8.9 R_A. The “background contamination” of the subridge mantle, by radiogenic helium derived from off-ridge hotspots, displays a maximum between 20 and 24°S. The HePb and HeSr isotope relations along the ridge indicate that the³He/⁴He ratios are lower for the hotspot sources of St. Helena, Tristan da Cunha and Gough than for the MORB source, consistent with direct measurements of³He/⁴He ratios in the island lavas. Details of the HeSrPb isotope systematics between 12 and 22°S are consistent with early, widespread dispersion of the St. Helena plume into the asthenosphere, probably during flattening of the plume head beneath the thick lithosphere prior to continental breakup. The geographical variation in theHe/Pbratio deduced from the isotope systematics suggests only minor degassing of the plume during this stage. Subsequently, it appears that the plume component reaching the mid-Atlantic ridge was partially outgassed of He during off-ridge hotspot volcanism and related melting activity.

Overall, the similar behavior of He and Pb isotopes along the ridge indicates that the respective mantle sources have evolved under conditions which produced related He and Pb isotope variations.     

Geochemistry of waters in the Valley of Ten Thousand Smokes region, Alaska

Meteoric waters from cold springs and streams outside of the 1912 eruptive deposits filling the Valley of Ten Thousand Smokes (VTTS) and in the upper parts of the two major rivers draining the 1912 deposits have similar chemical trends. Thermal springs issue in the mid-valley area along a 300-m lateral section of ash-flow tuff, and range in temperature from 21 to 29.8°C in early summer and from 15 to 17°C in mid-summer. Concentrations of major and minor chemical constituents in the thermal waters are nearly identical regardless of temperature. Waters in the downvalley parts of the rivers draining the 1912 deposits are mainly mixtures of cold meteoric waters and thermal waters of which the mid-valley thermal spring waters are representative. The weathering reactions of cold waters with the 1912 deposits appear to have stabilized and add only subordinate amounts of chemical constituents to the rivers relative to those contributed by the thermal waters. Isotopic data indicate that the mid-valley thermal spring waters are meteoric, but data is inconclusive regarding the heat source. The thermal waters could be either from a shallow part of a hydrothermal system beneath the 1912 vent region or from an incompletely cooled, welded tuff lens deep in the 1912 ash-flow sheet of the upper River Lethe area.Bicarbonate-sulfate waters resulting from interaction of near-surface waters and the cooling 1953–1968 southwest Trident plug issue from thermal springs south of Katmai Pass and near Mageik Creek, although the Mageik Creek spring waters are from a well-established, more deeply circulating hydrothermal system. Katmai caldera lake waters are a result of acid gases from vigorous drowned fumaroles dissolving in lake waters composed of snowmelt and precipitation.     

Data quality in applied geochemistry: the requirements, and how to achieve them

To obtain a reliable interpretation of geochemical data, we need to ensure that the measurements are of sufficiently high quality. To do this we need to be able to specify quantitatively the data quality that we need, and then to be certain that we are getting it. The latter can be achieved only if properly validated analytical methods are employed, and well-designed quality control schemes are used as a matter of routine.