Genesis of lavas of the Taupo Volcanic Zone, North Island, New Zealand

The Taupo Volcanic Zone forms part of the Taupo-Hikurangi subduction system, and comprises five volcanic centres: Tongariro, Taupo, Maroa, Okataina and Rotorua. Tongariro Volcanic Centre is formed almost entirely of andesite while the other four centres contain predominantly rhyolitic volcanics and later fissure eruptions of high-Al basalt. Estimated total volume of each lava type are as follows: 2 km³ of high-Al basalt (< 0.1%); 260 km³ of andesite (< 2.5%); 5 km³ of dacite (< 0.1%); > 10,000 km³ of rhyolite and ignimbrite (> 97.4%).The location of the andesites and vent alignments suggest a source from a subduction zone underlying the area. However, the lavas differ chemically from island-arc andesites such as those of Tonga; in particular by having higher contents of the alkali elements, light REE and Sr and Pb isotopes. This suggests some crustal contamination, and it is considered that this may occur beneath the wide accretionary prism of the subduction system. Amphibolite of the subduction zone will break down between 80 and 100 km and a partial melt will rise. A multi-stage process of magma genesis is then likely to occur. High-Al basalts are thought to be derived from partial melting of a garnet-free peridotite near the top of the mantle wedge overlying the subduction zone, locations of the vents controlled largely by faults within the crust. Rhyolites and ignimbrites were probably derived from partial melting of Mesozoic greywacke and argillite under the Taupo Volcanic Zone. Initial partial melting may have been due to hydration of the base of the crust; the “water” having come from dehydration of the downgoing slab. The partial melts would rise to form granodiorite plutons and final release of the magma to form rhyolites and ignimbrites was allowed because of extension within the Taupo graben.Dacites of the Bay of Plenty probably resulted from mixing of andesitic magma with small amounts of rhyolitic magma, but those on the eastern side of the Rotorua-Taupo area were more likely formed by a higher degree of partial melting of the Mesozoic greywacke-argillite basement. This may be due to intrusion of andesite magma on this side of the Taupo volcanic zone.     

Origin of dacites of Taupo Volcanic Zone, New Zealand

Dacites form a relatively small proportion of lavas in the Taupo Volcanic Zone, New Zealand (5km³), and occur mainly on the eastern side. In this paper their origin is considered in terms of three processes: (a) partial melting of crustal rocks; (b) fractional crystallisation of basalt and andesite; and (c) sub-surface mixing of basic and acid magma. Modelling techniques are used to calculate the most acceptable degree of fractional crystallisation and magma mixing to fit major-element data, and these values are used to compare calculated and observed trace-element values. The success or failure of the model is determined by the closeness of the two sets of values. For partial-melt models, trace-element values alone are calculated by the batch-melting equation.Results indicate that White Island dacite can best be modelled by fractional crystallisation; Manawahe by fractional crystallisation plus limited crustal contamination; Maungaongaonga by partial melting of Western Basement greywacke, and Tauhara by partial melting of this greywacke together with minor mixing with a more basic magma. Results from Parekauau and Horohoro indicate that these lavas are unlikely to have formed by any of the processes examined.     

Buchitic metagreywacke xenoliths from Mount Ngauruhoe,Taupo Volcanic Zone,New Zealand

Buchitic sedimentary xenoliths, a few centimetres to several decimetres diameter, occur in Recent andesite from Mount Ngauruhoe, Tongariro Volcanic Center, Taupo Volcanic Zone, New Zealand. Bulk chemistry and Sr isotope compositions of the xenoliths indicate that they are greywacke and argillite derived from Mesozoic Torlesse terrane basement that partly underlies the Taupo Volcanic Zone. The xenoliths contain up to 80% glass with quartz, apatite and zircon remaining as unmelted phases. Glasses within the xenoliths are peraluminous (A/CNK = 1.0 − 1.4), have high normative corundum (2–7%), appreciable FeO (2–4 wt.%), MgO (0.2–1.5 wt.%), TiO₂ (0.17–0.84 wt.%), relatively high normative An (1.0–5.3%), and do not represent S-type granitic melts. In the argillite the glass has higher amounts of AI₂O₃, FeO, MgO, CaO and K₂O, and has less SiO₂ and Na₂O than glass in the greywacke. Silica-rich glass (up to 80 wt.% SiO₂) surrounds partially melted quartz. Variable glass chemistry reflects the heterogeneous (layered) nature of the xenoliths. Cordierite (Mg/(Mg + Fe + Mn) = 0.78-0.58), orthopyroxene (En_43–56), Mg-rich ilmenite, rutile, pleonaste, V-Cr-Ti spinel, and pyrrhotite occur in the glass of the xenoliths. The dominant cordierite, orthopyroxene, spinel assemblage can be accounted for by disequilibrium breakdown reactions under low oxidation conditions < QFM) involving phengite and chlorite which are abundant in Torlesse greywacke and argillite cropping out along the eastern side of the Taupo Volcanic Zone. Comparison with glass compositions and phase relations of disequilibrium melting experiments on Torlesse greywacke and argillite indicates a minimum temperature of 775°C and a maximum pressure of 1.5 kbar for fusion of the xenoliths that underwent a rapid rate of heating at a depth of less than 5 km and a cooling period constrained by the time of quenching when they were erupted.     

Asymmetric,multiple-block collapse at Rotorua Caldera,Taupo Volcanic Zone,New Zealand

Calderas worldwide have been classified according to their dominant collapse styles, although there is a good deal of speculation about the processes involved. Recent laboratory experiments have tried to constrain these processes by modelling magma withdrawal and observing the effects on overlying materials. However, many other factors also contribute to final caldera morphology. Rotorua Caldera formed during the eruption of the Mamaku Ignimbrite. Collapse structure and evolution of Rotorua Caldera is interpreted based its geophysical response, geology and geomorphology, and the stratigraphy of the Mamaku Ignimbrite. Rotorua Caldera is situated at the edge of the extensional Taupo Volcanic Zone, in which major faults strike NE-SW. A second, less dominant fault set strikes NW-SE. These two fault sets have a strong influence on the morphology of Rotorua Caldera. No one style of collapse can be applied to Rotorua Caldera; it was formed during a single eruption, but subsided as many blocks and shows features of trapdoor, piecemeal and downsag types of collapse. Here Rotorua Caldera is described, according to its composition, activity and geometry, as a rhyolitic, single event, asymmetric, multiple-block, single locus collapse structure. The Mamaku Ignimbrite is the only ignimbrite to have erupted from Rotorua Caldera. Extracaldera thickness of the Mamaku Ignimbrite is up to 145 m, whereas inside the caldera it may be greater than 1 km thick. The Mamaku Ignimbrite can be separated into a basal tephra sequence and main ignimbrite sequence. The main ignimbrite sequence contains no observable flow unit boundaries but can be split into lower, middle and upper parts (LMI, mMI, uMI respectively) based on crystal content, welding, jointing, devitrification and vapour phase alteration. Juvenile clasts within the ignimbrite comprise three consanguineous silicic pumice types and andesitic fragments. Only the most evolved pumice type occurs in the basal tephra sequence. All three pumice types occur together throughout the main ignimbrite sequence, whereas the andesitic fragments are only present in uMI. Lithic lag breccias in uMI indicate a late stage of caldera collapse. Concentration of lithic fragments increases towards the middle of the ignimbrite, and may also reflect increased subsidence rate during an earlier stage. Collapse of Rotorua Caldera is thought to have occurred throughout the eruption of the main ignimbrite sequence of the Mamaku Ignimbrite, allowing simultaneous eruption of all the different pumice types and causing the abrupt transition from deposition of the basal tephra sequence to the main ignimbrite sequence.     

Taupo Volcanic Zone calc-alkaline tephras on the peralkaline Mayor Island volcano, New Zealand: identification and uses as marker horizons

Mayor Island is a peralkaline rhyolitic caldera volcano characterised by numerous, sector-confined pyroclastic deposits, together with lavas forming at least five composite shields. Correlation of sequences between sectors is difficult because of the scarcity of island-wide marker beds. However, eight distal calc-alkaline fall tephras (ca. 7.3 ¹⁴C ka to 64 ka) from Okataina and Taupo volcanic centres in the nearby Taupo Volcanic Zone (TVZ) have been identified on the island. These “foreign” TVZ tephras provide marker planes to correlate activity in different sectors of Mayor Island volcano, and refine an eruptive chronology. At least seventeen pyroclastic eruptions and fourteen lava-producing events (including multiple, shield-forming events) have occurred in the past ca. 64 ka. Age controls provided by the calc-alkaline tephras confirm the extremely local dispersal characteristics of many of the Mayor Island eruptives and show that K/Ar ages as young as 25–33 ka on obsidians with 4.2–4.4% K₂O are reliable.     

Glass-bearing plutonic fragments from ignimbrites of the Okataina caldera complex, Taupo Volcanic Zone, New Zealand: remnants of a partially molten intrusion associated with preceding eruptions

Glass-bearing plutonic fragments occur as rare accessory lithics within the ca. 64 ka Rotoiti and Earthquake Flat ignimbrites that were erupted from Okataina caldera complex, Taupo Volcanic Zone, New Zealand. Granitoid lithic fragments are only found in the Rotoiti ignimbrite and fall into two groups. Group 1 granitoids have textures consistent with a period of slow cooling followed by rapid quenching, and were excavated by the Rotoiti eruption from a single incompletely solidified magma body. Although isotopic ratios for the Group 1 granitoids are similar to the host ignimbrite, they are not cognate, having different chemistry, mineralogy, mineral chemistry and crystallisation history. It is more likely that they represent fragments of a separate incompletely solidified magma chamber that was intercepted by the erupting Rotoiti ignimbrite magma. Low LILE and high HFSE abundances favour a comagmatic link with the ca. 0.28 Ma Matahina ignimbrite and it is suggested they are derived from an isolated cupola of the Matahina magma chamber that remained at depth (between 3.5 and 5 kbar pressure) after eruption of the Matahina ignimbrite. Migration toward the surface probably accompanied development of the Rotoiti magma system in the upper crust. Most geochemical variation in Group 1 granitoids is related to the abundance of biotite, the concentration of which is controlled by differential shear. REE abundance is controlled by light REE-enriched accessory minerals preferentially included within biotite. Although Eu_n remains constant in the Group 1 granitoids, Eu/Eu^* varies systematically with (La/Yb)_n and is controlled by variations in Sm and Gd rather than in Eu. Group 2 granitoid fragments have a wide range of composition, comparable to many Okataina rhyolites, including those found as lithic fragments in the Rotoiti ignimbrite. Rare microdiorite fragments occur in both Rotoiti and Earthquake Flat ignimbrites and typically contain vesicular interstitial glass indicating that they were incompletely solidified prior to eruption. Those from the Rotoiti ignimbrite are comparable to the (>64 ka) Matahi basaltic tephra and probably represent part of the same magmatic event which generated the Matahi tephra.     

Variations in the chemical and isotopic composition of fluids discharged from the Taupo Volcanic Zone, New Zealand

Contents of H₂O, CO₂ and Cl in well discharges from six explored geothermal systems of the Taupo Volcanic Zone, New Zealand, point to the existence of two distinct source fluids. The fluid present in discharges from systems along the eastern boundary is characterised by high CO₂ contents, 1.6 ± 0.5 , at mole ratios of 3.9 ± 1.5. High (0.06) and (12) weight ratios in these waters suggest that all four constituents are derived from associated andesitic rock. Geothermal discharges in the western parts of the TVZ, dominated by rhyolitic magmatism, are characterised by low CO₂ contents, 0.12 ± 0.05 , and low (0.14 ± 0.1) ratios. Again, relative Cl, B, Li and Cs contents agree with those of this potential source rock. High and ratios in the east are typical of fluids affected by the addition of volatiles released from subducted marine sediments. For the western systems, these ratios resemble more closely those expected for mantle-derived volatiles. The isotopic compositions of all deep waters point to the presence of variable amounts of a magmatic component, some 14 ± 5% in the eastern and 6 ± 2% in the western systems. The observed variations are explained in terms of interaction of volatiles released from the subducted sediments with material of the mantle wedge to form a volatile-charged, high-alumina basalt. Its convective rise, in a direction opposite to that of the down-going slab, leads to high enrichment in volatiles of the magmas generated beneath the eastern parts of the TVZ and increases their ability to intrude the continental crust. Further fractional crystallisation and assimilation leads to the formation of volatile-rich andesitic melts, partly extruded to form the volcanoes of the andesitic arc, partly intruded to act as source rocks for the high-gas geothermal systems. Batches of high-alumina basalt, depleted in subducted volatiles, travel farther west to pond beneath a zone of crustal extension. Following extensive fractionation, highly siliceous melts, carrying predominantly mantle-type volatiles, rise beneath the western part of the TVZ to supply both heat and volatiles to the geothermal systems there.     

Pyroclastic phases of a rhyolitic dome-building eruption: Puketarata tuff ring,Taupo Volcanic Zone,New Zealand

The 14 ka Puketarata eruption of Maroa caldera in Taupo Volcanic Zone was a dome-related event in which the bulk of the 0.25 km³ of eruption products were emplaced as phreatomagmatic fall and surge deposits. A rhyolitic dike encountered shallow groundwater during emplacement along a NE-trending normal fault, leading to shallow-seated explosions characterised by low to moderate water/magma ratios. The eruption products consist of two lava domes, a proximal tuff ring, three phreatic collapse craters, and a widespread fall deposit. The pyroclastic deposits contain dominantly dense juvenile clasts and few foreign lithics, and relate to very shallow-level disruption of the growing dome and its feeder dike with relatively little involvement of country rock. The distal fall deposit, representing 88% of the eruption products is, despite its uniform appearance and apparently subplinian dispersal, a composite feature equivalent to numerous discrete proximal phreatomagmatic lapilli fall layers, each deposited from a short-lived eruption column. The Puketarata products are subdivided into four units related to successive phases of:(A) shallow lava intrusion and initial dome growth; (B) rapid growth and destruction of dome lobes; (C) slower, sustained dome growth and restriction of explosive disruption to the dome margins; and (D) post-dome withdrawal of magma and crater-collapse. Phase D was phreatic, phases A and C had moderate water: magma ratios, and phase B a low water: magma ratio. Dome extrusion was most rapid during phase B, but so was destruction, and hence dome growth was largely accomplished during phase C. The Puketarata eruption illustrates how vent geometry and the presence of groundwater may control the style of silicic volcanism. Early activity was dominated by these external influences and sustained dome growth only followed after effective exclusion of external water from newly emplaced magma.     

Mamaku Ignimbrite: a caldera-forming ignimbrite erupted from a compositionally zoned magma chamber in Taupo Volcanic Zone, New Zealand

Mamaku Ignimbrite was deposited during the formation of Rotorua Caldera, Taupo Volcanic Zone, New Zealand, 220–230 ka. Its outflow sheet forms a fan north, northwest and southwest of Rotorua, capping the Mamaku–Kaimai Plateau. Mamaku Ignimbrite can be divided into a partly phreatomagmatic basal sequence, and a main sequence which comprises lower, middle, and upper ignimbrite. The internal stratigraphy indicates that it was emplaced progressively from a pyroclastic density current of varying energy that became less particulate away from source. Gradational contacts between lower, middle, and upper ignimbrite are consistent with it being deposited during one eruptive event from the same source. Variations in lithic clast content and coexistence of different pumice types through the ignimbrite sequence indicate that caldera collapse occurred throughout the eruption, but particularly when middle Mamaku Ignimbrite was deposited and in the final stages of deposition of upper Mamaku Ignimbrite. Maximum lithic data and the location of lithic lag breccias in upper Mamaku Ignimbrite confirm Rotorua Caldera as the source. At least 120 m of geothermally altered intra-caldera Mamaku Ignimbrite occurs inside Rotorua Caldera. Pumice clasts in the Mamaku Ignimbrite are dacite to high-silica rhyolite and can be chemically divided into three types: high–silica rhyolite (type 1), rhyolite (type 2), and dacite (type 3). All are petrogenetically related and types 1 and 2 may be derived by up to 20% crystal fractionation from the type 3 dacite. All three types probably resided in a single, gradationally zoned magma chamber. Andesitic juvenile fragments are found only in upper Mamaku Ignimbrite and inferred to represent a discrete magma that was injected into the silicic chamber and is considered to have accumulated as a sill at the base of the magma chamber. The contrast in density between the andesitic and silicic magmas did not allow eruption of the andesitic fragments during the deposition of lower and middle Mamaku Ignimbrite. The advanced stage of caldera collapse, late in the main eruptive phase, created withdrawal dynamics that allowed andesitic magma to reach the surface as fragments within upper Mamaku Ignimbrite.     

Volcanic hazards of North Island, New Zealand-overview

In October 1980, a National Civil Defence Planning Committee on Volcanic Hazards was formed in New Zealand, and solicited reports on the likely areas and types of future eruptions, the risk to public safety, and the need for special precautions. Reports for eight volcanic centres were received, and made available to the authors. This paper summarises and quantifies the type and frequency of hazard, the public risk, and the possibilities for mitigation at the 7 main volcanic centres: Northland, Auckland, White Island, Okataina, Taupo, Tongariro, and Egmont. On the basis of Recent tephrostratigraphy, eruption probabilities up to 20% per century (but commonly 5%), and tephra volumes up to 100 km³ are credible.     

Petrogenesis of Tauhara Dacite (Taupo Volcanic Zone,New Zealand) - Evidence for magma mixing between high-alumina andesite and rhyolite

Tauhara dacites have petrographic, geochemical and isotopic characteristics which indicate an origin by magma mixing between andesite and rhyolite. Phenocrysts typically exhibit strong zoning near their rims, are resorbed or display fusion textures. Assemblages are not in equilibrium with host lavas and compositions are bimodal: plagioclase An_23–43 and An_66–91; orthopyroxene En_44–51 and En_69–79. Chemical and isotopic trends pass through the bulk compositions of high-alumina andresite and rhyolite which crop out in the vicinity of the dacite domes. Least squares mixing models indicate 40–75% of a rhyolite endmember mixed with andesite can generate the full range of dacite compositions. Subtle geochemical differences between domes suggest that magma mixing may have proceeded as three or more general episodes, each punctuated by several events. These episodes may have catalyzed some of the larger pyroclastic flow eruptions of Taupo Volcanic Zone in the past 50,000 years.     

Seismicity in the Rotorua and Kawerau geothermal systems,Taupo Volcanic Zone,New Zealand,based on improved velocity models and cross-correlation measurements

We analyze earthquakes occurring in and around the Rotorua and Kawerau geothermal systems, Taupo Volcanic Zone, New Zealand. The two data sets contain 504 and 1875 shallow (≤ 20 km deep) earthquakes, respectively, and span the 21 year period between 1984 and 2004. The arrival time data for these earthquakes are first used to calculate 1-D P- and S-wave seismic velocity models and accompanying station correction terms for both areas. In order to address the non-uniqueness of the joint hypocenter-velocity model estimation problem, we analyze suites of 1000 velocity models computed from random initial models. The final velocity models are well constrained, particularly at depths between 4 and 15 km, and consistent with the results obtained in previous seismic refraction studies of the central Taupo Volcanic Zone. Using a combination of cross-correlation-derived and catalog-based arrival times, we relocate subsets of the Rotorua and Kawerau data sets. In Rotorua, the relocated earthquakes cluster near the geothermally active parts of Rotorua City and beneath the Mount Ngongotaha rhyolite dome. Earthquake clusters and alignments reveal seismogenic structures in the mid-crust whose positions and geometries are consistent with previously published fault mechanisms and known near-surface faults. In Kawerau, the earthquakes within the geothermal field align along northeast-trending lineations, consistent with the predominant alignment of surface-mapped faults in the area.     

The Reporoa Caldera,Taupo Volcanic Zone: source of the Kaingaroa Ignimbrites

The Reporoa Caldera occupies the northern end of the Reporoa Depression, previously described as a tectonic fault-angle depression. Earlier confirmation of the topographic basin as a caldera had been hindered by the lack of an associated young pyroclastic flow deposit of large enough volume to have caused caldera collapse. New exposures on the eastern margin of the Reporoa basin reveal thick lithic lag breccias (>30 m) interbedded within the 0.24 Ma Kaingaroa Ignimbrites. These ignimbrites were previously attributed to the adjacent Okataina Volcanic Centre. Lag breccia thicknesses and maximum clast sizes decrease rapidly outward from the caldera rim, and discrete breccias are absent from ignimbrite sections more than 3 km from the rim. The lithic lag breccias, together with structural and geophysical evidence, confirm Reporoa Caldera as the source of the c. 100 km³ Kaingaroa Ignimbrites, adding another major rhyolitic volcanic centre to the seven previously recognized in the Taupo Volcanic Zone. Other, older, calderas may also be present in the Reporoa Depression.     

Crustal heat transfer in the Taupo Volcanic Zone (New Zealand): comparison with other volcanic arcs and explanatory heat source models

The Taupo Volcanic Zone (TVZ) is a 200-km-long volcanic arc segment which developed ≤2 Ma ago within the continental crust of the North Island of New Zealand and lies at the southern end of the much larger Tonga-Kermadec arc system. The total crustal heat transfer of the TVZ is at present c. 2600 MW/100 km, most of the heat being transferred by convective geothermal systems. The rate of transfer is anomalously high in comparison to that of other active arcs, and arguably the highest world wide for such a setting. Heat transfer of other active arcs appear to vary almost linearly with subduction speed (about 150 MW/100 km for 10 mm/yr). The mass rate of common type arc extrusions (basalts, andesites, dacites) also increases almost linearly with subduction speed. This allows separation of the TVZ heat transfer into a “normal” component, associated with extrusions and intrusions of andesites and dacites (about 600 MW/100 km), and an “anomalous” component of about 2000 MW/100 km, related to extrusions and intrusions of rhyolitic melts whose generation is not directly controlled by subduction processes.Rhyolitic melts in the TVZ are partial melts of dominantly crustal origin. Comparison with other arcs indicates that the long-term extrusion rate of TVZ rhyolites (about 400 kg/s per 100 km) is also the highest world wide for this setting. The occurrence of voluminous Quaternary rhyolitic pyroclastics is a rare phenomenon and appears to be associated with a few arc segments (TVZ, Sumatra, Kyushu) that undergo significant crustal deformation.Various models have been proposed to explain the phenomenon of the anomalously high heat transfer within the TVZ. Models which require only heat transfer from plumes and subcrustal melts, either ponded at the crust/mantle boundary or intruding a spreading crust, are not suitable because the associated heat transfer at the contact is too low by a factor 2 to explain the required transfer rate of about 0.8 W/m² representing the “anomalous” crustal heat component of the TVZ. Heat generation by focussed plastic deformation within the ductile lithosphere is an alternative mechanism to explain “endogenous crustal heating” which yields heating rates that are also too low by a factor of two, although important parameters (average yield strength of lithosphere and opening rate of the TVZ) are not well known. A further search for a suitable combination of heat source models is required.     

Lacustrine shore platforms at Lake Waikaremoana,North Island,New Zealand

This paper examines the morphology and processes governing the development of shore platforms at Lake Waikaremoana, North Island, New Zealand. Shore platforms at Lake Waikaremoana are recent features, and were formed when a new sequence of shoreline development was initiated, due to lowering of the lake by 5 m in 1946 for hydroelectric power development. Three predominant platform morphologies were identified around the lake. These include gently sloping platforms (c.1·5 to 3·9°), ramp platforms (c.6·8 to 9·2°), and concave ramp platforms (c.7·9 to 12°). Platform widths ranged from 11 to 31 m, with the gently sloping platforms characterized by the widest morphologies. Erosion rates were estimated using perched sandstone boulders and were found to range from 3·4 to 12·5 mm a^?1, with a mean erosion rate of 5·9 mm a^?1. Higher rates of erosion were identified at lower platform elevations, due to a greater frequency of wetting and drying cycles coincident with storm waves, while lower erosion rates were identified at higher elevations. Field evidence suggests that shore platforms at Lake Waikaremoana were likely initiated and continue to develop as a result of subaerial wetting and drying cycles. Waves, coincident with fluctuating lake levels, play an important role by removing the weathered material from the platforms, and appear to control the width of the platforms. A conceptual model of platform development is presented, and analogies are drawn between this model, and the formation of shore platforms in oceanic environments. Copyright © 2002 John Wiley & Sons, Ltd.     

The 10 ka multiple vent pyroclastic eruption sequence at Tongariro Volcanic Centre, Taupo Volcanic Zone, New Zealand: Part 2. Petrological insights into magma storage and transport during regional extension

The six eruption episodes of the 10 ka Pahoka–Mangamate (PM) sequence (see companion paper) occurred over a ?200–400-year period from a 15-km-long zone of multiple vents within the Tongariro Volcanic Centre (TgVC), located at the southern end of the Taupo Volcanic Zone (TVZ). Most TgVC eruptives are plagioclase-dominant pyroxene andesites and dacites, with strongly porphyritic textures indicating their derivation from magmas that ascended slowly and stagnated at shallow depths. In contrast, the PM pyroclastic eruptives show petrographic features (presence of phenocrystic and groundmass hornblende, and the coexistence of olivine and augite without plagioclase during crystallisation of phenocrysts and microphenocrysts) which suggest that their crystallisation occurred at depth. Depths exceeding 8 km are indicated for the dacitic magmas, and >20 km for the andesitic and basaltic andesitic magmas. Other petrographic features (aphyric nature, lack of reaction rims around hornblende, and the common occurrence of skeletal microphenocrystic to groundmass olivine in the andesites and basaltic andesites) suggest the PM magmas ascended rapidly immediately prior to their eruption, without any significant stagnation at shallow depths in the crust. The PM eruptives show three distinct linear trends in many oxide–oxide diagrams, suggesting geochemical division of the six episodes into three chronologically-sequential groups, early, middle and late. Disequilibrium features on a variety of scales (banded pumice, heterogeneous glassy matrix and presence of reversely zoned phenocrysts) suggest that each group contains the mixing products of two end-member magmas. Both of these end-member magmas are clearly different in each of the three groups, showing that the PM magma system was completely renewed at least three times during the eruption sequence. Minor compositional diversity within the eruptives of each group also allows the PM magmas to be distinguished in terms of their source vents. Because petrography suggests that the PM magmas did not stagnate at shallow levels during their ascent, the minor diversity in magmas from different vents indicates that magmas ascended from depth through separate conduits/dikes to erupt at different vents either simultaneously or sequentially. These unique modes of magma transport and eruption support the inferred simultaneous or sequential tapping of small separate magma bodies by regional rifting in the southern Taupo Volcanic Zone during the PM eruption sequence (see companion paper).     

Terrestrial heat flow in the North Island of New Zealand

Large variations in terrestrial heat flow from 21 to 209 mW/m² have been observed over the North Island, New Zealand. This is generally in good agreement with the pattern of existing geological and geophysical observations. A high heat flow zone with a value of 92 ±3 mW/m², which corresponds to melting temperatures near the base of the crust, is delineated in the northern part of the Taranaki Basin. In the rest of the island, heat flow appears to be low to normal, but some isolated high values are also found. Observed results are interpreted in terms of crust and mantle structure in a region of plate subduction.     

Geophysical imaging of buried volcanic structures within a continental back-arc basin: The Central Volcanic Region, North Island, New Zealand

Hidden beneath the ~ 2 km thick low-velocity volcaniclastics on the western margin of the Central Volcanic Region, North Island, New Zealand, are two structures that represent the early history of volcanic activity in a continental back-arc. These ~ 20 × 20 km structures, at Tokoroa and Mangakino, form an adjacent gravity high and low, respectively. Interpretations from seismic refraction arrivals and gravity modelling indicate the − 65 mgal Mangakino residual gravity anomaly can be modelled, in part, by two low-density bodies that reach depths of ~ 6.5 km, whereas the Tokoroa gravity anomaly is due to a higher density rock coming, at most, to within ~ 650 m of the surface. The Mangakino anomaly is interpreted to be due to the remnants of magma chambers that fed large ignimbrite eruptions from about 1.2 Ma. An andesite volcano or complex volcanic structure is the preferred interpretation for the Tokoroa gravity high. The size of the putative volcanic structure is comparable to the presently active Tongariro Volcanic Complex in the centre of North Island.