Hostname: page-component-745bb68f8f-kw2vx Total loading time: 0 Render date: 2025-02-06T09:43:31.841Z Has data issue: false hasContentIssue false

Formation of spinel-cordierite-plagioclase symplectites replacing andalusite in metapelitic migmatites of the Alvand aureole, Iran

Published online by Cambridge University Press:  22 October 2010

A. SAKI*
Affiliation:
Department of Geology, Faculty of Science, Shahid Chamran University, Ahvaz, 65355-141, Iran
Rights & Permissions [Opens in a new window]

Abstract

Spinel–cordierite–plagioclase symplectites partially replacing andalusite occur in the metapelitic migmatite rocks of the Alvand aureole within the Sanandaj–Sirjan metamorphic belt, Hamadan, Iran. The presence of melt shows that corona development occurred under partial melting conditions. Spinel is predicted to grow with cordierite at around 700°C. Exhaustion of the available SiO2 and/or separation of sillimanite/andalusite from SiO2-rich matrix domains by cordierite resulted in the formation of localized low-silica activity domains and thus triggered the growth of spinel in the rim of andalusite, the reaction Sil/And + Bt = Crd + Spl + Kfs + melt, as the most common reaction for the development of coronas in the metapelitic of Alvand aureole. The breakdown of garnet to plagioclase + sillimanite, dehydration melting and the formation of spinel–plagioclase symplectite could occur during heating or decompression; these textures are limited to the contact aureole in the studied area, so heating is perhaps the more likely explanation for formation of the symplectites in the metapelitic rocks of the Alvand aureole. The P–T diagram, inferred paths and zoning profiles of garnet do not account for the decompression history of the terrane.

Type
Original Article
Copyright
Copyright © Cambridge University Press 2010

1. Introduction

The mineral assemblage garnet–spinel–cordierite–feldspar–quartz–biotite–sillimanite is typical of metapelitic systems at high temperature (HT) and low (LP) to medium pressure (MP) (e.g. Hand et al. Reference Hand, Scrimgeour, Powell, Stüwe and Wilson1994; McDade & Harley, Reference McDade and Harley2001; White, Powell & Holland, Reference White, Powell and Holland2001; White, Powell & Clarke, Reference White, Powell and Clarke2002; Cenki, Kriegsman & Braun, Reference Cenki, Kriegsman and Braun2002; Ings & Owen, Reference Ings and Owen2002; Riesco et al. Reference Riesco, Stüwe, Reche and Martinez2004). Such P–T conditions are also propitious for the assemblage to be accompanied by melt, derived by partial melting of the rocks. A typical microstructural feature in these metapelites is the intergrowth of spinel and cordierite, which usually replaces garnet–sillimanite-bearing assemblages (e.g. Mezger, Chacko & Erdmer, Reference Mezger, Chacko and Erdmer2001; Topuz et al. Reference Topuz, Altherra, Kalta, Satirb, Wernera and Schwarza2004), and/or andalusite (Cesare et al. Reference Cesare, Marchesi, Hermann and Gomez-Pungnaire2003; Johnson et al. Reference Johnson, Brown, Gibson and Wing2004). Symplectites of spinel and cordierite replacing andalusite and/or sillimanite have been reported from many high-temperature regional metamorphic terranes and several contact aureoles, and commonly these symplectites are inferred to record decompression with reaction occurring within restricted silica-deficient compositional domains focused around aluminosilicate (e.g. Bucher-Nurminen & Droop, Reference Bucher-Nurminen and Droop1983; Droop, Reference Droop1989; Grant & Frost, Reference Grant and Frost1990; Clarke & Powell, Reference Clarke and Powell1991; Waters, Reference Waters1991; Carson et al. Reference Carson, Powell, Wilson and Dirks1997; Greenfield, Clarke & White, Reference Greenfield, Clarke and White1998; Pitra & de Waal, Reference Pitra and De Wall2001; Norlander et al. Reference Norlander, Whitney, Teyssier and Vanderhaeghe2002; White, Powell & Clarke, Reference White, Powell and Clarke2003). Spinel–cordierite assemblages may occur as coronitic and/or symplectite microstructures, related to particular kinetic conditions of a low chemical diffusion rate (Hand et al. Reference Hand, Dirks, Powell and Buick1992; Passchier & Trouw, Reference Passchier and Trouw1996; Vernon, Reference Vernon1996; Whittington, Harris & Baker, Reference Whittington, Harris, Baker, Treolar and O'Brien1998). The presence of such symplectites and the correct interpretation of their petrogenesis have consequences for understanding the tectonothermal evolution of segments of the crust, and in investigating equilibration volumes and diffusion-controlled fractionated domainal bulk compositions (e.g. Stüwe, Reference Stüwe1997; Marmo, Clarke & Powell, Reference Marmo, Clarke and Powell2002). Spinel–cordierite symplectites replacing andalusite are found locally in metapelites within the migmatite rocks in the Alvand aureole. However, the occurrence of such symplectites within the aureole of the Alvand Complex, and the absence of such symplectites within regional metamorphic rocks elsewhere in the Hamadan area, raises doubts about the general applicability of this model. In this paper, we investigate in detail migmatized andalusite- and symplectite-bearing metapelitic rocks from the Alvand aureole within the Sanandaj–Sirjan metamorphic belt, Hamadan.

2. Geological setting

The study area is a part of the 1500 km long Sanandaj–Sirjan metamorphic belt of the Zagros orogen of western Iran (Fig. 1a, b). The Zagros formed during subduction of a Neo-Tethyan seaway and subsequent oblique collision of Afro-Arabia (Gondwana) with the Iranian microcontinent in Late Cretaceous–early Tertiary times (Berberian & King, Reference Berberian and King1981; Alavi, Reference Alavi1994; Mohajjel & Fergusson, Reference Mohajjel and Fergusson2000). Crustal shortening associated with subduction and collision metamorphosed and deformed Late Palaeozoic and Mesozoic sedimentary rocks, including a sequence of pelitic, psammitic, mafic, calc-pelitic and calc-silicate rocks near Hamadan (Fig. 1b).

Figure 1. (a) Location of the Sanandaj–Sirjan belt in Iran. (b) Geological map of the area (modified from Baharifar et al. Reference Baharifar, Moinevaziri, Bellon and Pique2004), with sample locations and stratigraphic column of the area. (c) Generalized geological map of the study area near Hamadan. Isograds (dashed lines) and metamorphic zones are shown (modified from Sepahi, Whitney & Baharifar, Reference Sepahi, Whitney and Baharifar2004). Some faults are shown (heavy lines within unit 4); contacts between zones 4a–4b, 4a–5, 2–5, and possibly 4–7 are faults. Zones: 1 – biotite zone; 2 – garnet zone; 3 – andalusite zone (±fibrolite); 4 – staurolite zone (a = staurolite + andalusite ± fibrolite; b = staurolite + garnet, no andalusite); 5 – sillimanite–andalusite zone; 6 – sillimanite–Kfs zone; 7 – cordierite zone (a = cordierite with no new growth of Al2SiO5; b = cordierite + neoblastic andalusite + sillimanite); 8 – cordierite Kfs zone; 9 – spotted schist; 10 – spinel + orthopyroxene zone (migmatites).

During the Palaeozoic period, the Sanandaj–Sirjan Zone (SSZ) was part of northeast Gondwanaland, separated from the Eurasian plate by the Palaeo-Tethys ocean (Golonka, Reference Golonka2004; Sears, George & Winne, Reference Sears, George and Winne2005; Shahabpour, Reference Shahabpour2005, Ahmadi et al. Reference Ahmadi, Esmaeily, Valizadeh and Rahimpour-Bonab2007; Fazlnia et al. Reference Fazlnia, Schenk, Straaten and Mirmohammadi2009). From early Permian to early Triassic times, the SSZ was situated along the southern margin of the Eurasian plate, separated from northern Gondwanaland by the Neo-Tethyan ocean (Mohajjel, Fergusson & Sahandi, Reference Mohajjel, Fergusson and Sahandi2003; Khalaji et al. Reference Khalaji, Esmaeily, Valizadeh and Rahimpour-Bonab2007; Fazlnia et al. Reference Fazlnia, Schenk, Straaten and Mirmohammadi2009). During the Mesozoic, the oceanic crust of the Neo-Tethys was subducted beneath the Eurasian plate (Mohajjel, Fergusson & Sahandi, Reference Mohajjel, Fergusson and Sahandi2003; Golonka, Reference Golonka2004; Shahabpour, Reference Shahabpour2005; Ahmadi et al. Reference Ahmadi, Esmaeily, Valizadeh and Rahimpour-Bonab2007; Fazlnia et al. Reference Fazlnia, Schenk, Straaten and Mirmohammadi2009), and the SSZ occupied the position of a magmatic arc (Berberian & Berberian, Reference Berberian, Berberian, Gupta and Delany1981; Arvin et al. Reference Arvin, Pan, Dargahi, Malekzadeh and Babaei2007). The final closure of the Neo-Tethys and the collision of the Arabian and Eurasian plates took place during the Tertiary period (Berberian & Berberian, Reference Berberian, Berberian, Gupta and Delany1981), and during the same period, the Zagros Fold and Thrust Belt formed as part of the Alpine–Himalayan mountain chain, extending to about 2000 km from eastern Turkey to the Oman line in southern Iran (Agard et al. Reference Agard, Omrani, Joliver and Mouthereau2005; Mohajjel, Fergusson & Sahandi, Reference Mohajjel, Fergusson and Sahandi2003; Ahmadi et al. Reference Ahmadi, Esmaeily, Valizadeh and Rahimpour-Bonab2007; Omrani et al. Reference Omrani, Agard, Whitechurch, Benoit, Prouteau and Jolivet2008).

The Sanandaj–Sirjan zone contains abundant Mesozoic arc-related plutonic and plutono-metamorphic complexes, such as Alvand, Almogholagh, Samen, Ghorveh, Borojerd, Urumieh, Arak, Astaneh, Qori and Siah-Kuh.

The U–Pb dating results demonstrate that all granitoids were exclusively emplaced during the Jurassic period instead of being Cretaceous or younger in age as suggested previously (Shahbazi et al. Reference Shahbazi, Siebel, Pourmoafee, Ghorbani, Sepahi, Shang and Abedini2010). The pluton was assembled incrementally over about 10 Ma. Gabbros formed at 166.5±1.8 Ma, granites between 163.9±0.9 Ma and 161.7±0.6 Ma, and leucocratic granitoids between 154.4±1.3 and 153.3±2.7 Ma. Granites and leucocratic granitoids show some A-type affinity. It is concluded that the Alvand plutonic complex was generated in a continental-arc-related extensional regime during subduction of Neo-Tethyan oceanic crust beneath the SSZ (Shahbazi et al. Reference Shahbazi, Siebel, Pourmoafee, Ghorbani, Sepahi, Shang and Abedini2010).

The tectonic evolution of the Sanandaj–Sirjan belt involved continental arc magmatism followed by collision. Mafic to intermediate plutonic bodies (olivine gabbro, gabbro, gabbro-norite, diorite, quartz diorite and tonalite) are older than crustally derived granitic plutons in the region (Alvand Plutonic Complex). All metamorphic rocks in the area that were regionally metamorphosed prior to the contact metamorphism up to greenschist to amphibolite facies underwent contact metamorphism due to emplacement of the Alvand Complex. The plutons, including the granites, are commonly associated with contact aureoles defined by migmatites, hornfelsic textures and mineral assemblages that overprint earlier minerals and fabrics. Migmatitic rocks are located in the southwestern part of the Alvand Complex near the city of Toyserkan (Fig. 1b).

3. Field observations

In the field area near Hamadan, igneous and metapelitic rocks are the most abundant rock type and are interlayered with minor metabasaltic rocks (amphibole schist and amphibolite), metacarbonates and calc-silicate rocks.

Calc-silicate rocks have formed in the eastern part of this batholite. The peak metamorphic assemblage of calc-silicate rocks is garnet, diopside, vesuvinite, wollastonite and calcite. P–T results indicate that a contact metamorphism condition occurs in pyroxene hornfels facies in this part of the Hamadan crust (Saki & Pourkaseb, Reference Saki and Pourkaseb2010).

Critical metamorphic mineral assemblages in the metapelitic rocks indicate regional metamorphism under amphibolite-facies conditions. The dominant metamorphic assemblage of regional metapelitic rocks is Ms + Bt + Grt + St + And + Sil/Ky.

Contact metamorphism overprinted the older regional metamorphic assemblages. The main peak contact metamorphic assemblage is Grt + Bt+ And + Ms + Sil + Crd + Spl ± Opx.

Metapelitic rocks occur as slate, phyllite, mica schist, garnet schist, garnet–andalusite (±sillimanite or kyanite) schist, garnet–staurolite schist and garnet–sillimanite (±kyanite) schist, hornfelsic and migmatitic rocks (near the Alvand Plutonic Complex). These contact zone rocks include cordierite + K-feldspar (±andalusite, fibrolite) and garnet–staurolite (±kyanite) hornfels. In the Hamadan metamorphic sequence, cordierite and K-feldspar are only found near (<3 km from) the Alvand granite. The general trend of increasing metamorphic grade is towards the pluton, with a succession of isograds in metapelitic rocks: biotite-in, garnet-in, andalusite-in, staurolite-in, cordierite-in, K-feldspar-in and spinel-in. In some parts of the field area, the complete sequence occurs, but in other areas, some of the lower grade zones are missing (Fig. 1c).

Faulting has disrupted the sequences in places, and some apparent isograds may actually be faults. The most obvious example of disruption of the sequence by faulting is the fault-bounded panel of andalusite–sillimanite schists (zone 5) that occurs between the andalusite and staurolite zones (Fig. 1c).

As mentioned above, pelitic rocks, containing pale quartz-feldspar (granitoid) veins, appear in the upper Spinel Zone and Sillimanite Zone. Because these rocks are mixed rocks with quartz-feldspathic parts (leucosomes) and darker coloured parts, more enriched in ferromagnesian minerals, we will use the term ‘migmatite’ (Sederholm, Reference Sederholm1967) for them.

The migmatites are best developed on the southern side of the Alvand complex (Fig. 2). Migmatitic metamorphic zones occur at several localities that are in close proximity to batholithic contacts within the Alvand aureole pelites of the Hamadan region (Saki, Moazzen & Baharifar, unpub. data). The migmatite zone shown in Figure 1c is the most extensive of these and can be related to a very well-defined metamorphic gradient through the metapelite unit.

Figure 2. Field photographs of migmatites. A colour version of this figure is available at http://www.cambridge.org/journals/geo. (a, b) Migmatite in which smaller leucosomes in (a) and (b) appear to be crystallized from melt derived locally, whereas the larger leucosome likely records evidence of melt moving through the system. (c, d) Plane-polarized light view of a Spl + Crd + And corona in the mesosome. (e) Sil + Crd + Spl in the mesosome. (f) The garnet as inclusion in the cordierite. (g) Feldspar–plagioclase–quartz patches (Lu), with the melt apparently segregated into mesoscopic patches or veins. (h) Euhedral shape of the plagioclase phenocrysts against intersertal quartz in (g).

The leucosomes that developed during migmatization can be morphologically divided into: (1) concordant (stromatitic) and pervasive network leucosomes (Fig. 2a), (2) discordant dykes and pockets of pegmatites and aplites (Fig. 2b), and (3) agmatitic and randomly distributed leucosome patches of various sizes, which appear to represent the gathering of quartzo-feldspathic segregations into larger, granitic or pegmatitic tubes and dykes.

4. Petrography

Petrographic characteristics of the symplectite-bearing rocks are based on observations of several samples. The rocks are composed of three mesoscopic components: a cordierite-, K-feldspar, andalusite and biotite-rich matrix (mesosome); partial pseudomorphs in which resorbed andalusite is surrounded by spinel–cordierite symplectites; quartz-, K-feldspar- and/or plagioclase-rich leucosome that commonly shows a close spatial relationship with the pseudomorphs but which also occurs as pockets within the matrix.

4.a. Mesosomes

Mesosomes in patches and stromatic migmatites (patch migmatites formed by in situ melting) are massive (Fig. 2a, b). The matrix (mesosomes) contains cordierite, K-feldspar, biotite, andalusite and quartz with minor plagioclase, garnet and ilmenite. Cordierite forms randomly oriented subhedral grains up to 10 mm in length and has garnet inclusions (Fig. 2f). At least two microstructural variants of cordierite occur in these samples: (1) cordierite as a constituent of the coronas around andalusite (Fig. 2c, d); (2) cordierite crystals in the matrix, with inclusions of melt, biotite and sillimanite (Fig. 2e). Andalusites have equidimensional shape and resorbed rims, and are generally <20 mm in length. Locally, spinel has a thin rim of plagioclase. K-feldspar forms anhedral grains up to 2 mm across; cores are riddled with finely disseminated opaque inclusions and larger inclusions of quartz and biotite and exhibit a micro-perthitic microstructure.

4.b. Leucosomes

Migmatites in the Alvand aureole exhibit different leucosome shapes. The leucosomes exhibit mesoscopic (Fig. 2a, b) to microscopic (Fig. 2g) structure, are fine- to coarse-grained and show magmatic texture (Fig. 2h). They are leucotonalitic to granitic in composition, as some strings and layers consist of distinct portions that are either rich in Kfs or in Pl or are made up of Qtz + Grt + two feldspars. This seems to indicate that melt segregation was not completely efficient and that melt did not migrate over long distances; some migration did occur, however, as indicated by the discordant and concordant leucosomes and by the variable concentration of leucosomes at outcrop scale. Leucosomes are composed of quartz, plagioclase and K-feldspar in different proportions. Pl may be the only feldspar in some leucotonalitic types; Kfs exceeds 50 vol.% within some granitic bodies; Qtz, independently of the composition of the rocks, may exceed 60 vol.%. Garnet is present in leucosomes as coarse crystals up to 1 cm in diameter. The granitic bodies associated with metapelites are composed of large subhedral and embayed crystals of Kfs, more rarely Pl with interstitial fine- to medium-grained Qtz, Pl and Kfs and large euhedral crystals of garnet. Leucosomes preserve textures resulting from the presence of melt, such as euhedral phenocrysts/peritectic phases, textures that mimic solid–melt relations (quartz/feldspar films along grain boundaries, interstitial–xenomorphic quartz, and residual quartz in feldspar), melt–solid reaction textures and fractured residual grains (Saki & Baharifar, Reference Saki and Baharifar2010).

4.c. Pseudomorphs

The partial pseudomorphs of andalusite are up to 10 mm long and make up 5–10 vol.% of the rock. They comprise a variably resorbed andalusite surrounded by a rim of symplectite composed of spinel, cordierite and plagioclase that may be in excess of 1 mm wide (Fig. 2e, f). The symplectites are surrounded by a moat of pure cordierite, up to 1 mm wide, in which individual stubby cordierite grains are aligned perpendicular to the symplectite–andalusite interface. The symplectite rim and cordierite moat are surrounded by K-feldspar- and/or plagioclase-rich, quartz-poor leucosomes and/or by a coarse-grained biotite fringe that grades into K-feldspar-poor matrix (Fig. 2g, h).

4.d. Corona microstructure

The corona microstructure is composed of a partly to completely resorbed andalusite up to 20 mm in diameter which shows a corona of Spl + Crd + Pl that is only present between andalusite and the matrix. The corona displays a shell of subhedral spinel within a matrix of cordierite, K-feldspar and plagioclase. Spinel is in contact with the matrix or with andalusite (Fig. 2c, d). In a few cases, spinel in the corona shows a thin rim of plagioclase. The melt was present throughout the history of the andalusite formation and its subsequent partial replacement by coronitic phases.

5. Mineralogy

Major element compositions of biotite, garnet, spinel, plagioclase, cordierite and ilmenite in selected samples (migmatite rocks) were determined by wavelength-dispersive spectrometry using a Cameca SX100 microprobe. The operational conditions were 15 kV, 10–20 nA specimen current. The analytical spot diameter was set between 3 and 5 μm, keeping the same current conditions. Natural and synthetic standards were used for calibration. Ferric iron contents in ferromagnesian minerals are calculated using stoichiometric criteria (Droop, Reference Droop1989). Mineral abbreviations are from Kretz (Reference Kretz1983).

5.a. Garnet

Garnet compositions are listed in Table 1. In general, garnets are Fe-rich, with almandine content varying between 80 and 90%. Ti is not presented or occurs in very small amounts. Figure 3 shows the flat zoning profiles of garnet in the migmatite rocks (contact metamorphic rocks). Zoning profiles are from core to rim of the contact and regional garnets (Fig. 4).

Table 1. Representative mineral analyses of garnet and cordierite

Fe2+/Fe3+ ratio in garnet is calculated using stoichiometry.

Figure 3. Zoning profile from core to rim through the garnet in some samples of migmatite (mesosome; contact metamorphism), which shows a flat profile attributed to garnet crystallization in these rocks.

Figure 4. Zoning profile of the garnet in some samples from the regional metapelitic rocks, which shows prograde metamorphism (normal zoning) attributed to garnet crystallization in different regional metamorphic zones.

The pronounced normal zoning in the regional metamorphic garnet from some samples can be seen in Figure 4.

5.b. Cordierite

Cordierite compositions are listed in Table 1. The average EMP total of cordierite is >98.5. In the K-feldspar–cordierite, sillimanite–K-feldspar and spinel–cordierite zones, cordierite is homogeneous in thin-section. Cordierite in the muscovite–cordierite zone is completely pinitized. XMg of cordierite shows no systematic change with ascending metamorphic grade from the K-feldspar–cordierite zone to the spinel–cordierite zone (0.45–0.62).

5.c. Biotite

Table 2 includes the microprobe analyses of biotite in the studied samples. Biotite has high TiO2 contents, ranging from 3.02 to 3.25 wt%. Phlogopite–annite is the dominant constituent in the biotite. Classification after Robinson et al. (Reference Robinson, Hollocher, Tracy, Dietsch, Joester and Quarrier1982) indicates that peak metamorphism biotite belongs to the Sil–Kfs and Grt–Crd–Kfs zone type.

Table 2. Representative mineral analyses of biotite, spinel and ilmenite

For biotite, Fe total is considered as Fe2+.

5.d. Spinel

Spinel compositions are listed in Table 2. This is a hercynite-rich solid solution that has compositions in the range X(Her) 0.86–0.89.

5.e. Ilmenite

This occurs in all the zones as an accessory phase, not only in pelitic but also in siliceous rocks. Electron-microprobe analyses of ilmenite are provided in Table 2; the ilmenites have a rather constant composition, with XIlm = 0.93–0.96. The MgO content is <2.02 wt% in all the zones.

5.f. Feldspar

Electron-microprobe analyses of feldspar are provided in Table 3. The anorthite content (XAn) of plagioclase is variable; for example, Ca/(Na+K+Ca) is 0.33–0.88. The plagioclases are rich in anorthite. Plagioclases with high Ca content are used for barometry (see below).

Table 3. Representative mineral analyses of feldspar

6. Conventional geothermobarometry

The Hamadan metamorphic rocks have experienced multiple episodes of metamorphism driven by burial and heating during arc construction and collision, and these events are associated with local partial melting (at high grades, near the pluton) and infiltration of aqueous fluids. Determination of the peak P–T conditions and the paths for these various events was carried out using conventional thermobarometry methods such as cation exchange reaction thermometry, multiple equilibria calculations, GASP barometry and petrogenetic grids. Fe–Mg exchange between garnet and biotite thermometers using different calibrations was applied. Solution models of Helffrich & Wood (Reference Helffrich and Wood1989) and Ganguly & Saxena (Reference Ganguly and Saxena1984) were used for garnets. For pressure estimation, GASP (garnet–sillimanite–quartz–plagioclase; Koziol & Newton, Reference Koziol and Newton1988) barometers were employed. The outer parts of garnet along with biotites with highest Ti content were used to find the peak temperature, however; plagioclase does not show significant zoning. P–T estimation for peak metamorphic paragenesis gives a temperature of 650–700 °C and pressure of 2 kbar. Furthermore, using multiple equilibria (THERMOCALC program version 3.2; Holland & Powell, Reference Holland and Powell1998), temperature and pressure were calculated for the formation of the metapelitic rocks (Fig. 5). Table 4 shows the activity values that were used for calculating equilibrium reactions by THERMOCALC. The results are in good agreement with those from other methods. Activities of biotite, muscovite, garnet and plagioclase end-members were calculated using the AX program of Holland & Powell (Reference Holland and Powell1998). On the basis of these data, the studied rocks were at a depth of 6–12 km when they experienced a peak metamorphic temperature. This reflects a mean geothermal gradient of about 58°C/km, which is in the range of contact aureoles (Spear, Reference Spear1993). Therefore, peak metamorphism occurred under high-temperature and low-pressure conditions (LP/HT metamorphism) and reflects the high heat flow in this part of the crust. This suggests a genetic link between granite ascent from the deep crust and the high-temperature/low-pressure (contact metamorphism) type of metamorphism in the upper crust.

Figure 5. Representative results of multiple equilibria calculation by THERMOCALC (3.2) program; a pressure of ~2 kbar and temperature of ~750°C can be considered for peak metamorphism.

Table 4. Activity values used for calculating equilibrium reactions by THERMOCALC

py – pyrope; alm – almandine; ann – annite; phl – phlogopite; east – eastonite; crd – cordierite; fCrd – Fe-cordierite; An – anorthite; Ab – albite; Spl – spinel; herc – hercynite; mt – magnetite; ilm – ilmenite; Opx – orthopyroxene; Sil – sillimanite; Qtz – quartz

7. Discussion

The observations in the studied area are consistent with a thermal-tectonic history in which early regional metamorphism (greenschist to amphibolite facies) occurred during initial collision/contraction of continental margin basins and development of a continental margin arc and associated high geothermal gradient. Intermediate stages of contractional history and arc evolution were associated with intrusion of the Alvand Plutonic Complex, including the granite in the Hamadan area that created a low-pressure–high-temperature contact aureole (Sepahi, Whitney & Baharifar, Reference Sepahi, Whitney and Baharifar2004). Therefore, metamorphic rocks of the Hamadan have experienced multiple episodes of metamorphism driven by burial and heating during arc construction and collision, and these events are associated with local partial melting (at high grades, near the pluton).

The metamorphic conditions of the low-grade andalusite-bearing rocks (regional metamorphic rocks) are therefore not well defined by mineral assemblages, but peak temperatures and pressures were likely in the range of 375–500°C and 2–3 kbar, based on petrogenetic grid considerations, estimates of geothermal gradients, and calculations and inferences for rocks in higher grade zones. For example, garnet–staurolite schist in zone 4 records garnet–biotite Fe–Mg exchange temperatures of 520–570°C and a pressure of 3 kbar (A. A. Baharifar, unpub. M.Sc. thesis, Tarbiat Moallem Univ. Tehran, Iran, 1997).

The Alvand Plutonic Complex was emplaced at shallow crustal levels and is associated with a low-pressure–high-temperature contact aureole (A. A. Baharifar, unpub. M.Sc. thesis, Tarbiat Moallem Univ. Tehran, Iran, 1997; A. A. Sepahi, unpub. Ph.D. thesis, Tarbiat Moallem Univ. Tehran, Iran, 1999). Cordierite occurs as a texturally late phase near the pluton (<3 km) in rocks with a hornfels texture, and overprints an earlier andalusite ± sillimanite assemblage. Migmatitic rocks and the assemblages K-feldspar + sillimanite and spinel + plagioclase + cordierite also occur only near the pluton, implying that the highest temperatures were associated with intrusion of the granite. Mineral assemblages in hornfelsic rocks and the presence of granitic leucosomes suggest that the stability of muscovite in quartz-bearing rocks was exceeded, and that contact metamorphic P–T conditions were therefore in excess of 700°C at moderate to low pressures (~4 kbar). Contact metamorphism texturally overprinted regional metamorphic assemblages and fabrics. In some zones near the granite, breakdown of garnet to plagioclase + sillimanite and cordierite is replaced (pseudomorphed) by sillimanite + biotite, implying back reaction of cordierite perhaps during the post-migmatite/granite path that accounts for growth of late kyanite (Sepahi, Whitney & Baharifar, Reference Sepahi, Whitney and Baharifar2004).

The P–T diagram (Fig. 6) illustrates the implications of these data. Some rocks (ASKV, ASV, ASVH) record a prograde heating path that is clockwise on a P–T diagram. Other rocks that were buried and heated along apparently similar paths to ASV/H, or at slightly higher pressures (kyanite zone), experienced low-pressure–high-temperature contact metamorphism near the Alvand granite, reaching conditions sufficient to break down muscovite and initiate partial melting in pelitic rocks. The contact metamorphosed rocks, the granite itself, and some rocks that are beyond the contact aureole experienced kyanite zone conditions following intrusion of the granite. Whether these rocks were buried and heated (e.g. during contraction associated with collision and suturing) or whether the P–T path was dominated by cooling is not defined by our observations, garnet zoning and analyses. The P–T diagram, inferred paths and zoning profiles of garnet do not account for the decompression history of the terrane. The breakdown of garnet to plagioclase + sillimanite, dehydration melting and the formation of spinel–plagioclase symplectite on andalusite could occur during decompression or heating; these textures are limited to the contact aureole, so heating is perhaps the more likely explanation for these textures and formation of the symplectites in the metapelitic rocks of the Alvand aureole.

Figure 6. Pressure–temperature diagram (modified after Sepahi, Whitney & Baharifar, Reference Sepahi, Whitney and Baharifar2004) and representative picture of corona evolution in a P–T diagram. Al2SiO5 triple point from Holdaway (Reference Holdaway1971); other equilibria calculated using database and software of Berman (Reference Berman1991). P–T conditions and paths are shown for selected samples. Dsp – diaspore; Kln – kaolinite; Prl – pyrophyllite. The shaded field labelled ‘And’ (andalusite) shows the conditions at which the andalusite schists and andalusite–quartz veins may have equilibrated. The field labelled ‘St’ (staurolite) was calculated for the staurolite–garnet schist (SschH) using garnet–biotite thermometry, brackets from lower and higher grade rocks, and results of Baharifar (1997). The fields labelled ‘St + Ky’ and ‘St + Sil’ are partially defined by equilibria (e.g. the upper stability of staurolite + quartz), but the maximum pressure for St + Ky is not known.

Based on the phases observed in the coronas and on their development between andalusite and matrix (Fig. 2c, d), (exhaustion of the available SiO2 and/or separation of sillimanite/andalusite from SiO2-rich matrix domains by cordierite resulted in the formation of localized low-silica activity domains and thus triggered the growth of spinel in the rim of andalusite) and according to the quartz-absent KFMASH reaction, a first-order approximation of the prograde melting reaction is: Sil/And + Bt = Crd + Spl + Kfs + melt.

8. Conclusions

The microstructural and thermobarometric study of this work enabled us to form some conclusions regarding the spinel–cordierite–plagioclase symplectites replacing andalusite in the Hamadan pelites in the Alvand aureole:

(1) The presence of melt shows that spinel–cordierite–plagioclase symplectites replacing andalusite development occurred under partial melting conditions.

(2) Estimates of the temperature (700°C) and pressure (~2–4 kbar) for the peak metamorphism show that it occurred under LP/HT conditions; the geothermal gradient obtained for the peak metamorphism is 58°C km−1, which is in the upper range of contact aureole values.

(3) The reaction Sil/And + Bt = Crd + Spl + Kfs + melt is the most common reaction for the development of coronas in the metapelitic rocks of the Alvand aureole.

(4) The P–T diagram and inferred paths do not account for the decompression history of the terrane, so breakdown of garnet to plagioclase + sillimanite, dehydration melting and the formation of spinel–plagioclase symplectite are limited to the contact aureole; therefore, heating is perhaps the more likely explanation for formation of the symplectites in the metapelitic rocks of the Alvand aureole.

Acknowledgements

I thank Dr Baharifar for his generous help with electron microprobe analysis. The author is deeply grateful to Dr Mark Allen (Editor of Geological Magazine) and two anonymous reviewers for constructive comments.

References

Agard, P., Omrani, J., Joliver, L. & Mouthereau, F. 2005. Convergence history across Zagros (Iran): constraints from collisional and earlier deformation. International Journal of Earth Sciences 94, 401–19.CrossRefGoogle Scholar
Ahmadi, A. A., Esmaeily, D., Valizadeh, M. V. & Rahimpour-Bonab, H. 2007. Petrology and geochemistry of the granitoid complex of Boroujerd, Sanandaj-Sirjan Zone, Western Iran. Journal of Asian Earth Sciences 29, 859–77.Google Scholar
Alavi, M. 1994. Tectonics of the Zagros Orogenic belt of Iran: new data and interpretations. Tectonophysics 229, 211–38.CrossRefGoogle Scholar
Arvin, M., Pan, Y., Dargahi, S., Malekzadeh, A. & Babaei, A. 2007. Petrochemistry of the Siah-Kuh granitoid stock southwest of Kerman, Iran: implication for initiation of Neotethys subduction. Journal of Asian Earth Sciences 30, 474–89.CrossRefGoogle Scholar
Baharifar, A. A., Moinevaziri, H., Bellon, H. & Pique, A. 2004. The crystalline complexes of Hamadan (Sanandaj–Sirjan zone, western Iran): metasedimentary Mesozoic sequences affected by Late Cretaceous tectono-metamorphic and plutonic events, II. 40K–40Ar dating. Geoscience 336, 1443–52.CrossRefGoogle Scholar
Berberian, F. & Berberian, M. 1981. Tectono-plutonic episodes in Iran. In Zagros Hindukush, Himalaya Geodynamic Evolution (eds Gupta, H. K. & Delany, F. M.), pp. 532. Washington, DC: American Geophysical Union.CrossRefGoogle Scholar
Berberian, M. & King, G. C. 1981. Towards a paleogeography and tectonic evolution of Iran. Canadian Journal of Earth Sciences 18, 210–65.CrossRefGoogle Scholar
Berman, R. 1991. Thermobarometry using multi-equilibrium calculations: a new technique, with petrological applications. Canadian Mineralogist 29, 833–55.Google Scholar
Bucher-Nurminen, K. & Droop, G. 1983. The metamorphic evolution of garnet–cordierite–sillimanite–gneisses of the Gruf-Complex, Eastern Pennine Alps. Contributions to Mineralogy and Petrology 84, 215–27.CrossRefGoogle Scholar
Carson, C. J., Powell, R., Wilson, C. J. L. & Dirks, P. H. G. M. 1997. Partial melting during tectonic exhumation of a granulite terrane: an example from the Larsemann Hills, East Antarctica. Journal of Metamorphic Geology 15, 105–26.CrossRefGoogle Scholar
Cenki, B., Kriegsman, L. M. & Braun, I. 2002. Melt-producing and melt-consuming reactions in anatectic granulites: P–T evolution of the Achankovil cordierite gneisses, South India. Journal of Metamorphic Geology 20, 543–61.CrossRefGoogle Scholar
Cesare, B., Marchesi, C., Hermann, J. & Gomez-Pungnaire, M. T. 2003. Primary melt inclusions in andalusite from anatectic graphitic metapelites: implications for the position of the Al2SiO5 triple point. Journal of Metamorphic Geology 31, 573–6.Google Scholar
Clarke, G. L. & Powell, R. 1991. Decompressional coronas and symplectites in granulites of the Musgrave Complex, central Australia. Journal of Metamorphic Geology 9, 441–50.CrossRefGoogle Scholar
Droop, G. T. R. 1989. Reaction history of garnet-sapphirine granulites and conditions of Archaean high-pressure granulite-facies metamorphism in the Central Limpopo Mobile Belt, Zimbabwe. Journal of Metamorphic Geology 7, 383403.CrossRefGoogle Scholar
Fazlnia, A., Schenk, V., Straaten, F. & Mirmohammadi, M. 2009. Petrology, geochemistry, and geochronology of trondhjemites from the Qori Complex, Neyriz, Iran. Lithos 112, 413–33.CrossRefGoogle Scholar
Ganguly, J. & Saxena, S. 1984. Mixing properties of aluminosilicate garnets: constraints from natural and experimental data and applications to geothermobarometry. American Mineralogist 69, 8897.Google Scholar
Golonka, J. 2004. Plate tectonic evolution of the southern margin of Eurasia in the Mesozoic and Cenozoic. Tectonophysics 38, 235–73.CrossRefGoogle Scholar
Grant, J. A. & Frost, B. R. 1990. Contact metamorphism and partial melting of pelitic rocks in the aureole of the Laramie anorthosite complex, Morton Pass, Wyoming. American Journal of Science 290, 425–72.CrossRefGoogle Scholar
Greenfield, J. E., Clarke, G. L. & White, W. 1998. A sequence of partial melting reactions at Mount Stafford, central Australia. Journal of Metamorphic Geology 16, 363–78.CrossRefGoogle Scholar
Hand, M., Dirks, P. H. G. M., Powell, R. & Buick, I. S. 1992. How well established is isobaric cooling in Proterozoic orogenic belts? An example from the Arunta inlier, central Australia. Geology 20, 649–52.2.3.CO;2>CrossRefGoogle Scholar
Hand, M., Scrimgeour, I., Powell, R., Stüwe, K. & Wilson, C. J. L. 1994. Metapelitic granulites from Jetty Peninsula, east Antarctica: formation during a single event or by polymetamorphism? Journal of Metamorphic Geology 12, 557–73.CrossRefGoogle Scholar
Helffrich, G. & Wood, B. 1989. Subregular model for multicomponent solutions. American Mineralogist 74, 1016–22.Google Scholar
Holdaway, M. J. 1971. Stability of andalusite and the aluminum silicate phase diagram. American Journal of Science 271, 97131.CrossRefGoogle Scholar
Holland, T. J. B. & Powell, R. 1998. An internally consistent thermodynamic data set for phases of petrological interest. Journal of Metamorphic Geology 16, 309–43.CrossRefGoogle Scholar
Ings, S. J. & Owen, J. V. 2002. Decompressional reaction textures formed by isobaric heating: an example from the thermal aureole of the Taylor Brook Gabbro Complex, western Newfoundland. Mineralogical Magazine 66, 941–51.CrossRefGoogle Scholar
Johnson, T., Brown, M., Gibson, R. & Wing, B. 2004. Spinel–cordierite symplectites replacing andalusite: evidence for melt-assisted diapirism in the Bushveld Complex, South Africa. Journal of Metamorphic Geology 22, 529–45.CrossRefGoogle Scholar
Khalaji, A. A., Esmaeily, D., Valizadeh, M. V. & Rahimpour-Bonab, H. 2007. Petrology and geochemistry of the granitoid complex of Boroujerd, Sanandaj-Sirjan Zone, Western Iran. Journal of Asian Earth Sciences 29, 859–77.CrossRefGoogle Scholar
Koziol, A. M. & Newton, R. C. 1988. Redetermination of the anorthite breakdown reaction and improvement of the plagioclase–garnet–aluminosilicate–quartz geobarometer. American Mineralogist 73, 216–23.Google Scholar
Kretz, R. 1983. Symbols for rock-forming minerals. American Mineralogy 68, 277–9Google Scholar
Marmo, B. A., Clarke, G. L. & Powell, R. 2002. Fractionation of bulk rock composition due to porphyroblast growth; effects on eclogite facies mineral equilibria, Pam Peninsula, New Caledonia. Journal of Metamorphic Geology 20, 151–65.CrossRefGoogle Scholar
McDade, P. & Harley, S. L. 2001. A petrogenetic grid for aluminous granulite facies metapelites in the KFMASH system. Journal of Metamorphic Geology 19, 4559.CrossRefGoogle Scholar
Mezger, J. E., Chacko, T. & Erdmer, P. 2001. Metamorphism along a late Mesozoic accretionary continental margin: a case study from the northern Coast Belt of the North American Cordillera. Journal of Metamorphic Geology 19, 121–38.CrossRefGoogle Scholar
Mohajjel, M. & Fergusson, C. L. 2000. Dextral transpression in Late Cretaceous continental collision, Sandandaj–Sirjan Zone, western Iran. Journal of Structural Geology 22, 1125–39.CrossRefGoogle Scholar
Mohajjel, M., Fergusson, C. L. & Sahandi, M. R. 2003. Cretaceous–Tertiary convergence and continental collision, Sanandaj–Sirjan Zone, Western Iran. Journal of Asian Earth Sciences 21, 397412.CrossRefGoogle Scholar
Norlander, B. H., Whitney, D. L., Teyssier, C. & Vanderhaeghe, O. 2002. Partial melting and decompression of the Thor-Odin Dome, Shuswap metamorphic core complex, Canadian Cordillera. Lithos 61, 103–25.CrossRefGoogle Scholar
Omrani, J., Agard, Ph., Whitechurch, H., Benoit, M., Prouteau, G. & Jolivet, L. 2008. Arc magmatism and subduction history beneath the Zagros Mountains, Iran: a new report of adakites and geodynamic consequences. Lithos 106, 380–98.CrossRefGoogle Scholar
Passchier, C. W. & Trouw, R. A. J. 1996. Microtectonics. Heidelberg: Springer-Verlag, 289 pp.Google Scholar
Pitra, P. & De Wall, S. A. 2001. High-temperature, low-pressure metamorphism and development of prograde symplectites, Marble Hall Fragment, Bushveld Complex (South Africa). Journal of Metamorphic Geology 19, 311–25.Google Scholar
Riesco, M., Stüwe, K., Reche, J. & Martinez, J. 2004. Silica depleted melting of pelites. Petrogenetic grid and application to the Susqueda Aureole, Spain. Journal of Metamorphic Geology 22, 475–94.CrossRefGoogle Scholar
Robinson, P. R., Hollocher, K. T., Tracy, R. J. & Dietsch, C. W. 1982. High grade Acadian regional metamorphism in south-central Massachusetts. In NEIGC 74th Annual Meeting of the State Geological and Natural History Survey of Connecticut, guidebook for fieldtrips in Connecticut and South-Central Massachusetts (eds Joester, R. A. & Quarrier, S. S.), pp. 289340. Storrs: The University of Connecticut.Google Scholar
Saki, A. & Baharifar, A. A. 2010. Common melting reactions and their characteristics in the Alvand aureole metapelites, Hamadan. Iranian Jounal of Geosciences (in press).Google Scholar
Saki, A. & Pourkaseb, H. 2010. P–T Conditions and fluid composition for wollastonite-clinopyroxene-garnet-bearing calc-silicate rocks from contact aureole of the Alvand batholite, Hamadan. Shahid Chamran University Journal of Sciences (in press).Google Scholar
Sears, J. W., George, G. M. S. & Winne, J. C. 2005. Continental rift systems and anorogenic magmatism. Lithos 80, 147–54.CrossRefGoogle Scholar
Sederholm, J. J. 1967. Selected Works: Granite and Migmatites. Edinburgh: Oliver and Boyd.Google Scholar
Sepahi, A. A., Whitney, D. L. & Baharifar, A. A. 2004. Petrogenesis of And–Ky–Sil veins and host rocks, Sanandaj-Sirjan metamorphic belt, Hamadan, Iran. Journal of Metamorphic Geology 22 (2), 119–34.CrossRefGoogle Scholar
Shahabpour, J. 2005. Tectonic evolution of the orogenic belt in the region located between Kerman and Neyriz. Journal of Asian Earth Sciences 24, 405–17.CrossRefGoogle Scholar
Shahbazi, H., Siebel, W., Pourmoafee, M., Ghorbani, M., Sepahi, A. A., ShangC. K., Vousoughi C. K., Vousoughi. & Abedini, M. 2010. Geochemistry and U–Pb zircon geochronology of the Alvand plutonic complex in Sanandaj–Sirjan Zone (Iran): new evidence for Jurassic magmatism. Journal of Asian Earth Sciences, doi:10.1016/j.jseaes.2010.04.014, in press.CrossRefGoogle Scholar
Spear, F. S. 1993. Metamorphic Phase Equilibria and Pressure–Temperature–Time Paths. Monograph 1. Washington, DC: Mineralogical Society of America.Google Scholar
Stüwe, K. 1997. Effective bulk composition changes due to cooling: a model predicting complexities in retrograde reaction textures. Contributions to Mineralogy and Petrology 129, 4352.CrossRefGoogle Scholar
Topuz, G., Altherra, R., Kalta, A., Satirb, M., Wernera, O. & Schwarza, W. H. 2004. Aluminous granulites from the Pulur complex, NE Turkey: a case of partial melting, efficient melt extraction and crystallisation. Lithos 72, 183207.CrossRefGoogle Scholar
Vernon, R. H. 1996. Problems with inferring P–T–t paths in low-P granulite facies rocks. Journal of Metamorphic Geology 14, 143–53.CrossRefGoogle Scholar
Waters, D. J. 1991. Hercynite-quartz granulites: phase relations and implications for crustal processes. European Journal of Mineralogy 3, 367–86.CrossRefGoogle Scholar
White, R. W., Powell, R. & Clarke, G. L. 2002. The interpretation of reaction textures in Fe-rich metapelitic granulites of the Musgrave Block, central Australia: constraints from mineral equilibria calculations in the system K2O–FeO–MgOAl2O3–SiO2–H2O–TiO2–Fe2O3. Journal of Metamorphic Geology 20, 4155.CrossRefGoogle Scholar
White, R. W., Powell, R. & Clarke, G. L. 2003. Prograde metamorphic assemblage evolution during partial melting of metasedimentary rocks at low pressures: migmatites from Mt Stafford, Central Australia. Journal of Petrology 44, 1937–60.CrossRefGoogle Scholar
White, R. W., Powell, R. & Holland, T. J. B. 2001. Calculation of partial melting equilibria in the system Na2O–CaO–K2O–FeO–MgO–Al2O3–SiO2–H2O (NCKFMASH). Journal of Metamorphic Geology 19, 139–53.CrossRefGoogle Scholar
Whittington, A., Harris, N. & Baker, J. 1998. Low-pressure crustal anatexis. In What Drives Metamorphism and Metamorphic Reactions? (eds Treolar, P. J. & O'Brien, P. J.), pp. 183–98. Geological Society of London, Special Publication no. 138.Google Scholar
Figure 0

Figure 1. (a) Location of the Sanandaj–Sirjan belt in Iran. (b) Geological map of the area (modified from Baharifar et al. 2004), with sample locations and stratigraphic column of the area. (c) Generalized geological map of the study area near Hamadan. Isograds (dashed lines) and metamorphic zones are shown (modified from Sepahi, Whitney & Baharifar, 2004). Some faults are shown (heavy lines within unit 4); contacts between zones 4a–4b, 4a–5, 2–5, and possibly 4–7 are faults. Zones: 1 – biotite zone; 2 – garnet zone; 3 – andalusite zone (±fibrolite); 4 – staurolite zone (a = staurolite + andalusite ± fibrolite; b = staurolite + garnet, no andalusite); 5 – sillimanite–andalusite zone; 6 – sillimanite–Kfs zone; 7 – cordierite zone (a = cordierite with no new growth of Al2SiO5; b = cordierite + neoblastic andalusite + sillimanite); 8 – cordierite Kfs zone; 9 – spotted schist; 10 – spinel + orthopyroxene zone (migmatites).

Figure 1

Figure 2. Field photographs of migmatites. A colour version of this figure is available at http://www.cambridge.org/journals/geo. (a, b) Migmatite in which smaller leucosomes in (a) and (b) appear to be crystallized from melt derived locally, whereas the larger leucosome likely records evidence of melt moving through the system. (c, d) Plane-polarized light view of a Spl + Crd + And corona in the mesosome. (e) Sil + Crd + Spl in the mesosome. (f) The garnet as inclusion in the cordierite. (g) Feldspar–plagioclase–quartz patches (Lu), with the melt apparently segregated into mesoscopic patches or veins. (h) Euhedral shape of the plagioclase phenocrysts against intersertal quartz in (g).

Figure 2

Table 1. Representative mineral analyses of garnet and cordierite

Figure 3

Figure 3. Zoning profile from core to rim through the garnet in some samples of migmatite (mesosome; contact metamorphism), which shows a flat profile attributed to garnet crystallization in these rocks.

Figure 4

Figure 4. Zoning profile of the garnet in some samples from the regional metapelitic rocks, which shows prograde metamorphism (normal zoning) attributed to garnet crystallization in different regional metamorphic zones.

Figure 5

Table 2. Representative mineral analyses of biotite, spinel and ilmenite

Figure 6

Table 3. Representative mineral analyses of feldspar

Figure 7

Figure 5. Representative results of multiple equilibria calculation by THERMOCALC (3.2) program; a pressure of ~2 kbar and temperature of ~750°C can be considered for peak metamorphism.

Figure 8

Table 4. Activity values used for calculating equilibrium reactions by THERMOCALC

Figure 9

Figure 6. Pressure–temperature diagram (modified after Sepahi, Whitney & Baharifar, 2004) and representative picture of corona evolution in a P–T diagram. Al2SiO5 triple point from Holdaway (1971); other equilibria calculated using database and software of Berman (1991). P–T conditions and paths are shown for selected samples. Dsp – diaspore; Kln – kaolinite; Prl – pyrophyllite. The shaded field labelled ‘And’ (andalusite) shows the conditions at which the andalusite schists and andalusite–quartz veins may have equilibrated. The field labelled ‘St’ (staurolite) was calculated for the staurolite–garnet schist (SschH) using garnet–biotite thermometry, brackets from lower and higher grade rocks, and results of Baharifar (1997). The fields labelled ‘St + Ky’ and ‘St + Sil’ are partially defined by equilibria (e.g. the upper stability of staurolite + quartz), but the maximum pressure for St + Ky is not known.

Supplementary material: Image

Saki supplementary material

Colour figure.tif

Download Saki supplementary material(Image)
Image 10.8 MB