Prediction of rock alteration patterns: A potential tool in mineral exploration

Hydrothermal alteration of a quartz‐K‐feldspar rock is simulated numerically by coupling fluid flow and chemical reactions. Introduction of CO₂ gas generates an acidic fluid and produces secondary quartz, muscovite and/or pyrophyllite at constant temperature and pressure of 300°C and 200 MPa. The precipitation and/or dissolution of the secondary minerals is controlled by either mass‐action relations or rate laws. In our simulations the mass of the primary elements are conserved and the mass‐balance equations are solved sequentially using an implicit scheme in a finite‐element code. The pore‐fluid velocity is assumed to be constant. The change of rock volume due to the dissolution or precipitation of the minerals, which is directly related to their molar volume, is taken into account. Feedback into the rock porosity and the reaction rates is included in the model. The model produces zones of pyrophyllite quartz and muscovite due to the dissolution of K‐feldspar. Our model simulates, in a simplified way, the acid‐induced alteration assemblages observed in various guises in many significant mineral deposits. The particular aluminosilicate minerals produced in these experiments are associated with the gold deposits of the Witwatersrand Basin.     

Fractionation of hydrogen and oxygen isotopes of gypsum hydration water and assessment of its geochemical indications

Fundamental knowledge of the isotopic fractionation between the hydration water and the mother solution and whether the primary information recorded in hydration water can be preserved or not in deposits or mines have long been unclear. In order to calculate the accurate hydrogen and oxygen isotopic fractionation factors between gypsum hydration water and its mother solution with new methods, to understand the mechanism of fractionation and synthetically assess the record-keeping abilities of the isotopic composition of hydration water during the process of diagenesis after deposition, experiments on the hydrogen and oxygen isotopic compositions of gypsum hydration water and its mother solution at different isothermal temperatures from 5 to 50°C were systematically conducted. In addition, samples from two typical gypsum deposits formed in different environmental conditions were also determined. Results show that during gypsum crystallisation, both hydrogen and oxygen isotopes show significant fractionation between the hydration water and the mother solution. The calculated hydrogen isotopic fractionation factors are <1, while the oxygen isotopic fractionation factors are >1 at temperatures from 5 to 50°C. The fractionation factors show no functional relationships with temperature. Isotopic compositions of gypsum hydration water in arid lake sediments can be used to trace the source of water and primary deposit environmental information. However, the isotopic composition of the gypsum hydration water can easily be altered by dissolution and secondary precipitation of gypsum during later diagenesis, particularly in areas with humid climate and abundant groundwater. A very careful assessment on record-keeping abilities of the primary isotopic composition of hydration water in gypsum during later diagenesis must be considered before application.     

The problem of scaling in grid-related hydrological process modelling

This article shows some important scaling aspects that exist in hydrological process grid modelling with GIS, with respect to selecting the optimal grid width. Hydrologically relevant space parameter fields should often be regarded as self-similar fractals for which certain invariance characteristics are applied. These characteristics are especially important if there is an extensive change in the model grid width. The self-similarity generates simple scaling laws, which can be used for rescaling procedures. © 1997 John Wiley & Sons, Ltd.     

Evolution of localized folding for a thin elastic layer in a softening visco-elastic medium

Long compressed elastic struts on softening elastic foundations have a tendency to buckle locally. The same tendency is demonstrated here for the instantaneous response of elastic struts supported by visco-elastic media. A governing nonlinear partial differential equation is derived to describe the evolution of the localized form in time. Under the assumed constant end-shortening this is found to be approximated by a coupled set of seven ordinary differential (diffusion) equations. As the load drops to zero, the localized buckle pattern evolves towards the form of the single long wave, but remains aperiodic for all time. Three-dimensional plots show how this localized pattern changes over time.