Genesis of metasomatic sapphirine-corundum-spinel-bearing granulites in Sri Lanka

The goal of the present study is to understand the mechanism of mass transfer, the composition and the role of fluids during crustal metasomatism in high-temperature metamorphic terranes. A well constrained case study, a locality at Rupaha, Sri Lanka was selected. It is located in the Highland Complex of Sri Lanka, which represents a small, but important fragment of the super-continent Gondwana. Excellent exposures of ultramafic rocks, which are embedded in granulites, were found at 10 localities. These provide a unique background for understanding the metasomatic processes. The boundary between the ultramafic and the granulite rocks are lined with metasomatic reaction zones up to 50cm in width. Progressing from the ultramafics to the granulite host rock, three distinct zones with the following mineral assemblages can be distinguished: (1). phlogopite + spinel + sapphirine, (2). spinel + sapphirine and (3). corundum + biotite + plagioclase. In order to assess the P-T-t path, the peak metamorphism and the exhumation history were constrained using different thermobarometers, as well as a diffusion model of garnet zoning. A maximum temperature of 875 ± 20oC (Opx-Cpx thermometer) and at the peak pressure of 9.0 ± 0.1 kbar (Grt-Cpx-Pl-Qtz) was calculated for the silicic granulite. The ultramafic rocks recorded a peak temperature of 840 ± 70oC (Opx-Cpx thermometer) at 9 kbar. Coexisting spinel and sapphirine from the reaction zone yield a temperature of 820 ± 40oC. This is in agreement with the peak-temperatures recorded in the adjacent granulites and ultramafics rocks. The structural concordance of the ultramafic rocks with the siliceous granulite host rock further support the suggestion, that all units have experienced the same peak metamorphism. Diffusion modeling of retrograde zoning in garnets from mafic granulites suggests a three-step cooling history. A maximum cooling rate of 1oC/Ma is estimated during the initial stage of cooling, followed by a cooling rate of ~30oC/Ma. The outermost rims of garnet indicate a slightly slower cooling rate at about 10-15oC/Ma. The sequences of mineral zones, containing a variety of Al-rich, silica undersaturated minerals in the reaction zones separating the ultramafic rocks from the silica-rich rocks can be explained by a diffusion model. This involves the diffusion of Mg from ultramafic rocks across the layers, and K and Si diffuse in opposite direction. Chemical potential of Mg and Si generated continuous monotonic gradient, allowing steady state diffusional transport across the profile. The strong enrichment in Al, and the considerable loss of Si, during the formation of reaction bands can be inferred from isocon diagrams. Some Al was probably added to the reaction zones, while Si was lost. This is most likely due to fluids percolating parallel to the zones at the boundary of the rock units. This study has shown that not only pressure and temperature conditions but most importantly PH2O and the concentration of the chlorine and fluorine in aqueous fluids also control the mass transport in different geological environments.

Structural and strain analysis in the Late Paleozoic coastal accretionary wedge of central Chile

Abstract In this study structural and finite strain data are used to explore the tectonic evolution and the exhumation history of the Chilean accretionary wedge. The Chilean accretionary wedge is part of a Late Paleozoic subduction complex that developed during subduction of the Pacific plate underneath South America. The wedge is commonly subdivided into a structurally lower Western Series and an upper Eastern Series. This study shows the progressive development of structures and finite strain from the least deformed rocks in the eastern part of the Eastern Series of the accretionary wedge to higher grade schist of the Western Series at the Pacific coast. Furthermore, this study reports finite-strain data to quantify the contribution of vertical ductile shortening to exhumation. Vertical ductile shortening is, together with erosion and normal faulting, a process that can aid the exhumation of high-pressure rocks. In the east, structures are characterized by upright chevron folds of sedimentary layering which are associated with a penetrative axial-plane foliation, S1. As the F1 folds became slightly overturned to the west, S1 was folded about recumbent open F2 folds and an S2 axial-plane foliation developed. Near the contact between the Western and Eastern Series S2 represents a prominent subhorizontal transposition foliation. Towards the structural deepest units in the west the transposition foliation became progressively flat lying. Finite-strain data as obtained by Rf/Phi and PDS analysis in metagreywacke and X-ray texture goniometry in phyllosilicate-rich rocks show a smooth and gradual increase in strain magnitude from east to west. There are no evidences for normal faulting or significant structural breaks across the contact of Eastern and Western Series. The progressive structural and strain evolution between both series can be interpreted to reflect a continuous change in the mode of accretion in the subduction wedge. Before ~320-290 Ma the rocks of the Eastern Series were frontally accreted to the Andean margin. Frontal accretion caused horizontal shortening and upright folds and axial-plane foliations developed. At ~320-290 Ma the mode of accretion changed and the rocks of the Western Series were underplated below the Andean margin. This basal accretion caused a major change in the flow field within the wedge and gave rise to vertical shortening and the development of the penetrative subhorizontal transposition foliation. To estimate the amount that vertical ductile shortening contributed to the exhumation of both units finite strain is measured. The tensor average of absolute finite strain yield Sx=1.24, Sy=0.82 and Sz=0.57 implying an average vertical shortening of ca. 43%, which was compensated by volume loss. The finite strain data of the PDS measurements allow to calculate an average volume loss of 41%. A mass balance approximates that most of the solved material stays in the wedge and is precipitated in quartz veins. The average of relative finite strain is Sx=1.65, Sy=0.89 and Sz=0.59 indicating greater vertical shortening in the structurally deeper units. A simple model which integrates velocity gradients along a vertical flow path with a steady-state wedge is used to estimate the contribution of deformation to ductile thinning of the overburden during exhumation. The results show that vertical ductile shortening contributed 15-20% to exhumation. As no large-scale normal faults have been mapped the remaining 80-85% of exhumation must be due to erosion.