Cautions and Suggestions for Geochemical Sampling in Fractured Rock

Collecting water samples for geochemical analyses in open bedrock boreholes or in discrete intervals of boreholes intersected by multiple fractures is likely to yield ambiguous results for ground water chemistry because of the variability in the transmissivity, storativity, and hydraulic head of fractures intersecting the borehole. Interpreting chemical analyses of water samples collected in bedrock boreholes requires an understanding of the hydraulic conditions in the borehole under the ambient flow regime in the aquifer as well as during sampling. Pumping in open boreholes, regardless of the pumping rate and the location of the pump intake, first draws water from the borehole and then from fractures intersecting the borehole. The time at which the volumetric rate of water entering the borehole from fractures is approximately equal to the pumping rate can be identified by monitoring the logarithm of drawdown in the borehole as a function of the logarithm of time. Mixing of water entering the borehole from fractures with water in the borehole must be considered in estimating the time at which the pump discharge is representative of aquifer water. In boreholes intersected by multiple fractures, after the contribution from the borehole volume has diminished, the contribution of fractures to the pump discharge will be weighted according to their transmissivity, regardless of the location of the pump intake. This results in a flux-averaged concentration in the pump discharge that is biased by the chemical signature of those fractures with the highest transmissivity. Under conditions where the hydraulic head of fractures varies over the length of the borehole, open boreholes will be subject to ambient flow in the water column in the borehole. In some instances, the magnitude of the ambient flow may be similar to the designated pumping rate for collecting water samples for geochemical analyses. Under such conditions, the contributions to the pump discharge from individual fractures will be a function not only of the transmissivity of the fractures, but also of the distribution of hydraulic head in fractures intersecting the borehole. To reduce or eliminate the deleterious effects of conducting geochemical sampling in open boreholes, a straddle-packer apparatus that isolates a single fracture or a series of closely spaced fractures is recommended. It is also recommended that open boreholes be permanently outfitted with borehole packers or borehole liners in instances where maintaining the hydraulic and chemical stratification in the aquifer is of importance. In a field example, a comparison of results from sampling in an open borehole and in discrete intervals of the same borehole showed dramatic differences in the concentrations of chemical constituents in the water samples, even though chemical field parameters stabilized prior to both open borehole and discrete interval sampling.     

新形势下林区的就业机会与渠道

我国实施“天然林保护工程”以来,林区人民积极转产分流,寻找新的就业门路和生计,为走出困境、改善生活做出了不懈的努力。为了进一步开阔视野,拓展思路,借鉴他国的经验,我们选刊了“国际森林趋势组织”的几篇关于寻找新的林区就业渠道和森林环境服务市场开发方面的研究文章,以期对我国的林业改革与建设有所启示。(注:这些文章已经“森林趋势组织”的执行主管Mivhael Jenkins先生同意译成中文在本刊刊登,本刊作了部分删节和改动。原文及参考文献请参阅森林趋势组织的网页:http://www. forest-trends. org)     

The interpretation of reaction textures in Fe‐rich metapelitic granulites of the Musgrave Block,central Australia: constraints from mineral equilibria calculations in the system K2O–FeO–MgO–Al2O3–SiO2–H2O–TiO2–Fe2O3

Fe‐rich metapelitic granulites of the Musgrave Block, central Australia, contain several symplectic and coronal reaction textures that post‐date a peak S2 metamorphic assemblage involving garnet, sillimanite, spinel, ilmenite, K‐feldspar and quartz. The earliest reaction textures involve spinel‐ and quartz‐bearing symplectites that enclose garnet and to a lesser extent sillimanite. The symplectic spinel and quartz are in places separated by later garnet and/or sillimanite coronas. The metamorphic effects of a later, D3, event are restricted to zones of moderate to high strain where a metamorphic assemblage of garnet, sillimanite, K‐feldspar, magnetite, ilmenite, quartz and biotite is preserved. Quantitative mineral equilibria calculations in the system K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–Fe₂O₃ (KFMASHTO) using Thermocalc 3.0 and the accompanying internally consistent dataset provide important constraints on the influence of TiO₂ and Fe₂O₃ on biotite‐bearing and spinel‐bearing equilibria, respectively. Biotite‐bearing equilibria are shifted to higher temperatures and spinel‐bearing equilibria to higher pressures and lower temperatures in comparison to the equivalent equilibria in K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O (KFMASH). The sequence of reaction textures involving spinel is consistent with a D2 P–T path that involved a small amount of decompression followed predominantly by cooling within a single mineral assemblage stability field. Thus, the reaction textures reflect changes in modal proportions within an equilibrium assemblage rather than the crossing of a univariant reaction. The D3 metamorphic assemblage is consistent with lower temperatures than those inferred for D2.     

A mineral equilibria study of the hydrothermal alteration in mafic greenschist facies rocks at Kalgoorlie, Western Australia

The influx of a H₂O–CO₂‐dominated fluid into actinolite‐bearing metabasic rocks during greenschist facies metamorphism in the Kalgoorlie area of Western Australia resulted in a zoned alteration halo around inferred fluid conduits that contain gold mineralisation. The alteration halo is divided into two outer zones, the chlorite zone and the carbonate zone, and an inner pyrite zone adjacent to the inferred fluid conduits. Reaction between the fluid and the protolith resulted in the breakdown of actinolite and the development of chlorite, dolomite, calcite and siderite. In addition, rocks in the pyrite zone developed muscovite‐bearing assemblages as a consequence of the introduction of potassium by the fluid. Mineral equilibria calculations undertaken using the computer software thermocalc in the model system Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–CO₂ show that mineral assemblages in the outer zones of the alteration halo are consistent with equilibrium of the protoliths with a fluid of composition X_CO2 = CO₂/(CO₂ + H₂O) = 0.1–0.25 for temperatures of 315–320 °C. The inner zone of the alteration halo reflect equilibrium with a fluid of composition X_CO2≈ 0.25. Fluid‐rock buffering calculations show that the alteration halo is consistent with interaction with a single fluid composition and that the zoned structure of the halo reflects the volume of this fluid with which the rocks reacted. This fluid is likely to have also been the one responsible for the gold mineralisation at Kalgoorlie.     

The influence of thermochemical convection on the fixity of mantle plumes

A general feature of both isochemical and thermochemical studies of mantle convection is that horizontal plume velocities tend to be smaller than typical convective velocities, however, it is not clear which system leads to a greater fixity of mantle plumes. We perform two- and three-dimensional numerical calculations and compare both thermochemical and isochemical cases with similar convective vigor to determine whether presence of a dense component in the mantle can lead to smaller ratios of horizontal plume velocity to surface velocity. We investigate different viscosity and density contrasts between chemical components in the thermochemical calculations, and we perform isochemical calculations with both free-slip and no-slip bottom boundary conditions. We then compare both visually and quantitatively the results of the thermochemical and isochemical calculations to determine which leads to greater plume fixity. We find that horizontal plume velocities for thermochemical calculations are similar to those from isochemical calculations with no-slip bottom boundary conditions. In addition, we find that plumes tend to be more fixed for isochemical cases with free-slip bottom boundary conditions for two-dimensional calculations, however, in three dimensions, we find that plume fixity is similar to that observed in thermochemical calculations.     

Dust production and the release of iron oxides resulting from the aeolian abrasion of natural dune sands

Global dust trajectories indicate that signi?cant quantities of aeolian‐transported iron oxides originate in contemporary dryland areas. One potential source is the iron‐rich clay coatings that characterize many sand‐sized particles in desert dune?elds. This paper uses laboratory experiments to determine the rate at which these coatings can be removed from dune sands by aeolian abrasion. The coatings impart a red colour to the grains to which previous researchers have assigned variable geomorphological signi?cance. The quantities of iron removed during a 120 hour abrasion experiment are small (99 mg kg^?1) and dif?cult to detect by eye; however, high resolution spectroscopy clearly indicates that ferric oxides are released during abrasion and the re?ectance of the particles alters. One of the products of aeolian abrasion is ?ne particles (<10 µm diameter) with the potential for long distance transport. Copyright © 2004 John Wiley & Sons, Ltd.     

Coherence methods in mapping AVO anomalies and predicting P-wave and S-wave impedances

Filters for migrated offset substacks are designed by partial coherence analysis to predict ‘normal’ amplitude variation with offset (AVO) in an anomaly free area. The same prediction filters generate localized prediction errors when applied in an AVO‐anomalous interval. These prediction errors are quantitatively related to the AVO gradient anomalies in a background that is related to the minimum AVO anomaly detectable from the data. The prediction‐error section is thus used to define a reliability threshold for the identification of AVO anomalies. Coherence analysis also enables quality control of AVO analysis and inversion. For example, predictions that are non‐localized and/or do not show structural conformity may indicate spatial variations in amplitude–offset scaling, seismic wavelet or signal‐to‐noise (S/N) ratio content. Scaling and waveform variations can be identified from inspection of the prediction filters and their frequency responses. S/N ratios can be estimated via multiple coherence analysis. AVO inversion of seismic data is unstable if not constrained. However, the use of a constraint on the estimated parameters has the undesirable effect of introducing biases into the inverted results: an additional bias‐correction step is then needed to retrieve unbiased results. An alternative form of AVO inversion that avoids additional corrections is proposed. This inversion is also fast as it inverts only AVO anomalies. A spectral coherence matching technique is employed to transform a zero‐offset extrapolation or near‐offset substack into P‐wave impedance. The same technique is applied to the prediction‐error section obtained by means of partial coherence, in order to estimate S‐wave velocity to P‐wave velocity (V_S/V_P) ratios. Both techniques assume that accurate well ties, reliable density measurements and P‐wave and S‐wave velocity logs are available, and that impedance contrasts are not too strong. A full Zoeppritz inversion is required when impedance contrasts that are too high are encountered. An added assumption is made for the inversion to the V_S/V_P ratio, i.e. the Gassmann fluid‐substitution theory is valid within the reservoir area. One synthetic example and one real North Sea in‐line survey illustrate the application of the two coherence methods.     

Variation in nitrate and calcium as indicators of recharge pathways in Nolte Spring, PA

This study documents variations in calcium and nitrate concentrations that suggest changes in recharge pathways in a karst spring. The nitrate concentrations increased at the end of the growing season, showing the importance of the soil zone in the recharge pathway. The increase occurred over just a few days, which may be indicative of a change in contribution of baseflow in different seasons from deep to shallow groundwater. The calcium concentrations decreased several days after storm events. A change in the carbonate equilibrium is hypothesized because chloride was not diluted during these events. The decrease in calcium could be due to outgassing and calcite precipitation in the recharge area when older, higher ionic strength matrix water mixes with stormwater in open conduits. The use of geochemical indicators to better understand recharge pathways benefited from long-term monitoring and periods of daily sampling.