Skip to main content

Physicochemical conditions and timing of rodingite formation: evidence from rodingite-hosted fluid inclusions in the JM Asbestos mine, Asbestos, Québec

Abstract

Fluid inclusions and geological relationships indicate that rodingite formation in the Asbestos ophiolite, Québec, occurred in two, or possibly three, separate episodes during thrusting of the ophiolite onto the Laurentian margin, and that it involved three fluids. The first episode of rodingitization, which affected diorite, occurred at temperatures of between 290 and 360°C and pressures of 2.5 to 4.5 kbar, and the second episode, which affected granite and slate, occurred at temperatures of between 325 and 400°C and pressures less than 3 kbar. The fluids responsible for these episodes of alteration were moderately to strongly saline (~1.5 to 6.3 m eq. NaCl), rich in divalent cations and contained appreciable methane. A possible third episode of alteration is suggested by primary fluid inclusions in vesuvianite-rich bodies and secondary inclusions in other types of rodingite, with significantly lower trapping temperatures, salinity and methane content. The association of the aqueous fluids with hydrocarbon-rich fluids containing CH4 and higher order alkanes, but no CO2, suggests strongly that the former originated from the serpentinites. The similarities in the composition of the fluids in all rock types indicate that the ophiolite had already been thrust onto the slates when rodingitization occurred.

Introduction

Calc-silicate rocks associated with serpentinites are called rodingites, a name which was given originally by Bell et al. [1] to lime-rich, coarse- to fine-grained, gabbro-like rocks composed primarily of grossular or prehnite in the vicinity of the Roding River in the Dun Mountain ultramafic complex, New Zealand. Since the first description of rodingites by Bell et al. [1], these rocks have been studied by numerous authors in a variety of tectonic settings containing serpentinized peridotite [213]. For example, calc-silicate rocks are frequently found at the contacts between serpentinites and mafic or felsic intrusive bodies. Calc-silicate rocks also commonly form rinds around tectonic inclusions of these intrusions in serpentinites. Finally, calc-silicate rocks occur at the contacts between serpentinized peridotites and country rocks of diverse origin. In all cases, these calc-silicate rocks are metasomatic products of the alteration of the various rock types in contact with serpentinite, and, in some cases, of the serpentinite itself. As a result, the term rodingite is no longer restricted to altered gabbros but now refers to any rock type that contains calc-silicate-rich mineral assemblages and occurs in contact with serpentinite [2, 6, 8].

One of the reasons why rodingites have attracted considerable attention, is that their complex mineralogy provides important opportunities to use mineral stability relationships to determine the conditions under which serpentinization, or de-serpentinization, occurs [10, 11, 14]. However, rodingites also commonly host abundant primary fluid inclusions, which, in principle, can also be used for this purpose. Despite this, very few researchers have studied fluid inclusions in rodingites and fewer still have conducted the detailed microthermometric studies needed to establish the pressure-temperature conditions of serpentinization. Moreover, the temperature and pressures that have been reported vary widely not only among deposits but also within a single deposit. Thus, for example, O'Hanley et al. [11] reported a temperature of ~300°C and a pressure of < 800 bar for the Cassiar serpentinite, British Columbia, Dubinska et al. [13] temperatures of 270 to 300°C and a pressure of ~1 kbar for a metasomatic shell around rodingite from the Jordanów-Gogołów serpentinite massif, Poland, Mittwede and Schandl [15] temperatures of 350 to 500°C and pressures of 1 to 2.5 kbar for the Hammett Grove serpentinites, and Schandl and Mittwede [16] temperatures of 250 to 450°C and a pressure of ~3 kbar for the Acipayam serpentinites, Turkey.

In this paper, we present results from a detailed study of fluid inclusions in rodingitized slate, diorite and granite, and vesuvianite-rich veins from the JM Asbestos mine, Asbestos, Québec. The data presented here, in combination with published information on the tectonic evolution of the ophiolite complexes of the south-eastern Quebec Appalachians, help to constrain the pressure and temperature conditions during which serpentinization took place. Moreover, they establish the chemistry of the fluids responsible for rodingitization, and help clarify the timing relationships of serpentinization and rodingitization during obduction of the Asbestos ophiolite onto the Laurentian continental margin.

Geological setting

The JM Asbestos mine is located in serpentinized ultramafic rocks that form part of the Ordovician Asbestos ophiolite, which was obducted onto the margin of the Laurentian craton during the ~470–460 Ma Taconic orogeny [17, 18]. The ultramafic rocks at the base of the ophiolite comprise serpentinized harzburgite tectonites overlain by serpentinized cumulate dunites and pyroxenites [19] (Figure 1). These rocks are in fault contact with greenschist facies slates and greywackes that were originally deposited on the Laurentian continental margin [20, 21]. Dioritic-monzodioritic and later granodioritic-granitic rocks, many of which contain primary magmatic andalusite [22], were intruded into the basal ultramafic portion of the ophiolite prior to emplacement into its present tectonic setting. According to Laurent and Hébert [23], the dioritic intrusive rocks are consanguineous with the ophiolite and were emplaced into peridotites prior to serpentinization. By contrast, these authors [23] suggested, based on field evidence, that the later granitic rocks were emplaced in at least partially serpentinized peridotite. Concurrent with serpentinization and transportation of the ophiolite, the intrusive rocks were subjected to sodic-calcic and rodingitic hydrothermal alteration [22]. These two alteration facies also developed in the continentally derived metasedimentary units (slate) below the ophiolite during thrusting of the latter onto the Laurentian continental margin [22, 24]. The fluids responsible for rodingitization were trapped in a variety of rodingite minerals and are the principal focus of this paper.

Figure 1
figure 1

Geology of the Asbestos ophiolite in the area of the JM Asbestos mine, Asbestos, Québec (After maps by Johns-Manville Inc, Bédard et al. [18], Hébert [19]). The contact between the ophiolite and the Caldwell Group metasediments shows evidence of earlier thrusting and later normal faulting.

Rodingites

A large variety of rodingitized rocks, containing a complex array of mineral assemblages, were studied in the field, in hand samples and in thin sections. These comprise altered serpentinite, felsic intrusive dykes that cut the serpentinite, and slate that lies in thrust contact with serpentinites. A suite of samples representative of the various types of rodingites was described in detail in Normand [22] and comprises four samples of rodingitized diorite, five samples of rodingitized granite, three rodingites that replaced slate, and two samples of vein-like rodingite rich in vesuvianite (Table 1). Mineral identification and textural relationships were established using an optical binocular microscope and an electron microprobe that has scanning electron microscopic capability. Minerals characteristic of the rodingites comprise grossular and hydrogrossular, clinopyroxenes that belong to the diopside-hedenbergite series (simply referred to as clinopyroxene below), epidote-group minerals (largely zoisite), chlorite-group minerals, prehnite, vesuvianite and wollastonite. The essential characteristics of the different rodingite types are described below.

Table 1 Description of samples.

Rodingitized diorite

The thinner dykes of rodingitized diorite (< 15 cm) are generally composed of grossular-clinopyroxene-predominant assemblages, whereas the mineral assemblages in the thicker ones (> 2 m) also include zoisite and prehnite. The mineral assemblages observed in the pervasively altered diorites include grossular + clinopyroxene + chlorite, grossular + clinopyroxene, grossular + clinopyroxene + zoisite or grossular + clinopyroxene + zoisite + prehnite and zoisite + clinopyroxene. Veins cutting the pervasively rodingitized diorites contain highly variable proportions of the same phases present in the host rock (e.g., Dior-2b*; Table 1). Where they cut diorite, which has not been pervasively rodingitized, the veins measure up to 8 cm thick and are prehnite-, feldspar-(either albite or K-feldspar) or more rarely grossular-rich. A rim of prismatic to acicular clinopyroxene crystals is commonly observed on the walls of these veins, whereas the core is filled with coarse-grained Ca-Al-silicate minerals. In addition to these veins, there are also veinlets composed of dark brown vesuvianite and chlorite cutting grossular + clinopyroxene rodingites (Table 1; Dior-3).

Sample Dior-1 (Table 1) represents an example of the smaller, completely rodingitized diorite dykes. The dyke is 8 cm thick and shows blackwall margins that are in sharp contact with the host serpentinized harzburgites. This contact is not deformed. The rodingite is composed of fluid inclusion-bearing grossular and clinopyroxene, and minor interstitial proportions of chlorite-group minerals.

Sample Dior-2 (Table 1) was collected from a pervasively rodingitized margin of a 3.84 m thick diorite dyke. The serpentinites in contact with the dyke were sheared and replaced by antigorite. A thin, less than 2 cm thick, biotite-bearing blackwall zone of alteration occurs in the dyke at the contact with serpentinite. The diorite was completely rodingitized to a pinkish zoisite + diopside + garnet ± prehnite assemblage over a distance of 7 cm adjacent to the blackwall, grading sharply into strongly rodingitized diorite with relict biotite up to 18 cm from the contact with serpentinite. Fluid inclusions are common in zoisite and grossular.

Rodingitized granite

Two contrasting types of mineral assemblage are present in the two dykes of rodingitized granitic rock investigated in this study. The first dyke is a pegmatitic leucogranite in which the assemblage chlorite >> grossular + prehnite locally forms a thick zone of blackwall-type alteration near its contact with serpentinite (Table 1; Gran-1a), whereas the pervasively rodingitized granite is composed of a coarse-grained assemblage of grossular + clinopyroxene (Table 1; Gran-1b) which contains abundant fluid inclusions.

The second dyke (large block; provenance in the mine unknown) is a biotite granite which was mylonitized at the margins prior to alteration. The first four cm of the dyke from the contact with serpentinite were altered to chlorite + hydrogrossular + clinopyroxene (Table 1; Gran-2a), and are followed by a six cm wide zone composed of wollastonite + vesuvianite + clinopyroxene + hydrogrossular (Table 1; Gran-2b, -2c), and a 35 cm wide zone of bleached, albite and K-feldspar-dominant rock. Fluid inclusions suitable for microthermometric measurements were not identified in the samples.

Rodingitized slate

Rodingitization of slate was preceded by an episode of sodic-calcic alteration, which produced a halo of albite-tremolite-biotite rock extending up to two meters from the contact with serpentinite. Relict patches of this type of alteration are commonly found in the rodingitized slate, which forms a zone ~60 cm thick between albitized slate and a ~10 cm thick zone of blackwall-altered slate immediately adjacent to the serpentinites. The original sedimentary textures were preserved during albitization and rodingitization.

At one location, pervasive rodingitization of the slate produced a fine-grained calc-silicate mineral assemblage composed of 1) colorless hydrogrossular + zoisite + clinopyroxene ± prehnite ± K-feldspar at the contact with the blackwall (Table 1; Slate-1a) and 2) colorless hydrogrossular + prehnite + clinopyroxene ± K-feldspar further from this contact (Table 1; Slate-1b). The rodingitized slate adjacent to the blackwall is cut by numerous veins that contain fluid inclusion-rich orange grossular and clinopyroxene crystals measuring up to 5 mm in diameter, and lesser proportions of space-filling prehnite practically devoid of fluid inclusions. No fluid inclusions were observed in the fine-grained mineral assemblages that compose the pervasively altered slate.

At one other location, a 20 cm zone of slate in contact with blackwall was pervasively replaced by a porous assemblage of coarse-grained, fluid inclusion-bearing clinopyroxene and small proportions of prehnite (Table 1; Slate-2).

Vesuvianite-rich veins

Vesuvianite-rich veins are concentrated at the base of the cumulate dunite-pyroxenite unit near the transition into the tectonized harzburgite. They commonly display a very high vug porosity (visually estimated at up to 25%). Many of the veins have thick haloes of diopsidized serpentinite, which are cut by vesuvianite-rich veinlets. Mineral parageneses include white grossular + green and purple vesuvianite + clinopyroxene, vesuvianite + clinopyroxene, vesuvianite + wollastonite and vesuvianite + clinopyroxene + chlorite (see Table 1).

Fluid inclusions

The most representative examples of fluid inclusion-bearing rodingitized dykes and slate were subjected to a detailed fluid inclusion study. Twenty-four doubly polished thin sections (150–200 μm thick) were prepared from two samples of rodingitized diorite (Dior-1, -2), one sample of rodingitized granite (Gran-1b), two samples of rodingitized slate (Slate-1a, -2) and two samples of vesuvianite from vesuvianite-rich "veins" (Ves-1, -2). Selected inclusions were analysed microthermometrically using a U.S.G.S.-type heating-freezing stage mounted on a Leitz Sm-Lux Pol polarizing microscope. Calibration was performed between -56.6 and 374.1°C using synthetic fluid inclusion standards of water and carbon dioxide.

Fluid inclusion petrography

Two types of primary fluid inclusions were recognized from microscopic examination of the samples at 25°C. Type 1 inclusions consist of a single fluid phase and type 2 inclusions of two immiscible fluids (an aqueous fluid and a carbonic phase forming a bubble). A primary origin for these inclusions was inferred based on their distributions as three dimensional groups or occurrence as comparatively large, isolated inclusions. Where fluid inclusions formed planar arrays within or cutting across crystals, they were interpreted to be of secondary origin.

Type 1 inclusions measure up to 20 μm in diameter, commonly appear very dark under the microscope and have variable shapes (irregular, tubular and rectangular). They are abundant in sample Dior-2, where they occur as three dimensional clusters in zoisite and grossular. However, they were not found in sample Dior-1. Type 1 fluid inclusions also occur as three dimensional clusters in clinopyroxene in the rodingitized slate. They are rare in sample Slate-1a (an isolated group composed of four inclusions was observed) and common in sample Slate-2. Only two type 1 fluid inclusions were detected in vesuvianite from the vesuvianite-rich veins (Ves-1), where they formed an isolated group. The occurrence of type 1 inclusions in zoisite from sample Dior-2 is illustrated in Figure 2.

Figure 2
figure 2

Type 1 fluid inclusions in zoisite from sample Dior-2. (A): prismatic zoisite crystals showing fractures running perpendicularly to the elongation. (B): enlargement of inset in A. The black arrow points to an irregularly shaped and dark colored type 1 fluid inclusion which forms part of a three dimensional cluster. Type 2a fluid inclusions coexist with the type 1 inclusions in the same clusters (the white arrow points to a slightly out of focus example). (C): detail of the irregular type 1 fluid inclusion shown in B.

The two-phase type 2 inclusions occur in all the rodingite-forming minerals. They are generally prismatic in clinopyroxene (Figure 3B and 3E) and vesuvianite (Figure 3F), and oriented parallel to the c axis (and cleavage in clinopyroxene). In garnet and zoisite, they are prismatic, equant or irregularly shaped (Figure 3A and 3D). Their diameter varies between < 1 μm and ~50 μm (one tubular inclusion in garnet in sample Gran-1b measured 150 μm in length). The volume fraction occupied by the aqueous solution at room temperature in type 2 inclusions was visually estimated to be 84 ± 6% in rodingitized diorite, 86 ± 4% in rodingitized granite and slate, and approximately 95% in vesuvianite. Type 2 inclusions were considered primary when they were large and isolated, or occurred in tight three dimensional groups away from the contacts between grains. They were considered secondary when they occurred in planar arrays.

Figure 3
figure 3

A serie of photomicrographs showing the distribution, shape and interrelationships with host minerals and other inclusions of primary type 2 fluid inclusions. (A): Comparatively large type 1 and type 2b fluid inclusion coexisting in the same zoisite crystal (sample Dior-2). The type 2b inclusion is isolated. Two smaller type 1 inclusions are present near the lower left part of the one indicated by a black arrow, which together form a small group. (B): Rod-shaped fluid inclusions in a diopside crystal (sample Dior-1). The inclusions are randomly distributed and are oriented parallel to the c axis of the crystal. (C): Image showing acicular inclusions in grossular (Gr) from sample Gran-1b. Some of the inclusions consist of diopside needles and others are fluid filled. The occurrence of these inclusions is restricted to certain sectors within larger grossular crystals (black arrow). Type 2a fluid inclusions in grossular most commonly occur in these acicular clusters. Coarse-grained diopside crystals (Di) are abundant in the grossular and contain the same type of fluid inclusions. (D): Detail of an elongated type 2a fluid inclusion in grossular from sample Gran-1b. (E): Image showing diopside filling interstices between euhedral grossular crystals in veins from sample Slate-1a. A type 2a fluid inclusion in the diopside is shown by a black arrow. The grossular also contain type 2a fluid inclusions. Type 2a fluid inclusions in samples Slate-1a and Gran-1b showed the same behaviour during cryogenic experiments. (F): A type 1c fluid inclusion (black arrow) in vesuvianite sample Ves-2. This inclusion is oriented parallel to the c axis of the host crystal. Other needle-shaped inclusions can be seen oriented in the same way. They consist mostly of empty cavities.

Type 1 and type 2 inclusions coexist in the same crystals in grossular and zoisite from one of the rodingitized diorite samples (sample Dior-2; Figure 2B and 3A). In rodingitized slate, type 1 and type 2 fluid inclusions were observed in the same crystals only in sample Slate-2. The two type 1 fluid inclusions observed in vesuvianite (sample Ves-1) were in close spatial association with type 2 inclusions.

Microthermometry

Cryogenic experiments

The fluid inclusions displayed a variety of phase changes, which are summarized in Table 2. All type 1 inclusions in rodingite nucleated a bubble when they were cooled to temperatures below -90°C. One inclusion in sample Dior-2 contained solid, liquid and gas at a temperature of -184.4°C. On being heated, these inclusions homogenized by dissolution of the vapour bubble into the denser liquid phase. This occurred at temperatures varying between -123 and -89°C (Figure 4), suggesting that the type 1 inclusions consist dominantly of a CH4-rich fluid. The lowest homogenization temperatures were obtained for inclusions in sample Dior-2 (-119°C on average). Homogenization temperatures for type 1 inclusions in rodingitized slate (Slate-1a, -2) were between -120 and -89°C. The two type 1 inclusions in sample Ves-1 homogenized at -90°C.

Table 2 Microthermometric data.
Figure 4
figure 4

Histogram showing the distribution of total homogenization temperatures of type 1 fluid inclusions in rodingitized diorite and slate, and in vesuvianite from vesuvianite-rich veins.

Type 2 inclusions can be divided into three sub-types based on the melting behaviour of ice and clathrate, and on whether or not a vapour-liquid phase change was observed. Type 2a inclusions have low ice and clathrate melting temperatures. They occur in the rodingitized granite sample Gran-1b and in veins in the rodingitized slate sample Slate-1a. Type 2b inclusions have low to moderately low ice melting temperatures and moderate to elevated clathrate melting temperatures. They occur in the samples of rodingitized diorite and slate (Dior-1, -2; Slate-1a, -2). Type 2c inclusions form a sub-population of type 2b inclusions that are characterized by low density (see below), and occur in vesuvianite (Ves-1, -2).

After being cooled to a temperature of -196°C, type 2a inclusions contained ice, a greenish birefringent hydrate, a clathrate, and solid and vapour (and liquid?) hydrocarbons. The first phase change during heating was the melting of a solid or disappearance of liquid (difficult to distinguish) inside the inner fluid at temperatures between -196 (observed in two inclusions in sample Gran-1b) and -183°C (observed in one inclusion in sample Gran-1b). This suggests that the bubbles are composed of CH4 and another component representing a system with a triple point located at lower temperature than that of CH4 (triple point -182.5). Between -102 and -85°C, a vapour bubble homogenized inside the inner fluid (liquid), further suggesting that the latter is composed predominantly of CH4. Melting of the ice started between -74 and -42°C in sample Gran-1b, and between -76 and -48°C in sample Slate-1a (Figure 5). The ice finally melted between -27.3 and -23.4°C in sample Gran-1b and between -25.6 and -22.8°C in sample Slate-1a. Upon further heating, the clathrate melted at temperatures between -7.1 and -3.2°C in sample Gran-1b, and between -6.4 and -1.8°C in sample Slate-1a. The hydrate phase persisted to temperatures between -6.5 and 6.7°C in sample Slate-1a, and between -4.6 and 20.8°C in sample Gran-1b. A pinkish solid was observed in two inclusions in sample Slate-1a and melted at a temperature of between 0.2 and 0.3°C above the melting point of the hydrate phase (-5.6 and -5.4°C). The salinities calculated from last ice melting temperatures, in terms of equivalent NaCl and uncorrected for the presence of clathrate, vary between 24 and 27 wt. % (5.3–6.3 m).

Figure 5
figure 5

Histograms showing the distribution of the initial ice melting temperatures of type 2a and 2b fluid inclusions in rodingitized granite (Gran-1b), rodingitized slate (Slate-1a, -2), and rodingitized diorite (Dior-1, -2), and of type 2c fluid inclusions in vesuvianite from vesuvianite-rich veins (Ves-1, -2).

The concentration of initial ice melting temperatures around -52°C (Figure 5) and the occurrence of five inclusions with initial ice melting at temperatures of -67 to -76°C suggest strongly that the salts dissolved in type 2a fluid inclusions are predominantly NaCl and CaCl2. The eutectic temperature of the system NaCl-CaCl2-H2O is -52°C and there is a metastable eutectic or recrystallization event that takes place at around -70°C (see [2527]). The sequence of eutectic, final ice and hydrate melting temperatures of type 2a fluid inclusions in sample Slate-1a can be used to graphically estimate the NaCl/(NaCl+CaCl2) molar ratio of these inclusions using the data presented in Williams-Jones and Samson [28]. The bulk fluid composition of two type 2a inclusions, that displayed sequences of phase change consistent with the data of Williams-Jones and Samson [28], suggests NaCl/(NaCl+CaCl2) molecular ratios of 0.84 and 0.85, respectively. Calculation of the NaCl/(NaCl+CaCl2) molecular ratio of nine inclusions using CalcicBrine [29] suggests values between 0.76 and 0.97, with an average of 0.84.

Behaviour similar to that observed in type 2a inclusions was noted during heating of type 2b inclusions from -196 to 25°C. The main difference was in the temperature of melting of the clathrate, which varied from 0 to 13°C. Type 2b inclusions show a much wider range of clathrate and ice melting temperatures. However, these temperatures are positively correlated in type 2b inclusions, suggesting that the latter were trapped under similar pressure and temperature conditions. The relationship between the clathrate melting temperatures and the salinity in wt. % equivalent NaCl (based on final ice melting temperatures, not corrected for the presence of clathrate and other salts) of type 2 inclusions is illustrated in Figure 6. Curves for internal inclusion pressures of 25, 50, 100, 150 and 200 bar at the temperature of final clathrate melting for the system CH4-NaCl-H2O are also shown for comparison to illustrate the effects of temperature and pressure on the equilibrium fluid compositions. The curves were calculated using the program CURVES [30]. The data indicate that type 2b inclusions are denser than type 2a inclusions, and, consistent with their coexistence with high-density type 1 inclusions, that they were thus most likely trapped at higher pressure and/or lower temperature than type 2a inclusions. The lower final clathrate melting temperatures and the higher total homogenization temperatures of type 2a fluid inclusions in sample Gran-1b relative to the high salinity type 2b inclusions in sample Dior-2 is consistent with trapping conditions of a lower density fluid in sample Gran-1b at higher temperature. First ice melting temperatures of 42 type 2b inclusions ranged from -29 to -57°C, with most data being concentrated between -40 and -50°C suggesting the presence of a significant proportion of dissolved divalent cations.

Figure 6
figure 6

Diagrams showing the relationship between the final melting temperature of the clathrate and the salinity (uncorrected for the presence of clathrate) of type 2 fluid inclusions in rodingitized diorite (A; samples Dior-1, -2), rodingitized granite (B; sample Gran-1b), rodingitized slate (C; samples Slate-1b, -2), and in vesuvianite (D; sample Ves-1). Primary type 2 fluid inclusions form three distinct populations: 1) Low density type 2a fluid inclusions characterized by high salinity and low final clathrate melting temperatures, 2) high density type 2b fluid inclusions showing a wide range of salinities and elevated final clathrate melting temperatures, and 3) low density and low salinity type 2c fluid inclusions found in vesuvianite. Also shown on the diagrams are isobars at pressures of 25, 50, 100, 150 and 200 bar for the system CH4-NaCl-H2O calculated using program CURVES [30]. Filled symbols represent data for primary fluid inclusions and open symbols represent data for secondary fluid inclusions. Refer to the text for a detailed discussion.

Whereas in type 2a and 2b inclusions a vapour bubble formed in the inner fluid at low temperatures (< -82°C), this was not the case in type 2c inclusions. This suggests that the bubble in these latter inclusions contains a much lower proportion of CH4, i.e., that it is composed predominantly of water vapour, or that it consists of a low density CH4 phase. The ice began to melt between -42 and -35°C, and finally melted between -5.7 and -4.5°C. Clathrate was the last solid to melt in type 2c inclusions, and did so between 4.8 and 10.8°C. The clathrate was very difficult to observe, further suggesting that only a very small amount of CH4 was dissolved in the vapour bubble.

Heating experiments

Almost all inclusions ≥ ~10 μm in diameter decrepitated upon being heated above 100°C, whether or not they had been frozen previously, consistent with the presence of CH4-rich inclusions which have high internal pressures at room temperature and steep isochoric paths. To avoid problems related to leakage or stretching of inclusions, heating experiments were performed after freezing experiments. High temperature phase changes could only be successfully measured for very small inclusions, which in many cases were too small to see all the low temperature phase changes.

Complete homogenization of type 2a inclusions in sample Gran-1b occurred at temperatures between 288 and 385°C, and there was a tight clustering of homogenization temperatures at 344 ± 15°C (Figure 7). In sample Dior-1, the homogenization temperatures of type 2b inclusions varied between 277 and 326°C, and in sample Dior-2, they were between 286 and 333°C. Secondary inclusions in sample Dior-2 (cryogenic data similar to those of low salinity type 2b and 2c inclusions) homogenized at temperatures between 222 and 275°C. Homogenization temperatures of primary inclusions in rodingitized slate (Slate-1a) were concentrated between 290 and 344°C. Five secondary inclusions homogenized at temperatures between 230 and 269°C. Primary type 2c inclusions of prismatic habit in vesuvianite (Ves-2) homogenized at a mean temperature of 191°C.

Figure 7
figure 7

Histograms showing the distribution of the temperatures of total homogenization of type 2a and 2b fluid inclusions in rodingitized granite (Gran-1b), rodingitized slate (Slate-1a), and rodingitized diorite (Dior-1, -2), and of type 2c fluid inclusions in vesuvianite from vesuvianite-rich veins (Ves-2).

Gas chromatographic analyses

Gas chromatographic analyses of the inclusion fluids in samples Slate-1a, Gran-1b, Dior-2 and Ves-1 confirmed the presence of CH4. The latter formed the principal carbonic species (> 90 mol % excluding water) in the volatile mixtures. Hydrocarbons detected in addition to CH4 are C2H6, C3H8, n-C4H10 and n-C5H12. Carbon dioxide was not detected in any of the samples.

Small amounts of N2 or CO or Ar or H2 were also detected. However, co-elution of N2, CO, Ar, and H2 did not permit differentiation between these species by gas chromatographic analysis using He as the carrier gas. Results of gas chromatographic analyses for selected rodingite samples are presented in Table 3. Also included in this table are values for N2 assuming this species represents the unknown gas. The values for H2O reported in Table 3 should be considered qualitative as a good calibration curve for this species could not be obtained. The lower H2O/CH4 ratio in sample Dior-2 is consistent with the abundant occurrence of type 1 inclusions in the sample.

Table 3 Volatile content of selected samples measured by gas chromatography (in nanomoles).

Assuming that the unidentified gas is N2, X N 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqqGybawdaWgaaWcbaGaeeOta40aaSbaaWqaaiabikdaYaqabaaaleqaaaaa@305C@ in the bulk fluid (excluding H2O) would be ~0.05 and 0.07 for samples Dior-2 and Gran-1b, respectively, ~0.09 for sample Slate-1a, and up to ~0.35 for vesuvianite. The similarity of these values among the samples of rodingitized felsic intrusive rocks and slates suggests that the difference in the inferred internal pressure of type 2a and type 2b inclusions is not due to large variations in the proportion of the volatiles, but is real. However, the proportions of this unidentified species in the samples from vesuvianite-rich veins, which contain type 1 and type 2c inclusions, are significantly higher.

Details on the analytical procedures and results, and a discussion of the origin of the hydrocarbons will be presented elsewhere.

Table 4 Calculated type 2a and 2b fluid inclusion compositions and homogenization pressures.

V-X properties

Microthermometric and bulk gas chromatographic data indicate that type 1 inclusions in rodingitized diorite and slate are composed predominantly of CH4. Although a rim of aqueous solution could not be detected optically, this does not mean that it is not present in type 1 inclusions. Very low proportions of an aqueous solution in the form of a thin film on the walls of an inclusion may appear invisible. The above notwithstanding, the molar volume of type 1 inclusions was calculated from their homogenization temperature assuming, as a first order approximation, that they are composed of only one phase and that this phase is pure CH4. The calculations were performed using the program BULK [30, 31] and the equation of state for pure CH4 of Duan et al. [32, 33].

The mole fraction of components in type 2 fluid inclusions was estimated using the programs CURVES (version 12/02) and ICE (version 12/02), implemented in the package of programs CLATHRATES [30, 34]. Fugacity calculations based on clathrate equilibria and molar volumes of the aqueous solutions at the temperature of final clathrate melting were estimated using equations of state in Duan et al. [32, 33]. The program DENSITY could not be used as clathrate always formed in the fluid inclusions before a vapour bubble nucleated in the CH4 phase [30, 34]. The program ICE could be used only for fluid inclusions in vesuvianite containing a low density gas phase [30].

Owing to the large errors inherent in the visual estimation of the bubble volume fraction in fluid inclusions (especially if they are irregularly shaped), we used a method to estimate the volumetric properties that involved application of the programs CURVES and GEOFLUIDS model 1 (available at http://geotherm.ucsd.edu; [3237]). The method consisted of first determining the mole fractions of the inclusion components and their molar volumes for a range of liquid proportions relative to fluid inclusion volume using CURVES. Compositional data determined using CURVES were then used in GEOFLUIDS to obtain molar volumes for the inclusions at CH4 saturation (i.e., at the temperature of total homogenization). The correct composition of a specific inclusion corresponded to that for which the molar volumes calculated using the two programs corresponded. This method also allowed determination of the minimum trapping pressure. Lamb et al. [38] have shown that location of the solvus in the system CH4-H2O-NaCl calculated using GEOFLUIDS is consistent with their experimentally determined data at temperatures between 300 and 600°C, and pressures of 1 and 2 kbar [38, 39].

Program CURVES permits determination of inclusion composition provided the salt content is known. An assumption about the salt concentration in the aqueous phase was therefore necessary. For this, we used last ice melting data to obtain an equivalent wt.% NaCl. We realize that this method may overestimate the salinity as the presence of clathrate at final ice melting indicates that the aqueous phase is impoverished in water. In addition, first ice melting data indicate that the salts dissolved in the inclusions comprise a significant proportion of divalent cations. However, in the absence of additional compositional data, this was the best approach available to us. The only salt considered by GEOFLUIDS Model 1 is NaCl. Thus, the homogenization pressures determined using this program may be slightly overestimated.

We selected data from two representative type 2b fluid inclusions in sample Dior-2 for the calculation of fluid inclusion composition and total homogenisation pressure in rodingitized diorite. The first inclusion, with a final ice melting temperature of -14.2°C and a clathrate melting temperature of 9.1°C, is typical of a large proportion of the fluid inclusions. The volume fraction of aqueous fluid in this inclusion was visually estimated to be 86%. The second inclusion, with a final ice melting temperature of -22.4°C and a clathrate melting temperature of 3.7°C, is typical of the most saline inclusions in the sample. The volume fraction of aqueous fluid in this inclusion was also visually estimated to be 86%. Using the equation of Bodnar [40], the wt.% equivalent concentration of NaCl in the two fluid inclusions was calculated to be 18 and 24, respectively. The latter value slightly exceeds the range of salinities that can be reliably estimated with the equation (the eutectic temperature in the system NaCl-H2O is -21.2°C) but is considered to be a close approximation. The final homogenization temperature could not be obtained for these fluid inclusions. However, as indicated above, type 2b inclusions in the sample homogenized between 286 and 333°C and at an average temperature of 312°C. Applying the method described above, the volume fraction of aqueous fluid at final clathrate melting for the first inclusion was calculated to be between 84 and 80%, corresponding to total homogenization temperatures of 301 (the lowest temperature for which the program GEOFLUID is designed) and 333°C, respectively. For the second inclusion, the volume fraction of aqueous fluid at final clathrate melting was calculated to vary between 85 and 81%. The composition and molar volumes of the fluid of the two inclusions calculated using the method are reported in Table 4. The molar fraction of methane is estimated to vary between 0.033 and 0.051.

Volumetric and compositional values for type 2a fluid inclusions in the rodingitized granite sample Gran-1b and the rodingitized slate sample Slate-1a were estimated in the same manner. A representative type 2a fluid inclusion was selected from data collected for sample Gran-1b with a final ice melting temperature of -23.8°C and a clathrate melting temperature of -6.1°C. As stated above, use of the program GEOFLUID is limited to the system H2O-CH4-NaCl. A final ice melting temperature of -23.8°C in the system H2O-NaCl indicates salinities too elevated to be used in calculations using the program CURVES (where calculations are limited to fluids having salinities at or below that corresponding to the eutectic temperature of the H2O-NaCl system), Volumetric and compositional data for the type 2a inclusion considered where thus extrapolated from a set of data calculated at salinities corresponding to the eutectic temperature and below. Again, a complete data set with final ice and clathrate melting, and total homogenization temperatures could not be compiled. The final homogenization temperatures of type 2a inclusions in sample Gran-1b, however, form a tight grouping at 344 ± 15°C. Calculations of composition and volume were therefore conducted for temperatures between 330 and 360°C. The results are reported in Table 4. A data set including the temperature of final ice and clathrate melting, and the final homogenization temperature was obtained for two type 2a includions in sample Slate-1a. The final homogenization temperatures for these two inclusions were 300 and 322.5°C. Because the data set of homogenization temperatures has a mode at 319°C in this sample, calculations were performed only for the inclusion that homogenized at 322.5°C (Table 4). As is immediately apparent from inspection of Table 4, the calculated homogenization pressures are consistent with the lower density of fluid inclusions in samples Gran-1b and Slate-1a inferred above from clathrate melting temperatures and salinity relationships.

The behaviour of type 2c fluid inclusions in vesuvianite during cryogenic experiments indicates that these inclusions are of relatively low density and that the salt dissolved in the aqueous phase is predominantly NaCl. In the absence of detailed information on the nature of the salts dissolved in the fluid inclusions, we assumed, as a first order approximation, that the aqueous phase contained only NaCl. Based on estimates of the bubble volume fraction and ice and clathrate melting temperatures, the program ICE yields molar fractions of CH4, H2O and NaCl of 0.01, 0.97 and 0.02, respectively, assuming methane is the only gas present. Using these data, an homogenization pressure of less than 1000 bar was obtained using the equation of state of Duan et al. [41].

As can be observed in Figure 6, type 2b fluid inclusions have a wide range of salinities from less than 10 wt.% to approximately 25 wt.% eq. NaCl. A variation in salinity could be produced by progressive separation of a methane-rich fluid during cooling, consumption of water during serpentinization, or mixing between more and less saline fluids. Fractional phase separation would produce a distribution of data points moving (approximately) away from the CH4 apex, consumption of water due to serpentinization would produce a distribution of data points moving away from the H2O apex, and mixing of two aqueous fluids would produce a line parallel to the H2O-NaCl tie-line provided that both fluids contained the same proportion of methane. The observed distribution is most consistent with mixing accompanied by a minor loss of CH4 (Figure 8).

Figure 8
figure 8

Ternary diagram showing the compositions of type 2b fluid inclusions (solid circles) in the system H2O-NaCl-CH4. The compositions were calculated with the program CURVES [30] using visual estimates of the volume fraction of aqueous solution in the inclusions. Dashed lines represent the solvus in the system at various pressures and temperatures (calculated using GEOFLUIDS model 1 [32,33,35-37]) along a geotherm of ~0.12°C/bar. The trend in inclusion fluid compositions suggests mixing between two fluids of contrasting salinity. See text for a detailed discussion.

Discussion

Trapping conditions

The most probable conditions of trapping of the inclusion fluids is the P-T region between the homogenization pressures, isochores limiting maximum and minimum pressures, and, for rodingitized diorite, granite and slate, the upper thermal stability of prehnite, which provides a maximum estimate of the temperature. The isochores for type 2a inclusions in samples Gran-1b and Slate-1a, and for type 2b inclusions in sample Dior-2 were calculated using the program GEOFLUIDS Model 1 and bulk molar volumes were estimated as described above. Although the P-T conditions associated with the earlier phase of rodingitization of the slates could not be determined, the lower homogenization temperatures of type 1 inclusions in sample Slate-2 suggest that they were trapped at pressures comparable to those trapped in diorite.

None of the mineral assemblages in the samples studied define invariant or univariant reactions in P-T space, for which the position can be estimated by correcting activities of end-member components for the effects of solid solution. There are also no suitable mineralogical geothermometers for these rocks. However, the upper boundary of the field of stability of prehnite can be used to constrain maximum trapping temperatures, as expressed by the decomposition/dehydration reaction

5Prehnite = 3Quartz + 2Grossular + 2Zoisite (or Clinozoisite) + 4H2O

which has a negative slope in P-T space. The position of this univariant reaction was determined with the program package THERMOCALC v2.3 [42, 43] assuming pure phases. Because grossular and zoisite do not have end-member compositions, and quartz is not present in the rodingites, the above phase boundary provides a maximum P-T limit for the projection of the isochores.

The position of the isochores for type 2 fluid inclusions in samples Dior-2, Gran-1b and Slate-1a, and of reaction (1) are shown on the P-T diagrams in Figure 9A and 9B. Also shown on this figure are the positions of the isochores for type 1 inclusions in zoisite and grossular from sample Dior-2 and in clinopyroxene from sample Slate-1a. The position of the isochores for type 1 inclusions was calculated using program ISOC [30, 31] and the equation of state for pure CH4 of Duan et al. [32, 33]. Volume corrections were performed for the compressibility and expansion of the clinopyroxene hosting type 1 fluid inclusions in sample Slate-1a.

Figure 9
figure 9

Pressure-temperature diagrams showing calculated trapping conditions for type 2b fluid inclusions in sample Dior-2 (A) and type 2a fluid inclusions in samples Gran-1b and Slate-1a (B). (A): The solid lines marked 1 and 2 (in white circles) represent P-T conditions at saturation (CH4-rich fluid) for type 2b fluid inclusions with salinities of 18 and 24 wt.% eq. NaCl, respectively. The heavy long-dashed line represent the average isochore of type 1 fluid inclusions occurring in sample Dior-2, and the short-dashed lines form the ± 1σ envelope of the average isochore. Stars indicate the maximum P-T conditions of equilibration of calc-silicate layers from the Thetford Mines ophiolite calculated by Laird et al. [50]. Because quartz is absent from their reported assemblage, the actual conditions should plot at lower temperatures and pressures along the univariant reaction 5Pumpellyite = 5Zoisite + 2Prehnite + 2Grossular + Chlorite + 9H2O which approximately parallels the isochores. (B): The grey shaded area represents the calculated trapping conditions for fluid inclusions in granite sample Gran-1b. The filled dot represents the calculated total homogenization conditions for a representative fluid inclusion in rodingitized slate sample Slate-1a. The solid line connected to the filled dot represents the calculated isochore of the corresponding fluid inclusion. Dashed lines represent isochores of type 1 fluid inclusions in sample Slate-1a. The occurrence of prehnite that formed contemporaneously with other calc-silicates in the studied rodingitized dykes and slate imposes the limit on homogenization pressures.

As noted in the section discussing the gas chromatographic analyses, bulk fluid samples from rodingitized diorite, granite and slate contain an unidentified species that may form up to 9 mol % of the total volatiles excluding H2O if this species is N2. To test the potential error in estimation of trapping pressures introduced by assuming that only CH4 is present in type 1 inclusions, we used the program BULK and the equation of state of Bakker [44] to calculate the molar volume of various mixtures of CH4 with N2 and H2, and compared them with results obtained for pure CH4. A value of -118°C was selected for the homogenization temperature, which represents the average for type 1 inclusions in sample Dior-2. The molar volumes were then used to calculate isochores using program ISOC. Our calculations suggest that trapping pressures would be overestimated by about 80–90 bar between 300 and 350°C for a fluid mixture with molecular proportions of 0.9 CH4 and 0.1 N2.

Molecular hydrogen is a common constituent of aqueous and gas phases in the serpentinization environment, and is produced during oxidation of iron [4549]. Assuming the molecular proportions of CH4 and H2 in the fluid mixture trapped by type 1 inclusions were 0.9 and 0.1, respectively, the results of our calculations suggest that the trapping pressures would be underestimated by 85–110 bar between 300 and 350°C.

Potential errors involved in the calculation of the isochores in type 2a and type 2b fluid inclusions due to the presence of a gas other than CH4 were also evaluated. Using data presented in Table 4, new fluid compositions were calculated using program CURVES and the equation of state of Duan et al. [37]. For this purpose, a mixture in the non aqueous phase with molecular proportions of 0.9 CH4 and 0.1 N2 was used. For each measurement, the volume fraction of aqueous solution after final clathrate melting was kept the same as that determined for fluids for which CH4 was the only dissolved gas species. The results indicate that the pressures may be underestimated by a maximum of approximately 80 bar for the mixture considered. Unfortunately, the method could not be applied to evaluate the effect of mixing H2 with CH4 because suitable equations of state for modeling clathrate involving a mixture of these volatiles in saline aqueous inclusions are not available. However, by analogy to the calculations involving 0.9 CH4-0.1 N2 and 0.9 CH4-0.1 H2 mixtures for type 1 inclusions, and 0.9 CH4-0.1 N2 for type 2 inclusions, the change in bulk molar volume of type 2 inclusions produced by a progressive increase in the proportion of H2 is probably small. We consider that an error in the order of ± 100 bar is a reasonable approximation on the position of the isochores for type 1 and type 2 fluid inclusions.

Despite the compositional uncertainties considered above, it is reasonable to assume, based on the petrographic observations and the estimation of the position of the isochores, that type 1 and type 2b inclusions in Dior-2 were trapped in the same environment. For these two types of inclusions to have been trapped simultaneously, a process involving immiscibility must be considered. To verify this, we calculated the composition of the low density, CH4-rich fluid phase that should coexist with type 2a inclusions at saturated vapour pressure (Ph in Table 4) using GEOFLUIDS. The calculated compositions were then used to calculate the volumes of water and methane in the low density phase at room temperature using BULK. The results indicate that between 4 and 8 vol. % of aqueous solution should be present in type 1 inclusions for the investigated total homogenization temperature interval of between 301 and 333°C for type 2a inclusions. Such small proportions of aqueous fluid would form a film that may appear invisible in the dark coloured type 1 inclusions. The homogenization temperature calculated for the methane in the CH4-rich, low density fluid coexisting with the water-rich fluid having a composition corresponding to that calculated for the lower salinity type 2a fluid inclusion at 301°C, as reported in Table 4, is -117.5°C. This value agrees well with those measured for type 1 fluid inclusions in sample Dior-2. For this to be the case, the trapping pressure of type 1 and type 2a inclusions in this sample would be equal to the total homogenization pressure of type 2a inclusions, i.e., between ~2.5 and 3 kbar. In sample Dior-1, which does not contain type 1 fluid inclusions, the trapping pressure of type 2a inclusions would be above that of the total homogenization.

Isochores for type 1 inclusions in Slate-1a cross the calculated P-T field of type 2a inclusions in sample Gran-1b and the calculated isochore for type 2a inclusions in sample Slate-1a. These inclusions are considered to have been trapped in a different environment to that in which type 1 inclusions in rodingitized diorite were trapped.

Based on the data analysis presented above, we estimate that pressure conditions accompanying rodingitization of diorite (samples Dior-1 and Dior-2) were between 2.5 and 4.5 kbar, and temperatures were between the minimum total homogenization temperature (286°C) and 360°C, the temperature at which the isochore projected from the highest homogenization temperature crosses reaction (1). These conditions are similar to those estimated by Laird et al. [50] from a study of calc-silicate layers in serpentinized ultramafic north of Richmond, Québec. Based on thermobarometric calculations using TWEEQU [51] on the assemblage pumpellyite + zoisite + prehnite + grossular + chlorite (quartz absent), they proposed maximum pressures between 3 and 5 kbar and maximum temperatures between 300 and 400°C. Pressures accompanying the formation of veins cutting pervasively rodingitized slate (sample Slate-1a) and rodingitization of the granite (sample Gran-1b) are interpreted to have been between 0.6 and 3 kbar; the temperature was potentially as high as 410°C.

The relatively high temperature of rodingitization of the granite raises the possibility that deserpentinization may have occurred, albeit locally, by reaction of serpentine with brucite to form forsteritic olivine by the reaction

Chrysotile + Brucite = 2Forsterite + 3H2O

or that temperature was buffered by this reaction. This is suggested by the fact that this reaction boundary is crossed by the isochore for sample Slate-1a and the lower temperature isochore of the P-T field for sample Gran-1b, whereas the higher temperature isochore for sample Gran-1b lies entirely in the forsterite field (Figure 9B; note that for natural forsteritic olivine from ultramafic rocks, temperatures for reaction (2) are shifted by 10°C to the left; Moody [52]).

Because the fluid inclusions in the vesuvianite samples contain significant proportions of an unknown species, no attempt was made to estimate the entrapment pressure of type 2c inclusions. However, if the two type 1 inclusions observed in vesuvianite were trapped at the same time as the type 2c inclusions, their homogenization temperature of -90°C would suggest a minimum pressure of entrapment of ~1 kbar.

Evolution of mineral parageneses

The fluid inclusion and petrographic data indicate that there were probably three episodes of rodingitization. The earliest episode corresponded to rodingitization of diorite and probably of slate at high pressure, and was followed by a second episode of rodingitization during which there was rodingitization of granite and further rodingitization of slate. The vesuvianite-rich "veins" formed during the latest stages of ophiolite emplacement. The data thus define an early high pressure, and a later, lower pressure regime of alteration. The pressure and temperature conditions estimated from type 2b (and coexisting type 1) fluid inclusions are consistent with stability fields calculated using the GEO-CALC [5355] and THERMOCALC v2.3 [42, 43] software packages for the assemblages with which they occur in the system CMASH. The interpreted trapping temperatures of fluid inclusions in sample Gran-1b (340 – 400°C; supported by occurrence of the assemblage hydrogrossular + clinopyroxene + wollastonite + vesuvianite in sample Gran-2c, which has a lower thermal stability of ~400°C at 1 to 2 kbar [5658]) and the total homogenization temperatures of "vein" sample Ves-2 (Th 180–203°C) suggest that vesuvianite formed over a very wide range of temperatures.

Because the granites must have intruded the serpentinites before rodingitization of the slates, the lower temperatures of the vesuvianite-rich "veins" must indicate a very late stage of rodingitization, recrystallization and intense veining. If the unidentified gas species was N2, its presence in the vesuvianite-rich "veins" in proportions up to twice those in the other rodingites might indicate a progressive increase in the contribution of a fluid derived from sedimentary rocks of the Laurentian continental margin or of meteoric waters.

Spatial and temporal relationships

The emplacement history of the Québec Appalachian ophiolites is complex and a clearer picture of this history has emerged from recent geological and geochronological studies that are consistent with the obduction history of other ophiolites that evolved in similar tectonic settings throughout the world. It is now well-established, based on U-Pb zircon age dating of oceanic plagiogranites and granitic rocks of continental affinity intruded in the ophiolites, on Ar-Ar cooling age determinations of amphibolite soles at the basal thrust of the ophiolites, and on continental-margin metasediments in tectonic contact with the ophiolites, that the latter were obducted very soon after their formation and that the granites were emplaced during obduction on the Laurentian margin during the Ordovician [5961].

At the suggested time of intrusion of the felsic intrusive rocks through the base of the ophiolite sequences, between 475 and 465 Ma [59], the only rock unit that existed above the ophiolite was the Saint-Daniel Formation, which comprises, among other lithologies, pelagic sediments deposited on the seafloor [21, 62]. Southeast of Asbestos, this formation is < 6 km thick. The combined thickness of this sedimentary unit (assuming it represents the original thickness) and the underlying igneous crust (2 km) indicates that the minimum depth of emplacement of the felsic intrusive rocks in the Asbestos ophiolite was 8 km, corresponding to a lithostatic pressure of > 2.2–2.4 kbar, i.e., just below the range of pressures estimated from the fluid inclusion data in rodingitized diorite.

Our calculated P-T conditions for the trapping of fluids in rodingitized diorite are intersected by the hotter P-T-t paths estimated by Whitehead et al. [59] during obduction of the Thetford-Mines ophiolite and generation of granitic rocks by melting of rock units underlying the young and hot mantle rocks at the base of the thrust. We suggest that P-T conditions similar to those for the Thetford-Mines ophiolite developed during obduction of the Asbestos ophiolite and that the high density fluids in type 1 and type 2a inclusions were responsible for rodingitization of diorite during this time (Figure 10).

Figure 10
figure 10

Pressure-temperature diagram showing the calculated trapping conditions of the fluid inclusions in rodingites from the JM Asbestos mine and superimposed pressure-temperature-time (P-T-t) paths estimated by Whitehead et al. [59] for the footwall at the base of the Thetford-Mines ophiolite (dashed lines). Only the paths that cross the P-T conditions determined by Clague et al. [65] and Feininger [66] for the metamorphism in the footwall at the base of the Thetford-Mines ophiolite are reproduced. The P-T-t paths define a range of geotherms that transect the calculated trapping conditions of type 2b fluid inclusions in rodingitized diorite. Trapping conditions for type 2a fluid inclusions in samples Gran-1b and Slate-1a, and for two type 1 fluid inclusions in vesuvianite (isochore indicated by solid line labelled 'Ves'), lie below those geotherms. The diagram also show the limits of stability of forsterite and antigorite in the system MgO-SiO2-H2O, and prehnite, in part, calculated using THERMOCALC v2.3 [42,43] assuming pure phases. The wet pelite solidus is after Thompson and Algor [67] and the limit for the stability of andalusite after Richardson et al. [68]. The diagram suggests emplacement of andalusite granites at relatively shallow levels in the Asbestos ophiolite, and possibly through peridotites that were already serpentinized.

The similar trapping temperatures but lower trapping pressures determined for the fluids in the rodingitized granite and slate samples imply that rodingitization also occurred near the final stages of emplacement and deformation of the Asbestos ophiolite. The relatively low, greenschist grade metamorphism that affected the slates at the base of the Asbestos ophiolite, the fact that these slates are rodingitized, and the evidence for thrusting followed by normal movement along the basal contact (Alain Tremblay, pers. comm. 2005), all suggest that the saline and Ca-rich fluids in type 2b inclusions interacted with the host rocks later than was the case for the less saline fluids represented by type 2a inclusions. The kinematic indicators for both thrust and normal faulting (striations) are superimposed on the blackwall replacing the slates between the rodingite and the serpentinized ultramafic rocks, and on the serpentinites at the contact.

A precise age for this event can not be given. However, based on the geological evidence presented above and 40Ar/39Ar dating of muscovite in the Caldwell sediments at the base of the Asbestos ophiolite [61], we can state confidently that type 2b fluids were trapped after granite intrusion at 465–475 Ma., and between Taconic metamorphism at 461–468 Ma. and extensional deformation which followed during the Silurian at 430–410 Ma. [62, 63].

Conclusion

Petrographic and fluid inclusion data indicate that there were at least two, and probably three, episodes of rodingitization involving distinct fluids. In an early, relatively high-pressure (~2.5–4.5 kbar) episode, fluids of moderate to high salinity (~1.4 to 5.7 m eq. NaCl) rodingitized diorite at temperatures of 290–360°C. A second episode of rodingitization at overlapping temperatures (325–400°C), but lower pressure, is recorded by primary, high salinity (~5.3 to 6.3 m eq. NaCl) fluid inclusions in veins cutting rodingitized slate and in rodingitized granite. Primary, low-salinity (1.5 m eq. NaCl), low temperature (< 200°C) fluid inclusions in vesuvianite-rich, highly porous mineral assemblages were the latest to form. Gas chromatographic analyses of the inclusion fluids in the rodingites indicate that all contain a carbonic component, in which CH4 is the dominant species among mixtures containing short-chain alkanes. This suggests that all the fluids interacted with serpentinite in Fischer-Tropsch and related hydrocarbon synthesizing reactions. It is interesting to note that type 2a and type 2b fluid inclusions represent the most concentrated brines yet reported from rodingites (cf, Compagnoni et al. [64]). Another interesting feature of the fluids is that they present a rare example of first ice melting temperatures consistent with divalent cation-rich compositions (the only other example of this that we were able to find was for zircon-hosted fluid inclusions in a metasomatic shell around rodingite from the Jordanów-Gogołów serpentinite massif, Poland [13]). Fluids responsible for rodingitization are thus not necessarily Ca-rich.

The exact nature of the inclusion fluids in the rodingites can not be resolved until more detailed analyses of their dissolved salts are performed and C-O-H stable isotopic data are obtained. Nevertheless, some general observations can be made taking into account the geological environment in which the Asbestos ophiolite evolved. The most important information concerning the fluid inclusions is that type 1, type 2a and type 2b inclusion fluids occur both in the rodingitized slates and in the rodingitized felsic dykes. Based on this, we conclude that the rodingites formed during thrusting of the Asbestos ophiolite onto the Laurentian margin. Thus, we rule out the involvement in rodingitization of seawater modified in a hydrothermal system driven by mafic magmatism in or near an oceanic zone of extension. For the same reason we also rule out dehydration of deeply subducted oceanic crust and sediments as a source for the rodingitizing fluid. Potential candidates for rodingitizing fluids represented by the moderately to strongly saline type 2b and 2a fluid inclusions of the Asbestos ophiolite are metamorphic waters, including evolved seawater, and evolved meteoric and connate waters.

References

  1. Bell JM, Clarke EdeC, Marshall P: The geology of the Dun Mountain subdivision, Nelson. 1911, New Zealand Geological Survey Bulletin 12 new series

    Google Scholar 

  2. Coleman RG: New Zealand serpentinites and associated metasomatic rocks. 1966, New Zealand Geological Survey Bulletin 76

    Google Scholar 

  3. Coleman RG: Low-temperature reaction zones and alpine ultramafic rocks of California, Oregon and Washington. 1967, United States Geological Survey Bulletin 1247

    Google Scholar 

  4. Honnorez J, Kirst P: Petrology of rodingites from the equatorial Mid-Atlantic fracture zones and their geotectonic significance. Contributions to Mineralogy and Petrology. 1975, 49: 233-257. 10.1007/BF00376590.

    Article  Google Scholar 

  5. Anhaeusser CR: Rodingite occurrences in some Archean ultramafic complexes in the Barberton Mountain Land, South Africa. Precambrian Research. 1979, 8: 49-76. 10.1016/0301-9268(79)90038-X.

    Article  Google Scholar 

  6. Wares RP, Martin RF: Rodingitization of granite and serpentinite in the Jeffrey mine, Asbestos, Quebec. Canadian Mineralogist. 1980, 18: 231-240.

    Google Scholar 

  7. Barriga F, Fyfe WS: Development of rodingite in basaltic rocks in serpentinites, East Liguria, Italy. Contributions to Mineralogy and Petrology. 1983, 84: 146-151. 10.1007/BF00371281.

    Article  Google Scholar 

  8. Katoh T, Niida K: Rodingites from the Kamuikotan, tectonic belt, Hokkaido. Journal of the Faculty of Science, Hokkaido University, Serie IV. 1983, 20: 151-169.

    Google Scholar 

  9. Schandl E, O'Hanley DS, Wicks FJ: Rodingites in serpentinized ultramafic rocks of the Abitibi greenstone belt, Ontario. Canadian Mineralogist. 1989, 27: 579-591.

    Google Scholar 

  10. Schandl E, O'Hanley DS, Wicks FJ, Kyser TK: Fluid inclusions in rodingite: a geothermometer for serpentinization. Economic Geology. 1990, 85: 1273-1276.

    Article  Google Scholar 

  11. O'Hanley DS, Schandl ES, Wicks FJ: The origin of rodingites from Cassiar, British Columbia, and their use to estimate T and P(H2O) during serpentinization. Geochimica et Cosmochimica Acta. 1992, 56: 97-108. 10.1016/0016-7037(92)90119-4.

    Article  Google Scholar 

  12. Dubiñska E: Rodingites of the eastern part of Jordanów-Gogołów serpentinite massif, Lower Silesia, Poland. Canadian Mineralogist. 1995, 33: 3585-3608.

    Google Scholar 

  13. Dubiñska E, Bylina P, Kozlowski A, Dörr W, Nejbert K, Schastok J, Kulicki C: U-Pb dating of serpentinization: hydrothermal zircon from metasomatic rodingite shell (Sudetic ophiolite, SW Poland). Chemical Geology. 2004, 203: 183-203. 10.1016/j.chemgeo.2003.10.005.

    Article  Google Scholar 

  14. Frost BR: Contact metamorphism of serpentinite, chloritic blackwall and rodingite at Paddy-Go-Easy Pass, Central Cascades, Washington. Journal of Petrology. 1975, 16: 272-313.

    Article  Google Scholar 

  15. Mittwede SK, Schandl ES: Rodingites from the southern Appalachian Piedmont, South Carolina, USA. European Journal of Mineralogy. 1992, 4: 7-16.

    Article  Google Scholar 

  16. Schandl ES, Mittwede SK: Evolution of the Acipayam (Denizli, Turkey) rodingites. International Geology Review. 2001, 43: 611-623.

    Article  Google Scholar 

  17. Hébert R, Bédard JH: Architecture and petrogenesis of arc ophiolites: example from the ophiolitic complexes of southern Quebec. Geological Association of Canada, Mineralogical Association of Canada, Association des Géologues et Géophysiciens du Québec, International Association of Hydrogeologists, Canadian Geophysical Union Joint Annual Meeting, Québec, Québec, Canada, Fieldtrip A4 Guidebook. 1988

    Google Scholar 

  18. Bédard JH, Shroetter J-M, Pagé P, Tremblay A, Bécu V: An overview of the geology and Cr-PGE potential of the southern Québec ophiolite belt. Mineral Resources of Canada: A Synthesis of Major Deposit-types, District Metallogeny, the Evolution of Geological Provinces, and Exploration Methods (Goodfellow WD, ed.): Geological Association of Canada, Mineral deposits Division, Special Publication x.

  19. Hébert R: Étude pétrologique des roches ophiolitiques d'Asbestos et du Mont Ham (Ham Sud). 1980, M. Sc. Thesis, Université Laval, Québec

    Google Scholar 

  20. Cousineau PA: Le Groupe de Caldwell et le domaine océanique entre Saint-Joseph-de-Beauce et Sainte-Sabine. Ministère de l'Énergie, des Mines et des Ressources du Québec, MM 87-02. 1990

    Google Scholar 

  21. Pinet N, Tremblay A: Tectonic evolution of the Quebec-Maine Appalachians: From oceanic spreading to obduction and collision in the Northern Appalachians. American Journal of Science. 1995, 295: 173-200.

    Article  Google Scholar 

  22. Normand C: Experimental and field investigations of serpentinization and rodingitization. 2001, Ph.D. thesis, McGill University, Montreal

    Google Scholar 

  23. Laurent R, Hébert Y: Paragenesis of serpentine assemblages in harzburgite tectonite and dunite cumulate from the Quebec Appalachians. Canadian Mineralogist. 1979, 17: 857-869.

    Google Scholar 

  24. De A: Post-ultramafic dykes and associated rocks in the asbestos mining areas in the Eastern Townships, Quebec. 1968, Geological Survey of Canada, Report of Activities, Paper 68-1, Part A

    Book  Google Scholar 

  25. Davis DW, Lowenstein TK, Spencer RJ: Melting behavior of fluid inclusions in laboratory-grown halite crystals in the systems NaCl-H2O, NaCl-KCl-H2O, NaCl-MgCl2-H2O and NaCl-CaCl2-H2O. Geochimica et Cosmochimica Acta. 1990, 54: 591-601. 10.1016/0016-7037(90)90355-O.

    Article  Google Scholar 

  26. Samson IM, Walker RT: Cryogenic Raman spectroscopic studies in the system NaCl-CaCl2-H2O and implications for low-temperature phase behaviour in aqueous fluid inclusions. Canadian Mineralogist. 2000, 38: 35-43. 10.2113/gscanmin.38.1.35.

    Article  Google Scholar 

  27. Bakker RJ: Raman spectra of fluid and crystal mixtures in the systems H2O, H2O-NaCl and H2O-MgCl2 at low temperatures: applications to fluid-inclusion research. Canadian Mineralogist. 2004, 42: 1283-1314. 10.2113/gscanmin.42.5.1283.

    Article  Google Scholar 

  28. Williams-Jones AE, Samson IM: Theoretical estimation of halite solubility in the system NaCl-CaCl2-H2O: applications to fluid inclusions. Canadian Mineralogist. 1990, 28: 299-304.

    Google Scholar 

  29. Naden J: CalcicBrine: A Microsoft Excel 5.0 add-in for calculating salinities from microthermometric data in the system NaCl-CaCl2-H2O. 1996, PACROFI VI, Madison, USA, 97-98.

    Google Scholar 

  30. Bakker RJ, Brown P: Computer modelling in fluid inclusion research. Fluid Inclusions: Analysis and Interpretation (Samson I, Anderson A, Marshall D, eds.), Mineralogical Association of Canada Short Course. 2003, 32: 175-212.

    Google Scholar 

  31. Bakker RJ: Package FLUIDS 1. New computer programs for the analysis of fluid inclusion data and for modelling bulk fluid properties. Chemical Geology. 2003, 194: 3-23. 10.1016/S0009-2541(02)00268-1.

    Article  Google Scholar 

  32. Duan Z, Møller N, Weare J: An equation of state for the CH4-CO2-H2O system: I. Pure systems from 0 to 8000°C and 0 to 8000 bar. Geochimica et Cosmochimica Acta. 1992, 56: 2605-2617. 10.1016/0016-7037(92)90347-L.

    Article  Google Scholar 

  33. Duan Z, Møller N, Weare J: An equation of state for the CH4-CO2-H2O system: II. Mixtures from 50 to 1000°C and from 0 to 1000 bar. Geochimica et Cosmochimica Acta. 1992, 56: 2619-2631. 10.1016/0016-7037(92)90348-M.

    Article  Google Scholar 

  34. Bakker RJ: Clathrates: Computer programs to calculate fluid inclusion V-X properties using clathrate melting temperatures. Computers & Geosciences. 1997, 23: 1-18. 10.1016/S0098-3004(96)00073-8.

    Article  Google Scholar 

  35. Duan Z, Møller N, Weare J: Equation of state for the NaCl-H2O-CO2 system: Prediction of phase equilibria and volumetric properties. Geochimica et Cosmochimica Acta. 1995, 59: 2869-2882. 10.1016/0016-7037(95)00182-4.

    Article  Google Scholar 

  36. Duan Z, Møller N, Weare J: Accurate prediction of the thermodynamic properties of fluids in the system H2O-CO2-CH4-N2 up to 2000 K and 100 kbar from a corresponding states/one fluid equation of state. Geochimica et Cosmochimica Acta. 1996, 64: 1069-1075. 10.1016/S0016-7037(99)00368-3.

    Article  Google Scholar 

  37. Duan Z, Møller N, Weare J: A general equation of state for supercritical fluid mixtures and molecular dynamics simulation of mixture PVTX properties. Geochimica et Cosmochimica Acta. 1996, 60: 1209-1216. 10.1016/0016-7037(96)00004-X.

    Article  Google Scholar 

  38. Lamb WM, McShane CJ, Popp RK: Phase relations in the CH4-H2O-NaCl system at 2 kbar, 300 to 600°C as determined using synthetic fluid inclusions. Geochimica et Cosmochimica Acta. 2002, 66: 3971-3986. 10.1016/S0016-7037(02)00963-8.

    Article  Google Scholar 

  39. Lamb WM, Popp RK, Boockoff LA: The determination of phase relations in the CH4-H2O-NaCl system at 1 kbar, 400 to 600°C using synthetic fluid inclusions. Geochimica et Cosmochimica Acta. 1996, 60: 1885-1897. 10.1016/0016-7037(96)00070-1.

    Article  Google Scholar 

  40. Bodnar RJ: Revised equation and table for determining the freezing point depression of H2O-NaCl solutions. Geochimica et Cosmochimica Acta. 1993, 57: 683-684. 10.1016/0016-7037(93)90378-A.

    Article  Google Scholar 

  41. Duan Z, Møller N, Greenberg J, Weare J: The prediction of methane solubility in natural waters to high ionic strength from 0 to 250°C and from 0 to 1600 bar. Geochimica et Cosmochimica Acta. 1992, 56: 1451-1460. 10.1016/0016-7037(92)90215-5.

    Article  Google Scholar 

  42. Powell R, Holland TJB: An internally consistent thermodynamic dataset with uncertainties and correlations: 3. Applications to geobarometry, worked examples and a computer program. Journal of Metamorphic Geology. 1988, 6: 173-204. 10.1111/j.1525-1314.1988.tb00415.x.

    Article  Google Scholar 

  43. Holland TJB, Powell R: An enlarged and updated internally consistent thermodynamic dataset with uncertainties and correlations: the system K2O-Na2O-CaO-MgO-MnO-FeO-Fe2O3-Al2O3-TiO2-SiO2-C-H2-O2. Journal of Metamorphic Geology. 1990, 8: 89-124. 10.1111/j.1525-1314.1990.tb00458.x.

    Article  Google Scholar 

  44. Bakker RJ: Adaptation of the Bowers and Helgeson (1983) equation of state to the H2O-CO2-CH4-N2-NaCl system. Chemical Geology. 1999, 154: 225-236. 10.1016/S0009-2541(98)00133-8.

    Article  Google Scholar 

  45. Thayer TP: Serpentinization considered as a constant volume metasomatic process. American Mineralogist. 1966, 51: 685-710.

    Google Scholar 

  46. Neal C, Stanger G: Hydrogen generation from mantle source rocks in Oman. Earth and Planetary Science Letters. 1983, 66: 315-320. 10.1016/0012-821X(83)90144-9.

    Article  Google Scholar 

  47. Abrajano TA, Sturchio NC, Kennedy BM, Lyon GL, Muehlenbachs K, Bohlke JK: Geochemistry of reduced gas related to serpentinization of the Zambales ophiolite, Philippines. Applied Geochemistry. 1990, 5: 625-630. 10.1016/0883-2927(90)90060-I.

    Article  Google Scholar 

  48. Peretti A, Dubessy J, Mullis J, Frost BR, Trommsdorff V: Highly reducing conditions during Alpine metamorphism of the Malenco peridotite (Sondrio, northern Italy) indicated by mineral paragenesis and H2 in fluid inclusions. Contributions to Mineralogy and Petrology. 1992, 112: 329-340. 10.1007/BF00310464.

    Article  Google Scholar 

  49. Charlou JL, Donval JP, Fouquet Y, Jean Baptiste P, Holm N: Geochemistry of high H2 and CH4 vent fluids issuing from ultramafic rocks at the Rainbow hydrothermal field (36°14'N, MAR). Chemical Geology. 2002, 191: 345-359. 10.1016/S0009-2541(02)00134-1.

    Article  Google Scholar 

  50. Laird J, Trzcienski WE, Bothner WA: High pressure, Taconian, and subsequent polymetamorphism of Southern Quebec and Northern Vermont. Geological Society of America Field Trip Guidebook for the Northeastern United States Volume 2. Department of Geology and Geography, University of Massachusetts Contribution 67, O 1-O 32. 1993

    Google Scholar 

  51. Berman RG: Thermobarometry using multiequilibrium calculations: a new technique with petrologic applications. Canadian Mineralogist. 1991, 29: 833-855.

    Google Scholar 

  52. Moody JB: An experimental study on the serpentinization of iron-bearing olivines. Canadian Mineralogist. 1976, 14: 462-478.

    Google Scholar 

  53. Berman RG, Brown TH, Perkins EH: GEØ-CALC: Software for calculation and display of P-T-X phase diagrams. American Mineralogist. 1987, 72: 861-862.

    Google Scholar 

  54. Brown TH, Berman RG, Perkins EH: GEØ-CALC: a software package for rapid calculation of stable pressure – temperature – activity phase diagrams. Geological Society of America, Abstracts with Programs. 1987, 19: 603-

    Google Scholar 

  55. Brown TH, Berman RG, Perkins EH: GEØ-CALC: a software package for calculation and display of pressure – temperature – composition phase diagrams using an IBM or compatible personal computer. Computers and Geoscience. 1988, 14: 279-289. 10.1016/0098-3004(88)90061-1.

    Article  Google Scholar 

  56. Hochella MF, Liou JG, Keskinen MJ, Kim HS: Synthesis and stability relations of magnesium idocrase. Economic Geology. 1982, 77: 798-808.

    Article  Google Scholar 

  57. Valley JW, Peacor DR, Bowman JR, Essene EJ, Allard MJ: Crystal chemistry of a Mg-vesuvianite and implications of phase equilibria in the system CaO-MgO-Al2O3-SiO2-H2O-CO2. Journal of Metamorphic Geology. 1985, 3: 137-153. 10.1111/j.1525-1314.1985.tb00311.x.

    Article  Google Scholar 

  58. Hoisch TD: The solid solution chemistry of vesuvianite. Contributions to Mineralogy and Petrology. 1985, 89: 205-214. 10.1007/BF00379454.

    Article  Google Scholar 

  59. Whitehead J, Dunning GR, Spray JG: U-Pb geochronology and origin of granitoid rocks in the Thetford Mines ophiolite, Canadian Appalachians. Geological Society of America Bulletin. 2000, 112: 915-928. 10.1130/0016-7606(2000)112<0915:UPGAOO>2.3.CO;2.

    Article  Google Scholar 

  60. Whitehead J, Reynolds PH, Spray JG: The sub-ophiolitic metamorphic rocks of the Québec Appalachians. Journal of Geodynamics. 1995, 19: 325-350. 10.1016/0264-3707(94)00021-M.

    Article  Google Scholar 

  61. Whitehead J, Reynolds PH, Spray JG: 40Ar/39Ar age constraints on Taconian and Acadian events in the Quebec Appalachians. Geology. 1996, 24: 359-362. 10.1130/0091-7613(1996)024<0359:AAACOT>2.3.CO;2.

    Article  Google Scholar 

  62. Schroetter J-M, Tremblay A, Bedard JH, Villeneuve ME: Syncollisional basin development in the Appalachian Orogen; the Saint-Daniel Melange, southern Quebec, Canada. Geological Society of America Bulletin. 2006, 118: 109-125. 10.1130/B25779.1.

    Article  Google Scholar 

  63. Tremblay A, Castonguay S: The structural evolution of the Laurentian margin revisited (southern Quebec): implications for the Salinian Orogeny and Appalachian successor basins. Geology. 2002, 30: 79-82. 10.1130/0091-7613(2002)030<0079:SEOTLM>2.0.CO;2.

    Article  Google Scholar 

  64. Compagnoni A, Conte R, Facchinelli A, Orione P, Rossetti P, Trossarelli C: The rodingite of Bellecombe (Val D'Aosta, Italy): a petrographic and fluid inclusion study. 17th general meeting, International Mineralogical Association, August 9–14, Toronto, Canada, Programme and Abstracts, A77L. 1998

    Google Scholar 

  65. Clague D, Rubin J, Brackett R: The age and origin of the of the garnet amphibolite underlying the Thetford mines ophiolite, Quebec. Canadian Journal of Earth Sciences. 1981, 18: 469-486.

    Article  Google Scholar 

  66. Feininger T: Amphibolite associated with the Thetford Mines ophiolite complex at Belmina Ridge, Quebec. Canadian Journal of Earth Sciences. 1981, 18: 1878-1892.

    Article  Google Scholar 

  67. Thompson AB, Algor JR: Model systems for anatexis of politic rocks I. Theory of modeling reactions in the system KAlO2-NaAlO2-Al2O3-SiO2-H2O. Contributions to Mineralogy and Petrology. 1977, 63: 247-269. 10.1007/BF00375575.

    Article  Google Scholar 

  68. Richardson SW, Gilbert MC, Bell PM: Experimental determination of kyanite-andalusite and andalusite-sillimanite equilibria; the aluminium silicate triple point. America Journal of Science. 1969, 267: 259-272.

    Article  Google Scholar 

Download references

Acknowledgements

Prokopis Kranidiotis is thanked for collecting part of the fluid inclusion data set. The comments of two reviewers helped to improve this paper. This research was funded by NSERC and FCAR scholarships to CN and NSERC and FCAR research grants to AEW-J.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Charles Normand.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Normand, C., Williams-Jones, A.E. Physicochemical conditions and timing of rodingite formation: evidence from rodingite-hosted fluid inclusions in the JM Asbestos mine, Asbestos, Québec. Geochem Trans 8, 11 (2007). https://doi.org/10.1186/1467-4866-8-11

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1467-4866-8-11

Keywords