Abstract Granitic rocks represent a ubiquitous component of upper continental crust but their origin remains highly controversial. This controversy stems from the fact that the granites may result from fractionation of mantle-derived basaltic magmas or partial melting of different crustal protoliths at contrasting pressure–temperature conditions, either water-fluxed or fluid-absent. Consequently, many different mechanisms have been proposed to explain the compositional variability of granites ranging from whole igneous suites down to mineral scale. This Special Publication presents an overview of the state of the art and envisages future avenues towards a better understanding of granite petrogenesis.

Abstract Granites ( sensu lato ) come in many types and flavours, defining distinct magmatic series/suites/types. A good classification not only gives generally accepted and understandable names to similar rocks but also links the bulk chemical composition to the stoichiometry of the constituent minerals and, potentially, also to the likely source, magmatic evolution and tectonic setting. The ‘ideal’ granitoid classification should be based on chemical criteria amenable to an objective treatment. Statistical analysis helps to identify the most discriminant variables. The key properties are (1) acidity/maficity, (2) alkalinity (balance of Na + K v. Ca), (3) aluminosity (balance of Al v. Ca, Na and K), (4) Fe/Mg balance and (5) Na/K balance and K contents at the given SiO 2 level. These are used by successful classifications, e.g. the I/S dichotomy is based mainly on aluminosity, while the Frost et al. (2001 ; ‘A geochemical classification for granitic rocks', Journal of Petrology , 42 , 2033–2048, https://doi.org/10.1093/petrology/42.11.2033 ) classification includes all but Na/K. Even though it is commonplace to use weight percentages of oxides, we suggest that a better strategy is to employ simple atomic parameters (e.g. millications-based) that can be directly linked to modal proportions and compositions/crystal structure of individual rock-forming minerals. This facilitates a petrological interpretation, which, in turn, can be related to petrogenesis and, ultimately, to likely tectonic setting(s).

Abstract The classical S–I–A-type granites from the Lachlan Orogen, SE Australia, formed as a tectonic end-member of the accretionary orogenic spectrum, the Paleozoic Tasmanides. The sequence of S- to I- to A-type granite is repeated at least three times. All the granites are syn-extensional, formed in a dominantly back-arc setting behind a single, stepwise-retreating arc system between 530 and 230 Ma. Peralkaline granites are rare. Systematic S–I–A progressions indicate the progressive dilution of an old crustal component as magmatism evolved from arc (S-type) to proximal back-arc (I-type) to distal back-arc (A-type) magmatism. The alkaline and peralkaline A-type Younger granites of Nigeria were generally hotter and drier than the Lachlan A-type granites and were emplaced into an anhydrous Precambrian basement during intermittent intracontinental rifting. This geodynamic environment contrasts with the distal back-arc setting of the Lachlan A-type granites, where magmatism migrated rapidly across the orogen. Tectonic discrimination diagrams are inappropriate for the Lachlan granites, placing them in the wrong settings. Only the peralkaline Narraburra suite of the Lachlan Orogen fits the genuine ‘within-plate’ setting of the Nigerian A-type granites. Such discrimination diagrams require re-evaluation in the light of an improved modern understanding of tectonic processes, particularly the role of extensional tectonism and its geodynamic drivers.

Abstract The origin of large I-type batholiths remains a disputed topic. One model states that I-type granites form by partial melting of older crustal lithologies (amphibolites or intermediate igneous rocks). In another view, granites are trapped rhyolitic liquids occurring at the end of fractionation trends defining a basalt–andesite–dacite–rhyolite series. This paper explores the thermal implications of both scenarios, using a heat balance model that abstracts the heat production and consumption during crustal melting. Heat is consumed by melting and by losses through the surface (conductive or advective, as a result of eruption). It is supplied as a basal conductive heat flux, as internal heat production or as advective heat carried by an influx of hot basalt into the crust. Using this abstract approach, it is possible to explore the role different parameters play in the balance of granites formed by differentiation of basalts or by crustal melting. Two end-member situations appear equally favourable to generating large volumes of granites: (1) short-lived environments dominated by high basaltic flux, where granites result mostly from basalt differentiation; and (2) long-lived systems with no or minimal basalt flux, with granites resulting chiefly from crustal melting.

Abstract New laboratory experiments using granulite xenoliths support a dual origin for I-type granites as primary and secondary. Primary I-type granites represent fractionated liquids from intermediate magma systems of broadly andesitic composition. Fluid-fluxed melting of igneous rocks that resided in the continental crust generates secondary I-type granites. The former are directly related to subduction, with Cordilleran batholiths as the most characteristic examples. Experiments with lower crust granulite sources, in the presence of water, show that amphibole is formed by a water-fluxed peritectic rehydration melting reaction . Entrainment of only 10% of restites composed of amphibole, pyroxene, plagioclase and magnetite, is sufficient to account for discrepancies in aluminium saturation index and maficity in secondary I-type granites. Lower crust granulite xenoliths, attached to a sanukitoid containing 6 wt% water, have been used in two-layer capsules to test fluid-fluxed melting reactions as the origin of secondary I-type granites. It is proposed that sanukitoid magmas act as water donors that trigger extensive melting of the lower crust, giving rise to granodioritic liquids. Because primary granites are related to coeval subduction, and secondary ones are crustal melts from older subduction-related rocks, the distinction between both I-types is essential in tectonic reconstructions of ancient orogenic belts.

Abstract Peraluminous granites and trondhjemites make up small plutonic bodies intruded into high-grade paragneisses in the Peloritani Mountains, marking the beginning of late Variscan granitoid magmatism in southernmost Italy. The granites range from low-Ca monzogranites to alkali feldspar granites, while the trondhjemites vary from trondhjemites s.s . to low-Ca trondhjemites. Relatively high radiogenic ( 87 Sr/ 86 Sr) i ratios (mostly from 0.7073 to 0.7125) and negative ε Nd values (mostly from −5.66 to −8.73) point to crustal sources for all the granitoids. Major and trace element compositions indicate an absence of genetic relationships between the trondhjemites s.s . and the low-Ca granitoids, but possible relationships between the low-Ca trondhjemites and the granites. All of the studied granitoids have near-pure melts compositions, consistent with H 2 O-fluxed and dehydration melting of metasediments for the trondhjemites and the granites, respectively. However, the unusual compositions of the low-Ca trondhjemites and microstructural evidence in these rocks for pervasive subsolidus replacement of magmatic feldspars by secondary sodic plagioclase indicate that they were derived instead from metasomatic alteration of the granites. Thus, water may be involved in the production of trondhjemites in two different ways, driving water-fluxed melting in the magma source and driving alkali metasomatism at the sites of granite emplacement in the upper crust.

Abstract Leucogranites are a characteristic feature of collisional orogens. Their generation is intimately related to crustal thickening and the active deformation and metamorphism of metapelites. Data from Proterozoic to present day orogenic belts show that collisional leucogranites (CLGs) are peraluminous, with muscovite, biotite and tourmaline as characteristic minerals. Isotopic ratios uniquely identify the metapelitic sequences in which CLGs occur as sources. Organic material in pelitic sources results in f O 2 in CLGs that is usually below the fayalite–magnetite–quartz buffer. Most CLGs form under vapour-poor conditions with melting involving a peritectic breakdown of muscovite. The low concentrations of Mg, Fe and Ti that characterize CLGs are largely related to biotite–melt equilibria in the source rocks. Concentrations of Zr, Th and rare earth elements are lower than expected from zircon and monazite saturation models because these minerals often remain enclosed in residual biotite during melting. Melting involving muscovite may limit the temperatures achieved in the source regions. A lack of nearby mantle heat sources in thick collisional orogens has led to thermal models for the generation of CLGs that involve flux melting, or large amounts of radiogenic heat generation, or decompression melting or shear heating, the last one emphasizing the link of leucogranites and their sources to crustal-scale shear zone systems.

Abstract Modern quantitative phase equilibria modelling allows the calculation of the stable phase assemblage of a rock system given its pressure, temperature and bulk composition. A new software tool (Rcrust) has been developed that allows the modelling of points in pressure–temperature–bulk composition space in which bulk compositional changes can be passed from point to point as the system evolves. This new methodology enables quantitative process-oriented investigation of the evolution of rocks. Procedures are outlined here for using this tool to model: (1) the control of the water content of a subsolidus system based on available pore space; (2) the triggering of melt loss events when a critical melt volume threshold is exceeded, while allowing a portion of melt retention; (3) the entrainment of crystals during segregation and ascent of granitic magmas from its source; (4) the modification of the composition of granite magmas owing to fractional crystallization; and (5) the progressive availability (through dissolution) of slow diffusing species and their control of the effective bulk composition of a system. These cases collectively illustrate thermodynamically constrained methods for modelling systems that involve mass transfer.

Abstract The ability of Rcrust software to conduct path-dependent phase equilibrium modelling with automated changing bulk compositions allows for a phase equilibrium approach to investigate an array of source controls for their effect on the bulk compositions of melts produced by sequential melting events. The following source controls of the rock system are considered: (1) initial magnesium and iron content; (2) initial sodium and calcium content; (3) pressure–temperature path followed by the system; and (4) threshold by which melt extractions in the system are triggered. These source controls are investigated in a water-restricted system and a water-in-excess system. The permutation of these cases resulted in 128 different modelled pressure–temperature bulk composition paths investigating the melting of an average pelite composition. The resultant melt compositions are compared to that of a natural granite dataset and provide a good fit for the incompatible elements Na 2 O and K 2 O with the allowance that granites most likely form as magmas consisting of melt and ferromagnesian-rich crystals. The fluid state of the system is shown to have the strongest control on melt compositions, with the pressure–temperature path having subordinate control on the volume and composition of melts produced.

Abstract The Buddusò Pluton in NE Sardinia (Italy) is a normally zoned intrusion composed of three units with chemical composition ranging from hornblende-bearing tonalites (SiO 2 ∼ 65 wt%) to leucocratic monzogranites (SiO 2 ∼ 76 wt%). Zircon crystals in the pluton are dated at 292.2 ± 0.7 Ma and have ε Hf values ranging from −4 to −8, with no systematic differences observed between the units. The pluton, which is isotopically homogeneous at the whole-rock scale in terms of Sr and Nd isotopes, shows textural evidence indicating local crystal–melt segregation. In this paper, we have implemented a novel approach based on path-dependent phase-equilibria modelling to test the hypothesis that the internal chemical variability of the pluton was generated by crystallization differentiation of a homogeneous parental magma. Our modelling indicates that this hypothesis is valid if the mechanism by which this occurs is compaction in a rheologically locked crystal-rich magma and if the separation occurs at 0.3 GPa from a tonalitic magma with water content >2 wt%. Finally, a subset of the magmatic enclaves in the pluton are considered to be autoliths, formed by the disruption of the compacted crystal mush and interaction between these cumulates and the felsic melt.

Abstract Whole-rock geochemistry represents a powerful tool in deciphering petrogenesis of magmatic suites, including granitoids, which can be used to formulate and test hypotheses qualitatively and often also quantitatively. Typically, it can rule out impossible/improbable scenarios and further constrain the process inferred on geological and petrological grounds. With the current explosion of high-precision data, both newly acquired and retrieved from extensive databases, the whole-rock geochemistry-based petrogenetic modelling of igneous rocks will gain further importance. Especially promising is its combination with thermodynamic modelling into a single, coherent and comprehensive software, using the R and Python languages.

Granites ( sensu lato ) represent the dominant rock-type forming the upper–middle continental crust but their origin remains a matter of long-standing controversy. The granites may result from fractionation of mantle-derived basaltic magmas, or partial melting of different crustal protoliths at contrasting P–T conditions, either water-fluxed or fluid-absent. Consequently, many different mechanisms have been proposed to explain the compositional variability of granites ranging from whole igneous suites down to mineral scale. This book presents an overview of the state of the art, and envisages future avenues towards a better understanding of granite petrogenesis. The volume focuses on the following topics: compositional variability of granitic rocks generated in contrasting geodynamic settings during the Proterozoic to Phanerozoic Periods; main permissible mechanisms producing subduction-related granites; crustal anatexis of different protoliths and the role of water in granite petrogenesis; and new theoretical and analytical tools available for modelling whole-rock geochemistry in order to decipher the sources and evolution of granitic suites.

Combined field and geophysical data show that plutons from the Bega Batholith are elongate, meridional, wedge-shaped bodies which intruded during a period of regional east-west extension in the Palaeozoic eastern Lachlan orogen, eastern Australia. Plutons within the core of the batholith have intruded coeval, syn-rift sediments and co-magmatic volcanics. The batholith is bound by high-temperature, dip-slip faults, and contains several major NE- trending transtensional faults which were active during batholith construction. In the central part of the batholith, the Kameruka pluton is an asymmetric, eastward-thickening, wedge-shaped body with the base exposed as the western contact, which is characterised by abundant, shallow-dipping schlieren migmatites which contain recumbent folds and extensional shear bands. A shallow (<30°), east-dipping, primary magmatic layering in the Kameruka pluton steepens progressively westward, where it becomes conformable to the east-dipping basal migmatites. The systematic steepening of the layering is comparable to sedimentary units formed during floor depression in syn-rift settings. The present authors suggest that the wedge-shaped plutons of the Bega Batholith are the deeper, plutonic expression of a hot, active rift. The batholith was fed and sustained by injection of magma through sub-vertical dykes. Displacement along syn-magmatic, NE- trending faults suggests up to 25 km of arc-perpendicular extension during batholith construction. The inferred tectonic setting for batholith emplacement is a continental back-arc, where modern half-extension rates of 20–40 mm yr −1 are not unusual, and are sufficient to emplace the entire batholith in ∼ 1Ma. This structural model provides a mechanism for the emplacement of some wedge-shaped plutons and is one solution to the ‘room problem’ of batholith emplacement.

Abstract A global three-dimensional chemistry-transport model has been applied to study the tropospheric sulphur cycle, and in particular the volcanic component. The model is in general agreement with previous studies of the global S budget. We find that volcanic emissions constitute 10% of the present-day global SO 2 source to the atmosphere, but form 26% of the SO 2 burden, and 14% of the sulphate aerosol burden. Two previous modelling studies suggested that the volcanic fraction of sulphate was 18% and 35%, from sources representing 7% and 14%, respectively, of the global total SO 2 emission. The results are dependent upon various assumptions about volcanic emissions (magnitude, geographical location, altitude), the global distribution of oxidants, and the physical processes of dry and wet deposition. Because of this dependence upon poorly constrained parameters, it is unclear which modelling study is closest to the truth.

Field, petrographic and geochemical evidence from the K-feldspar megacrystic Kameruka pluton, Lachlan Fold Belt, southeastern Australia, suggests that complex, multi-component, mafic microgranular enclaves (MME) are produced by two-stage hybridisation processes. Stage 1 mixing occurs in composite dykes below the pluton, as mafic and silicic melts ascend through shared conduits. Pillows formed in these conduits are homogeneous, fine-to medium-grained stage 1 MME, which typically range from basaltic to granitic compositions that plot as a sublinear array on Marker diagrams. Stage 2 hybridisation occurs in the magma chamber when the composite dykes mix with the resident magma as synplutonic dykes. The stage 2 hybrids also form linear chemical arrays and range from basaltic to granodioritic compositions, the latter resembling the more mafic phases of the pluton. Stage 2 MME are distinguished from stage 1 types by the presence of K-feldspar xenocrysts and a more heterogeneous nature: they commonly contain stage 1 enclaves. Subsequent disaggregation and dispersal of stage 2 hybrid synplutonic dykes within the magma chamber produces a diverse array of multi-component MME. Field evidencefor conduit mixing is consistent with published analogue experimental studies, which show that hybrid thermo-mechanical boundary layers (TMBL) develop between mafic and silicic liquids in conduits. A mechanical mixing model is developed, suggesting that the TMBL expands and interacts with the adjacent contrasting melts during flow, producing an increasing compositional range of hybrids with time that are mafic in the axial zone, grading to felsic in the peripheral zones in the conduit. Declining flow rates in the dyke and cooling of the TMBL zones produce a pillowing sequence progressing from mafic to felsic, which explains the general observation of more MME in more silicic hosts. The property of granitic magmas to undergo transient brittle failure in seismic regimes allows analogies with fractured solids to be drawn. The fracture network in granitic magmas consists of through-going ‘backbone’ mafic and silicic ± composite dykes, and smaller ‘dangling’ granitic dykes locally generated in the magma chamber. Stage 1 hybrids form in composite backbone dykes and stage 2 hybrids form where they intersect dangling dykes in the magma chamber. With subsequent shear stress recovery, the host magma chamber reverts to a visco-plastic material capable of flow, resulting in disaggregation and dispersal of these complex, hybrid synplutonic dykes, and a vast array of double and multicomponent enclaves potentially develop in the pluton.

Abstract Granulite facies anatexis ( T ≈ 900 °C) in the Wuluma Hills region of the Arunta Inlier was synchronous with deformation. During D3 contractional deformation strain was partitioned into S3 shear zones, which alternate with lower strain domains containing F3 fold hinges. Subsequent D4 deformation was minor and in part extensional. Leucosomes in the S3 shear zones are principally veins oriented parallel, or subparallel, to the pervasive S3 foliation. Leucosomes in the F3 hinge domains are more complex, and occur parallel to anisotropy due to lithological layering, the pre-existing S1/2 foliation, S3 and fold axial planes (F3 and F4). Some leucosomes (generally high Na 2 O, low K 2 O and Rb/Sr) record melt migration paths, and other sites of melt accumulation. All the migmatites are residual and lost melt when deformation forced melt from matrix grain boundaries, through a network of small lensoid channelways to accumulation sites in fold hinges, there larger batches of magma developed. Leucosomes in accumulation sites develop a schlieric or diatexitic appearance because inflowing melt eroded the host rocks. Later increments of D3 contractional strain overpressured the accumulated granitic magma and it migrated again to other (more stable) low pressure sites through veins generally oriented parallel to S3. Magma/melt movement stopped when the solidus was reached, or the magma reached a structurally stable site (e.g. pluton).

The paradox of Lachlan Fold Belt (LFB) granitoids is that although contrasted chemical types (S- and I-types) imply melting of distinct crustal sources, the simple Nd–Sr–Pb–O isotopic arrays indicate a continuum, suggesting mixing of magmatic components. The paradox is resolved by the recognition that the previously inferred, isotopically primitive end-member is itself a crust–mantle mix, so that three general source components, mantle, lower crust and middle crust, comprise the granitoids. Based on Nd isotopic evidence, mantle-derived basaltic magmas melted and mixed with Neoproterozoic–Cambrian, arc–backarc-type material to produce primitive I-type, parental granitoid magmas in the lower–middle crust. Ordovician metasediment, locally underthrust to mid-crustal levels, was remobilised under the elevated geotherms and is most clearly recognised as diatexite in the Cooma complex, but it also exists as gneissic enclaves in S-type granites. The diatexite mixed with the hybrid I-type magmas to produce the parental S-type magmas. Unique parent magma compositions of individual granite suites reflect variations within any or all of the three major source components, or between the mixing proportions. For example, chemical tie-lines between Cooma diatexite and mafic I-type Jindabyne suite magma encompass almost all mafic S-type granites of the vast Bullenbalong supersuite, consistently in the proportion Jindabyne: Cooma, 30:70. The modelling shows that LFB S-type magmas are heavily contaminated I-type magmas, produced by large-scale mixing of hot I-type material with lower temperature diatexite in the middle crust. The model implies a genetic link between migmatite and pluton-scale, crustally derived (S-type) granites. Given the chemical and isotopic contrasts of the crustally derived source components, and their typically unequal proportions in the magmas, it is not surprising that the LFB granitoids are so distinctive and have been categorised as S- and I-type. The sublinear chemical trends of the granitoid suites are considered to be secondary effects associated with crystal fractionation of unique parental magmas that were formed by three-component mixing. The model obviates the necessity for multiple underplating events and Proterozoic continental basement, in accordance with the observed tectono-stratigraphy of the Lachlan Fold Belt.