Elasto-plastic analysis of jointed rocks using discrete continuum and equivalent continuum approaches

Results from elasto-plastic numerical simulations of jointed rocks using both the equivalent continuum and discrete continuum approaches are presented, and are compared with experimental measurements. Initially triaxial compression tests on different types of rocks with wide variation in the uniaxial compressive strength are simulated using both the approaches and the results are compared. The applicability and relative merits and limitations of both the approaches for the simulation of jointed rocks are discussed. It is observed that both the approaches are reasonably good in predicting the real response. However, the equivalent continuum approach has predicted somewhat higher stiffness values at low strains. Considering the modelling effort involved in case of discrete continuum approach, for problems with complex geometry, it is suggested that a proper equivalent continuum model can be used, without compromising much on the accuracy of the results. Then the numerical analysis of a tunnel in Japan is taken up using the continuum approach. The deformations predicted are compared well against the field measurements and the predictions from discontinuum analysis. (C) 2012 Elsevier Ltd. All rights reserved.

Ultrahigh-temperature metamorphism from an unusual corundum plus orthopyroxene intergrowth bearing Al-Mg granulite from the Southern Marginal Zone, Limpopo Complex, South Africa

The Southern Marginal Zone of the Limpopo Complex is composed of granite-greenstone cratonic rocks reworked by a Neoarchean high-grade tectono-metamorphic event. Petrographic and mineral chemical characterization of an Al-Mg granulite from this zone is presented here. The granulite has a gneissic fabric with distinct Al-rich and Si-rich layers, with the former preserving the unusual lamellar (random and regular subparallel) intergrowths of corundum and symplectic intergrowth of spinel with orthopyroxene. The Al-rich layer preserves mineral assemblages such as rutile with orthopyroxene + sillimanite +/- A quartz, Al-rich orthopyroxene (similar to 11 wt%), spinel + quartz, and corundum in possible equilibrium with quartz, while the Si-rich layer preserves antiperthites and orthopyroxene + sillimanite +/- A quartz, all considered diagnostic of ultrahigh-temperature metamorphism. Application of Al-in-opx thermometry, ternary feldspar thermometry and construction of suitable pressure-temperature phase diagrams, compositional and model proportion isopleth results indicate P-T conditions as high as similar to 1,050-1,100 A degrees C, and similar to 10-12 kbars for the Al-Mg granulite. Our report of ultrahigh-temperature conditions is significant considering that the very high temperature was reached during decompression of an otherwise high-pressure granulite complex (clockwise P-T path), whereas most other ultrahigh-temperature granulites are linked to magma underplating at the base of the crust (counterclockwise P-T path).

High-temperature, high-pressure granulites (retrogressed eclogites) in the central region of the Lewisian, NW Scotland: Crustal-scale subduction in the Neoarchaean

Eclogites and associated high-pressure (HP) rocks in collisional and accretionary orogenic belts preserve a record of subduction and exhumation, and provide a key constraint on the tectonic evolution of the continents. Most eclogites that formed at high pressures but low temperatures at > 10-11 kbar and 450-650 degrees C can be interpreted as a result of subduction of cold oceanic lithosphere. A new class of high-temperature (HT) eclogites that formed above 900 degrees C and at 14 to 30 kbar occurs in the deep continental crust, but their geodynamic significance and processes of formation are poorly understood. Here we show that Neoarchaean mafic-ultramafic complexes in the central granulite facies region of the Lewisian in NW Scotland contain HP/HT garnet-bearing granulites (retrogressed eclogites), gabbros, Iherzolites, and websterites, and that the HP granulites have garnets that contain inclusions of omphacite. From thermodynamic modeling and compositional isopleths we calculate that peak eclogite-facies metamorphism took place at 24-22 kbar and 1060-1040 degrees C. The geochemical signature of one (G-21) of the samples shows a strong depletion of Eu indicating magma fractionation at a crustal level. The Sm-Nd isochron ages of HP phases record different cooling ages of ca. 2480 and 2330 Ma. We suggest that the layered mafic-ultramafic complexes, which may have formed in an oceanic environment, were subducted to eclogite depths, and exhumed as HP garnet-bearing orogenic peridotites. The layered complexes were engulfed by widespread orthogneisses of tonalite-trondhjemite-granodiorite (TTG) composition with granulite facies assemblages. We propose two possible tectonic models: (1) the fact that the relicts of eclogitic complexes are so widespread in the Scourian can be taken as evidence that a >90 km x 40 km-size slab of continental crust containing mafic-ultramafic complexes was subducted to at least 70 km depth in the late Archaean. During exhumation the gneiss protoliths were retrogressed to granulite facies assemblages, but the mafic-ultramafic rocks resisted retrogression. (2) The layered complexes of mafic and ultramafic rocks were subducted to eclogite-facies depths and during exhumation under crustal conditions they were intruded by the orthogneiss protoliths (TTG) that were metamorphosed in the granulite facies. Apart from poorly defined UHP metamorphic rocks in Norway, the retrogressed eclogites in the central granulite/retrogressed eclogite facies Lewisian region, NW Scotland have the highest crustal pressures so far reported for Archaean rocks, and demonstrate that lithospheric subduction was transporting crustal rocks to HP depths in the Neoarchaean. (C) 2012 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.

Numerical simulations of CO2 migration during charnockite genesis

Charnockite is considered to be generated either through the dehydration of granitic magma by CO2 purging or by solid-state dehydration through CO2 metasomatism during granulite facies metamorphism. To understand the extent of dehydration, CO2 migration is quantitatively modeled in silicate melt and metasomatic fluid as a function of temperature, H2O wt%, pressure, basal CO2 flux and dynamic viscosity. Numerical simulations show that CO2 advection through porous and permeable high-grade metamorphic rocks can generate dehydrated patches close to the CO2 flow path, as illustrated by the occurrences of ``incipient charnockites.'' CO2 reaction-front velocity constrained by field observations is 0.69 km/m.y., a reasonable value, which matches well with other studies. On the other hand, temperature, rate of cooling, and basal CO2 flux are the critical parameters affecting CO2 diffusion through a silicate melt. CO2 diffusion through silicate melt can only occur at temperature greater than 840 degrees C and during slow cooling (<= 3.7 x 10(-5) degrees C/yr), features that are typical of magma emplacement in the lower crust. Stalling of CO2 fluxing at similar to 840 degrees C explains why some deep-level plutons contain both hydrous and anhydrous (charnockitic) mineral assemblages. CO2 diffusion through silicate melt is virtually insensitive to pressure. Addition of CO2 basal flux facilitates episodic dehydrated melt migration by generating fracture pathways.

Chromium isotope variations (delta Cr-53/52) in mantle-derived sources and their weathering products: Implications for environmental studies and the evolution of delta Cr-53/52 in the Earth's mantle over geologic time

Here we report chromium isotope compositions, expressed as delta Cr-53/ 52 in per mil (&) relative to NIST 979, measured in selected Cr-rich minerals and rocks formed by the primary magmatic as well as the secondary metamorphic and weathering processes. The main objectives of this study were: (i) to further constrain the isotope composition of the Earth's mantle Cr inventory and its possible variation during geological history, based on the analysis of globally distributed and stratigraphically constrained mantle-derived chromites; and (ii) to investigate the magnitude and systematics of Cr isotope fractionation during oxidative weathering and secondary alteration (i. e., hydration, serpentinization) of the magmatic Cr sources. Specifically, we analyzed delta Cr-53/ 52 in a set of globally distributed mantle-derived chromites (FeMgCr2O4, n = 30) collected from various locations in Europe, Asia, Africa and South America, and our results confirm that a chromite-hosted Earth's mantle Cr inventory is uniform at - 0.079 +/- 0.129& (2SD), which we named here as a ` canonical' mantle d 53/ 52 Cr signature. Furthermore our dataset of stratigraphically constrained chromites, whose crystallization ages cover most of the Earth's geological history, indicate that the bulk Cr isotope composition of the chromite-hosted mantle inventory has remained uniform, within about +/- 0.100&, since at least the Early Archean times (similar to 3500 million years ago, Ma). To investigate the systematics of Cr isotope fractionation associated with alteration processes we analyzed a number of secondary Cr-rich minerals and variably altered ultramafic rocks (i. e., serpentinized harzburgites, lherzolites) that revealed large positive delta Cr-53/ 52 anomalies that are systematically shifted to higher values with an increasing degree of alteration and serpentinization. The degree of aqueous alteration and serpentinization was quantified by the abundances of fluid-mobile (Rb, K) elements, and by the Loss On Ignition (LOI) parameter, which determines the amount of structurally bound water (OH/ H2O) present in secondary hydrated minerals like serpentine. Overall, we observed that altered ultramafic rocks that yielded the highest LOI values, and the lowest amounts of fluid mobile elements, also yielded the heaviest delta Cr-53/ 52 signatures. Therefore, we conclude that secondary alteration (i.e., hydration, serpentinization) of ultramafic rocks in near-surface oxidative environments tend to shift the bulk Cr isotope composition of the weathered products to isotopically heavier values, pointing to a dynamic redox cycling of Cr in the Earth's crustal and near-surface environments. Hence, if validated by future

A Rodinian suture in western India: New insights on India-Madagascar correlations

We report detailed evidence for a new paleo-suture zone (the Kumta suture) on the western margin of southern India. The c. 15-km-wide, westward dipping suture zone contains garnet-biotite, fuchsite-haematite, chlorite-quartz, quartz-phengite schists, biotite augen gneiss, marble and amphibolite. The isochemical phase diagram estimations and the high-Si phengite composition of quartz-phengite schist suggest a near-peak condition of c. 18 kbar at c. 550 degrees C, followed by near-isothermal decompression. The detrital SHRIMP U-Pb zircon ages from quartz-phengite schist give four age populations ranging from 3280 to 2993 Ma. Phengite from quartz-phengite schist and biotite from garnet-biotite schist have K-Ar metamorphic ages of ca. 1326 and ca. 1385 Ma respectively. Electron microprobe-CHIME ages of in situ zircons in quartz-phengite schist (ca. 3750 Ma and ca. 1697 Ma) are consistent with the above results. The Bondla ultramafic-gabbro complex in the west of the Kumta suture compositionally represents an arc with K-Ar biotite ages from gabbro in the range 1644-1536 Ma. On the eastern side of the suture are weakly deformed and unmetamorphosed shallow westward-dipping sedimentary rocks of the Sirsi shelf, which has the following upward stratigraphy: pebbly quartzite/sandstone, turbidite, magnetite iron formation, and limestone; farther east the lower lying quartzite has an unconformable contact with ca. 2571 Ma quartzo-feldspathic gneisses of the Dharwar block with a ca. 1733 Ma biotite cooling age. To the west of the suture is a c. 60-km-wide Karwar block mainly consisting of tonalite-trondhjemite-granodiorite (TTG) and amphibolite. The TTGs have U-Pb zircon magmatic ages of ca. 3200 Ma with a rare inherited core age of ca. 3601 Ma. The K-Ar biotite cooling age from the TTGs (1746 Ma and 1796 Ma) and amphibolite (ca. 1697 Ma) represents late-stage uplift. Integration of geological, structural and geochronological data from western India and eastern Madagascar suggest diachronous ocean closure during the amalgamation of Rodinia; in the north at around ca. 1380 Ma, and a progression toward the south until ca. 750 Ma. Satellite imagery based regional structural lineaments suggests that the Betsimisaraka suture continues into western India as the Kumta suture and possibly farther south toward a suture in the Coorg area, representing in total a c. 1000 km long Rodinian suture. (C) 2013 Elsevier B.V. All rights reserved.

A Short Review on Lu-Hf Isotope System in Zircon: Implications for Crustal Evolution

Zircon has been recognized as the unaltered part of the Earth's history which preserves nearly 4 billion year record of earth's evolution. Zircon preserves igneous and metamorphic processes during its formation and remains unaffected by sedimentary processes and crustal recycling. U-Pb and Lu-Hf in zircon work as geochronometer and geochemical tracer respectively. Zircon provide valuable information about the source composition of the rocks and the intrinsic details of an unseen crust-mantle processes. The world wide data of U-Pb and Lu-Hf isotope systems in zircon reveal crustal evolution through geological history. Moreover, the U-Pb age pattern of zircons show distinct peaks attributed to preservation of crustal rocks or mountain building during supercontinent assembly. The histogram of continental crust preservation shows that nearly one-third of continental crust was formed during the Archean, almost 20% was formed during Paleoproterozoic and 14% in last 400 Ma.