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.

Large amplitude folding in finely layered viscoelastic rock structures

We analyze folding phenomena in finely layered viscoelastic rock. Fine is meant in the sense that the thickness of each layer is considerably smaller than characteristic structural dimensions. For this purpose we derive constitutive relations and apply a computational simulation scheme (a finite-element based particle advection scheme; see MORESI et al., 2001) suitable for problems involving very large deformations of layered viscous and viscoelastic rocks. An algorithm for the time integration of the governing equations as well as details of the finite-element implementation is also given. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered plate. The domain is compressed along the layer axis by prescribing velocities along the sides. First, for the viscous limit we consider the response to a series of harmonic perturbations of the director orientation. The Fourier spectra of the initial folding velocity are compared for different viscosity ratios. Turning to the nonlinear regime we analyze viscoelastic folding histories up to 40% shortening. The effect of layering manifests itself in that appreciable buckling instabilities are obtained at much lower viscosity ratios (1:10) as is required for the buckling of isotropic plates (1:500). The wavelength induced by the initial harmonic perturbation of the director orientation seems to be persistent. In the section of the parameter space considered here elasticity seems to delay or inhibit the occurrence of a second, larger wavelength. Finally, in a linear instability analysis we undertake a brief excursion into the potential role of couple stresses on the folding process. The linear instability analysis also provides insight into the expected modes of deformation at the onset of instability, and the different regimes of behavior one might expect to observe.

Combined chemical separation of Lu, Hf, Sm, Nd, and REEs from a single rock digest: Precise and accurate isotope Determinations of Lu-Hf and Sm-Nd using multicollector-ICPMS

A combined procedure for separating Lu, Hf, Sm, Nd, and rare earth elements (REEs) from a single sample digest is presented. The procedure consists of the following five steps: (1) sample dissolution via sodium peroxide sintering; (2) separation of the high field strength elements from the REEs and other matrix elements by a HF-free anion-exchange column procedure; (3) purification of Hf on a cation-exchange resin; (4) separation of REEs from other matrix elements by cation exchange; (5) Lu, Sm, and Nd separation from the other REEs by reversed-phase ion chromatography. Analytical reproducibilities of Sm-Nd and Lu-Hf isotope systematics are demonstrated for standard solutions and international rock reference materials. Results show overall good reproducibilities for Sm-Nd systematics independent of the rock type analyzed. For the Lu-Hf systematics, the reproducibility of the parent/daughter ratio is much better for JB-1 (basalt) than for two analyzed felsic crustal rocks (DR-N and an Archaean granitoid). It is demonstrated that this poorer reproducibility of the Lu/Hf ratio is truly caused by sample heterogeneity; thus, results are geologically reasonable.