Geochronology of weathering and landscape evolution, Dugald River valley, NW Queensland, Australia

40Ar/39Ar laser incremental heating analyses of individual grains of supergene jarosite, alunite, and cryptomelane from weathering profiles in the Dugald River area, Queensland, Australia, show a strong positive correlation between a sample’s age and its elevation. We analyzed 125 grains extracted from 35 hand specimens collected from weathering profiles at 11 sites located at 3 distinct elevations. The highest elevation profile hosts the oldest supergene minerals, whereas progressively younger samples occur at lower positions in the landscape. The highest elevation sampling sites (three sites), located on top of an elongated mesa (255 to 275 m elevation), yield ages in the 16 to 12 Ma range. Samples from an intermediate elevation site (225 to 230 m elevation) yield ages in the 6 to 4 Ma range. Samples collected at the lowest elevation sites (200 to 220 m elevation) yield ages in the 2.2 to 0.8 Ma interval. Grains of supergene alunite, jarosite, and cryptomelane analyzed from individual single hand specimens yield reproducible results, confirming the suitability of these minerals to 40Ar/39Ar geochronology. Multiple samples collected from the same site also yield reproducible results, indicating that the ages measured are true precipitation ages for the samples analyzed. Different sites, up to 3 km apart, sampled from weathering profiles at the same elevation again yield reproducible results. The consistency of results confirms that 40Ar/39Ar geochronology of supergene jarosite, alunite, and cryptomelane yields ages of formation of weathering profiles, providing a reliable numerical basis for differentiating and correlating these profiles. The age versus elevation relationship obtained suggest that the stepped landscapes in the Dugald River area record a progressive downward migration of a relatively flat weathering front. The steps in the landscape result from differential erosion of previously weathered bedrock displaying different susceptibility to weathering and contrasting resistance to erosion. Combined, the age versus elevation relationships measured yield a weathering rate of 3.8 m. Myr−1 (for the past 15 Ma) if a descending subhorizontal weathering front is assumed. The results also permit the calculation of the erosion rate of the more easily weathered and eroded lithologies, assuming an initially flat landscape as proposed in models of episodic landscape development. The average erosion rate for the past 15 Ma is 3.3 m. Myr−1, consistent with erosion rates obtained by cosmogenic isotope studies in the region.

Sedimentologic, petrographic, and sulfur isotope constraints on fine-grained pyrite formation at Mount Isa Mine and environs, Northwest Queensland, Australia

Fine-grained pyrite is the earliest generation of pyrite and the most abundant sulfide within the Urquhart Shale at Mount Isa, northwest Queensland. The pyrite is intimately interbanded with ore-grade Pb-Zn miner alization at the Mount Isa mine but is also abundant north and south of the mine at several stratigraphic horizons within the Urquhart Shale. Detailed sedimentologic, petrographic, and sulfur isotope studies of the Urquhart Shale, mostly north of the mine, reveal that the fine-grained pyrite (delta(34)S = -3.3 to +26.3 parts per thousand) formed by thermochemical sulfate reduction during diagenesis. The sulfate source was local sulfate evaporites, pseudo morphs of which are present throughout the Urquhart Shale (i.e., gypsum, anhydrite, and barite). Deep-burial diagenetic replacement of these evaporites resulted in sulfate-bearing ground waters which migrated parallel to bedding. Fine-grained pyrite formed where these fluids infiltrated and then interacted with carbon-rich laminated siltstones. Comparison of the sulfur isotope systematics of fine-grained pyrite and spatially associated base metal sulfides from the Mount Isa Pb-Zn and Cu orebodies indicates a common sulfur source of ultimately marine origin for all sulfide types. Different sulfur isotope ratio distributions for the various sulfides are the result of contrasting formation mechanisms and/or depositional conditions rather than differing sulfur sources. The sulfur isotope systematics of the base metal and associated iron sulfide generations are consistent with mineralization by reduced hydrothermal fluids, perhaps generated by bulk reduction of evaporite-sourced sulfate-bearing waters generated deeper in the Mount Isa Group, the sedimentary sequence which contains the Urquhart Shale. The available sulfur isotope data from the Mount Isa orebodies are consistent with either a chemically and thermally zoned, evolving Cu-Pb-Zn system, or discrete Cu and Pb-Zn mineralizing events linked by a common sulfur source.

Vascular smooth muscle and arterial calcification

Smooth muscle cultures can calcify under certain circumstances. As a model system these cultures therefore provide information on why calcification occurs in atherosclerotic plaques. Whether all smooth muscle cells (under certain conditions), or only specific populations, can produce this mineralization has not been resolved. Demer's group has cloned calcifying vascular cells from subcultured bovine aorta and studied them in detail. They have speculated on whether the cells are smooth muscle which have altered in phenotype, or whether they are derived from a stem cell population within the artery wall. The article argues that while the normal process of smooth muscle phenotypic modulation seen in arterial repair could account for the observations, this view may be two simplistic considering the complex nature of the artery wall. Certainly there is evidence for heterogeneity of smooth muscle cells in the artery wall and recent evidence suggests that stem cells can circulate in the blood and repopulate tissues. Further studies are required to resolve the important question as to the origin of cells which produce mineralization in atheroma.

Finite element modelling of dissipative structures for nonequilibrium chemical reactions in fluid-saturated porous media

We use the finite element method to model and predict the dissipative structures of chemical species for a nonequilibrium chemical reaction system in a fluid-saturated porous medium. In particular, we explore the conditions under which dissipative structures of the species may exist in the Brusselator type of nonequilibrium chemical reaction. Since this is the first time the finite element method and related strategies have been used to study the chemical instability problems in a fluid-saturated porous medium, it is essential to validate the method and strategies before they are put into application. For this purpose, we have rigorously derived the analytical solutions for dissipative structures of chemical species in a benchmark problem, which geometrically is a square. Comparison of the numerical solutions with the analytical ones demonstrates that the proposed numerical method and strategy are robust enough to solve chemical instability problems in a fluid-saturated porous medium. Finally, the related numerical results from two application examples indicate that both the regime and the magnitude of pore-fluid flow have significant effects on the nature of the dissipative structures that developed for a nonequilibrium chemical reaction system in a fluid-saturated porous medium. The motivation for this study is that self-organization under conditions of pore-fluid flow in a porous medium is a potential mechanism of the orebody formation and mineralization in the upper crust of the Earth. (C) 2000 Elsevier Science S.A. All rights reserved.

Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

Finite element modelling of three-dimensional convection problems in fluid-saturated porous media heated from below

We use the finite element method to model three-dimensional convective pore-fluid flow in fluid-saturated porous media when they are heated from below. In particular, we employ the particle-tracking technique to mimic the trajectories of particles in three-dimensional fluid flow problems. The related numerical results demonstrated that: (1) The progressive asymptotic approach procedure, which was previously developed for the finite element modelling of two-dimensional convective pore-fluid flow problems, is equally applicable to the finite element modelling of three-dimensional convective pore-fluid flow in fluid-saturated porous media heated from below. (2) The perturbation of gravity at different planes has a significant effect on the pattern of three-dimensional convective pore-fluid flow and therefore, may influence the pattern of orebody formation and mineralization in three-dimensional hydrothermal systems. Copyright (C) 2001 John Wiley & Sons, Ltd.

Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

Finite element modelling of rock alteration and metamorphic process in hydrothermal systems

We use the finite element method to simulate the rock alteration and metamorphic process in hydrothermal systems. In particular, we consider the fluid-rock interaction problems in pore-fluid saturated porous rocks. Since the fluid rock interaction takes place at the contact interface between the pore-fluid and solid minerals, it is governed by the chemical reaction which usually takes place very slowly at this contact interface, from the geochemical point of view. Due to the relative slowness of the rate of the chemical reaction to the velocity of the pore-fluid flow in the hydrothermal system to be considered, there exists a retardation zone, in which the conventional static theory in geochemistry does not hold true. Since this issue is often overlooked by some purely numerical modellers, it is emphasized in this paper. The related results from a typical rock alteration and metamorphic problem in a hydrothermal system have shown not only the detailed rock alteration and metamorphic process, but also the size of the retardation zone in the hydrothermal system. Copyright (C) 2001 John Wiley & Sons, Ltd.

Finite element modelling of heat transfer through permeable cracks in hydrothermal systems with upward throughflow

We use the finite element method to model the heat transfer phenomenon through permeable cracks in hydrothermal systems with upward throughflow. Since the finite element method is an approximate numerical method, the method must be validated before it is used to soh,e any new, kind of problem. However, the analytical solution, which can be used to validate the finite element method and other numerical methods, is rather limited in the literature, especially, for the problem considered here. Keeping this in mind, we have derived analytical solutions for the temperature distribution along the vertical axis of a crack in a fluid-saturated porous layer. After the finite element method is validated by comparing the numerical solution with the analytical solution for the same benchmark problem, it is used to investigate the pore-fluid flow and heat transfer in layered hydrothermal systems with vertical permeable cracks. The related analytical and numerical results have demonstrated that vertical cracks are effective and efficient members to transfer heat energy from the bottom section to the top section in hydrothermal systems with upward throughflow.

Computer simulations of coupled problems in geological and geochemical systems

The finite element method is used to simulate coupled problems, which describe the related physical and chemical processes of ore body formation and mineralization, in geological and geochemical systems. The main purpose of this paper is to illustrate some simulation results for different types of modelling problems in pore-fluid saturated rock masses. The aims of the simulation results presented in this paper are: (1) getting a better understanding of the processes and mechanisms of ore body formation and mineralization in the upper crust of the Earth; (2) demonstrating the usefulness and applicability of the finite element method in dealing with a wide range of coupled problems in geological and geochemical systems; (3) qualitatively establishing a set of showcase problems, against which any numerical method and computer package can be reasonably validated. (C) 2002 Published by Elsevier Science B.V.

Analysis of steady-state heat transfer through mid-crustal vertical cracks with upward throughflow in hydrothermal systems

We conduct a theoretical analysis of steady-state heat transfer problems through mid-crustal vertical cracks with upward throughflow in hydrothermal systems. In particular, we derive analytical solutions for both the far field and near field of the system. In order to investigate the contribution of the forced advection to the total temperature of the system, two concepts, namely the critical Peclet number and the critical permeability of the system, have been presented and discussed in this paper. The analytical solution for the far field of the system indicates that if the pore-fluid pressure gradient in the crust is lithostatic, the critical permeability of the system can be used to determine whether or not the contribution of the forced advection to the total temperature of the system is negligible. Otherwise, the critical Peclet number should be used. For a crust of moderate thickness, the critical permeability is of the order of magnitude of 10(-20) m(2), under which heat conduction is the overwhelming mechanism to transfer heat energy, even though the pore-fluid pressure gradient in the crust is lithostatic. Furthermore, the lower bound analytical solution for the near field of the system demonstrates that the permeable vertical cracks in the middle crust can efficiently transfer heat energy from the lower crust to the upper crust of the Earth. Copyright (C) 2002 John Wiley Sons, Ltd.

Effects of hot intrusions on pore-fluid flow and heat transfer in fluid-saturated rocks

We use the finite element method to solve coupled problems between pore-fluid flow and heat transfer in fluid-saturated porous rocks. In particular, we investigate the effects of both the hot pluton intrusion and topographically driven horizontal flow on the distributions of the pore-flow velocity and temperature in large-scale hydrothermal systems. Since general mineralization patterns are strongly dependent on distributions of both the pore-fluid velocity and temperature fields, the modern mineralization theory has been used to predict the general mineralization patterns in several realistic hydrothermal systems. The related numerical results have demonstrated that: (1) The existence of a hot intrusion can cause an increase in the maximum value of the pore-fluid velocity in the hydrothermal system. (2) The permeability of an intruded pluton is one of the sensitive parameters to control the pore-fluid flow, heat transfer and ore body formation in hydrothermal systems. (3) The maximum value of the pore-fluid velocity increases when the bottom temperature of the hydrothermal system is increased. (4) The topographically driven flow has significant effects on the pore-fluid flow, temperature distribution and precipitation pattern of minerals in hydrothermal systems. (5) The size of the computational domain may have some effects on the pore-fluid flow and heat transfer, indicating that the size of a hydrothermal system may affect the pore-fluid flow and heat transfer within the system. (C) 2003 Elsevier Science B.V. All rights reserved.

An equivalent algorithm for simulating thermal effects of magma intrusion problems in porous rocks

An equivalent algorithm is proposed to simulate thermal effects of the magma intrusion in geological systems, which are composed of porous rocks. Based on the physical and mathematical equivalence, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with a physically equivalent heat source. From the analysis of an ideal solidification model, the physically equivalent heat source has been determined in this paper. The major advantage in using the proposed equivalent algorithm is that the fixed finite element mesh with a variable integration time step can be employed to simulate the thermal effect of the intruded magma solidification using the conventional finite element method. The related numerical results have demonstrated the correctness and usefulness of the proposed equivalent algorithm for simulating the thermal effect of the intruded magma solidification in geological systems. (C) 2003 Elsevier B.V. All rights reserved.

Theoretical and numerical analyses of convective instability in porous media with temperature-dependent viscosity

Exact analytical solutions of the critical Rayleigh numbers have been obtained for a hydrothermal system consisting of a horizontal porous layer with temperature-dependent viscosity. The boundary conditions considered are constant temperature and zero vertical Darcy velocity at both the top and bottom of the layer. Not only can the derived analytical solutions be readily used to examine the effect of the temperature-dependent viscosity on the temperature-gradient driven convective flow, but also they can be used to validate the numerical methods such as the finite-element method and finite-difference method for dealing with the same kind of problem. The related analytical and numerical results demonstrated that the temperature-dependent viscosity destabilizes the temperature-gradient driven convective flow and therefore, may affect the ore body formation and mineralization in the upper crust of the Earth. Copyright (C) 2003 John Wiley Sons, Ltd.