A review of the metallogeny and tectonics of the Lachlan Orogen

Future mineral exploration within eastern Australia will be enhanced by resolving the tectonic evolution of the Lachlan Orogen to establish the spatial and temporal terrane distribution of the various mineral deposits. The Lachlan Orogen, from north-eastern Tasmania through to central and eastern New South Wales, is host to a number of major mineral deposit styles—including orogenic gold (e.g. Stawell, Ballarat, Bendigo), volcanic-hosted massive sulphide (e.g. Woodlawn, Currawong), sediment-hosted Cu–Au (e.g. Cobar Basin deposits), porphyry Au–Cu (e.g. Cadia, Parkes, Cowal) and granite-related Sn (e.g. Ardlethan, Beechworth). Each of these mineral deposit styles is a sensitive and diagnostic indicator of the prevalent tectonic environment during their formation. In this review, we briefly summarise the deposit- to large-scale factors that define the diverse metallogenic evolution of the Lachlan Orogen. This overview is intended to “set the scene” for subsequent specialist papers published in this thematic issue on the metallogeny and tectonics of the Lachlan Orogen in south-east Australia.     

Unravelling the evolution of an Alpine to post-glacially active fault in the Swiss Alps

The Gemmi fault is a prominent NW–SE striking lineament that crosses the Gemmi Pass in the central Swiss Alps. A multidisciplinary investigation of this structure that included geological mapping, joint profiling, cathodoluminescence and scanning electron microscopy, stable isotope measurements, luminescence- and U-TH-dating, 3D ground penetrating radar (GPR) surveying and trenching reveals a history of fault movements from the Miocene to the Holocene. The main fault zone comprises a 0.5–3 m thick calcite cataclasite formed during several cycles of veining and brittle deformation. Displaced Cretaceous rock layers show an apparent dextral slip of 10 m along the fault.A detailed study of a small sediment-filled depression that crosses the fault provides evidence for a post-glacial reactivation of the fault. A trench excavated across the fault exposed a Late-Glacial-age loess layer and late Holocene colluvial-like slope-wash deposits that showed evidence for fault displacement of a few centimeters, indicating a recent strike-slip reactivation of the fault. Focal mechanisms of recent instrumentally recorded earthquakes are consistent with our findings that show that the fault at the Gemmi Pass, together with other parallel faults in this area, may be reactivated in today's stress field. Taking together all the observations of its ancient and recent activity, the Gemmi fault can be viewed as a window through geological space and time.     

Isotopic and noble gas study of Chalk groundwater in the London Basin, England

The chemical and isotopic composition of groundwater from 52 sites in the London (U.K.) area was determined as part of a project aimed at assessing the spatial variation in the age of Chalk groundwater, and in determining the relationship between fracture and matrix groundwater in this dual porosity system.Systematic changes in groundwater chemistry take place in the downgradient direction in response to several chemical processes. These processes include early concentration by evaporation and congruent dissolution of calcite followed by widespread incongruent dissolution and ion exchange in addition to local oxidation-reduction reactions, gypsum dissolution and saline intrusion. As a result of the above processes, Chalk groundwater follows an evolutionary path from Ca bicarbonate type to Na bicarbonate type.The age of Chalk groundwater was modelled using¹⁴C, δ¹³C,³H, δ²H and δ¹⁸0. There is a general increase in the groundwater age in a downgradient direction with the oldest water found in N central areas of the basin. Groundwater in the unconfined zones and in areas S of the Greenwich fault is almost entirely of unevolved, modem composition. Carbon-14 modelling suggests that Chalk groundwater in the S basin is generally less than 10000 a old while that in the north is generally between 10000 and 25000 a old. The presence of³H in concentrations of up to 7 TU in groundwater which yields ages of several 1000 a, however, indicates that mechanisms exist for the rapid introduction of recent groundwater to the confined aquifer. Results of palaeorecharge temperature determinations using δ²H, δ¹⁸0 and noble gas analytical results suggest that significant Devensian recharge did indeed occur in the aquifer.A model of the development of the Chalk recognizes that it is a classic dual porosity aquifer in which groundwater flow occurs predominantly in the fracture system. The upper 50 m of the aquifer was flushed with fresh water during the 2–3 × 10⁶ a of the Quaternary and therefore meteoric water largely replaced the Tertiary and Cretaceous marine water that previously saturated the system. Most processes which control the chemistry of the groundwater occur in the matrix where the surface area is exceptionally high. Although fracture flow dominates the flow regime, diffusion from the matrix into the fracture porosity controls the chemistry of Chalk groundwater.     

Effects of diffusional modification of garnet growth zoning on P-T path calculations

The effect of intragranular diffusion on chemical zoning in garnet and on P-T paths calculated from that zoning was evaluated using a numerical model of multicomponent diffusion in combination with simulations of garnet growth. Syn-and post-growth diffusion of Mg–Fe–Mn–Ca species in garnet was calculated for a model pelitic assemblage over a range of temperatures from 485 to 635°C. Compositions from zoned garnet, as modified by diffusion, hypothetical inclusions of plagioclase within garnet and matrix phases were used to reconstruct pressure-temperature (P-T) paths from isobaric and polybaric model histories. P-T path calculations, based on numerical simulations conducted over an input isobaric heating path that reached peak temperatures between 585 and 635°C, show that relaxation of garnet compositional gradients by diffusion can induce modest to appreciable curvature in the inferred paths. Retrieved paths also indicated somewhat smaller overall temperature changes relative to the actual temperature difference of the input path. The magnitude of these distortions is shown to depend upon the heating and cooling rate and garnet crystal size as well as the actual peak temperature condition. The effect of diffusion on path trajectories in simulations with thermal histories that also included cooling were comparable to heating-only models that reached peak temperatures approximately 15–30°C higher. Compositions of garnets with radii less than 1 mm, that reached actual peak temperatures of 605°C along temperature-time histories characteristic of regional metamorphism, experienced sufficient diffusional relaxation to introduce errors of hundreds of bars to in excess of one kilobar in path trajectories. Path distortions were significant at heating/cooling rates less than 10°C/Ma, but rapidly diminished for rates faster than this. In polybaric simulations diffusion effects were least noticeable when the actual pressure-temperature conditions changed in a clockwise sense (i.e., convex to higher P and higher T), but apprecciable modification was seen in path models that underwent counterclockwise changes in P and T. Reequilibration of garnet rim compositions occurred during cooling on all paths, and temperature maxima obtained from garnet-biotite geothermometry underestimated actual peak conditions by 40 to 70°C. Calculations suggest that P-T path trajectories calculated from garnets of at least 1 mm size, and that experienced actual thermal maxima below 585°C, are not likely to be distorted by diffusional effects during regional metamorphism. However, P-T path reconstructions based on garnet zonation with smaller grains or higher temperatures may lead to misinterpretation of crystallization history. The partitioning record of peak metamorphic temperatures may be destroyed by diffusional reequilibration of garnet rim compositions during cooling, seriously complicating the task of quantitatively estimating diffusion effects on path calculations.     

Tectonic and sedimentary basin evolution of the eastern Bangong–Nujiang zone (Tibet): a Reading cycle

The ophiolite-bearing Bangong-Nujiang zone (BNZ) traversing central Tibet from east to west separates the Qiangtang block in the north from the Lhasa block in the south. Their stratigraphic development indicates that both blocks once formed a continuous continental platform until the Late Triassic. Following Late Paleozoic-Triassic rifting, ocean crust formed between both blocks during the Late Triassic creating the Dongqiao-Naqu basin (DNB) among other basins (Yu et al. 1991). The analysis of the rift flank sequences reveals that rifting was dominated by transtension. The basin was shortened by post-Mid-Cretaceous transpression. Thus, the overall basin evolution represents a Reading cycle despite some active margin processes which gave this cycle a special imprint. Major basin parts were preserved despite transpressional shortening suggesting that the eastern BNZ represents a remnant basin. Our understanding of the DNB solves the prior problem of viewing the BNZ as a Mid-Late Jurassic collisional suture although typical collision-related deformation, thickening, mountain building, as well as related molasse formation are lacking. Our model also explains the scattered linear ophiolite distribution by local transpression of remnant oceanic basin floor without having to consider problematic long range ophiolite thrusting.     

Carbon Dioxide in igneous petrogenesis: I

A number of experimental CO₂ solubility data for silicate and aluminosilicate melts at a variety of P- T conditions are consistent with solution of CO₂ in the melt by polymer condensation reactions such as SiO _4(m ^4? +CO_2(v)+Si _n O _3n+1(m) ⁽²ⁿ⁺¹⁾ ?Si _n+1O _3n+4(m) ^(2n+4)? +CO _3(m ^)2? . For various metalsilicate systems the relative solubility of CO₂ should depend markedly on the relative Gibbs free change of reaction. Experimental solubility data for the systems Li₂O-SiO₂, Na₂O-SiO₂, K₂O-SiO₂, CaO-SiO₂, MgO-SiO₂ and other aluminosilicate melts are in complete accord with predictions based on Gibbs Free energies of model polycondesation reactions. A rigorous thermodynamic treatment of published P- T-wt.% CO₂ solubility data for a number of mineral and natural melts suggests that for the reaction CO_2(m) ? CO_2(v)

1. CO₂-melt mixing may be considered ideal (i.e., { \(a_{{\text{CO}}_{\text{2}} }^m = X_{{\text{CO}}_{\text{2}} }^m \) );
2. \(\bar V_{{\text{CO}}_{\text{2}} }^m \) , the partial molal volume of CO₂ in the melt, is approximately equal to 30 cm³ mole^?1 and independent of P and T;
3. Δ C _p ⁰ is approximately equal to zero in the T range 1,400° to 1,650 °C and
4. enthalpies and entropies of the dissolution reaction depend on the ratio of network modifiers to network builders in the melt. Analytic expressions which relate the CO₂ content of a melt to P, T, and \(f_{{\text{CO}}_{\text{2}} } \) for andesite, tholeiite and olivine melilite melts of the form

$$\ln X_{{\text{CO}}_{\text{2}} }^m = \ln f_{{\text{CO}}_{\text{2}} } - \frac{A}{T} - B - \frac{C}{T}(P - 1)$$ have been determined. Regression parameters are (A, B, C): andesite (3.419, 11.164, 0.408), tholeiite (14.040, 5.440,0.393), melilite (9.226, 7.860, 0.352). The solubility equations are believed to be accurate in the range 3<P<30 kbar and 1,100°<T<1,650 °C. A series of CO₂ isopleth diagrams for a wide range of T and P are drawn for andesitic, tholeiitic and alkalic melts.     

Uncertainties of geochemical modeling during CO<Subscript>2</Subscript> sequestration applying batch equilibrium calculations

One of the uncertainties in the field of carbon dioxide capture and storage (CCS) is caused by the parameterization of geochemical models. The application of geochemical models contributes significantly to calculate the fate of the CO₂ after its injection. The choice of the thermodynamic database used, the selection of the secondary mineral assemblage as well as the option to calculate pressure dependent equilibrium constants influence the CO₂ trapping potential and trapping mechanism. Scenario analyses were conducted applying a geochemical batch equilibrium model for a virtual CO₂ injection into a saline Keuper aquifer. The amount of CO₂ which could be trapped in the formation water and in the form of carbonates was calculated using the model code PHREEQC. Thereby, four thermodynamic datasets were used to calculate the thermodynamic equilibria. Furthermore, the equilibrium constants were re-calculated with the code SUPCRT92, which also applied a pressure correction to the equilibrium constants. Varying the thermodynamic database caused a range of 61% in the amount of trapped CO₂ calculated. Simultaneously, the assemblage of secondary minerals was varied, and the potential secondary minerals dawsonite and K-mica were included in several scenarios. The selection of the secondary mineral assemblage caused a range of 74% in the calculated amount of trapped CO₂. Correcting the equilibrium constants with respect to a pressure of 125 bars had an influence of 11% on the amount of trapped CO₂. This illustrates the need for incorporating sensitivity analyses into reaction pathway modeling.     

The effect of ionic interactions on the oxidation of metals in natural waters

The effect of ionic interactions of the major components of natural waters on the oxidation of Cu(I) and Fe(II) has been examined. The various ion pairs of these metals have been shown to have different rates of oxidation. For Fe(II), the chloride and sulfate ion pairs are not easily oxidized. The measured decrease in the rate constant at a fixed pH in chloride and sulfate solutions agrees very well with the values predicted. The effect of pH (6 to 8) on the oxidation of Fe(II) in water and seawater have been shown to follow the rate equation

$\text{-d in [Fe(II)]/dt = k}{}_1\text{β}{}_1\text{α}{}_{\mathrm{Fe}}\text{/[H}{}^{\mathrm{+}}\text{] + k}{}_2\text{β}{}_2\text{α}{}_{\mathrm{Fe}}\text{/[H}{}^{\mathrm{+}}\text{]}{}^2$

where k₁ and k₂ are the pseudo first order rate constants, β₁ and β₂ are the hydrolysis constants for Fe(OH)⁺ and Fe(OH)⁰. The value of α_FE is the fraction of free Fe²⁺. The value of k₁ (2.0 ±0.5 min^?1) in water and seawater are similar within experimental error. The value of k₂ (1.2 × 10⁵ min^?1) in seawater is 28% of its value in water in reasonable agreement with predictions using an ion pairing model.For the oxidation of Cu(I) a rate equation of the form

$\text{?d ln [Cu(I)]/dt = k}{}_0\text{α}{}_{\mathrm{Cu}}\text{+ k}{}_1\text{β}{}_1\text{α}{}_{\mathrm{Cu}}\text{[Cl]}$

was found where k₀ (14.1 sec^?1) and k₁ (3.9 sec^?1) are the pseudo first order rate constants for the oxidation of Cu⁺ and CuCl⁰, β₁ is the formation constant for CuCl⁰ and α_Cu is the fraction of free Cu⁺. Thus, unlike the results for Fe(II), Cu(I) chloride complexes have measurable rates of oxidation.