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 CO2 gas generates an acidic fluid and produces secondary quartz, muscovite and/or pyrophyllite at constant temperature and pressure of 300 degrees 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.

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.

A numerical study of flexural buckling of foliated rock slopes

The occurrence of foliated rock masses is common in mining environment. Methods employing continuum approximation in describing the deformation of such rock masses possess a clear advantage over methods where each rock layer and each inter-layer interface (joint) is explicitly modelled. In devising such a continuum model it is imperative that moment (couple) stresses and internal rotations associated with the bending of the rock layers be properly incorporated in the model formulation. Such an approach will lead to a Cosserat-type theory. In the present model, the behaviour of the intact rock layer is assumed to be linearly elastic and the joints are assumed to be elastic-perfectly plastic. Condition of slip at the interfaces are determined by a Mohr-Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformations. The model is incorporated into the finite element program AFENA and validated against an analytical solution of elementary buckling problems of a layered medium under gravity loading. A design chart suitable for assessing the stability of slopes in foliated rock masses against flexural buckling failure has been developed. The design chart is easy to use and provides a quick estimate of critical loading factors for slopes in foliated rock masses. It is shown that the model based on Euler's buckling theory as proposed by Cavers (Rock Mechanics and Rock Engineering 1981; 14:87-104) substantially overestimates the critical heights for a vertical slope and underestimates the same for sub-vertical slopes. Copyright (C) 2001 John Wiley & Sons, Ltd.

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.

Coaxial flattening at deep levels of orogenic belts: evidence from blueschists and eclogites on Syros and Sifnos (Cyclades, Greece)

This work presents new Structural data from a high-pressure/low-temperature (HP/LT) metamorphic terrane exposed on the islands of Syros and Sifnos (Cyclades, Greece). The structure and the metamorphism of a relatively coherent HP/LT rock section were studied in order to elucidate how strain was accommodated at deep crustal levels during the formation and exhumation of HP/LT rocks. At least three deformation phases associated with eclogite- and blueschist-facies conditions (P = 8-15 kbar; T = 400-550 degreesC) were recognised. The earliest deformation fabric (S1), preserved as inclusion trails within garnet porphyroblasts, is aligned to define a sub-vertical schistosity (at present orientation), which is frequently orthogonal to the flat matrix schistosity (S2), and may indicate that deep crustal thickening involved upright folding. The currently dominant fabric in the HP rock section, S2, is Usually moderately dipping and locally contains NW-trending glaucophane lineations, symmetric pressure-shadows and eclogitic boudins. The symmetric structures associated with this fabric seem to indicate coaxial vertical thinning, although the existence of non-coaxial structures out of the study area cannot be excluded. Glaucophane-bearing shear bands (S3), with top-to-NW sense of shearing, locally crosscut the earlier structures. The latest recognised fabric (D4) is scarce and often absent within the HP rocks. It is associated with top-to-NE kinematic criteria that formed at greenschist-facies conditions (P = 4-7 kbar; T = 400-450 degreesC). Based on these observations, it is suggested that partitioning of strain occurred at different crustal levels and at different times. Deep crustal deformation was governed by thickening via upright folding followed by coaxial vertical thinning, whereas non-coaxial shearing occurred when the rocks were already exhumed to relatively shallow crustal levels. The earliest fabrics (D1 to D3) pertain to Alpine orogenesis and possibly to syn-orogenic extension, whereas the latest correspond to whole-crust back-are extension. (C) 2002 Elsevier Science Ltd. All rights reserved.

A director theory for visco-elastic folding instabilities in multilayered rock

A model for finely layered visco-elastic rock proposed by us in previous papers is revisited and generalized to include couple stresses. We begin with an outline of the governing equations for the standard continuum case and apply a computational simulation scheme suitable for problems involving very large deformations. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered beam under compression. We analyse folding up to 40% shortening. The standard continuum solution becomes unstable for extreme values of the shear/normal viscosity ratio. The instability is a consequence of the neglect of the bending stiffness/viscosity in the standard continuum model. We suggest considering these effects within the framework of a couple stress theory. Couple stress theories involve second order spatial derivatives of the velocities/displacements in the virtual work principle. To avoid C-1 continuity in the finite element formulation we introduce the spin of the cross sections of the individual layers as an independent variable and enforce equality to the spin of the unit normal vector to the layers (-the director of the layer system-) by means of a penalty method. We illustrate the convergence of the penalty method by means of numerical solutions of simple shears of an infinite layer for increasing values of the penalty parameter. For the shear problem we present solutions assuming that the internal layering is oriented orthogonal to the surfaces of the shear layer initially. For high values of the ratio of the normal-to the shear viscosity the deformation concentrates in thin bands around to the layer surfaces. The effect of couple stresses on the evolution of folds in layered structures is also investigated. (C) 2002 Elsevier Science Ltd. All rights reserved.

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.

Environmental influence on the stress corrosion cracking of rock bolts

In order to understand rock bolt Stress Corrosion Cracking (SCC), a series of experiments have been performed in Linearly Increasing Stress Test (LIST) apparatus. One series of experiments determined the threshold stress of various bolt metallurgies (900 MPa for Steel A, and 800 MPa for Steel B and C). The high values of threshold stress suggest that SCC begins in rock bolts when they are sheared by moving rock strata. Typical crack velocity values have been measured to be 2.5 x 10(-8) m s(-1), indicating that there is not much benefit for rock bolt steel of higher fracture toughness. Another series of experiments were performed to understand the environmental conditions causing SCC of steel A and galvanised Steel A rock bolt steel. SCC only occurred for environmental conditions for which produce hydrogen on the sample surface, leading to hydrogen embrittlement and SCC. Fracture surfaces of LIST samples failed by SCC were found to display the same fracture regions as fracture surfaces of rock bolts failed in service by SCC: Tearing Topography Surface (TTS), Corrugated Irregular Surface (CIS), quasi Micro Void Coalescence (qMVC) and Fast Fracture Surface (FFS). Water chemistry analysis were carried out on samples collected from various Australian mines in order to compare laboratory electrolyte conditions to those found in underground mines.

Royal Australian and New Zealand College of Psychiatrists clinical practice guidelines for the treatment of schizophrenia and related disorders

Background: The Royal Australian and New Zealand College of Psychiatrists is co-ordinating the development of clinical practice guidelines (CPGs) in psychiatry, funded under the National Mental Health Strategy (Australia) and the New Zealand Health Funding Authority. This paper presents CPGs for schizophrenia and related disorders. Over the past decade schizophrenia has become more treatable than ever before. A new generation of drug therapies, a renaissance of psychological and psychosocial interventions and a first generation of reform within the specialist mental health system have combined to create an evidence-based climate of realistic optimism. Progressive neuroscientific advances hold out the strong possibility of more definitive biological treatments in the near future. However, this improved potential for better outcomes and quality of life for people with schizophrenia has not been translated into reality in Australia. The efficacy-effectiveness gap is wider for schizophrenia than any other serious medical disorder. Therapeutic nihilism, under-resourcing of services and a stalling of the service reform process, poor morale within specialist mental health services, a lack of broad-based recovery and life support programs, and a climate of tenacious stigma and consequent lack of concern for people with schizophrenia are the contributory causes for this failure to effectively treat. These guidelines therefore tackle only one element in the endeavour to reduce the impact of schizophrenia. They distil the current evidence-base and make recommendations based on the best available knowledge. Method: A comprehensive literature review (1990-2003) was conducted, including all Cochrane schizophrenia reviews and all relevant meta-analyses, and a number of recent international clinical practice guidelines were consulted. A series of drafts were refined by the expert committee and enhanced through a bi-national consultation process. Treatment recommendations: This guideline provides evidence-based recommendations for the management of schizophrenia by treatment type and by phase of illness. The essential features of the guidelines are: (i) Early detection and comprehensive treatment of first episode cases is a priority since the psychosocial and possibly the biological impact of illness can be minimized and outcome improved. An optimistic attitude on the part of health professionals is an essential ingredient from the outset and across all phases of illness. (ii) Comprehensive and sustained intervention should be assured during the initial 3-5 years following diagnosis since course of illness is strongly influenced by what occurs in this 'critical period'. Patients should not have to 'prove chronicity' before they gain consistent access and tenure to specialist mental health services. (iii) Antipsychotic medication is the cornerstone of treatment. These medicines have improved in quality and tolerability, yet should be used cautiously and in a more targeted manner than in the past. The treatment of choice for most patients is now the novel antipsychotic medications because of their superior tolerability and, in particular, the reduced risk of tardive dyskinesia. This is particularly so for the first episode patient where, due to superior tolerability, novel agents are the first, second and third line choice. These novel agents are nevertheless associated with potentially serious medium to long-term side-effects of their own for which patients must be carefully monitored. Conventional antipsychotic medications in low dosage may still have a role in a small proportion of patients, where there has been full remission and good tolerability; however, the indications are shrinking progressively. These principles are now accepted in most developed countries. (vi) Clozapine should be used early in the course, as soon as treatment resistance to at least two antipsychotics has been demonstrated. This usually means incomplete remission of positive symptomatology, but clozapine may also be considered where there are pervasive negative symptoms or significant or persistent suicidal risk is present. (v) Comprehensive psychosocial interventions should be routinely available to all patients and their families, and provided by appropriately trained mental health professionals with time to devote to the task. This includes family interventions, cognitive-behaviour therapy, vocational rehabilitation and other forms of therapy, especially for comorbid conditions, such as substance abuse, depression and anxiety. (vi) The social and cultural environment of people with schizophrenia is an essential arena for intervention. Adequate shelter, financial security, access to meaningful social roles and availability of social support are essential components of recovery and quality of life. (vii) Interventions should be carefully tailored to phase and stage of illness, and to gender and cultural background. (viii) Genuine involvement of consumers and relatives in service development and provision should be standard. (ix) Maintenance of good physical health and prevention and early treatment of serious medical illness has been seriously neglected in the management of schizophrenia, and results in premature death and widespread morbidity. Quality of medical care for people with schizophrenia should be equivalent to the general community standard. (x) General practitioners (GPs)s should always be closely involved in the care of people with schizophrenia. However, this should be truly shared care, and sole care by a GP with minimal or no special Optimal treatment of schizophrenia requires a multidisciplinary team approach with a consultant psychiatrist centrally involved.

Material influence on the stress corrosion cracking of rock bolts

Rock bolt stress corrosion cracking (SCC) has been investigated using the linearly increasing stress test (LIST). One series of experiments determined the threshold stress of various bolt metallurgies (900 MPa for 1355AXRC, and 800 MPa for MAC and MA840B steels). The high values of threshold stress suggest that SCC begins in rock bolts when they are sheared by moving rock strata. SCC only occurred for environmental conditions which produce hydrogen on the sample surface, leading to hydrogen embrittlement and SCC. Different threshold potentials were determined for a range of metallurgies. Cold work was shown to increase the resistance of the steel to SCC. Rock bolt rib geometry does not have a direct impact on the SCC resistance properties of the bolt, although the process by which the ribs are produced can introduce tensile stresses into the bolt which lower its resistance to SCC. (C) 2004 Elsevier Ltd. All rights reserved.

Incidence and distribution of viruses of Taro (Colocasia esculenta) in Pacific Island countries

Four viruses have been reported from taro; Dasheen mosaic virus (DsMV), Taro bacilliform virus (TaBV) and two putative rhabdoviruses, Colocasia bobone disease virus (CBDV) and Taro vein chlorosis virus (TaVCV). A fifth virus, tentatively named Taro reovirus (TaRV), has also been recently identified. The distribution of these viruses throughout the Pacific Islands, and the symptoms associated with their infection, are unknown in many cases due to a lack of sensitive diagnostic tests. We have used recently developed PCR-based diagnostic tests to survey taro growing in 11 Pacific Island countries for the presence of known viruses. DsMV and TaBV were widespread, whereas TaVCV and TaRV were more restricted in their distribution. CBDV was restricted to PNG and Solomon Islands and was always associated with the two most serious viral diseases of taro; alomae disease and bobone disease, but the causal agent of these two diseases remains unclear.