The Palu Metamorphic Complex, NW Sulawesi, Indonesia: Origin and evolution of a young metamorphic terrane with links to Gondwana and Sundaland

The Palu Metamorphic Complex (PMC) is exposed in a late Cenozoic orogenic belt in NW Sulawesi, Indonesia. It is a composite terrane comprising a gneiss unit of Gondwana origin, a schist unit composed of meta-sediments deposited along the SE Sundaland margin in the Late Cretaceous and Early Tertiary, and one or more slivers of amphibolite with oceanic crust characteristics. The gneiss unit forms part of the West Sulawesi block underlying the northern and central sections of the Western Sulawesi Province. The presence of Late Triassic granitoids and recycled Proterozoic zircons in this unit combined with its isotopic signature suggests that the West Sulawesi block has its origin in the New Guinea margin from which it rifted in the late Mesozoic. It docked with Sundaland sometime during the Late Cretaceous. U–Th–Pb dating results for monazite suggest that another continental fragment may have collided with the Sundaland margin in the earliest Miocene. High-pressure (HP) and ultrahigh-pressure (UHP) rocks (granulite, peridotite, eclogite) are found as tectonic slices within the PMC, mostly along the Palu–Koro Fault Zone, a major strike-slip fault that cuts the complex. Mineralogical and textural features suggest that some of these rocks resided at depths of 60–120 km during a part of their histories. Thermochronological data (U–Th–Pb zircon and 40Ar/39Ar) from the metamorphic rocks indicate a latest Miocene to mid-Pliocene metamorphic event, which was accompanied by widespread granitoid magmatism and took place in an extensional tectonic setting. It caused recrystallization of, and new overgrowths on, pre-existing zircon crystals, and produced andalusite–cordierite–sillimanite–staurolite assemblages in pelitic protoliths, indicating HT–LP (Buchan-type) metamorphism. The PMC was exhumed as a core complex at moderate rates (c. 0.7–1.0 mm/yr) accompanied by rapid cooling in the Plio-Pleistocene. Some of the UHP rocks were transported to the surface at significantly higher rates (⩾16 mm/yr). The results of our study do not support recent plate tectonic reconstructions that propose a NW Australia margin origin for the West Sulawesi block (e.g. Hall et al., 2009).

In search of a hidden long-term isolated sub-chondritic 142Nd/144Nd reservoir in the deep mantle : implications for the Nd isotope systematics of the Earth

Here we search for evidence of the existence of a sub-chondritic 142Nd/144Nd reservoir that balances the Nd isotope chemistry of the Earth relative to chondrites. If present, it may reside in the source region of deeply sourced mantle plume material. We suggest that lavas from Hawai’i with coupled elevations in 186Os/188Os and 187Os/188Os, from Iceland that represent mixing of upper mantle and lower mantle components, and from Gough with sub-chondritic 143Nd/144Nd and high 207Pb/206Pb, are favorable samples that could reflect mantle sources that have interacted with an Early-Enriched Reservoir (EER) with sub-chondritic 142Nd/144Nd. High-precision Nd isotope analyses of basalts from Hawai’i, Iceland and Gough demonstrate no discernable 142Nd/144Nd deviation from terrestrial standards. These data are consistent with previous high-precision Nd isotope analysis of recent mantle-derived samples and demonstrate that no mantle-derived material to date provides evidence for the existence of an EER in the mantle. We then evaluate mass balance in the Earth with respect to both 142Nd/144Nd and 143Nd/144Nd. The Nd isotope systematics of EERs are modeled for different sizes and timing of formation relative to ε143Nd estimates of the reservoirs in the μ142Nd = 0 Earth, where μ142Nd is ((measured 142Nd/144Nd/terrestrial standard 142Nd/144Nd)−1 * 10−6) and the μ142Nd = 0 Earth is the proportion of the silicate Earth with 142Nd/144Nd indistinguishable from the terrestrial standard. The models indicate that it is not possible to balance the Earth with respect to both 142Nd/144Nd and 143Nd/144Nd unless the μ142Nd = 0 Earth has a ε143Nd within error of the present-day Depleted Mid-ocean ridge basalt Mantle source (DMM). The 4567 Myr age 142Nd–143Nd isochron for the Earth intersects μ142Nd = 0 at ε143Nd of +8 ± 2 providing a minimum ε143Nd for the μ142Nd = 0 Earth. The high ε143Nd of the μ142Nd = 0 Earth is confirmed by the Nd isotope systematics of Archean mantle-derived rocks that consistently have positive ε143Nd. If the EER formed early after solar system formation (0–70 Ma) continental crust and DMM can be complementary reservoirs with respect to Nd isotopes, with no requirement for significant additional reservoirs. If the EER formed after 70 Ma then the μ142Nd = 0 Earth must have a bulk ε143Nd more radiogenic than DMM and additional high ε143Nd material is required to balance the Nd isotope systematics of the Earth.

Basin analysis and economic geology of the Northern Mount Isa Basin

The Mount Isa Basin is a new concept used to describe the area of Palaeo- to Mesoproterozoic rocks south of the Murphy Inlier and inappropriately described presently as the Mount Isa Inlier. The new basin concept presented in this thesis allows for the characterisation of basin-wide structural deformation, correlation of mineralisation with particular lithostratigraphic and seismic stratigraphic packages, and the recognition of areas with petroleum exploration potential. The northern depositional margin of the Mount Isa Basin is the metamorphic, intrusive and volcanic complex here referred to as the Murphy Inlier (not the "Murphy Tectonic Ridge"). The eastern, southern and western boundaries of the basin are obscured by younger basins (Carpentaria, Eromanga and Georgina Basins). The Murphy Inlier rocks comprise the seismic basement to the Mount Isa Basin sequence. Evidence for the continuity of the Mount Isa Basin with the McArthur Basin to the northwest and the Willyama Block (Basin) at Broken Hill to the south is presented. These areas combined with several other areas of similar age are believed to have comprised the Carpentarian Superbasin (new term). The application of seismic exploration within Authority to Prospect (ATP) 423P at the northern margin of the basin was critical to the recognition and definition of the Mount Isa Basin. The Mount Isa Basin is structurally analogous to the Palaeozoic Arkoma Basin of Illinois and Arkansas in southern USA but, as with all basins it contains unique characteristics, a function of its individual development history. The Mount Isa Basin evolved in a manner similar to many well described, Phanerozoic plate tectonic driven basins. A full Wilson Cycle is recognised and a plate tectonic model proposed. The northern Mount Isa Basin is defined as the Proterozoic basin area northwest of the Mount Gordon Fault. Deposition in the northern Mount Isa Basin began with a rift sequence of volcaniclastic sediments followed by a passive margin drift phase comprising mostly carbonate rocks. Following the rift and drift phases, major north-south compression produced east-west thrusting in the south of the basin inverting the older sequences. This compression produced an asymmetric epi- or intra-cratonic clastic dominated peripheral foreland basin provenanced in the south and thinning markedly to a stable platform area (the Murphy Inlier) in the north. The fmal major deformation comprised east-west compression producing north-south aligned faults that are particularly prominent at Mount Isa. Potential field studies of the northern Mount Isa Basin, principally using magnetic data (and to a lesser extent gravity data, satellite images and aerial photographs) exhibit remarkable correlation with the reflection seismic data. The potential field data contributed significantly to the unravelling of the northern Mount Isa Basin architecture and deformation. Structurally, the Mount Isa Basin consists of three distinct regions. From the north to the south they are the Bowthorn Block, the Riversleigh Fold Zone and the Cloncurry Orogen (new names). The Bowthom Block, which is located between the Elizabeth Creek Thrust Zone and the Murphy Inlier, consists of an asymmetric wedge of volcanic, carbonate and clastic rocks. It ranges from over 10 000 m stratigraphic thickness in the south to less than 2000 min the north. The Bowthorn Block is relatively undeformed: however, it contains a series of reverse faults trending east-west that are interpreted from seismic data to be down-to-the-north normal faults that have been reactivated as thrusts. The Riversleigh Fold Zone is a folded and faulted region south of the Bowthorn Block, comprising much of the area formerly referred to as the Lawn Hill Platform. The Cloncurry Orogen consists of the area and sequences equivalent to the former Mount Isa Orogen. The name Cloncurry Orogen clearly distinguishes this area from the wider concept of the Mount Isa Basin. The South Nicholson Group and its probable correlatives, the Pilpah Sandstone and Quamby Conglomerate, comprise a later phase of now largely eroded deposits within the Mount Isa Basin. The name South Nicholson Basin is now outmoded as this terminology only applied to the South Nicholson Group unlike the original broader definition in Brown et al. (1968). Cored slimhole stratigraphic and mineral wells drilled by Amoco, Esso, Elf Aquitaine and Carpentaria Exploration prior to 1986, penetrated much of the stratigraphy and intersected both minor oil and gas shows plus excellent potential source rocks. The raw data were reinterpreted and augmented with seismic stratigraphy and source rock data from resampled mineral and petroleum stratigraphic exploration wells for this study. Since 1986, Comalco Aluminium Limited, as operator of a joint venture with Monument Resources Australia Limited and Bridge Oil Limited, recorded approximately 1000 km of reflection seismic data within the basin and drilled one conventional stratigraphic petroleum well, Beamesbrook-1. This work was the first reflection seismic and first conventional petroleum test of the northern Mount Isa Basin. When incorporated into the newly developed foreland basin and maturity models, a grass roots petroleum exploration play was recognised and this led to the present thesis. The Mount Isa Basin was seen to contain excellent source rocks coupled with potential reservoirs and all of the other essential aspects of a conventional petroleum exploration play. This play, although high risk, was commensurate with the enormous and totally untested petroleum potential of the basin. The basin was assessed for hydrocarbons in 1992 with three conventional exploration wells, Desert Creek-1, Argyle Creek-1 and Egilabria-1. These wells also tested and confrrmed the proposed basin model. No commercially viable oil or gas was encountered although evidence of its former existence was found. In addition to the petroleum exploration, indeed as a consequence of it, the association of the extensive base metal and other mineralisation in the Mount Isa Basin with hydrocarbons could not be overlooked. A comprehensive analysis of the available data suggests a link between the migration and possible generation or destruction of hydrocarbons and metal bearing fluids. Consequently, base metal exploration based on hydrocarbon exploration concepts is probably. the most effective technique in such basins. The metal-hydrocarbon-sedimentary basin-plate tectonic association (analogous to Phanerozoic models) is a compelling outcome of this work on the Palaeo- to Mesoproterozoic Mount lsa Basin. Petroleum within the Bowthom Block was apparently destroyed by hot brines that produced many ore deposits elsewhere in the basin.

The largest volcanic eruptions on Earth

Large igneous provinces (LIPs) are sites of the most frequently recurring, largest volume basaltic and silicic eruptions in Earth history. These large-volume (N1000 km3 dense rock equivalent) and large-magnitude (NM8) eruptions produce areally extensive (104–105 km2) basaltic lava flow fields and silicic ignimbrites that are the main building blocks of LIPs. Available information on the largest eruptive units are primarily from the Columbia River and Deccan provinces for the dimensions of flood basalt eruptions, and the Paraná–Etendeka and Afro-Arabian provinces for the silicic ignimbrite eruptions. In addition, three large-volume (675– 2000 km3) silicic lava flows have also been mapped out in the Proterozoic Gawler Range province (Australia), an interpreted LIP remnant. Magma volumes of N1000 km3 have also been emplaced as high-level basaltic and rhyolitic sills in LIPs. The data sets indicate comparable eruption magnitudes between the basaltic and silicic eruptions, but due to considerable volumes residing as co-ignimbrite ash deposits, the current volume constraints for the silicic ignimbrite eruptions may be considerably underestimated. Magma composition thus appears to be no barrier to the volume of magma emitted during an individual eruption. Despite this general similarity in magnitude, flood basaltic and silicic eruptions are very different in terms of eruption style, duration, intensity, vent configuration, and emplacement style. Flood basaltic eruptions are dominantly effusive and Hawaiian–Strombolian in style, with magma discharge rates of ~106–108 kg s−1 and eruption durations estimated at years to tens of years that emplace dominantly compound pahoehoe lava flow fields. Effusive and fissural eruptions have also emplaced some large-volume silicic lavas, but discharge rates are unknown, and may be up to an order of magnitude greater than those of flood basalt lava eruptions for emplacement to be on realistic time scales (b10 years). Most silicic eruptions, however, are moderately to highly explosive, producing co-current pyroclastic fountains (rarely Plinian) with discharge rates of 109– 1011 kg s−1 that emplace welded to rheomorphic ignimbrites. At present, durations for the large-magnitude silicic eruptions are unconstrained; at discharge rates of 109 kg s−1, equivalent to the peak of the 1991 Mt Pinatubo eruption, the largest silicic eruptions would take many months to evacuate N5000 km3 of magma. The generally simple deposit structure is more suggestive of short-duration (hours to days) and high intensity (~1011 kg s−1) eruptions, perhaps with hiatuses in some cases. These extreme discharge rates would be facilitated by multiple point, fissure and/or ring fracture venting of magma. Eruption frequencies are much elevated for large-magnitude eruptions of both magma types during LIP-forming episodes. However, in basaltdominated provinces (continental and ocean basin flood basalt provinces, oceanic plateaus, volcanic rifted margins), large magnitude (NM8) basaltic eruptions have much shorter recurrence intervals of 103–104 years, whereas similar magnitude silicic eruptions may have recurrence intervals of up to 105 years. The Paraná– Etendeka province was the site of at least nine NM8 silicic eruptions over an ~1 Myr period at ~132 Ma; a similar eruption frequency, although with a fewer number of silicic eruptions is also observed for the Afro- Arabian Province. The huge volumes of basaltic and silicic magma erupted in quick succession during LIP events raises several unresolved issues in terms of locus of magma generation and storage (if any) in the crust prior to eruption, and paths and rates of ascent from magma reservoirs to the surface.

New insights into Crustal Contributions to Large-volume Rhyolite Generation in the Mid-Tertiary Sierra Madre Occidental Province, Mexico, revealed by U Pb Geochronology

Voluminous (≥3·9 × 105 km3), prolonged (∼18 Myr) explosive silicic volcanism makes the mid-Tertiary Sierra Madre Occidental province of Mexico one of the largest intact silicic volcanic provinces known. Previous models have proposed an assimilation–fractional crystallization origin for the rhyolites involving closed-system fractional crystallization from crustally contaminated andesitic parental magmas, with <20% crustal contributions. The lack of isotopic variation among the lower crustal xenoliths inferred to represent the crustal contaminants and coeval Sierra Madre Occidental rhyolite and basaltic andesite to andesite volcanic rocks has constrained interpretations for larger crustal contributions. Here, we use zircon age populations as probes to assess crustal involvement in Sierra Madre Occidental silicic magmatism. Laser ablation-inductively coupled plasma-mass spectrometry analyses of zircons from rhyolitic ignimbrites from the northeastern and southwestern sectors of the province yield U–Pb ages that show significant age discrepancies of 1–4 Myr compared with previously determined K/Ar and 40Ar/39Ar ages from the same ignimbrites; the age differences are greater than the errors attributable to analytical uncertainty. Zircon xenocrysts with new overgrowths in the Late Eocene to earliest Oligocene rhyolite ignimbrites from the northeastern sector provide direct evidence for some involvement of Proterozoic crustal materials, and, potentially more importantly, the derivation of zircon from Mesozoic and Eocene age, isotopically primitive, subduction-related igneous basement. The youngest rhyolitic ignimbrites from the southwestern sector show even stronger evidence for inheritance in the age spectra, but lack old inherited zircon (i.e. Eocene or older). Instead, these Early Miocene ignimbrites are dominated by antecrystic zircons, representing >33 to ∼100% of the dated population; most antecrysts range in age between ∼20 and 32 Ma. A sub-population of the antecrystic zircons is chemically distinct in terms of their high U (>1000 ppm to 1·3 wt %) and heavy REE contents; these are not present in the Oligocene ignimbrites in the northeastern sector of the Sierra Madre Occidental. The combination of antecryst zircon U–Pb ages and chemistry suggests that much of the zircon in the youngest rhyolites was derived by remelting of partially molten to solidified igneous rocks formed during preceding phases of Sierra Madre Occidental volcanism. Strong Zr undersaturation, and estimations for very rapid dissolution rates of entrained zircons, preclude coeval mafic magmas being parental to the rhyolite magmas by a process of lower crustal assimilation followed by closed-system crystal fractionation as interpreted in previous studies of the Sierra Madre Occidental rhyolites. Mafic magmas were more probably important in providing a long-lived heat and material flux into the crust, resulting in the remelting and recycling of older crust and newly formed igneous materials related to Sierra Madre Occidental magmatism.

Petrology and geochemistry of the Quaternary Caldera-forming, Phonolitic Granadilla Eruption, Tenerife (Canary Islands)

The Granadilla eruption at 600 ka was one of the largest phonolitic explosive eruptions from the Las Cañadas volcano on Tenerife, producing a classical plinian eruptive sequence of a widespread pumice fall deposit overlain by an ignimbrite. The eruption resulted in a major phase of caldera collapse that probably destroyed the shallow-level magma chamber system. Granadilla pumices contain a diverse phenocryst assemblage of alkali feldspar + biotite + sodian diopside to aegirine–augite + titanomagnetite + ilmenite + nosean/haüyne + titanite + apatite; alkali feldspar is the dominant phenocryst and biotite is the main ferromagnesian phase. Kaersutite and partially resorbed plagioclase (oligoclase to sodic andesine) are present in some eruptive units, particularly in pumice erupted during the early plinian phase, and in the Granadilla ignimbrite at the top of the sequence. Associated with the kaersutite and plagioclase are small clots of microlitic plagioclase and kaersutite interpreted as quenched blebs of tephriphonolitic magma within the phonolite pumice. The Granadilla Member has previously been recognized as an example of reverse-then-normal compositional zonation, where the zonation is primarily expressed in terms of substantial variations in trace element abundances with limited major element variation (cryptic zonation). Evidence for cryptic zonation is also provided by the chemistry of the phenocryst phases, and corresponding changes in intensive parameters (e.g. T, f O2, f  H2O). Geothermometry estimates indicate that the main body of phonolite magma had a temperature gradient from 860 °C to ∼790 °C, with hotter magma (≥900 °C) tapped at the onset and terminal phases of the eruption. The reverse-then-normal chemical and thermal zonation reflects the initial tapping of a partially hybridized magma (mixing of phonolite and tephriphonolite), followed by the more sequential tapping of a zoned and relatively large body of highly evolved phonolite at a new vent and during the main plinian phase. This suggests that the different magma types within the main holding chamber could have been laterally juxtaposed, as well as in a density-stratified arrangement. Correlations between the presence of mixed phenocryst populations (i.e. presence of plagioclase and kaersutite) and coarser pumice fall layers suggest that increased eruption vigour led to the tapping of hybridized and/or less evolved magma probably from greater depths in the chamber. New oxygen isotope data for glass and mineral separates preclude syn-eruptive interaction between the vesiculating magma and hydrothermal fluids as the cause of the Sr isotope disequilibrium identified previously for the deposit. Enrichment in radiogenic Sr in the pumice glass has more likely been due to low-temperature exchange with meteoric water that was enriched in 87Sr by sea spray, which may be a common process affecting porous and glassy pyroclastic deposits on oceanic islands.

Reconstructing the event stratigraphy from the complex structural - stratigraphic architecture of an Archaean volcanic – intrusive – sedimentary succession : the Boorara Domain, Eastern Goldfields Superterrane, Western Australia

Two main deformational phases are recognised in the Archaean Boorara Domain of the Kalgoorlie Terrane, Eastern Goldfields Superterrane, Yilgarn Craton, Western Australia, primarily involving southover- north thrust faulting that repeated and thickened the stratigraphy, followed by east northeast – west-southwest shortening that resulted in macroscale folding of the greenstone lithologies. The domain preserves mid-greenschist facies metamorphic grade, with an increase to lower amphibolite metamorphic grade towards the north of the region. As a result of the deformation and metamorphism, individual stratigraphic horizons are difficult to trace continuously throughout the entire domain. Volcanological and sedimentological textures and structures, primary lithological contacts, petrography and geochemistry have been used to correlate lithofacies between faultbounded structural blocks. The correlated stratigraphic sequence for the Boorara Domain comprises quartzo-feldspathic turbidite packages, overlain by high-Mg tholeiitic basalt (lower basalt), coherent and clastic dacite facies, intrusive and extrusive komatiite units, an overlying komatiitic basalt unit (upper basalt), and at the stratigraphic top of the sequence, volcaniclastic quartz-rich turbidites. Reconstruction of the stratigraphy and consideration of emplacement dynamics has allowed reconstruction of the emplacement history and setting of the preserved sequence. This involves a felsic, mafic and ultramafic magmatic system emplaced as high-level intrusions, with localised emergent volcanic centres, into a submarine basin in which active sedimentation was occurring.

Subliquidus glass-glass immiscibility along the albite-diopside join

A calorimetric study has shown that glasses along the albite-diopside join in the system albiteanorthite-diopside have positive enthalpies of mixing. Thermodynamic calculations based on these data describe a nearly symmetric, metastable, subliquidus irascibility gap along the join with a critical temperature at 910 K. The existence of the miscibility gap was tested experimentally by annealing an Ab50Di50 glass at 748 K and 823 K. Annealed glasses were examined by optical microscopy and by scanning and transmission electron microscopy. The glasses showed morphological and chemical features consistent with unmixing of two glass phases. The apparent mechanism of phase separation involves initial spinodal decomposition followed by coarsening to produce 0.1 μm–0.3 μm spherical glass phases.

Micas in experimentally shocked gneiss

Biotites and muscovites from a gneiss have been experimentally shocked between 18 and 70 GPa using powder-propellant guns at NASA Johnson Space Center and at the California Institute of Technology. This study shows that shock in biotite and muscovite can produce homogeneous and devolatilized glasses within microseconds. Shock-deformed micas display fracturing, kinking, and complex extinction patterns over the entire pressure range investigated. However, these deformation features are not a sensitive pressure indicator. Localized melting of micas begins at 33 GPa and goes to completion at 70 GPa. Melted biotite and muscovite are optically opaque, but show extensive microvesiculation and flow when observed with the SEM. Electron diffraction confirms that biotite and muscovite have transformed to a glass. The distribution of vesicles in shock-vitrified mica shows escape of volatiles within the short duration of the shock experiment. Experimentally shocked biotite and muscovite undergo congruent melting. Compositions of the glasses are similar to the unshocked micas except for volatiles (H2O loss and K loss). These unusual glasses derived from mica may be quenched by rapid cooling conditions during the shock experiment. Based on these results, the extremely low H2O content of tektites may be reconciled with a terrestrial origin by impact. Release of volatiles in shock-melted micas affects the melting behavior of coexisting dry silicates during the short duration of the shock experiment. Transportation and escape of volatiles released from shock-melted micas may provide plausible mechanisms for the origin of protoatmospheres on terrestrial planets, hydrothermal activity on phyllosilicate-rich meteorite parent bodies, and fluid entrapment in meteorites.

Integrated TEM, XRD and electron-microprobe investigation of mixed-layer chlorite smectite from the Point-Sal ophiolite, California

Five basalt samples from the Point Sal ophiolite, California, were examined using HRTEM and AEM in order to compare observations with interpretations of XRD patterns and microprobe analyses. XRD data from ethylene-glycol-saturated samples indicate the following percentages of chlorite in mixed-layer chlorite-smectite identified for each specimen: (i) L2036 almost-equal-to 50%, (ii) L2035 almost-equal-to 70 and 20%, (iii) 1A-13 almost-equal-to 70%, (iv) 1B-42 almost-equal-to 70%, and (v) 1B-55 = 100%. Detailed electron microprobe analyses show that 'chlorite' analyses with high Si, K, Na and Ca contents are the result of interlayering with smectite-like layers. The Fe/(Fe + Mg) ratios of mixed-layer phyllosilicates from Point Sal samples are influenced by the bulk rock composition, not by the percentage of chlorite nor the structure of the phyllosilicate. Measurements of lattice-fringe images indicate that both smectite and chlorite layers are present in the Point Sal samples in abundances similar to those predicted with XRD techniques and that regular alternation of chlorite and smectite occurs at the unit-cell scale. Both 10- and 14-angstrom layers were recorded with HRTEM and interpreted to be smectite and chlorite, respectively. Regular alternation of chlorite and smectite (24-angstrom periodicity) occurs in upper lava samples L2036 and 1A-13, and lower lava sample 1B-42 for as many as seven alternations per crystallite with local layer mistakes. Sample L2035 shows disordered alternation of chlorite and smectite, with juxtaposition of smectite-like layers, suggesting that randomly interlayered chlorite (< 0.5)-smectite exists. Images of lower lava sample 1B-55 show predominantly 14-angstrom layers. Units of 24 angstrom tend to cluster in what may otherwise appear to be disordered mixtures, suggesting the existence of a corrensite end-member having thermodynamic significance.