Gold ore-forming fluids of the Tanami region, Northern Australia

Fluid inclusion studies have been carried out on major gold deposits and prospects in the Tanami region to determine the compositions of the associated fluids and the processes responsible for gold mineralization. Pre-ore, milky quartz veins contain only two-phase aqueous inclusions with salinities ≤19 wt% NaCl eq. and homogenization temperatures that range from 110 to 410°C. In contrast, the ore-bearing veins typically contain low to moderate salinity (<14 wt% NaCl eq.), H₂O + CO₂ ± CH₄ ± N₂-bearing fluids. The CO₂-bearing inclusions coexist with two-phase aqueous inclusions that exhibit a wider range of salinities (≤21 wt% NaCl eq.). Post-ore quartz and carbonate veins contain mainly two-phase aqueous inclusions, with a last generation of aqueous inclusions being very CaCl₂-rich. Salinities range from 7 to 33 wt% NaCl eq. and homogenization temperatures vary from 62 to 312°C. Gold deposits in the Tanami region are hosted by carbonaceous or iron-rich sedimentary rocks and/or mafic rocks. They formed over a range of depths at temperatures from 200 to 430°C. The Groundrush deposit formed at the greatest temperatures and depths (260–430°C and ≤11 km), whereas deposits in the Tanami goldfield formed at the lowest temperatures (≥200°C) and at the shallowest depths (1.5–5.6 km). There is also evidence in the Tanami goldfield for late-stage isothermal mixing with higher salinity (≤21 wt% NaCl eq.) fluids at temperatures between 100 and 200°C. Other deposits (e.g., The Granites, Callie, and Coyote) formed at intermediate depths and at temperatures ranging from 240 to 360°C. All ore fluids contained CO₂ ± N₂ ± CH₄, with the more deeply formed deposits being enriched in CH₄ and higher level deposits being enriched in CO₂. Fluids from deposits hosted mainly by sedimentary rocks generally contained appreciable quantities of N₂. The one exception is the Tanami goldfield, where the quartz veins were dominated by aqueous inclusions with rare CO₂-bearing inclusions. Calculated δ ¹⁸O values for the ore fluids range from 3.8 to 8.5‰ and the corresponding δD values range from −89 to −37‰. Measured δ ¹³C values from CO₂ extracted from fluid inclusions ranged from −5.1 to −8.4‰. These data indicate a magmatic or mixed magmatic/metamorphic source for the ore fluids in the Tanami region. Interpretation of the fluid inclusion, alteration, and structural data suggests that mineralization may have occurred via a number of processes. Gold occurs in veins associated with brittle fracturing and other dilational structures, but in the larger deposits, there is also an association with iron-rich rocks or carbonaceous sediments, suggesting that both structural and chemical controls are important. The major mineralization process appears to be boiling/effervescence of a gas-rich fluid, which leads to partitioning of H₂S into the vapor phase resulting in gold precipitation. However, some deposits also show evidence of desulfidation by fluid–rock interaction and/or reduction of the ore-fluid by fluid mixing. These latter processes are generally more prevalent in the higher crustal-level deposits.     

Accelerated initial stiffness schemes for elastoplasticity

Iterative methods for the solution of non‐linear finite element equations are generally based on variants of the Newton–Raphson method. When they are stable, full Newton–Raphson schemes usually converge rapidly but may be expensive for some types of problems (for example, when the tangent stiffness matrix is unsymmetric). Initial stiffness schemes, on the other hand, are extremely robust but may require large numbers of iterations for cases where the plastic zone is extensive. In most geomechanics applications it is generally preferable to use a tangent stiffness scheme, but there are situations in which initial stiffness schemes are very useful. These situations include problems where a nonassociated flow rule is used or where the zone of plastic yielding is highly localized. This paper surveys the performance of several single‐parameter techniques for accelerating the convergence of the initial stiffness scheme. Some simple but effective modifications to these procedures are also proposed. In particular, a modified version of Thomas' acceleration scheme is developed which has a good rate of convergence. Previously published results on the performance of various acceleration algorithms for initial stiffness iteration are rare and have been restricted to relatively simple yield criteria and simple problems. In this study, detailed numerical results are presented for the expansion of a thick cylinder, the collapse of a rigid strip footing, and the failure of a vertical cut. These analyses use the Mohr–Coulomb and Tresca yield criteria which are popular in soil mechanics. Copyright © 2000 John Wiley & Sons, Ltd.     

Influence of wave and sediment dynamics on cordgrass (<Emphasis Type="Italic">Spartina anglica</Emphasis>) growth and sediment accumulation on an exposed intertidal flat

Isolated patches ofSpartina anglica (cordgrass) at two sites on a wave-exposed mid-intertidal flat of the 340 km² Manukau Harbor (Auckland, New Zealand) have developed very differently since being planted in the mid-1970s. Although the two sites are only 0.5 km apart and at the same intertidal elevation,Spartina patches at the easternmost site (site 1) have as much as an order of magnitude higher biomass and accumulated sediment volume thanSpartina patches at site 2. A field experiment was conducted to characterize waves and associated sediment dynamics at each site, which might explain whySpartina patches at the two sites have developed so differently over the past 25 yr or so. Suspended sediments were measured and wave characteristics were inferred from subsurface pressure data measured for 5 wk at bothSpartina sites and at an intermediate location. Bed-orbital speeds and frictional wave-energy dissipation were consistently lower at the easternmost site with the largerSpartina patches. The west-to-east reduction in wave energy is due to the spatial arrangement of theSpartina sites relative to the predominant wind fetches. The wave-energy gradient is maintained by tidal-cycle variations in fetch and bed friction and results in a west-to-east reduction in sand suspension. Silt, which is largely resuspended under southwest winds, is redeposited in the low wave-energy conditions in and around the larger site 1Spartina patches. Shell accumulation bySpartina patches at site 1 occurs infrequently, during southwest winds >10 m s⁻¹ and water depths >0.7 m, when waves are least attenuated by bed friction. Large between-site differences in the growth of and sediment accumulation by theSpartina patch are consistent with the observed wave-energy gradient. The resulting spatial patterns of silt, sand, and shell resuspension and deposition directly influence the rate of sediment accumulation bySpartina patches and the composition of accumulated sediment on this wave-exposed intertidal flat.     

Raman microspectroscopy: A non-destructive tool for routine calibration of apatite crystallographic structure for fission-track analyses

Compositional control on the annealing kinetics of fission-tracks (FT) in apatite requires routine measurement of sample grain composition. However, for practical reasons the bulk composition of analysed grains is not routinely measured and instead grain chlorine content or etch-pit dimensions are used to characterise a samples annealing behaviour. A more desirable approach is to measure crystallographic parameters (i.e. unit cell dimension) of a grain as these represent the summed effect of all substitutions and crystal defects. We show how Raman microspectrometry can be used as a routine non-destructive tool to obtain rapid measurement of the crystallographic structure of apatite grains etched for FT analysis. Variations of unit cell parameter a are found to correspond to a systematic variation of Raman shift in the range of 452–440 cm^− 1 for measurements made on c-parallel sections of apatite where the direction of the polarized incident beam is parallel to the c axis.     

Using streamflow characteristics to explore permafrost thawing in northern Swedish catchments

The recent and rapid warming of the Arctic leads to thawing of permafrost, which influences and changes subsurface water-flow systems in such landscapes. This study explores the utility of catchments as “sentinels of change” by considering long-term discharge data from 17 stations on unregulated rivers in northern Sweden and analyzing trends in annual minimum discharge and recession flow characteristics. For the catchments considered, the annual minimum discharge has increased significantly (based on the Mann Kendall test at a 95 % confidence level) in nine of the catchments and decreased significantly in one catchment. Considering changes in recession-flow characteristics, seven catchments showed significant trends consistent with permafrost thawing while two catchments showed significant trends in the opposite direction. These results are mechanistically consistent with generic physically based modeling studies and the geological setting, as the catchments considered span the spatial limit of permafrost extent. This study illuminates the potential for using hydrologic observations to monitor changes in catchment-scale permafrost. Further, this opens the door for research to isolate the mechanisms behind the different trends observed and to gauge their ability to reflect actual permafrost conditions at the catchment scale.     

A latest Pleistocene and Holocene glacial history and paleoclimate reconstruction at Three Sisters and Broken Top Volcanoes, Oregon, U.S.A.

At least three sets of moraines mark distinct glacial stands since the last glacial maximum (LGM) in the Three Sisters region of the Oregon Cascade Range. The oldest stand predates 8.1 ka (defined here as post-LGM), followed by a second between ∼ 2 and 8 ka (Neoglacial) and a third from the Little Ice Age (LIA) advance of the last 300 years. The post-LGM equilibrium line altitudes were 260 ± 100 m lower than that of modern glaciers, requiring 23 ± 9% increased winter snowfall and 1.4 ± 0.5°C cooler summer temperatures than at present. The LIA advance had equilibrium line altitudes 110 ± 40 m lower than at present, implying 10 ± 4% greater winter snowfall and 0.6 ± 0.2°C cooler summer temperatures.     

Elemental and isotopic fractionation of Type B calcium-, aluminum-rich inclusions: experiments, theoretical considerations, and constraints on their thermal evolution

Experiments exposing Type B calcium-, aluminum-rich inclusion (CAI)-like melts at high temperatures to high vacuum or reducing hydrogen-rich gas mixtures were used to determine the rates and consequences of elemental and isotopic fractionation by evaporation. Silicon and magnesium were found to evaporate much faster than calcium and aluminum, and the resulting residual liquid trajectories in composition space are reproduced via a thermodynamic model for the saturation vapor pressure of the evaporating species. Isotopic fractionations associated with evaporation were measured for magnesium. The resulting relationship between fraction of magnesium lost and enrichment of the residue in the heavy isotopes of magnesium follows a Rayleigh fractionation curve with a fractionation factor that is close to, but not exactly, the theoretically expected value. The rate of evaporation is found to be a strong function of temperature, oxygen fugacity, and melt composition, which can be understood and modeled in terms of the dependence of the saturation vapor pressures on these variables. The relationship between evaporation rate, which we measure, and calculated saturation vapor involves empirical evaporation coefficients that we find to be significantly less than one (∼0.1). Analytical and numerical models are used to characterize how diffusion in both the melt and in the surrounding gas affects evaporation rates and the degree of chemical and isotopic fractionation. The experimental data and theoretical considerations are combined to give a parameterization of the rates and consequences of evaporation of Type B CAI-like liquids, which is then used to translate the measured isotopic fractionation of Type B CAIs into constraints on their thermal history. Cooling rates of the order of 10°C per hour are indicated.     

Carbonatization and mobility of Ti,Y, and Zr in Ascot Formation metabasalts,SE Quebec

The middle Ordovician Ascot Formation of southeastern Quebec consists of greenschist facies metamorphosed silicic to mafic pyroclastic rocks and lava flows and associated metasediments. Chemical analyses of lavas reveal a preponderance of metarhyolites and metabasalts, together with some porphyritic rocks with intermediate SiO₂ contents. The metabasalts exhibit wide ranges in concentrations of TiO₂ (0.25–2.0 wt.°), Y (9–46 ppm), and Zr (5–135 ppm). The extent of the ranges, and unusual interelement ratios, suggest that the concentrations of these normally immobile elements have been affected by secondary processes.There is a strong correlation between trace-element concentrations and the degree of carbonatization of the metabasalts. Low carbonate rocks are severely depleted in Ti, Y, and Zr whereas high carbonate rocks are depleted in Y and Zr and enriched in Ti. The differing movement of Ti can be explained in terms of variable chemical potential due to the various carbonatization reactions affecting titaniferous phases. Overall mobility of these generally immobile/rd elements is attributed to high CO₂ levels in the fluid phase during metamorphism.Extrapolation of the two alteration trends to a common origin enables one to infer primary concentrations of the trace-elements. Primary inter-element ratios arrived at in this way are compatible with an island-arc origin for the Ascot Formation although TiO₂ concentrations are a little high (1.5 wt/%).     

Lithological and spatial controls on the distribution of quartz veins in andesite- and rhyolite-hosted epithermal Au–Ag deposits of the Hauraki Goldfield, New Zealand

Vein distributions in line samples from four epithermal Au–Ag deposits of the Hauraki Goldfield were logged and quantified by vein spacing, vein density, vein thickness and percentage of vein extension. One deposit is hosted in andesite lavas (Martha Hill), one in andesite lavas and dacite porphyry, dacitic tuffs and pyroclastic breccias (Golden Cross), and two in rhyolite lavas and rhyolitic tuffs with minor andesite lavas or andesite dikes (Ohui and Wharekirauponga). The vein systems in these deposits form fault-controlled arrays of extensional veins. Vein spacing distributions are non-fractal over two to three orders of magnitude (1 mm to 5 m), and therefore fractal dimension statistics are not applicable. The coefficient of variation (C_v) of vein spacing was used as a measure of the degree of vein clustering. Rock type has a marked influence on vein spacing distributions, with veining in rhyolite lava having lower average thickness and percentage extension, but a generally higher degree of vein clustering compared with veining in andesite lava in the same deposit. Vein spacing distributions in well-jointed lithologies, mainly andesite lava, have C_v values (0.8–1.2) that are indicative of anticlustered to weakly clustered patterns, particularly in the vein stockwork of the upper part of the Golden Cross deposit. These C_v values are consistent with field observations that joints are a major control on vein spacing. In the poorly jointed dacitic and rhyolitic rocks, the veins are weakly to strongly clustered as shown by higher C_v values (1.2–2.4), and are commonly associated with normal faults. Overall, andesite lava and dacite porphyry and pyroclastics host thicker and more persistent veins than rhyolite lava and tuff. These larger veins contain significant volumes of high-grade gold mineralisation. The higher chemical reactivity to hydrothermal fluids of andesite and dacite compared with rhyolite may have aided propagation and thickening of the veins in andesite-hosted deposits. Within an individual epithermal deposit, location close to thick veins, representing major fluid conduits, commonly overrides the effect of different lithologies. Sites that are deeper and located within or adjacent to major vein structures have higher average vein thickness, percentage extension and degree of vein clustering. Systematic collection and analysis of vein spacing, thickness and density data can be used to define trends that are useful in the exploration of gold-bearing epithermal vein deposits. Received: 25 August 1998 / Accepted: 23 December 1999