Skip Navigation


Journal of Petrology Advance Access originally published online on January 21, 2005
Journal of Petrology 2005 46(5):893-919; doi:10.1093/petrology/egi004
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow All Versions of this Article:
46/5/893    most recent
egi004v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (3)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by SLAGSTAD, T.
Right arrow Articles by CULSHAW, N. G.
Right arrow Search for Related Content
GeoRef
Right arrow GeoRef Citation
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2005. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oupjournals.org

Formation, Crystallization, and Migration of Melt in the Mid-orogenic Crust: Muskoka Domain Migmatites, Grenville Province, Ontario

T. SLAGSTAD*, R. A. JAMIESON and N. G. CULSHAW

DEPARTMENT OF EARTH SCIENCES, DALHOUSIE UNIVERSITY, HALIFAX, N.S., CANADA, B3H 3J5

RECEIVED APRIL 30, 2004; ACCEPTED DECEMBER 1, 2004


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING
 STUDY AREA
 MUSKOKA MIGMATITES
 PETROGRAPHY AND GEOCHEMISTRY
 MIGMATITE PETROGENESIS
 MELT MIGRATION
 CONCLUSIONS
 REFERENCES
 
Migmatitic orthogneisses in the Muskoka domain, southwestern Grenville Province, Ontario, formed during the Ottawan stage (c. 1080–1050 Ma) of the Grenvillian orogeny. Stromatic migmatites are volumetrically dominant, comprising granodioritic gneisses with 2–5 cm thick granitic leucosomes, locally rimmed by thin melanosomes, that constitute 20–30 vol. %, and locally 40–50 vol. %, of the outcrops. Patch migmatites in dioritic gneisses form large (>10 m) pinch-and-swell structures within the stromatic migmatites, and consist of decimetre-scale, irregular patches of granitic leucosome, surrounded by medium-grained hornblende–plagioclase melanosomes interpreted as restite. The patches connect to larger networks of zoned pegmatite dykes. Petrographic and geochemical evidence suggests that the patch leucosomes formed by 20–40% fluid-present, equilibrium melting of the dioritic gneiss, followed by feldspar-dominated crystallization. The dyke networks may have resulted from hydraulic fracturing, probably when the melts reached water saturation during crystallization. Field and geochemical data from the stromatic migmatites suggest a similar petrogenesis to the patch migmatites, but with significant additions of externally derived melts, indicating that they acted as conduits for melts derived from deeper structural levels within the orogen. We hypothesize that the Muskoka domain represents a transfer zone for melts migrating to higher structural levels during Grenvillian deformation.

KEY WORDS: migmatite geochemistry; partial melting; melt crystallization; melt transport; Grenville orogen


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING
 STUDY AREA
 MUSKOKA MIGMATITES
 PETROGRAPHY AND GEOCHEMISTRY
 MIGMATITE PETROGENESIS
 MELT MIGRATION
 CONCLUSIONS
 REFERENCES
 
Partial melting in the middle to lower crust may take place in response to dehydration of hydrous minerals such as muscovite, biotite, and amphibole (e.g. Thompson, 1982Go; Nédélec et al., 1993Go) or influx of externally derived hydrous fluids (e.g. Brown, 1979Go; Montel et al., 1992Go). In many cases, geochemical evidence suggests that the melts may undergo magmatic differentiation, either in situ or while migrating from the site of melting (e.g. Sawyer, 1987Go; Sawyer et al., 1999Go). Most migmatites form in regions of orogenically thickened crust and a number of field and experimental studies have indicated a significant potential for melt migration during deformation (e.g. Collins & Sawyer, 1996Go; Brown & Solar, 1998Go; Sawyer et al., 1999Go; Barraud et al., 2004Go). Partial melting and melt migration are important in the compositional differentiation of the continental crust (e.g. Fyfe, 1973Go; Vielzeuf et al., 1990Go) and may influence orogenic evolution by lowering crustal strength (e.g. Hollister & Crawford, 1986Go; Davidson et al., 1994Go; Beaumont et al., 2001Go; Jamieson et al., 2002Go). Understanding how melts form, evolve, and move during orogenesis is therefore fundamental for understanding the geochemical structure of the continental crust (e.g. Sandiford & McLaren, 2002Go), the formation of middle to shallow crustal granitoid plutons and batholiths (e.g. Collins & Sawyer, 1996Go; Solar & Brown, 2001Go), and structural evolution during orogenesis (e.g. Hollister, 1993Go; Beaumont et al., 2004Go).

The Muskoka domain in the Central Gneiss Belt, SW Grenville Province, Ontario (Fig. 1), is characterized by leucosome-rich migmatitic orthogneisses (Timmermann et al., 2002Go; Slagstad, 2003Go) that formed at mid-crustal levels at the peak of Grenvillian (Ottawan) contractional deformation and high-grade metamorphism (Timmermann et al., 1997Go; Slagstad et al., 2004bGo). The Muskoka domain, therefore, presents a natural laboratory in which the effects of partial melting in the middle orogenic crust can be studied. Despite several recent studies (Timmermann, 1998Go; McMullen, 1999Go; Timmermann et al., 2002Go), a number of questions regarding the petrogenesis of these rocks remain. Unresolved problems include the role of fluid during partial melting (e.g. Timmermann, 1998Go; Timmermann et al., 2002Go), whether the abundant leucosomes represent in situ or externally derived melts, the extent to which leucosome compositions were modified by fractionation, and the possible relationships between melting, melt migration, and orogenic evolution. The purpose of this paper is to present field, petrographic, and geochemical data from the Muskoka domain migmatites in an attempt to answer these questions, and thereby contribute to a better understanding of the causes and consequences of partial melting in the Grenville orogen and, by analogy, other large collisional orogens.



View larger version (47K):
[in this window]
[in a new window]
 
Fig. 1. (a) Location and geology of the study area. Bold barbed lines indicate major thrust boundaries and bold non-barbed lines indicate other important structural contacts. Fine lines represent domain boundaries after Davidson (1984)Go and Davidson & van Breemen (1988)Go, with minor modifications by Culshaw et al. (1997)Go. (b) Crustal-scale cross-section from the Grenville Front to the Central Metasedimentary Belt, modified after Culshaw et al. (1997)Go; based on geological and structural data, discussed by Culshaw et al. (1997)Go, and seismic reflection profiles AGL-31, 32, 33 and GLIMPCE-J (Green et al., 1988Go; White et al., 1994Go). Composite Arc Belt defined by Carr et al. (2000)Go. GFTZ, Grenville Front Tectonic Zone; SSZ, Shawanaga Shear Zone; PSSZ, Parry Sound Shear Zone; UR, Upper Rosseau domain; LR, Lower Rosseau domain; UGH, Upper Go Home domain; LGH, Lower Go Home domain; CMBBZ, Central Metasedimentary Belt boundary thrust zone; CMB, Central Metasedimentary Belt; PSD, Parry Sound domain. In cross-section: SH, Shawanaga domain; MR-MS, Moon River–Muskoka domains; GH-RA, Go Home–Rosseau domains.

 

    GEOLOGICAL SETTING
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING
 STUDY AREA
 MUSKOKA MIGMATITES
 PETROGRAPHY AND GEOCHEMISTRY
 MIGMATITE PETROGENESIS
 MELT MIGRATION
 CONCLUSIONS
 REFERENCES
 
The Central Gneiss Belt (Wynne-Edwards, 1972Go) of the Grenville Province, Ontario, exposes large tracts of upper amphibolite- to granulite-facies migmatitic gneisses of dominantly Laurentian affinity, formed by continent–arc(s) and/or continent–continent collision(s) along the southeastern margin of Laurentia between 1200 and 1000 Ma (e.g. Rivers, 1997Go; Carr et al., 2000Go). In Ontario (Fig. 1), the Central Gneiss Belt comprises pre-1350 Ma high-grade gneisses, formed along the active southeastern margin of the pre-Grenvillian Laurentian craton (Culshaw & Dostal, 1997Go; Rivers, 1997Go; Rivers & Corrigan, 2000Go; Culshaw & Dostal, 2002Go; Slagstad et al., 2004aGo), and assembled during multiple stages of NW-directed convergence (Culshaw et al., 1997Go). The earliest stage of convergence, termed the ‘Elzevirian orogeny’ (c. 1200–1160 Ma; Moore & Thompson, 1980Go), involved assembly of arc terranes offshore from, or at the distal southeastern edge of, Laurentia (e.g. Carr et al., 2000Go). In the study area, the main collisional event, termed the ‘Ottawan orogeny’ (Moore & Thompson, 1980Go), began some time after 1120 Ma and resulted in emplacement of post-1400 Ma arc and back-arc assemblages of the Central Metasedimentary Belt over the Central Gneiss Belt along a crustal-scale shear zone termed the Central Metasedimentary Belt boundary thrust zone (Hanmer & McEachern, 1992Go). The Ottawan orogeny was associated with foreland-propagating convergence and widespread upper amphibolite- to granulite-facies metamorphism and partial melting at 1085–1035 Ma (e.g. Culshaw et al., 1997Go; Timmermann et al., 1997Go; McMullen, 1999Go; Carr et al., 2000Go; Wodicka et al., 2000Go; Slagstad, 2003Go). This orogenic phase was followed by regionally significant ductile extension at c. 1020 Ma (Ketchum et al., 1998Go), and there is evidence for both earlier and later extensional deformation in some parts of the region (van der Pluijm & Carlson, 1989Go; Timmermann et al., 2002Go). A final stage of orogen-wide convergence affected the Grenville Front Tectonic Zone at c. 1000 Ma (e.g. Haggart et al., 1993Go; Krogh, 1994Go).


    STUDY AREA
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING
 STUDY AREA
 MUSKOKA MIGMATITES
 PETROGRAPHY AND GEOCHEMISTRY
 MIGMATITE PETROGENESIS
 MELT MIGRATION
 CONCLUSIONS
 REFERENCES
 
The Muskoka domain constitutes the uppermost structural level of the Central Gneiss Belt (Culshaw et al., 1983Go), in the immediate footwall to the Central Metasedimentary Belt boundary thrust zone (Figs 1 and 2). East- to SE-dipping, grey migmatitic orthogneisses with protolith compositions ranging from gabbro to granite, with granodiorite predominant (Timmermann, 1998Go; Slagstad et al., 2004aGo), and major and trace element compositions suggest formation in a continental magmatic arc (McMullen, 1999Go; Slagstad et al., 2004aGo). Radiometric (U–Pb zircon) ages from the grey gneisses range from c. 1480 to 1430 Ma, but most fall in a narrow range between c. 1450 and 1460 Ma (Timmermann et al., 1997Go; Nadeau & van Breemen, 1998Go; McMullen, 1999Go). The Muskoka domain is bounded above by the Central Metasedimentary Belt boundary thrust zone, which contains evidence of thrusting between c. 1080 and 1050 Ma (van Breemen & Hanmer, 1986Go; McEachern & van Breemen, 1993Go; Burr & Carr, 1994Go). Adjacent to the Muskoka domain, the Central Metasedimentary Belt boundary thrust zone includes reworked migmatites similar to those in its immediate footwall. Peak metamorphic conditions of 750–850°C at 0·9–1·15 GPa (Fig. 3) were attained at c. 1080–1060 Ma (Timmermann, 1998Go; Timmermann et al., 2002Go; S. Gagné & J. Hawken, unpublished data, 2002). Timmermann et al. (1997)Go dated a foliation-parallel leucosome in the Muskoka domain at 1064 ± 18 Ma (U–Pb zircon, TIMS), and this age has recently been corroborated by sensitive high-resolution ion microprobe (SHRIMP) dating (1067 ± 9 Ma; Slagstad et al., 2004bGo).



View larger version (93K):
[in this window]
[in a new window]
 
Fig. 2. Geological map of the Muskoka domain showing dominant lithological units and the locations of the investigated outcrops. UTM coordinates (UTM zone 17) of the investigated outcrops: o/c M03083 (Parkersville): 5002 400N, 632 700E; o/c M30072 [GenBank] (Bonnie Lake): 4995 500N, 636 500E; o/c M2707 (Wren Lake): 5005 430N, 668 550E.

 


View larger version (25K):
[in this window]
[in a new window]
 
Fig. 3. PT diagram summarizing some relevant melting reactions and results from experimental melting studies. For clarity, divariant reactions are illustrated using a single line. Continuous lines denote fluid-absent melting reactions, dashed lines denote fluid-present reactions. (1a) and (1b) illustrate the reaction Qtz + Pl ± Kfs + Fl = M in the presence of a pervasive H2O-rich (aH2O = 1) and CO2-rich (aH2O = 0·1) fluid, respectively (Stevens & Clemens, 1993Go). (2) Water-saturated melting of tonalite Bt + Pl + Qtz + Fl (water-rich) = Hbl + M (Büsch et al., 1974Go). (3) Partial melting of a biotite–plagioclase–quartz gneiss (Gardien et al., 2000Go). (3a) Biotite-out with no H2O added to the system. At pressures between 1 and 2 GPa and no H2O added, <4 vol. % melt was produced at T = 900°C. (3b) Biotite-out with >2 wt % H2O added. At pressures between 1 and 2 GPa and 4 wt % H2O added, between 50 and 60 vol. % melt was produced at T as low as 800°C. (4) Fluid-absent melting Bt + Pl + Qtz = Opx + Kfs + M (Stevens & Clemens, 1993Go). (5) Fluid-absent biotite dehydration melting of a tonalite at 1 GPa (Rutter & Wyllie, 1988Go). Between 825 and 900°C (indicated by double-pointed arrow) the biotite melting reaction produced 20% melt. (6) Partial melting of dacite at 1 GPa (Conrad et al., 1988Go). Boxes with vertical (6a), NE–SW diagonal (6b), and NW–SE diagonal (6c) lines represent the solidus at 1 GPa with aH2O = 0·75, 0·5, and 0·25, respectively. (7) Beginning of biotite dehydration melting in a biotite gneiss according to the reaction Bt + Pl + Qtz = M + Opx + Ox (Patiño Douce & Beard, 1995Go). (8a) Bt + Qtz + Kfs = Opx + M (reaction [V] of Peterson & Newton (1989)Go), and (8b) Bt + Qtz + Kfs + Fl (water-rich) = M (reaction [O] of Peterson & Newton (1989)Go). (9) Wet basalt solidus (Green, 1982Go). Boxes A (shaded) and B (white) indicate range of PT conditions in the Muskoka domain determined by Timmermann et al. (2002)Go and S. Gagné & J. Hawken (unpublished data, 2002), respectively. Mineral abbreviations are after Kretz (1983)Go; M, melt; Fl, fluid; Ox, oxides.

 

    MUSKOKA MIGMATITES
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING
 STUDY AREA
 MUSKOKA MIGMATITES
 PETROGRAPHY AND GEOCHEMISTRY
 MIGMATITE PETROGENESIS
 MELT MIGRATION
 CONCLUSIONS
 REFERENCES
 
We distinguish between two types of migmatite, referred to as stromatic and patch migmatite, based on leucosome morphology and type of host rock. The stromatic migmatites typically have light grey, quartz monzodioritic to granodioritic mesosomes, and are by far the dominant type of migmatite in the Muskoka domain. The patch migmatites have dioritic mesosomes, and are typically found as outcrop-scale bodies within the stromatic migmatites. Because of their generally lower strain, these rocks provide some constraints on melt generation and transport processes indeterminable from the stromatic migmatites.

Three representative outcrops were investigated in detail (Fig. 2). The first two outcrops are located south of Parkersville (o/c M03083) and Bonnie Lake (o/c M30072 [GenBank] ), in the western part of the Muskoka domain, and display stromatic and patch migmatites, respectively. The third outcrop is located at Wren Lake (o/c M2707), in the eastern part of the Muskoka domain, and comprises both stromatic and patch migmatites separated by sharp, foliation-parallel contacts. Most of the investigated migmatites preserve pre-migmatization compositional banding and foliation and can be classified as metatexite (Brown, 1973Go). The terms ‘leucosome’, ‘mesosome’, and ‘melanosome’ are used below as defined by Ashworth (1985)Go.

Patch migmatites
The patch migmatites are morphologically more varied than the stromatic migmatites and range from small, irregular patches, <10–20 cm across, to dyke networks with individual dykes 30–40 cm across and 3–4 m long (Fig. 4a). For clarity, we refer to these morphological end-members as patches and dykes, respectively, but emphasize that the field, petrographic, and geochemical data strongly suggest that they are part of a continuum. Most of the patches contain coarse-grained granitic to pegmatitic-textured, strongly recrystallized leucosome, containing subhedral to euhedral, poikiloblastic hornblende up to 5 cm long (Fig. 4b and c). The inclusions are dominantly plagioclase, with minor quartz and biotite. The mineralogy of the patches varies from hornblende-dominated, with seam-like networks of granitic material, to leucosome-dominated, with little or no hornblende. The dykes are pegmatitic and typically zoned, with quartz-rich cores surrounded by K-feldspar-rich pegmatite (Fig. 4e), resembling typical pegmatite zonation and suggesting crystallization from relatively water-rich melts (e.g. Jahns & Burnham, 1969Go; Symmes & Ferry, 1995Go). Large hornblende crystals, morphologically resembling those in the patches, are common along the dyke margins (Fig. 4d) or form local accumulations within the dykes. Commonly, patches and dykes are connected by thin seams of granitic material oriented subparallel to the foliation in the dioritic mesosome. Immediately adjacent to the patches, and locally to the dykes, the host rock is enriched in subhedral, medium-grained (c. 0·5 cm) hornblende (Fig. 4c). These zones are interpreted as melanosome and referred to as patch melanosome below. Despite clear petrographic and geochemical differences, discussed below, the patch melanosomes are difficult to distinguish from the dioritic mesosome in the field. The distribution of melanosome is, therefore, unclear; in particular, some of the dykes appear to be locally intrusive and may lack associated melanosomes. Leucosome abundances are difficult to estimate in the patch migmatites because of their irregular distribution, but appear to be lower than in the stromatic migmatites.



View larger version (136K):
[in this window]
[in a new window]
 
Fig. 4. (a) Network of dykes and patches in a patch migmatite; hammer for scale c. 30 cm long. (b) Hornblende-rich patch comprising subequal amounts of granitic leucosome and coarse hornblende; coin for scale c. 2·8 cm across. (c) Close-up of hornblende-rich patch showing coarse, subhedral, poikiloblastic hornblende and leucosome. (d) Close-up of pegmatite dyke with hornblende porphyroblasts, similar to those depicted in (b) and (c), concentrated along the margins; coin for scale c. 2·8 cm across. (e) Quartz-rich centre of zoned pegmatite dyke; coin for scale c. 2·4 cm across. (f) Foliation in the stromatic migmatites wrapping around a body of patch migmatite, forming a pinch-and-swell structure; field of view is c. 5 m across.

 
In outcrops containing both stromatic and patch migmatites (e.g. Wren Lake outcrop), the foliation in the former can be seen wrapping around the latter, forming a pinch-and-swell structure (Fig. 4f). This observation suggests significant rheological contrast between the two types of migmatite. Unlike in the stromatic migmatites described below, there is no evidence that the location of the patch leucosomes is structurally controlled. We ascribe this difference to the higher strength of the dioritic host rocks relative to the surrounding granodioritic orthogneisses, which may have shielded the patch migmatites from the effects of deformation (see also Sawyer, 1999Go). None of the patches or dykes were observed to cut the patch–stromatic migmatite contact, thus the relative age of stromatic and patch migmatites is not known.

Stromatic migmatite
The stromatic migmatites are characterized by a layered structure made up of 0·5–5 cm thick layers of medium- to coarse-grained, strongly recrystallized granitic leucosome, typically a few metres long. The granodioritic host rock is commonly banded with centimetre-scale layers of dark grey diorite (Fig. 5a). The leucosomes are generally rimmed by thin, discontinuous melanosomes, referred to as concordant melanosome below. Leucosomes in the stromatic migmatites at Parkersville and Wren Lake range from concordant through slightly to strongly discordant to the host-rock foliation. Only concordant to slightly discordant leucosomes were sampled for geochemistry, and, for convenience, these leucosomes are referred to here as concordant leucosomes. The concordant leucosomes constitute 20–30 vol. %, and rarely up to 40–50 vol. %, of the stromatic migmatites, measured on a series of vertical cross-sections through the outcrops. The amount of leucosome varies strongly within an outcrop. At the Parkersville outcrop, for example, leucosome abundances range from <5 to c. 50 vol. % over a distance of a few tens of metres. We could not, however, identify any gradual variation between leucosome-poor and -rich parts of the outcrop. Locally, where leucosome proportions reach 40–50 vol. %, the stromatic migmatites contain rotated ‘blocks’ of migmatitic gneiss ‘floating’ in leucosome (Fig. 5c). The concordant leucosomes are variably folded and, in places, younger leucosomes can be seen cutting older, deformed leucosomes (Fig. 5a); the observed 20–50 vol. % leucosome is therefore cumulative, and is unlikely to represent the amount of melt present at any one time, nor does it represent the fraction of melt (F) generated by melting of the mesosome. Locally, the melts appear to have migrated along foliation planes and accumulated within shear bands or in boudin necks (Fig. 5b). These observations suggest that melting and/or melt movement in the rocks was coeval with deformation (e.g. Collins & Sawyer, 1996Go; Timmermann et al., 2002Go).



View larger version (91K):
[in this window]
[in a new window]
 
Fig. 5. (a) Concordant leucosome in stromatic migmatite. Arrows indicate locations where a discordant relationship is observed between different generations of leucosome; lens cap for scale c. 5 cm across. (b) Small leucosome patch in boudin neck; pen for scale c. 13 cm long. (c) Part of outcrop with a very large proportion of leucosome, containing small, rotated blocks of ‘normal’ stromatic migmatite; hammer for scale c. 30 cm long. (d) Discordant leucosome in small shear bands; field of view is c. 1 m across (vertical view).

 
The very high leucosome proportions observed in the stromatic migmatites (e.g. in the Parkersville outcrop) are inconsistent with in situ partial melting, both because granodioritic gneisses are unlikely to yield such high melt proportions and because complementary residual rocks are missing. Although the concordant melanosomes may represent residues after partial melting, the small volume of melanosome (typically only 1–2 mm thick) cannot account for the volume of leucosome in the rocks. A rigorous mass balance calculation is precluded by difficulties in estimating the composition of the concordant melanosome and the interpretation (below) that most of the leucosomes do not represent melt compositions.

The contrast in geometry between the patch and stromatic migmatites could result from differences in the timing of leucosome formation (i.e. patch migmatites formed post-kinematically) and/or from local differences in strain during deformation.Geochronological data from the stromatic migmatites indicate that the concordant leucosomes crystallized at c. 1065 Ma (Timmermann et al., 1997Go; Slagstad et al., 2004bGo), and that high-grade metamorphism and deformation persisted until c. 1050 Ma; thus it is possible that the patch migmatites formed later in the metamorphic history of the Muskoka domain and therefore experienced less deformation.The available field and geochronological data neither support nor reject this hypothesis. In contrast, there is good field evidence to suggest a significant rheological contrast between the two types of migmatite and we believe that this contrast may have been the main factor controlling the geometry of the leucosomes. Although the exact mechanism by which stromatic structures develop is obscure, work by Robin (1979)Go and Brown et al. (1995)Go suggests that they can develop in response to applied differential stress acting on a strongly anisotropic protolith, such as the banded, granodioritic host rock to the stromatic migmatites. In contrast, the dioritic host rock to the patch migmatites was stronger and more homogeneous, possibly resulting in a stress field that was closer to isotropic, favouring the observed patch geometry. At present, this is our preferred hypothesis for the geometrical differences between the patch and stromatic migmatites.


    PETROGRAPHY AND GEOCHEMISTRY
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING
 STUDY AREA
 MUSKOKA MIGMATITES
 PETROGRAPHY AND GEOCHEMISTRY
 MIGMATITE PETROGENESIS
 MELT MIGRATION
 CONCLUSIONS
 REFERENCES
 
Analytical procedures
Major and trace elements were determined on fused glass beads and pressed powder pellets, respectively, using standard X-ray fluorescence (XRF) techniques at St. Mary's University, Halifax, Nova Scotia. Rare earth elements (REE), Hf, and Th were determined by inductively coupled plasma-mass spectrometry (ICP-MS) using a Na2O2 sintering technique at Memorial University, St. John's, Newfoundland. Dostal et al. (1986)Go and Longerich et al. (1990)Go discussed the analytical procedures, accuracy, and precision of the XRF and ICP-MS analyses, respectively. The geochemical data are presented in Tables 1 and 2; modes of representative samples are presented in Table 3. Mineral abbreviations used in the text and figures are from Kretz (1983)Go.


View this table:
[in this window]
[in a new window]
 
Table 1: Geochemistry of the stromatic migmatites

 

View this table:
[in this window]
[in a new window]
 
Table 2: Geochemistry of the patch migmatites

 

View this table:
[in this window]
[in a new window]
 
Table 3: Modes of representative migmatite samples

 
Patch migmatite
The dioritic mesosome is fine grained, with a granoblastic to weakly foliated texture, dominantly comprising hornblende and plagioclase with minor biotite and quartz. Biotite laths, typically 1–3 mm long, define a weak foliation. Apatite and zircon are the main accessory minerals; allanite was not observed.

The patches are characterized by a mix of medium- to coarse-grained granitic leucosome (patch leucosome) with large (1–5 cm long), randomly oriented, euhedral to subhedral, poikiloblastic or skeletal hornblende (Figs. 4b, c and 6d). Plagioclase inclusions are abundant and typically constitute 30–40 vol. % of the poikiloblastic hornblendes, along with minor quartz and biotite inclusions. The plagioclase inclusions are typically irregular and many are oriented roughly along the cleavage of the hornblende. We interpret this texture to indicate that the hornblendes crystallized in the presence of melt, and their random orientation suggests minimal post-crystallization strain.



View larger version (121K):
[in this window]
[in a new window]
 
Fig. 6. (a) ‘Melt’ films along K-feldspar grain boundaries in granodioritic mesosome, from Parkersville outcrop. Similar ‘melt’ films have been observed in the dioritic mesosome. Scale bar represents 0·1 mm. (b) Coarse hornblende with rounded inclusions of plagioclase, biotite, and quartz from the concordant melanosome, from Parkersville outcrop. Coarse biotite on the right-hand side of the photograph is apparently in equilibrium with the hornblende. Scale bar represents 1 mm. (c) Small carbonate grain associated with biotite in granodioritic mesosome at Parkersville outcrop. Scale bar represents 0·5 mm. (d) Coarse, subhedral hornblende porphyroblast (dark) with abundant plagioclase inclusions in hornblende-rich patch, from Wren Lake outcrop. The hornblende is optically continuous. Scale bar represents 5 mm.

 
The patches grade into melanosome comprising medium-grained (~0·5 cm) hornblende (60–70 vol. %), plagioclase (30–40 vol. %), and minor biotite. This type of melanosome is referred to as patch melanosome to distinguish it from the concordant melanosome in the stromatic migmatites. Coarse hornblende in the patch melanosome typically contains rounded inclusions of plagioclase, biotite, and quartz. Biotite typically forms dark red to nearly opaque, corroded grains, in contrast to the fine biotite laths in the dioritic mesosome.

The dykes are pegmatitic but otherwise petrographically similar to the patch leucosomes. Along their margins, the dykes contain significant amounts of hornblende; some of this hornblende is texturally similar to the large, poikiloblastic hornblende in the patches, and some is similar to hornblende in the patch melanosome.

The dioritic mesosome and patch melanosome have overlapping major element compositions, but on average, the melanosome is lower in SiO2, Al2O3, and Na2O, and higher in TiO2, Fe2O3, MgO, and CaO (Fig. 7). The K2O content is similar in the mesosome and melanosome. The patch leucosomes and pegmatite dykes are enriched in SiO2 and K2O relative to the mesosome and melanosome, but depleted in the other major elements. Apart from patch leucosome sample M30072 [GenBank] -3-L2, the leucosomes and dykes have K2O >5·9 wt %, unusually high for granitic melts. This is illustrated in the normative Q–Ab–Or diagram (Fig. 7h), where the patch leucosomes and dykes are displaced towards the Or apex relative to cotectic melts. In contrast, sample M30072 [GenBank] -3-L2 plots closer to the eutectic composition for melts formed under water-undersaturated conditions corresponding to aH2O ~0·5 at c. 1 GPa (Johannes & Holtz, 1990Go). Sample M30072 [GenBank] -3-L2 is enriched in Al2O3, CaO, and Na2O, and depleted in K2O relative to the other patch leucosomes, indicating a higher plagioclase/K-feldspar ratio in this sample than in the other leucosomes.



View larger version (24K):
[in this window]
[in a new window]
 
Fig. 7. Harker variation diagrams and normative Q–Ab–Or ternary diagram, showing the compositions of patch migmatites. {blacksquare}, dioritic mesosome; {square}, patch leucosome; {triangleup}, dykes; x, patch melanosome.

 
Complete separation of melanosome from leucosome prior to analysis was generally achieved for the leucosome samples, whereas the sampled melanosome probably contains several percent leucosome. The effect of ‘contamination’ by 5–10 wt % leucosome is relatively minor for most major elements (<10–15%), but not for K2O, which increases by several tens of percent. The measured K2O content of the patch melanosome is, therefore, likely to be too high. The patch leucosomes are low in REE compared with the melanosomes (see below), thus ‘contamination’ by leucosome does not affect REE contents of the melanosome significantly.

The patch leucosomes and dykes are enriched in trace elements associated with feldspar (Ba, Rb, Sr) but depleted in other trace elements relative to the dioritic mesosome and patch melanosome. The patch leucosomes and dykes show a negative correlation between K2O and both Na2O and CaO, but a positive correlation between K2O and Rb/Sr (Fig. 8c), indicating that crystallization of feldspar was a major factor in the compositional variation of the leucosomes. The REE correlate with P2O5, Zr, Hf, Y, and Nb indicating that accessory apatite and zircon exerted a significant control on the behaviour of the REE. The (La/Yb)N ratios vary little within the patch migmatites, ranging from about five in the dioritic mesosome and patch melanosome to c. 10 in the patch leucosomes; the dykes have significantly higher (La/Yb)N ratios ranging from nine to 90. The Eu anomalies [(Eu/Eu*)N] of the patch leucosomes are inversely correlated with total REE (Fig. 8), increasing from 0·75 in the most REE-rich sample to 1·79 in the sample with lowest total REE. In the dykes, (Eu/Eu*)N ratios increase from 1·37 to 19·67 with decreasing total REE content. Despite significant major element differences between patch leucosome sample M30072 [GenBank] -3-L2 and the other patch leucosomes/dykes, the trace element contents are similar and fall along the same trends.



View larger version (39K):
[in this window]
[in a new window]
 
Fig. 8. Major and trace element bivariate diagrams for the patch leucosomes and dykes, and chondrite-normalized REE patterns for the patch migmatites comparing the different migmatite components. Normalized to values of Sun & McDonough (1989)Go. Symbols as in Fig. 7.

 
Stromatic migmatites
The concordant leucosomes and granodioritic mesosomes were sampled by serial cutting of large slabs. The concordant melanosomes were too thin to be sampled for geochemical analysis. Timmermann (1998)Go and Timmermann et al. (2002)Go described the petrography of the stromatic migmatites in the Muskoka domain in detail, and the following description therefore focuses on features pertinent to their petrogenesis.

The granodioritic to quartz monzodioritic mesosome (for simplicity referred to as granodioritic) is fine to medium grained and moderately foliated. Minor, interstitial ‘melt films’ (see Sawyer, 2001Go) along some grain boundaries (Fig. 6a) and rare, small (<1 cm) quartz–feldspar patches, interpreted to represent former melt patches, suggest some melting in the mesosome. However, there is no petrographic or geochemical evidence that significant volumes of melt were extracted from the mesosome. The proportion of melt films in the granodioritic mesosome, estimated by point counting as <0·5{per thousand}, is volumetrically insignificant. Apatite, zircon, and allanite are common accessory phases.

The concordant leucosomes are medium to coarse grained and granitic, with c. 1–5 vol. % hornblende porphyroblasts and flakes of biotite. The hornblende in the leucosomes is texturally similar to that found in the associated concordant melanosomes, and probably represents residual hornblende that was mechanically incorporated into the melt.

The concordant melanosomes are generally <5–10 mm thick and comprise hornblende, plagioclase, and biotite with accessory apatite, zircon, allanite, and titanite. Much of the hornblende in the concordant melanosome is texturally similar to hornblende in the granodioritic mesosome; however, a significant portion is coarser grained with rounded inclusions of plagioclase, quartz, and biotite (Fig. 6b); that is, similar to the patch melanosomes. Similar coarse hornblende is absent from the granodioritic mesosome. Coarse plagioclase with inclusions of biotite and quartz is common in the concordant melanosomes. Biotite in the concordant melanosomes is compositionally similar to biotite in the granodioritic mesosome, but is coarser and forms clusters intergrown and apparently in equilibrium with hornblende (Fig. 6b). However, the PT conditions experienced by these rocks (Timmermann et al., 2002Go; S. Gagné & J. Hawken, unpublished data, 2002) suggest that biotite was unlikely to have been stable during peak metamorphism (Fig. 3); it seems more likely that the biotite in the concordant melanosomes grew during retrograde reaction with the melt as the melt crystallized (see Waters, 2001Go).

The granodioritic mesosome ranges from 65 to 68 wt % SiO2. Samples from the two investigated outcrops at Parkersville and Wren Lake have similar major (Fig. 9) and trace element contents. The granodioritic mesosome from Wren Lake is slightly higher in total REE with more pronounced negative Eu anomalies [(Eu/Eu*)N = 0·40–0·45] than mesosome from Parkersville [(Eu/Eu*)N = 0·58], but the REE patterns are similarly fractionated [(La/Yb)N = 5–12] (Fig. 10).



View larger version (21K):
[in this window]
[in a new window]
 
Fig. 9. Harker variation diagrams and normative Q–Ab–Or ternary diagram, showing the compositions of stromatic migmatites. For granodioritic mesosome: •, Parkersville outcrop; {diamondsuit}, Wren Lake outcrop; for concordant leucosome: {circ}, Parkersville outcrop; {diamond}, Wren Lake outcrop.

 


View larger version (32K):
[in this window]
[in a new window]
 
Fig. 10. Chondrite-normalized REE patterns for the stromatic migmatites comparing the various migmatite components, and bivariate trace element diagrams. Normalized to values of Sun & McDonough (1989)Go. Symbols as in Fig. 9.

 
The concordant leucosomes contain 70–73 wt % SiO2; although clustered, most major elements decrease with increasing SiO2 (Fig. 9) whereas K2O and most trace elements show little systematic variation. As for the patch leucosomes, K2O is negatively correlated with CaO and Na2O, and positively correlated with Rb/Sr, suggesting crystallization of feldspar. In the normative Q–Ab–Or diagram (Fig. 9h) the concordant leucosomes plot close to the cotectic line separating the primary fields for quartz and K-feldspar, similar to the patch leucosomes. The REE patterns (Fig. 10b) are moderately fractionated, with (La/Yb)N increasing from eight to 24 with increasing SiO2 (Fig. 10c). Total REE and Eu anomalies [(Eu/Eu*)N = 0·7–1·3] are negatively correlated (Fig. 10d), but neither are correlated with SiO2.

Relative to the granodioritic mesosome, the concordant leucosomes are depleted in all major elements except SiO2 and K2O. Barium, Rb, and Sr are enriched in the leucosomes. Yttrium, Zr, Hf, and Nb are lower in the leucosomes, and Th is similar in leucosomes and granodioritic mesosomes. Total REE contents are lower, REE patterns more fractionated, and Eu anomalies are more positive in the leucosomes than in the mesosomes (Fig. 10a and b).


    MIGMATITE PETROGENESIS
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING
 STUDY AREA
 MUSKOKA MIGMATITES
 PETROGRAPHY AND GEOCHEMISTRY
 MIGMATITE PETROGENESIS
 MELT MIGRATION
 CONCLUSIONS
 REFERENCES
 
Evidence for fluid-present partial melting
Feldspars, quartz, hornblende, and biotite in variable proportions dominate the mineralogy of the mesosome, melanosome, and leucosome in both types of migmatite. Despite the relatively simple mineralogy, the textures described above provide some clues to the melting reactions in these rocks. In both the patch and concordant melanosome, coarse-grained hornblende with numerous rounded inclusions of plagioclase, biotite, and quartz suggest that melting took place by the reaction

(1)
Our interpreted melting reaction is comparable with that suggested by Timmermann et al. (2002)Go, who proposed fluid-absent biotite-dehydration melting according to the reaction

(2)
Our reasons for invoking fluid-present partial melting, and the constraints on fluid composition, are explained below. In the following discussion, the term ‘fluid’ is used to describe a volatile phase (e.g. H2O), whereas ‘melt’ is used to describe silicate melt ± entrained crystals.

Petrographic evidence suggests that partial melting in both stromatic and patch migmatite was associated with growth of new hornblende. Experimental studies (e.g. Büsch et al., 1974Go; Naney, 1983Go; Gardien et al., 2000Go) have shown that crystallization and growth of hornblende during partial melting of a biotite–plagioclase–quartz assemblage requires at least 3–4 wt % H2O in the melt. New hornblende can form by biotite-dehydration melting under fluid-absent conditions (e.g. Skjerlie & Johnston, 1996Go); however, in this case, the breakdown of biotite to form hornblende produces <5 vol. % melt at 850°C. Larger volumes of melt are produced at higher temperatures (>900°C), where hornblende reacts to form melt and residual pyroxene (e.g. Wolf & Wyllie, 1994Go); however, this reaction is not consistent with observations from the Muskoka migmatites.

If the granodioritic and dioritic mesosomes in the Muskoka migmatites represent the protoliths of the concordant and patch leucosomes, respectively, the relatively low abundance of biotite (<10 vol. %) is insufficient to produce significant melt volumes. For example, given 10 vol. % biotite with 3·9 wt % H2O, and reasonable estimates of density (3·0 for biotite; 2·7 for source rock), up to 0·4 wt % H2O is released by biotite breakdown. According to Clemens & Vielzeuf (1987)Go, this amount of H2O can produce c. 10 vol. % melt in a quartzofeldspathic source rock at 1 GPa. This is an upper estimate because the calculation assumes that all the fluid liberated by dehydration enters the melt (i.e. the residue is anhydrous); in contrast, partial melting in the Muskoka domain involved growth of a new hydrous phase (hornblende) that would have reduced the amount of fluid available to the melt. The estimated 10% melting is significantly less than the degree of partial melting estimated from the geochemical data (discussed below), thus fluid-absent, biotite-dehydration melting appears unlikely to have produced large amounts of melt in the Muskoka domain. Therefore, although small amounts of melt are likely to have formed by dehydration melting of biotite at T >800°C, we conclude that partial melting dominantly took place in the presence of an H2O-bearing fluid [reactions (3b) and (6c) in Fig. 3]. A similar model was invoked by McMullen (1999)Go for the eastern Muskoka domain, and has been proposed to account for partial melting of similar lithologies elsewhere (e.g. Lappin & Hollister, 1980Go; Kenah & Hollister, 1983Go; McLellan, 1988Go; Mogk, 1992Go; Nédélec et al., 1993Go; Sawyer, 1998Go; Escuder Viruete, 1999Go).

The composition of the fluid is impossible to constrain without isotopic and/or fluid inclusion data. Phase equilibrium evidence from high-grade metamorphic rocks suggests that pure aqueous fluids are generally not present in excess [see Stevens & Clemens (1993)Go and references therein], thus, aH2O was probably less than unity during migmatization. This is also consistent with the lack of leucosome in metabasite enclaves in the area [reaction (9) in Fig. 3]. Values of aH2O <1·0 can be achieved in systems that contain insufficient H2O-rich fluid to saturate the melt, or that contain a mixed-volatile fluid (e.g. H2O + CO2) (Clemens & Watkins, 2001Go). In the Muskoka domain, several observations suggest the presence of a mixed H2O + CO2 fluid during melting. The migmatites locally contain small amounts of carbonate associated with skeletal biotite and biotite–hornblende intergrowths that are interpreted to have been involved in melting (Fig. 6c). Diffuse granulite patches, also containing small amounts of carbonate (T. Slagstad, unpublished data, 2002), are common in parts of the Muskoka domain, and were interpreted by Timmermann et al. (2002)Go to have formed by influx of low-aH2O, CO2-rich fluids coeval with formation of the migmatites. Infiltration of CO2-rich fluid was also invoked by Pattison (1991)Go to explain the development of light-coloured ‘net veins’ rich in pyroxene, plagioclase, and locally quartz, in metagabbro enclaves in the adjacent Seguin subdomain (Fig. 2). Garnet porphyroblasts in pelitic gneisses in the underlying Algonquin domain (Fig. 2; north of Dorset) contain CO2-rich fluid inclusions interpreted to have formed as the result of partitioning of H2O into coexisting melt (Layman, 2003Go).

We speculate that influx of a mixed H2O–CO2 fluid during high-grade Grenvillian metamorphism can account for both extensive partial melting (e.g. Giorgetti et al., 1996Go) and the local development of granulite-facies assemblages. According to this model, H2O would have been partitioned into the melts, whereas most of the CO2 would have remained in the rocks (see Clemens, 1993Go), forming the patch granulites and small accumulations of carbonate. During cooling and melt crystallization, back-reaction between H2O and orthopyroxene could have formed the hornblende-rimmed orthopyroxenes typical of the patch granulites (Timmermann et al., 2002Go). Detailed isotopic and fluid inclusion analyses from both the migmatites and the patch granulites are needed to determine the role and nature of fluids during high-grade metamorphism in these rocks.

Petrogenetic constraints from field and petrographic data
The field and petrographic observations from the patch migmatites provide important clues to their petrogenesis. The patch leucosomes are irregular, randomly distributed, and have associated melanosomes, consistent with in situ partial melting of dioritic mesosome (see McLellan, 1988Go). The partial melting reaction apparently produced a solid residue comprising 60–70 vol. % hornblende and 30–40 vol. % plagioclase as indicated by the petrographic observations. In addition, both patch and concordant melanosomes contain a significant proportion of biotite that we interpret to have formed during melt crystallization.

Thin seams of granitic material typically connect the patches and dykes and could represent former melt conduits along which granitic melts migrated from the patches. If this is true, the patch leucosomes may represent feldspar-rich cumulates, probably with some trapped, evolved melt. Melt migration probably took place by porous flow over distances of a few centimetres from melanosome to patch, and by channelled flow over distances of up to a few metres from patch to dyke. The zoning in the dykes, with feldspar-rich rims and quartz-rich centres, is similar to that observed in many granite–pegmatite systems (e.g. Jahns & Burnham, 1969Go; Simmons et al., 1987Go), suggesting that feldspar crystallized from the melts. The pegmatitic texture indicates that the melts were water rich, and the network of granitic seams and pegmatite dykes could indicate that the magma and/or fluid pressure was great enough to cause fracturing of the host rock (e.g. Clemens & Mawer, 1992Go). Fluid-present melting is associated with a negative volume change, thus hydraulic fracturing is unlikely to have taken place during partial melting of the rocks. A more likely hypothesis, based on evidence for crystallization of hydrous melts (discussed further below), is that the magma reached water saturation, resulting in the formation of a free fluid phase, and thus a positive volume change, which caused fracturing of the surrounding rocks. Clear evidence of hydraulic fracturing (see Davidson et al., 1994Go) is, however, missing from these rocks, and the process is considered only in passing in the petrogenetic model outlined below.

The high strain in the stromatic migmatites obscures many field and petrographic clues to their petrogenesis. Locally, however, the stromatic migmatites preserve features suggesting a comparable evolution to the patch migmatites. Irregularly distributed, thin melanosomes are common, suggesting that partial melting produced a hornblende + plagioclase-dominated residue like that associated with the patch migmatites. Coarse hornblende crystals, similar to the poikiloblastic hornblende in the patch migmatites, are locally present. Cross-cutting leucosome dykes are relatively common (Fig. 5d) and may have formed by hydraulic fracturing during melt crystallization, as proposed for the patch migmatites.

Petrogenetic constraints of patch migmatites from trace element (REE) modelling
Partial melting
Trace element concentrations of model melts were calculated assuming equilibrium batch melting, in which the melt continually reacts and re-equilibrates with the solid residue at the site of melting until it is able to escape as a single ‘batch’. The concentration of trace element i in the model melt () is given by the equation

(3)
(Shaw, 1970Go), where ; is the bulk partition coefficient of element i in the residual solid, is the partition coefficient of element i between melt and mineral a, XRS is the weight proportion of mineral a in the residual solid, is the concentration of element i in the original solid, and F is the degree of partial melting. The concentration of element i in the complementary residue () is obtained by mass balance, such that

(4)

Crystallization
The trace element distribution during melt crystallization is modelled assuming equilibrium crystallization, meaning that the crystallized solids (cumulate phases) remain in equilibrium with the melt throughout its crystallization history. We prefer the equilibrium crystallization model over the Rayleigh fractional crystallization model, in which only the crystal surfaces are in equilibrium with the melt, because a wealth of metamorphic, textural, and geochronological data from the Central Gneiss Belt (Culshaw et al., 1997Go; Wodicka et al., 2000Go) suggest protracted high-grade metamorphic conditions and slow regional cooling, which would have facilitated equilibration between melt and cumulate phases (e.g. Sawyer, 1987Go).

The concentration of trace element i in the evolved melt () is given by the equation

(5)
(Allègre & Minster, 1978Go; Hanson, 1978Go), where ; is the bulk partition coefficient of element i in the phases crystallized, is the partition coefficient of element i between melt and mineral a, XCP is the weight proportion of mineral a crystallized, is the concentration of element i in the primary melt, and F is the fraction of melt remaining. The concentration of element i in the complementary cumulate () is obtained by mass balance, such that

(6)
The trace element modelling is performed in two stages. Stage 1 models partial melting of the dioritic mesosome leaving a hornblende + plagioclase-dominated residue, producing a ‘primary’ melt (20–40% melting). Stage 2 models feldspar-dominated equilibrium crystallization in the patches. Partition coefficients and sources are given in Table 4.


View this table:
[in this window]
[in a new window]
 
Table 4: Partition coefficients used in partial melting modelling

 
As is nearly always the case with trace element modelling, the models are not unique, and several alternative models involving different phase proportions are possible. Here, we use the petrographic data to constrain the phase proportions in the models. This approach effectively tests whether the petrogenetic model devised from the petrographic data is consistent with the observed geochemistry of the rocks. Another well-known problem in applying trace element modelling to granitic (sensu lato) melts is the common presence of accessory phases such as apatite, zircon, allanite, and monazite. These phases are extremely rich in certain trace elements [e.g. light REE (LREE) in allanite and monazite, heavy REE (HREE) in zircon] and in many cases can obscure the effects produced by the major phases during partial melting and crystallization (e.g. Bea, 1996Go, and references therein). Zircon, apatite, and allanite are common accessory phases in the rocks studied here, particularly in the melanosomes, and therefore should be included in the models. Apatite has only moderate partition coefficients and does not differentiate the REE significantly, whereas zircon and allanite strongly fractionate the HREE and LREE, respectively. The partition coefficient for Eu in allanite is an order of magnitude smaller than that of the other middle REE (MREE) (Table 4); thus, fractionation of allanite influences the Eu anomaly in addition to the LREE. The effects of the accessory phases on the model results presented here are discussed in detail in the next section.

Stage 1: equilibrium melting of dioritic mesosome
The first stage of the model involves 20–40% equilibrium batch melting of the average dioritic mesosome. The model restite is hornblende dominated, as indicated by the petrographic data, comprising 63 wt % hornblende and 36 wt % plagioclase, with trace amounts of zircon, apatite, and allanite (see Fig. 11a and caption for details). The model melts are depleted in REE relative to the source and have negligible Eu anomalies ((Eu/Eu*)N = 0·91 and 0·87 for 20 and 40% partial melting, respectively), resembling the REE pattern for sample M30072 [GenBank] -3-L2 with a (Eu/Eu*)N of 0·89.



View larger version (34K):
[in this window]
[in a new window]
 
Fig. 11. Results of REE modelling on the patch migmatites using the equations presented in the text and partition coefficients from Table 4. (a) Partial melting of average dioritic mesosome (source) leaving a hornblende + plagioclase-dominated residue. Residue comprises (wt %) hornblende (63), plagioclase (36), zircon (0·0078), apatite (0·0053), and allanite (0·0027). Model melts and residues are calculated for 20 and 40% partial melting (i.e. F = 0·2 and 0·4, respectively). Symbols as in Fig. 7. (b) Model melts formed by 40% partial melting, similar to (a), but with varying proportions of accessory phases. R, reference line ‘Model melt F = 0·4’ in (a); 1, model melt with no accessory phases in the residue; 2, model melt with no apatite and –25% zircon (0·0059 wt %) and allanite (0·0020 wt %) relative to reference line (R); 3, model melt with 2 wt % apatite and +25% zircon (0·0096 wt %) and allanite (0·0034 wt %) relative to reference line (R). Reducing or increasing the proportions of allanite and zircon by more than 25% relative to the reference line yields model REE patterns that diverge significantly from those observed in the patch leucosomes. (c) Crystallization of a ‘primary’ melt produced by 40% partial melting in (a). Crystallizing phases are (in wt %) K-feldspar (60), plagioclase (20), and quartz (20). The modelling shows that the composition of the patch leucosomes and dykes is compatible with a cumulate, produced by 40% crystallization, mixed with various proportions (0–50 wt %) of evolved melt.

 
Apatite has a negligible impact on the model result; varying the amount of apatite in the residue from zero to 0·02 wt % produces a 15% variation in the total REE and <10% variation in the (La/Yb)N ratio. Apatite was, nevertheless, included because the petrographic data suggest that it is a common accessory phase. The geochemical data from the single patch melanosome analysis indicate c. 0·001 wt % apatite in the residue. Zircon effectively controls the amount of Tm, Yb, and Lu, and to a lesser extent Ho and Er, in the model melt; allanite controls the amount of the LREE and the size of the Eu anomaly. The proportion of zircon in the residue is estimated from the geochemical data to be c. 0·015 wt %. However, because zircon typically occurs as inclusions in refractory phases shielding it from the melt (e.g. Watson et al., 1989Go) or because of limited solubility in the melt (e.g. Watson & Harrison, 1983Go), this is a maximum estimate of how much zircon participated in the melting reaction (e.g. Watson & Harrison, 1983Go; e.g. Watson et al., 1989Go). Here, the proportions of zircon and allanite in the model residue were selected such that the REE pattern of the model melt resembled that observed for the patch leucosomes; that is, moderately negatively sloping LREE pattern and flat HREE pattern. The weight fraction of zircon participating in the melting reaction (0·0078 wt %) is about half that estimated from the geochemical data, suggesting that about half of the zircon present in the residue (patch melanosome) did not participate in the melting reaction. The amount of Zr in the leucosomes is relatively high, and calculations by Slagstad (2003)Go suggest that the melts were saturated in Zr. Thus, it is likely that limited solubility of Zr, possibly in addition to shielding by refractory phases, inhibited further dissolution of zircon. If the weight proportions of zircon and allanite participating in the melting reaction vary by >25% relative to those used in the partial melting model, the resulting model REE patterns (but not REE abundances) deviate significantly from the observed REE patterns of the patch leucosomes (Fig. 11b).

The REE contents of the model restite are similar to those of the patch melanosome, consistent with the latter representing the residue after melt extraction. The model melt produced by 40% partial melting was used as the ‘primary’ melt in stage 2.

Stage 2: crystallization in the patches
Stage 2 models equilibrium crystallization of the ‘primary’ melt produced in stage 1, crystallizing 60% K-feldspar, 20% plagioclase, and 20% quartz. The petrographic data suggest that K-feldspar is the main fractionating phase, but the K-feldspar/plagioclase ratio is unconstrained. Because both feldspars have small and essentially similar REE partition coefficients, changing the proportion of K-feldspar to plagioclase has a negligible effect on the model result. The interpretation that K-feldspar was a crystallizing phase implies that the partial melt generated in Stage 1 was relatively potassic. This is difficult to account for in the partial melting reaction proposed above, and it is possible that the infiltrating fluid introduced potassium to the melt system, as observed in other migmatites (e.g. Olsen, 1984Go). In contrast, the solubility of REE in water-rich fluids is generally considered low, thus the REE contents of the patch migmatites were probably little affected by the infiltrating fluid. The results of the crystallization modelling suggest that most of the patch leucosomes and dykes are feldspar cumulates formed by extreme degrees of crystallization (95–99%). An alternative, and probably more realistic, interpretation is that the leucosomes and dykes represent cumulates, formed by lower degrees of crystallization than above, with some trapped evolved melt. Figure 11c shows the REE pattern obtained by mixing a cumulate derived by 40% crystallization with a varying weight fraction (0–50 wt %) of trapped, evolved melt. In this case, the majority of patch leucosomes and several of the dyke samples can be interpreted as cumulates mixed with 20–50 wt % evolved melt.

Petrogenesis of the stromatic migmatites
Despite the high strain, the stromatic migmatites preserve petrographic and geochemical features suggesting that they evolved in much the same way as the less deformed patch migmatites. In particular, the concordant leucosomes are rimmed by thin, discontinuous melanosomes with textures resembling those observed in the patch melanosomes. The observations are consistent with fluid-present melting of plagioclase, biotite, and quartz, leaving a hornblende + plagioclase-dominated residue, and producing a granitic melt. The REE patterns of the concordant leucosomes display a more positive Eu anomaly with decreasing total REE, similar to the patch leucosomes, indicating that feldspar-dominated crystallization was important in their evolution. In general, the concordant leucosomes have less positive Eu anomalies than the patch leucosomes. A likely explanation for this difference is that the granodioritic mesosome, a likely source composition for some of the concordant leucosomes, has a more negative Eu anomaly than the dioritic mesosome, which would result in a ‘primary’ melt with a more negative Eu anomaly, assuming a similar melting reaction. Thus, despite their less positive Eu anomalies, we interpret the concordant leucosomes to represent crystal-rich melts or cumulates with some trapped, evolved melt; that is, similar to the patch leucosomes. This interpretation suggests that melts were continuously drained from the rocks, and may explain why the rocks did not develop into diatexites despite having present-day integrated leucosome proportions of up to 40–50 vol. % (see also Guernina & Sawyer, 2003Go).

The field, petrographic, and geochemical evidence from the stromatic migmatites suggests that in addition to in situ partial melting, a significant proportion of the leucosome present in the investigated outcrops was derived from external sources. There is no petrographic or geochemical evidence that the granodioritic mesosomes reported here (Table 1), or leucosome-poor grey gneisses throughout the Muskoka domain (Slagstad et al., 2004aGo), are residual after large-scale melt extraction. This suggests that any externally derived melts must have migrated from structurally deeper, or more interior, parts of the orogen than those now exposed in the Muskoka domain, implying minimum melt transport distances of several kilometres.


    MELT MIGRATION
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING
 STUDY AREA
 MUSKOKA MIGMATITES
 PETROGRAPHY AND GEOCHEMISTRY
 MIGMATITE PETROGENESIS
 MELT MIGRATION
 CONCLUSIONS
 REFERENCES
 
Several mechanisms for melt transfer in the crust have been proposed, including diapirism (Bateman, 1984Go), dyking (Petford et al., 1993aGo, 1993bGo), fracture propagation (Clemens & Mawer, 1992Go; Davidson et al., 1994Go), and migration through faults or shear zones (Hollister & Crawford, 1986Go; Brown, 1994Go). In a study of the Mount Hay region of Australia, Collins & Sawyer (1996)Go proposed that melt migrated through a distributed, structurally controlled network of conduits developed during compressional, non-coaxial deformation. The driving force for melt transfer was interpreted to be fluid-pressure-driven dilatancy resulting from hydraulic fracturing, which formed a pervasive interconnected network through which the melts flowed. Although the geometry of the dykes in the patch migmatites is consistent with hydraulic fracturing, evidence of a similar process in the stromatic migmatites is less obvious. In these rocks fracturing may have been parallel to the foliation, controlled by the magnitude of difference in tensile strength normal and parallel to the schistosity, the orientation of the principal stress axes relative to the schistosity, and the differential stress applied to the rocks (see Wickham, 1987Go). However, local discordant, leucosome-filled veins in the stromatic migmatites (Fig. 5d) indicate that some steep fractures formed, possibly in structurally favourable sites or late in the deformation history.

Field observations from the stromatic migmatites suggest that melt was present during deformation, consistent with geochronological data (Timmermann et al., 1997Go; Slagstad et al., 2004bGo) showing that leucosomes formed at c. 1070–1050 Ma, during the main stage of the Ottawan orogeny. Based on the available data, we conclude that melt migration was driven by differential pressure associated with syn-melt deformation on the local to regional scale (see Sawyer, 1991Go, 1994Go; Brown, 1994Go; Collins & Sawyer, 1996Go). Melt transfer may have been dominantly vertical (towards the upper crust), with originally steep melt-transfer networks subsequently deformed into their present orientations, or melts may have migrated laterally along the sub-horizontal leucosome networks observed today. The prevalence of structures with low to moderate dips, even in relatively low-strain regions, is more consistent with the latter possibility. This implies lateral melt migration, and thus lateral pressure gradients, on a regional scale. If so, the Muskoka domain may represent the remnants of a regional-scale, melt-transfer zone developed at the mid-crustal level of an orogen large enough to have developed a significant differential pressure between its core and its flanks.

Despite indications that a large volume of melt migrated through the Muskoka domain, no map-scale Ottawan granite plutons that might correspond to melt collection zones have been identified in or near the study area. The absence of upper-crustal rocks in the Central Gneiss Belt means that the corresponding plutons, if they existed, may have been removed by erosion. This possibility is suggested by the abundance of Ottawan-age detrital zircons in Neoproterozoic sedimentary rocks from western North America (e.g. Ross et al., 1992Go; Rainbird et al., 1997Go) and in modern Appalachian rivers (Eriksson et al., 2003Go).


    CONCLUSIONS
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING
 STUDY AREA
 MUSKOKA MIGMATITES
 PETROGRAPHY AND GEOCHEMISTRY
 MIGMATITE PETROGENESIS
 MELT MIGRATION
 CONCLUSIONS
 REFERENCES
 
(1) We recognize two distinct types of migmatite in the Muskoka domain, patch and stromatic migmatites, based on leucosome geometry and host-rock characteristics.

(2) The present-day morphology of migmatites in the Muskoka domain appears to be controlled by differences in strain. The geometry of leucosomes in the patch migmatites and the observation that deformed stromatic migmatites wrap around smaller, less deformed bodies of patch migmatite suggest that the dioritic rocks hosting the patch migmatites were stronger than the surrounding granodioritic orthogneisses, shielding the leucosomes from the effects of deformation.

(3) Both types of migmatite probably formed by fluid-present melting of plagioclase, biotite, and quartz, producing a hornblende + plagioclase-dominated residue and granitic melt. Crystallization of feldspar and quartz produced granitic cumulates, represented by most of the studied leucosomes. The fate of the complementary evolved melts could not be determined.

(4) Field and geochemical data from the stromatic migmatites suggest significant melt contributions from external sources, and indicate that the Muskoka domain may have acted as a mid-crustal melt transfer zone during Grenvillian (Ottawan) contractional deformation.


    ACKNOWLEDGEMENTS
 
Barrie Clarke and Ed Sawyer provided very helpful comments on an earlier version of the paper, resulting in significant improvements. Pat Bogutyn provided able assistance in the field, and together with Kat Eisnor did much of the sample preparation. Gordon Brown made the thin sections. Pierre Barbey, Ed Sawyer, and Serge Fourcade are thanked for insightful reviews that improved and clarified the paper; Nick Arndt provided helpful editorial comments. T.S. was supported by a Dalhousie University Killam Scholarship and Ph.D. grant NFR-138552/432 from the Norwegian Research Council (NFR). Field and analytical work was funded by NSERC Discovery Grants to R.A.J. and N.G.C. Writing of this paper was carried out while T.S. was a post-doctoral fellow at the Geological Survey of Norway (NGU).


* Corresponding author. Present address: Geological Survey of Norway (NGU), Leiv Erikssons v. 39, 7491 Trondheim, Norway. Telephone: +47 73 90 42 29. Fax: +47 73 92 16 20. E-mail: Trond.Slagstad{at}ngu.no


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING
 STUDY AREA
 MUSKOKA MIGMATITES
 PETROGRAPHY AND GEOCHEMISTRY
 MIGMATITE PETROGENESIS
 MELT MIGRATION
 CONCLUSIONS
 REFERENCES
 
Allègre, C. J. & Minster, J. F. (1978). Quantitative models of trace element behavior in magmatic processes. Earth and Planetary Science Letters 38, 1–25.[CrossRef][Web of Science]

Arth, J. G. (1976). Behaviour of trace elements during magmatic processes—a summary of theoretical models and their applications. Journal of Research of the US Geological Survey 4, 41–47.[Web of Science]

Arth, J. G. & Barker, F. (1976). Rare earth partitioning between hornblende and dacitic liquid and implications for the genesis of trondhjemitic–tonalitic magmas. Geology 4, 534–536.[Abstract/Free Full Text]

Arth, J. G. & Hanson, G. N. (1975). Geochemistry and origin of the early Precambrian crust of northeastern Minnesota. Geochimica et Cosmochimica Acta 39, 325–362.[CrossRef][Web of Science]

Ashworth, J. R. (1985). Introduction. In: Ashworth, J. R. (ed.) Migmatites. Glasgow: Blackie, pp. 1–35.

Barbey, P., Bertrand, J.-M., Angoua, S. & Dautel, D. (1989). Petrology and U/Pb geochronology of the Telohat migmatites, Aleksod, Central Hoggar, Algeria. Contributions to Mineralogy and Petrology 101, 207–219.[CrossRef][Web of Science]

Barraud, J., Gardien, V., Allemand, P. & Grandjean, P. (2004). Analogue models of melt-flow networks in folding migmatites. Journal of Structural Geology 26, 307–324.[CrossRef][Web of Science]

Bateman, R. (1984). On the role of diapirism in the segregation, ascent and final emplacement of granitoid magmas. Tectonophysics 110, 211–231.[CrossRef][Web of Science]

Bea, F. (1996). Residence of REE, Y, Th, and U in granites and crustal protoliths: implications for the chemistry of crustal melts. Journal of Petrology 37, 521–552.[Abstract/Free Full Text]

Beaumont, C., Jamieson, R. A., Nguyen, M. H. & Lee, B. (2001). Himalayan tectonics explained by extrusion of a low-viscosity crustal channel coupled to focused surface denudation. Nature 414, 738–742.[CrossRef][Medline]

Beaumont, C., Jamieson, R. A., Nguyen, M. H. & Medvedev, S. (2004). Crustal channel flows: 1. Numerical models with application to the tectonics of the Himalayan–Tibetan orogen. Journal of Geophysical Research 109, B06406, doi: 10.1029/2003JB002809.[CrossRef]

Brown, M. (1973). The definition of metatexis, diatexis and migmatite. Proceedings of the Geologists' Association 84, 371–382.

Brown, M. (1979). The petrogenesis of the St. Malo migmatite belt, Armorican Massif, France, with particular reference to the diatexites. Neues Jahrbuch für Mineralogie, Abhandlungen 135, 48–74.

Brown, M. (1994). The generation, segregation, ascent and emplacement of granite magma: the migmatite-to-crustally-derived granite connection in thickened orogens. Earth-Science Reviews 36, 83–130.

Brown, M. & Solar, G. S. (1998). Granite ascent and emplacement during contractional deformation in convergent orogens. Journal of Structural Geology 20, 1365–1393.[CrossRef][Web of Science]

Brown, M., Averkin, Y. A., McLellan, E. L. & Sawyer, E. W. (1995). Melt segregation in migmatites. Journal of Geophysical Research 100, 15655–15679.[CrossRef][Web of Science]

Burr, J. L. & Carr, S. D. (1994). Structural geometry and U–Pb geochronology near Lithoprobe seismic Line 32, Western Central Metasedimentary Belt, Grenville Province, Ontario: new results from high precision U–Pb geochronology. In: Paper presented at Lithoprobe Abitibi–Grenville Project, Workshop Report, pp. 59–62.

Büsch, W., Schneider, G. & Mehnert, K. R. (1974). Initial melting at grain boundaries. Part II: melting in rocks of granodioritic, quartz dioritic, and tonalitic composition. Neues Jahrbuch für Mineralogie, Monatshefte 8, 345–370.

Carr, S. D., Easton, R. M., Jamieson, R. A. & Culshaw, N. G. (2000). Geologic transect across the Grenville orogen of Ontario and New York. Canadian Journal of Earth Sciences 37, 193–216.[CrossRef][Web of Science]

Clemens, J. D. (1993). Experimental evidence against CO2-promoted deep crustal melting. Nature 363, 336–338.[CrossRef]

Clemens, J. D. & Mawer, C. K. (1992). Granitic magma transport by fracture propagation. Tectonophysics 204, 339–360.[CrossRef][Web of Science]

Clemens, J. D. & Vielzeuf, D. (1987). Constraints on melting and magma production in the crust. Earth and Planetary Science Letters 86, 287–306.[CrossRef][Web of Science]

Clemens, J. D. & Watkins, J. M. (2001). The fluid regime of high-temperature metamorphism during granitoid magma genesis. Contributions to Mineralogy and Petrology 140, 600–606.[Web of Science]

Collins, W. J. & Sawyer, E. W. (1996). Pervasive granitoid magma transfer through the lower–middle crust during non-coaxial compressional deformation. Journal of Metamorphic Geology 14, 565–579.[CrossRef][Web of Science]

Conrad, W. K., Nicholls, I. A. & Wall, V. J. (1988). Water-saturated and -undersaturated melting of metaluminous and peraluminous crustal compositions at 10 kb: evidence for the origin of silicic magmas in the Taupo Volcanic Zone, New Zealand, and other occurrences. Journal of Petrology 29, 765–803.[Abstract/Free Full Text]

Culshaw, N. G. & Dostal, J. (1997). Sand Bay gneiss association, Grenville Province, Ontario: a Grenvillian rift- (and -drift) assemblage stranded in the Central Gneiss Belt? Precambrian Research 85, 97–113.[CrossRef][Web of Science]

Culshaw, N. G. & Dostal, J. (2002). Amphibolites of the Shawanaga domain, Central Gneiss Belt, Grenville Province, Ontario: tectonic setting and implications for relations between the Central Gneiss Belt and Midcontinental USA. Precambrian Research 113, 65–85.[CrossRef][Web of Science]

Culshaw, N. G., Davidson, A. & Nadeau, L. (1983). Structural subdivisions of the Grenville Province in the Parry Sound–Algonquin region, Ontario. Current Research, Part B, Geological Survey of Canada, Paper 83-1B, 243–252.

Culshaw, N. G., Jamieson, R. A., Ketchum, J. W. F., Wodicka, N., Corrigan, D. & Reynolds, P. H. (1997). Transect across the northwestern Grenville orogen, Georgian Bay, Ontario: polystage convergence and extension in the lower orogenic crust. Tectonics 16, 966–982.[CrossRef][Web of Science]

Davidson, A. (1984). Identification of ductile shear zones in the southwestern Grenville Province of the Canadian Shield. In: Kröner, A. & Greiling, E. (eds), Precambrian Tectonics Illustrated. Stuttgart: Schweizerbart, pp. 263–279.

Davidson, A. & van Breemen, O. (1988). Baddeleyite–zircon relationships in coronitic metagabbro, Grenville Province, Ontario: implications for geochronology. Contributions to Mineralogy and Petrology 100, 291–299.[CrossRef][Web of Science]

Davidson, C., Schmid, S. M. & Hollister, L. S. (1994). Role of melt during deformation in the deep crust. Terra Nova 6, 133–142.[Web of Science]

Dostal, J., Baragar, W. R. A. & Dupuy, C. (1986). Petrogenesis of the Natkusiak continental basalts, Victoria Island, Northwest Territories, Canada. Canadian Journal of Earth Sciences 23, 622–632.[Web of Science]

Eriksson, K. A., Campbell, I. H., Palin, M. & Allen, C. M. (2003). Predominance of Grenvillian magmatism recorded in detrital zircons from modern Appalachian rivers. Journal of Geology 111, 707–717.[CrossRef][Web of Science]

Escuder Viruete, J. (1999). Hornblende-bearing leucosome development during syn-orogenic crustal extension in the Tormes Gneiss Dome, NW Iberian Massif, Spain. Lithos 46, 751–772.[CrossRef][Web of Science]

Fyfe, W. S. (1973). The granulite facies, partial melting and the Archean crust. Philosophical Transactions of the Royal Society of London, Series A 273, 457–461.[CrossRef]

Gardien, V., Thompson, A. B. & Ulmer, P. (2000). Melting of biotite + plagioclase + quartz gneisses: the role of H2O in the stability of amphibole. Journal of Petrology 41, 651–666.[Abstract/Free Full Text]

Giorgetti, G., Frezzotti, M.-L.-E., Palmeri, R. & Burke, E. A. J. (1996). Role of fluids in migmatites: CO2–H2O fluid inclusions in leucosomes from the Deep Freeze Range migmatites (Terra Nova Bay, Antarctica). Journal of Metamorphic Geology 14, 307–317.[CrossRef][Web of Science]

Green, A. G., Milkereit, B., Davidson, A., Spencer, C., Hutchinson, W. F., Cannon, W. F., Lee, M. W., Agena, W. F., Behrendt, J. C. & Hinze, W. J. (1988). Crustal structure of the Grenville Front and adjacent terranes. Geology 16, 788–792.[Abstract/Free Full Text]

Green, T. H. (1982). Anatexis of mafic crust and high pressure crystallization of andesite. In: Thorpe, R. S. (ed.) Andesites. London: John Wiley, pp. 465–487.

Guernina, S. & Sawyer, E. W. (2003). Large-scale melt-depletion in granulite terranes: an example from the Archean Ashuanipi Subprovince of Quebec. Journal of Metamorphic Geology 21, 181–201.[CrossRef][Web of Science]

Haggart, M. J., Jamieson, R. A., Reynolds, P. H., Krogh, T. E., Beaumont, C. & Culshaw, N. G. (1993). Last gasp of the Grenville orogeny—thermochronology of the Grenville Front Tectonic Zone near Killarney, Ontario. Journal of Geology 101, 575–589.[Web of Science]

Hanmer, S. & McEachern, S. J. (1992). Kinematical and rheological evolution of a crustal-scale ductile thrust zone, Central Metasedimentary Belt, Grenville orogen, Ontario. Canadian Journal of Earth Sciences 29, 1779–1790.[Web of Science]

Hanson, G. N. (1978). The application of trace elements to the petrogenesis of igneous rocks of granitic composition. Earth and Planetary Science Letters 38, 26–43.[CrossRef][Web of Science]

Hollister, L. S. (1993). The role of melt in the uplift and exhumation of orogenic belts. Chemical Geology 108, 31–48.[CrossRef][Web of Science]

Hollister, L. S. & Crawford, M. L. (1986). Melt-enhanced deformation: a major tectonic process. Geology 14, 558–561.[Abstract/Free Full Text]

Jahns, R. H. & Burnham, C. W. (1969). Experimental studies of pegmatite genesis: I. A model for the derivation and crystallization of granitic pegmatites. Economic Geology 64, 843–864.[Abstract/Free Full Text]

Jamieson, R. A., Beaumont, C., Nguyen, M. H. & Lee, B. (2002). Interaction of metamorphism, deformation and exhumation in large convergent orogens. Journal of Metamorphic Geology 20, 9–24.[CrossRef][Web of Science]

Johannes, W. & Holtz, F. (1990). Formation and composition of H2O-undersaturated granitic melts. In: Ashworth, J. R. & Brown, M. (eds) High-temperature Metamorphism and Crustal Anatexis. Mineralogical Society Series 2, 87–104.

Kenah, C. & Hollister, L. S. (1983). Anatexis in the Central Gneiss Complex, British Columbia. In: Atherton, M. P. & Gribble, C. D. (eds), Migmatites, Melting and Metamorphism. Nantwich: Shiva, pp. 142–162.

Ketchum, J. W. F., Heaman, L. M., Krogh, T. E., Culshaw, N. G. & Jamieson, R. A. (1998). Timing and thermal influence of late orogenic extension in the lower crust: a U–Pb geochronological study from the southwest Grenville orogen, Canada. Precambrian Research 89, 25–45.[CrossRef][Web of Science]

Kretz, R. (1983). Symbols for rock-forming minerals. American Mineralogist 68, 277–279.[Abstract]

Krogh, T. E. (1994). Precise U–Pb ages for Grenvillian and pre-Grenvillian thrusting of Proterozoic and Archean metamorphic assemblages in the Grenville Front Tectonic Zone, Canada. Tectonics 13, 963–982.[CrossRef][Web of Science]

Lappin, A. R. & Hollister, L. S. (1980). Partial melting in the Central Gneiss Complex near Prince Rupert, British Columbia. American Journal of Science 280, 518–545.[Abstract/Free Full Text]

Layman, A. J. (2003). Composition and significance of fluid inclusions in garnet, Algonquin Domain, Central Gneiss Belt, Grenville Province, Ontario. B.Sc.(Hons) thesis, Dalhousie University, Halifax.

Longerich, H. P., Jenner, G. A., Fryer, B. J. & Jackson, S. E. (1990). Inductively coupled plasma-mass spectrometric analysis of geological samples: a critical evaluation based on case studies. Chemical Geology 83, 105–118.[CrossRef][Web of Science]

Martin, H. (1987). Petrogenesis of Archean trondhjemites, tonalites, and granodiorites from eastern Finland: major and trace element geochemistry. Journal of Petrology 28, 921–953.[Abstract/Free Full Text]

McEachern, S. J. & van Breemen, O. (1993). Age of deformation within the Central Metasedimentary Belt boundary thrust zone, southwest Grenville orogen: constraints on the collision of the Mid-Proterozoic Elzevir terrane. Canadian Journal of Earth Sciences 30, 1155–1165.[Web of Science]

McLellan, E. L. (1988). Migmatite structures in the Central Gneiss Complex, Boca de Quadra, Alaska. Journal of Metamorphic Geology 6, 517–542.[Web of Science]

McMullen, S. M. (1999). Tectonic evolution of the Bark Lake area, eastern Central Gneiss Belt, Ontario Grenville: constraints from geology, geochemistry and U–Pb geochronology. M.Sc. thesis, Carleton University, 175 pp.

Mogk, D. W. (1992). Ductile shearing and migmatization at mid-crustal levels in an Archean high-grade gneiss belt, northern Gallatin Range, Montana, USA. Journal of Metamorphic Geology 10, 427–438.[Web of Science]

Montel, J. M., Marignac, C., Barbey, P. & Pichavant, M. (1992). Thermobarometry and granite genesis: the Hercynian low-P, high-T Velay anatectic dome (French Massif Central). Journal of Metamorphic Geology 10, 1–15.[Web of Science]

Moore, J. M. & Thompson, P. H. (1980). The Flinton group: a late Precambrian metasedimentary succession in the Grenville Province of eastern Ontario. Canadian Journal of Earth Sciences 17, 1685–1707.[Web of Science]

Nadeau, L. & van Breemen, O. (1998). Plutonic ages and tectonic setting of the Algonquin and Muskoka allochthons, Central Gneiss Belt, Grenville Province, Ontario. Canadian Journal of Earth Sciences 35, 1423–1438.[CrossRef][Web of Science]

Naney, M. T. (1983). Phase equilibria of rock-forming ferromagnesian silicates in granitic systems. American Journal of Science 238, 993–1033.

Nédélec, A., Minyem, D. & Barbey, P. (1993). High-P–high-T anatexis of Archean tonalitic grey gneisses: the Eseka migmatites, Cameroon. Precambrian Research 62, 191–205.[CrossRef][Web of Science]

Olsen, S. N. (1984). Mass-balance and mass-transfer in migmatites from the Colorado Front Range. Contributions to Mineralogy and Petrology 85, 30–44.[CrossRef][Web of Science]

Patiño Douce, A. E. & Beard, J. S. (1995). Dehydration-melting of biotite gneiss and quartz amphibolite from 3 to 15 kbar. Journal of Petrology 36, 707–738.[Abstract/Free Full Text]

Pattison, D. R. M. (1991). Infiltration-driven dehydration and anatexis in granulite facies metagabbro, Grenville Province, Ontario, Canada. Journal of Metamorphic Geology 9, 315–332.[Web of Science]

Peterson, J. W. & Newton, R. C. (1989). Reversed experiments on biotite–quartz–feldspar melting in the system KMASH: implications for crustal anatexis. Journal of Geology 97, 465–485.[Web of Science]

Petford, N., Kerr, R. C. & Lister, J. R. (1993a). Dike transport of granitoid magmas. Geology 21, 845–848.[Abstract/Free Full Text]

Petford, N., Lister, J. R. & Kerr, R. C. (1993b). The ascent of felsic magmas in dykes. Lithos 32, 161–168.

Rainbird, R. H., McNicoll, V. J., Thériault, R. J., Heaman, L. M., Abbott, J. G., Long, D. G. F. & Thorkelson, D. J. (1997). Pan-continental river system draining Grenville orogen recorded by U–Pb and Sm–Nd geochronology of Neoproterozoic quartzarenites and mudrocks, northwestern Canada. Journal of Geology 105, 1–17.[Web of Science]

Rivers, T. (1997). Lithotectonic elements of the Grenville Province: review and tectonic implications. Precambrian Research 86, 117–154.[CrossRef][Web of Science]

Rivers, T. & Corrigan, D. (2000). Convergent margin on southeastern Laurentia during the Mesoproterozoic: tectonic implications. Canadian Journal of Earth Sciences 37, 359–383.[CrossRef][Web of Science]

Robin, P.-Y. F. (1979). Theory of metamorphic segregation and related processes. Geochimica et Cosmochimica Acta 43, 1587–1600.[CrossRef][Web of Science]

Ross, G. M., Parrish, R. R. & Winston, D. (1992). Provenance and U–Pb geochronology of the Mesoproterozoic Belt Supergroup (northwestern United States): implications for the age of deposition and pre-Panthalassa plate reconstructions. Earth and Planetary Science Letters 113, 57–76.[CrossRef][Web of Science]

Rutter, M. J. & Wyllie, P. J. (1988). Melting of vapour-absent tonalite at 10 kbar to simulate dehydration-melting in the deep crust. Nature 331, 159–160.[CrossRef]

Sandiford, M. & McLaren, S. (2002). Tectonic feedback and the ordering of heat producing elements within the continental lithosphere. Earth and Planetary Science Letters 204, 133–150.[CrossRef][Web of Science]

Sawyer, E. W. (1987). The role of partial melting and fractional crystallization in determining discordant leucosome compositions. Journal of Petrology 28, 445–473.[Abstract/Free Full Text]

Sawyer, E. W. (1991). Disequilibrium melting and the rate of melt–residuum separation during migmatization of mafic rocks from the Grenville Front, Quebec. Journal of Petrology 32, 701–738.[Abstract/Free Full Text]

Sawyer, E. W. (1994). Melt segregation in the continental crust. Geology 22, 1019–1022.[Abstract/Free Full Text]

Sawyer, E. W. (1998). Formation and evolution of granite magmas during crustal reworking: the significance of diatexites. Journal of Petrology 39, 1147–1167.[CrossRef][Web of Science]

Sawyer, E. W. (1999). Criteria for the recognition of partial melting. Physics and Chemistry of the Earth (A) 24, 269–279.

Sawyer, E. W. (2001). Melt segregation in the continental crust: distribution and movement of melt in anatectic rocks. Journal of Metamorphic Geology 19, 291–309.[CrossRef][Web of Science]

Sawyer, E. W., Dombrowski, C. & Collins, W. J. (1999). Movement of melt during synchronous regional deformation and granulite-facies anatexis, an example from the Wuluma Hills, central Australia. In: Castro, A., Fernandez, C. & Vigneresse, J. L. (eds) Understanding Granites: Integrating New and Classical Techniques. Geological Society, London, Special Publications 168, 221–237.

Shaw, D. M. (1970). Trace element fractionation during anatexis. Geochimica et Cosmochimica Acta 34, 237–243.[CrossRef][Web of Science]

Simmons, W. B., Lee, M. T. & Brewster, R. H. (1987). Geochemistry and evolution of the South Platte granite–pegmatite system, Jefferson County, Colorado. Geochimica et Cosmochimica Acta 51, 455–471.[CrossRef][Web of Science]

Skjerlie, K. P. & Johnston, A. D. (1996). Vapour-absent melting from 10 to 20 kbar of crustal rocks that contain multiple hydrous phases: implications for anatexis in the deep to very deep continental crust and active continental margins. Journal of Petrology 37, 661–691.[Abstract/Free Full Text]

Slagstad, T. (2003). Muskoka and Shawanaga domains, Central Gneiss Belt, Grenville Province, Ontario: geochemical and geochronological constraints on pre-Grenvillian and Grenvillian geological evolution. Ph.D. thesis, Dalhousie University, Halifax.

Slagstad, T., Culshaw, N. G., Jamieson, R. A. & Ketchum, J. W. F. (2004a). Early Mesoproterozoic tectonic history of the southwestern Grenville Province, Ontario: constraints from geochemistry and geochronology of high-grade gneisses. In: Tollo, R. P., Corriveau, L., McLelland, J. & Bartholomew, M. J. (eds) Proterozoic Tectonic Evolution of the Grenville Orogen in North America. Geological Society of America, Memoir 197, 209–241.

Slagstad, T., Hamilton, M. A., Jamieson, R. A. & Culshaw, N. G. (2004b). Timing and duration of melting in the mid orogenic crust: constraints from U–Pb (SHRIMP) data, Muskoka and Shawanaga domains, Grenville Province, Ontario. Canadian Journal of Earth Sciences 41, 1339–1365.[CrossRef][Web of Science]

Solar, G. S. & Brown, M. (2001). Petrogenesis of migmatites in Maine, USA: possible source of peraluminous leucogranite in plutons. Journal of Petrology 42, 789–823.[Abstract/Free Full Text]

Stevens, G. & Clemens, J. D. (1993). Fluid-absent melting and the roles of fluids in the lithosphere: a slanted summary? Chemical Geology 108, 1–17.[CrossRef][Web of Science]

Sun, S.-S. & McDonough, W. F. (1989). Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes. In: Saunders, A. D. & Norry, M. J. (eds) Magmatism in the Ocean Basins. Geological Society, London, Special Publications 42, 313–345.

Symmes, G. H. & Ferry, J. M. (1995). Metamorphism, fluid flow and partial melting in pelitic rocks from the Onawa contact aureole, central Maine, USA. Journal of Petrology 36, 587–612.[Abstract/Free Full Text]

Thompson, A. B. (1982). Dehydration melting of pelitic rocks and the generation of H2O-undersaturated granitic liquids. American Journal of Science 282, 1567–1595.[Abstract/Free Full Text]

Timmermann, H. (1998). Geology, metamorphism, and U–Pb geochronology in the Central Gneiss Belt between Huntsville and Haliburton, southwestern Grenville Province, Ontario. Ph.D. thesis, Dalhousie University, Halifax, 410 pp.

Timmermann, H., Parrish, R. R., Jamieson, R. A. & Culshaw, N. G. (1997). Time of metamorphism beneath the Central Metasedimentary Belt boundary thrust zone, Grenville orogen, Ontario: accretion at 1080 Ma? Canadian Journal of Earth Sciences 34, 1023–1029.[Web of Science]

Timmermann, H., Jamieson, R. A., Parrish, R. R. & Culshaw, N. G. (2002). Coeval migmatites and granulites, Muskoka domain, southwestern Grenville Province. Canadian Journal of Earth Sciences 39, 239–258.[CrossRef][Web of Science]

van Breemen, O. & Hanmer, S. (1986). Zircon morphology and U–Pb geochronology in active shear zones: studies of syntectonic intrusions along the northwest boundary of the Central Metasedimentary Belt, Grenville Province, Ontario. Current Research, Part B, Geological Survey of Canada, Paper 86-1B, 775–784.

van der Pluijm, B. A. & Carlson, K. A. (1989). Extension in the Central Metasedimentary Belt of the Ontario Grenville: timing and tectonic significance. Geology 17, 161–164.[Abstract/Free Full Text]

Vielzeuf, D., Clemens, J. D., Pin, C. & Moinet, E. (1990). Granites, granulites and crustal differentiation. In: Vielzeuf, D. & Vidal, P. (eds), Granulites and Crustal Evolution. Dordrecht: Kluwer Academic, pp. 551–568.

Waters, D. J. (2001). The significance of prograde and retrograde quartz-bearing intergrowth microstructures in partially melted granulite-facies rocks. Lithos 56, 97–110.[CrossRef][Web of Science]

Watson, E. B. & Harrison, T. M. (1983). Zircon saturation revisited: temperature and composition effects in a variety of crustal magma types. Earth and Planetary Science Letters 64, 295–304.[CrossRef][Web of Science]

Watson, E. B., Vicenzi, E. P. & Rapp, R. P. (1989). Inclusion/host relations involving accessory minerals in high-grade metamorphic and anatectic rocks. Contributions to Mineralogy and Petrology 101, 220–231.[CrossRef][Web of Science]

White, D. J., Easton, R. M., Culshaw, N. G., Milkereit, B., Forsyth, D. A., Carr, S. D., Green, A. G. & Davidson, A. (1994). Seismic images of the Grenville Orogen in Ontario. Canadian Journal of Earth Sciences 31, 293–307.[Web of Science]

Wickham, S. M. (1987). The segregation and emplacement of granitic magmas. Journal of the Geological Society, London 144, 281–297.[Abstract/Free Full Text]

Wodicka, N., Ketchum, J. W. F. & Jamieson, R. A. (2000). Grenvillian metamorphism of monocyclic rocks, Georgian Bay, Ontario, Canada: implications for convergence history. Canadian Mineralogist 38, 471–510.[CrossRef][Web of Science]

Wolf, M. B. & Wyllie, P. J. (1994). Dehydration-melting of amphibolite at 10 kbar: the effects of temperature and time. Contributions to Mineralogy and Petrology 115, 369–383.[CrossRef][Web of Science]

Wynne-Edwards, H. R. (1972). The Grenville Province. In: Price, R. A. & Douglas, R. J. W. (eds) Variations in Tectonic Styles in Canada. Geological Association of Canada Special Paper 11, 263–334.


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
Journal of the Geological SocietyHome page
M. Brown
Crustal melting and melt extraction, ascent and emplacement in orogens: mechanisms and consequences
Journal of the Geological Society, July 1, 2007; 164(4): 709 - 730.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow All Versions of this Article:
46/5/893    most recent
egi004v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (3)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by SLAGSTAD, T.
Right arrow Articles by CULSHAW, N. G.
Right arrow Search for Related Content
GeoRef
Right arrow GeoRef Citation
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?