Longevity and stability of cratonic lithosphere: Insights from numerical simulations of coupled mantle convection and continental tectonics

[1] The physical conditions required to provide for the tectonic stability of cratonic crust and for the relative longevity of deep cratonic lithosphere within a dynamic, convecting mantle are explored through a suite of numerical simulations. The simulations allow chemically distinct continents to reside within the upper thermal boundary layer of a thermally convecting mantle layer. A rheologic formulation, which models both brittle and ductile behavior, is incorporated to allow for plate-like behavior and the associated subduction of oceanic lithosphere. Several mechanisms that may stabilize cratons are considered. The two most often invoked mechanisms, chemical buoyancy and/or high viscosity of cratonic root material, are found to be relatively ineffective if cratons come into contact with subduction zones. High root viscosity can provide for stability and longevity but only within a thick root limit in which the thickness of chemically distinct, high-viscosity cratonic lithosphere exceeds the thickness of old oceanic lithosphere by at least a factor of 2. This end-member implies a very thick mechanical lithosphere for cratons. A high brittle yield stress for cratonic lithosphere as a whole, relative to oceanic lithosphere, is found to be an effective and robust means for providing stability and lithospheric longevity. This mode does not require exceedingly deep strength within cratons. A high yield stress for only the crustal or mantle component of the cratonic lithosphere is found to be less effective as detachment zones can then form at the crust-mantle interface which decreases the longevity potential of cratonic roots. The degree of yield stress variations between cratonic and oceanic lithosphere required for stability and longevity can be decreased if cratons are bordered by continental lithosphere that has a relatively low yield stress, i.e., mobile belts. Simulations that combine all the mechanisms can lead to crustal stability and deep root longevity for model cratons over several mantle overturn times, but the dominant stabilizing factor remains a relatively high brittle yield stress for cratonic lithosphere.

Importance of water for Archaean granitoid petrology: a comparative study of TTG and potassic granitoids from Barberton Mountain Land, South Africa

A new model for Archaean granitoid magmatism is presented which reconciles the most important geochemical similarities and differences between tonalite-trondhjemite-granodiorite (TTG) and potassic granitoids. Trace element abundances reveal a strong arc magmatism signature in all studied granitoids from Barberton Mountain Land. Characteristic features include HFSE depletion as well as distinct enrichment peaks of fluid-sensitive trace elements such as Pb in N-MORB normalisation, clearly indicating that all studied granitoids are derived from refertilised mantle above subduction zones. We envisage hydrous basaltic (s.l.) melts as parental liquids, which underwent extensive fractional crystallisation. Distinctive residual cumulates evolved depending on initial water content. High-H2O melts crystallised garnet/amphibole together with pyroxenes and minor plagioclase, but no olivine. This fractionation path ultimately led to TTG-like melts. Less hydrous basaltic melts also crystallised garnet/amphibole, but the lower compatible element content indicates that olivine was also a liquidus phase. Pronounced negative Eu-anomalies of the granitic melts, correlating with Na, Ca and Al, indicate plagioclase to be of major importance. In the context of our model, the post-Archaean disappearance of TTG and concomitant preponderance of granites (s.l.), therefore, is explained with secular decrease of aqueous fluid transport into subduction zones and/or efficiency of deep fluid release from slabs.

Anisotropic viscous models of large-deformation Mohr-Coulomb failure

We have developed a way to represent Mohr-Coulomb failure within a mantle-convection fluid dynamics code. We use a viscous model of deformation with an orthotropic viscoplasticity (a different viscosity is used for pure shear to that used for simple shear) to define a prefered plane for slip to occur given the local stress field. The simple-shear viscosity and the deformation can then be iterated to ensure that the yield criterion is always satisfied. We again assume the Boussinesq approximation, neglecting any effect of dilatancy on the stress field. An additional criterion is required to ensure that deformation occurs along the plane aligned with maximum shear strain-rate rather than the perpendicular plane, which is formally equivalent in any symmetric formulation. We also allow for strain-weakening of the material. The material can remember both the accumulated failure history and the direction of failure. We have included this capacity in a Lagrangian-integration-point finite element code and show a number of examples of extension and compression of a crustal block with a Mohr-Coulomb failure criterion. The formulation itself is general and applies to 2- and 3-dimensional problems.

Subduction of the Nazca Ridge and the Inca Plateau: Insights into the formation of ore deposits in Peru

A large number of ore deposits that formed in the Peruvian Andes during the Miocene (15-5 Ma) are related to the subduction of the Nazea plate beneath the South American plate. Here we show that the spatial and temporal distribution of these deposits correspond with the arrival of relatively buoyant topographic anomalies, namely the Nazca Ridge in central Peru and the now-consumed Inca Plateau in northern Peru, at the subduction zone. Plate reconstruction shows a rapid metallogenic response to the arrival of the topographic anomalies at the subduction trench. This is indicated by clusters of ore deposits situated within the proximity of the laterally migrating zones of ridge subduction. It is accordingly suggested that tectonic changes associated with impingement of the aseismic ridge into the subduction zone may trigger the formation of ore deposits in metallogenically fertile suprasubduction environments. (c) 2005 Elsevier B.V All rights reserved.

Partitioning of deformation within a subduction channel during exhumation of high-pressure rocks: a case study from the Western Alps

The metamorphic belt of the Western Alps was subjected to widespread extensional tectonism at the end of the Eocene (ca. 45-35 Ma). Extension was accommodated by hinterland-directed movements along gently inclined extensional shear zones, which facilitated rapid exhumation of high-pressure and ultra-high-pressure rocks. This deformation resulted in a normal metamorphic sequence. Extension in the inner parts of the Western Alps was coeval with shortening at the front of the belt (foreland-directed thrusts), which took place during decompression, and emplaced higher grade metamorphic units over lower grade metamorphic rocks, thus forming an inverse metamorphic sequence. Two mechanisms for this extensional episode are discussed: (1) collapse of an overthickened lithosphere, and (2) internal readjustments within the orogenic wedge due to subduction channel dynamics. We favour the latter mechanism because it can account for the development of the observed inverse and normal metamorphic sequences along foreland-directed thrusts and hinterland-directed detachments, respectively. This hypothesis is supported by published structural, metamorphic and geochronological data from four geological transects through the Western Alps. This study also emphasizes the importance of post-shearing deformation (e.g. horizontal buckling versus vertical flattening), which can modify the distribution of hinterland- and foreland-directed shear zones in orogenic belts. (c) 2006 Elsevier Ltd. All rights reserved.

Local adaptation and species segregation in two mussel (Mytilus edulis X M. trossulus) hybrid zones

Few marine hybrid zones have been studied extensively, the major exception being the hybrid zone between the mussels Mytilus edulis and M. galloprovincialis in southwestern Europe. Here, we focus on two less studied hybrid zones that also involve Mytilus spp.; M. edulis and M. trossulus are sympatric and hybridize on both western and eastern coasts of the Atlantic Ocean. We review the dynamics of hybridization in these two hybrid zones and evaluate the role of local adaptation for maintaining species boundaries. In Scandinavia, hybridization and gene introgression is so extensive that no individuals with pure M. trossulus genotypes have been found. However, M. trossulus alleles are maintained at high frequencies in the extremely low salinity Baltic Sea for some allozyme genes. A synthesis of reciprocal transplantation experiments between different salinity regimes shows that unlinked Gpi and Pgm alleles change frequency following transplantation, such that post-transplantation allelic composition resembles native populations found in the same salinity. These experiments provide strong evidence for salinity adaptation at Gpi and Pgm (or genes linked to them). In the Canadian Maritimes, pure M. edulis and M. trossulus individuals are abundant, and limited data suggest that M. edulis predominates in low salinity and sheltered conditions, whereas M. trossulus are more abundant on the wave-exposed open coasts. We suggest that these conflicting patterns of species segregation are, in part, caused by local adaptation of Scandinavian M. trossulus to the extremely low salinity Baltic Sea environment.

Convective instability of 3-D fluid-saturated geological fault zones heated from below

We conduct a theoretical analysis to investigate the convective instability of 3-D fluid-saturated geological fault zones when they are heated uniformly from below. In particular, we have derived exact analytical solutions for the critical Rayleigh numbers of different convective flow structures. Using these critical Rayleigh numbers, three interesting convective flow structures have been identified in a geological fault zone system. It has been recognized that the critical Rayleigh numbers of the system have a minimum value only for the fault zone of infinite length, in which the corresponding convective flow structure is a 2-D slender-circle flow. However, if the length of the fault zone is finite, the convective flow in the system must be 3-D. Even if the length of the fault zone is infinite, since the minimum critical Rayleigh number for the 2-D slender-circle flow structure is so close to that for the 3-D convective flow structure, the system may have almost the same chance to pick up the 3-D convective flow structures. Also, because the convection modes are so close for the 3-D convective flow structures, the convective flow may evolve into the 3-D finger-like structures, especially for the case of the fault thickness to height ratio approaching zero. This understanding demonstrates the beautiful aspects of the present analytical solution for the convective instability of 3-D geological fault zones, because the present analytical solution is valid for any value of the ratio of the fault height to thickness. Using the present analytical solution, the conditions, under which different convective flow structures may take place, can be easily determined.

Double diffusion-driven convective instability of three-dimensional fluid-saturated geological fault zones heated from below

We conduct a theoretical analysis to investigate the double diffusion-driven convective instability of three-dimensional fluid-saturated geological fault zones when they are heated uniformly from below. The fault zone is assumed to be more permeable than its surrounding rocks. In particular, we have derived exact analytical solutions to the total critical Rayleigh numbers of the double diffusion-driven convective flow. Using the corresponding total critical Rayleigh numbers, the double diffusion-driven convective instability of a fluid-saturated three-dimensional geological fault zone system has been investigated. The related theoretical analysis demonstrates that: (1) The relative higher concentration of the chemical species at the top of the three-dimensional geological fault zone system can destabilize the convective flow of the system, while the relative lower concentration of the chemical species at the top of the three-dimensional geological fault zone system can stabilize the convective flow of the system. (2) The double diffusion-driven convective flow modes of the three-dimensional geological fault zone system are very close each other and therefore, the system may have the similar chance to pick up different double diffusion-driven convective flow modes, especially in the case of the fault thickness to height ratio approaching 0. (3) The significant influence of the chemical species diffusion on the convective instability of the three-dimensional geological fault zone system implies that the seawater intrusion into the surface of the Earth is a potential mechanism to trigger the convective flow in the shallow three-dimensional geological fault zone system.

Numerical simulation of double-diffusion driven convective flow and rock alteration in three-dimensional fluid-saturated geological fault zones

Numerical methods are used to simulate the double-diffusion driven convective pore-fluid flow and rock alteration in three-dimensional fluid-saturated geological fault zones. The double diffusion is caused by a combination of both the positive upward temperature gradient and the positive downward salinity concentration gradient within a three-dimensional fluid-saturated geological fault zone, which is assumed to be more permeable than its surrounding rocks. In order to ensure the physical meaningfulness of the obtained numerical solutions, the numerical method used in this study is validated by a benchmark problem, for which the analytical solution to the critical Rayleigh number of the system is available. The theoretical value of the critical Rayleigh number of a three-dimensional fluid-saturated geological fault zone system can be used to judge whether or not the double-diffusion driven convective pore-fluid flow can take place within the system. After the possibility of triggering the double-diffusion driven convective pore-fluid flow is theoretically validated for the numerical model of a three-dimensional fluid-saturated geological fault zone system, the corresponding numerical solutions for the convective flow and temperature are directly coupled with a geochemical system. Through the numerical simulation of the coupled system between the convective fluid flow, heat transfer, mass transport and chemical reactions, we have investigated the effect of the double-diffusion driven convective pore-fluid flow on the rock alteration, which is the direct consequence of mineral redistribution due to its dissolution, transportation and precipitation, within the three-dimensional fluid-saturated geological fault zone system. (c) 2005 Elsevier B.V. All rights reserved.

Establishment of native ecosystems after mining - Australian experience across diverse biogeographic zones

The Australian minerals industry, which is dominated by coal, gold, bauxite, iron ore, base metals and mineral sand operations, is widely scattered across a continent which has a wide range of climatic zones ranging from moist temperate in the south through hot deserts in the centre to moist tropical in the north. There is an emphasis at most mines on establishing native ecosystems after mining, and technologies have had to be developed to ensure successful establishment and stability of these ecosystems under often adverse climatic conditions. This paper describes some of the innovative practices used to establish native ecosystenms in bauxite, mineral sand and coal operations across diverse biogeographic zones. Additionally, brief reference is made to an ecosystem function analysis, which has been developed to assess the success of establishment of these ecosystems. (C) 2001 Elsevier Science B.V. All rights reserved.

Delineation of sulphide ore-zones by borehole radar tomography at Hellyer Mine, Australia

Velocity and absorption tomograms are the two most common forms of presentation of radar tomographic data. However, mining personnel, geophysicists included, are often unfamiliar with radar velocity and absorption. In this paper, general formulae are introduced, relating velocity and attenuation coefficient to conductivity and dielectric constant. The formulae are valid for lossy media as well as high-resistivity materials. The transformation of velocity and absorption to conductivity and dielectric constant is illustrated via application of the formulae to radar tomograms from the Hellyer zinc-lead-silver mine, Tasmania, Australia. The resulting conductivity and dielectric constant tomograms constructed at Hellyer demonstrated the potential of radar tomography to delineate sulphide ore zones. (C) 2001 Elsevier Science B.V. All rights reserved.

Location of innervation zones of sternocleidomastoid and scalene muscles - a basis for clinical and research electromyography applications

Objectives: Advances in surface electromyography (sEMG) techniques provide a clear indication that refinement of electrode location relative to innervation zones (IZ) is required in order to optimise the accuracy, relevance and repeatability of the sEMG signals. The aim of this study was to identify the IZ for the sternocleidomastoid and anterior scalene muscles to provide guidelines for electrode positioning for future clinical and research applications. Methods: Eleven volunteer subjects participated in this study. Myoelectric signals were detected from the sternal and clavicular heads of the stemocleidomastoid and the anterior scalene muscles bilaterally using a linear array of 8 electrodes during isometric cervical flexion contractions. The signals were reviewed and the IZ(s) were identified, marked on the subjects' skin and measurements were obtained relative to selected anatomical landmarks. Results: The position of the IZ lay consistently around the mid-point or in the superior portion of the muscles studied. Conclusions: Results suggest that electrodes should be positioned over the lower portion of the muscle and not the mid-point, which has been commonly used in previous studies. Recommendations for sensor placement on these muscles should assist investigators and clinicians to ensure improved validity in future sEMG applications. (C) 2002 Elsevier Science Ireland Ltd. All rights reserved.