Abstract Raman spectroscopy is widely used to identify mineral and fluid inclusions in host crystals, as well as to calculate pressure-temperature (P-T) conditions with mineral inclusion elastic thermobarometry, for example quartz-in-garnet barometry (QuiG) and zircon-in-garnet thermometry (ZiG). For thermobarometric applications, P-T precision and accuracy depend crucially on the reproducibility of Raman peak position measurements. In this study, we monitored long-term instrument stability and varied analytical parameters to quantify peak position reproducibility for Raman spectra from quartz and zircon inclusions and reference crystals. Our ultimate goal was to determine the reproducibility of calculated inclusion pressures (“Pinc”) and entrapment pressures (“Ptrap”) or temperatures (“Ttrap”) by quantifying diverse analytical errors, as well as to identify optimal measurement conditions and provide a baseline for interlaboratory comparisons. Most tests emphasized 442 nm (blue) and 532 nm (green) laser sources, although repeated analysis of a quartz inclusion in garnet additionally used a 632.8 nm (red) laser. Power density was varied from <1 to >100 mW and acquisition time from 3 to 270s. A correction is proposed to suppress interference on the ~206 cm–1 peak in quartz spectra by a broad nearby (~220 cm–1) peak in garnet spectra. Rapid peak drift up to 1 cm–1/h occurred after powering the laser source, followed by minimal drift (<0.2 cm–1/h) for several hours thereafter. However, abrupt shifts in peak positions as large as 2–3 cm–1 sometimes occurred within periods of minutes, commonly either positively or negatively correlated to changes in room temperature. An external Hg-emission line (fluorescent light) can be observed in spectra collected with the green laser and shows highly correlated but attenuated directional shifts compared to quartz and zircon peaks. Varying power density and acquisition time did not affect Raman peak positions of either quartz or zircon grains, possibly because power densities at the levels of inclusions were low. However, some zircon inclusions were damaged at higher power levels of the blue laser source, likely because of laser-induced heating. Using a combination of 1, 2, or 3 peak positions for the ~128, ~206, and ~464 cm–1 peaks in quartz to calculate Pinc and Ptrap showed that use of the blue laser source results in the most reproducible Ptrap values for all methods (0.59 to 0.68 GPa at an assumed temperature of 450 °C), with precisions for a single method as small as ±0.03 GPa (2σ). Using the green and red lasers, some methods of calculating Ptrap produce nearly identical estimates as the blue laser with similarly good precision (±0.02 GPa for green laser, ±0.03 GPa for red laser). However, using 1- and 2-peak methods to calculate Ptrap can yield values that range from 0.52 ± 0.06 to 0.93 ± 0.16 GPa for the green laser, and 0.53 ± 0.08 GPa to 1.00 ± 0.45 GPa for the red laser. Semiquantitative calculations for zircon, assuming a typical error of ±0.25 cm–1 in the position of the ~1008 cm–1 peak, imply reproducibility in temperature (at an assumed pressure) of approximately ±65 °C. For optimal applications to elastic thermobarometry, analysts should: (1) delay data collection approximately one hour after laser startup, or leave lasers on; (2) collect a Hg-emission line simultaneously with Raman spectra when using a green laser to correct for externally induced shifts in peak positions; (3) correct for garnet interference on the quartz 206 cm–1 peak; and either (4a) use a short wavelength (blue) laser for quartz and zircon crystals for P-T calculations, but use very low-laser power (<12 mW) to avoid overheating and damage or (4b) use either the intermediate wavelength (green; quartz and zircon) or long wavelength (red; zircon) laser for P-T calculations, but restrict calculations to specific methods. Implementation of our recommendations should optimize reproducibility for elastic geothermobarometry, especially QuiG barometry and ZiG thermometry.
more »
« less
Regional Quartz Inclusion Barometry and Comparison with Conventional Thermobarometry and Intersecting Isopleths from the Connecticut Valley Trough, Vermont and Massachusetts, USA
Abstract A comparative analysis of Raman shifts of quartz inclusions in garnet was made along two traverses across the Connecticut Valley Trough (CVT) in western New England, USA, to examine the regional trends of quartz inclusion in garnet (QuiG) Raman barometry pressure results and to compare this method with conventional thermobarometry and the method of intersecting garnet core isopleths. Overall, Raman shifts of quartz inclusions ranged from 1·2 to 3·5 cm–1 over all field areas and displayed a south to north decrease, matching the overall decrease in mapped metamorphic grade. Raman shifts of quartz inclusions typically did not show systematic variation with respect to their radial position within a garnet crystal, and indicate that garnet probably grew at nearly isothermal and isobaric pressure–temperature (P–T) conditions. The P–T conditions inferred from conventional thermobarometry were in the range of ∼500–575 °C and ∼7·4–10·3 kbar over the sample suite and are in good agreement with previous published thermobarometry throughout the CVT. These P–T results are broadly consistent with QuiG barometry and also suggest that garnet grew isothermally and isobarically at near peak P–T conditions. However, P–T conditions and P–T paths inferred using either garnet core thermobarometry or garnet core intersecting isopleths yield results that are internally inconsistent and generally disagree with the pressure results from QuiG barometry. Garnet core isopleth intersections consistently plotted between the nominal garnet-in curve on mineral assemblage diagrams and the P–T conditions constrained by QuiG isomekes for the majority of the sample suite. Additionally, most samples’ P–T results from QuiG barometry and rim thermobarometry show marked disagreement from those derived from garnet core thermobarometry, compared with the minority that showed agreement within uncertainty. Pressures calculated from QuiG barometry ranged from 8·5 to 9·5 kbar along the traverses in western Massachusetts (MA) and central Vermont (VT) and from 6·5 to 7·5 kbar in northern VT indicating an increase in peak burial of 3–6 km from north to south. Along the western end of the central VT traverse, there are differences in measured Raman shifts and inferred peak pressures of up to 1 kbar across the Richardson Memorial Contact (RMC), indicating a possible fault contact with minor post-peak metamorphic shortening of up to ∼3 km. In contrast, along an east–west traverse in the vicinity of the Goshen Dome, MA, there was little observed variation in Raman shifts across the contact. By contrast, QuiG barometry clearly indicates significant discontinuities in peak pressure east of the Strafford Dome in central VT. This supports the interpretation that post-peak metamorphic shortening was necessary to juxtapose upper staurolite–kyanite zone rocks next to lower garnet zone pelites. Overall, it is concluded that garnet core thermobarometry and garnet core isopleths may provide unreliable results for the P–T conditions of garnet nucleation and inferred P–T paths during garnet growth unless independently verified. The consistency of QuiG results with rim thermobarometry indicates that peak metamorphic conditions previously reported for the CVT using garnet rim thermobarometry are robust and that variation in QuiG barometry results is a valuable tool to analyze structural features within a metamorphic terrane.
more »
« less
- NSF-PAR ID:
- 10394626
- Date Published:
- Journal Name:
- Journal of Petrology
- Volume:
- 61
- Issue:
- 8
- ISSN:
- 0022-3530
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract Metamorphic rocks from the Connecticut Valley Trough (CVT), Vermont, and Massachusetts, have been examined using quartz‐in‐garnet (QuiG) and conventional thermobarometry, thermodynamic reaction modelling, diffusion modelling, and40Ar/39Ar thermochronology to constrain their
P–T–t paths during Acadian metamorphism and subsequent exhumation. Numerous samples, collected in the vicinity of the Acadian domes, contain garnet porphyroblasts that display cloudy zones characterized by numerous fluid inclusions and modified garnet compositions associated with the replacement of the original garnet by biotite±muscovite±plagioclase±quartz±lowX grs/enrichedX sps. QuiG and conventional thermobarometry constrain both the conditions of garnet nucleation and peakP–T conditions to have occurred at ~0.85–1.05 GPa, ~550–600°C. Most notably, QuiG barometry was performed on inclusions adjacent to these reaction zones in conjunction with Gibbs method reaction modelling to reveal that these dissolution–reprecipitation reactions occurred during nearly isothermal decompression from the peakP–T conditions to around ~0.3 GPa, 550°C. Diffusion modelling reveals that the Mn zoning profiles created during garnet resorption that accompanied decompression formed in less thanc . 3 Ma, which constrains the tectonic exhumation to have occurred at 8–10 mm/year. Subsequent cooling to 500°C occurred rapidly at a rate of 100°C/Ma, followed by slower cooling reaching 1.7°C /Ma by the mid Carboniferous. This is the first reported example of QuiG barometry revealing a multi‐stage metamorphic history and highlights the utility of this method for unravelling complex metamorphic terranes. -
Abstract Garnet–kyanite–staurolite assemblages with large, late porphyroblasts of amphibole form garbenschists in Ordovician volcaniclastic rocks lying immediately south of the Pearya terrane on northernmost Ellesmere Island, Canada. The schist, which together with carbonate olistoliths makes up the Petersen Bay Assemblage (PBA), displays a series of parallel isograds that mark an increase in metamorphic grade over a distance of 10 km towards the contact with Pearya; however, a steep, brittle Cenozoic strike-slip fault with an unknown amount displacement disturbs the earlier accretionary relationship. The late amphibole growth, probably due to fluid ingress, is clear evidence of disequilibrium conditions in the garbenschist. In order to recover the P–T history of the schists, we construct isochemical phase equilibrium models for a nearby garnet–mica schist that escaped the fluid event and compare the results to quartz inclusion in garnet (QuiG) barometry for a garbenschist and the metapelitic garnet schist. Quartz inclusions are confined to garnet cores and the QuiG results, combined with Ti-in-biotite and garnet–biotite thermometry, delineate a prograde path from 480 to 600°C and 0.7 to 0.9 GPa. This path agrees with growth zoning in garnet deduced from X-ray maps of the spessartine component in garnet. The peak conditions obtained from pseudosection modelling using effective bulk composition and the intersection of garnet rim with matrix biotite and white mica isopleths in the metapelite are 665°C at ≤0.85 GPa. Three generations of monazite (I, II and III) were identified by textural characterization, geochemical composition (REE and Y concentrations) and U–Pb ages measured by ion microprobe. Monazite I occurs in the matrix and as inclusions in garnet rims and grew at peak P–T conditions at 397 ± 2 Ma (2σ) from the breakdown of allanite. Monazite II forms overgrowths on matrix Monazite I grains that are oriented parallel to the main schistosity and yield ages of 385 ± 2 Ma. Monazite III, found only in the garbenschist, is 374 ± 6 Ma, which is interpreted as the time of amphibole growth during fluid infiltration at lower temperature and pressure on a clockwise P–T path that remained in the kyanite stability field. These results point to a relatively short (≈12 Myr) Barrovian metamorphic event that affected the schists of the PBA. An obvious heat source is lacking in the adjacent Pearya terrane, but we speculate it was large Devonian plutons—similar to the 390 ± 10 Ma Cape Woods granite located 40 km across strike from the fault—that have been excised by strike-slip. Arc fragments that are correlative to the PBA are low grade; they never saw the heat and were not directly involved in Pearya accretion.more » « less
-
more » « less
Abstract Mineral inclusions are ubiquitous in metamorphic rocks and elastic models for host‐inclusion pairs have become frequently used tools for investigating pressure–temperature (
P–T ) conditions of mineral entrapment. Inclusions can retain remnant pressures () that are relatable to their entrapment P–T conditions using an isotropic elastic model andP–T–V equations of state for host and inclusion minerals. Elastic models are used to constrainP–T curves, known as isomekes, which represent the possible inclusion entrapment conditions. However, isomekes require a temperature estimate for use as a thermobarometer. Previous studies obtained temperature estimates from thermometric methods external of the host‐inclusion system. In this study, we present the firstP–T estimates of quartz inclusion entrapment by integrating the quartz‐in‐garnet elastic model with titanium concentration measurements of inclusions and a Ti‐in‐quartz solubility model (QuiG‐TiQ). QuiG‐TiQ was used to determine entrapmentP–T conditions of quartz inclusions in garnet from a quartzofeldspathic gneiss from Goodenough Island, part of the (ultra)high‐pressure terrane of Papua New Guinea. Raman spectroscopic measurements of the 128, 206, and 464 cm−1bands of quartz were used to calculate inclusion pressures using hydrostatic pressure calibrations (), a volume strain calculation ( ), and elastic tensor calculation ( ), that account for deviatoric stress. values calculated from the 128, 206, and 464 cm−1bands’ hydrostatic calibrations are significantly different from one another with values of 1.8 ± 0.1, 2.0 ± 0.1, and 2.5 ± 0.1 kbar, respectively. We quantified elastic anisotropy using the 128, 206 and 464 cm−1Raman band frequencies of quartz inclusions and stRAinMAN software (Angel, Murri, Mihailova, & Alvaro, 2019, 234 :129–140). The amount of elastic anisotropy in quartz inclusions varied by ~230%. A subset of inclusions with nearly isotropic strains gives an averageand of 2.5 ± 0.2 and 2.6 ± 0.2 kbar, respectively. Depending on the sign and magnitude, inclusions with large anisotropic strains respectively overestimate or underestimate inclusion pressures and are significantly different (<3.8 kbar) from the inclusions that have nearly isotropic strains. Titanium concentrations were measured in quartz inclusions exposed at the surface of the garnet. The average Ti‐in‐quartz isopleth (19 ± 1 ppm [2 σ ]) intersects the average QuiG isomeke at 10.2 ± 0.3 kbar and 601 ± 6°C, which are interpreted as theP–T conditions of quartzofeldspathic gneiss garnet growth and entrapment of quartz inclusions. TheP–T intersection point of QuiG and Ti‐in‐quartz univariant curves represents mechanical and chemical equilibrium during crystallization of garnet, quartz, and rutile. These three minerals are common in many bulk rock compositions that crystallize over a wide range ofP–T conditions thus permitting application of QuiG‐TiQ to many metamorphic rocks. -
Abstract. Retrograde metamorphic rocks provide key insights into the pressure–temperature (P–T) evolution of exhumed material, and resultant P–T constraints have direct implications for the mechanical and thermal conditions of subduction interfaces. However, constraining P–T conditions of retrograde metamorphic rocks has historically been challenging and has resulted in debate about the conditions experienced by these rocks. In this work, we combine elastic thermobarometry with oxygen isotope thermometry to quantify the P–T evolution of retrograde metamorphic rocks of the Cycladic Blueschist Unit (CBU), an exhumed subduction complex exposed on Syros, Greece. We employ quartz-in-garnet and quartz-in-epidote barometry to constrain pressures of garnet and epidote growth near peak subduction conditions and during exhumation, respectively. Oxygen isotope thermometry of quartz and calcite within boudin necks was used to estimate temperatures during exhumation and to refine pressure estimates. Three distinct pressure groups are related to different metamorphic events and fabrics: high-pressure garnet growth at ∼1.4–1.7 GPa between 500–550 ∘C, retrograde epidote growth at ∼1.3–1.5 GPa between 400–500 ∘C, and a second stage of retrograde epidote growth at ∼1.0 GPa and 400 ∘C. These results are consistent with different stages of deformation inferred from field and microstructural observations, recording prograde subduction to blueschist–eclogite facies and subsequent retrogression under blueschist–greenschist facies conditions. Our new results indicate that the CBU experienced cooling during decompression after reaching maximum high-pressure–low-temperature conditions. These P–T conditions and structural observations are consistent with exhumation and cooling within the subduction channel in proximity to the refrigerating subducting plate, prior to Miocene core-complex formation. This study also illustrates the potential of using elastic thermobarometry in combination with structural and microstructural constraints, to better understand the P–T-deformation conditions of retrograde mineral growth in high-pressure–low-temperature (HP/LT) metamorphic terranes.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/petrology/egaa076
