No description available.

A continuum of rhyolite compositions has been observed throughout the Taupo Volcanic Zone (TVZ) over the past 550 kyr. reflecting changes in the ƒH2O, ƒO₂, and P-T conditions in a lower crustal 'hot-zone' (10-30 km) where these evolved melts are generated by crystal fractionation of successively intruded basaltic magmas. The rhyolite compositional continuum is bound by two distinct end-member types: R1 is characterized by hydrous minerals (hornblende ± biotite), low FeO*/MgO (calc-alkaline series), low MREE, Y, and Zr, and high Sr; and R2 is characterized by anhydrous minerals (orthopyroxene ± clinopyroxene), high FeO*/MgO (tholeiitic series), high MREE, Y, and Zr, and low Sr. Slab-derived aqueous fluid components (Ba, Cl) correlate well with oxygen fugacity, and other well defined characteristics of silicic magmas in the Taupo Volcanic Zone (TVZ) between a cold-wet-oxidizing magma type (R1: amphibole ± biotite; high Sr, low Zr and FeO*/MgO, depleted MREE) and a hot-dry-reducing magma type (R2: orthopyroxene ± clinopyroxene; low Sr, high Zr, and FeO*/MgO, less depleted MREE). Oxygen fugacity was obtained from analysis of Fe-Ti oxides and ranges between -0.039 to +2.054 log units (ΔQFM; where QFM = quartz + fayalite + magnetite buffer) and is positively correlated with the bulk-rock Ba/La ratio, indicating that slab-derived fluid is the oxidizing agent in the rhyolites. Chlorine contents in hornblende also correlate with the bulk-rock Ba/La ratio. Hence, high fluid-flux typically correlates with the R1 and low fluid-flux with R2 rhyolite magma types. A geochemical evolution and distribution can be tracked in time and space throughout the central region of the TVZ from 550 ka to present and has revealed two distinct magmatic cycles that vary in length. The first cycle included widespread R1 type magmatism across the central TVZ beginning ca. 550 ka and was directly associated with previously unreported dome-building and ignimbrite-forming volcanism, and led to a voluminous (>3000 km³) ignimbrite 'flare-up' between ca. 340 and 240 ka. These magmas also display the highest K₂O and Pb isotopic compositions compared to those erupted more recently, and is consistent with a peak in slab-derived sediment input. The second cycle began roughly 180 ka, erupting ca. 800 km³ of magma, and continues to the present. The duration, rate, and composition of melt production within these cycles appears to be governed by the flux of fluid/sediment released from the subducting slab, while the distribution of melts may be governed more by extension along the central rift axis. The Matahina Ignimbrite (~160 km³ rhyolite magma; 330 ka) was deposited during a caldera-forming eruption from the Okataina Volcanic Centre, TVZ. The outflow sheet is distributed primarily from the northeast to southeast and consists of a basal plinian fall member and three ash-flow members. Pumice clasts are separated into three groups defined by differences in bulk geochemistry and mineral contents: high CaO, MgO, Fe₂O₃T, TiO₂, and low Al₂O₃, +hornblende (A2), low CaO, MgO, Fe2O3T, TiO2, ±hornblende (A1), and a subset to A1, which has high-K, +biotite (B). Two types of crystal-rich mafic clasts were also deposited during the final stages of the eruption. The distinct A and B rhyolite magma types are petrogenetically related to corresponding type A and B andesitic magma by up to 50% crystal fractionation under varying ƒO₂-ƒH₂O conditions. Further variations in the low- to high-silica rhyolites can be accounted for by up to 25% crystal fractionation, again under distinct ƒO₂-ƒH₂O conditions. Reconstruction of the P-T-ƒO₂-ƒ’H₂O conditions of the andesite to rhyolite magmas are consistent with the existence of a compositional and thermal gradient prior to the eruption. Magma mingling/mixing between the basalt to andesite and main compositionally zoned rhyolitic magma occurred during caldera-collapse, modifying the least-evolved rhyolite at the bottom of the reservoir and effectively destroying the pre-eruptive gradients. A detailed examination of the diverse range of calcic-amphibole compositions from the ca. 330 ka Matahina eruption (ca. 160 km³ rhyolitic magma) of the Okataina Volcanic Complex, Taupo Volcanic Zone, including crystal-rich basalt to dacite pumice from post-collapse deposits, reveals several pre- and syn-eruption magmatic processes. (1) Amphibole phenocrysts in the basaltic-andesite and andesite crystallized at the highest pressures and temperatures (P: up to 0.6±0.06 GPa and T: up to 950°C), equivalent to mid-crustal depths (13-22 km). Inter- and intra-crystalline compositions range from Ti-magnesiohornblende → Ti-tschermakite → tschermakite → magnesiohornblende and some display gradual decreases in T from core to rim, both consistent with magma differentiation by cooling at depth. (2) The largest amphibole crystals from the basaltic-andesite to andesite display several core to rim increases in T (up to 70°C), indicating new hotter magma periodically fluxed the crystal mush. (3) The dominant population of amphibole (magnesiohornblende) from the rhyolite is small and bladed and crystallized at low P-T conditions (P: 0.3 GPa, T: 765°C), equivalent to the eruptive P-T conditions. Amphibole (tschermakite-magnesiohornblende) from the dacitic and low-silica rhyolitic pumice form two distinct populations, which nucleated at two different T (High: 820°C and Low: 750°C). These compositional variations, governed primarily by differences in T conditions during crystal growth, record the mixing of two distinct amphibole populations that approached a thermal equilibrium at the eruptive T. Therefore, the diversity in amphibole compositions can be reconciled as an exchange of crystals+liquid between the basaltic-andesite to dacite from the mid-crust and rhyolite from the upper-crust, which quenched against one another, modifying the dacite to low-silica rhyolite compositions as the eruption progressed.

Sixteen percent of New Zealand’s power comes from geothermal sources which are primarily located within the Taupo Volcanic Zone (TVZ). The TVZ hosts twenty three geothermal fields, seven of which are currently utilised for power generation. Ngatamariki Geothermal Field is the latest geothermal power generation site in New Zealand, located approximately 15 km north of Taupo. This was the location of interest in this project, with testing performed on a range of materials to ascertain the physical properties and microstructure of reservoir rocks. The effect of burial diagenesis on the physical properties was also investigated. Samples of reservoir rocks were taken from the Tahorakuri Formation and Ngatamariki Intrusive Complex from a range of wells and depths (1354-3284 mbgl). The samples were divided into four broad lithologies: volcaniclastic lithic tuff, primary tuff, welded ignimbrite and tonalite. From the supplied samples twenty one small cylinders (~40-50mm x 20-25mm) were prepared and subjected to the following analyses: dual weight porosity, triple weight porosity, dry density, ultrasonic velocity (saturated and dry) and permeability (over a range of confining pressures). Thin sections impregnated with an epoxy fluorescent dye were created from offcuts of each cylinder and were analysed using polarised light microscopy and quantitative fluorescent light microstructural microscopy. The variety of physical testing allowed characterisation of the physical properties of reservoir rocks within the Ngatamariki Geothermal Field. Special attention was given to the petrological and mineralogical fabrics and their relation to porosity and matrix permeability. It was found that the pore structures (microfractures or vesicles) had a large influence on the physical properties. Microfractured samples were associated with low porosity and permeability, while the vesicular samples were associated with high porosity and permeability. The microfractured samples showed progressively lower permeability with increased confining pressure whereas samples with a vesicular microstructure showed little response to increased confining pressure. An overall trend of decreasing porosity and permeability with increasing density and sonic velocity was observed with depth, however large fluctuations with depth indicate this trend may be uncertain. The large variations correlate with changes in lithology suggest that the lithology is the primary control of the physical properties with burial diagenesis being a subsidiary factor. This project has established a relationship between the microstructure and permeability, with vesicular samples showing high permeability and little response increased confining pressure. The effects of burial diagenesis on the physical properties are subsidiary to the observed variations in lithology. The implications of these results suggest deep drilling in the Tahorakuri Formation may reveal unexploited porosity and permeability at depth.

The central segment of the Taupo Volcanic Zone (TVZ) is one of the world’s most productive areas of silicic volcanism and geothermal activity. Rhyolites largely predominate the eruptive output in the central TVZ, with only minor basalts, andesites and dacites. The rhyolites show diversity in composition, and form a compositional continuum between two end-member types (R1 and R2), as suggested in previous studies. In this thesis I present results from a quartz- (and rare plagioclase-) hosted melt inclusions study, focussing on the volatile concentration (i.e. H2O, Cl, F, CO2) and their relative distribution between R1 and R2 rhyolites. The main objective is to add further constraints on the magmatic systems with regard to their contribution to the hydrothermal systems in the central TVZ. A comparative study between R1 and R2 melt inclusions show distinct volatile, fluid-mobile, and highly incompatible element compositions. Differences in the bulk volatile concentration of the parental magmas (i.e. basalts intruding the lower crust) are suggested to be at the origin of these volatile disparities. Further analysis on the volatile exsolution of R1 and R2 melts lead to the observation that the two rhyolite types exsolve a volatile phase at different stages in their magmatic history. From Cl and H2O concentrations, it is suggested that R1 magmas exsolve a vapour phase first, whereas R2 rhyolites more likely exsolve a hydrosaline fluid phase. These results have considerable implications for the magmatic contribution into the hydrothermal systems in the central TVZ, as differences in the composition of the resulting volatile phase may be expected. The hydrothermal systems in the central TVZ are subdivided into two groups based on their gas and fluid chemistry; and the current model suggests that there are two distinct contributions: a typical ‘arc’ system, with geochemical affinity with andesitic fluids, located along the eastern margin of the TVZ, and a typical ‘rift’ system, with geochemical affinity with rhyolitic/basaltic fluids, located along the central and/or western region of the TVZ. The addition of the new data on the rhyolitic melt inclusions, leads to a re-evaluation of the magmatic contribution into the hydrothermal systems, with a particular focus on B and Cl. The results indicate a more diverse variety of contributions to the meteoric water in the hydrothermal systems, and also show that the east-west distribution of ‘arc’ and ‘rift’ fluids is not a viable model for the central TVZ. This work emphasises that melt inclusion data and their volatile degassing history cannot be underestimated when characterising and quantifying the magmatic component in hydrothermal fluids. The melt inclusion data also provide further insight into the pre-eruptive magmatic plumbing systems and are particularly important from a hazard perspective. Included in the thesis is a detailed petrological analysis of rhyolite melt inclusions across the central TVZ and an interpretation that large silicic magma systems (in the TVZ) are typically comprised of multiple batches of magma emplaced at some of the shallowest depths on Earth. Tectonic activity is suggested to play an important role in triggering large caldera-forming eruptions as the evacuation of one magma batch could cause a regional-scale readjustment that is sufficient enough to trigger and allow simultaneous eruption of an adjacent melt batch.

A continuum of rhyolite compositions has been observed throughout the Taupo Volcanic Zone (TVZ) over the past 550 kyr. reflecting changes in the ƒH2O, ƒO₂, and P-T conditions in a lower crustal 'hot-zone' (10-30 km) where these evolved melts are generated by crystal fractionation of successively intruded basaltic magmas. The rhyolite compositional continuum is bound by two distinct end-member types: R1 is characterized by hydrous minerals (hornblende ± biotite), low FeO*/MgO (calc-alkaline series), low MREE, Y, and Zr, and high Sr; and R2 is characterized by anhydrous minerals (orthopyroxene ± clinopyroxene), high FeO*/MgO (tholeiitic series), high MREE, Y, and Zr, and low Sr. Slab-derived aqueous fluid components (Ba, Cl) correlate well with oxygen fugacity, and other well defined characteristics of silicic magmas in the Taupo Volcanic Zone (TVZ) between a cold-wet-oxidizing magma type (R1: amphibole ± biotite; high Sr, low Zr and FeO*/MgO, depleted MREE) and a hot-dry-reducing magma type (R2: orthopyroxene ± clinopyroxene; low Sr, high Zr, and FeO*/MgO, less depleted MREE). Oxygen fugacity was obtained from analysis of Fe-Ti oxides and ranges between -0.039 to +2.054 log units (ΔQFM; where QFM = quartz + fayalite + magnetite buffer) and is positively correlated with the bulk-rock Ba/La ratio, indicating that slab-derived fluid is the oxidizing agent in the rhyolites. Chlorine contents in hornblende also correlate with the bulk-rock Ba/La ratio. Hence, high fluid-flux typically correlates with the R1 and low fluid-flux with R2 rhyolite magma types. A geochemical evolution and distribution can be tracked in time and space throughout the central region of the TVZ from 550 ka to present and has revealed two distinct magmatic cycles that vary in length. The first cycle included widespread R1 type magmatism across the central TVZ beginning ca. 550 ka and was directly associated with previously unreported dome-building and ignimbrite-forming volcanism, and led to a voluminous (>3000 km³) ignimbrite 'flare-up' between ca. 340 and 240 ka. These magmas also display the highest K₂O and Pb isotopic compositions compared to those erupted more recently, and is consistent with a peak in slab-derived sediment input. The second cycle began roughly 180 ka, erupting ca. 800 km³ of magma, and continues to the present. The duration, rate, and composition of melt production within these cycles appears to be governed by the flux of fluid/sediment released from the subducting slab, while the distribution of melts may be governed more by extension along the central rift axis. The Matahina Ignimbrite (~160 km³ rhyolite magma; 330 ka) was deposited during a caldera-forming eruption from the Okataina Volcanic Centre, TVZ. The outflow sheet is distributed primarily from the northeast to southeast and consists of a basal plinian fall member and three ash-flow members. Pumice clasts are separated into three groups defined by differences in bulk geochemistry and mineral contents: high CaO, MgO, Fe₂O₃T, TiO₂, and low Al₂O₃, +hornblende (A2), low CaO, MgO, Fe2O3T, TiO2, ±hornblende (A1), and a subset to A1, which has high-K, +biotite (B). Two types of crystal-rich mafic clasts were also deposited during the final stages of the eruption. The distinct A and B rhyolite magma types are petrogenetically related to corresponding type A and B andesitic magma by up to 50% crystal fractionation under varying ƒO₂-ƒH₂O conditions. Further variations in the low- to high-silica rhyolites can be accounted for by up to 25% crystal fractionation, again under distinct ƒO₂-ƒH₂O conditions. Reconstruction of the P-T-ƒO₂-ƒ’H₂O conditions of the andesite to rhyolite magmas are consistent with the existence of a compositional and thermal gradient prior to the eruption. Magma mingling/mixing between the basalt to andesite and main compositionally zoned rhyolitic magma occurred during caldera-collapse, modifying the least-evolved rhyolite at the bottom of the reservoir and effectively destroying the pre-eruptive gradients. A detailed examination of the diverse range of calcic-amphibole compositions from the ca. 330 ka Matahina eruption (ca. 160 km³ rhyolitic magma) of the Okataina Volcanic Complex, Taupo Volcanic Zone, including crystal-rich basalt to dacite pumice from post-collapse deposits, reveals several pre- and syn-eruption magmatic processes. (1) Amphibole phenocrysts in the basaltic-andesite and andesite crystallized at the highest pressures and temperatures (P: up to 0.6±0.06 GPa and T: up to 950°C), equivalent to mid-crustal depths (13-22 km). Inter- and intra-crystalline compositions range from Ti-magnesiohornblende → Ti-tschermakite → tschermakite → magnesiohornblende and some display gradual decreases in T from core to rim, both consistent with magma differentiation by cooling at depth. (2) The largest amphibole crystals from the basaltic-andesite to andesite display several core to rim increases in T (up to 70°C), indicating new hotter magma periodically fluxed the crystal mush. (3) The dominant population of amphibole (magnesiohornblende) from the rhyolite is small and bladed and crystallized at low P-T conditions (P: 0.3 GPa, T: 765°C), equivalent to the eruptive P-T conditions. Amphibole (tschermakite-magnesiohornblende) from the dacitic and low-silica rhyolitic pumice form two distinct populations, which nucleated at two different T (High: 820°C and Low: 750°C). These compositional variations, governed primarily by differences in T conditions during crystal growth, record the mixing of two distinct amphibole populations that approached a thermal equilibrium at the eruptive T. Therefore, the diversity in amphibole compositions can be reconciled as an exchange of crystals+liquid between the basaltic-andesite to dacite from the mid-crust and rhyolite from the upper-crust, which quenched against one another, modifying the dacite to low-silica rhyolite compositions as the eruption progressed.

The Ngatamariki Geothermal Field is located 20 km north of Taupo in the Taupo Volcanic Zone and has a boundary of 12 km² as delineated by magneto-telluric surveys (Urzua 2008). Rhyolitic deposits, derived from the Maroa Volcanic Centre, dominate the geology of the area with the 186 AD (Wilson et al. 2009) Taupo pumice mantling stream valleys in the area. The majority of thermal features at Ngatamariki are located along the Orakonui Stream on the western boundary of the field; the stream area is dominated by a 50x30 m geothermal pool filling a hydrothermal eruption crater. This crater was formed during a hydrothermal eruption in 1948, with a subsequent eruption in April 2005. Orakei Korako is located 7 km north of Ngatamariki and has one of the largest collections of thermal features in New Zealand. The geology at Orakei Korako is similar to Ngatamariki, but the area is dominated by a series of south-west trending normal faults which create sinter terraces on the eastern bank of Lake Ohakuri. Water samples from springs and wells at Ngatamariki and Orakei Korako were taken to assess the nature of both fields. Spring waters at Ngatamariki have chloride contents of 56 to 647 mg/l with deep waters from wells ranging from 1183 to 1574 mg/l. This variation is caused by mixing of deep waters with a steam heated groundwater, above clay caps within the reservoir. Stable isotopic results (δ¹⁸O and δD) suggest that reservoir waters are meteoric waters mixed with magmatic (andesitic) water at Ngatamariki. Reservoir water chemistry at Orakei Korako exhibits low chloride contents, which is anomalous in the Taupo Volcanic Zone. Chloride content in well and spring waters is similar ranging from 546 to 147 mg/l, due to mixing of reservoir fluids with a ‘hot water’ diluent at depth. Isotopic compositions of spring waters suggest that they are meteoric waters which mix with magmatic (rhyolitic) water, more enriched in δ¹⁸O and δD than ‘andesitic’ water. Relationships between major ion concentrations and known subsurface geology suggest there is no hydraulic connection between the two fields.

The Tahorakuri Formation was introduced as a stratigraphic term to simplify the sometimes complex and inconsistent naming conventions in subsurface deposits within the geothermal fields of the central Taupo Volcanic Zone (TVZ). It consists of all volcaniclastic and sedimentary deposits between the ~350 ka Whakamaru-group ignimbrites and the greywacke basement that cannot be correlated with known ignimbrites. As such, it represents a long period in which relatively little is known about the volcano-tectonic history of the TVZ. The thesis focuses on the Tahorakuri Formation at Ngatamariki and Rotokawa geothermal fields and the implications for the volcano-tectonic evolution of the TVZ. Drill cuttings from wells NM5 and NM6 are re-examined, and new U-Pb zircon dates from the Tahorakuri Formation are presented and implications discussed. Potassium feldspars identified in the drill cuttings from NM5 were examined by Raman spectroscopy and electron microprobe (EMP) analysis. Although petrographically many of the feldspars appear similar to sanidine, a primary volcanic mineral phase, this showed them to be adularia which formed during hydrothermal alteration. Raman spectroscopy was found to be ideal for analysing a large number of grains quickly, with the spectral peak at ~140 cm⁻¹ being particularly useful for identifying adularia as it is absent in sanidine. EMP analysis was found to be somewhat slower, but definitively identified the feldspars as adularia, with typical potassium-rich compositions of Or₉₄-Or₉₉. U-Pb dating shows that the Tahorakuri Formation formed over a very long time, with pyroclastic deposits ranging from 1.89 - 0.70 Ma. This was followed by a period with little or no explosive volcanism until ~0.35 Ma during which sediments were deposited at Ngatamariki. The periods at ~1.9 Ma and ~0.9 Ma were particularly active phases of pyroclastic deposition, with the second phase likely correlating with the Akatarewa ignimbrite. The oldest deposits overlie a large andesitic composite cone volcano. Significant subsidence of the andesite must have preceded emplacement of the silicic deposits, indicating that rifting within the central TVZ may have started earlier than previously thought. While the origin of the deposits is uncertain, the distribution of the oldest deposits outcropping at the surface, as well as the likely early initiation of rifting, would suggest a source within the TVZ is likely.

Ngauruhoe is an active basaltic andesite to andesite composite cone volcano at the southern end of the Tongariro volcanic complex, and most recently erupted in 1954-55 and 1974-75. These eruptions constructed the inner crater of Ngauruhoe, largely composed of 1954-55 deposits, which are the basis of this study. The inner crater stratigraphy, exposed on the southern wall, is divided into seven lithostratigraphic units (A to G), while the northern stratigraphy is obscured by the inward collapse of the crater rim. The units are, from oldest to youngest: Unit A, (17.5 m thick), a densely agglutinated spatter deposit with sharp clast outlines; Unit B, (11.2 m) a thick scoria lapilli deposit with local agglutination and scattered spatter bombs up to 1 m in length; Unit C, (6.4 m thick) a clastogenic lava deposit with lateral variations in agglutination; and Unit D, (10 m thick) a scoria lapilli with varying local agglutination. The overlying Unit E (15 cm thick) is a fine ash fallout bed that represents the final vulcanian phase of the 1954-55 eruption. Unit F is a series of six lapilli and ash beds that represent the early vulcanian episode of the 1974-75 eruption. The uppermost Unit G (averaging 10 m thick) is a densely agglutinated spatter deposit that represents the later strombolian phase of the 1974-75 eruption. Units A-D juvenile clasts are porphyritic, with phenocrysts of plagioclase, orthopyroxene, clinopyroxene, minor olivine, within a microlitic glassy groundmass. Quartzose and greywacke xenoliths are common in most units, and are derived from the underlying basement. The 1954-55 and 1974-75 eruptions are a product of a short-lived, continental arc medium-K calc-alkaline magma. The magma originated from the mantle, then filtered through the crust, undergoing assimilation and fractionation, and evolving to basaltic andesite and andesite compositions. The magma body stagnated in shallow reservoirs where it underwent further crustal assimilation and fractionation of plagioclase and olivine, and homogenisation through magma mixing. Prior to the 1954-55 eruption a more primitive magma body was incorporated into the melt. The melt homogenised and fed both the 1954-55 and 1974-75 eruptions, with a residence time of at least 20 years. The 1954-55 eruption produced alternating basaltic andesite and andesite strombolian activity and more intense fire fountaining, erupting scoria and spatter that built up the bulk of the inner crater. A period of relative quiescence allowed the formation of a cooled, solid cap rock that resulted in the accumulation of pressure due to volatile exsolution and bubble coalescence. The fracturing of the cap rock then resulted in a vulcanian eruption, depositing a thin layer of fine ash and ballistic blocks. The 1974-75 eruption commenced with the rupturing of the near-solid cap rock from the 1954-55 eruption in an explosive vulcanian blast, the result of decompressional volatile exsolution and bubble coalescence, and possible magma-water interaction. The eruption later changed to strombolian style, producing a clastogenic lava that partially flowed back into the crater.