New petrological, geochemical, and geochronological perspectives on andesite-dacite magma genesis at Ruapehu volcano, New Zealand

Time–composition relationships in eruptive sequences at composite volcanoes can show how the ongoing intrusion of magmas progressively affects the lithosphere at continental convergent margins. Here, new whole-rock and microanalytical major and trace element data from andesite-dacite lava flows are integrated with previous studies and existing isotopic data, and placed within the framework of a high-resolution chronostratigraphy for Ruapehu volcano (southern Taupo Volcanic Zone, New Zealand). The geochemical evolution of lavas erupted over the ∼200 kyr lifetime of the exposed edifice reflects variable degrees of fractionation and systematic changes in the type of crustal assimilation in the Ruapehu magma system. Lavas erupted from ∼200–150 ka have previously been distinguished from those erupted <150 ka based on Sr-Nd isotopic characteristics, which indicate that the oldest lavas were sourced from magmas that assimilated oceanic crust. Such source rocks underlie the regionally widespread Mesozoic meta-sedimentary greywacke-argillite basement, which was conversely assimilated by <150 ka magmas. New results from this work reveal that since 150 ka, an upper limit of magma differentiation occurred from ∼50–35 ka. High K2O (∼6 wt%) and Rb contents (∼270 ppm) in melt inclusions, interstitial glass, and glass from in situ quenched melts of partially fused crustal xenoliths are reported for andesite-dacite lavas erupted during this period. In addition to crystal fractionation, selective partial melting and assimilation of K- and Rb-rich mineral phases (e.g., biotite, K-feldspar) that are significant components of the meta-sedimentary basement rocks is inferred to explain these geochemical characteristics. These processes coincided also with the effusion of high-MgO andesitedacite lavas that display petrological evidence for mixing between andesite-dacite and more mafic magmas. An influx of hotter mafic magma into the system explains why the extent of crustal assimilation recorded by Ruapehu lavas peaked during the ∼50–35 ka eruptive period. From 26 ka to the present, andesite lavas have reverted to more mafic compositions with less potassic melt inclusion and whole-rock compositions when compared to the ∼50–35 ka lavas. We suggest that the younger lavas assimilated less-enriched melts because fertile phases had been preferentially extracted from the crustal column during earlier magmatism. This scenario of bottom-up heating of the lithosphere and exhaustion of fertile phases due to the progressive intrusion of magma explains the geochemical evolution of Ruapehu lavas. This model may be applicable to other long-lived composite volcanoes of the circum-Pacific continental arcs

The Last Glacial Maximum in the central North Island, New Zealand: palaeoclimate inferences from glacier modelling

Abstract. Quantitative palaeoclimate reconstructions provide data for evaluating the mechanisms of past, natural climate variability. Geometries of former mountain glaciers constrained by moraine mapping afford the opportunity to reconstruct palaeoclimate, due to the close relationship between ice extent and local climate. In this study, we present results from a series of experiments using a 2D coupled energy-balance/ice-flow model that investigate the palaeoclimate significance of Last Glacial Maximum moraines within nine catchments in central North Island, New Zealand. We find that the former ice limits can be simulated when present day temperatures are reduced by between 4 °C and 7 °C, when precipitation remains unchanged from present. The spread in the results between the nine catchments is likely to represent the combination of chronological and model uncertainties. The temperature decrease required to simulate the former glaciers falls in the range of 5.1 °C and 6.3 °C for the majority of catchments targeted, which represents our best estimate of the peak temperature anomaly in central North Island, New Zealand during the Last Glacial Maximum. A decrease in precipitation, as suggested by proxy evidence and climate models, of up to 25 % from present, increases the magnitude of the required temperature changes by up to 0.8 °C. Glacier model experiments using reconstructed topographies that exclude the volume of post-glacial (&lt;15 ka) volcanism, generally increased the magnitude of cooling required to simulate the former ice limits by up to 0.5 °C. Our palaeotemperature estimates expand the spatial coverage of proxy-based quantitative palaeoclimate reconstructions in New Zealand, and are consistent with independent, proximal temperature reconstructions from fossil pollen assemblages, as well as similar glacier modelling reconstructions from central Southern Alps. </jats:p

Surface rupture of multiple crustal faults in the 2016 Mw 7.8 Kaikōura, New Zealand, earthquake

Multiple (>20 >20 ) crustal faults ruptured to the ground surface and seafloor in the 14 November 2016 M w Mw 7.8 Kaikōura earthquake, and many have been documented in detail, providing an opportunity to understand the factors controlling multifault ruptures, including the role of the subduction interface. We present a summary of the surface ruptures, as well as previous knowledge including paleoseismic data, and use these data and a 3D geological model to calculate cumulative geological moment magnitudes (M G w MwG ) and seismic moments for comparison with those from geophysical datasets. The earthquake ruptured faults with a wide range of orientations, sense of movement, slip rates, and recurrence intervals, and crossed a tectonic domain boundary, the Hope fault. The maximum net surface displacement was ∼12  m ∼12  m on the Kekerengu and the Papatea faults, and average displacements for the major faults were 0.7–1.5 m south of the Hope fault, and 5.5–6.4 m to the north. M G w MwG using two different methods are M G w MwG 7.7 +0.3 −0.2 7.7−0.2+0.3 and the seismic moment is 33%–67% of geophysical datasets. However, these are minimum values and a best estimate M G w MwG incorporating probable larger slip at depth, a 20 km seismogenic depth, and likely listric geometry is M G w MwG 7.8±0.2 7.8±0.2 , suggests ≤32% ≤32% of the moment may be attributed to slip on the subduction interface and/or a midcrustal detachment. Likely factors contributing to multifault rupture in the Kaikōura earthquake include (1) the presence of the subduction interface, (2) physical linkages between faults, (3) rupture of geologically immature faults in the south, and (4) inherited geological structure. The estimated recurrence interval for the Kaikōura earthquake is ≥5,000–10,000  yrs ≥5,000–10,000  yrs , and so it is a relatively rare event. Nevertheless, these findings support the need for continued advances in seismic hazard modeling to ensure that they incorporate multifault ruptures that cross tectonic domain boundaries

Neogene Evolution of the Pacific - Australia Plate Boundary Zone in NE Marlborough, South Island, New Zealand

Six new palaeomagnetic localities in NE Marlborough, sampled from Late Cretaceous - Early Tertiary Amuri Formation and Middle Miocene Waima Formation, all yield clockwise declination anomalies of 100 - 150 degrees. Similarity in the magnitude of all new declination anomalies and integration of these results with previous data implies that clockwise vertical-axis rotation of this magnitude affected the entire palaeomagnetically sampled part of NE Marlborough (an area of ~700sq. km) after ~18 Ma. Previous palaeomagnetic sampling constrains this rotation to have occurred before ~7 Ma. The regional nature of this rotation implies that crustal-scale vertical-axis rotations were a fundamental process in the Miocene evolution of the Pacific - Australia plate boundary in NE South Island. The Flags Creek Fault System (FCFS) is a fold-and-thrust belt that formed in marine conditions above a subduction complex that developed as the Pacific - Australia plate boundary propagated through Marlborough in the Early Miocene. Thin-skinned fault offset accommodated at least 20 km of horizontal shortening across a leading-edge imbricate fan. Mesoscopic structures in the deformed belt indicate thrust vergence to the southeast. The palaeomagnetically-determined regional clockwise vertical axis rotation of ~100 degrees must be undone in order to evaluate this direction in the contemporary geographic framework of the thrust belt. Therefore the original transport direction of the thrust sheets in the FCFS was to the NE, in accordance with NE-SW plate motion vector between the Pacific and Australian plates during the Early Miocene. The two new palaeomagnetic localities that are within ~3 km of the active dextral strike-slip Kekerengu Fault have the highest clockwise declination anomalies (up to 150 degrees). Detailed structural mapping suggests that the eastern ends of the FCFS are similarly clockwise-rotated, by an extra 45 degrees relative to the regional average, to become south-vergent in proximity to the Kekerengu Fault. This structural evidence implies the presence of a zone of Plio-Pleistocene dextral shear and vertical-axis rotation within 2-3 km of the Kekerengu Fault. Local clockwise vertical-axis rotations of up to 50 degrees are inferred to have accrued in this zone, and to have been superimposed on the older, regional. ~100 degrees Miocene clockwise vertical-axis rotation. The Late Quaternary stratigraphy of fluvial terraces in NE Marlborough has been revised by the measurement of five new optically stimulated luminescence (OSL) dates on loess. This new stratigraphy suggests that the latest aggradation surface in the Awatere Valley (the Starborough-1 terrace) is, at least locally, ~9 ka old, several thousand years younger than the previous 16 ka thermoluminescence age for the same site. This new surface abandonment age implies that terrace-building events in NE Marlborough lasted well after the last glacial maximum (~17 ka). The timing of terrace aggradation in this peri-glacial region is compared with oxygen isotope data. Downstream transport of glacially derived sediment at the time of maximum deglaciation/warming is concluded to be the primary influence on the aggradation of major fill terraces in coastal NE Marlborough. This interpretation is generally applicable to peri-glacial central New Zealand. Patterns of contemporary uplift and directions of landscape tilting have been analysed by assessing the rates of stream incision and by the evolution of drainage networks over a wide tract of NE Marlborough that includes the termination of the dextral strike-slip Clarence Fault. Relative elevations of differentially aged terraces suggests an increase in rates of incision over the last ~10 ka. Uplift is highest in the area immediately surrounding the fault tip and is generally high where Torlesse basement rocks are exposed. Independently derived directions of Late Quaternary tilting of the landscape display a similar pattern of relative uplift in a broad dome to the north and west of the fault tip. This pattern of uplift suggests dissipation of strike-slip motion at the Clarence Fault tip into a dome-shaped fold accommodating: 1) crustal thickening (uplift) and 2) up to 44 degrees of vertical-axis rotation of a ~40 km2 crustal block, relative to more inland domains, into which the fault terminates. The distribution of incision rates is compared with the pattern of crustal thickening predicted by elastic models of strike-slip fault tips. The observed pattern and spatial extent of uplift generally conforms with the distribution of thickening predicted by the models, although the rate of incision/uplift over the last ~120 ka has been variable. These differences may be due to variability in the strike-slip rate of the Clarence Fault, superimposition of the regional uplift rate or to interaction with nearby fault structures not accounted for in the models

Volcanic records of the Laschamp geomagnetic excursion from Mt Ruapehu, New Zealand

We present palaeodirectional records of the Laschamp geomagnetic excursion from lavas on Mt Ruapehu, New Zealand. Fourteen lava flows on the northwestern and southern flanks of Mt Ruapehu, with 40Ar/39Ar weighted mean plateau ages that range from 46.3±2.0 to 39.9±1.4 ka, were studied. The youngest and older flows carry a normal polarity magnetization; however, six flows, dated between 46.3±2.0 and 42.7±1.8 ka, record excursional directions. Three of these flows record southerly palaeomagnetic declinations and negative inclinations that agree well with a published Laschamp record from the Auckland Volcanic Field (AVF). Together, the AVF and Mt Ruapehu lavas currently represent the only volcanic records of the Laschamp excursion outside the Chaîne des Puys region, France. Thus, they make an important contribution to the global set of Laschamp excursion records. Virtual geomagnetic pole (VGP) groups for the New Zealand and French records early in the excursion are compatible with a dipole-dominated field that rotated to an equatorial orientation while simultaneously decaying in strength. In contrast, younger excursional flows from France and New Zealand yield separate VGP groups, which suggest either that the field had a nondipolar morphology in this later phase, or that the VGP groups were not synchronous. 40Ar/39Ar ages for the Mt Ruapehu record are on average slightly older than published northern hemisphere ages and from the relative palaeointensity minimum in the GLOPIS sedimentary stack. Although few individual ages differ significantly at the 2σ level, the spread suggests an overall excursion duration that is longer than the currently accepted 1500 years. This age spread may result from excess Ar in magmas at the time of the eruption biasing the results to slightly older ages, or from non-synchronous excursional field behaviour at near-antipodal locations, or, possibly, a precursory phase prior to the main excursio

Chemical and isotopic changes induced by pyrometamorphism in metasedimentary xenoliths at Tongariro volcano, New Zealand

Andesites erupted from Tongariro volcano, North Island, New Zealand contain feldspathic and quartzose xenoliths derived from basement rocks. New major oxide, trace element and Sr-Nd-Pb isotopic data indicate that both the Waipapa and Kaweka (meta)sedimentary terranes are represented in erupted xenoliths, rather than only the Kaweka terrane as previously thought. Xenolith mineral assemblages differ from their likely source materials, notably through the lack of white mica, illite, chlorite and quartz, which is reflected in contrasting chemical and isotopic compositions. Major and trace element data indicate that most xenoliths underwent bulk mass decreases of about 50% when pyrometamorphosed at temperatures of ~800–980 °C, similar to typical Tongariro magma temperatures of ~800–1000 °C. Bulk Eu concentrations were retained (in restitic plagioclase); however, other rare earth elements are commonly lower in xenoliths than in protoliths. In xenoliths, the 143Nd/144Nd ratios of protoliths were also retained, which indicates that xenoliths were derived from the Kaweka and Waipapa terranes in subequal amounts. Reductions in 87Sr/86Sr ratios by up to 0.003 in xenoliths, relative to their likely protoliths, were accompanied by decreases in Rb/Sr ratios from 0.1–0.8 down to <0.1, reflecting the dissolution of hydrous, Rb-rich minerals (white mica ± illite) with their radiogenic isotopic ingrowths liberated into surrounding andesitic magmas. Varied amounts of U/Pb, Th/Pb and Th/U fractionation demonstrably occurred between xenoliths and protoliths, but these are challenging to correlate with Pb isotopic fractionation that also occurred. One xenolith contains a vein of clinopyroxene, calcic plagioclase, silicic glass and graphite that formed when quartz + calcite veins were pyrometamorphosed. The vein-bearing xenolith possesses unusual chemical and isotopic features, which include a negative Ce anomaly, LaN/YbN ~ 2 and high 143Nd/144Nd (0.51284), which are also reported for xenolithic material erupted at neighbouring Ruapehu volcano

A high resolution 40Ar/39Ar lava chronology and edifice construction history for Tongariro volcano, New Zealand

Detailed mapping and geochronological investigations of edifice-forming materials reconstruct the growth history of Tongariro volcano, New Zealand, and subdivide the edifice into thirty six distinct units which are organised into twelve formations and constituent members. Twenty nine new 40Ar/39Ar age determinations, along with published K/Ar ages combined with volume estimates, petrographic observations and rock chemistry provide an integrated history of the volcano's growth through edifice-forming lavas and pyroclastic deposits. The oldest lava (512 ± 59 ka, 2 s.d.) is a small inlier of basaltic-andesite on Tongariro's NW sector that may reflect a nearly buried independent volcano. The next oldest material that can be confidently attributed to a Tongariro source is 304 ± 11 ka andesite, incorporated as boulders in late Pleistocene ejecta from the Tama Lakes area. In-situ lavas at Tongariro date from 230 ka to present, including numerous flows erupted during glacial periods and building the edifice unevenly due to emplacement against valley-filling ice bodies. Tongariro has a total edifice volume of ~90 km3, 19 km3 of which is represented by exposed map units, with glacial deposits amounting to <1 km3. The ring plain volume immediately adjacent to the volcano contains ~60 km3 of material. Sequential eruptive records, from 230 ka to present, reveal an irregular cyclicity in MgO concentrations over ~10–70 kyr intervals. During these cycles, rapid (≤10 kyr) increases in MgO concentrations to ≥5–9 wt% are inferred to reflect episodes of enhanced mafic magma replenishment, with maxima at ~230, ~160, ~117, ~88, ~56, ~35, ~17.5 ka and during the Holocene, which are each followed by gradual declines to ~2–5 wt%. Field evidence, including extensive moraines and U-shaped valleys, and lava textures, implies repeated occupation of valleys on Tongariro by major glaciers and possibly ice caps. During periods of major ice coverage, which generally correlate with global cold climate/glacial Marine Isotope Stages, edifice-building rates on Tongariro were only 17–26% of those during warmer climatic periods. Because the changes in edifice-building rates do not coincide with changes in the magmatic system, these contrasts are inferred to reflect a preservation bias whereby materials erupted onto ice were contemporaneously (or subsequently, as ice masses melted) conveyed to the ring plain as debris rather than building the edifice. Although the Tongariro edifice is smaller than that of neighbouring Ruapehu (~150 km3), the exposed edifice materials on Tongariro record a longer and more complex growth history. The wider geographic distribution of <50 ka vent locations at Tongariro reflects greater rifting rates than at Ruapehu

Ruapehu and Tongariro stratovolcanoes: a review of current understanding

Ruapehu (150 km3 cone, 150 km3 ring-plain) and Tongariro (90 km3 cone, 60 km3 ring-plain) are iconic stratovolcanoes, formed since ∼230 and ∼350 ka, respectively, in the southern Taupo Volcanic Zone and Taupo Rift. These volcanoes rest on Mesozoic metasedimentary basement with local intervening Miocene sediments. Both volcanoes have complex growth histories, closely linked to the presence or absence of glacial ice that controlled the distribution and preservation of lavas. Ruapehu cone-building vents are focused into a short NNE-separated pair, whereas Tongariro vents are more widely distributed along that trend, the differences reflecting local rifting rates and faulting intensities. Both volcanoes have erupted basaltic andesite to dacite (53–66 wt.% silica), but mostly plagioclase-two pyroxene andesites from storage zones at 5–10 km depth. Erupted compositions contain evidence for magma mixing and interaction with basement rocks. Each volcano has an independent magmatic system and a growth history related to long-term (>104 years) cycles of mantle-derived magma supply, unrelated to glacial/interglacial cycles. Historic eruptions at both volcanoes are compositionally diverse, reflecting small, dispersed magma sources. Both volcanoes often show signs of volcanic unrest and have erupted with a wide range of styles and associated hazards, most recently in 2007 (Ruapehu) and 2012 (Tongariro)