Hello Guest User, You are visiting this website from a computer with an IP address of 172.69.59.95 with the name of '?' since Thu Apr 18, 2024 at 3:10:18 PM PT for approx. 0 minutes now.
The Sheahan Diamond Literature Reference Compilation - Scientific and Media Articles based on Major Keyword - Geothermobarometry
The Sheahan Diamond Literature Reference Compilation is compiled by Patricia Sheahan who publishes on a monthly basis a list of new scientific articles related to diamonds as well as media coverage and corporate announcements called the Sheahan Diamond Literature Service that is distributed as a free pdf to a list of followers. Pat has kindly agreed to allow her work to be made available as an online digital resource at Kaiser Research Online so that a broader community interested in diamonds and related geology can benefit. The references are for personal use information purposes only; when available a link is provided to an online location where the full article can be accessed or purchased directly. Reproduction of this compilation in part or in whole without permission from the Sheahan Diamond Literature Service is strictly prohibited. Return to Diamond Keyword Index
Sheahan Diamond Literature Reference Compilation - Scientific Articles by Author for all years
Each article reference in the SDLRC is tagged with one or more key words assigned by Pat Sheahan to highlight the main topics of the article. In an effort to make it easier for users to track down articles related to a specific topic, KRO has extracted these key words and developed a list of major key words presented in this Key Word Index to which individual key words used in the article reference have been assigned. In most of the individual Key Word Reports the references are in crhonological order, though in some such as Deposits the order is first by key word and then chronological. Only articles classified as "technical" (mainly scientific journal articles) and "media" (independent media articles) are included in the Key Word Index. References that were added in the most recent monthly update are highlighted in yellow.
Geothermobarometry is the science of measuring the pressure and temperature history of metamorphic and igneous rocks through examination of the chemical composition of individual minerals and the assemblages of different minerals. Geobarometry focuses on establishing the pressure conditions associated with mineral formation whereas Geothermometry focuses on the temperature conditions. The SDLRC also uses Thermobarometry and Thermometry as key words. This topic is relevant to diamonds because it enables the establishment of unique correlations between diamonds and other minerals based on pressure-temperature regimes.
The Effect of Analytical and Experimental Errors on Temperature and Pressure Estimates Based on Analyses of Pyroxenes From Garnet Lherzolite Xenoliths in Kimberlites.
An Ultrapotassic Basaltic Suite from the Central Sierra Nevada, California: a Study of the Mineralogy, Petrology, Geochemistry and Isotopic Composition.
Ph.d. Thesis, University California, Santa Barbara., 100P.
United States, California, West Coast
Basanite, Whole Rock Geochemistry, Isotope, Geothermometry
A major change in the thermal state of the earth at the Archean Proterozoic boundary: consequences for the nature and preser-vation of continental lithosphere
Journal of Petrology, Special Volume 1988- Oceanic and Continental, pp. 39-52
high pressure experimental calibration of the olivine ortho pyroxene spinel oxygen geobarometer-implications for the oxidation state of the upper mantle
Contributions to Mineralogy and Petrology, Vol. 107, No. 1, pp. 27-40
Heterogeneity in the thermal state of the lower crust and upper mantle beneath eastern Australia
Australian Society of Exploration Geophysicists and Geological Society of Australia, 8th. Exploration Conference in the Bulletin., Vol. 22, No. 2, June pp. 295-298
Mantle xenoliths from the Quaternary Pali-Aike volcanic field of southernmost South America: implications for the accretion of Phanerozoic continentallithosphere
Proceedings of Fifth International Kimberlite Conference held Araxa June 1991, Servico Geologico do Brasil (CPRM) Special, pp. 395-397
Peru, South America
Petrochemical, mineral chemistry, geothermometry, Basalts
Mineral chem. of silicate and oxide phases from fertile peridotite equilibrated with a C-O-H fluid phase- a low fO2 dat a set- evaluation of mineralbarometers, therM.
Proceedings of Fifth International Kimberlite Conference held Araxa June 1991, Servico Geologico do Brasil (CPRM) Special, pp. 417-419
Geodynamic evolution and thermal history of the central Flin Flon Domain, Trans Hudson Orogen: constraints from Structural development 40 Ar 39 Ar and stable isotope
Baker, J., Chazot, G., Menzies, M.A., Thirlwall, M.
Lithospheric mantle beneath Arabia: a Pan-African protolith modified by the Afar and older plumes, rather than a source for continental flood volcanism?
Geological Society of America Special Paper, No. 362, pp. 65-80.
New pieces to the Archean terrane jigsaw puzzle in the Nuuk region, southern West Greenland: steps in transforming a simple insight into a complex regional tecton thermal model.
Journal of the Geological Society, Vol. 162, 1, pp. 147-162.
Mass transport mechanism between the upper and lower mantle in numerical simulations of thermochemical mantle convection with multicomponent phase changes.
Earth and Planetary Science Letters, Vol. 230, 1-2, pp. 11-27.
Mantle potential temperature at Hawaii, Iceland, and the mid-Ocean Ridge system, as inferred from olivine phenocrysts: evidence thermally driven mantle plumes
Mantle potential temperature at Hawaii, Iceland, and the mid-Ocean Ridge system, as inferred from olivine phenocrysts: evidence thermally driven mantle plumes
The Trans Hudson Orogen of North America and the Himalayan Karakoram Tibetan Orogen of Asia: structural and thermal evolution of the lower and upper plates.
GAC Annual Meeting Halifax May 15-19, Abstract 1p.
Electrical signatures due to thermal anomalies along mobile belts reactivated by the trail and outburst of mantle plume: evidences from the Indian subcontinent.
Andreoli, M.A.G., Hart, R.J., Ashwal, L.D., Coetzee, H.
Correlations between U, Th content and metamorphic grade in the Western Namaqualand Belt, South Africa: with implications for radioactive heating of the crust.
Chalapathi Rao, N.V., Burgess, R., Anand, M., Mainkar, D.
Evidence for a Phanerozoic (478 Ma) Diamondiferous kimberlite emplacement epoch in the Indian Shield from 40 Ar/ 39Ar dating of the Kodomali kimberlite: implications ....
Geological Society of India, Bangalore November Meeting Group Discussion on Kimberlites and Related Rocks India, Abstract p. 103-106.
India, Bastar Craton, Rodinia
Tectonics - Kodomali, Pan African , Geothermometry
Multistage exhumation and juxaposition of lower continental crust in the western Canadian Shield: linking high resolution U Pb and 40 Ar / 39 Ar thermochronometry with pressure temperature deformation paths.
Internally consistent thermodynamic dat a set for dense hydrous magnesium silicates up to 35 GPa, 1600 degree C: implications for water circulation in deep mantle.
Physics of the Earth and Planetary Interiors, Vol. 156, 1-2, pp. 89-107.
Constraints on the coupled thermal evolution of the Earth's core and mantle, the age of the inner core and the origin of the 186 Os 188 Os core signal in plume..
Earth and Planetary Science Letters, In press - available
United States, Hawaii
Geothermometry - potassium, not specific to diamonds
Trans Hudson Orogen of North America and Himalaya Karakoram Tibetan Orogen of Asia: structural and thermal characteristics of the lower and upper plates.
A new method to simulate convection with strongly temperature and pressure dependent viscosity in a spherical shell: applications to the Earth's mantle.
Physics of the Earth and Planetary Interiors, in press available
Electrical signatures due to thermal anomalies along mobile belts reactivated by the trail and outburst of mantle plume: evidence from the Indian subcontinent.
Journal of Applied Geophysics, Vol. 58, 4, April, pp. 313-320.
Hu, S., Raza, A., Min, K., Kohn, B.P., Reiners, Ketcham, Wang, Gleadow
Late Mesozoic and Cenozoic thermotectonic evolution along a transect from the north Chin a craton through the Qinling orogen into the Yangtze craton, central.
Density profiles of oceanic slabs and surrounding mantle: integrated thermodynamic and thermal modeling, and implications for the fate of slabs at the 660 km
Physics of the Earth and Planetary Interiors, Vol. 172, 3-4, pp. 257-267.
Keepler, H., Dubrovinsky, L.S., Narygina, O., Kantor, I.
Optical absorption and radioactive thermal conductivity silicate perovskite to 125 Gpa at high pressures, silicate perovskite, abundant in Earth's mantle....
Science, Vol. 322, 5907 Dec. 5, pp. 1529-1531.
Mantle
Geothermometry Radioactive heat important in deep Earth
Van Summeren, J.R., Vandenberg, A.P., Van der Hilst, R.D.
Upwellings from a deep mantle reservoir filtered at the 660 km phase transition in thermochemical convection models and implications for intra-plate volcanism.
Physics of the Earth and Planetary Interiors, Vol. 172, 3-4, pp. 210-224.
Geothermobarometry for ultramafic assemblages from the Emeishan large igneous province, southwest Chin a and the Nikos and Zulu kimberlites, Nunavut, Canada.
GAC/MAC/AGU Meeting held May 23-27 Toronto, Abstract only
Zozulya, D.R., O'Brien, H., Peltonen, P., Lehtonen, M.
Thermobarometry of mantle derived garnets and pyroxenes of Kola region ( NW Russia): lithosphere composition, thermal regime and diamond prospectivity.
Bulletin of the Geological Society of Finland, Vol. 81, pp. 143-158.
Zozulya, D.R., O'Brien, H., Peltonen, P., Lehtonen, M.
Thermobarometry of mantle derived garnets and pyroxenes of Kola region ( NW Russia): lithosphere composition, thermal regime and diamond prospectivity.
Bulletin of the Geological Society of Finland, Vol. 81, pp. 143-158.
Strain heating as a mechanism for partial melting and ultrahigh temperature metamorphism in convergent orogens: implications of temperature dependent thermal
Journal of Geophysical Research, Vol. 115, B 12 B12417
Schmadicke, E., Okrusch, M., Rupprecht-Gutpwski, P., Will, T.M.
Garnet pyroxenite, eclogite and alkremite xenoliths from the off-craton Gibeon kimberlite field, Namibia: a window into the upper mantle of Rehoboth Terrane.
Time dependent convection models of mantle thermal structure constrained by seismic tomography and geodynamics: implications for mantle plume dynamics and CMB heat flow.
Geophysical Journal International, Vol. 190, 2, pp. 785-815.
An experimental study of Fe oxidation states in garnet and clinopyroxene as a function of temperature in the system CaO FeO Fe2O3 MgO Al2O3 SiO2: implications for garnet-clinopyroxene geothermometry.
Contributions to Mineralogy and Petrology, Vol. 164, 4, pp. 623-639.
Geological mapping of the Francistown area in northeastern Botswana by surface temperature and spectral emissivity information derived from advanced spaceborn thermal emission and reflection radiometer (ASTER) thermal infrared data.
Multiple reaction oxygen barometry for mantle peridotite: an internally consistent thermodynamic model for reactions and garnet solid-solutions, with applications to the oxidation state of lithospheric mantle.
Volcanic and Magmatic Studies Group meeting, Poster Held Jan. 6-8. See minsoc website
Heat flow dat a in the Four Corners area suggest Neogene crustal warming resulting from partial lithosphere replacement in the Colorado Plateau interior, southwest USA.
Geological Society of America Bulletin, Vol. 126, pp. 1084-1092.
Journal of African Earth Sciences, Vol. 100, pp. 20-41.
Africa, Namibia
Geothermometry
Abstract: Intracontinental deformation accommodated along major lithospheric scale shear zone systems and within associated extensional basins has been well documented within West, Central and East Africa during the Late Cretaceous. The nature of this deformation has been established by studies of the tectonic architecture of sedimentary basins preserved in this part of Africa. In southern Africa, where the post break-up history has been dominated by major erosion, little evidence for post-break-up tectonics has been preserved in the onshore geology. Here we present the results of 38 new apatite fission track analyses from the Damara region of northern Namibia and integrate these new data with our previous results that were focused on specific regions or sections only to comprehensively document the thermo-tectonic history of this region since continental break-up in the Early Cretaceous. The apatite fission track ages range from 449 ± 20 Ma to 59 ± 3 Ma, with mean confined track lengths between 14.61 ± 0.1 ?m (SD 0.95 ?m) to 10.83 ± 0.33 ?m (SD 2.84 ?m). The youngest ages (c. 80–60 Ma) yield the longest mean track lengths, and combined with their spatial distribution, indicate major cooling during the latest Cretaceous. A simple numerical thermal model is used to demonstrate that this cooling is consistent with the combined effects of heating caused by magmatic underplating, related to the Paraná-Etendeka continental flood volcanism associated with rifting and the opening of the South Atlantic, and enhanced erosion caused by major reactivation of major lithospheric structures within southern Africa during a key period of plate kinematic change that occurred in the South Atlantic and SW Indian ocean basins between 87 and 56 Ma. This phase of intraplate tectonism in northern Namibia, focused in discrete structurally defined zones, is coeval with similar phases elsewhere in Africa and suggests some form of trans-continental linkage between these lithospheric zones.
Contributions to Mineralogy and Petrology, Vol. 169, 21p.
Mantle
Geobarometry - eclogites
Abstract: Thermodynamic parameters have been calibrated for a geobarometer suitable for use on eclogitic mantle xenoliths. The barometer is based on the incorporation of tetrahedrally coordinated aluminum in clinopyroxene coexisting with garnet and has been calibrated using the results of piston cylinder and multi-anvil experiments performed between pressures of 3 and 7 GPa and temperatures from 1,200 to 1,550 °C. Starting materials were hydrous and anhydrous synthetic mixtures of basaltic bulk compositions that yielded homogeneous bimineralic garnet-clinopyroxene phase assemblages. The experimental data set was expanded by employing results from previous experimental studies conducted in eclogitic systems, which widened the range of applicable conditions and compositions. The calibration reproduces experimental pressures of bimineralic eclogite assemblages, in addition to SiO2-saturated and kyanite-bearing eclogites, to within 0.4 GPa at the 95 % confidence interval. The barometer was then used to examine equilibration pressures recorded by natural mantle eclogites from various xenolith locations covering a wide pressure, temperature, and compositional range.
Thesis, Doctoral Beyreuther Graduiertenschule fur Mathematik und Naturwissenschaften IN: ENGLISH, 222p. Available pdf
Mantle
Geobarometry
Abstract: Eclogite rocks, composed mainly of garnet and clinopyroxene, form principally as a metamorphic product of oceanic crust as it undergoes subduction. The equilibrium between garnet and clinopyroxene is sensitive to temperature and pressure, therefore eclogitic outcrops and xenoliths can reveal important information on the chemical and mineralogical processes occurring during such episodes. This is particularly the case for lithospheric eclogitic xenoliths from Archean cratons, which can potential reveal information on ancient (> 1 Ga) subduction events. To obtain information on the depth of origin of eclogitic xenoliths the first project of this thesis was designed to establish a thermodynamically grounded geobarometer, which is based on the incorporation of tetrahedrally coordinated aluminum in clinopyroxene coexisting with garnet: 2/3 Ca3Al2Si3O12 + 1/3 Mg3Al2Si3O12 = CaAl2SiO6 + CaMgSi2O6 Grossular Pyrope CaTs Diopside The reaction was calibrated against high-pressure and high-temperature experiments carried out in the multi-anvil and piston-cylinder apparatus between pressures of 3 and 7 GPa and temperatures from 1200 to 1550 °C. Starting materials were hydrous and anhydrous synthetic mixtures of basaltic bulk compositions that yielded homogeneous bimineralic garnet-clinopyroxene phase assemblages. The experimental data set was expanded by employing results from previous experimental studies conducted in eclogitic systems, which widened the range of applicable conditions and compositions. The calibration reproduces experimental pressures of bimineralic eclogite assemblages, in addition to SiO2-saturated and kyanite-bearing eclogites, to within 0.4 GPa at the 95 % confidence interval. The barometer was then used to examine equilibration pressures recorded by natural mantle eclogites from various xenolith locations covering a wide pressure, temperature, and compositional range. The results seem to indicate that many eclogite xenoliths fall towards the hotter side of the range of geothermal temperatures displayed by peridotitic xenoliths from the same localities. The second project calibrates the composition of majoritic garnets, which contain excess silicon substituted onto the octahedrally coordinated garnet site, coexisting with clinopyroxene as a function of pressure, temperature and bulk composition. Majorite substitution for a given bulk composition increases with pressure, and its proportion can in principal be used as a geobarometer. Single majoritic garnet crystals are found as inclusions in diamonds, which are generally used to support a sublithospheric origin in the deeper upper mantle or transition zone. The chemical compositions of such inclusions indicate that they have formed from a number of different lithologies, including mafic, ultramafic and pyroxenitic. These inclusions give important insight into the environment where diamonds crystallize and the evolution of deep subducted crustal material. The empirical barometer studied here is based on the three major majoritic substitutions: 2Al3+ = Mg2+ + Si4+ (Mj), Mg2+ + Al3+ = Na1+ + Si4+ (NaSi), Mg2+ + Al3+ = Na1+ + Ti4+ (NaTi), and the secondary effect of chromium on the stability of Mj. The barometer was calibrated against experiments conducted in the multi-anvil apparatus between pressures of 6 and 16 GPa and temperatures of 1000 to 1400 °C. In order to expand the applicability to a wide range of compositions experiments were performed in three different mafic compositions and in one pyroxenitic composition. Moreover, existing experimental data in mafic and ultramafic systems, including mid-ocean ridge basalts, kimberlite, komatiite and peridotite bulk compositions were included in the calibration covering pressures from 6 to 20 GPa and temperatures from 900 to 2200 °C. Applying the geobarometer to natural majoritic diamond inclusions reveals clearly that eclogitic and pyroxenitic inclusions form dominantly at conditions near the top of the transition zone (300-400 km). Peridotitic inclusions, however, have formed generally at shallower depths ~200 km within the Earth’s upper mantle. This may reflect the differences expected for the oxygen fugacity between peridotitic and eclogitic/pyroxenitic rocks. If diamonds form through the reduction of carbonate or CO2 bearing melts then most peridotitic rocks would be already within the diamond stability field by depths of 200 km. At greater pressures carbon would remain immobilized as diamond and there is no mechanism by which new diamonds can form. Eclogitic rocks formed by the subduction of oceanic crust, however, should form an intrinsically more oxidized environment that remains within the carbonate stability field to much higher pressures. The diamond stability field would be eventually reached, however, due to either the effect of pressure on controlling Fe3+/Fe2+ equilibria or due to partial melting, which would preferentially remove ferric iron and lower the Fe3+/?Fe of the residue, on which the oxygen fugacity is mainly dependent. In fact deep partial melting as slabs heat up may be the mechanism by which both the local oxygen fugacity is lowered, carbon is mobilized and pyroxenite rocks are formed, with the latter occurring through reaction with the surrounding peridotite. The third project was focused on the elasticity of garnet solid solutions formed from eclogitic compositions. Garnet is an important mineral because it is a major phase in the upper mantle and dominates mafic rocks at these conditions. The elastic behavior of garnet solid solutions plays a role in the interpretation of seismic data and is important for estimating the density contrast between subducting slabs and the surrounding mantle. High-precision single-crystal X-ray diffraction measurements at high pressure have been conducted on three different ternary garnet solid solutions with varying chemical composition in order to examine possible non-ideal contributions to the volume and compressibility of garnet solid solutions. Furthermore one experiment has been conducted at high pressures and high temperatures to examine the effect of temperature on the elasticity of complex garnet solid solutions. Experimental results reveal that the concentration of the almandine (Fe3Al2Si3O12) component has a significant effect on the elasticity, whereby 10 - 20 mole% of almandine is sufficient to overprint the previously reported minima of the bulk modulus along the pyrope grossular join. It has also been shown that minor compositional variations of Ca and Mg within the Mg Fe Ca garnet ternary are not resolvable within the analytical errors. Therefore, the two eclogitic samples have similar bulk moduli within the analytical uncertainties. In contrast to previous studies, no evidence was found that garnets have a K’ significantly different from 4. The high-temperature experiment revealed that the relatively small fraction of almandine in a solid solution increased the softening of garnet with temperature. Finally, the experimental volumes and calculated densities have been compared to the self-consistent thermodynamic model of Stixrude and Lithgow-Bertelloni (2005, 2011). The comparison clearly reveals that volume and elastic properties cannot be linearly interpolated as a function of composition. Moreover, it has been shown that the excess properties vary not only as a function of composition and pressure, but also as a function of temperature. The final project is focused on the partitioning of fluorine (F) between garnet, clinopyroxene, and silicate melt within eclogitic compositions. Fluorine is the most abundant halogen on Earth and plays an important role in the formation of ultrapotassic lithologies, i.e. lamproites that contain several wt.% F, in contrast to the average lithospheric mantle that contains only tens of µg/g F. The cycling and partitioning behavior of fluorine in the Earth’s mantle are not well understood. High-pressure experiments have been conducted in the multi-anvil apparatus to obtain mineral-melt partition coefficients between garnet, clinopyroxene, and coexisting silicate melt of fluorine in a mafic system under conditions of the Earth’s upper mantle. The results show that mafic crust can host significantly more fluorine than the surrounding ultramafic mantle, due to the much higher proportion of clinopyroxene and its high fluorine partition coefficient of D_F^(clinopyroxene/melt)= 0.057 - 0.074. Combining the fluorine partitioning data with water partitioning data yields a plausible process to generate lamproitic magmas with a high F/H2O ratio. The enrichment of fluorine relative to H2O is triggered by multiple episodes of small degree melting which deplete the residual more in H2O than in fluorine, caused by the approximately three times smaller mineral-melt partition coefficients of H2O.
Geological Society of America Special Paper, No. 514, pp. SPE514-07.
Mantle
Geothermometry
Abstract: Calculations of mantle convection generally use constant rates of internal heating and time invariant core-mantle boundary temperature. In contrast parameterized convection calculations, sometimes called thermal history calculations, allow these properties to vary with time but only provide a single average temperature for the entire mantle. Here I consider 3D spherical convection calculations that run for the age of the Earth with heat producing elements that decrease with time, a cooling core boundary condition, and a mobile lid. The calculations begin with a moderately hot initial temperature, consistent with a relatively short accretion time for the formation of the planet. I find that the choice of a mobile or stagnant lid has the most significant effect on the average temperature as a function of time in the models. However the choice of mobile versus stagnant lid has less of an effect on the distribution of hot and cold anomalies within the mantle, or planform. I find the same low-degree (one upwelling or two upwelling) temperature structures in the mobile lid calculations that have previously been found in stagnant lid calculations. While having less of an effect on the mean mantle temperature, the viscosity of the asthenosphere has a profound effect on the pattern of temperature anomalies, even in the deep mantle. If the asthenosphere is weaker than the upper mantle by more than an order of magnitude, then the low-degree (one or two giant upwellings) pattern of temperature anomalies results. If the asthenosphere is less than an order of magnitude weaker than the upper mantle, then the pattern of temperature anomalies has narrow cylindrical upwellings and cold down going sheets. The low-degree pattern of temperature anomalies is more consistent with the plate model than the plume model (Foulger, 2007).
Geochemistry, Geophysics, Geosystems: G3, Vol. 16, 10, pp. 3400-3413.
Mantle
Geothermometry
Abstract: We investigate the influence of tectonic mode on the thermochemical evolution of simulated mantle convection coupled to a parameterized core cooling model. The tectonic mode is controlled by varying the friction coefficient for brittle behavior, producing the three tectonic modes: mobile lid (plate tectonics), stagnant lid, and episodic lid. The resulting compositional structure of the deep mantle is strongly dependent on tectonic mode, with episodic lid resulting in a thick layer of subducted basalt in the deep mantle, whereas mobile lid produces only isolated piles and stagnant lid no basaltic layering. The tectonic mode is established early on, with subduction initiating at around 60 Myr from the initial state in mobile and episodic cases, triggered by the arrival of plumes at the base of the lithosphere. Crustal production assists subduction initiation, increasing the critical friction coefficient. The tectonic mode has a strong effect on core evolution via its influence on deep mantle structure; episodic cases in which a thick layer of basalt builds up experience less core heat flow and cooling and a failed geodynamo. Thus, a continuous mobile-lid mode existing from early times matches Earth's mantle structure and core evolution better than an episodic mode characterized by large-scale flushing (overturn) events.
Geochimica et Cosmochimica Acta, in press available
Geothermometry
Abstract: The use of wet chemistry and X-ray absorption near edge structure (XANES) spectroscopy to determine the oxidation state of Fe in submarine glasses and olivine-hosted melt inclusions has provided important new insights into the global systematics of Fe3+/?Fe in mid-ocean ridge basalts (MORB) [1, 2]. Because MORB are aggregates of near-fractional partial melts formed by decompression melting of variably depleted peridotite, it is difficult to judge the extent to which they directly reflect the oxidation state of the oceanic upper mantle. To provide a theoretical framework within which to interpret Fe3+/?Fe in MORB, I have developed a model that describes the behavior of Fe3+/?Fe during spinel lherzolite partial melting in a system closed to oxygen. Modeling is carried out by calculating the Fe3+/?Fe of olivine using the point defect model of [3], and determining Fe3+/?Fe of the bulk peridotite from mineral-mineral partitioning. The inter-mineral Fe3+/Fe2+ exchange coefficients are derived from Mössbauer data on natural spinel peridotites, and are parameterized in terms of oxygen fugacity, temperature, and the Fe content of the olivine. The Fe3+/?Fe of the melt is determined by combining mass-balance with an equation relating the Fe3+/?Fe of the melt to the fugacity of oxygen [4]. Spinel lherzolite partial melting is modeled after [5]. Modeling results indicate that oxygen fugacity does not follow the fayalite-magnetite-quartz (FMQ) buffer during partial melting. For isobaric partial melting, the system becomes reduced relative to FMQ with increasing extent of melting. This results from an increase in the FMQ buffer with increasing temperature, whereas oxygen fugacity in the peridotite remains nearly constant. Conversely, during polybaric partial melting the oxidation state of the residual peridotite increases relative to FMQ. The effective partition coefficient for Fe3+is larger than previously thought, so that a redox couple with S is not required to explain its compatibility during partial melting.
Contributions to Mineralogy and Petrology, Vol. 171, pp. 51-
Mantle
Geobarometry
Abstract: The stability of the majorite component in garnet has been experimentally investigated at high pressure and high temperature, focusing on the effect of bulk composition and temperature. High-pressure experiments were performed in a multi-anvil apparatus, at pressures ranging from 6 to 14.5 GPa, and temperatures between 1400 and 1700 °C. Experiments were performed in a range of bulk compositions in the system SiO2-Al2O3-Cr2O3-CaO-MgO with varying Cr/(Cr + Al) ratios. The majorite content of garnet gradually increases with pressure, and the composition of the garnet, specifically the Cr/(Cr + Al) ratio, exerts a significant effect on the majorite substitution. We found no significant effect of temperature. We use the experimental results in combination with the literature data to derive two empirical geobarometers, which can be used to determine the equilibration pressure of natural majoritic garnets of peridotitic and eclogitic bulk compositions. The barometer for peridotitic compositions is
P=?77.1+27.6×Si+1.67×Cr
And the barometer for eclogitic compositions is
P=?29.6+11.8×Si+7.81×Na+4.49×Ca.
Contributions to Mineralogy and Petrology, Vol. 171, 7, 14p.
China
Geothermometry
Abstract: The large equilibrium Mg isotope fractionation between clinopyroxene and garnet observed in eclogites makes it a potential high-precision geothermometer, but calibration of this thermometer by natural samples is still limited. Here, we report Mg isotopic compositions of eclogite whole rocks as well as Mg and O isotopic compositions of clinopyroxene and garnet separates from 16 eclogites that formed at different temperatures from the Dabie orogen, China. The whole-rock ?26Mg values vary from ?1.20 to +0.10 ‰. Among them, 11 samples display limited ?26Mg variations from ?0.36 to ?0.17 ‰, similar to those of their protoliths. The mineral separates exhibit very different ?26Mg values, from ?0.39 to +0.39 ‰ for clinopyroxenes and from ?1.94 to ?0.81 ‰ for garnets. The clinopyroxene -garnet Mg isotope fractionation (?26Mgclinopyroxene -garnet = ?26Mgclinopyroxene -?26Mggarnet) varies from 1.05 to 2.15 ‰. The clinopyroxene -garnet O isotope fractionation (?18Oclinopyroxene -garnet = ?18Oclinopyroxene -?18Ogarnet) varies from ?1.01 to +0.98 ‰. Equilibrium Mg isotope fractionation between clinopyroxene and garnet in the investigated samples is selected based on both the ?26Mgclinopyroxene versus ?26Mggarnet plot and the state of O isotope equilibrium between clinopyroxene and garnet. The equilibrium ?26Mgclinopyroxene -garnet and corresponding temperature data obtained in this study, together with those available so far in literatures for natural eclogites, are used to calibrate the clinopyroxene -garnet Mg isotope thermometer. This yields a function of ?26Mgclinopyroxene -garnet = (0.99 ± 0.06) × 106/T 2, where T is temperature in Kelvin. The refined function not only provides the best empirically calibrated clinopyroxene -garnet Mg isotope thermometer for precise constraints of temperatures of clinopyroxene- and garnet-bearing rocks, but also has potential applications in high-temperature Mg isotope geochemistry.
Abstract: We present major and trace element data on coexisting garnet and clinopyroxene from experiments carried out between 1.3 and 10 GPa and 970 and 1400 °C. We demonstrate that the lattice strain model, which was developed for applications to mineral-melt partitioning, can be adapted to garnet-clinopyroxene partitioning. Using new and published experimental data we develop a geothermometer for coexisting garnet and clinopyroxene using the concentration of rare earth elements (REE). The thermometer, which is based on an extension of the lattice strain model, exploits the tendency of minerals at elevated temperatures to be less discriminating against cations that are too large or too small for lattice sites. The extent of discrimination against misfit cations is also related to the apparent elasticity of the lattice site on which substitution occurs, in this case the greater stiffness of the dodecahedral X-site in garnet compared with the eightfold M2-site in clinopyroxene. We demonstrate that the ratio of REE in clinopyroxene to that in coexisting garnet is particularly sensitive to temperature. We present a method whereby knowledge of the major and REE chemistry of garnet and clinopyroxene can be used to solve for the equilibrium temperature. The method is applicable to any scenario in which the two minerals are in equilibrium, both above and below the solidus, and where the mole fraction of grossular in garnet is less than 0.4. Our method, which can be widely applied to both peridotitic and eclogitic paragenesis with particular potential for diamond exploration studies, has the advantage over commonly used Fe-Mg exchange thermometers in having a higher closure temperature because of slow interdiffusion of REE. The uncertainty in the calculated temperatures, based on the experimental data set, is less than ±80 °C.
Abstract: New thermodynamic data for skiagite garnet (Fe3Fe23+Si3O12) are derived from experimental phase-equilibrium data that extend to 10 GPa and are applied to oxybarometry of mantle peridotites using a revised six-component garnet mixing model. Skiagite is more stable by 12 kJ mol-1 than in a previous calibration of the equilibrium 2 skiagite = 4 fayalite + ferrosilite + O2, and this leads to calculated oxygen fugacities that are higher (more oxidized) by around 1-1•5 logfO2 units. A new calculation method and computer program incorporates four independent oxybarometers (including 2 pyrope + 2 andradite + 2 ferrosilite = 2 grossular + 4 fayalite + 3 enstatite + O2) for use on natural peridotite samples to yield optimum logfO2 estimates by the method of least squares. These estimates should be more robust than those based on any single barometer and allow assessment of possible disequilibrium in assemblages. A new set of independent oxybarometers for spinel-bearing peridotites is also presented here, including a new reaction 2 magnetite + 3 enstatite = 3 fayalite + 3 forsterite + O2. These recalibrations combined with internally consistent PT determinations for published analyses of mantle peridotites with analysed Fe2O3 data for garnets, from both cratonic (Kaapvaal, Siberia and Slave) and circumcratonic (Baikal Rift) regions, provide revised estimates of oxidation state in the lithospheric mantle. Estimates of logfO2 for spinel assemblages are more reduced than those based on earlier calibrations, whereas garnet-bearing assemblages are more oxidized. Importantly, this lessens considerably the difference between garnet and spinel oxybarometry that was observed with previous published calibrations.
Abstract: Heat flow and heat production data complement seismic information and provide strong constraints on crustal composition, thickness and evolution. They have helped understand the nature of the Mohorovicic discontinuity and the variations in seismic velocities below the Moho. Notably, heat flow studies have delineated the vertical distribution of heat producing elements throughout the crust and in the upper most mantle lithosphere. Analysis of global data sets on heat flow and crustal thickness demonstrate that there is no correlation between these two variables. This is due to the large spatial variations in crustal composition and heat production that exist within a single geological province. For a given crustal thickness, the Moho temperature varies within a wide range (? 300 K) depending on surface heat flux and crustal heat production. Thus one cannot use generic models based on a “type” crustal column to calculate crustal geotherms. In stable regions, lower crustal temperatures depend on the amount and vertical distribution of heat producing elements in the crust. These temperatures determine the conditions of crustal stability and impose a limit on the maximum thickness of a stabilized crust.
Abstract: The relationship between plate- and plume-tectonics is considered in view of the growth and breakdown of supercontinents, active rifting, the formation of passive volcanic-type continental margins, and the origin of time-progressive volcanic chains on oceanic and continental plates. The mantle wind phenomenon is described, as well as its effect on plume morphology and anisotropy of the ambient mantle. The interaction of plumes and mid-ocean ridges is discussed. The principles and problems of plume activity analysis in subduction- and collision-related foldbelts are considered and illustrated with examples.
Abstract: Trapped-charge dating methods including luminescence and electron spin resonance dating have high potential as low temperature (< 100 °C) thermochronometers. Despite an early proof of concept almost 60 years ago, it is only in the past two decades that thermoluminescence (TL), electron-spin-resonance (ESR), and optically stimulated luminescence (OSL), have begun to gain momentum in geological thermochronometry and thermometry applications. Here we review the physics of trapped-charge dating, the studies that led to its development and its first applications for deriving palaeo-temperatures and/or continuous cooling histories. Analytical protocols, which enable the derivation of sample specific kinetic parameters over laboratory timescales, are also described. The key limitation of trapped-charge thermochronometry is signal saturation, which sets an upper limit of its application to < 1 Ma, thus restricting it to rapidly exhuming terrains (> 200 °C Ma? 1), or elevated-temperature underground settings (> 30 °C). Despite this limitation, trapped-charge thermochronometry comprises a diverse suite of versatile methods, and we explore potential future applications and research directions.
Earth and Planetary Science Letters, Vol. 461, pp. 30-39.
Mantle
Majorite, geobarometry
Abstract: Most diamonds form in the Earth's lithosphere but a small proportion contain Si-rich majoritic garnet inclusions that indicate formation in the deeper mantle. The compositions of syngenetic garnet inclusions can potential yield information on both the depth and mantle lithology in which the diamonds formed. Pressure dependent changes in garnet compositions have been calibrated using the results of experiments conducted in a multi-anvil apparatus at pressures between 6 and 16 GPa and temperatures of 1000 to 1400?°C. Using the results of these experiments a barometer was formulated based on an empirical parameterisation of the two major majoritic substitutions, referred to as majorite (Maj ; Al3+=Mg2++Si4+Al3+=Mg2++Si4+), and Na-majorite (Na-Maj ; Mg2++Al3+=Na++Si4+Mg2++Al3+=Na++Si4+). Moreover, previously published experimental garnet compositions from basaltic, kimberlite, komatiite and peridotite bulk compositions were included in the calibration, which consequently covers pressures from 6 to 20 GPa and temperatures from 900 to 2100?°C. Experimental pressures are reproduced over these conditions with a standard deviation of 0.86 GPa. The barometer is used to determine equilibration pressures of approximately 500 reported garnet inclusions in diamonds from a range of localities. As the majority of these inclusions are proposed to be syngenetic this allows a detailed picture of diamond formation depths and associated source rocks to be established using inclusion chemistry. Geographic differences in diamond source rocks are mapped within the sub-lithospheric mantle to over 500 km depth. Continuous diamond formation occurs over this depth range within lithologies with eclogitic affinities but also in lithologies that appear transitional between eclogitic and peridotitic bulk compositions, with an affinity to pyroxenites. The geographic differences between eclogitic and pyroxenitic diamond source rocks are rationalised in terms of diamond formation within downwelling and upwelling regimes respectively. Macroscopic diamond formation in rocks with pyroxenite compositions are likely facilitated in the deep mantle by higher average oxidation states and low mineral H2OH2O solubility compared to the surrounding mantle, which aid the mobility of C-O-H volatile species. The apparent lack of inclusions with a peridotite affinity may result from generally low oxygen fugacities in such lithologies, which reduces carbon mobility, and the lack of a suitable oxidising agent to allow diamonds to form from CH4. This glimpse of deep carbon cycle processes implies that heterogeneities in the carbon content, redox state and chemical composition of the mantle may be strongly coupled.
Earth and Planetary Science Letters, Vol. 465, pp. 59-69.
Mantle
Geothermometry
Abstract: Crustal foundering is an important mechanism in the differentiation and recycling of continental crust. Nevertheless, little is known about the dynamics of the lower crust, the temporal scale of foundering and its role in the dynamics of active margins and orogens. This particularly applies to active settings where the lower crust is typically still buried and direct access is not possible. Crustal xenoliths derived from mantle depth in the Pamir provide a unique exception to this. The rocks are well-preserved and comprise a diverse set of lithologies, many of which re-equilibrated at high-pressure conditions before being erupted in their ultrapotassic host lavas. In this study, we explore the petrological and chronological record of eclogite and felsic granulite xenoliths. We utilized accessory minerals - zircon, monazite and rutile - for coupled in-situ trace-element analysis and U-(Th-)Pb chronology by laser-ablation (split-stream) inductively coupled plasma mass spectrometry. Each integrated analysis was done on single mineral zones and was performed in-situ in thin section to maintain textural context and the ability to interpret the data in this framework. Rutile thermo-chronology exclusively reflects eruption (View the MathML source11.17±0.06Ma), which demonstrates the reliability of the U-Pb rutile thermo-chronometer and its ability to date magmatic processes. Conversely, zircon and monazite reveal a series of discrete age clusters between 55-11 Ma, with the youngest being identical to the age of eruption. Matching age populations between samples, despite a lack of overlapping ages for different chronometers within samples, exhibit the effectiveness of our multi-mineral approach. The REE systematics and age data for zircon and monazite, and Ti-in-zircon data together track the history of the rocks at a million-year resolution. The data reveal that the rocks resided at 30-40 km depth along a stable continental geotherm at 720-750?°C until 24-20 Ma, and were subsequently melted, densified, and buried to 80-90 km depth - 20 km deeper than the present-day Moho - at View the MathML source930±35°C. The material descended rapidly, accelerating from 0.9-1.7 mm?yr?1 to 4.7-5.8 mm?yr?1 within 10-12 Myr, and continued descending after reaching mantle depth at 14-13 Ma. The data reflect the foundering of differentiated deep-crustal fragments (2.9-3.5 g?cm?3) into a metasomatized and less dense mantle wedge. Through our new approach in constraining the burial history of rocks, we provided the first time-resolved record of this crustal-recycling process. Foundering introduced vestiges of old evolved crust into the mantle wedge over a relatively short period (c. 10 Myr). The recycling process could explain the variability in the degree of crustal contamination of mantle-derived magmatic rocks in the Pamir and neighboring Tibet during the Cenozoic without requiring a change in plate dynamics or source region.
Earth and Planetary Science Letters, Vol. 465, pp. 29-37.
Mantle
Geothermometry
Abstract: (Mg,?Fe)O ferropericlase (Fp) is one of the important minerals comprising Earth's lower mantle, and its thermal conductivity could be strongly influenced by the iron content and its spin state. We examined the lattice thermal conductivity of (Mg,?Fe)O Fp containing 19 mol% iron up to 111 GPa and 300 K by means of the pulsed light heating thermoreflectance technique in a diamond anvil cell. We confirmed a strong reduction in the lattice thermal conductivity of Fp due to iron substitution as reported in previous studies. Our results also show that iron spin crossover in Fp reduces its lattice thermal conductivity as well as its radiative conduction. We also measured the electrical conductivity of an identical Fp sample up to 140 GPa and 2730 K, and found that Fp remained an insulator throughout the experimental conditions, indicating the electronic thermal conduction in Fp is negligible. Because of the effects of strong iron impurity scattering and spin crossover, the total thermal conductivity of Fp at the core-mantle boundary conditions is much smaller than that of bridgmanite (Bdg). Our findings indicate that Bdg (and post-perovskite) is the best heat conductor in the Earth's lower mantle, and distribution of iron and its valence state among the lower mantle minerals are key factors to control the lower mantle thermal conductivity.
Abstract: Decompression of hot mantle rock upwelling beneath oceanic spreading centers causes it to exceed the melting point (solidus), producing magmas that ascend to form basaltic crust ~6 to 7 kilometers thick. The oceanic upper mantle contains ~50 to 200 micrograms per gram of water (H2O) dissolved in nominally anhydrous minerals, which -relative to its low concentration-has a disproportionate effect on the solidus that has not been quantified experimentally. Here, we present results from an experimental determination of the peridotite solidus containing known amounts of dissolved hydrogen. Our data reveal that the H2O-undersaturated peridotite solidus is hotter than previously thought. Reconciling geophysical observations of the melting regime beneath the East Pacific Rise with our experimental results requires that existing estimates for the oceanic upper mantle potential temperature be adjusted upward by about 60°C.
Geochemistry, Geophysics, Geosystems: G3, Vol. 18, pp. 872-888.
Mantle
geothermometry
Abstract: Chemical composition of mafic magmas is a critical indicator of physicochemical conditions, such as pressure, temperature, and fluid availability, accompanying melt production in the mantle and its evolution in the continental or oceanic lithosphere. Recovering this information has fundamental implications in constraining the thermal state of the mantle and the physics of mantle convection throughout the Earth's history. Here a statistical approach is applied to a geochemical database of about 22,000 samples from the mafic magma record. Potential temperatures (Tps) of the mantle derived from this database, assuming melting by adiabatic decompression and a Ti-dependent (Fe2O3/TiO2?=?0.5) or constant redox condition (Fe2+/?Fe?=?0.9 or 0.8) in the magmatic source, are thought to be representative of different thermal “horizons” (or thermal heterogeneities) in the ambient mantle, ranging in depth from a shallow sublithospheric mantle (Tp minima) to a lower thermal boundary layer (Tp maxima). The difference of temperature (?Tp) observed between Tp maxima and minima did not change significantly with time (?170°C). Conversely, a progressive but limited cooling of ?150°C is proposed since ?2.5 Gyr for the Earth's ambient mantle, which falls in the lower limit proposed by Herzberg et al. [2010] (?150-250°C hotter than today). Cooling of the ambient mantle after 2.5 Ga is preceded by a high-temperature plateau evolution and a transition from dominant plumes to a plate tectonics geodynamic regime, suggesting that subductions stabilized temperatures in the Archaean mantle that was in warming mode at that time.
Abstract: Detailed petrography, microstructure, and geochemistry of garnet pyroxenite xenoliths in Holocene basanite tuffs from maars at Lakes Bullenmerri and Gnotuk (western Victoria, southeastern Australia) have been used to track their igneous and metamorphic history, enabling the reconstruction of the thermal-tectonic evolution of the lithospheric mantle. The exsolution of orthopyroxene and garnet and rare spinel, plagioclase, and ilmenite from complex clinopyroxene megacrysts suggests that the xenoliths originally were clinopyroxene-dominant cumulates associated with minor garnet, orthopyroxene, or spinel. The compositions of exsolved phases and their host clinopyroxene were reintegrated using measured modal proportions to show that the primary clinopyroxene was enriched in Al2O3 (5.53-13.63 wt%) and crystallized at ~1300-1500 °C and 16-30 kbar. These cumulates then underwent extensive exsolution, recrystallization, and reaction during cooling, and finally equilibrated at ~950-1100 °C and 12-18 kbar before entrainment in the basanites. Rare earth element (REE) thermobarometry of garnets and coexisting clinopyroxenes preserves evidence of an intermediate stage (1032 °C and 21 kbar). These results imply that the protoliths of the garnet pyroxenite formed at a range of depths from ~50 to 100 km, and then during or shortly after cooling, they were tectonically emplaced to higher levels (~40-60 km; i.e., uplifted by at least 10-20 km) along the prevailing geotherm. This uplift may have been connected with lithosphere-scale faulting during the Paleozoic orogeny, or during Mesozoic-Cenozoic rifting of eastern Australia.
Abstract: The previously unstudied cooling and exhumation history of mid-crustal rocks exposed on southeastern Baffin Island (Canada) provides new insights into the post-orogenic evolution of the Paleoproterozoic Trans-Hudson Orogen (THO). New 40Ar/39Ar step-heat analyses of biotite, muscovite and phlogopite and core-to-rim intra-grain 40Ar/39Ar analyses of muscovite have a range of apparent ages compatible with slow regional cooling following peak metamorphism. Twenty-nine amphibolite- to granulite-facies rocks were dated using the 40Ar/39Ar step-heating laser (CO2) method. 40Ar/39Ar spot analyses were performed across muscovite grains from three samples using an ultraviolet (UV) laser to investigate intra-grain 40Ar/39Ar age variations. Step-heating apparent ages range from ca. 1788–1622 Ma for biotite, 1720–1630 Ma for phlogopite and 1729–1657 Ma for muscovite. UV spot 40Ar/39Ar analyses in the three muscovite grains range from ca. 1661–1640 Ma, 1675–1645 Ma and 1680–1652 Ma, with core-to-rim apparent age gradients of 20–30 Myr. Previous studies resolved peak metamorphism in this region to between ca. 1860 and 1820 Ma and identified late- to post-THO zircon and monazite populations at ca. 1800–1750 Ma. Numerical diffusion models for Ar in muscovite were conducted to test different Proterozoic cooling and exhumation scenarios. Comparisons with our 40Ar/39Ar ages attest to cooling rates of ~ 1–2 °C/Myr following peak metamorphism and ~ 1.5–2.5 °C/Myr after ca. 1740 Ma. Anomalously old apparent 40Ar/39Ar ages, in cases equivalent to U–Pb zircon rim and monazite ages, likely result from incorporation of excess Ar. The results suggest that mid-crustal rocks on southeastern Baffin Island remained hotter than ~ 420–450 °C for ~ 150–200 Myr after peak metamorphism, with subsequent slow cooling and denudation rates that are typical of Proterozoic orogens. The apparent absence of orogenic collapse implies that, despite high temperatures and estimated maximum crustal thicknesses comparable to those of large, hot orogens, the THO remained gravitationally stable during its terminal phase.
Abstract: Kimberlites are mantle-derived ultramafic rocks preserved in volcanic and sub-volcanic edifices and are the main primary source of diamonds. The temperatures of formation, transport, eruption and deposition remain poorly constrained despite their importance for understanding the petrological and thermodynamic properties of kimberlite magmas and styles of volcanic eruption. Here, we present measured values of Colour Alteration Indices (CAI) for conodonts recovered from 76 Paleozoic carbonate xenoliths found within 11 pipes from the Chidliak kimberlite field on Baffin Island, Nunavut, Canada. The dataset comprises the largest range of CAI values (1.5 to 8) and the highest CAI values reported to date for kimberlite-hosted xenoliths. Thermal models for cooling of the Chidliak kimberlite pipes and synchronous heating of conodont-bearing xenoliths indicate time windows of 10–20 000 h and, for these short time windows, the measured CAI values indicate heating of the xenoliths to temperatures of 225 to >925 ?C. We equate these temperatures with the minimum temperatures of the conduit-filling kimberlite deposit (i.e. emplacement temperature, TE). The majority of the xenoliths record CAI values of between 5 and 6.5 suggesting heating of xenoliths to temperatures of 460 ?C–735 ?C. The highest CAI values are consistent with being heated to 700 ?C–925 ?C and establish the minimum conditions for welding or formation of clastogenic kimberlite deposits. Lastly, we use TE variations within and between individual pipes, in conjunction with the geology of the conduit-filling deposits, to constrain the styles of explosive volcanic eruption.
Contrasting thermal structure, melt depletion and metasomatism of mantle lithosphere beneath two Proterozoic terranes west of the Kaapvaal craton, southern Africa.
Journal of Metamorphic Geology, Vol. 35, 6, pp. 631-661.
Mantle
geothermometry
Abstract: The Proterozoic belts that occur along the margins of the West Australian Craton, as well as those in intraplate settings, generally share similar geological histories that suggest a common plate-margin driver for orogeny. However, the thermal drivers for intraplate orogenesis are more poorly understood. The Mutherbukin Tectonic Event records a protracted period of Mesoproterozoic reworking of the Capricorn Orogen and offers significant insight into both the tectonic drivers and heat sources of long-lived intraplate orogens. Mineral assemblages and tectonic fabrics related to this event occur within a 50 km-wide fault-bound corridor in the central part of the Gascoyne Province in Western Australia. This zone preserves a crustal profile, with greenschist facies rocks in the north grading to upper amphibolite facies rocks in the south. The P–T–t evolution of 13 samples from 10 localities across the Mutherbukin Zone is investigated using phase equilibria modelling integrated with in situ U–Pb monazite and zircon geochronology. Garnet chemistry from selected samples is used to further refine the P–T history and shows that the dominant events recorded in this zone are prolonged D1 transpression between c. 1,320 and 1,270 Ma, followed by D2 transtension from c. 1,210 to 1,170 Ma. Peak metamorphic conditions in the mid-crust reached >650°C and 4.4–7 kbar at c. 1,210–1,200 Ma. Most samples record a single clockwise P–T evolution during this event, although some samples might have experienced multiple perturbations. The heat source for metamorphism was primarily conductive heating of radiogenic mid- and upper crust, derived from earlier crustal differentiation events. This crust was thickened during D1 transpression, although the thermal effects persisted longer than the deformation event. Peak metamorphism was terminated by D2 transtension at c. 1,210 Ma, with subsequent cooling driven by thinning of the radiogenic crust. The coincidence of a sedimentary basin acting as a thermal lid and a highly radiogenic mid-crustal batholith restricted to the Mutherbukin Zone accounts for reworking being confined to a discrete crustal corridor. Our results show that radiogenic regions in the shallow to mid crust can elevate the thermal gradient and localize deformation, causing the crust to be more responsive to far-field stresses. The Mutherbukin Tectonic Event in the Capricorn Orogen was synchronous with numerous Mesoproterozoic events around the West Australian Craton, suggesting that thick cratonic roots play an important role in propagating stresses generated at distant plate boundaries.
Abstract: The long-term cooling of Earth's mantle is recorded in the declining temperature and volume of its volcanic outpourings over time. However, analyses of 89-million-year-old lavas from Costa Rica suggest that extremely hot mantle still lurks below.
Journal of Geophysical Research: Solid Earth, Vol. 122, 8, pp. 5992-6013.
Mantle
geothermometry
Abstract: Cratons form the stable core roots of the continental crust. Despite long-term stability, cratons have failed in the past. Cratonic destruction (e.g., North Atlantic Craton) due to chemical rejuvenation at the base of the lithosphere remains poorly constrained numerically. We use 2-D petrological-thermomechanical models to assess cratonic rifting characteristics and mantle CO2 degassing in the presence of a carbonated subcontinental lithospheric mantle (SCLM). We test two tectonothermal SCLM compositions: Archon (depleted) and Tecton (fertilized) using 2 CO2 wt % in the bulk composition to represent a metasomatized SCLM. We parameterize cratonic breakup via extensional duration (7-12 Ma; full breakup), tectonothermal age, TMoho (300-600°C), and crustal rheology. The two compositions with metasomatized SCLMs share similar rifting features and decarbonation trends during initial extension. However, we show long-term (>67 Ma) stability differences due to lithospheric density contrasts between SCLM compositions. The Tecton model shows convective removal and thinning of the metasomatized SCLM during failed rifting. The Archon composition remained stable, highlighting the primary role for SCLM density even when metasomatized at its base. In the short-term, three failed rifting characteristics emerge: failed rifting without decarbonation, failed rifting with decarbonation, and semifailed rifting with dry asthenospheric melting and decarbonation. Decarbonation trends were greatest in the failed rifts, reaching peak fluxes of 94 × 104 kg m?3. Increased TMoho did not alter the effects of rifting or decarbonation. Lastly, we show mantle regions where decarbonation, mantle melting in the presence of carbonate, and preservation of carbonated mantle occur during rifting.
Abstract: Zircon- and bulk-rock Zr-based thermometric parameters have become fundamental to petrogenetic models of magmatism, from which broader geochronological and tectonic implications are being made. In particular, petrogenetic models have become increasingly reliant on Ti concentration in zircon geothermometry (TZircTi) and zircon saturation temperature (TZircsat). A feature of many of these studies is an implicit assumption that all zircons present in the host igneous rock are autocrystic, that is, crystallised from the surrounding melt. However, it has long been recognised that zircons present in an igneous rock can be inherited either from the surrounding country rock or source region (xenocrysts), or from earlier phases of magmatism or the magmatic plumbing system (antecrysts). Distinguishing these different origins for zircon crystals or domains within crystals is not straightforward. Here, we first review the utility and reliability of zircon-based thermometers for petrogenetic studies and show that TZircsat is a theoretical temperature and cannot be used to constrain magmatic or partial melting temperatures. It is a dynamic variable that changes during magma crystallisation, and essentially increases as fractional crystallisation proceeds, whereas true magmatic temperatures (TMagma) decrease. Generally, in Temperature-SiO2 space, the cross-over point of these two temperatures is magmatic system dependent, and also affected by the type of calibration used for the TZircsat calculations. Consequently, each magmatic system needs to be evaluated independently to assess the validity and usefulness of TZircsat. A fundamental conclusion of TZircsat and TMagma relationships assessed here is that new zircon generally only crystallises in silicic (granitic/rhyolitic) melt compositions, and thus autocrystic zircons should not be assumed to be present in igneous rocks with bulk compositions < 64 wt% SiO2, although inherited and minor zircons crystallising from late-stage differentiated melt pockets can be present. This highlights the importance of discriminating autocrystic from inherited zircons in igneous rocks. We then review techniques available to discriminate autocrystic from inherited zircons, and propose a new methodology to assist in the identification of autocrystic zircons for emplacement age determination and separate evaluation of inherited zircon components. The approach uses two strands of data: 1) zircon data such as zircon morphologies, textures, compositions and U-Pb ages, and 2) whole-rock data, in particular SiO2 and coupled geothermometry (TZircsat and TMagma) to estimate whether the magma was zircon-saturated or undersaturated. To test this new protocol, we use as examples, several Phanerozoic granitic rocks intersected by drilling in Queensland where contextual information is limited, and show how antecrystic and xenocrystic zircons and monazites can be distinguished. In contrast, where zircons are metamict (for example, high U and Th-rich zircons), much of the ability to discriminate is impacted because such zircons have suffered Pb loss and have modified compositions (e.g., higher TZircTi). We recommend an integrated approach incorporating whole-rock chemistry, independent geothermometric constraints, zircon composition, textures and ages obtained by routine cathodoluminescence and LA-ICP-MS or ion microprobe analysis to provide increased confidence for the discrimination of inherited zircons from autocrystic zircons and determination of the emplacement age.
Abstract: Elastic geothermobarometry is a method of determining metamorphic conditions from the excess pressures exhibited by mineral inclusions trapped inside host minerals. An exact solution to the problem of combining non-linear Equations of State (EoS) with the elastic relaxation problem for elastically isotropic spherical host-inclusion systems without any approximations of linear elasticity is presented. The solution is encoded into a Windows GUI program EosFit-Pinc. The program performs host-inclusion calculations for spherical inclusions in elastically isotropic systems with full P-V-T EoS for both phases, with a wide variety of EoS types. The EoS values of any minerals can be loaded into the program for calculations. EosFit-Pinc calculates the isomeke of possible entrapment conditions from the pressure of an inclusion measured when the host is at any external pressure and temperature (including room conditions), and it can calculate final inclusion pressures from known entrapment conditions. It also calculates isomekes and isochors of the two phases.
Abstract: The fragmentary Phanerozoic geological record of the anomalously elevated Zimbabwe Craton makes reconstructing its history difficult using conventional field methods. Here we constrain the cryptic Phanerozoic evolution of the Zimbabwe Craton using a spatially extensive apatite (U-Th-Sm)/He (AHe), apatite fission track (AFT), and zircon (U-Th)/He (ZHe) data set. Joint thermal history modeling reveals that the region experienced two cooling episodes inferred to be the denudational response to surface uplift. The first and most significant protracted denudation period was triggered by stress transmission from the adjacent ~750-500 Ma Pan-African orogenesis during the amalgamation of Gondwana. The spatial extent of this rejuvenation signature, encompassing the current broad topographic high, could indicate the possible longevity of an ancient topographic feature. The ZHe data reveal a second, minor denudation phase which began in the Paleogene and removed a kilometer-scale Karoo cover from the craton. Within our data set, the majority of ZHe ages are younger than their corresponding AHe and AFT ages, even at relatively low eU. This unexpectedly recurrent age “inversion” suggests that in certain environments, moderately, as well as extremely, damaged zircons have the potential to act as ultra-low-temperature thermochronometers. Thermal history modeling results reveal that the zircon radiation damage accumulation and annealing model (ZRDAAM) frequently overpredicts the ZHe age. However, the opposite is true for extremely damaged zircons where the ZHe and AHe data are also seemingly incompatible. This suggests that modification of the ZRDAAM may be required for moderate to extreme damage levels.
Abstract: Zircon- and bulk-rock Zr-based thermometric parameters have become fundamental to petrogenetic models of magmatism, from which broader geochronological and tectonic implications are being made. In particular, petrogenetic models have become increasingly reliant on Ti concentration in zircon geothermometry (TZircTi) and zircon saturation temperature (TZircsat). A feature of many of these studies is an implicit assumption that all zircons present in the host igneous rock are autocrystic, that is, crystallised from the surrounding melt. However, it has long been recognised that zircons present in an igneous rock can be inherited either from the surrounding country rock or source region (xenocrysts), or from earlier phases of magmatism or the magmatic plumbing system (antecrysts). Distinguishing these different origins for zircon crystals or domains within crystals is not straightforward. Here, we first review the utility and reliability of zircon-based thermometers for petrogenetic studies and show that TZircsat is a theoretical temperature and cannot be used to constrain magmatic or partial melting temperatures. It is a dynamic variable that changes during magma crystallisation, and essentially increases as fractional crystallisation proceeds, whereas true magmatic temperatures (TMagma) decrease. Generally, in Temperature-SiO2 space, the cross-over point of these two temperatures is magmatic system dependent, and also affected by the type of calibration used for the TZircsat calculations. Consequently, each magmatic system needs to be evaluated independently to assess the validity and usefulness of TZircsat. A fundamental conclusion of TZircsat and TMagma relationships assessed here is that new zircon generally only crystallises in silicic (granitic/rhyolitic) melt compositions, and thus autocrystic zircons should not be assumed to be present in igneous rocks with bulk compositions < 64 wt% SiO2, although inherited and minor zircons crystallising from late-stage differentiated melt pockets can be present. This highlights the importance of discriminating autocrystic from inherited zircons in igneous rocks. We then review techniques available to discriminate autocrystic from inherited zircons, and propose a new methodology to assist in the identification of autocrystic zircons for emplacement age determination and separate evaluation of inherited zircon components. The approach uses two strands of data: 1) zircon data such as zircon morphologies, textures, compositions and U-Pb ages, and 2) whole-rock data, in particular SiO2 and coupled geothermometry (TZircsat and TMagma) to estimate whether the magma was zircon-saturated or undersaturated. To test this new protocol, we use as examples, several Phanerozoic granitic rocks intersected by drilling in Queensland where contextual information is limited, and show how antecrystic and xenocrystic zircons and monazites can be distinguished. In contrast, where zircons are metamict (for example, high U and Th-rich zircons), much of the ability to discriminate is impacted because such zircons have suffered Pb loss and have modified compositions (e.g., higher TZircTi). We recommend an integrated approach incorporating whole-rock chemistry, independent geothermometric constraints, zircon composition, textures and ages obtained by routine cathodoluminescence and LA-ICP-MS or ion microprobe analysis to provide increased confidence for the discrimination of inherited zircons from autocrystic zircons and determination of the emplacement age.
Geochimica et Cosmochimica Acta, Vol. 222, Feb 1, pp. 421-435.
Mantle
geobarometry
Abstract: We have performed an experimental cross calibration of a suite of mineral equilibria within mantle rock bulk compositions that are commonly used in geobarometry to determine the equilibration depths of upper mantle assemblages. Multiple barometers were compared simultaneously in experimental runs, where the pressure was determined using in-situ measurements of the unit cell volumes of MgO, NaCl, Re and h-BN between 3.6 and 10.4?GPa, and 1250 and 1500?°C. The experiments were performed in a large volume press (LVPs) in combination with synchrotron X-ray diffraction. Noble metal capsules drilled with multiple sample chambers were loaded with a range of bulk compositions representative of peridotite, eclogite and pyroxenite lithologies. By this approach, we simultaneously calibrated the geobarometers applicable to different mantle lithologies under identical and well determined pressure and temperature conditions. We identified discrepancies between the calculated and experimental pressures for which we propose simple linear or constant correction factors to some of the previously published barometric equations. As a result, we establish internally-consistent cross-calibrations for a number of garnet-orthopyroxene, garnet-clinopyroxene, Ca-Tschermaks-in-clinopyroxene and majorite geobarometers.
Abstract: “Super-deep” diamonds are thought to have a sub-lithospheric origin (i.e., below ~300 km depth) because some of the mineral phases entrapped within them as inclusions are considered to be the products of retrograde transformation from lower-mantle or transition-zone precursors. CaSiO3-walstromite, the most abundant Ca-bearing mineral inclusion found in super-deep diamonds, is believed to derive from CaSiO3-perovskite, which is stable only below ~600 km depth, although its real depth of origin is controversial. The remnant pressure (Pinc) retained by an inclusion, combined with the thermoelastic parameters of the mineral inclusion and the diamond host, allows calculation of the entrapment pressure of the diamond-inclusion pair. Raman spectroscopy, together with X-ray diffraction, is the most commonly used method for measuring the Pinc without damaging the diamond host. In the present study we provide, for the first time, a calibration curve to determine the Pinc of a CaSiO3-walstromite inclusion by means of Raman spectroscopy without breaking the diamond. To do so, we performed high-pressure micro-Raman investigations on a CaSiO3-walstromite crystal under hydrostatic stress conditions within a diamond-anvil cell. We additionally calculated the Raman spectrum of CaSiO3-walstromite by ab initio methods both under hydrostatic and non-hydrostatic stress conditions to avoid misinterpretation of the results caused by the possible presence of deviatoric stresses causing anomalous shift of CaSiO3-walstromite Raman peaks. Last, we applied single-inclusion elastic barometry to estimate the minimum entrapment pressure of a CaSiO3-walstromite inclusion trapped in a natural diamond, which is ~9 GPa (~260 km) at 1800 K. These results suggest that the diamond investigated is certainly sub-lithospheric and endorse the hypothesis that the presence of CaSiO3-walstromite is a strong indication of super-deep origin.
Abstract: New versions of the universal Jd-Di exchange clinopyroxene barometer for peridotites, pyroxenites and eclogites, and also garnet barometer for eclogites and peridotites were developed. They were checked using large experimental data sets for eclogitic (?530) and peridotitic systems (>650). The precision of the universal Cpx barometer for peridotites based on Jd-Di exchange is close to Cr-Tschermakite method produced by Nimis and Taylor (2000). Cpx barometer was transformed by the substitution of major multiplier for KD by the equations dependent from Al-Na-Fe. Obtained equation in combination with the thermometer of Nimis and Taylor (2000) allow to reconstruct position of the magma feeder systems of the alkali basaltic magma within the mantle diapirs in modern platforms like in Vitim plateau and other Southern Siberia localities and several localities worldwide showing good agreement of pressure ranges for black and green suites. These equations allow construct PTX diagrams for the kimberlite localities in Siberia and worldwide calculating simultaneously the PT parameters for different groups of mantle rocks. They give very good results for the concentrates from kimberlite lamproites and placers with mantle minerals. They are useful for PT estimates for diamond inclusions. The positions of eclogite groups in mantle sections are similar to those determined with new Gar–Cpx barometer produced by C. Beyer et al. (2015). The Fe rich eclogites commonly trace the boundary between the lower upper parts of subcontinental lithospheric mantle (SCLM) at 3–4 GPa marking pyroxenite eclogites layer. Ca-rich eclogites and especially grospydites in SCLM beneath Precambrian kimberlites occurs near pyroxenite layer but in younger mantle sections they became common in the lower parts. The diamondiferous Mg Cr-less group eclogites referring to the ancient island arc complexes are also common in the middle part of mantle sections and near 5–6 GPa. Commonly eclogites in lower apart of mantle sections are remelted and trace the high temperature convective branch. The Mg- and Fe-rich pyroxenites also show the extending in pressure trends which suggest the anatexic melting under the influence of volatiles or under the interaction with plums.
Geochimica et Cosmochimica Acta, Vol. 225, pp. 1-16.
Technology
geobarometry
Abstract: Majoritic garnet [(Ca, Mg, Fe2+)3(Fe3+, Al, Si)2(SiO4)3] is one of the predominant and important constituents of upper mantle peridotite and ultra-deep subducted slabs. Majoritic substitution in garnet depends on pressure, and it has been used to estimate the formation pressure of natural majoritic garnet. Ferric iron (Fe3+) substitution occurs in natural majoritic garnets from mantle diamonds and shocked meteorites. However, available majorite geobarometers were developed without considering the effect of Fe3+ substitution in the structure. In this study, we systematically synthesized Fe3+- bearing majoritic garnets from 6.5?GPa to 15?GPa to evaluate the effect of Fe3+ on the majorite geobarometer. The Fe3+ contents of synthetic majoritic garnets were analyzed using the "Flank method" with the electron probe microanalyzer (EPMA). The results were compared with those based on the charge balance calculations. From the known synthesis pressures and measured Fe3+ contents, we developed a new majorite geobarometer for Fe3+-bearing majoritic garnets. Our results show that the existing majorite geobarometer, which does not take into account the Fe3+ substitution, could underestimate the formation pressure of majoritic garnets, especially for samples with a high majoritic component.
Abstract: Using low?temperature thermochronology on apatite and zircon crystals, we show that the western Reguibat Shield, located in the northern part of the West African Craton, experienced significant cooling and heating events between Jurassic and present times. The obtained apatite fission track ages range between 49 and 102 Ma with mean track lengths varying between 11.6 and 13.3 ?m and Dpar values between 1.69 and 3.08 ?m. Zircon fission track analysis yielded two ages of 159 and 118 Ma. Apatite (U-Th)/He uncorrected single?grain ages range between 76 and 95 Ma. Thermal inverse modelling indicates that the Reguibat Shield was exhumed during the Early Cretaceous, Late Cretaceous, Palaeocene-Eocene and Quaternary. These exhumation events were coeval with regional tectonic and geodynamic events, and were probably driven by a combined effect of plate tectonics and mantle dynamics.
Earth Planetary Science Letters, Vol. 491, pp. 148-159.
Africa, Australia, Canada, Europe
geothermometry
Abstract: The long-term stability of cratons has been attributed to low temperatures and depletion in iron and water, which decrease density and increase viscosity. However, steady-state thermal models based on heat flow and xenolith constraints systematically overpredict the seismic velocity-depth gradients in cratonic lithospheric mantle. Here we invert for the 1-D thermal structure and a depth distribution of metasomatic minerals that fit average Rayleigh-wave dispersion curves for the Archean Kaapvaal, Yilgarn and Slave cratons and the Proterozoic Baltic Shield below Finland. To match the seismic profiles, we need a significant amount of hydrous and/or carbonate minerals in the shallow lithospheric mantle, starting between the Moho and 70 km depth and extending down to at least 100-150 km. The metasomatic component can consist of 0.5-1 wt% water bound in amphibole, antigorite and chlorite, ?0.2 wt% water plus potassium to form phlogopite, or ?5 wt% CO2 plus Ca for carbonate, or a combination of these. Lithospheric temperatures that fit the seismic data are consistent with heat flow constraints, but most are lower than those inferred from xenolith geothermobarometry. The dispersion data require differences in Moho heat flux between individual cratons, and sublithospheric mantle temperatures that are 100-200?°C less beneath Yilgarn, Slave and Finland than beneath Kaapvaal. Significant upward-increasing metasomatism by water and CO2-rich fluids is not only a plausible mechanism to explain the average seismic structure of cratonic lithosphere but such metasomatism may also lead to the formation of mid-lithospheric discontinuities and would contribute to the positive chemical buoyancy of cratonic roots.
Physics of the Earth and Planetary Interiors, Vol. 277, 1, pp. 10-29.
Mantle
Geothermometry
Abstract: The dynamics of Earth’s lowermost mantle exert significant control on the formation of mantle plumes and the core-mantle boundary (CMB) heat flux. However, it is not clear if and how the variation of CMB heat flux and mantle plume activity are related. Here, we perform geodynamic model experiments that show how temporal variations in CMB heat flux and pulses of mantle plumes are related to morphologic changes of the thermochemical piles of large-scale compositional heterogeneities in Earth’s lowermost mantle, represented by the large low shear velocity provinces (LLSVPs). We find good correlation between the morphologic changes of the thermochemical piles and the time variation of CMB heat flux. The morphology of the thermochemical piles is significantly altered during the initiation and ascent of strong mantle plumes, and the changes in pile morphology cause variations in the local and the total CMB heat flux. Our modeling results indicate that plume-induced episodic variations of CMB heat flux link geomagnetic superchrons to pulses of surface volcanism, although the relative timing of these two phenomena remains problematic. We also find that the density distribution in thermochemical piles is heterogeneous, and that the piles are denser on average than the surrounding mantle when both thermal and chemical effects are included.
Mineralogy and Petrology, 10.1007/ s00710-018- 0577-8, 19p.
Africa
metasomatism, subduction, geobarometry
Abstract: The Kaapvaal craton and its underlying mantle is probably one of the best studied Archean entity in the world. Despite that, discussion is still vivid on important aspects. A major debate over the last few decades is the depth of melting that generated the mantle nuclei of cratons. Our new evaluation of melting parameters in peridotite residues shows that the Cr2O3/Al2O3 ratio is the most useful pressure sensitive melting barometer. It irrevocably constrains the pressure of melting (melt separation) to less than 2 GPa with olivine (ol), orthopyroxene (opx) and spinel (sp) as residual phases. Garnet (grt) grows at increasing pressure during lithosphere thickening and subduction via the reaction opx?+?sp ? grt?+?ol. The time of partial melting is constrained by Re-depletion model ages (TRD) mainly to the Archean (Pearson and Wittig 2008). However, only 3% of the ages are older than 3.1 Ga while crustal ages lie mainly between 3.1 to 2.8 Ga for the W- and 3.7 to 2.8 Ga for the E-block. Many TRD-ages are probably falsified by metasomatism and the main partial melting period was older than 3.1 Ga. Also, Nd- and Hf- model ages of peridotitic lithologies from the W-block are 3.2 to 3.6 Ga old. The corresponding very negative ?Nd (?40) and ?Hf values (?65) signal the presence of subducted crustal components in these old mantle portions. Subducted components diversify the mantle in its chemistry and thermal structure. Adjustment towards a stable configuration occurs by fluid transfer, metasomatism, partial melting and heat transfer. Ages of metasomatism from the Lu-Hf isotope system are 3.2 Ga (Lace), 2.9 Ga (Roberts Victor) and 2.62 Ga (Finsch) coinciding with the collision of cratonic blocks, the growth of diamonds, metamorphism of eclogites and of Ventersdoorp magmatism. The cratonic lithosphere was stabilized thermally by the end of the Archean and cooled since then with a rate of 0.07 °C/Ma.
Abstract: Some seismic models derived from tomographic studies indicate elevated shear?wave velocities (?4.7 km/s) around 120?150 km depth in cratonic lithospheric mantle. These velocities are higher than those of cratonic peridotites, even assuming a cold cratonic geotherm (i.e., 35 mW/m2 surface heat flux) and accounting for compositional heterogeneity in cratonic peridotite xenoliths and the effects of anelasticity. We reviewed various geophysical and petrologic constraints on the nature of cratonic roots (seismic velocities, lithology/mineralogy, electrical conductivity, and gravity) and explored a range of permissible rock and mineral assemblages that can explain the high seismic velocities. These constraints suggest that diamond and eclogite are the most likely high?Vs candidates to explain the observed velocities, but matching the high shear?wave velocities requires either a large proportion of eclogite (>50 vol.%) or the presence of up to 3 vol.% diamond, with the exact values depending on peridotite and eclogite compositions and the geotherm. Both of these estimates are higher than predicted by observations made on natural samples from kimberlites. However, a combination of ?20 vol.% eclogite and ~2 vol.% diamond may account for high shear?wave velocities, in proportions consistent with multiple geophysical observables, data from natural samples, and within mass balance constraints for global carbon. Our results further show that cratonic thermal structure need not be significantly cooler than determined from xenolith thermobarometry.
Annual Review of Earth and Planetary Sciences, Vol. 46, pp. 47-66.
Mantle
geothermometry
Abstract:
The thermal conductivity of iron alloys at high pressures and temperatures is a critical parameter in governing (a) the present-day heat flow out of Earth's core, (b) the inferred age of Earth's inner core, and (c) the thermal evolution of Earth's core and lowermost mantle. It is, however, one of the least well-constrained important geophysical parameters, with current estimates for end-member iron under core-mantle boundary conditions varying by about a factor of 6. Here, the current state of calculations, measurements, and inferences that constrain thermal conductivity at core conditions are reviewed. The applicability of the Wiedemann-Franz law, commonly used to convert electrical resistivity data to thermal conductivity data, is probed: Here, whether the constant of proportionality, the Lorenz number, is constant at extreme conditions is of vital importance. Electron-electron inelastic scattering and increases in Fermi-liquid-like behavior may cause uncertainties in thermal conductivities derived from both first-principles-associated calculations and electrical conductivity measurements. Additional uncertainties include the role of alloying constituents and local magnetic moments of iron in modulating the thermal conductivity. Thus, uncertainties in thermal conductivity remain pervasive, and hence a broad range of core heat flows and inner core ages appear to remain plausible.
Abstract: Empirical and experimental calibration of single element solubility thermometers, such as Zr-in-rutile, Zr-in-titanite, Ti-in-zircon, and Ti-in-quartz, within the past 13 years has greatly expanded our ability to assess the pressure and temperature conditions of individual minerals associated with specific textures in metamorphic rocks. Combined with advances in in situ techniques for analyzing trace concentrations, this has led to an increase in the combined use of single element thermometers, geochronometers, and isotope ratios, often simultaneously, in metamorphic minerals. Here we review the calibration and application of single element thermometers at the pressure and temperature conditions of interest in metamorphic rocks. We discuss to what extent accessory phase equilibrium and trace element equilibrium are attained in metamorphic systems, and the thermodynamic and kinetic framework within which trace element equilibrium is assessed. As an example, we present a comprehensive study of trace element distribution during rutile replacement by titanite in rocks that experienced high-temperature amphibolite-facies overprinting and those that underwent low-temperature blueschist-facies overprinting from a variety of subduction-related terranes worldwide. We find that trace element distributions approach equilibrium partition coefficients in rocks from amphibolite-facies overprinted terranes, whereas trace element distributions do not approach equilibrium in rocks that experienced blueschist-facies overprinting. We caution that single element thermometers that rely upon slow-diffusing high field strength elements should not be applied to rocks equilibrated at <600 °C unless attainment of trace element equilibrium can be demonstrated.
Earth and Planetary Science Letters, Vol. 500, pp. 86-96.
Mantle
geothermometry
Abstract: The seismically-observed large low shear velocity provinces in the Earth's lowermost mantle have been hypothesized to be caused by thermochemical piles of compositionally distinct, more-primitive material which may be remnants of Earth's early differentiation. However, one critical question is how the Earth's thermal evolution is affected by the long-term presence of the large-scale compositional heterogeneity in the lowermost mantle. Here, we perform geodynamical calculations to investigate the time evolution of the morphology of large-scale compositional heterogeneity and its influence on the Earth's long-term thermal evolution. Our results show that a global layer of intrinsically dense material in the lowermost mantle significantly suppresses the CMB heat flux, which leads to faster cooling of the background mantle relative to an isochemical mantle. As the background mantle cools, the intrinsically dense material is gradually pushed into isolated thermochemical piles by cold downwellings. The size of the piles also decreases with time due to entraining of pile material into the background mantle. The morphologic change of the accumulations of intrinsic dense material eventually causes a gradual increase of CMB heat flux, which significantly reduces the cooling rate of Earth's mantle.
The utility of clinopyroxene in diamond exploration.
2018 Yellowknife Geoscience Forum , p. 13. abstract
Global
thermobarometry
Abstract: Clinopyroxene single-crystal thermo-barometry is an essential tool in the identification and evaluation of prospective kimberlites. The paleogeothermal gradient preserved by clinopyroxene xenocrysts elucidates the thermal structure of the underlying lithospheric mantle; indicates the depth to and thickness of the “diamond window”. The widely used clinopyroxene thermometer-barometer of Nimis and Taylor (2000) requires that clinopyroxene equilibrated with both garnet and orthopyroxene. With the rare exception of wehrlites, equilibration with orthopyroxene is nearly a given for the majority of chrome-diopside clinopyroxene xenocrysts. Demonstrating equilibration with garnet, however, is a major obstacle for clinopyroxene-based thermobarometry. The most commonly used method for clinopyroxene discrimination is an Al2O3-Cr2O3 diagram proposed by Ramsay and Thompkins in 1994 supplemented with an additional MgO-Al2O3 from Nimis (1998) and an additional 1-dimensional filter based on chemical composition. Despite the aggressiveness of the filtering method, single-clinopyroxene pressure-temperature results have large scatter that can obscure the true paleogeothermal gradient. This is especially true of areas where the lithospheric mantle has undergone chemical modification by melt/fluid influx. Using a database of clinopyroxenes derived from kimberlite-borne mantle-derived lherzolites, we have developed a simple and effective discrimination plot that identifies clinopyroxene from garnet lherzolites and simultaneously removes clinopyroxene from metaosomatised peridotites. Calculated paleogeothermal gradients from clinopyroxene xenocrysts cut across model conductive geotherms which can complicate the interpretation of thermobarometry data. Grütter (2009) presented a solution to the problem by way of relative reference geotherms. He used xenocryst data from three Canadian locations with different thermal structures as references in comparison to the dataset under investigation. Taking a cue from this earlier work, we have developed a new set of relative reference geotherms that are based on single-clinopyroxene thermobarometry data for xenoliths from well-characterized regions - Somerset Island, Kaapvaal on-craton, and the Central Slave. A simple linear fit through the data produces sub-parallel clinopyroxene reference geotherms that are simpler to use and easier to visualize compared to the xenocrysts reference geotherms. Using these two new and simple tools will greatly help maximize the utility of clinopyroxene data in large exploration databases.
Doklady Earth Sciences, Vol. 482, 2, pp. 1317-1319.
Mantle
geothermometry
Abstract: It is known that the ?-? parameters of diamond-bearing kimberlite xenoliths correspond to subductive paleogeotherms lying between the 36 and 41 mW/m2 conductive models. There are some studies showing the correlation of diamond ability with oxygen fugacity and the fluid composition of mantle xenoliths.The most diamondiferous samples correspond to the water compositions of the calculated O-H-C fluid with a minimum atomic carbon content in it. From the calculations it follows that the fluid carbon atomic content increases with a temperature increase and with the pressure decreasing. The most minor C contents have the 35 mW/m2 conductive model in comparison with the 40 and 45 mW/m2 models. As a result, it is possible to conclude that the low temperature fields (less than 1100°C) of the “cold” geotherms have the highest diamondiferous ability.
Earth and Planetary Science Letters, Vol. 490, 1, pp. 161-169.
Mantle
geothermometry
Abstract: In this study, we have determined the phase boundary between Mg0.735Fe0.21Al0.07Si0.965O3-Bm and PPv and the thermal equations of state of both phases up to 202 GPa and 2600 K using synchrotron X-ray diffraction in laser heated diamond anvil cells. Our experimental results have shown that the combined effect of Fe and Al produces a wide two-phase coexistence region with a thickness of 26 GPa (410 km) at 2200 K, and addition of Fe lowers the onset transition pressure to 98 GPa at 2000 K, consistent with previous experimental results. Furthermore, addition of Fe was noted to reduce the density (?) and bulk sound velocity () contrasts across the Bm-PPv phase transition, which is in contrast to the effect of Al. Using the obtained phase diagram and thermal equations of state of Bm and PPv, we have also examined the effect of composition variations on the ? and
profiles of the lowermost mantle. Our modeling results have shown that the pyrolitic lowermost mantle should be highly heterogeneous in composition and temperature laterally to match the observed variations in the depth and seismic signatures of the D? discontinuity. Normal mantle in a pyrolitic composition with ?10% Fe and Al in Bm and PPv will lack clear seismic signature of the D? discontinuity because the broad phase boundary could smooth the velocity contrast between Bm and PPv. On the other hand, Fe-enriched regions close to the cold slabs may show a seismic signature with a change in the velocity slope of the D? discontinuity, consistent with recent seismic observations beneath the eastern Alaska. Only regions depleted in Fe and Al near the cold slabs would show a sharp change in velocity. Fe in such regions could be removed to the outer core by strong core-mantle interactions or partitions together with Al to the high-pressure phases in the subduction mid ocean ridge basalts. Our results thus have profound implication for the composition of the lowermost mantle.
Abstract: Thermal history models, that have been used to understand the geological history of Earth, are now being coupled to climate models to map conditions that allow planets to maintain surface water over geologic time - a criteria considered crucial for life. However, the lack of intrinsic uncertainty assessment has blurred guidelines for how thermal history models can be used toward this end. A model, as a representation of something real, is not expected to be complete. Unmodeled effects are assumed to be small enough that the model maintains utility for the issue(s) it was designed to address. The degree to which this holds depends on how unmodeled factors affect the certainty of model predictions. We quantify this intrinsic uncertainty for several parameterized thermal history models (a widely used subclass of planetary models). Single perturbation analysis is used to determine the reactance time of different models. This provides a metric for how long it takes low amplitude, unmodeled effects to decay or grow. Reactance time is shown to scale inversely with the strength of the dominant feedback (negative or positive) within a model. A perturbed physics analysis is then used to determine uncertainty shadows for model outputs. This provides probability distributions for model predictions and tests the structural stability of a model. That is, do model predictions remain qualitatively similar, and within assumed model limits, in the face of intrinsic uncertainty. Once intrinsic uncertainty is accounted for, model outputs/predictions and comparisons to observational data should be treated in a probabilistic way.
Earth and Planetary Science Letters, Vol. 505, pp. 65-75.
Mantle
geothermometry
Abstract: It is estimated that around three quarters of Earth's first generation continental crust had been produced by the end of the Archaean Eon, 2.5 billion years ago. This ancient continental crust is mostly composed of variably deformed and metamorphosed magmatic rocks of the tonalite-trondhjemite-granodiorite (TTG) suite that formed by partial melting of hydrated mafic rocks. However, the geodynamic regime under which TTG magmas formed is a matter of ongoing debate. Using a filtered global geochemical dataset of 563 samples with ages ranging from the Eoarchaean to Neoarchaean (4.0-2.5 Ga), we interrogate the bulk rock major oxide and trace element composition of TTGs to assess evidence for secular change. Despite a high degree of scatter in the data, the concentrations or ratios of several key major oxides and trace elements show statistically significant trends that indicate maxima, minima and/or transitions in the interval 3.3-3.0 Ga. Importantly, a change point analysis of K2O/Na2O, Sr/Y and LaN/YbN demonstrates a statistically significant (>99% confidence) change during this 300 Ma period. These shifts may be linked to a fundamental change in geodynamic regime around the peak in upper mantle temperatures from one dominated by non-uniformitarian, deformable stagnant lid processes to another dominated by the emergence of global mobile lid or plate tectonic processes by the end of the Archaean. A notable change is also evident at 2.8-2.7 Ga that coincides with a major jump in the rate of survival of metamorphic rocks with contrasting thermal gradients, which may relate to the emergence of more potassic continental arc magmas and an increased preservation potential during collisional orogenesis. In many cases, the chemical composition of TTGs shows an increasing spread through the Archaean, reflecting the irreversible differentiation of the lithosphere.
Geochemistry, Geophysics, Geosystems, Vol. 20, 2, pp. 673-687.
Mantle
Thermometry
Abstract: We reevaluate the effects of slab age, speed, and dip on slab temperature with numerical models. The thermal parameter ? = t v sin ?, where t is age, v is speed, and ? is angle, is traditionally used as an indicator of slab temperature. However, we find that an empirically derived quantity, in which slab temperature T ? log (t?av?b) , is more accurate at depths <120 km, with the constants a and b depending on position within the slab. Shallower than the decoupling depth (~70-80 km), a~1 and b~0, that is, temperature is dependent on slab age alone. This has important implications for the early devolatilization of slabs in the hottest (youngest) cases and for shallow slab seismicity. At subarc depths (~100 km), within the slab mantle, a~1 and b~0 again. However, for the slab crust, now a~0.5 and b~1, that is, speed has the dominant effect. This is important when considering the generation of arc magmatism, and in particular, slab melting and the generation of slab?derived melange diapirs. Moving deeper into the Earth, the original thermal parameter performs well as a temperature indicator, initially in the core of the slab (the region of interest for deep water cycling). Finally, varying the decoupling depth between 40 and 100 km has a dominant effect on slab temperatures down to 140?km depth, but only within the slab crust. Slab mantle temperature remains primarily dependent on age.
Contributions to Mineralogy and Petrology, Vol. 173, pp. 39-
Global
FTIR
Abstract: Platelets are one of the most common defects occurring in natural diamonds but their behaviour has not previously been well understood. Recent technical advances, and a much improved understanding of the correct interpretation of the main infrared (IR) feature associated with platelets (Speich et al. 2017), facilitated a systematic study of platelets in 40 natural diamonds. Three different types of platelet behaviour were identified here. Regular diamonds show linear correlations between both B-centre concentrations and platelet density and also between platelet size and platelet density. Irregular diamonds display reduced platelet density due to platelet breakdown, anomalously large or small platelets and a larger platelet size distribution. These features are indicative of high mantle storage temperatures. Finally, a previously unreported category of subregular diamonds is defined. These diamonds experienced low mantle residence temperatures and show smaller than expected platelets. Combining the systematic variation in platelet density with temperatures of mantle storage, determined by nitrogen aggregation, we can demonstrate that platelet degradation proceeds at a predictable rate. Thus, in platelet-bearing diamonds where N aggregation is complete, an estimate of annealing temperature can now be made for the first time.
Abstract: We examine the partial melting and the cooling history of a ~5?km section of mantle lithosphere preserved in the Donjek massif, part of a Permian ophiolite in the northern Cordillera of Yukon, Canada. The mantle rocks are depleted spinel harzburgite containing <3% clinopyroxene displaying steep rare-earth element (REE) chondrite-normalized profiles and low (Gd/Yb)n (0.02 to 0.07) compared to most other ophiolites. The REE patterns of clinopyroxene can be modeled as 16-20% partial melts of typical depleted mid-ocean ridge (MOR) mantle. The REE exchange between coexisting ortho- and clinopyroxene preserves temperatures (TREE) of 1150-1360?°C, some of the highest values recorded in ophiolites and abyssal peridotites, and show a positive correlation with CaMg exchange (solvus) temperatures (TBKN) of 900-970?°C. The harzburgite represents lithosphere formed at an initial melting temperature of ~ 1350?°C that cooled at rate of 10?1 to 10?4?°C/year as deduced by TREE values with cation diffusion and grain size data. The TREE temperatures and cooling rates for the Donjek massif show a regular systematic variation with depth from the crust-mantle transition along a trend similar to the Samail ophiolite of Oman, consistent with conductive heat transfer beneath a cool lower crust. High near-solidus temperatures and the cooling rates in the massif were a consequence of rapid obduction against oceanic crust along either a transform or low angle detachment soon after melt extraction. Final emplacement of the ophiolite as klippen on underlying continental crust occurred ~ 40?m.y. later.
Geochemistry, Geophysics, Geosystems, Vol. 20, 3, pp. 1358-1386.
Mantle
geothermometry
Abstract: Continents influence the mantle's convective wavelength and the heat flow escaping from the planet's surface. Over the last few decades, many numerical and analytical studies have contributed to the debate about whether the continents can warm up the subcontinental mantle or not and if they do, then to what extent? However, a consensus regarding the exact nature and magnitude of this correlation between continents and elevated temperatures in the subcontinental mantle remains to be achieved. By conducting a systematic parameter study using 2?D global mantle convection simulations with mobile continents, we provide qualitative and quantitative observations on the nature of this correlation. In our incompressible and compressible convection models, we observe the general processes of downwellings bringing cold material into the mantle along continental margins and a subsequent buildup of warm thermal anomalies underneath the continents. We compute the amplitude and degree of this correlation using spectral decomposition of the temperature and composition fields. The dominant degree of correlation evolves with time and changes with continental configuration. Using simple empirical fits, we observe that this correlation decreases with increasing core temperature, number of continents, internal heating, or decreasing reference viscosity. We also report simple regressions of the time dependence of this correlation. Additionally, we show that decompression melting as a result of a mantle upwelling or small?scale sublithospheric convection leads to voluminous volcanism. The emplacement of this dense basalt?eclogite material breaks the continents apart and destroys the correlation.
Geophysical Research Letters, Vol. 46, 7, pp. 3652-3662.
Mantle
geothermometry
Abstract: High mountains in Norway have long puzzled scientists because it is challenging to explain their existence. Numerous explanations have been proposed including processes deep inside the Earth. Our results show that these processes must be located above 410?km depth. This observation is critical for the ongoing debate on the cause of the enigmatic mountains in Scandinavia. New data acquired between 2012 and 2017 by the collaborative ScanArray project between European institutions allow mapping of the mantle transition zone—the deepest layer possibly involved in the mountain support. We show that the mantle transition zone boundaries beneath Fennoscandia are close to reference depths and the zone has a standard thickness. As the depths to these boundaries are sensitive to temperatures, this indicates that the mantle transition zone in this area is unaffected by any ongoing deep process. Therefore, the explanation for the high topography in Norway must be found above the mantle transition zone. This study provides the first map of the mantle transition zone below Fennoscandia, which will be valuable for any further global studies of the mantle transition zone.
Geochimica et Cosmochimica Acta, Vol. 253, pp. 290-306.
Mantle
geothermometry
Abstract: Carbon isotope exchange between carbon-bearing high temperature phases records the carbon (re-) processing in the Earth's interior, where the vast majority of global carbon is stored. Redox reactions between carbonate phases and elemental carbon govern the mobility of carbon, which then can be traced by its isotopes. We determined the carbon isotope fractionation factor between graphite and a Na2CO3-CaCO3 melt at 900-1500?°C and 1?GPa; The failure to isotopically equilibrate preexisting graphite led us to synthesize graphite anew from organic material during the melting of the carbonate mixture. Graphite growth proceeds by (1) decomposition of organic material into globular amorphous carbon, (2) restructuring into nano-crystalline graphite, and (3) recrystallization into hexagonal graphite flakes. Each transition is accompanied by carbon isotope exchange with the carbonate melt. High-temperature (1200-1500?°C) equilibrium isotope fractionation with type (3) graphite can be described by (temperature T in K). As the experiments do not yield equilibrated bulk graphite at lower temperatures, we combined the ?1200?°C experimental data with those derived from upper amphibolite and lower granulite facies carbonate-graphite pairs (Kitchen and Valley, 1995; Valley and O'Neil, 1981). This yields the general fractionation function usable as a geothermometer for solid or liquid carbonate at ?600?°C. Similar to previous observations, lower-temperature experiments (?1100?°C) deviate from equilibrium. By comparing our results to diffusion and growth rates in graphite, we show that at ?1100?°C carbon diffusion is slower than graphite growth, hence equilibrium surface isotope effects govern isotope fractionation between graphite and carbonate melt and determine the isotopic composition of newly formed graphite. The competition between diffusive isotope exchange and growth rates requires a more careful interpretation of isotope zoning in graphite and diamond. Based on graphite crystallization rates and bulk isotope equilibration, a minimum diffusivity of Dgraphite?=?2?×?10?17 m2s?1 for T?>?1150?°C is required. This value is significantly higher than calculated from experimental carbon self-diffusion constants (?1.6?×?10?29?m2?s?1) but in good agreement with the value calculated for mono-vacancy migration (?2.8?×?10?16?m2?s?1).
Abstract: The majority of continental crust formed during the hotter Archean was composed of Tonalite-Trondhjemite-Granodiorite (TTG) rocks. In contrast to the present-day loci of crust formation around subduction zones and intra-plate tectonic settings, TTGs are formed when hydrated basalt melts at garnet-amphibolite, granulite or eclogite facies conditions. Generating continental crust requires a two step differentiation process. Basaltic magma is extracted from the pyrolytic mantle, is hydrated, and then partially melts to form continental crust. Here, we parameterise the melt production and melt extraction processes and show self-consistent generation of primordial continental crust using evolutionary thermochemical mantle convection models. To study the growth of TTG and the geodynamic regime of early Earth, we systematically vary the ratio of intrusive (plutonic) and eruptive (volcanic) magmatism, initial core temperature, and internal friction coefficient. As the amount of TTG that can be extracted from the basalt (or basalt-to-TTG production efficiency) is not known, we also test two different values in our simulations, thereby limiting TTG mass to 10% or 50% of basalt mass. For simulations with lower basalt-to-TTG production efficiency, the volume of TTG crust produced is in agreement with net crustal growth models but overall crustal (basaltic and TTG) composition stays more mafic than expected from geochemical data. With higher production efficiency, abundant TTG crust is produced, with a production rate far exceeding typical net crustal growth models but the felsic to mafic crustal ratio follows the expected trend. These modelling results indicate that (i) early Earth exhibited a "plutonic squishy lid" or vertical-tectonics geodynamic regime, (ii) present-day slab-driven subduction was not necessary for the production of early continental crust, and (iii) the Archean Earth was dominated by intrusive magmatism as opposed to "heat-pipe" eruptive magmatism.
Abstract: arth’s magnetic field is generated by turbulent motion in its fluid outer core. Although the bulk of the outer core is vigorously convecting and well mixed, some seismic, geomagnetic and geodynamic evidence suggests that a global stably stratified layer exists at the top of Earth’s core. Such a layer would strongly influence thermal, chemical and momentum exchange across the core-mantle boundary and thus have important implications for the dynamics and evolution of the core. Here we argue that the relevant scenario is not global stratification, but rather regional stratification arising solely from the lateral variations in heat flux at the core-mantle boundary. Using our extensive suite of numerical simulations of the dynamics of the fluid core with heterogeneous core-mantle boundary heat flux, we predict that thermal regional inversion layers extend hundreds of kilometres into the core under anomalously hot regions of the lowermost mantle. Although the majority of the outermost core remains actively convecting, sufficiently large and strong regional inversion layers produce a one-dimensional temperature profile that mimics a globally stratified layer below the core-mantle boundary—an apparent thermal stratification despite the average heat flux across the core-mantle boundary being strongly superadiabatic.
Abstract: Sulphides are the most common type of inclusions found in diamonds and are widely used to determine the timing and lithology of diamond formation. Typical inclusions are monosulfide solid solutions (MSS) in the Fe-Ni-S system with minor amounts of Cu, Co and Mo. Previous experimental studies show that oxygen partitions into sulphide melts but most importantly measurements of natural sulphide inclusions indeed show measureable oxygen concentrations. If the parameters that control sulphide oxygen concentration can be determined then they could be potentially used to understand formation conditions of diamonds. We performed a series of high pressure (3-11 GPa) and high temperature (1573-1973 K) experiments in order to parameterize the oxygen content in sulphides in equilibrium with a mantle peridotite assemblage relevant to diamond formation. Multi-anvil experiments were carried out in graphite capsules and a peridotite silicate composition was equilibrated with molten FeS for at least 5 hrs. Run products that contained mantle silicate minerals and quenched sulphide melts were analysed using the electron microprobe. In some cases Ir was added in sufficient quantities to saturate the sulphides and form an Fe-Ir alloy from which the oxygen fugacity could be accurately determined. We measured up to 16 weight % of FeO in our experimental sulphide melts at mantle conditions. Moreover, the content of oxygen in the sulphide is found to be not controlled by fO2 or fS2, which is in disagreement with previous experimental studies conducted at ambient pressure conditions. The experiments indicate that the oxygen concentration is mainly controlled by the FeO activity in coexisting silicate phases and the temperature. In order to fit the data and to account for the observed FeO dependence, we developed a thermodynamic model using an end-member equilibrium between olivine, pyroxene and FeO in the sulphide melt. Using this relationship with measurements of oxygen in natural sulphide inclusions in diamonds reveals temperatures for lithospheric diamond formation in the range of 1140 – 1410 ºC.
Abstract: In the Archaean, global surface heat flow was substantially higher than today because of greater internal radiogenic heat production and primordial heat content within the Earth. Nonetheless, the lithospheric roots of Archaean cratons were apparently surprisingly cool, recording similarly low ambient temperatures to those inferred today, allowing e.g. for the stabilisation of diamond. This finding is seemingly in conflict with a generally ‘hotter’ Archaean mantle, as is widely postulated, but the paradox could be explained if the sub-cratonic lithospheric mantle was substantially thicker in the Archaean than today. Here, we report a re-investigation of the thermal structure of the Archaean Kalahari lithosphere using published and unpublished petrological data of diamond inclusions indicated to be of Archaean age. Our thermobarometric calculations agree with earlier findings that the Archaean cratonic mantle root was surprisingly cool. Importantly, the shape of the inclusion-derived P-T array deviates from the modern geotherm recorded by peridotite xenoliths. Specifically, diamond inclusions define a systematically steeper geothermal gradient than is observed in cratonic xenoliths. We find that Archaean diamond inclusion and modern xenolith P-T data cannot be reconciled by a single steady-state geotherm. The P-T conditions recorded in diamond inclusions are incompatible with the current characteristically low present-day heat-production of the overlying crust. Instead, the steeper geotherm implies high heat production in the crust during diamond formation and the distinctive geothermal gradient recorded in the studied diamond inclusions could reflect ancient mantle conditions. We modelled a suite of ‘fossil’ geotherms, with increased radiogenic heat production within the crust during the Archaean. Solutions providing very good fits with the diamond inclusion geotherm all require that the Archaean lithosphere must have extended to far greater depths than is preserved today. The required depth ranges from ~ 300 km to ~ 450 km depth, for a modern (~ 1350°C) and a significantly hotter (~ 1600°C) mantle potential temperature, respectively. In either case, it is clear that the Kalahari lithosphere must have experienced significant (at least 100 km) basal erosion subsequent to its formation.
Abstract: Thermomechanical models of mantle convection and melting in an inferred hotter Archean Earth show the emergence of pressure-temperature (P-T) regimes that resemble present-day plate tectonic environments yet developed within a non-plate tectonics regime. The models’ P-T gradients are compatible with those inferred from evolving tonalite-trondhjemite-granodiorite series rocks and the paired metamorphic belt record, supporting the feasibility of divergent and convergent tectonics within a mobilized, yet laterally continuous, lithospheric lid. “Hot” P-T gradients of 10-20 °C km-1 form along asymmetric lithospheric drips, then migrate to areas of deep lithospheric downwelling within ?300-500 m.y., where they are overprinted by high-pressure warm and, later, cold geothermal signatures, up to ?8 °C km-1. Comparisons with the crustal production and reworking record suggest that this regime emerged in the Hadean.
Abstract: Mineral inclusions entrapped in other minerals may record the local stresses at the moment of their entrapment in the deep Earth. When rocks are exhumed to the surface of the Earth, residual stresses and strains may still be preserved in the inclusion. If measured and interpreted correctly through elastic geobarometry, they give us invaluable information on the pressures (P) and temperatures (T) of metamorphism. Current estimates of P and T of entrapment rely on simplified models that assumes that the inclusion is spherical and embedded in an infinite host, and that their elastic properties are isotropic. We report a new method for elastic geobarometry for anisotropic inclusions in quasi-isotropic hosts. The change of strain in the inclusion is modelled with the axial equations of state of the host and the inclusion. Their elastic interaction is accounted for by introducing a 4th rank tensor, the relaxation tensor, that can be evaluated numerically for any symmetry of the host and the inclusion and for any geometry of the system. This approach can be used to predict the residual strain/stress state developed in an inclusion after exhumation from known entrapment conditions, or to estimate the entrapment conditions from the residual strain measured in real inclusions. In general, anisotropic strain and stress states are developed in non-cubic mineral inclusions such as quartz and zircon, with deviatoric stresses typically limited to few kbars. For garnet hosts, the effect of the mutual crystallographic orientation between the host and the inclusion on the residual strain and stress is negligible when the inclusion is spherical and isolated. Assuming external hydrostatic conditions, our results suggest that the isotropic and the new anisotropic models give estimations of entrapment conditions within 2%.
Geophysical Research Letters, Vol. 46, 20, pp. 11065-110670.
Mantle
geothermometry
Abstract: Earth's magnetic field is produced by a dynamo in the core that requires motion of the fluid Fe alloy. Both thermal convection, arising from the transport of heat in excess of conducted heat, and compositional convection, arising from light element exsolution at the freezing inner core boundary, are suggested as energy sources. The contribution of thermal convection (possibly ranging from nothing to significant) depends on thermal conductivity of the outer core. Our experimental measurements of electrical resistivity of solid and liquid Fe at high pressures show that resistivity is constant along the pressure?dependent melting boundary of Fe. Using our derived thermal conductivity value at the inner core (freezing) boundary, we calculate the heat conducted in the liquid outer core and find that thermal convection is needed to carry additional heat through the outer core to match the heat extracted through the core?mantle boundary.
Abstract: The heat flux across the core-mantle boundary (QCMB) is the key parameter to understand the Earth/s thermal history and evolution. Mineralogical constraints of the QCMB require deciphering contributions of the lattice and radiative components to the thermal conductivity at high pressure and temperature in lower mantle phases with depth-dependent composition. Here we determine the radiative conductivity (krad) of a realistic lower mantle (pyrolite) in situ using an ultra-bright light probe and fast time-resolved spectroscopic techniques in laser-heated diamond anvil cells. We find that the mantle opacity increases critically upon heating to ~3000 K at 40-135 GPa, resulting in an unexpectedly low radiative conductivity decreasing with depth from ~0.8 W/m/K at 1000 km to ~0.35 W/m/K at the CMB, the latter being ~30 times smaller than the estimated lattice thermal conductivity at such conditions. Thus, radiative heat transport is blocked due to an increased optical absorption in the hot lower mantle resulting in a moderate CMB heat flow of ~8.5 TW, at odds with present estimates based on the mantle and core dynamics. This moderate rate of core cooling implies an inner core age of about 1 Gy and is compatible with both thermally- and compositionally-driven ancient geodynamo.
Earth and Planetary Science Letters, Vol. 536, 116161 7p. Pdf
Mantle
geothermometry
Abstract: Knowledge of thermal conductivity of mantle minerals is crucial for understanding heat transport from the Earth's core to mantle. At the pressure-temperature conditions of the Earth's core-mantle boundary, calculations of lattice thermal conductivity based on atomistic models have determined values ranging from 1 to 14 W/m/K for bridgmanite and bridgmanite-rich mineral assemblages. Previous studies have been performed at room temperature up to the pressures of the core-mantle boundary, but correcting these to geotherm temperatures may introduce large errors. Here we present the first measurements of lattice thermal conductivity of mantle minerals up to pressures and temperatures near the base of the mantle, 120 GPa and 2500 K. We use a combination of continuous and pulsed laser heating in a diamond anvil cell to measure the lattice thermal conductivity of pyrolite, the assemblage of minerals expected to make up the lower mantle. We find a value of W/m/K at 80 GPa and 2000 to 2500 K and 5.9 W/m/K at 124 GPa and 2000 to 3000 K. These values rule out the highest calculations of thermal conductivity of the Earth's mid-lower mantle (i.e. W/m/K at 80 GPa), and are consistent with both the high and low calculations of thermal conductivity near the base of the lower mantle.
Precambrian Research, doi.org/10.1016/j.precamres.2020.105703 in press available 80p. Pdf
Canada
geothermometry
Abstract: The northern Canadian Shield is comprised of multiple Archean cratons that were sutured by the late Paleoproterozoic to form the Canadian component of supercontinent Nuna. More than 2000 combined K-Ar and 40Ar/39Ar cooling ages from across the region reveal a stark contrast in upper and lower plate thermal responses to Nuna-forming events, with the Churchill Province in particular revealing near complete thermal reworking during the late Paleoproterozoic. We review the detailed cooling history for five regions that span the Churchill Province and Trans-Hudson orogen (THO): Thelon Tectonic Zone, South Rae, Reindeer Zone, South Hall Peninsula, and the Cape Smith Belt. The cooling patterns across Churchill Province are revealed in two >1500 km transects. At the plate scale, Churchill’s cooling history is dominated by THO accretionary and collisional events, during which it formed the upper plate. Cooling ages generally young from west to east across both southern and central Churchill, and latest cooling in the THO is 50 myr older in southernmost Churchill (Reindeer Zone) compared to eastern Churchill (Hall Peninsula), indicating diachronous thermal equilibration across 2000 km strike length of the THO. Churchill exhibits relatively high post-terminal THO cooling rates of ~4 °C/myr, which support other geological evidence for widespread rapid exhumation of the THO upper plate following terminal collision, potentially in response to lithospheric delamination.
Earth and Planetary Science Letters, Vol. 537, 116171 12p. Pdf
Mantle
geothermometry
Abstract: A better understanding of the Earth's compositional structure is needed to place the geochemical record of surface rocks into the context of Earth accretion and evolution. Cosmochemical constraints imply that lower-mantle rocks may be enriched in silica relative to upper-mantle pyrolite, whereas geophysical observations support whole-mantle convection and mixing. To resolve this discrepancy, it has been suggested that subducted mid-ocean ridge basalt (MORB) segregates from subducted harzburgite to accumulate in the mantle transition zone (MTZ) and/or the lower mantle. However, the key parameters that control basalt segregation and accumulation remain poorly constrained. Here, we use global-scale 2D thermochemical convection models to investigate the influence of mantle-viscosity profile, planetary-tectonic style and bulk composition on the evolution and distribution of mantle heterogeneity. Our models robustly predict that, for all cases with Earth-like tectonics, a basalt-enriched reservoir is formed in the MTZ, and a harzburgite-enriched reservoir is sustained at 660?800 km depth, despite ongoing whole-mantle circulation. The enhancement of basalt and harzburgite in and beneath the MTZ, respectively, are laterally variable, ranging from ?30% to 50% basalt fraction, and from ?40% to 80% harzburgite enrichment relative to pyrolite. Models also predict an accumulation of basalt near the core mantle boundary (CMB) as thermochemical piles, as well as moderate enhancement of most of the lower mantle by basalt. While the accumulation of basalt in the MTZ does not strongly depend on the mantle-viscosity profile (explained by a balance between basalt delivery by plumes and removal by slabs at the given MTZ capacity), that of the lowermost mantle does: lower-mantle viscosity directly controls the efficiency of basalt segregation (and entrainment) near the CMB; upper-mantle viscosity has an indirect effect through controlling slab thickness. Finally, the composition of the bulk-silicate Earth may be shifted relative to that of upper-mantle pyrolite, if indeed significant reservoirs of basalt exist in the MTZ and lower mantle.
Earth and Planetary Science Letters, Vol. 536, 116161, 11p. Pdf
Mantle
geothermometry
Abstract: Knowledge of thermal conductivity of mantle minerals is crucial for understanding heat transport from the Earth's core to mantle. At the pressure-temperature conditions of the Earth's core-mantle boundary, calculations of lattice thermal conductivity based on atomistic models have determined values ranging from 1 to 14 W/m/K for bridgmanite and bridgmanite-rich mineral assemblages. Previous studies have been performed at room temperature up to the pressures of the core-mantle boundary, but correcting these to geotherm temperatures may introduce large errors. Here we present the first measurements of lattice thermal conductivity of mantle minerals up to pressures and temperatures near the base of the mantle, 120 GPa and 2500 K. We use a combination of continuous and pulsed laser heating in a diamond anvil cell to measure the lattice thermal conductivity of pyrolite, the assemblage of minerals expected to make up the lower mantle. We find a value of W/m/K at 80 GPa and 2000 to 2500 K and 5.9 W/m/K at 124 GPa and 2000 to 3000 K. These values rule out the highest calculations of thermal conductivity of the Earth's mid-lower mantle (i.e. W/m/K at 80 GPa), and are consistent with both the high and low calculations of thermal conductivity near the base of the lower mantle.
Rendiconti Lincei. Scienze Fisiche e Naturali *** in Eng., doi.org/10.1007/ s12210-020-00897-8 9p. Pdf
Mantle
geothermometry
Abstract: Natural diamonds and their mineral inclusions represent the deepest regions of our planet. Diamonds form between about 120/130 km in the upper mantle and possibly up to 1000 km depth in the lower mantle, and they can transport mineral inclusions, which are fragments directly from regions that are inaccessible to geologists. Diamond-inclusion system is a very precious geological object not only due to the depth of provenance in the mantle but also because of the diamond age, which ranges between 3.6 and 0.09 Ga providing information over a very wide evolution time of the Earth. It is evident that the determination of the depth of formation of the diamond-inclusion system is one of the crucial aspects to retrieve fundamental geological information about the deep mantle. However, the determination of such depth is not trivial at all and different approaches could be adopted; one of the most promising is represented by the so-called "elastic geobarometry". In this review, I will focus on elastic geobarometry and on the role that anisotropy has on the determination of the depth of diamond formation. The work will also provide a short overview of the most common approaches used to retrieve the depth of diamond formation.
Abstract: Paleo-temperature data indicates that the Earth's mantle did not cool at a constant rate over geologic time. Post magma ocean cooling was slow with an onset of more rapid mantle cooling between 2.5 and 3.0 Gyr. We explore the hypothesis that this multi-stage cooling is a result of deep water cycling coupled to thermal mantle convection. As warm mantle ascends, producing melt, the mantle is dehydrated. This tends to stiffens the mantle, which slows convective vigor causing mantle heating. At the same time, an increase in temperature tends to lower mantle viscosity which acts to increase convective vigor. If these two tendencies are in balance, then mantle cooling can be weak. If the balance is broken, by a switch to a net rehydration of the mantle, then the mantle can cool more rapidly. We use coupled water cycling and mantle convection models to test the viability of this hypothesis. We test models with different parameterizations to allow for variable degrees of plate margin strength. We also perform a layered uncertainty analysis on all the models to account for input, parameter, and structural model uncertainties. Within model and data uncertainty, the hypothesis that deep water cycling, together with a combination of plate strength and mantle viscosity resisting mantle overturn, can account for paleo data constraints on mantle cooling.
Proceedinds of the National Academy of Sciences, Vol. 117, 19, 9p. Pdf
United States, Colorado
geothermometry
Abstract: The Great Unconformity marks a major gap in the continental geological record, separating Precambrian basement from Phanerozoic sedimentary rocks. However, the timing, magnitude, spatial heterogeneity, and causes of the erosional event(s) and/or depositional hiatus that lead to its development are unknown. We present field relationships from the 1.07-Ga Pikes Peak batholith in Colorado that constrain the position of Cryogenian and Cambrian paleosurfaces below the Great Unconformity. Tavakaiv sandstone injectites with an age of ?676 ± 26 Ma cut Pikes Peak granite. Injection of quartzose sediment in bulbous bodies indicates near-surface conditions during emplacement. Fractured, weathered wall rock around Tavakaiv bodies and intensely altered basement fragments within unweathered injectites imply still earlier regolith development. These observations provide evidence that the granite was exhumed and resided at the surface prior to sand injection, likely before the 717-Ma Sturtian glaciation for the climate appropriate for regolith formation over an extensive region of the paleolandscape. The 510-Ma Sawatch sandstone directly overlies Tavakaiv-injected Pikes granite and drapes over core stones in Pikes regolith, consistent with limited erosion between 717 and 510 Ma. Zircon (U-Th)/He dates for basement below the Great Unconformity are 975 to 46 Ma and are consistent with exhumation by 717 Ma. Our results provide evidence that most erosion below the Great Unconformity in Colorado occurred before the first Neoproterozoic Snowball Earth and therefore cannot be a product of glacial erosion. We propose that multiple Great Unconformities developed diachronously and represent regional tectonic features rather than a synchronous global phenomenon.
Abstract: Raman spectroscopy is widely applied in metamorphic petrology and offers many opportunities for geological and tectonic research. Minimal sample preparation preserves sample integrity and microtextural information, while use with confocal microscopes allows spatial resolution down to the micrometer level. Raman spectroscopy clearly distinguishes mineral polymorphs, providing crucial constraints on metamorphic conditions, particularly ultrahigh-pressure conditions. Raman spectroscopy can also be used to monitor the structure of carbonaceous material in metamorphic rocks. Changes in structure are temperature-sensitive, so Raman spectroscopy of carbonaceous material is widely used for thermometry. Raman spectroscopy can also detect and quantify strain in micro-inclusions, offering new barometers that can be applied to understand metamorphic and tectonic processes without any assumptions about chemical equilibrium.
Earth and Planetary Science Letters, Vol. 542, 116317 16p. Pdf
Mantle
plumes, geothermometry
Abstract: Mantle convection is the principal mechanism by which heat is transferred from the deep Earth to the surface. Cold subducting slabs sink into the mantle and steadily warm, whilst upwelling plumes carry heat to the base of lithospheric plates where it can subsequently escape by conduction. Accurate estimation of the total heat carried by these plumes is important for understanding geodynamic processes and Earth's thermal budget. Existing estimates, based upon swell geometries and velocities of overriding plates, yield a global heat flux of ?2 TW and indicate that plumes play only a minor role in heat transfer. Here, we revisit the Icelandic and Hawaiian plumes to show that their individual flux estimates are likely to be incorrect due to the assumption that buoyancy is mainly produced within the lithosphere and therefore translates at plate velocities. We develop an alternative methodology that depends upon swell volume, is independent of plate velocities, and allows both for decay of buoyancy through time and for differential motion between asthenospheric buoyancy and the overlying plate. Reanalysis of the Icelandic and Hawaiian swells yields buoyancy fluxes of Mg s?1 and Mg s?1, respectively. Both swells are used to calibrate a buoyancy decay timescale of ?45 Myr for the new volumetric approach, which enables buoyancy fluxes to be estimated for a global inventory of 53 swells. Estimates from magmatic hotspots yield a cumulative lower bound on global plume flux of 2 TW, which increases to 6 TW if amagmatic swells are also included and if all buoyancy is assumed to be thermal in origin. Our results suggest that upwelling plumes play a significant role in the transfer of heat into the uppermost mantle.
Physics of the Earth and Planetary Interiors, Vol. 306, 106509, 18p. Pdf
Mantle
geothermometry
Abstract: Thermal structure of the lithosphere exerts a primary control on its strength and density and thereby its dynamic evolution as the outer thermal and mechanic boundary layer of the convecting mantle. This contribution focuses on continental lithosphere. We review constraints on thermal conductivity and heat production, geophysical and geochemical/petrological constraints on thermal structure of the continental lithosphere, as well as steady-state and non-steady state 1D thermal models and their applicability. Commonly used geotherm families that assume that crustal heat production contributes an approximately constant fraction of 25-40% to surface heat flow reproduce the global spread of temperatures and thermal thicknesses of the lithosphere below continents. However, we find that global variations in seismic thickness of continental lithosphere and seismically estimated variations in Moho temperature below the US are more compatible with models where upper crustal heat production is 2-3 times higher than lower crustal heat production (consistent with rock estimates) and the contribution of effective crustal heat production to thermal structure (i.e. estimated by describing thermal structure with steady-state geotherms) varies systematically from 40 to 60% in tectonically stable low surface heat flow regions to 20% or lower in higher heat flow tectonically active regions. The low effective heat production in tectonically active regions is likely partly the expression of a non-steady thermal state and advective heat transport.
Abstract: The thermobarometric analysis of inclusions in lithospheric diamonds indicates that they originated from a wide range of depths, with a global mode at ca. 170±15 km [1]. Studies based on diamond depth distribution at global scale, however, cannot clarify if this mode reflects a real concentration of diamonds, preferential sampling of materials from this level by rising kimberlites, or even a statistical distribution within the hard limits imposed by diamond stability, lithosphere thickness, and mantle adiabat under typical cratonic thermal regimes. We addressed this problem by comparing depth distributions for peridotitic diamonds from the three localities that have been the most prolific for diamond geobarometry (Cullinan, Kimberley and Voorspoed, South Africa) with those of mantle xenocrysts from the same kimberlite sources. P-T estimates indicate that the diamonds were formed at T higher, equal or lower than the ambient geotherm. They may record old mantle thermal regimes or local thermal perturbations related to infiltration of parent fluids or melts. Nonetheless, the diamonds show similar depth distributions for different localities, with a distinct mode at ?175 ?? 10 km. The similarity of these distributions with that calculated for peridotitic diamonds worldwide, as well as the lack of systematic correlation with kimberlite sampling efficiency as recorded by mantle xenocrysts, suggests that this mode has genetic significance. Based on observed depth distributions at both local and global scale and on thermodynamic modeling of COH fluids, diamond-forming processes are predicted to become less efficient with decreasing depth from at least ?160 km. In addition, diamond endowment near the base of the lithosphere may be negatively affected by infiltration of carbon-undersaturated melts. Considering the poor correlation between diamond and xenocryst depth distributions in single kimberlites or kimberlite clusters, even limited xenocryst records from diamond favorable depths (especially the 160-190 km interval) may correspond to significant diamond potential.
Physics of the Earth and Planetary Interiors, Vol. 305, 106457, 17p. Pdf
Mantle
geothermometry
Abstract: The role of heat flow coming from the core is often overlooked or underestimated in simple models of Earth's thermal evolution. Throughout most of Earth's history, the mantle must have been extracting from the core at least the amount of heat that is required to operate the geodynamo. In view of recent laboratory measurements and theoretical calculations indicating a higher thermal conductivity of iron than previously thought, the above constraint has important implications for the thermal history of the Earth's mantle. In this work we construct a paramaterized mantle convection model that treats both the top and the core-mantle boundary heat fluxes according to the boundary layer theory, or alternatively employs the model of Labrosse (2015) to compute the thermal evolution of the Earth's core. We show that the core is likely to provide all the missing heat that is necessary in order to avoid the so-called “thermal catastrophe” of the mantle. Moreover, by analyzing the mutual feedback between the core and the mantle, we provide the necessary ingredients for obtaining thermal histories that are consistent with the petrological record and have reasonable initial conditions. These include a sufficiently high viscosity contrast between the lower and upper mantle, whose exact value is sensitive to the activation energy that governs the temperature dependence of the viscosity.
Abstract: In the two decades since Subduction: Top to Bottom was published in 1996, improved analytical and numerical thermal-petrologic models of subduction zones have been constructed and evaluated against new seismological and geological observations. Advances in thermal modeling include a range of new approaches to incorporating shear (frictional, viscous) heating along the subduction interface and to simulating induced flow in the mantle wedge. Forearc heat-flux measurements constrain the apparent coefficient of friction (??) along the plate interface to 0.1, but the extent to which ?? may vary between subduction zones remains challenging to discern owing to scatter in the heat-flux measurements and uncertainties in the magnitude and distribution of radiogenic heat production in the overriding crust. Flow in the mantle wedge and the resulting thermal structure depend on the rheology of variably hydrated mantle rocks and the depth at which the subducting slab becomes coupled to the overlying mantle wedge. Advances in petrologic modeling include the incorporation of sophisticated thermodynamic software packages into thermal models and the prediction of seismic velocities from mineralogic and petrologic models. Current thermal-petrologic models show very good agreement between the predicted location of metamorphic dehydration reactions and observed intermediate-depth earthquakes, and between the predicted location of the basalt-to-eclogite transition in subducting oceanic crust and observed landward-dipping, low-seismic-velocity layers. Exhumed high-pressure, low-temperature metamorphic rocks provide insight into subduction-zone temperatures, but important thermal parameters (e.g., convergence rate) are not well constrained, and metamorphic rocks exposed at the surface today may reflect relatively warm conditions in the past associated with subduction initiation or ridge subduction. We can anticipate additional advances in our understanding of subduction zones as a result of further testing of model predictions against geologic and geophysical observations, and of evaluating the importance of advective processes, such as diapirism and subduction-channel flow, that are not captured in hybrid kinematic-dynamic models of subduction zones but are observed in fully dynamical models under certain conditions.
Journal of the Geological Society, Vol. 177, pp. 784-798.
Africa, Madagascar
geothermometry
Abstract: Madagascar occupied an important place in the amalgamation of Gondwana and preserves a record of several Neoproterozoic events that are linked to orogenesis of the East African Orogen. In this study, we integrate remote sensing, field data and thermochronology to unravel complex deformation in the Ikalamavony and Itremo domains of central Madagascar. The deformation sequence comprises a gneissic foliation (S1), followed by south- to SW-directed, tight to isoclinal, recumbent folding (D2). These are overprinted by north-trending upright folds that formed during an approximately east-west shortening event (D3). Together these produced type 1 and type 2 fold interference patterns throughout the Itremo and Ikalamavony domains. We show that the Itremo and Ikalamavony domains were deformed together in the same orogenic system, which we interpret as the c. 630 Ma collision of Azania with Africa along the Vohibory Suture in southwestern Madagascar. In eastern Madagascar, deformation is syn- to post-550 Ma, and probably formed in response to final closure of the Mozambique Ocean along the Betsimisaraka Suture that amalgamated Madagascar with the Dharwar Craton of India. Apatite U-Pb and novel laser ablation triple quadrupole inductively coupled plasma mass spectrometry (LA-QQQ-ICP-MS) muscovite and biotite Rb-Sr thermochronology indicates that much of central Madagascar cooled through c. 500°C at c. 500 Ma.
Vestnik Otdelenia nauk o Zemle RAN *** ENG, Vol. 3, doi:10.2205/2011NZ000138, 5p. Pdf * note date
Mantle
geobarometry, peridotites
Abstract: Original monomineral thermobarometers for mantle peridotites for clinopyroxene, garnet, chromite and ilmenites for the mantle peridotites were statistically calibrated on the PT estimates for mantle peridotites [Ashchepkov et al., 2010] were tested using the mineral phases obtained in high pressure experiments with the natural peridotites (380 runs) [Brey et al.,1990; 2008 etc] and eclogites (240 runs)[Dasgupta et al., 2006 etc]. In the original program of that written on FORTRAN are assembled the most reliable methods of mineral thermometers (45) and barometers (36) and oxybarometers (5), including original monomineral and methods [Ashchepkov, 2003 Ashchepkov et al., 2008; 2009; 2010; 2011] for the mantle peridotites bases on the compositions of on clinopyroxene, garnet, chromite and ilmenite. Program reads the text files, which converted from Excel. Original data include standard silicate compositions for 12 components in standard order. The text file includes 15 columns of 8 symbols. The first is file name which is the same for all the minerals in the association. The second is indicator symbol for phases. E- enstatite, D - diopside, O-olivine, S-spinel, G- garnet, I -ilmenite, A- amphibole, F - phlogopite, P-plagioclase, L- liquid, R- bulk rock. Then follow oxides: SiO2, TiO2, Al2O3, Cr2O3, FeO, MnO, MgO, CaO, Na2O, K2O, NiO, V2O3 written with 2-3 decimals. The last column may contain description of the mineral or association up to 64 symbols. Monomineral methods use calculated values for Fe#Ol or Fe#Cpx. The input from console includes file name (8 symbols) (A8), then amount of PT pairs of numbers thermometers and barometers (2I2) and one for FO2 method. Program allows input of the iteration numbers (to 25 by default). It allow to choose whether to use the calculated Fe3+ for the minerals and also. It is possible also to put fixed values of T and P (default 1000o C and 40 kbar ).
Eclogitic geotherms of the Rio de la Plata craton archon-core: Estancia Trementina and Puentesino, Dpto. Of Concepccion-Paraguay. Compared to two large diamond deposits Argyle ( lamproitic ) and Orapa ( kimberlitic.
Abstract: QUIDDIT is a free Python software-package designed to process Fourier Transform Infrared (FTIR) spectra of diamonds automatically and efficiently. Core capabilities include baseline correction, determination of nitrogen concentration, nitrogen aggregation state and model temperature and fitting of both the 3107 cm-1 and platelet (B’) peaks. These capabilities have allowed the authors to study platelet defects and their relationship to nitrogen aggregation in previous studies. Data visualisation, vital to interpreting and evaluating results, is another key component of the software. QUIDDIT can be applied to single spectra as well as linescan and 2-dimensional map data. Recently, additional features such as manual platelet peak and nitrogen fitting, custom batch peak fitting and two-stage aggregation modelling were made available. QUIDDIT has been used successfully for natural diamonds containing aggregated forms of nitrogen in the past and has since been adapted for the study of diamonds containing C-centres as well.
Abstract: We have performed an analysis of the cases of synchronism in th egrowth temperature in local zones of diamond crystals and the concentration of hydrogen in them.The considered cases were observed by the authors and fined out in the iterature. Possible causes of the simbatic change in the crystal growth temperature and the concentration of hydrogen in it are considered.The determination of the temperature change over the zones was carried out on the basis of local FTIR spectroscopy from the ratio of the nitrogen concentration in the form of defects in the crystal structure of A and B1, and size the B2 defects.The change in the hydrogen concentration in various zones of diamond crystals was estimated from the 3107cm-1 band of the hydrogen-containing defect. It is shown that in the analyzed cases the concentration of hydrogen in diamond is determined mainly by its content in the growth medium.We accept the obtained results as evidence of the participation of hydrogen in the heat transfer in mantle mineral-forming systems.
Progress in Earth and Planetary Science, doi.org./10.1186/ s40645-020-00374-8 17p. Pdf
Mantle
geophysics, geothermometry
Abstract: In this review, I provide the current status and future prospects for the coupled core-mantle evolution and specifically summarize the constraints arising from geomagnetism and paleomagnetism on the long-term secular variations of the geomagnetic field. The heat flow across the core-mantle boundary (CMB) is essential for determining the best-fit scenario that explains the observational data of geomagnetic secular variations (e.g., onset timing of the inner core growth, geomagnetic polarity reversals, and westward drift) and should include the various origins of the heterogeneous structures in the deep mantle that have affected the heat transfer across the core-mantle boundary for billions of years. The coupled core-mantle evolution model can potentially explain the onset timing of the inner core and its influence on the long-term geomagnetic secular variations, but it is still controversial among modeling approaches on the core energetics because the paleomagnetic data contains various uncertainties. Additionally, with the coupled core-mantle evolution model in geodynamo simulations, the frequency of the geomagnetic polarity reversals can be explained with the time variations of the heat flow across the CMB. Additionally, the effects of the stable region in the outermost outer core to the magnetic evolution are also crucial but there would be still uncertain for their feasibility. However, despite this progress in understanding the observational data for geomagnetic secular variations, there are several unresolved issues that should be addressed in future investigations: (1) initial conditions—starting with the solidification of the global magma ocean with the onset timing of plate tectonics and geodynamo actions and (2) planetary habitability—how the dynamics of the Earth’s deep interior affects the long-term surface environment change that has been maintained in the Earth’s multisphere coupled system.
Earth and Planetary Letters, Vol. 550, 116549, 13p.
Global, United States, Wyoming, Canada, Northwest Territories, Europe, Baltic, India
geothermometry
Abstract: Cratonic lithosphere is believed to have been chemically buoyant and mechanically resistant to destruction over billions of years. Yet the absence of cratonic roots at some Archean terrains casts doubt on the craton stability and longevity on a global scale. As unique mantle-derived melts at ancient continents, silica-poor, kimberlitic melts are ideal tools to constrain the temporal variation of lithosphere thickness and the processes affecting the lithosphere root. However, no reliable thermobarometer exists to date for strongly silica-undersaturated, mantle-derived melts. Here we develop a new thermobarometer for silica-poor, CO2-rich melts using high-temperature, high-pressure experimental data. Our barometer is calibrated based on a new observation of pressure-dependent variation of Al2O3 in partial melts saturated with garnet and olivine, while our thermometer is calibrated based on the well-known olivine-melt Mg-exchange. For applications to natural magmas, we also establish a correction scheme to estimate their primary melt compositions. Applying this liquid-based thermobarometer to the estimated primary melt compositions for a global kimberlite dataset, we show that the equilibration depths between primary kimberlite melts and mantle peridotites indicate a decrease of up to ?150 km in cratonic lithosphere thickness globally during the past ?2 Gyr. Together with the temporal coupling between global kimberlite frequency and cold subduction flux since ?2 Gyr ago, our results imply a causal link between lithosphere thinning and supply of CO2-rich melts enhanced by deep subduction of carbonated oceanic crusts. While hibernating at the lithosphere root, these melts chemically metasomatize and rheologically weaken the rigid lithosphere and consequently facilitate destruction through convective removal in the ambient mantle or thermo-magmatic erosion during mantle plume activities.
Abstract: The mantle beneath the Western Dharwar Craton of the Indian shield comprises a suite of refractory and fertile peridotites and mafic granulites. Detailed petrographic studies coupled with new mineral analysis and geothermobarometric estimations permit to decipher the thermal architecture and get an insight into the evolution of this ancient craton. The refractory rocks are coarse grained harzburgites/dunites, whereas the more fertile ones are at times, porphyroclastic lherzolites. Both show a similar range of equilibration temperatures and pressures indicating intermixing between the two at various levels. The peridotites contain undeformed interstitial REE-enriched clinopyroxene, phlogopite, apatite and carbonates recording post-kinematic modal and cryptic metasomatic events in the Precambrian cratonic lithosphere. Xenoliths of mafic granulite contain layers of clinopyroxenite which also vein the granulite. The P-T range of the granulites overlaps that of the ultramafic rocks. This study in combination with previous investigations reveals a distinct change in the thermal architecture of the craton from a warm/hot geotherm in the Proterozoic to a highly perturbed, still hotter geotherm of the Palaeocene. The Cenozoic thermotectonic rifting episodes heated, refertilized and thinned the bulk of the cratonic lithosphere beneath the Western Dharwar Craton, which has witnessed the most re-activation among cratons of the Indian shield. The waning of the Deccan Traps volcanism in Palaeocene time saw the reworking of ancient cratonic lithosphere and its replacement by non-cratonic, juvenile mantle and magmatic accretions, indicated by compound xenoliths. Differing petrological and geochemical characteristics of refractory xenoliths and fertile lherzolites serve to constrain the relative timing and composition of non-cratonic lithosphere. By the end of the Palaeocene the Western Dharwar Craton was characterised by a thermal high, an attenuated continental lithosphere (60-80 km), and a thin crust (<10- ~ 21 km), reflecting the decratonization of at least the western part of the Western Dharwar Craton.
Abstract: At the Rio de la Plata Craton archon-core environment were inferred, based on 1D Vs profiles (on 208 numbers of points), of the peridotitic geotherms. Values for the archon-core environment, it was estimated 38.5 to 40 mW/m2 in its central northern portion and southern portion and in its edges/southern portion 40 to 42 mW/m2. Geotherm values that allowed estimate LAB between 243 to 237 km depth (northern portion) and 225 to 213 km depth (southern portion). The same 1D Vs information allowed recognizing for this geothermal environment the depth of the graphite-to-diamond phase transition, finding that it is located at ~135 km. depth. So, projecting 70-90 Km. (southern portion) to 102-108 km. (northern portion) thickness of the “diamond window” for the Rio de la Plata craton archon-core. "Diamond window" thickness very close to those of the Kalahari archon craton where the highest grade of diamond deposit is the Kimberley with 200 cpht. Thus, it is estimated for eventual diamond deposit, in the Río de la Plata craton core, are quite similar to Kimberley diamond deposits could be also expected in the archon-core of Río de la Plata craton.
Eclogitic geotherms of the Rio de la Plata craton archon-core: Estancia Trementina and Puentesino, Dpto. Of Concepion - Paraguay. Compared to two large diamond deposits Argyle ( lamproitic) and Orapa ( Kimberlitic).
Contributions to Mineralogy and Petrology, dor.org/10.1007/ s00410-020-01750-9 27p. Pdf
Mantle
Geothermometry
Abstract: The interplay between stress and chemical processes is a fundamental aspect of how rocks evolve, relevant for understanding fracturing due to metamorphic volume change, deformation by pressure solution and diffusion creep, and the effects of stress on mineral reactions in crust and mantle. There is no agreed microscale theory for how stress and chemistry interact, so here I review support from eight different types of the experiment for a relationship between stress and chemistry which is specific to individual interfaces: (chemical potential)?=?(Helmholtz free energy)?+?(normal stress at interface)?×?(molar volume). The experiments encompass temperatures from -100 to 1300 degrees C and pressures from 1 bar to 1.8 GPa. The equation applies to boundaries with fluid and to incoherent solid-solid boundaries. It is broadly in accord with experiments that describe the behaviours of free and stressed crystal faces next to solutions, that document flow laws for pressure solution and diffusion creep, that address polymorphic transformations under stress, and that investigate volume changes in solid-state reactions. The accord is not in all cases quantitative, but the equation is still used to assist the explanation. An implication is that the chemical potential varies depending on the interface, so there is no unique driving force for reaction in stressed systems. Instead, the overall evolution will be determined by combinations of reaction pathways and kinetic factors. The equation described here should be a foundation for grain-scale models, which are a prerequisite for predicting larger scale Earth behaviour when stress and chemical processes interact. It is relevant for all depths in the Earth from the uppermost crust (pressure solution in basin compaction, creep on faults), reactive fluid flow systems (serpentinisation), the deeper crust (orogenic metamorphism), the upper mantle (diffusion creep), the transition zone (phase changes in stressed subducting slabs) to the lower mantle and core mantle boundary (diffusion creep).
Abstract: The Earth’s cratons are traditionally regarded as tectonically stable cores that were episodically buried by thin sedimentary covers. Cratonic crust in southern Finland holds seven post-1.7 Ga tiered unconformities, with remnants of former sedimentary covers. We use the geometries of the tiered unconformities, along with previously dated impact structures and kimberlite and carbonatite pipes, to reconstruct the erosion and burial history of the craton and to derive estimates of depths of erosion in basement and former sedimentary rocks. The close vertical spacing (<200 m) of the unconformities and the survival of small (D ? 5 km) Neoproterozoic and Early Palaeozoic impact structures indicate minor later erosion. Average erosion rates (<2.5 m/Ma) in basement and cover are amongst the lowest reported on Earth. Ultra-slow erosion has allowed the persistence in basement fractures of Phanerozoic fracture coatings and Palaeogene groundwater and microbiomes. Maximum thicknesses of foreland basin sediments in Finland during the Sveconorwegian and Caledonide orogenies are estimated as ~1.0 km and <0.68-1.0 km, respectively. Estimated losses of sedimentary cover derived from apatite fission track thermochronology are higher by factors of at least 2 to 4. A dynamic epeirogenic history of the craton in Finland, with kilometre-scale burial and exhumation, proposed in recent thermochronological models is not supported by other geological proxies. Ultra-slow erosion rates in southern Finland reflect long term tectonic stability and burial of the craton surface for a total of ~1.0 Ga beneath generally thin sedimentary cover.
Contributions to Mineralogy and Petrology, Vol. 176, 10.1007/s0041 0-020-01768-z 21p. Pdf
Mantle
geothermometry
Abstract: The pressure dependence of the exchange of Cr between clinopyroxene and garnet in peridotite is applicable as a geobarometer for mantle-derived Cr-diopside xenocrysts and xenoliths. The most widely used calibration (Nimis and Taylor Contrib Miner Petrol 139: 541-554, 2000; herein NT00) performs well at pressures below 4.5 GPa, but has been shown to consistently underestimate pressures above 4.5 GPa. We have experimentally re-examined this exchange reaction over an extended pressure, temperature, and compositional range using multi-anvil, belt, and piston cylinder apparatuses. Twenty-nine experiments were completed between 3-7 GPa, and 1100-1400 °C in a variety of compositionally complex lherzolitic systems. These experiments are used in conjunction with several published experimental datasets to present a modified calibration of the widely-used NT00 Cr-in-clinopyroxene (Cr-in-cpx) single crystal geobarometer. Our updated calibration calculates P (GPa) as a function of T (K), CaCr Tschermak activity in clinopyroxene (acpxCaCrTs), and Cr/(Cr?+?Al) (Cr#) in clinopyroxene. Rearranging experimental results into a 2n polynomial using multiple linear regression found the following expression for pressure: P(GPa)=11.03+(?T(K) ln(acpxCaCrTs)×0.001088)+(1.526×ln(Cr#cpxT(K))) where Cr#cpx=(CrCr+Al), acpxCaCrTs=Cr?0.81?Cr#cpx?(Na+K), with all mineral components calculated assuming six oxygen anions per formula unit in clinopyroxene. Temperature (K) may be calculated through a variety of geothermometers, however, we recommend the NT00 single crystal, enstatite-in-clinopyroxene (en-in-cpx) geothermometer. The pressure uncertainty of our updated calibration has been propagated by incorporating all analytical and experimental uncertainties. We have found that pressure estimates below 4 GPa, between 4-6 GPa and above 6 GPa have associated uncertainties of 0.31, 0.35, and 0.41 GPa, respectively. Pressures calculated using our calibration of the Cr-in-cpx geobarometer are in good agreement between 2-7 GPa, and 900-1400 °C with those estimated from widely-used two-phase geobarometers based on the solubility of alumina in orthopyroxene coexisting with garnet. Application of our updated calibration to suites of well-equilibrated garnet lherzolite and garnet pyroxenite xenoliths and xenocrysts from the Diavik-Ekati kimberlite and the Argyle lamproite pipes confirm the accuracy and precision of our modified geobarometer, and show that PT estimates using our revised geobarometer result in systematically steeper paleogeotherms and higher estimates of the lithosphere?asthenosphere boundary compared with the original NT00 calibration.
Contributions to Mineralogy and Petrology, 176, 16p. Pdf
Mantle
geothermobarometry
Abstract: The temperature-dependent exchange of Ni and Mg between garnet and olivine in mantle peridotite is an important geothermometer for determining temperature variations in the upper mantle and the diamond potential of kimberlites. Existing calibrations of the Ni-in-garnet geothermometer show considerable differences in estimated temperature above and below 1100 °C hindering its confident application. In this study, we present the results from new synthesis experiments conducted on a piston cylinder apparatus at 2.25-4.5 GPa and 1100-1325 °C. Our experimental approach was to equilibrate a Ni-free Cr-pyrope-rich garnet starting mixture made from sintered oxides with natural olivine capsules (Niolv ? 3000 ppm) to produce an experimental charge comprised entirely of peridotitic pyrope garnet with trace abundances of Ni (10-100 s of ppm). Experimental runs products were analysed by wave-length dispersive electron probe microanalysis (EPMA) and laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS). We use the partition coefficient for the distribution of Ni between our garnet experimental charge and the olivine capsule (lnDNigrt/olv;NigrtNiolv), the Ca mole fraction in garnet (XCagrt; Ca/(Ca?+?Fe?+?Mg)), and the Cr mole fraction in garnet (XCrgrt; Cr/(Cr?+?Al)) to develop a new formulation of the Ni-in-garnet geothermometer that performs more reliably on experimental and natural datasets than existing calibrations. Our updated Ni-in-garnet geothermometer is defined here as: T(?C)=?8254.568((XCagrt×3.023)+(XCrgrt×2.307)+(lnDNigrtolv?2.639))?273±55 where DNigrt/olv=NigrtNiolv, Ni is in ppm, XCagrt = Ca/(Ca?+?Fe?+?Mg) in garnet, and XCrgrt= Cr/(Cr?+?Al) in garnet. Our updated Ni-in-garnet geothermometer can be applied to garnet peridotite xenoliths or monomineralic garnet xenocrysts derived from disaggregation of a peridotite source. Our calibration can be used as a single grain geothermometer by assuming an average mantle olivine Ni concentration of 3000 ppm. To maximise the reliability of temperature estimates made from our Ni-in-garnet geothermometer, we provide users with a data quality protocol method which can be applied to all garnet EPMA and LA-ICP-MS analyses prior to Ni-in-garnet geothermometry. The temperature uncertainty of our updated calibration has been rigorously propagated by incorporating all analytical and experimental uncertainties. We have found that our Ni-in-garnet temperature estimates have a maximum associated uncertainty of ± 55 °C. The improved performance of our updated calibration is demonstrated through its application to previously published experimental datasets and on natural, well-characterised garnet peridotite xenoliths from a variety of published datasets, including the diamondiferous Diavik and Ekati kimberlite pipes from the Lac de Gras kimberlite field, Canada. Our new calibration better aligns temperature estimates using the Ni-in-garnet geothermometer with those estimated by the widely used (Nimis and Taylor, Contrib Mineral Petrol 139:541-554, 2000) enstatite-in-clinopyroxene geothermometer, and confirms an improvement in performance of the new calibration relative to existing versions of the Ni-in-garnet geothermometer.
Abstract: Pressure and temperature estimates of rocks provide the fundamental data for the investigation of many geological processes such as subduction and exhumation, and yet their determination remains extremely challenging (Tajcmanova et al. 2020). A wide variety of methods are constantly being developed to tackle the ambitious objective of pinpointing the geological history of rocks through the many complex processes often interacting with one another at depth in our planet. Analytical advances are being pushed to the limit of conventional methods, allowing information preserved by mineral, fluid, and solid inclusions to be used for high spatial resolution determinations that can unravel a large variety of processes occurring at the micro- to the nano-scale. Among these, chemical geothermobarometry that is often challenging in many rock types due to alteration processes, chemical re-equilibration, diffusion, and kinetic limitations has been increasingly coupled with elastic geothermobarometry (e.g., Anzolini et al. 2019; Gonzalez et al. 2019). Elastic geothermobarometry of host-inclusion systems, in paper Mazzucchelli et al. 2021, this issue, is a new and complementary non-destructive method (see Fig. 1 for an example) to determine the pressures (P) and temperatures (T) of inclusion entrapment (i.e., the P-T conditions attained by rocks and minerals at depth in the Earth) from the remnant stress or strain measured in inclusions still trapped in their host mineral at room conditions (e.g., Nestola et al. 2011; Howell et al. 2012; Alvaro et al. 2020).
Abstract: The lithosphere and upper mantle of South America is investigated using multiple data sets, including the topography, crustal structure, regional seismic tomography, gravity, and mineral physics. These data are jointly inverted to estimate variations in temperature, density and composition in the lithospheric and sub-lithospheric upper mantle to a depth of 325 km. Our results show significant variations in lithospheric properties, including thick, depleted roots beneath large parts of the Amazon, São Francisco, and Paranapanema Cratons. However, portions of some cratons, such as the western Guyana Shield, lack a depleted root. We hypothesize that these regions either never developed a depleted root, or that the root was rejuvenated by lithospheric processes.
Geodynamics & Tectonophysics, Vol. 11, pp. 75-87. pdf
Russia, Siberia
geothermometry
Abstract: We present the first results of fission-track dating of apatite monofractions from two rock samples taken from the Southern carbonatite massif of the world’s largest alkaline ultrabasic Guli pluton (~250 Ma), located within the Maymecha-Kotuy region of the Siberain Traps. Based on the apatite fission-track data and computer modeling, we propose two alternative model of the Guli pluton's tectonothermal history. The models suggest (1) rapid post-magmatic cooling of the studied rocks in hypabyssal conditions at depth about 1.5 km, or (2) their burial under a 2-3 km thick volcano-sedimentary cover and reheating above 110°C, followed by uplift and exhumation ca. 218 Ma.
Abstract: Earthquakes occurring below ?300 km, especially in the mantle transition zone are some of the strongest events experienced on Earth. Deep earthquakes, whose nature and cause are poorly known, occur with regularity and are a deep and prominent result of plate tectonics. We model the paths of subducting slabs to relate pressure-temperature conditions to the experimentally determined mineralogies of the slab crust and mantle. We present a synthesis of mantle minerals included in diamonds derived from same depths as the deep earthquakes to show that fluids exist there. We show that decarbonization/melting reactions in the slab crust and dehydration reactions in the slab mantle can provide fluids to the earthquake generation regions, suggesting that fluids cause or are related to deep earthquakes.
Abstract: A colour-changing garnet exhibits the "alexandrite effect", whereby its colour changes from green in the presence of daylight to purplish red under incandescent light. This study examines this species of garnets as well as the causes of the colour change by using infrared and ultraviolet visible (UV-Vis) spectroscopy. The infrared spectra show that the colour-changing garnets in this paper belong to the solid solution of pyrope-spessartine type. CIE1931 XYZ colour matching functions are used to calculate the colour parameters influencing garnet colour-changing under different light sources. The UV-Vis spectra show two zones of transmittance, in the red region at 650-700 nm and the blue-green region at 460-510 nm. As they exhibit the same capacity to transmit light, the colour of the gem is determined by the external light source. The absorption bands of Cr3+ and V3+ at 574 nm in the UV-Vis spectra are the main cause of the change in colour. With the increase in the area of peak absorption, the differences in the chroma and colour of the garnet gradually increase in daylight and incandescent light, and it exhibits a more prominent colour-changing effect.
Abstract: Crystalline basement rocks of southwestern Montana have been subjected to multiple tectonothermal events since ?3.3 Ga: the Paleoproterozoic Big Sky/Great Falls orogeny, Mesoproterozoic extension associated with Belt-Purcell basin formation, Neoproterozoic extension related to Rodinia rifting, and the late Phanerozoic Sevier-Laramide orogeny. We investigated the long-term (>1 Ga), low-temperature (erosion/burial within 10 km of the surface) thermal histories of these tectonic events with zircon and apatite (U-Th)/He thermochronology. Data were collected across nine sample localities (n = 55 zircon and n = 26 apatite aliquots) in the northern and southern Madison ranges, the Blacktail-Snowcrest arch, and the Tobacco Root uplift. Our zircon (U-Th)/He data show negative trends between single aliquot date and effective uranium (a radiation damage proxy), which we interpreted with a thermal history model that considers the damage-He diffusivity relationship in zircon. Our model results for these basement ranges show substantial cooling from temperatures above 400°C to near surface conditions between 800 and 510 Ma. Subsequent Phanerozoic exhumation culminated by ?75 Ma. Late Phanerozoic cooling is coincident with along-strike Sevier belt thin-skinned thrusting in southeastern Idaho, and older than exhumation in basement-involved uplifts of the Wyoming Laramide province. Our long-term, low-temperature thermal record for these southwestern Montana basement ranges shows that: (a) these basement blocks have experienced multiple episodes of upper crustal exhumation and burial since Archean time, possibly influencing Phanerozoic thrust architecture and (b) the late Phanerozoic thick-skinned thrusting recorded by these rocks is among the earliest thermochronologic records of Laramide basement-involved shortening and was concomitant with Sevier belt thin-skinned thrusting.
Abstract: Petrological data indicate that upper mantle and mantle plume temperatures diverged 2.5 billion years ago. This has been interpreted as plate tectonics initiating at 2.5 Ga with Earth operating as a single plate planet before then. We take an Occam’s razor view that the continuous operation of plate tectonics can explain the divergence. We validate this hypothesis by comparing petrological data to results from mixed heating mantle convection models in a plate tectonic mode of mantle cooling. The comparison shows that the data are consistent with plate tectonics operating over geologic history.