Selective concentration of cesium in analcime during hydrothermal alteration,Yellowstone National Park,Wyoming

Chemical and mineralogical studies of fresh and hydrothermally altered rhyolitic material in Upper and Lower Geyser Basins, Yellowstone National Park, show that all the altered rocks are enriched in Cs and that Cs is selectively concentrated in analcime. The Cs content of unaltered rhyolite lava flows, including those from which the altered sediments are derived, ranges from 2.5 to 7.6 ppm. The Cs content of analcime-bearing altered sedimentary rocks is as high as 3000 ppm, and in clinoptilolite-bearing altered sedimentary rocks Cs content is as high as 180 ppm. Altered rhyolite lava flows which were initially vitrophyres, now contain up to 250 ppm Cs, and those which were crystallized prior to hydrothermal alteration contain up to 14 ppm. Mineral concentrates of analcime contain as much as 4700 ppm Cs. The Cs must have been incorporated into the analcime structure during crystallization, rather than by later cation substitution, because analcime does not readily exchange Cs. The $\text{Cs}\text{Cl}$ of the fluids circulating through the hydrothermal system varies, suggesting that Cs is not always a conservative ion and that Cs is lost from upflowing thermal waters due to water-rock interaction resulting in crystallization of Cs-bearing analcime.The source of Cs for Cs enrichment of the altered rocks is from leaching of rhyolitic rocks underlying the geyser basins, and from the top of the silicic magma chamber that underlies the area.Analcime is an important natural Cs sink, and the high Cs concentrations reported here may prove to be an important indicator of the environment of analcime crystallization.     

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.     

The value of a physically based model versus an empirical approach in the prediction of ephemeral gully erosion for loess-derived soils

A data set on soil losses and controlling factors for 58 ephemeral gullies has been collected in the Belgian loess belt from March 1997 to March 1999. Of the observed ephemeral gullies, 32 developed at the end of winter or in early spring (winter gullies) and 26 ephemeral gullies developed during summer (summer gullies). The assessed data have been used to test the physically based Ephemeral Gully Erosion Model (EGEM) and to compare its performance with the value of simple topographical and morphological indices in the prediction of ephemeral gully erosion.Analysis shows that EGEM is not capable of predicting ephemeral gully cross-sections well. Although conditions for input parameter assessment were ideal, some parameters such as channel erodibility, critical flow shear stress and local rainfall depth showed great uncertainty. Rather than revealing EGEM's inability of predicting ephemeral gully erosion, this analysis stresses the problematic nature of physically based models, since they often require input parameters that are not available or can hardly be obtained.With respect to the value of simple topographical and morphological indices in predicting ephemeral gully erosion, this study shows that for winter gullies and summer gullies, respectively, over 80% and about 75% of the variation in ephemeral gully volume can be explained when ephemeral gully length is known. Moreover, when previously collected data for ephemeral gullies in two Mediterranean study areas and the data for summer gullies formed in the Belgian loess belt are pooled, it appears that one single length (L)–volume (V) relation exists (V=0.048 L^1.29; R²=0.91). These findings imply that predicting ephemeral gully length is a valuable alternative for the prediction of ephemeral gully volume. A simple procedure to predict ephemeral gully length based on topographical thresholds is presented here. Secondly, the empirical length–volume relation can also be used to convert ephemeral gully length data extracted from aerial photos into ephemeral gully volumes.     

Latent heat effects of the major mantle phase transitions on low-angle subduction

Very low to zero shallow dip angles are observed at several moderately young subduction zones with an active trenchward moving overriding plate. We have investigated the effects of latent heat for this situation, where mantle material is pushed through the major mantle phase transitions during shallow low-angle subduction below the overriding plate. The significance of the buoyancy forces, arising from the latent heat effects, on the dynamics of the shallowly subducting slab is examined by numerical modeling. When a 32-Ma-old slab is overridden with 2.5 cm/yr by a continent, flat subduction occurs with a 4–5 cm/yr convergence rate. When latent heat is included in the model, forced downwellings cause a thermal anomaly and consequently thermal and phase buoyancy forces. Under these circumstances, the flat slab segment subducts horizontally about 350 km further and for about 11 Ma longer than in the case without latent heat, before it breaks through the 400-km phase transition. The style of subduction strongly depends on the mantle rheology: increasing the mantle viscosity by one order of magnitude can change the style of subduction from steep to shallow. Similarly, an overriding velocity of less than 1 cm/yr leads to steep subduction, which gradually changes to flat subduction when increasing the overriding velocity. However, these model parameters do not change the aforementioned effect of the latent heat, provided that low-angle subduction occurs. In all models latent heat resulted in a substantial increase of the flat slab length by 300–400 km. Varying the olivine–spinel transition Clapeyron slope γ from 1 to 6 MPa/K reveals a roughly linear relation between γ and the horizontal length of the slab. Based on these results, we conclude that buoyancy forces due to latent heat of phase transitions play an important role in low-angle subduction below an overriding plate.     

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.     

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.