Very Low- to Low-grade Metamorphic Processes

JOURNAL OF PETROLOGY
VOLUME 54
NUMBER 9
PAGES 1849^1874
2013
doi:10.1093/petrology/egt033
Very Low- to Low-grade Metamorphic Processes
Related to the Collisional Assembly of Avalonia
in SE Cape Breton Island (Nova Scotia, Canada)
ARNE P. WILLNER1,2*, HANS-JOACHIM MASSONNE1,
SANDRA M. BARR3 AND CHRIS E. WHITE4
1
INSTITUT FU«R MINERALOGIE UND KRISTALLCHEMIE, UNIVERSITA«T STUTTGART, AZENBERGSTR. 18, 70174
STUTTGART, GERMANY
2
INSTITUT FU«R GEOLOGIE, MINERALOGIE UND GEOPHYSIK, RUHR-UNIVERSITA«T, 44780 BOCHUM, GERMANY
3
DEPARTMENT OF EARTH AND ENVIRONMENTAL SCIENCE, ACADIA UNIVERSITY, WOLFVILLE, NS, B4P 2R6, CANADA
4
NOVA SCOTIA DEPARTMENT OF NATURAL RESOURCES, PO BOX 698, HALIFAX, NS, B3J 2T9, CANADA
RECEIVED MAY 2, 2012; ACCEPTED MAY 13, 2013
ADVANCE ACCESS PUBLICATION JULY 3, 2013
To better understand geodynamic processes related to the assembly
of various belts of volcanic^sedimentary rocks in the Avalonian
Mira terrane of SE Cape Breton Island (protolith ages 680^560
Ma) we investigated metamorphic processes in 20 white-mica-bearing mafic and felsic metavolcanic rocks. The felsic metavolcanic
rocks are foliated and partly sheared with blastoporphyritic relict
fabric and contain the assemblage phengite^epidote^chlorite^
albite^K-feldspar^quartz^titanite stilpnomelane calcite ilmenite. In contrast, many mafic metavolcanic rocks are relatively undeformed and display relict porphyritic and amygdaloidal textures.
They contain the assemblages epidote^chlorite^albite^quartz^titanite phengite pumpellyite prehnite calcite K-feldspar and
actinolite^epidote^chlorite^albite^quartz^titanite
phengite calcite K-feldspar. Heterogeneous metamorphic overprinting is indicated by local relicts of magmatic clinopyroxene, magnetite and plagioclase. Metamorphic minerals formed by local
precipitation in clusters and are due to continuous nucleation of very
low-grade phases during pulses of variably pervasive fluids, which
were released during intense dehydration at 250^3008C. Nucleaction
rate dominated over growth rate at these conditions. Potassic white
mica in both mafic and felsic rocks is mostly phengite with a wide compositional range (3·11^3·41 Si a.p.f.u.). Maximum Si contents are typically between 3·30 and 3·41a.p.f.u. P^T pseudosections were
calculated for the range 200^4008C and 1^7 kbar. The peak metamorphic assemblages occupy P^T fields consistent with the position of
*Corresponding author.Telephone: þ49 71168581218. Fax: þ49 7116858
1222. E-mail: [email protected]
isopleths for corresponding maximum Si contents in white mica. Peak
P^Tconditions in the Mira terrane samples lie within a narrow range
of 3·5 0·4 kbar and 280 308C in samples representing all of the
assembled volcanic^sedimentary belts. The derived peak metamorphic
conditions suggest syn-collisional burial to 11^14 km depth and a low
metamorphic geotherm of 20^258C km^1. Under these conditions subsequent strike-slip deformation is attributed to the final assembly of
magmatic arc slices to form the crust of the Mira terrane.
very low-grade metamorphism; prehnite^actinolite
facies; phengite; P^T pseudosections; mineral growth kinetics;
collision; Avalonia
KEY WORDS:
I N T RO D U C T I O N
Partial pressure^temperature (P^T) paths are commonly
derived from a series of equilibria preserved in metamorphic rocks, taking advantage of the compositional
variation of the phases in the prevalent assemblage. This
approach is indispensable for understanding tectonics in
three dimensions, particularly in exhumed metamorphic
areas. Although metamorphic rocks of very low to low
grade (subgreenschist to greenschist facies) are widespread
on Earth, the metamorphic evolution of these rocks has
The Author 2013. Published by Oxford University Press. All
rights reserved. For Permissions, please e-mail: journals.permissions@
oup.com
JOURNAL OF PETROLOGY
VOLUME 54
not been intensively studied by equilibrium thermodynamics, in contrast to medium- and high-grade metamorphic rocks. Experiments that limited the P^T stability
of very low- to low-grade minerals (e.g. Liou, 1971a, 1971b;
Nitsch, 1971; Schiffman & Liou, 1980) have allowed P^Testimates for rocks containing those minerals. For instance,
Schreyer & Abraham (1978), Cho et al. (1986) and Cho &
Liou (1987) used these stability limits to estimate P^Tconditions during metamorphism of mafic to felsic rocks. In
addition, the experimentally determined stability of very
low- to low-grade metamorphic minerals was used to construct petrogenetic grids for low-temperature metamorphism (e.g. Brown, 1975; Liou et al., 1985, 1987; Frey et al.,
1991; Banno, 1998), which have also been used to constrain
P^T conditions (e.g. Robinson et al., 2005; McMullin
et al., 2010). With the extraction of thermodynamic data on
metamorphic minerals in general, datasets have been created (e.g. Berman, 1988; Holland & Powell, 1998) that also
include minerals of the very low- to low-grade metamorphic realm. Such datasets have been applied to estimate the P^T conditions of very low- to low-grade
metamorphic rocks. For instance, Massonne (1995) reconstructed a P^T path for mafic and felsic rocks of the
Taunus Mountains, central Europe, which had experienced peak temperatures at 3008C. This path was based
on the chemical variability of minerals, such as potassic
white mica and Na- and Ca-amphiboles in these rocks
and the augmentation of the thermodynamic dataset by
data for these minerals. Further examples of the thermodynamic treatment of rocks of various lithologies from the
very low- to low-grade metamorphic realm to construct
P^T paths were given by Dingeldey et al. (1997), Vidal &
Parra (2000), Willner et al. (2000, 2009), Parra et al. (2002),
Willner (2005), Vidal et al. (2006) and Jolivet et al. (2010).
After numerous P^T pseudosections were successfully
applied to medium- and high-grade metamorphic rocks,
Massonne & Willner (2008) and Massonne (2010) augmented the thermodynamic dataset by Holland & Powell
(1998, updated 2003) to calculate P^T pseudosections also
for very low-grade metamorphic rocks including metasedimentary rocks. Subsequently Willner et al. (2009), Cruciani
et al. (2011) and Kryza et al. (2011) estimated P^T paths
for various low-temperature metamorphic rocks on the
basis of P^T pseudosections. Applications of pseudosections at low grade so far have included a wide range of
lithologies such as mafic rocks, calcareous and non-calcareous metagreywacke, and metatrachyte. This approach
has the additional advantage that the P^T pseudosections
can also be used to understand phase relationships around
the derived peak P^Tconditions.
McMullin et al. (2010) studied very low- and low-grade
mafic metamorphic rocks in Neoproterozoic magmatic
arcs in southeastern Cape Breton Island, Nova Scotia,
part of the West Avalonia microplate. They showed by
NUMBER 9
SEPTEMBER 2013
chemography and calculation of some multivariant reactions that equilibria exist locally at very low-grade conditions and relatively low pressures in the transitional field
between prehnite^actinolite and greenschist facies. Here,
we broaden these studies by applying P^T pseudosections
at 20 mafic and felsic metavolcanic samples, including an
assessment of the growth behaviour of minerals at very
low grade. We also partly focus on white mica, which is
common in all rock types, but which was not included in
the study by McMullin et al. (2010). White mica in the samples is generally phengite, a mineral typically associated
with intermediate to high pressures. Hence, it is a challenge to show its compatibility with a relatively low-pressure environment. Our aims are to (1) provide better
understanding of processes at very low metamorphic temperatures where studies are challenged by fine grain size
and high chemical variability in the minerals, (2) improve
knowledge of the kinetics of metamorphic reactions at
very low grade, (3) show how transient equilibrium conditions are realized at very low grade, (4) show the potential
application of P^T pseudosections for geothermobarometry at very low- and low-grade conditions and (5) contribute to a better knowledge of geodynamics during
and after collisional processes. Results of dating of the
metamorphic processes and related deformation will be
presented elsewhere.
GEOLOGIC A L S ET T I NG
The Mira terrane of southeastern Cape Breton Island
(Fig. 1) consists mainly of Neoproterozoic to Cambrian
rocks typical of Avalonia (Barr & Raeside, 1989). Differences in late Neoproterozoic magmatic evolution led to
the recognition of five now-juxtaposed volcanic^sedimentary^plutonic belts: (1) Stirling; (2) East Bay Hills; (3) Coxheath Hills; (4) Sporting Mountain; (5) Coastal belt
(Barr, 1993; Bevier et al., 1993; Barr et al., 1996). The Stirling
belt contains metasedimentary and metavolcanic rocks
(Stirling Group) with a protolith age (U/Pb zircon) of 681
þ6/^2 Ma (Bevier et al., 1993), dioritic and granodioritic
plutons of the Chisholm Brook suite at 620 þ3/^2 Ma (U/
Pb zircon) and Devonian granitoid plutons (Barr & Macdonald, 1992; Fig. 1). The protoliths of the Stirling Group
were predominantly andesitic to basaltic lapilli tuffs, litharenites, and siltstones, with subordinate basaltic lava flows
and breccias, rhyolitic lapilli tuffs, and rhyolite porphyry,
intruded by comagmatic gabbroic dykes and sills and
including a massive sulphide deposit (Barr, 1993; Macdonald & Barr, 1993; Barr et al., 1996). Based on the chemical
composition of its volcanic and plutonic rocks, it probably
formed within a calc-alkaline volcanic arc, but under
marine conditions as evidenced by abundant epiclastic turbiditic deposits (Barr, 1993; Macdonald & Barr, 1993; Barr
et al., 1996).
1850
WILLNER et al.
VERY LOW- TO LOW-GRADE METAMORPHIC PROCESSES
Fig. 1. Simplified geological map of southeastern Cape Breton Island after Barr et al. (1996), Giles et al. (2010) and McMullin et al. (2010). Sample
locations are indicated. The inset map shows the location of the Mira terrane within Avalonia and the terrane assemblage of the northern
Appalachian orogen after Hibbard et al. (2006). A, Avalonia; G, Ganderia; M, Meguma; L, Laurentia; PL, peri-Laurentian margin.
The East Bay Hills, Coxheath Hills, and Sporting
Mountain belts (Fig. 1) consist mainly of subaerial metavolcanic and metavolcaniclastic rocks of basaltic, andesitic
and rhyolitic composition (East Bay Hills, Coxheath, and
Pringle Mountain groups) with protolith ages of c. 620 Ma
(U/Pb zircon; Bevier et al., 1993; Barr et al., 1996; White
et al., 2003). These belts also include c. 620 Ma dioritic to
granitic plutons considered to be comagmatic with the
host volcanic rocks (Barr et al., 1990; Bevier et al., 1993).
The three belts are considered to represent parts of a
single continental margin volcanic-arc complex (Barr
et al., 1996).
The Fourchu and Main-a'-Dieu groups form the Coastal
belt of the Mira terrane (Barr et al., 1996). The Fourchu
Group primarily consists of metamorphosed dacitic lapilli
and ash tuffs as well as rare lava flows with minor basaltic,
andesitic, and rhyolitic tuff and flow units. A tuffaceous
unit and a granitic pluton yielded igneous crystallization
ages of 574 1 Ma and 574 3 Ma, respectively (Bevier
et al., 1993). In contrast, the Main-a'-Dieu Group consists
dominantly of slightly metamorphosed tuffaceous sedimentary and epiclastic rocks with subordinate basaltic
and rhyolitic flows and tuffs with a maximum depositional
age of c. 563 Ma based on U^Pb dating of a rhyolite flow
(Bevier et al., 1993); its minimum age is constrained by
overlying lower Cambrian sedimentary rocks (Barr et al.,
1996). The two units of the Coastal belt are considered to
represent calc-alkaline volcanism in a continental margin
arc (Fourchu Group) and slightly younger intra-arc extension (Main-a'-Dieu Group). Both units are intruded by
Devonian plutons with widespread contact metamorphic
aureoles in both Neoproterozoic rocks and overlying
Cambrian units (Barr et al., 1996).
The Mira terrane locally displays subvertical folds of
variable scale (centimetre to kilometre; here summarized
as D1) and vergence that trend NE^SW and moderately
plunge to the NE. These folds, collectively assigned to D1,
vary in intensity and scale as a result of contrasts in rock
types within and among the volcanic-sedimentary-plutonic
belts (Macdonald & Barr, 1993; Barr et al., 1996). A widely
spaced and weak cleavage (S1) is associated with this deformation and is interpreted to represent an axial-planar
foliation. Deformation occurred prior to emplacement of
the Devonian plutons and after deposition of the
Cambrian sedimentary rocks; folding is probably polyphase and of different ages. The folded belts are dissected
by prominent shear zones (D2) of similar NE^SW strike,
particularly along the boundaries between the belts and
more penetratively in the Coastal and East Bay Hills
belts. Foliation (S2) in these shear zones, which have
1851
JOURNAL OF PETROLOGY
VOLUME 54
widths at a 100 m scale, dips steeply to the NW or SE. The
shear zones are predominantly sinistral strike-slip, characterized by conspicuous strain gradients at metre scale
from undeformed to mylonitic (Macdonald & Barr, 1993;
Barr et al., 1996), and were probably related to juxtaposition
of the various magmatic arcs. Folded (D1) strata of the
Stirling Group are intruded by the c. 620 Ma Chisholm
Brook suite, implying that some folding in that belt at
least is of late Neoproterozoic age (Barr et al., 1996).
Furthermore, McMullin et al. (2010) noted that many pyroclastic rocks contain schistose clasts. These observations
and weak local unconformities at the base of the Cambrian
section show that at least part of the folding occurred in
the late Neoproterozoic. The Cambrian rocks show no
penetrative deformation, but broad open folds. Some sandstone units contain detrital white mica with 40Ar/39Ar
ages of 630^550 Ma (Reynolds et al., 2009).
In the only published study of the metamorphism in the
Mira terrane, McMullin et al. (2010) estimated regional
metamorphic conditions transitional between the prehnite^actinolite and greenschist facies. Metamorphic conditions of 250^3508C and 53 kbar were deduced, and
ascribed to burial in a volcanic-arc setting. Rocks of both
facies occur in all belts and in close proximity to each
other. In the vicinity of plutons, amphibolite facies is attained locally as a result of contact metamorphism. Potter
et al. (2008a, 2008b) and Petts et al. (2012) reported a pervasive low-dO18 anomaly throughout the Mira terrane and
also in other parts of Avalonia. They interpreted this
anomaly as a result of pre-metamorphic hydrothermal
alteration during Late Neoproterozoic transtension of
Avalonia.
Devonian granites intruded at 378 þ5/^1 Ma (U/Pb
zircon; Barr & Macdonald, 1992; Bevier et al., 1993) to a
high level in the crust with some associated porphyry and
skarn-type mineralization. These plutons are probably
related to collision, with either Ganderia to the north or
the Meguma terrane to the south (Fig. 1 inset).
M ET HODS
Sampling
Diverse lithologies in all the belts of the Mira terrane were
sampled. Based on thin section studies of these samples
augmented by samples from the collection of McMullin
et al. (2010), we selected 20 samples of mafic and felsic
rocks for further petrological investigation, with preference
for those containing relatively large white mica grains
and representing the various belts of the Mira terrane
(Fig. 1). Samples potentially affected by contact metamorphism around Devonian plutons were avoided. The
qualitative modal compositions of the 20 selected samples
are given in Table 1.
NUMBER 9
SEPTEMBER 2013
Analytical methods
Mineral compositions were obtained using a CAMECA
SX 100 electron microprobe with five wavelength-dispersive systems at Universita«t Stuttgart, Germany. Operating
conditions were an acceleration voltage of 15 kV, a beam
current of 10 nA, 20 s counting time per element each on
the peak and the background, and preferably a slightly defocused beam of 5 mm to avoid loss of alkalis in mica and
amphibole. The following standards were used: natural
wollastonite (Si, Ca), synthetic periclase (Mg), synthetic
corundum (Al), synthetic rutile (Ti), natural hematite
(Fe), natural albite (Na), natural orthoclase (K), natural
rhodonite (Mn) and natural baryte (Ba). The PAP correction procedure provided by Cameca was used for matrix
corrections. Representative analyses and structural formulae of minerals used for P^T calculations, together with
the calculation procedure of the structural formulae, are
presented in Table 2. Detection limits given in Table 2
refer to the operating conditions used. Further data are
available as Supplementary Data, which can be downloaded from http://www.petrology.oxfordjournals.org and
upon request to the first author.
Whole-rock analyses for major elements were made on
fused glass disks obtained by melting rock powder and
Spectromelt in the proportion 1:9. The disks were analysed
using a Philips PW2400 X-ray fluorescence spectrometer
at Universita«t Stuttgart. CO2 and H2O were determined
by IR spectroscopy using a Leco RC-412 C/H/H2O analyzer. The analytical data used for the calculation of the
P^T pseudosections are given in Table 3.
Calculation methods
For a detailed assessment of phase relations, geothermobarometric constraints and equilibrium conditions in the
samples, we calculated P^T pseudosections in the system
K2O^Na2O^CaO^FeO^O2^MgO^Al2O3^TiO2^SiO2^
H2O at 200^4008C and 1^7 kbar using the PERPLE__X
software package (Connolly, 1990, 2005; version of August
2006 downloaded from www.perplex.ethz.ch). The
thermodynamic data of Holland & Powell (1998, updated
2003) for minerals and aqueous fluid were applied, with
the addition of end-member data for Fe2þ- and Fe3þ-pumpellyite, Fe2þ- and Mg-stilpnomelane, actinolite and magnesioriebeckite (Massonne & Willner, 2008). Calculations
were performed using the solid-solution models of Powell
& Holland (1999) and Holland & Powell (2003) for white
mica, epidote, chlorite and biotite and those of Massonne
& Willner (2008) for amphibole, sodic clinopyroxene,
pumpellyite and stilpnomelane. In relatively Na- and
Fe-rich bulk-rock compositions, preliminary calculations
indicated that end-member riebeckite may appear over a
wide P^T range. In these cases we used a modified solidsolution model for Na-amphibole as introduced by Kryza
et al. (2011). The clinopyroxene model is that of Holland &
1852
1853
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
Wm
x
x
x
x
x
x
x
x
x
(x)
x
x
Kf
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
Ab
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
Qz
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
x
Tt
x
x
x
x
x
x
St
x
x
x
x
x
Pr
x
x
x
Cc
x
x
x
x
x
x
x
x
x
x
Mag
x
x
x
Cp
x
cl
pc, po
pc, po
po
cl
po
pc
cl
po
pc
pc
po, am
po, am
po, am
pc
po, am
pc, po
pc
pc, am
po
fabric
x
Relic
Pl
An410
fy
fy
fy
fy
f
f
f
fy
fy
f
fy
nf
nf
nf
f
nf
nf
f
f
nf
Deformation
3·14–3·32
3·20–3·41
3·31–3·43
3·24–3·32
3·20–3·35
3·20–3·32
3·28–3·36
3·34–3·38
3·24–3·31
3·35–3·38
3·27–3·37
3·16–3·41
3·11–3·34
3·30–3·35
3·20–3·33
3·34
3·21–3·30
3·35–3·41
(a.p.f.u.)
Si in Wm
250–310
265–280
270–295
280–290
260–275
265–285
270–300
255–310
275–300
240–290
270–285
270–295
265–300
280–310
280–290
290–310
265–280
(8C)
Peak T
3·2–3·6
3·9
3·1–3·5
3·7–3·8
3·3–4·7
3·4–3·7
3·2–3·8
3·2–4·0
3·4–3·8
2·2–3·8
2·8–3·4
3·1–3·8
3·2–4·8
43·0–3·3
3·2–3·4
3·3–3·9
3·1–3·8
(kbar)
Peak P
Stirling Group
East Bay Hills Group
East Bay Hills Group
East Bay Hills Group
Fourchu Group
Fourchu Group
Fourchu Group
Fourchu Group
Fourchu Group
East Bay Hills Group
Fourchu Group
Stirling Group
Main-à-Dieu Group
Sporting Mtn. Group
East Bay Hills Group
Main-à-Dieu Group
Stirling Group
Fourchu Group
Fourchu Group
East Bay Hills Group
Lithostratigraphy
Mineral abbreviations as in Fig. 6 and Ilm, ilmenite; Mag, magnetite. f, foliated; nf, not foliated; po, porphyritic; pc, pyroclastic; am, amygdaloid; cl, compositional
layering; (x) calculated but not observed phase.
*Sample from McMullin et al. (2010).
yFoliation related to late sinistral strike slip.
x
10Ca44
x
10Ca39
x
x
10Ca42
x
x
10Ca34
10Ca37
(x)
x
x
10Ca33
x
x
x
x
x
x
M84-41b
x
x
10Ca11
x
x
x
x
x
x
x
x
x
x
M84-41c
x
x
x
x
M83-60
(x)
x
x
x
M83-57
(x)
(x)
FA8785*
Felsic rocks
(x)
FC881614*
AS467*
x
x
(x)
10Ca40
(x)
(x)
10Ca23
x
x
x
x
10Ca13
x
(x)
x
x
M84-42
M84-41a
(x)
M83-59
Mafic rocks
Ch
Ilm
Ep
Am
Pu
Primary phases
Metamorphic phases
Table 1: Assemblages in the studied rock types of the Mira terrane and their peak conditions
WILLNER et al.
VERY LOW- TO LOW-GRADE METAMORPHIC PROCESSES
1854
0·037
1·888
1·933
4·000
19·996
Na
K
Sum
OH
O*
4·100
Sum
0·030
0·417
Mg
0·008
0·020
Mn
Ba
0·519
Fe2þ
19·923
4·000
1·795
1·741
0·041
0·013
0·026
4·100
0·642
0·017
0·071
0·698
0·014
2·659
8·000
1·598
6·403
98·1
4·34
9·88
0·15
19·973
4·000
1·945
1·901
0·025
0·000
0·019
4·100
0·365
0·020
0·564
0·261
0·006
2·885
8·000
1·377
6·623
100·7
4·41
11·0
0·10
19·969
4·000
1·938
1·904
0·015
0·002
0·017
4·100
0·352
0·016
0·445
0·361
0·003
2·924
8·000
1·511
6·490
99·7
4·37
10·9
0·06
0·32
50·07
0·36
50·07
19·927
4·000
1·818
1·774
0·020
0·024
0·007
4·082
0·542
0·005
0·498
0·000
0·011
3·027
8·000
1·213
6·787
100·2
4·48
10·4
0·08
0·12
19·943
4·000
1·870
1·843
0·019
0·008
0·008
4·100
0·591
0·012
0·253
0·426
0·013
2·806
8·000
1·465
6·535
98·6
4·37
10·5
0·07
0·15
50·07
50·14
50·14
0·17
2·89
2·71
19·977
4·000
1·928
1·815
0·088
0·010
0·015
4·083
0·524
0·000
0·319
0·000
0·003
3·238
8·000
1·410
6·590
101·0
4·52
10·7
0·34
0·29
0·07
50·14
2·65
2·87
19·989
4·000
1·961
1·920
0·042
0·000
0·009
4·060
0·383
0·000
0·476
0·000
0·001
3·199
8·000
1·320
6·680
100·5
4·47
11·2
0·16
0·16
50·07
50·14
1·92
4·24
20·021
4·000
1·993
1·962
0·023
0·008
0·021
4·092
0·132
0·000
0·368
0·000
0·002
3·590
8·000
1·777
6·223
99·9
4·43
11·4
0·09
0·39
50·07
50·14
0·66
3·25
19·995
4·000
1·991
1·955
0·019
0·004
0·013
4·100
0·713
0·000
0·428
0·067
0·004
2·889
8·000
1·179
6·821
101·0
4·48
11·4
0·07
0·25
50·07
50·14
3·57
3·82
0·66
19·965
4·000
1·929
1·884
0·026
0·009
0·010
4·100
0·748
0·005
0·285
0·325
0·007
2·731
8·000
1·288
6·712
101·3
4·49
11·1
0·10
0·18
0·07
50·14
3·76
2·55
3·23
2·393
AlVI
XMg
OH
Ca
K
Na
Sum
Mg
Fe
Mn
0·648
16·000
0·012
0·022
0·000
11·884
6·084
3·302
0·101
0·004
2·214
AlIV
Ti
5·787
100·0
11·8
0·09
50·05
50·07
0·59
20·1
19·5
Si
Sum
H2O*
K2O
Na2O
CaO
MnO
MgO
FeO
50·05
0·485
16·000
0·010
0·011
0·002
11·707
4·249
4·520
0·218
0·002
2·718
2·170
5·831
98·9
0·661
16·000
0·046
0·012
0·008
11·912
6·318
3·239
0·075
0·002
2·278
2·218
5·782
100·0
11·8
0·05
50·05
0·21
0·44
20·9
19·1
50·05
18·8
28·6
M8359
(continued)
0·552
16·000
0·018
0·014
0·000
12·061
5·343
4·338
0·078
0·003
2·299
2·475
5·525
99·0
11·4
0·05
11·3
50·05
50·05
0·08
0·44
17·0
24·6
50·05
50·05
50·07
1·22
13·5
25·6
50·05
19·2
26·2
M8442
NUMBER 9
Ca
0·013
0·015
Fe3þ
3·117
AlVI
0·09
0·47
0·14
0·17
1·72
2·20
4·13
TiO2
19·6
27·6
M8441b
VOLUME 54
Ti
8·000
1·366
Sum
Al
IV
6·635
4·38
H2O*
Si
10·8
K2O
99·1
0·14
Na2O
Sum
0·06
0·56
CaO
BaO
0·15
1·80
4·45
0·00
0·07
19·3
28·5
0·17
3·12
3·88
3·50
50·05
Al2O3
MnO
SiO2
2·04
4·96
2·56
50·05
25·5
50·3
MgO
0·61
6·71
50·05
25·8
50·9
4·52
50·05
33·7
46·0
0·14
0·12
28·6
49·8
FeO*
0·11
50·05
29·8
49·7
Fe2O3*
0·06
26·4
47·6
0·13
26·9
50·7
27·4
47·3
0·12
26·6
48·7
TiO2
M83-59
26·1
M83-59
46·3
AS467
27·8
AS467
48·4
M84-42
SiO2
10Ca39
Al2O3
10Ca39
10Ca13
M84-41b
10Ca13
10Ca13
M84-41b
Chlorite
White mica
Table 2: Representative mineral analyses
JOURNAL OF PETROLOGY
SEPTEMBER 2013
50·08
23·4
50·08
23·2
1·89
MgO
CaO
H2O*
3·000
2·221
0·002
0·005
0·020
0·769
3·018
1·967
1·000
11·992
Si
Al
Ti
Mg
Mn
Fe
Sum
Ca
OH
O*
100·1
50·05
50·05
TiO2
Sum
0·14
0·33
1855
12·210
1·000
2·037
3·116
0·795
0·009
0·007
0·001
2·304
3·000
99·6
1·85
12·048
1·000
1·965
3·051
0·836
0·053
0·000
0·012
2·150
3·000
99·8
1·87
22·9
50·08
0·20
0·87
13·8
12·106
1·000
2·002
3·067
0·932
0·008
0·000
0·004
2·124
3·000
100·7
1·87
23·3
50·08
0·06
0·13
15·4
0·06
Orthoclase
Anorthite
Albite
Sum
K
Na
Ba
0·199
4·235
95·567
1·015
0·002
0·970
0·000
0·043
0·003
Fe3þ
Ca
1·046
2·949
Al
Si
0·341
0·658
99·001
0·984
0·003
0·974
0·000
0·007
0·000
1·009
2·996
100·0
50·05
100·3
K2O
Sum
11·5
50·08
0·14
11·4
50·08
0·92
50·14
50·13
Na2O
BaO
CaO
FeO
19·6
68·7
10Ca11
20·3
67·5
Mn2O3
13·0
22·5
Al2O3
12·9
SiO2
Fe2O3
22·7
37·4
23·8
37·4
37·0
37·9
SiO2
Al2O3
24·1
M8441b
M8359
M8442
10Ca13
M8441b
Feldspar
Epidote
Table 2: Continued
0·803
34·157
64·975
1·024
0·008
0·666
0·001
0·350
0·012
1·310
2·666
99·3
0·14
7·63
50·08
7·26
0·31
24·7
59·3
M8442
98·108
0·115
1·081
0·967
0·949
0·010
0·007
0·001
0·023
1·073
2·931
99·7
16·0
0·12
Sum
OH
O*
Ca
Mn
Fe
Al
Ti
Si
Sum
H2O*
CaO
MnO
Fe2O3
50·07
0·37
Al2O3
TiO2
SiO2
0·58
19·6
63·1
10Ca33
4·974
0·102
4·872
0·906
0·001
0·076
0·128
0·819
0·988
99·5
0·47
26·0
50·14
3·11
3·33
33·5
30·4
M8359
Titanite
4·951
0·113
4·838
1·000
0·003
0·024
0·202
0·777
0·992
99·3
0·52
28·8
50·14
1·00
5·28
31·9
30·6
M8441b
0·894
0·757
0·019
2·000
6·001
Mg
Ca
Na
Sum
O*
0·007
0·194
Fe2þ
Mn
0·019
0·089
Fe3þ
AlVI
Ti
2·000
0·021
Sum
1·872
0·129
AlIV
100·1
0·26
19·1
16·2
0·23
6·28
3·19
0·67
3·43
50·6
Si
Sum
Na2O
CaO
MgO
MnO
FeO*
Fe2O3*
TiO2
Al2O3
SiO2
M8339
(continued)
6·000
2·000
0·010
0·793
0·920
0·007
0·208
0·031
0·009
0·023
2·000
0·061
1·939
99·0
0·13
19·9
16·6
0·22
6·69
1·11
0·30
1·91
52·1
10Ca23
Clinopyroxene
WILLNER et al.
VERY LOW- TO LOW-GRADE METAMORPHIC PROCESSES
24·3
17·4
8·25
0·41
Al2O3
Fe2O3*
FeO*
MnO
0·000
0·227
2·155
1·133
1856
0·012
1·702
0·057
6·000
0·010
0·008
0·033
0·678
0·728
21·052
6·000
0·655
0·600
Ti
Mg
Mn
Sum
Ca
Ba
Na
K
Sum
O*
OH
XFe3þ
XMg
XFe
XMg
OH
Sum
K
Na
Ca
Sum
Mg
0·942
0·630
2·000
2·058
0·033
0·151
1·874
13·000
2·976
1·646
Fe2þ
0·051
0·054
Fe3þ
Mn
0·005
0·267
0·713
0·292
XMg
21·771
O*
7·000
3·980
0·022
3·954
1·000
0·278
0·032
0·690
5·279
1·541
0·007
3·731
XFe3þ
OH
Sum
Na
Ca
Sum
Fe2þ
Mn
Mg
Sum
Fe3þ
Ti
Al
6·191
99·3
6·13
0·07
21·6
2·70
0·22
1·94
12·0
0·05
0·126
OH
Ca
2·000
1·974
1·004
Fe3þ
Sum
4·000
0·878
Sum
AlVI
3·010
0·990
Si
98·8
4·29
26·4
2·39
22·7
43·1
AlIV
Sum
H2O*
CaO
Fe2O3
Al2O3
SiO2
*Calculated.
White mica: the proportion of cations is based on 42 negative charges neglecting the interlayer cations; the sum of octahedrally coordinated cations is set at 4·1 to allow for an
estimation of Fe3þ. Chlorite: cations based on 56 negative charges; H2O calculated on the basis of OH ¼ 16. Epidote: proportions of cations are based on normalization of Si to three
cations. Feldspar: normalization on the basis of 16 negative charges. Titanite: sum of cations ¼ 3; OH ¼ (Al þ Fe) 0·5; O ¼ [( positive charges) – OH] 0·5. Clinopyroxene: normalization to four cations to calculate Fe3þ. Stilpnomelane: the proportion of cations is based on 47·375 negative charges, neglecting K þ Na; the proportion of Fe3þ is estimated assuming
15 cations; H2O calculated on the basis of OH ¼ 6. Amphibole: the proportion of cations is based on the sum of cations ¼ 13 except for Ca, Na and K for estimation of Fe3þ and on 46
negative charges. Pumpellyite: cations are based on 49 negative charges; H2O is calculated on the basis of OH ¼ 4; estimation of Fe2þ by assuming half of X-position filled with divalent
cations. Prehnite: the structural formula is based on 20 negative charges and OH ¼ 4.
6·000
21·020
0·657
0·034
AlVI
Ti
8·000
0·264
AlIV
Sum
7·736
Si
Si
Sum
H2O*
Na2O
CaO
MgO
MnO
FeO*
Fe2O3*
TiO2
18·5
36·2
10Ca13
Prehnite
NUMBER 9
0·614
0·007
0·002
5·290
0·622
0·782
2·05
98·8
H2O*
Sum
0·41
0·53
0·17
12·0
13·7
13·5
0·49
0·05
MnO
Na2O
K 2O
CaO
MgO
FeO*
Fe2O3*
TiO2
Al2O3
SiO2
10Ca23
Pumpellyite
VOLUME 54
0·002
2·664
1·220
Fe2þ
9·000
0·430
Fe3þ
AlIV
8·570
98·9
9·000
3·763
Si
0·941
5·237
Sum
0·11
5·68
Sum
98·3
H2O*
2·00
0·17
AlVI
0·12
5·48
BaO
0·10
3·23
Na2O
K2O
3·31
50·07
6·95
50·07
MgO
CaO
4·63
20·1
8·83
3·08
52·9
0·10
TiO2
Al2O3
54·1
31·9
SiO2
SiO2
M8442
M8441b
10Ca44
0·02
Amphibole
Stilpnomelane
Table 2: Continued
JOURNAL OF PETROLOGY
SEPTEMBER 2013
1857
0·90
0·85
15·16
6·57
4·51
Al2O3
FeO þ MnO
MgO
2·80
100·0
0·04
8·96
100·0
H2O
CO2
Sum
100·0
6·96
0·04
1·56
6·96
0·04
100·0
*Corrected for CaO content in apatite.
100·0
8·96
0·04
0·58
7·95
5·96
0·04
2·75
2·32
2·03
8·17
3·52
K2O
4·56
3·68
O2
8·24
2·56
5·84
5·89
9·33
17·72
0·90
44·3
2·40
9·02
9·94
14·95
0·68
46·7
100·5
Na2O
4·19
8·02
18·38
0·71
51·9
4·42
101·1
CaO*
5·96
9·73
15·37
47·7
53·6
SiO2
TiO2
Simplified compositions
3·50
99·8
4·56
3·22
101·4
2·97
0·26
2·76
100·6
0·13
0·27
2·33
8·31
5·91
0·26
10·1
17·8
0·91
44·5
Sum
0·13
1·60
3·60
8·53
9·22
0·30
10·9
15·3
0·70
47·7
H2O
0·19
0·61
3·77
4·83
4·29
0·16
8·93
18·8
0·73
53·1
2·81
0·17
2·71
8·97
6·31
0·24
11·2
16·3
0·95
50·5
10CA40
CO2
P2O5
2·03
5·30
MgO
K2O
0·14
MnO
7·03
7·88
Fe2O3
2·42
15·3
Al2O3
Na2O
0·88
TiO2
CaO
56·3
SiO2
Whole-rock analyses
10Ca23
100·0
5·96
0·04
2·46
8·92
4·29
9·29
16·51
0·54
52·0
100·6
4·60
0·03
0·05
2·48
9·05
4·33
0·16
10·2
16·7
0·54
52·5
FA8785
100·0
9·96
0·04
1·89
2·58
8·04
3·68
7·12
19·33
0·77
46·6
100·1
3·00
0·10
1·80
2·50
9·80
5·20
0·20
8·28
18·7
0·80
49·7
AS467
100·0
7·96
0·04
0·44
2·19
10·12
6·83
9·25
16·52
0·81
45·8
99·2
1·60
0·15
0·46
2·29
10·8
7·15
0·31
10·4
17·3
0·85
48·0
FC881614
100·0
6·96
0·04
1·09
4·91
1·93
1·73
3·83
14·33
0·56
64·6
100·5
1·91
0·11
1·15
5·16
2·18
1·82
0·17
4·28
15·0
0·59
67·8
100·0
6·96
0·04
1·10
5·36
0·96
0·39
1·99
12·78
0·36
70·1
98·8
0·15
0·05
1·17
5·66
1·09
0·41
0·17
2·15
13·5
0·38
74·0
M84-41B
M83-57
10CA13
M83-59
M84-42
Felsic rocks
Mafic rocks
100·0
6·96
0·04
1·79
4·04
1·53
0·72
2·65
13·03
0·39
68·9
98·7
0·50
0·06
1·88
4·24
1·68
0·76
0·09
2·99
13·7
0·41
72·2
M84-41C
100·0
6·96
0·04
0·90
5·10
1·41
1·01
2·93
12·87
0·38
68·4
98·6
0·25
0·08
0·95
5·35
1·58
1·06
0·11
3·29
13·5
0·40
71·8
10CA11
Table 3: Whole-rock analyses and simplified compositions for the calculation of P^T pseudosections (Figs 5^10)
100·0
6·96
0·04
0·40
4·49
2·93
0·33
2·29
11·57
0·30
70·7
99·4
0·72
0·04
0·43
4·74
3·14
0·35
0·13
2·54
12·2
0·31
74·6
10CA33
100·0
6·96
0·04
1·35
3·79
2·90
1·71
5·23
14·28
0·60
63·1
100·6
2·19
0·12
1·41
3·96
3·19
1·79
0·18
5·88
14·9
0·63
66·0
10CA34
100·0
1·21
5·96
0·04
3·05
3·07
2·04
1·25
3·42
14·54
0·55
64·9
99·8
1·88
1·25
0·11
3·16
3·17
2·26
1·30
0·07
3·92
15·0
0·57
67·0
10CA37
100·0
0·57
5·96
0·04
3·27
3·82
1·37
0·70
2·07
13·99
0·29
67·9
98·2
0·59
0·06
3·41
3·98
1·50
0·73
0·06
2·33
14·6
0·30
70·7
10CA39
100·0
6·96
0·04
0·72
4·85
1·08
0·19
1·10
10·87
0·14
74·1
98·6
0·01
0·77
5·14
1·16
0·20
0·06
1·23
11·5
0·14
78·4
10CA44
WILLNER et al.
VERY LOW- TO LOW-GRADE METAMORPHIC PROCESSES
JOURNAL OF PETROLOGY
VOLUME 54
Powell (1996) supplemented by an acmite component,
which is commonly enhanced at very low-grade conditions. Albite, K-feldspar, prehnite, quartz, titanite, H2O,
ilmenite and magnetite were considered as pure phases.
Pseudosections were calculated for eight mafic and nine
felsic samples (Table 1). The major element analyses were
simplified to the above 10-component system, correcting
Ca for apatite and adding MnO to FeO. Water contents
were increased to excess water conditions that are considered to have prevailed during peak P^T conditions.
Finally analyses were normalized to 100% (Table 3). For
samples where calcite was present in considerable amounts,
CO2 was added to the 10-component system. However,
Massonne (2010) showed that even fluids released from
marly limestone at low temperature contain very little
CO2. Hence CO2 can largely be neglected as a fluid component at temperatures54008C. Oxygen contents were arbitrarily chosen (0·04 wt %) to account for trivalent iron
in epidote and pumpellyite. Magmatic relict phases occur
only in a few samples, and do not exceed 5%. Because
of their low abundances, the whole-rock chemical data
were not adjusted to accommodate their presence.
R E S U LT S
Petrographic characteristics and mineral
chemistry
The protoliths of the selected samples include both lava
flows and various pyroclastic rocks. Metamorphic mineral
assemblages and fabrics are similar in samples from all of
the volcanic^sedimentary belts of the Mira terrane, and
hence we do not distinguish samples with respect to location in the descriptive sections.
Felsic rocks
The felsic rocks are sheared and foliated but retain recognizable relict fabrics such as porphyritic and/or pyroclastic
textures, compositional layering, and rare relict flow-banding. Locally, deformation has entirely obliterated the
primary fabric. The common assemblage is phengite^epidote^chlorite^albite^K-feldspar^quartz^titanite stilpnomelane calcite amphibole ilmenite. Plagioclase and
K-feldspar phenocrysts and clasts (0·1^2 mm size) are conspicuous relicts and commonly rotated between phyllosilicate bands (s-clasts) showing pressure shadows, some
pressure solution and abundant brittle deformation.
Plagioclase composition is always albite, which pseudomorphed the magmatic relicts or grew within the matrix.
Quartz occurs in places as primary phenocrysts with corrosion embayments (0·5^1mm), but mostly in the matrix
or as fissure fillings, where it is invariably recrystallized
to a polygonal fabric, especially in pressure shadows
(0·01^0·05 mm; Fig. 2a and b).
White mica is entirely of metamorphic origin, finegrained (0·005^0·03 mm) and commonly oriented parallel
NUMBER 9
SEPTEMBER 2013
to the predominant foliation (Fig. 2a and b). It typically
forms unoriented irregular clusters about 0·1^0·2 mm in
diameter mostly across feldspar grains, but also in the
matrix. White mica composition is mostly phengite [i.e.
3·2 Si a.p.f.u. following the proposal of Massonne &
Schreyer (1986)]. Phengite is an intermediate member of a
solid-solution series between the end-members muscovite
and (Mg, Fe)-aluminoceladonite (e.g. Rieder et al., 1998).
In our samples, it has a wide variation in Si content
(3·13^3·43 Si a.p.f.u.), but similar ranges in all of the
samples. Some compositions show wide deviation from
the ideal Tschermak’s substitution 2Al ¼ Si þ (Mg þ Fe2þ)
owing to highly variable Fe3þ-substitution (calculated
Fe3þ 0·0^0·7 a.p.f.u.; Fig. 3a) and to a minor extent to
the di/trioctahedral substitution 2Al þ œ ¼ 3(Mg, Fe2þ).
Consistent with very low metamorphic grade, Na contents are low (0·005^0·045 a.p.f.u.), as are the contents
of Ba (0·016 a.p.f.u.), Ti (0·018 a.p.f.u.) and Mn
(0·014 a.p.f.u.). Mica grains with unusually high Na contents 40·1a.p.f.u. and corresponding elevated Si contents
are interpreted as analyses that included both white mica
and adjacent albite as a result of the very fine grain size
and hence are discarded as artefacts. XMg [¼ Mg/
(Mg þ Fe2þ)] varies considerably between and within
samples (0·4^0·9); Si and XMg appear to be slightly anticorrelated within each sample (Fig. 3b). Detailed X-ray
element distribution maps (Fig. 4) of typical small white
mica grains show strongly irregular zoning with both
increasing and decreasing Si contents from inner portions
to the rim. Also, neighbouring grains can have very different compositions. The absence of consistent compositional
zoning patterns in white mica is typical for the very low
metamorphic grade samples and has also been observed
in all other coexisting solid-solution phases.
Chlorite is a minor phase in the felsic rocks. It shows
little compositional variation within samples, but more
variation between samples (Si 5·5^6·9 a.p.f.u.; XMg
0·30^0·74; Ti 50·02, locally 0·12^0·19; Mn 0·08^0·35). It is
mostly intergrown with oriented white mica of similar
grain size, but also forms clusters together with titanite
(Fig. 2a and b).
Epidote grains form conspicuous clusters 0·2^0·5 mm in
diameter with single grain sizes of 0·01^0·05 mm. The clusters are slightly elongate parallel to the foliation. Contents
of Fe3þ in epidote vary considerably between 0·5 and
1·0 a.p.f.u. Mn (50·12 a.p.f.u.) and Ti (50·07 a.p.f.u.) are
minor constituents.
Stilpnomelane commonly occurs as mostly unoriented
grains (0·01^0·05 mm in length) in small clusters
0·2^0·5 mm in diameter (Fig. 2a and b). In some samples
stilpnomelane occurs in the dominant foliation defined by
oriented phengite and chlorite. Characteristically, stilpnomelane compositions indicate a strong and highly variable
reduction of the interlayer cations and a wide variation in
1858
WILLNER et al.
VERY LOW- TO LOW-GRADE METAMORPHIC PROCESSES
Fig. 2. Photomicrographs of (a, b) felsic sample M84-41B with clusters of chlorite in albite and interstitial clusters of stilpnomelane (St) and
phengite and (c, d) mafic sample 10Ca13 showing clusters of prehnite (Pr), phengite (Ph) and epidote (Ep) within albite (Ab) and interstitial
chlorite (Ch) clusters. Scale bars represent 0·5 mm. þP, crossed polars; llP, parallel polars.
3þ
XFe
, here defined as Fe3þ/(Fe3þ þ Fe2þ), between 0·0
and 0·6 and XMg [¼ Mg/(Mg þ Fe2þ)] between 0·2 and
0·6. Most of this variation is assumed to be due to oxidation processes related to post-metamorphic alteration.
Other components include Mn (0·05^1·0 a.p.f.u.), Ti
(50·02 a.p.f.u.) and Ba (50·02 a.p.f.u.).
Titanite is a common accessory phase that occurs in
clusters 0·05^0·2 mm in diameter with single grain lengths
of 0·01^0·03 mm. Grains are elongate in laminae parallel
to the foliation. Owing to a minor vuagnatite substitution
[Ti4þ þ O2^ ¼ (Al,Fe)3þ þ (OH,F)] contents of Al
(0·12^0·27 a.p.f.u.) and Fe3þ (0·02^0·05 a.p.f.u.) are notable.
Ilmenite, presumably of magmatic origin, is present in
some samples as a second Ti phase.
Mafic rocks
In contrast to the felsic rocks, many mafic rocks are undeformed and preserve relict porphyritic, subophitic and
amygdaloidal textures. Mafic rocks with relict pyroclastic
textures show weak foliation defined by chlorite and
white mica, and alignment of small feldspar clasts and
lithoclasts. As described by McMullin et al. (2010) two
principal assemblages occur: (1) epidote^chlorite^albite^
quartz^titanite phengite pumpellyite prehnite calcite K-feldspar indicative of the prehnite^actinolite
facies; (2) actinolite^epidote^chlorite^albite^quartz^titanite phengite pumpellyite calcite K-feldspar indicative of the greenschist facies.
Heterogeneous overprinting is partly indicated by local
relicts of magmatic clinopyroxene and plagioclase.
Although plagioclase composition is generally albite,
some rare local patches of more calcic plagioclase (An14^
34) occur in phenocrysts and are probably relict. Augite
makes up to 5 vol. % in some of the samples. Its composition is variable both within and between samples (diopside0·36^0·56, hedenbergite0·06^0·22, orthopyroxene0·19^0·34,
acmite0·01^0·11, Tschermak component0·03^0·14) and also
includes minor Ti (0·004^0·025 a.p.f.u.) and Mn (0·004^
0·014 a.p.f.u.). Rare relict amygdules contain quartz, epidote and/or chlorite. Pyroclastic samples contain some
lithic clasts, which are schistose or mylonitic, indicating
pre-existing deformation.
1859
JOURNAL OF PETROLOGY
VOLUME 54
Fig. 3. Variation of Si content (a.p.f.u.) with (a) Al and (b) XMg [¼
Mg/(Mg þ Fe2þ)] in white mica of mafic and felsic samples. The line
of the ideal Tschermak’s substitution is indicated in (a).
White mica mainly forms 0·1^0·5 mm clusters of 0·01^
0·05 mm grains within plagioclase (Fig. 2c and d), but is
also disseminated in the matrix. Composition is mostly
phengite (i.e. Si43·2 a.p.f.u.) and its range is very similar
to that observed in felsic rocks (Si 3·11^3·41a.p.f.u.; calculated Fe3þ 0·00^0·76 a.p.f.u.; Na 0·01^0·05 a.p.f.u.; Ba
0·01a.p.f.u. with local 0·02 a.p.f.u.; Mn 0·015 a.p.f.u.;
Ti 0·00^0·01a.p.f.u.; XMg 0·17^0·93; Fig. 3a). Slight negative correlation is evident between Si and XMg [¼ Mg/
(Mg þ Fe2þ)] in each sample (Fig. 3b).
Ubiquitous chlorite with a grain size similar to that of
white mica is mostly disseminated in the matrix, and is
oriented in foliated samples. It also forms small clusters
(Fig. 2c and d) and fills fissures. It has similar contents of
Si (5·5^6·6 a.p.f.u.), Ti (50·06, locally 0·14^0·24 a.p.f.u.)
and Mn (0·04^0·21a.p.f.u.) to chlorite in felsic rocks, but
generally has higher XMg (0·44^0·71). Like white mica, it
displays weak negative correlation between Si and XMg.
Albite and K-feldspar were observed as part of the metamorphic assemblage in several mafic samples, and quartz
occurs in all samples as small pods of recrystallized
NUMBER 9
SEPTEMBER 2013
crystals in the matrix or in fissures. Amphibole is a minor
component, occurring as very small disseminated crystals
(0·01^0·02 mm). It is actinolite (XMg 0·6) and contains
small amounts of Na (0·09^0·12 a.p.f.u.).
Epidote grains (0·05^0·1mm in size) commonly occur in
conspicuous clusters (0·5^2 mm in size) or are disseminated in the matrix and in relict plagioclase phenocrysts
(Fig. 2c and d). Epidote is similar in composition to that
in felsic rocks (Fe3þ 0·52^1·0 a.p.f.u., Mn 50·05 a.p.f.u., Ti
50·03 a.p.f.u., and Mg 50·12 a.p.f.u.).
Pumpellyite forms clusters similar to those of epidote,
but is less common. Compositional variation of pumpel3þ
[¼ Fe3þ/(Fe3þ þ Al)]
lyite in sample 10Ca23 is XFe
0·27^0·34 and XMg [¼ Mg/(Mg þ Fe2þ)] 0·52^0·71 with
rare exceptional XMg down to 0·23 or up to 0·86, presumably as a result of locally variable oxidation conditions.
According to McMullin et al. (2010) XMg varies in the
range 0·30^0·77 between samples.
Prehnite is present mostly in clusters within former
plagioclase (Fig. 2c and d) or fills fissures. It contains notable amounts of Fe3þ (0·10^0·13 a.p.f.u.), as also noted by
McMullin et al. (2010). Prehnite was observed in coexistence with pumpellyite in only three samples [FC881614,
FA8785 and AS467, all taken from the sample set of
McMullin et al. (2010)].
Titanite commonly forms diffuse clusters elongate parallel to the foliation. Low contents of Al (0·06^0·21a.p.f.u.)
and Fe3þ (0·04^0·08 a.p.f.u.) indicate a minor vuagnatite
substitution. Also present are ilmenite and magnetite,
which may be relict magmatic phases.
Phase relationships deduced by P^T
pseudosections
P^T pseudosections show the extent of P^T fields of specific index minerals varying with whole-rock composition.
Pseudosections calculated for four mafic and felsic rock
compositions from the Mira terrane were reduced to the
P^T fields of critical very low- and low-grade phases
(Fig. 5). The P^T conditions are between the P^T fields of
laumontite and lawsonite on the low-temperature side and
biotite at the high-temperature side (Fig. 5a and b), because none of these minerals was observed. This field
approximately corresponds to the temperature range of
250^3008C. Notably, considerable similarity exists between
the P^T fields of the three index minerals in mafic and
felsic compositions. Pumpellyite has a smaller P^T field
in the felsic samples of this study (Fig. 5c), whereas
both prehnite and pumpellyite occur in the mafic samples
(Fig. 5d). The P^T fields of both minerals have very restricted overlap with a P^T maximum at 3·5 kbar and
2908C. Although we calculated the pseudosections
using only pure prehnite, the slight Fe3þ-substitution
observed in prehnite would have only a minor effect
(change in temperature up to 58C calculated for multivariant reactions). The observation of a restricted
1860
WILLNER et al.
VERY LOW- TO LOW-GRADE METAMORPHIC PROCESSES
Fig. 4. Maps of relative X-ray intensity distribution of Al, Si, Fe, Mg in white mica from very low-grade felsic samples M84-41 and 10Ca39.
Quantitative spot compositions as a.p.f.u. of Al, Si, Fe and Mg are shown. It should be noted that some increase or decrease of elements is due
to grain boundary effects.
prehnite^pumpellyite field has also been made elsewhere
by Frey et al. (1991). Owing to its limited stability field the
mineral pair should not be used to indicate its own facies,
but should be regarded as a subfacies within the prehnite^
actinolite facies, as also concluded by McMullin et al.
(2010).
Stilpnomelane and K-feldspar are also index minerals in
the studied mafic and felsic samples (Fig. 5e and f), but
with different extension of the P^T fields: stilpnomelane
appears in mafic and felsic rocks up to 250^3008C (at 2^6
kbar), whereas K-feldspar occurs only below 3·2 kbar in
mafic rocks, but also at higher pressure in very low-grade
felsic rocks. The appearance of epidote above 250^3008C
and the occasional replacement of white mica by K-feldspar at low pressure and temperature appears common in
both mafic and felsic rocks from the Mira terrane. In
general, the deduced P^T fields of the very low- and
low-grade metamorphic phases concur with earlier petrogenetic grids, which were constructed using approximated
mineral compositions (e.g. Brown, 1975; Liou et al., 1985,
1987; Frey et al., 1991; Banno, 1998).
Metamorphic conditions derived with
P^T pseudosections
Peak P^T conditions are estimated from P^T pseudosections by the position of isopleths representing the maximum observed contents of certain elements in one or two
phases in the field of the observed peak metamorphic
assemblage. For very low- and low-grade conditions, Si
contents of white mica, a common phase in the Mira terrane samples, are most variable. Hence in Fig. 6 isopleths
of Si contents are inserted into P^T pseudosections
of felsic sample 10Ca11 (Fig. 6a) and mafic sample
11Ca13 (Fig. 6b). Although Si contents often increase
1861
JOURNAL OF PETROLOGY
VOLUME 54
NUMBER 9
SEPTEMBER 2013
Fig. 5. P^T fields of minerals extracted from calculated pseudosections for the whole-rock compositions of mafic samples (a, c, e, g) and felsic
samples (b, d, f, h).
1862
WILLNER et al.
VERY LOW- TO LOW-GRADE METAMORPHIC PROCESSES
Fig. 6. P^T pseudosections calculated for the whole-rock compositions of (a) felsic sample 10Ca11 and (b) mafic sample 10Ca13 with isopleths
of Si a.p.f.u. in white mica (dashed lines). Shading indicates the degree of variance as shown. The highest and lowest Si contents are indicated
as bold dashed lines within the field of the observed assemblage (hatched). White dotted lines in (a) indicate XMg [¼ Mg/(Mg þ Fe)] for chlorite. Abbreviations of minerals and end-member components in this paper are: Ab, albite; Am, amphibole; Bt, biotite; CA, calcic amphibole;
Cc, calcite; Ch, chlorite; Cp, clinopyroxene; Ep, epidote; Gt, garnet; Kf, potassic feldspar; Lm, laumontite; Lw, lawsonite; Mt, magnetite; NA,
sodic amphibole; Pa, paragonite; Ph, phengite; Pl, plagioclase; Pr, prehnite; Pu, pumpellyite; Qz, quartz; St, stilpnomelane; Tt, titanite; V,
H2O as hydrous fluid; Wm, potassic white mica; Wk, wairakite.
1863
JOURNAL OF PETROLOGY
VOLUME 54
with pressure, the form of the Si isopleths is irregular in
both felsic and mafic samples and may be used as a
geothermobarometer in combination with the P^T field of
the observed metamorphic assemblage.
Thus, in felsic sample 10Ca11 (Fig. 6a) peak P^T conditions are reached at 265^2858C and 3·4^3·7 kbar with
growth of phengite with 3·36 Si a.p.f.u. in an assemblage
with chlorite, epidote, potassic feldspar, albite, titanite
and quartz. The isopleth for the lowest Si content measured in this sample (3·28 a.p.f.u.) plots in the same P^T
field and indicates retrograde conditions. Also, XMg [¼
Mg/(Mg þ Fe)] isopleths for chlorite inserted in Fig. 6a
show little variation around 0·40, similar to the observed
values and to other pseudosections.
Similar peak P^T conditions at 280^2908C and 3·2^3·4
kbar are indicated for mafic sample 10Ca13 (Fig. 6b) as
shown by the isopleth of maximum Si content in phengite
(3·33 a.p.f.u.) in an assemblage with epidote, prehnite, potassic feldspar, chlorite, albite, titanite and quartz. The isopleth for the lowest Si content of 3·20 a.p.f.u. also plots in
the same P^T field, representing retrograde white mica.
However, two limitations are indicated for this example.
First, ilmenite and magnetite were observed in the
NUMBER 9
SEPTEMBER 2013
sample, but do not occur in the calculated P^T pseudosections, as is also the case for pseudosections calculated for
the other three mafic samples. This indicates that these
accessory phases most probably represent primary
phases. On the other hand, a small amount of K-feldspar
(5 vol. %) was calculated for the likely field of the peak
assemblage, but not observed in the sample. However,
K-feldspar was observed in other mafic samples and small
amounts of K-feldspar in a fine-grained rock may have
been overlooked, or could occur in parts of the sample not
included in the thin section.
The P^T pseudosection for mafic sample AS467 that contains pumpellyite and prehnite in an assemblage with phengite, epidote, chlorite, albite and titanite is shown in Fig. 7.
Again, the very limited field for coexistence of pumpellyite
and prehnite is apparent. The peak metamorphic assemblage occurs at the upper P^T stability range for coexisting
prehnite and pumpellyite at 270^2958C and 3·2^3·8 kbar.
The isopleth for coexisting white mica with the highest Si
content (3·34 a.p.f.u.) almost coincides with this small field.
The field for coexisting prehnite and pumpellyite thus
limits the maximum P^Tconditions even in samples where
no white mica is present, such as sample FA8785 (Table 1).
Fig. 7. P^T pseudosection calculated for the whole-rock composition of mafic sample AS467 containing pumpellyite and prehnite.
Abbreviations and shading as in Fig. 6. The narrow P^T field of assemblages with coexisting prehnite and pumpellyite should be noted (6, 7,
21, 26, 19). Dashed lines indicate the Si contents in white mica. The maximum content of 3·34 a.p.f.u. is almost identical to the field of the peak
metamorphic assemblage (19; hatched). White dotted lines are XMg ¼ Mg/(Mg þ Fe2þ) in pumpellyite. Mineral abbreviations are as in Fig. 6.
1864
WILLNER et al.
VERY LOW- TO LOW-GRADE METAMORPHIC PROCESSES
In sample AS467, with white mica Si contents lower than
3·34 a.p.f.u. but mostly above 3·21a.p.f.u. (except for one
analysis with 3·11a.p.f.u.) probably grew mainly during
retrograde conditions. However, no retrograde K-feldspar
was observed in the sample, which would be expected in accessory amounts according to the calculation and hence
might have been missed in petrographic work. Calculated
XMg isopleths for pumpellyite range from 0·58 to 0·64
(Fig. 7), compatible with measured XMg values in this
sample (0·51 0·10; McMullin et al., 2010).
Mafic sample 10Ca23 lacks white mica but has the
critical assemblage pumpellyite þ epidote (Fig. 8). This assemblage is restricted to a narrow temperature range of
280^3108C at pressures exceeding 3·0^3·3 kbar. The calculated range of XMg ¼ 0·68^0·74 for pumpellyite (Fig. 8)
corresponds well to the observed XMg of 0·52^0·71 (Table
2). Both actinolite and clinopyroxene are part of the calculated assemblage in amounts less than 10 vol. %. No actinolite was observed in the sample, and may have been
overlooked in the extremely fine-grained matrix. Some
relict clinopyroxene is present; it was previously argued by
McMullin et al. (2010) that this clinopyroxene could also
be part of the metamorphic assemblage, an inference
supported by its presence in the calculated assemblage.
The calculated clinopyroxene composition for sample
10Ca23 at 2808C and 3·0 kbar is diopside0·63,
hedenbergite0·30, orthopyroxene0·01, and acmite0·06, relatively similar to the observed compositional range (see
above).
Two special cases where prograde metamorphic phases
are still preserved are shown in Fig. 9. Sample M84-42 contains a typical greenschist assemblage (Fig. 9a). Peak P^T
conditions are 290^3208C and 3·3^3·9 kbar based on the
highest Si contents of 3·30 a.p.f.u. in the P^T field of the
observed assemblage of phengite, actinolite, epidote, chlorite, albite, quartz and titanite. This sample contains ilmenite rather than titanite and, as in sample 10Ca13, the
ilmenite is interpreted as a relict magmatic phase.
The range of maximum temperature for this sample (290^
3208C) is only slightly higher than the range of 270^2958C
for the typical very low-grade assemblage with pumpellyite
and prehnite in sample AS467 (Fig.7). This slight difference
in peak temperature may be real, but is also within limits
of error. Na contents in actinolite increase steadily towards
the transition to Na-amphibole at 5^7 kbar and 250^
3008C (Fig. 9a). It is important that in this rock accessory
Fig. 8. P^T pseudosection calculated for the whole-rock composition of mafic sample 10Ca23 containing pumpellyite and epidote. The narrow
P^T fields of assemblages with coexisting prehnite and pumpellyite should be noted. Abbreviations and shading as in Fig. 6. Dashed lines are
XMg ¼ Mg/(Mg þ Fe2þ) in pumpellyite. Mineral abbreviations are as in Fig. 6.
1865
JOURNAL OF PETROLOGY
VOLUME 54
NUMBER 9
SEPTEMBER 2013
Fig. 9. P^T pseudosections calculated for the whole-rock composition of (a) mafic sample M84-42 with a typical greenschist-facies assemblage
(hatched) and (b) felsic sample M84-41B with prograde growth of stilpnomelane. Abbreviations and shading as in Fig. 6. Isopleths of Si contents
in white micas are the black dashed lines. White dotted lines indicate Na contents in actinolite in (a) and XMg [¼ Mg/(Mg þ Fe)] in stilpnomelane in (b). Mineral abbreviations are as in Fig. 6.
stilpnomelane, pumpellyite and K-feldspar were observed.
These minerals grew as a prograde assemblage in field 17
of Fig. 9a of the calculated pseudosection at 2·3^3·3 kbar
and 230^2508C. This finding implies a late prograde P^T
path with little pressure increase. The P^T field for white
mica calculated for sample M84-42 is restricted to pressures
43 kbar and Si contents exceeding 3·28 a.p.f.u. However,
white mica with Si contents as low as 3·2 a.p.f.u. is present,
and probably formed under retrograde conditions during
restricted water access (e.g. Cruciani et al., 2011) and a reacting chemical composition that deviated slightly from the
whole-rock composition.
1866
WILLNER et al.
VERY LOW- TO LOW-GRADE METAMORPHIC PROCESSES
The highest Si content (3·31a.p.f.u.) in phengite in felsic
sample M84-41B occurs in two fields with different
observed assemblages, occurring with stilpnomelane at
higher pressure and with potassic feldspar at higher temperature (Fig. 9b). Traces of sodic amphibole are also calculated in both fields, although this mineral was not observed
in thin section. In this case the possible range for peak pressure appears large. However, stilpnomelane is preserved as
a prograde phase formed around 2508C and the prograde
P^T path should pass near the isopleth with Si 3·31a.p.f.u.
up to maximum possible temperatures of 3108C. The isopleth with the lowest Si contents of 3·24 a.p.f.u. also cuts
the second field (at higher temperature), indicating retrograde growth of white mica below 3 kbar and 3008C. The
range of calculated values of 0·20^0·28 for XMg for stilpnomelane, shown as isopleths in Fig. 9b, fits well with observed
values of 0·21^0·28 (Table 2).
In a way similar to that illustrated in Figs 6^9, peak P^T
conditions were estimated for seven additional mafic and
nine additional felsic samples using maximum Si contents
in phengite (Table 1). All results strongly overlap and show
no consistent differences between samples. Remarkably
uniform peak metamorphic conditions within the range of
280 308C and 3·5 0·4 kbar are indicated for the entire
Mira terrane. Owing to the considerable range in wholerock compositions, assemblages of the prehnite^actinolite
and greenschist facies can overlap in this P^T range, but
slight local differences in temperature and/or pressure are
also possible. Such differences lie within the magnitude of
general errors of the geothermobarometric methods (0·5
kbar, 208C at 1s level). The thermometric approach
of McMullin et al. (2010), calculating multivariant reactions with mineral pairs, yielded an approximately similar
range of temperatures (244^3078C). No difference in
microfabric (grain size, lack of zonation, clustering,
degree of recrystallization) was identified between samples
in which subgreenschist- and greenschist-facies assemblages occur.
Isopleths of water content bound to solids also can be
extracted from the calculated P^T pseudosections. Calculated isopleths of water content (in wt %) of representative
mafic and felsic samples show that water was not released
continuously during metamorphism (Fig. 10). Most water
was released during breakdown of OH-bearing minerals
at very low-grade conditions. In the temperature range of
250^3008C, in which most metamorphic reactions took
place in these samples (Fig. 4), 3 wt % water in mafic
sample M84-42 and 1·5 wt % water in felsic sample
10Ca34 (Fig. 5) was released by prograde mineral reactions. This is equivalent to 40% and 25%, respectively,
of water bound to minerals at 2008C. In both cases water
release is most pronounced at around 2508C. This observation has important implications for metamorphic processes
at very low grade as discussed in the next section.
Fig. 10. Isopleths of weight per cent H2OSolid extracted from P^T
pseudosections of (a) mafic sample M84-42 and (b) felsic sample
10Ca34.
DISCUSSION
Kinetic and fluid control of metamorphic
reactions at very low-grade conditions
Before applying any geothermobarometric method the heterogeneity of the composition of the minerals at a specific
metamorphic grade has to be understood as well as the
availability of fluids (Putnis & Austrheim, 2010). This
understanding is as important as the attainment of mineral
equilibria. Prior to this work, McMullin et al. (2010)
showed with the aid of conventional chemographic projections that within single samples from the Mira terrane,
series of mineral pairs exist at thin section scale that vary
systematically (e.g. in their XMg) suggesting changing temperatures. In a similar way, the systematic variation of
XMg with changing Si content in white mica and chlorite
in the samples of this study can also be interpreted to reflect changing P^T conditions during mineral growth
within single samples.
1867
JOURNAL OF PETROLOGY
VOLUME 54
However, at very low-grade conditions, calculations of
mineral equilibria using multivariant reactions (e.g.Vidal &
Parra, 2000; Willner, 2005; Willner et al., 2009; Jolivet et al.,
2010) have been carried out using compositions of single
pairs of grains in close contact. Such an approach results in a
cluster of points around a P^T path.This application of equilibrium thermodynamics is based on the specific reaction behaviour at very low metamorphic grade, as follows.
Reactions at very low grade are governed (1) by grain
size and (2) by the availability of water, resulting in different mineral compositions pretending non-equilibrium conditions. However, this circumstance allows us to define
locally a series of equilibria during changing P^T conditions (Vidal et al., 2006). The very small grain size combined with the strong compositional variability of the
reaction products as observed, for instance, for white mica
(Fig. 4) reflects high nucleation and low growth rates. As
a consequence nucleation of new grains with a different
composition occurs during changing P^T conditions
rather than only further growth of the existing grains
(Fig. 11). In this way clusters are formed that can recrystallize under low-grade conditions, resulting in larger grains
grown at the expense of smaller grains. In addition, a different composition can form at the rim of larger grains
compared with the core composition. The phenomenon of
increasing grain size with metamorphism is well known
from clastic sedimentary rocks, where the fine-grained
clastic matrix vanishes by recrystallization at temperatures
of 3208C, when grain size increases by an order of magnitude (Brix et al., 2002).
Nevertheless, this specific kinetic behaviour at very lowgrade conditions favouring the crystallization of new
grains during changing P^Tconditions also results in prograde phases remaining metastably at peak P^Tconditions
and during retrograde mineral equilibration. Thus, prograde and retrograde minerals in general cannot be distinguished by their composition in the very low-grade
regime owing to a lack of the zoning that is characteristic
at higher grade (e.g. Kryza et al., 2011).
The reason for nucleation rates exceeding growth rates
is not yet clear. In rocks of very fine grain size close to contact metamorphic aureoles, Roselle et al. (1997) observed
an exponential increase in nucleation rates relative to
growth rates owing to a high reaction affinity, which is a
measure of the departure from equilibrium (excess of free
energy) or the degree of overstepping. Lasaga (1998) proposed that in the case of supersaturation the G of nucleation increases, with an increasing number of embryo
nuclei forming clusters until a maximum with a critical
cluster size is reached and nuclei become stable. The free
energy of nucleation decreases with formation of more
stable nuclei. At very low grade, such supersaturation
most probably prevails owing to the sluggishness of reactions and fluid transport at low temperature.
NUMBER 9
SEPTEMBER 2013
Fig. 11. Sketch showing small phengite crystals (grey) in a thin section that grew in a matrix (light grey) of quartz, feldspar and other
minerals during a water pulse at very low-grade conditions (stage I).
Subsequently at stage II, another water pulse still at very low-grade
conditions caused further growth of phengite with a different composition (dark grey). Stage III represents metamorphism at low-grade
conditions at which small phengite crystals were dissolved and larger
phengite grains grew instead, resulting in a more or less concentric
zoning pattern. During all three stages I^III deformation did not
take place.
Nucleation of new grains is activated by the availability
of water, but partly also by deformation. The apparently
incomplete and variable consumption of protolith minerals
is due to variable access to water. Hence, a variable degree
of preservation of primary minerals is a direct indicator
of a variable reaction progress. Highly variable oxidation
states mirrored, for example, by strongly variable Fe3þ
contents in white mica, reflect locally varying oxidation
1868
WILLNER et al.
VERY LOW- TO LOW-GRADE METAMORPHIC PROCESSES
conditions during pulses of fluid flux. The clustering of reaction products that dominates the metamorphic fabric in
all rock types is a prominent indication of water transport,
at least at thin-section scale. Nucleation of specific clusters
appears to be bound to specific precursor minerals such as
plagioclase, where white mica clusters form as a result of
the high local concentration of Al, or clusters of epidote
form as a result of a high local concentration of Ca. Water
drives the metamorphic reactions and forms abundant
OH-bearing minerals such as white mica, amphibole,
chlorite, pumpellyite, prehnite, stilpnomelane and epidote.
The clustering is controlled by the size and distribution of
precursor phases and hence distances between reactants
and products vary considerably. This influences the local
reaction progress. Putnis & Austrheim (2010) highlighted
the role of hydrous fluids during any metamorphic reaction (except at very high temperature) involving dissolution, material transport and reprecipitation. The extent
of material transport is variable and hence the difference
between metasomatism and metamorphism is rather semantic or a matter of scale. Also, clustering is a metasomatic phenomenon at small scale, resulting in its
characteristic reduction of number of phases. Despite the
local precipitation in clusters and the preferential nucleation of new minerals with differing composition during
changing metamorphic conditions at low grade, the metamorphic fluid causes material transport over a thin section
scale. Hence, the fluid is in contact with protolith minerals
and all phases continuously nucleated during changing
conditions; that is, with the full chemical reservoir represented by the whole-rock composition. Hence at any stage
of the reaction progress, water is interconnecting reactants
and products, ensuring transient equilibrium conditions.
Water at low temperature may be external water liberated by dehydration of neighbouring rocks during
prograde metamorphism. Prior to burial, considerable
amounts of water can be bound in alteration minerals
(e.g. illite, smectite, laumontite) formed during near-surface hydrothermal activity. No such phases were observed
in the rocks of the Mira terrane, because they would not
have survived dehydration, which is particularly intense
at very low grade.
The intensity of dehydration, particularly at very low
and low grade, is demonstrated by the isopleths for water
bound to solids extracted from the calculated results from
the P^T pseudosections (Fig. 10). By far the most water is
released between 250 and 3008C depending on pressure.
This observation was also made for various rock types
such as metapelite, metagreywacke and calcareous rocks
by Massonne & Willner (2008) and Massonne (2010), indicating that this phenomenon specifically occurs at very
low-grade conditions and low pressures. Massonne &
Willner (2008) and Massonne (2010) also showed that the
massive water release at very low grade facilitates
deformation owing to lowering of rock viscosity.
Deformation, in turn, facilitates the reaction progress.
In summary, a continuous adaptation of metamorphic
assemblages to the changing P^Tconditions takes place in
the higher temperature part of the very low-grade metamorphic regime. Each step of the reaction sequence is still
preserved in single samples.
Considerations about the P^T path
It has been shown above that all of the studied samples
from the Mira terrane experienced continuous nucleation
and growth of metamorphic minerals during changing P^
Tconditions within a limited P^T space in the higher temperature part of the very low-grade regime (at approximately 250^3008C and 2^4 kbar). Pressure conditions of
3·5 0·4 kbar were calculated for the Mira terrane samples on the basis of the Si contents of potassic white mica.
As no Si contents higher than 3·41a.p.f.u. were observed,
we conclude that these conditions represent the peak
pressure during the metamorphism. Temperatures of
280 308C were calculated, and we assume that the metamorphic evolution from peak P conditions followed an
exhumation path at these temperatures and, thus, a clockwise P^T trajectory, but we cannot deduce if peak temperatures were reached at the peak pressure, or later.
As no difference was observed in the range of results between samples of different composition (mafic and felsic)
or between samples from a single magmatic belt, we infer
a fairly homogeneous metamorphic evolution for the
entire terrane; that is, the P^T path was similar throughout the Mira terrane, within limits. A partial P^T path
for a single sample at very low-grade conditions can be
defined by a series of local equilibria calculated with multivariant reactions using coexisting mineral compositions
(Vidal & Parra, 2000; Willner, 2005; Jolivet et al., 2010).
On the other hand, the resulting cloud of P^T points
should lie within the P^T space defined by the lowest and
highest Si contents observed in white mica within the field
of the observed prograde, peak and retrograde assemblages; that is, the P^T space where metamorphic reactions continuously occurred to produce white mica of
different compositions. Such fields may have different extents in different samples, even from the same outcrop
area, because reactions at very low-grade conditions
strongly depend on local and temporary water access as
discussed above, as well as on rock composition.
Nevertheless, the white mica P^T fields of all samples
should have an optimal overlap not only at peak metamorphic conditions, but also along parts of the pro- and
retrograde P^T path, if a similar P^T history for all
rocks is assumed.
The entire range of measured phengite compositions
within the P^T fields of the observed (prograde, peak
and retrograde) assemblages are plotted on Fig. 12. For
better resolution separate plots are made for mafic and
1869
JOURNAL OF PETROLOGY
VOLUME 54
NUMBER 9
SEPTEMBER 2013
The prograde metamorphic path can be corroborated
by two fields with prograde stilpnomelane as discussed
above; that is, the field of stilpnomelane in an assemblage
with actinolite, pumpellyite, chlorite and K-feldspar in
mafic sample M84-42 and in an assemblage with phengite,
chlorite and epidote in felsic sample M84-41B. The prograde path between these fields and the peak P^T conditions constrained at 280 308C and 3·5 0·4 kbar can be
described as a heating path with very little pressure increase. By contrast, the subsequent early retrograde P^T
path is characterized by decompression with slight cooling.
It is constrained towards higher temperature by the upper
temperature limit of the retrograde white mica fields in
mafic samples 10Ca13 and AS467 (Fig. 10a) and towards
the lower temperature end by the lower temperature limit
of the retrograde white mica field for felsic sample 10Ca34
(Fig. 10b).
P^T evolution in a collisional setting
Fig. 12. P^T fields corresponding to the observed ranges of white
mica compositions within observed prograde, peak and retrograde assemblages in (a) mafic rocks and (b) felsic rocks. An inferred average
P^T path is indicated by the grey dashed line, which was deduced
from the optimal overlap of these P^T fields. Hatched fields represent
those of prograde stilpnomelane in assemblage (1) with actinolite,
pumpellyite, chlorite and K-feldspar in mafic sample M84-42 and
(2) with phengite, chlorite and epidote in felsic sample M84-41B. The
P^T path passes through the derived average peak conditions at
280 308C and 3·5 0·4 kbar.
felsic samples. Similar results are shown in samples from
all of the magmatic belts of the Mira terrane. This P^T
space, where all white mica (and all other metamorphic
phases) formed, is at 240^3108C; that is, that part where
most water is internally released in both mafic and felsic
rocks and hence is available for chemical reactions (see
Fig. 5). An inferred average P^T path can be inserted in
the area of optimal overlap.
The above derived peak metamorphic conditions suggest
burial to depths of 11^14 km (calculated with 2·8 g cm3
as mean crustal density) and a low transient metamorphic
gradient of 20^258C km1. Such a low geothermal gradient
at peak metamorphic conditions is incompatible with extensional settings, which are characterized by high gradients owing to compaction of isotherms. A relationship
with the collisional assembly of crustal segments represented by the magmatic belts appears more likely.
The maximum depth is comparable with those of
exhumed deeper parts of foreland thrust-and-fold belts
(e.g. Fielitz & Mansy, 1999); that is, of a crustal section
that was only moderately thickened by a collisional event.
No relicts of hydrothermal alteration at or near the surface
during the active magmatic arc stage are preserved.
Nevertheless, such processes probably happened, as evidenced by the pervasive low-dO18 anomaly throughout
the Mira terrane and other parts of West Avalonia (Potter
et al., 2008a, 2008b; Petts et al., 2012). If earlier metamorphic
phases had been formed at 200^2508C, they were erased
by subsequent events. On the other hand, the derived
metamorphic gradient appears relatively high for deeper
levels of a stacked crust. This may be due either to advective heating related to syn-collisional intrusions or to some
near-isobaric heating during thermal relaxation close to
maximum burial, as is typical for frontal accretionary
prisms (Willner et al., 2009). The Devonian granites,
which intruded at 378 þ5/^1 Ma (U/Pb zircon; Barr &
Macdonald, 1992; Bevier et al., 1993), could represent parts
of such syn-collisional magmatism. The timing of this
plutonism, and perhaps the regional metamorphism and
deformation in the Mira terrane, coincides with events in
the now-adjacent Meguma terrane to the south. The
Meguma terrane is characterized by high-temperature^
low-pressure metamorphism, collisional deformation and
substantial syn-collisional plutonism termed the
1870
WILLNER et al.
VERY LOW- TO LOW-GRADE METAMORPHIC PROCESSES
Neoacadian orogeny (400^365 Ma; White & Barr, 2012).
However, the part of the Meguma terrane now juxtaposed
against Mira terrane was situated much farther to the
east at that time, as terrane shuffling and horizontal translations continued by dextral strike-slip motion well into
the Carboniferous (Murphy et al., 2011).
The growth of white mica can be correlated with specific prograde deformation events in the Mira terrane.
The same pattern of white mica compositions (and the resulting metamorphic overprint) is observed in rocks that
are dominated by a foliation related to early D1 folding
(axial plane foliation that may include several generations
of folding) and in rocks, which are dominated by non-coaxial deformation D2 of a later sinistral strike-slip event.
The same pattern prevails also in apparently undeformed
rocks. These observations suggest that an early folding
and thrusting event D1 was mainly responsible for burial
of the rocks, whereas the later strike-slip event (D2)
occurred when the rocks now exposed were close to their
maximum depth of burial. Hence, maximum peak metamorphic conditions correspond to moderate crustal
thickening during a compressive deformational episode.
This episode was followed by a period of strike-slip deformation that was probably responsible for the dispersal
of the Avalonia microplate and final emplacement of the
single slices of magmatic arcs of the Mira terrane; that is,
the final assembly of the crust of the former microplate.
Our findings thus confirm suggestions by Barr et al. (1996)
that the various belts of the Mira terrane represent amalgamated volcanic arcs. It is likely that this model applies
throughout Avalonia, including the ‘type area’ in eastern
Newfoundland, which, like Mira, consists of belts of Neoproterozoic volcanic, sedimentary, and plutonic rocks of
different ages (e.g. O’Brien et al., 1996; Barr & Kerr, 1997;
Sparkes et al., 2005).
The strike-slip event may have been related to the
Acadian collision between Avalonia and the Ganderian
microplate to the north, and/or with the Meguma terrane to the south. Prior burial and D1 deformation
may have occurred during the same orogeny. The late prograde heating with little pressure increase that characterizes the estimated average P^T path (Fig. 12) may be due
to thermal relaxation along a horizontal particle path that
could be inferred for long-term propagation of duplex formation. Alternatively, some heating might have been
related to syn-collisional plutonism. The characteristic
retrograde decompression path with slight cooling might
be due to continuous cooling by newly underthrust rock
units.
CONC LUSIONS
The use of Si isopleths for potassic white mica in pseudosections for geothermobarometry has several advantages
over the calculation of multivariant reactions. The latter
approach uses the composition of minerals in mutual contact. However, although material transport between reactants and products is on the millimetre to centimetre scale
and fluids are in contact with the entire chemical reservoir
of the rocks, precipitation is localized in clusters of small
grains with differing compositions grown during a prevailing high nucleation rate at low temperature. This results
in the frequent neighbourhood of minerals with compositions that are not in equilibrium. Such mineral pairs are
usually discarded because of their strong deviation from
realistic P^T conditions. More importantly, thermometry
in calculated equilibria generally depends on the Fe^Mg
exchange between minerals (e.g. chlorite and white mica).
However, available datasets and activity models used to
calculate these exchange reactions, even when successfully
tested over a wide P^T range, can fail when applied to
other P^Tconditions and rock types. These two shortcomings inherent in the calculation of equilibria with multivariant reactions are overcome by the calculation of
pseudosections using the compositional variation of a
single phase within a calculated P^T field of an observed
assemblage.
The application of P^T pseudosections for geothermobarometry even at very low grade has several advantages:
(1) no coexisting mineral pair with assumed equilibrium
compositions has to be selected; (2) the pseudosections
can be used to study phase relationships; (3) peak P^Tconditions reached can be quantitatively derived with reasonable accuracy; (4) the P^T space, where metamorphic
reactions occurred along a partial P^T trajectory, can be
reconstructed; the overlap of such areas derived from several samples refines the P^T path; (5) important factors
driving metamorphic reactions such as release of water
along the partial P^T path can be evaluated.
However, for the successful application of pseudosections
at very low grade it is essential to use appropriate solid
solutions for all minerals stable at these conditions (e.g.
stilpnomelane, pumpellyite) and not end-member compositions (see Fagereng & Cooper, 2010). The activity models
for white mica, pumpellyite, chlorite, amphibole, epidote
and stilpnomelane, which we employed here, result in reasonable phase relations and compositional contouring of
the P^T pseudosections. Moreover, the specific reaction
kinetics at very low grade, where small grains are steadily
produced with changing composition, implies that the
whole-rock composition does not change during reactions,
as at higher grade with pronounced mineral zoning where
the composition of the cores is shielded from the reacting
rims of minerals.
In the Mira terrane, a region with a relatively uniform
metamorphic overprint, it has been shown for various
lithologies that phengite with Si contents of up to
3·41a.p.f.u. is compatible with formation within a range of
relatively low peak P^T conditions of 3·5 0·4 kbar and
1871
JOURNAL OF PETROLOGY
VOLUME 54
280 308C. In the very low-grade rocks of southeastern
Cape Breton Island most water was released at about
250^3008C and drove nucleation of metamorphic minerals
during changing P^T conditions. Kinetic factors at very
low-grade conditions, such as nucleation rates exceeding
growth rate, led to preservation of a series of mineral
pairs of different compositions in close proximity at thin
section scale, resulting from transient equilibration during
changing P^T conditions. Mineral formation at very low
grade depends on local access of internally generated or
externally derived water. It should be noted that at very
low-grade conditions fabrics related to the initiation of
metamorphic reactions are particularly well preserved.
Almost all metamorphic rocks pass through this grade
during their evolution.
The pronounced availability of water at very low-grade
conditions also lowers the competency and weakens the
rocks, facilitating deformation, which in turn also triggers
mineral formation. A compressional event characterized
by the folding of volcanic and sedimentary strata of the
volcanic arcs was responsible for burial to about 11^14 km
depth. At peak metamorphic conditions a transient metamorphic field gradient of 20^258C km1 prevailed, which
was probably related to collision. Stacking by deformation
similar to that of foreland fold-and-thrust belts can be predicted. The compressional deformation was followed by a
strike-slip event that partly overprinted rocks at depth,
causing the dispersal of the microplate and the final amalgamation of slices of former magmatic arcs to the crust
constituting the Mira terrane as part of the microplate of
Avalonia. The derived metamorphic conditions describe
the thermal state during this collisional assembly. Evidence
for the timing of the metamorphic processes and related
deformation will be presented in a forthcoming paper.
AC K N O W L E D G E M E N T S
David McMullin helped with insightful comments on earlier drafts of this paper. Critical reviews by J. Allaz, J.
Schumacher, C. van Staal and an anonymous reviewer as
well as careful editorial handling by R. Giere¤ improved
the paper substantially. We thank all of them for their
contributions.
F U N DI NG
This project was financed by Deutsche Forschungsgemeinschaft (grants Ma1126-27 and Wi847-9) to H.-J.M.
and A.P.W. Geological mapping by S.M.B. and C.E.W. in
southeastern Cape Breton Island during which some of
the studied samples were collected was funded by the Geological Survey of Canada through the 1984^1989 Canada^
Nova Scotia Mineral Development Agreement and the
1990^1992 Canada^Nova Scotia Cooperation Agreement,
NUMBER 9
SEPTEMBER 2013
as well as by research grants to S.M.B. from the Natural
Sciences and Engineering Research Council of Canada.
S U P P L E M E N TA RY DATA
Supplementary data for this paper are available at Journal
of Petrology online.
R E F E R E NC E S
Banno, S. (1998). Pumpellyite^actinolite facies of the Sanbagawa
metamorphism. Journal of Metamorphic Geology 16, 117^128.
Barr, S. M. (1993). Geochemistry and tectonic setting of late
Precambrian volcanic and plutonic rocks in southeastern Cape
Breton Island, Nova Scotia. Canadian Journal of Earth Sciences 30,
1147^1154.
Barr, S. M. & Kerr, A. (1997). Late Precambrian plutons in the Avalon
terrane of New Brunswick, Nova Scotia, and Newfoundland. In:
Sinha, A. K., Whalen, J. B. & Hogan, J. P. (eds) The Nature of
Magmatism in the Appalachian Orogen. Geological Society of America,
Memoirs 191, 45^74.
Barr, S. M. & Macdonald, A. S. (1992). Devonian plutons in southeastern Cape Breton Island, Nova Scotia. Atlantic Geology 28, 101^113.
Barr, S. M. & Raeside, R. P. (1989). Tectonostratigraphic terranes in
Cape Breton Island, Nova Scotia. Implications for the configuration of the northern Appalachian Orogen. Geology 17, 822^825.
Barr, S. M., White, C. E. & Macdonald, A. S. (1996). Stratigraphy, tectonic setting, and geologic history of Late Precambrian volcanic^sedimentary^plutonic belts in southeastern Cape Breton Island, Nova Scotia.
Geological Survey of Canada Bulletin 468, 84 p.
Barr, S. M., Dunning, G. R., Raeside, R. P. & Jamieson, R. A. (1990).
Contrasting U^Pb ages from plutons in the Bras d’Or and Mira
terranes of Cape Breton Island, Nova Scotia. Canadian Journal of
Earth Sciences 27, 1200^1208.
Berman, R. G. (1988). Internally-consistent thermodynamic data
for minerals in the system Na2O^K2O^CaO^MgO^FeO^
Fe2O3^Al2O3^SiO2^TiO2^H2O^CO2. Journal of Petrology 29,
445^522.
Bevier, M. L., Barr, S. M., White, C. E. & Macdonald, A. S. (1993).
U^Pb geochronologic constraints on the volcanic evolution of the
Mira (Avalon) terrane, southeastern Cape Breton Island, Nova
Scotia. CanadianJournal of Earth Sciences 30, 1^10.
Brix, M. R., Sto«ckhert, B., Seidel, E., Theye, T., Thomson, S. N. &
Ku«ster, M. (2002). Thermobarometric data from a fossil zircon partial annealing zone in high pressure^low temperature rocks of eastern and central Crete, Greece. Tectonophysics 349, 309^326.
Brown, E. H. (1975). A petrogenetic grid for reactions producing biotite and other Al^Fe^Mg-silicates in the greenschist facies. Journal
of Petrology 16, 258^271.
Cho, M. & Liou, J. G. (1987). Prehnite^pumpellyite to greenschist
facies transition in the Karmutsen mafic rocks, Vancouver Island,
B.C. Journal of Petrology 28, 417^443.
Cho, M., Liou, J. G. & Maruyama, S. (1986). Transition from the zeolite to prehnite^pumpellyite facies in the Karmutsen mafic rocks,
Vancouver Island, British Columbia. Journal of Petrology 27,
467^494.
Connolly, J. A. D. (1990). Multivariable phase diagrams; an algorithm
based on generalized thermodynamics. American Journal of Science
290, 666^718.
Connolly, J. A. D. (2005). Computation of phase equilibria by linear
programming: A tool for geodynamic modeling and its application
1872
WILLNER et al.
VERY LOW- TO LOW-GRADE METAMORPHIC PROCESSES
to subduction zone decarbonation. Earth and Planetary Science Letters
236, 524^541.
Cruciani, G., Franceschelli, M. I. & Massonne, H.-J. (2011). Low-temperature metamorphic evolution of a pre-Variscan gabbro: a case
study from the Palaeozoic basement of northwest Sardinia, Italy.
Mineralogical Magazine 75, 2793^2812.
Dingeldey, R. D., Dallmeyer, F., Koller, H.-J. & Massonne, H.-J.
(1997). P^T^t history of the Lower Austroalpine nappe complex in
the ‘Tarntaler Berge’ NW of the Tauern Window: implications for
the geotectonic evolution of the central Eastern Alps. Contributions
to Mineralogy and Petrology 129, 1^19.
Fagereng, —. & Cooper, A. F. (2010). The metamorphic history of rocks
buried, accreted and exhumed in an accretionary prism: an example from the Otago Schist, New Zealand. Journal of Metamorphic
Geology 28, 935^954.
Fielitz, W. & Mansy, J.-L. (1999). Pre- and synorogenic burial
metamorphism in the Ardennes and neighbouring areas (Rhenohercynian zone, central European Variscides). Tectonophysics 309,
227^256.
Frey, M., De Capitani, C. & Liou, J. G. (1991). A new petrogenetic
grid for low-grade mafic rocks. Journal of Metamorphic Geology 9,
497^509.
Giles, P. S., Naylor, R. D., Teniere, P. J., White, C. E., Barr, S. M.,
DeMont, G. J. & Force, E. R. (2010). Bedrock geology map of the
Port Hawkesbury area, Part of NTS Sheets 11F/06, 11F/07, 11F/10,
11F/11 and 11F/15, Inverness, Richmond, Guysborough, and
Antigonish Counties, Nova Scotia (scale 1:50 000)Nova Scotia
Department of Natural Resources, Mineral Resources Branch,
Open File Map Me 2010-006.
Hibbard, J. P., van Staal, C. R., Rankin, D. & Williams, H. (2006).
Lithotectonic map of the Appalachian orogen (north), Canada^
United States of America. Geological Survey of Canada, Map
2041A, scale 1:1500 000, 1 sheet.
Holland, T. J. B. & Powell, R. (1996). Thermodynamics of order^disorder in minerals. 2. Symmetric formalism applied to solid solutions.
American Mineralogist 81, 1425^1437.
Holland, T. J. B. & Powell, R. (1998). An internally consistent thermodynamic data set for phases of petrological interest. Journal of
Metamorphic Geology 16, 309^343.
Holland, T. J. B. & Powell, R. (2003). Activity^composition relations
for phases in petrological calculations: an asymmetric multicomponent formulation. Contributions to Mineralogy and Petrology 145,
492^501.
Jolivet, L., Trotet, F., Monie¤, P., Vidal, O., Goffe¤, B., Labrousse, L.,
Agard, P. & Ghorbal, B. (2010). Along-strike variations of P^Tconditions in accretionary wedges and syn-orogenic extension, the
HP^LT Phyllite^Quartzite Nappe in Crete and the Peloponnese.
Tectonophysics 480, 133^148.
Kryza, R., Willner, A. P., Massonne, H.-J., Muszyn¤ski, A. &
Schertl, H.-P. (2011). Blueschist-facies metamorphism in the
Kaczawa Mountains (Sudetes, SW Poland) of the CentralEuropean Variscides: P^Tconstraints by a jadeite-bearing metatrachyte. Mineralogical Magazine 75, 241^263.
Lasaga, A. C. (1998). KineticTheory in the Earth Sciences. Princeton Series in
Geochemistry. Princeton, NJ: Princeton University Press, 811 p.
Liou, J. G. (1971a). P^T stabilities of laumontite, wairakite, lawsonite,
and related minerals in the system CaAl2Si2O8^SiO2^H2O.
Journal of Petrology 12, 379^411.
Liou, J. G. (1971b). Synthesis and stability relations of prehnite,
Ca2Al2Si3O12(OH)2. American Mineralogist 56, 507^531.
Liou, J. G., Maruyama, S. & Cho, M. (1985). Phase equilibria and
mineral parageneses of mafic rocks in low-grade metamorphism.
Mineralogical Magazine 49, 321^333.
Liou, J. G., Maruyama, S. & Cho, M. (1987). Very low-grade metamorphism of volcanic and volcaniclastic rocksçmineral assemblages and mineral facies. In: Frey, M. (ed.) Low Temperature
Metamorphism. Glasgow: Blackie, pp. 59^113.
Macdonald, A. S. & Barr, S. M. (1993). Geological setting and depositional environment of the Stirling Group of southeastern Cape
Breton Island, Nova Scotia. Atlantic Geology 29, 137^147.
Massonne, H.-J. (1995). P^T evolution of metavolcanics from the
southern Taunus mountains. In: Dallmeyer, R. D., Franke, W. &
Weber, K. (eds) Pre-Permian Geology of Central and Eastern Europe.
Berlin: Springer, pp. 132^137.
Massonne, H.-J. (2010). Phase relations and dehydration behaviour of
calcareous sediments at very-low to low grade metamorphic conditions. Periodico di Mineralogia 79, 21^43.
Massonne, H.-J. & Schreyer, W. (1986). High-pressure syntheses and Xray properties of white micas in the system K2O^MgO^Al2O3^
SiO2^H2O. NeuesJahrbuch fu«r Mineralogie, Abhandlungen 153, 177^215.
Massonne, H.-J. & Willner, A. P. (2008). Phase relations and dehydration behaviour of psammopelite and mid-ocean ridge basalt at
very low-grade to low-grade metamorphic conditions. European
Journal of Mineralogy 20, 867^879.
McMullin, D. W. A., Barr, S. M. & Raeside, R. P. (2010). Very lowand low-grade metamorphism of mafic volcanic rocks of the Mira
terrane (Avalonia), southeastern Cape Breton Island, Nova Scotia.
Atlantic Geology 46, 95^126.
Murphy, J. B., Waldron, J. W. F., Kontak, D. J., Pe-Piper, G. &
Piper, D. J. W. (2011). Minas Fault Zone: Late Paleozoic history of
an intra-continental orogenic transform fault in the Canadian
Appalachians. Journal of Structural Geology 33, 312^328.
Nitsch, K. H. (1971). Stabilita«tsbeziehungen yon Prehnit- und
Pumpellyit-haltigen Paragenesen. Contributions to Mineralogy and
Petrology 30, 240^260.
O’Brien, S. J., O’Brien, B. H., Dunning, G. R. & Tucker, R. D. (1996).
Late Neoproterozoic Avalonian and related peri-Gondwanan
rocks of the Newfoundland Appalachians. In: Nance, R. D. &
Thompson, M. (eds) Avalonian and Related Peri-Gondwanan
Terranes of the Circum-North Atlantic. Geological Society of America,
Special Papers 304, 9^27.
Parra, T., Vidal, O. & Jolivet, L. (2002). Relation between the intensity
of deformation and retrogression in blueschist metapelites of Tinos
Island (Greece) evidenced by chlorite^mica local equilibria. Lithos
63, 41^66.
Petts, D., Longstaffe, F. J., Potter, J., Barr, S. M. & White, C. E. (2012).
Regional hydrothermal alteration and 18O-depletion of the ca. 620
Ma Huntington Mountain pluton and related rocks, Cape Breton
Island, Canada. Atlantic Geology 48, 54^69.
Potter, J., Longstaffe, F. J. & Barr, S. M. (2008a). Regional 18O-depletion of Neoproterozoic igneous rocks from Avalonia, Cape Breton
Island and southern New Brunswick, Canada. Geological Society of
America Bulletin 120, 347^367.
Potter, J., Longstaffe, F. J., Barr, S. M., Thompson, M. D. & White, C.
E. (2008b). Altering Avalonia: oxygen isotopes and terrane distinction in the Appalachian peri-Gondwanan realm. Canadian Journal
of Earth Sciences 45, 815^825.
Powell, R. & Holland, T. (1999). Relating formulations of the thermodynamics of mineral solid solutions: Activity modeling of pyroxenes, amphiboles, and micas. American Mineralogist 84, 1^14.
Putnis, A. & Austrheim, H. (2010). Fluid-induced processes: metasomatism and metamorphism. Geofluids 10, 254^269.
Reynolds, P. H., Barr, S. M. & White, C. E. (2009). Provenance of detrital muscovite in Cambrian Avalonia of Maritime Canada:
40
Ar/39Ar ages and chemical compositions. Canadian Journal of
Earth Sciences 46, 169^180.
1873
JOURNAL OF PETROLOGY
VOLUME 54
Rieder, M., Cavazini, G., D’Yakonov, Y. S., Frank-Kamenetskii, V. A.,
Guggenheim, S., Koval, P. V., Mu«ller, G., Neiva, A. M. R.,
Radoslovich, E. W., Robert, J.-L., Sassi, F. P., Takeda, H., Weiss, Z.
& Wones, D. R. (1998). Nomenclature of the micas. Canadian
Mineralogist 36, 905^912.
Robinson, D., Bevins, R. E. & Rubinstein, N. (2005). Subgreenschist
facies metamorphism of mafic rocks from the Precordillera terrane
of western Argentina: constraints on the later stages of accretion
onto Gondwana. EuropeanJournal of Mineralogy 17, 441^452.
Roselle, G. T., Baumgartner, L. P. & Chapman, J. A. (1997).
Nucleation-dominated crystallisation of forsterite in the Ubehebe
contact aureole, California. Geology 25, 823^826.
Schiffman, P. & Liou, J. G. (1980). Synthesis and stability relations of
Mg^Al pumpellyite, Ca4Al5MgSi6O21(OH)7. Journal of Petrology
21, 441^474.
Schreyer, W. & Abraham, K. (1978). Prehnite/chlorite and actinolite/
epidote bearing mineral assemblages in the metamorphic igneous
rocks of La Helle and Challes, Venn^Stavelot Massif, Belgium.
Annales de la Socie¤ te¤ Ge¤ ologique de Belgique 101, 227^241.
Sparkes, G. W., O’Brien, S. J., Dunning, G. R. & Dube, B. (2005). U^Pb
geochronological constraints on the timing of magmatism, epithermal alteration, and low-sulphidation gold mineralization, eastern
Avalon Zone, Newfoundland, Current Research (2005). Newfoundland
and Labrador Department of Natural Resources, Report 05-1, 115^130.
Vidal, O. & Parra, T. (2000). Exhumation paths of high pressure
metapelites obtained from local equilibria for chlorite^phengite assemblages. Geological Journal 35, 139^161.
NUMBER 9
SEPTEMBER 2013
Vidal, O., De Andrade, V., Lewin, E., Munoz, M., Parra, T. &
Pascarelli, S. (2006). P^T^deformation^Fe3þ/Fe2þ mapping at thin
section scale and comparison with XANES mapping: application
to a garnet-bearing metapelite from the Sambagawa metamorphic
belt (Japan). Journal of Metamorphic Geology 24, 669^683.
White, C. E. & Barr, S. M. (2012). The new Meguma: stratigraphy, metamorphism, paleontology and provenance. Field Trip B-5, Geological
Association of Canada^Mineralogical Association of Canada Field Trip
Guide, 68 p.
White, C. E., Barr, S. M. & Ketchum, J. W. F. (2003). New age controls on rock units in pre-Carboniferous basement blocks in southwestern Cape Breton Island and adjacent mainland Nova Scotia.
In: MacDonald, D. R. (ed.) Nova Scotia Department of Natural
Resources, Minerals and Energy Branch, Report of Activities 2002. Report
ME 2003-1, pp. 163^178.
Willner, A. P. (2005). Pressure^temperature evolution of an
Upper Paleozoic paired metamorphic belt in Central Chile
(348^35830’S). Journal of Petrology 46, 1805^1833.
Willner, A. P., Herve, F. & Massonne, H.-J. (2000). Mineral chemistry
and pressure^temperature evolution of two contrasting levels
within an accretionary complex in the Chonos Metamorphic
Complex, Southern Chile. Journal of Petrology 41, 309^330.
Willner, A. P., Sepu¤lveda, F. A., Herve¤, F., Massonne, H.-J. &
Sudo, M. (2009). Conditions and timing of pumpellyite^actinolite
facies metamorphism in the Early Mesozoic frontal accretionary
prism of the Madre de Dios Archipelago (50820’S; southern
Chile). Journal of Petrology 50, 2127^2155.
1874