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.

Finite element modelling of reactive fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins

We present the finite element simulations of reactive mineral carrying fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins. In particular we explore the mixing of reactive sulfide and sulfate fluids and the relevant patterns of mineralization for Load, zinc and iron minerals in the regime of temperature-gradient-driven convective flow. Since the mineralization and ore body formation may last quite a long period of time in a hydrothermal basin, it is commonly assumed that, in the geochemistry, the solutions of minerals are in an equilibrium state or near an equilibrium state. Therefore, the mineralization rate of a particular kind of mineral can be expressed as the product of the pore-fluid velocity and the equilibrium concentration of this particular kind of mineral Using the present mineralization rate of a mineral, the potential of the modern mineralization theory is illustrated by means of finite element studies related to reactive mineral-carrying fluids mixing problems in materially homogeneous and inhomogeneous porous rock basins.

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.

Experimental study of the deposition of combustion aerosols in the human respiratory tract

The total deposition of environmental tobacco smoke (ETS), diesel and petrol smoke in the respiratory tract of 14 non-smokers between the ages of 20 and 30 was determined experimentally. A scanning mobility particle sizer (SMPS) measuring a size range of 0.016-0.626 mu m was used to characterise the inhaled and exhaled aerosol during relaxed nasal breathing over a period of 10 min. The ETS, diesel, and petrol particles had average count median diameter (and geometric standard deviation) of 0.183 mu m (1.7), 0.125 mu m (1.7), and 0.069 mu m (1.7), respectively. The average total number deposition of ETS was 36% (standard deviation 10%), of diesel smoke 30% (standard deviation 9%), and of petrol smoke 41% (standard deviation 8%). The analysis of the deposition patterns as a function of particle size for the three aerosols in each individual showed that there is a significant difference between each aerosol for a majority of individuals (12 out of 14). This is an important result as it indicates that differences persist regardless of inter-subject variability. (C) 2005 Elsevier Ltd. All rights reserved.

The geochemistry of late Archaean microbial carbonate: Implications for ocean chemistry and continental erosion history

Trace element concentrations and combined Sr- and Nd-isotope compositions were determined on stromatolitic carbonates (microbialites) from the 2.52 Ga Campbellrand carbonate platform (South Africa). Shale-normalised rare earth element and yttrium patterns of the ancient samples are similar to those of modern seawater in having positive La and Y anomalies and in being depleted in light rare earth elements. In contrast to modem seawater (and microbialite proxies), the 2.52 Ga samples lack a negative Ce anomaly but possess a positive Eu anomaly. These latter trace element characteristics are interpreted to reflect anoxic deep ocean waters where, unlike today, hydrothermal Fe input was not oxidised, and scavenged and rare earth elements were not coprecipitated with Fe-oxyhydroxides. The persistence of a positive Eu anomaly in relatively shallow Campbellrand platform waters indicates a dramatic reversal from hydrothermally dominated (Archaean) to continental erosion-dominated (Phanerozoic) rare earth element flux ratio. The dominant hydrothermal input is also expressed in the initial Sr- and Nd-isotope ratios. There is collinear variation in Sr-Nd systematics, which range from primitive values (Sr-87/Sr-86 of 0.702386 and epsilon (Nd) of +2.1) to more evolved crustal ratios. Mixing calculations show that the range in trace element ratios (e.g., Y/Ho) and initial isotope ratios is not a result of contamination by trapped sediment, but that the chemical band isotopic variation reflects carbonate deposition in an environment where different water masses mixed. Calculated Nd flux ratios yield a hydrothermal input into the 2.52 Ga oceans one order of magnitude larger than continental input. Such a change in flux ratio most likely required substantially reduced continental inputs, which could, in turn, reflect a plate tectonic causation (e.g., reduced topography or expansion of epicontinental seas). Copyright (C) 2001 Elsevier Science Ltd.