J. V. Rowland

Evidence for focused magmatic accretion at segment centers from lateral dike injections captured beneath the Red Sea rift in Afar

Continental breakup occurs through repeated episodes of mechanical stretching and dike injection within discrete, narrow rift segments. However, the time and length scales of the dike intrusions, along with the source regions of melt within continental and oceanic rifts, are poorly constrained. We present measurements of spatial and temporal variability in deformation from the currently active 60-km-long Dabbahu segment of the Red Sea rift in Afar, using satellite radar, global positioning system, and seismicity data sets, that capture emplacement of two ~10-km-long, ~1–2-m-wide dike intrusions in June and July 2006. Our observations show that the majority of strain is accommodated by dikes that propagate laterally over ~4–5 h time scales along the rift axis and are sourced from a reservoir in the middle to lower crust, or upper mantle, beneath the center of the rift segment. New intrusions during the ongoing rifting episode in Afar show that the injection of lateral dikes fed from magma reservoirs beneath rift segment centers is a key component in creating and maintaining regular along-axis rift segmentation during the final stages of continental breakup. Our observations also provide evidence that the focused magmatic accretion at segment centers observed in slow-spreading mid-ocean ridges occurs prior to the onset of seafloor spreading.

Complex multifault rupture during the 2016 Mw 7.8 Kaikōura earthquake, New Zealand

An earthquake with a dozen faults The 2016 moment magnitude (Mw) 7.8 Kaikōura earthquake was one of the largest ever to hit New Zealand. Hamling et al. show with a new slip model that it was an incredibly complex event. Unlike most earthquakes, multiple faults ruptured to generate the ground shaking. A remarkable 12 faults ruptured overall, with the rupture jumping between faults located up to 15 km away from each other. The earthquake should motivate rethinking of certain seismic hazard models, which do not presently allow for this unusual complex rupture pattern. Science, this issue p. eaam7194 At least 12 faults spaced up to 15 kilometers apart ruptured during the magnitude 7.8 Kaikōura earthquake. INTRODUCTION On 14 November 2016 (local time), northeastern South Island of New Zealand was struck by a major moment magnitude (Mw) 7.8 earthquake. The Kaikōura earthquake was the most powerful experienced in the region in more than 150 years. The whole of New Zealand reported shaking, with widespread damage across much of northern South Island and in the capital city, Wellington. The earthquake straddled two distinct seismotectonic domains, breaking multiple faults in the contractional North Canterbury fault zone and the dominantly strike-slip Marlborough fault system. RATIONALE Earthquakes are conceptually thought to occur along a single fault. Although this is often the case, the need to account for multiple segment ruptures challenges seismic hazard assessments and potential maximum earthquake magnitudes. Field observations from many past earthquakes and numerical models suggest that a rupture will halt if it has to step over a distance as small as 5 km to continue on a different fault. The Kaikōura earthquake’s complexity defies many conventional assumptions about the degree to which earthquake ruptures are controlled by fault segmentation and provides additional motivation to rethink these issues in seismic hazard models. RESULTS Field observations, in conjunction with interferometric synthetic aperture radar (InSAR), Global Positioning System (GPS), and seismology data, reveal the Kaikōura earthquake to be one of the most complex earthquakes ever recorded with modern instrumental techniques. The rupture propagated northward for more than 170 km along both mapped and unmapped faults before continuing offshore at the island’s northeastern extent. A tsunami of up to 3 m in height was detected at Kaikōura and at three other tide gauges along the east coast of both the North and South Islands. Geodetic and geological field observations reveal surface ruptures along at least 12 major crustal faults and extensive uplift along much of the coastline. Surface displacements measured by GPS and satellite radar data show horizontal offsets of ~6 m. In addition, a fault-bounded block (the Papatea block) was uplifted by up to 8 m and translated south by 4 to 5 m. Modeling suggests that some of the faults slipped by more than 20 m, at depths of 10 to 15 km, with surface slip of ~10 m consistent with field observations of offset roads and fences. Although we can explain most of the deformation by crustal faulting alone, global moment tensors show a larger thrust component, indicating that the earthquake also involved some slip along the southern end of the Hikurangi subduction interface, which lies ~20 km beneath Kaikōura. Including this as a fault source in the inversion suggests that up to 4 m of predominantly reverse slip may have occurred on the subduction zone beneath the crustal faults, contributing ~10 to 30% of the total moment. CONCLUSION Although the unusual multifault rupture observed in the Kaikōura earthquake may be partly related to the geometrically complex nature of the faults in this region, this event emphasizes the importance of reevaluating how rupture scenarios are defined for seismic hazard models in plate boundary zones worldwide. Observed ground deformation from the 2016 Kaikōura, New Zealand, earthquake. (A and B) Photos showing the coastal uplift of 2 to 3 m associated with the Papatea block [labeled in (C)]. The inset in (A) shows an aerial view of New Zealand. Red lines denote the location of known active faults. The black box indicates the Marlborough fault system

Interconnected sills and inclined sheet intrusions control shallow magma transport in the Ferrar large igneous province, Antarctica

Field observations and structural data from intrusive complexes at Allan Hills and Terra Cotta Mountain, South Victoria Land, Antarctica, demonstrate that interconnected sills and inclined sheets transported magma through the shallow subsurface. These sills and sheets represent the upper-crustal (top 4 km) plumbing system of the 183 Ma Ferrar large igneous province. The sheets are short in length (

Evolution of the intra-arc Taupo-Reporoa Basin within the Taupo Volcanic Zone of New Zealand

The spatial and temporal distributions of volcaniclastic deposits in arc-related basins reflect a complex interplay between tectonic, volcanic, and magmatic processes that is typically difficult to unravel. We take advantage of comprehensive geothermal drill hole stratigraphic records within the Taupo-Reporoa Basin (TRB), and integrate them with new 40 Ar/ 39 Ar age determinations, existing age data, and new mapping to develop a four-dimensional model of basin evolution in the central Taupo Volcanic Zone (TVZ), New Zealand. Here, exceptional rhyolitic productivity and high rates of extensional tectonism have resulted in the formation of at least eight calderas and two subparallel, northeast-trending rift basins, each of which is currently subsiding at 3 to 4 mm/yr: the Taupo fault belt (TFB) to the northwest and the TRB to the southeast (the main subject of this paper). The basins are separated in the northeast by a high-standing, fault-controlled range termed the Paeroa block, which is the focus of mapping for this study, and in the southwest by an along strike alignment of smaller scale faults and an associated region of lower relief. Stratigraphic age constraints within the Paeroa block indicate that a single basin (∼120 km long by 60 km wide) existed within the central TVZ until 339 ± 5 ka (Paeroa Subgroup eruption age), and it is inferred to have drained to the west through a narrow and deep constriction, the present-day Ongaroto Gorge. Stratigraphic evidence and field relationships imply that development of the Paeroa block occurred within 58 ± 26 k.y. of Paeroa Subgroup emplacement, but in two stages. The northern Paeroa block underwent uplift and associated tilting first, followed by the southern Paeroa block. Elevations (>500 m above sea level) of lacustrine sediments within the southern Paeroa block are consistent with elevations of rhyolite lavas in the Ongaroto Gorge, the outlet to the paleolake in which these sediments were deposited, and indicate that the Paeroa block has remained relatively stable since development. East of the Paeroa block, stratigraphic relationships show that movement along the Kaingaroa Fault zone, the eastern boundary of the central TVZ, is associated with volcano-tectonic events. Stratigraphic and age data are consistent with rapid formation of the paired TRB and TFB at 339 ± 5 ka, and indicate that gradual, secular rifting is punctuated by volcano-tectonic episodes from time to time. Both processes influence basin evolution.

Oblique rifting along the central and southern Kermadec Arc front (30°–36°S), SW Pacific

We present a series of eleven structural maps of the central Kermadec arc front between 30°S and 35°50S interpreted from multibeam bathymetric data. These maps show the structural morphology of eighteen submarine arc volcanoes proximal to the Havre backarc, providing a critical window into the tectonomagmatic processes that operate at divergent arcs transitional to seafloor spreading. Each arc volcano, unaffected by subaerial processes, shows the structural imprint of oblique backarc rifting, providing an exceptional opportunity to evaluate tectonic influences on edifice construction and destruction. Mapped seafloor structures include faults, volcanic ridges, aligned vents and cones, cones aligned with faults, elongate calderas, and structurally controlled sector collapse. In addition to the structural maps, we present the results of quantitative analyses of fault orientation and growth behaviors. We demonstrate that >50% of identified faults strike NE-SW between 030° and 060°, with a modal value of 047°, evidence in support of an extension axis dextrally oblique to the NNE trend of the rift. Further to this, minimum and maximum values for the extension axis are calculated as 137° and 150°, perpendicular to the modal strikes for all identified faults and for faults within the most oblique surveyed backarc rift domain, respectively. These values confirm moderate obliquity of opening along the Havre rift, with the SHT opening slightly less obliquely than the NHT, consistent with accommodation of a lesser partitioned component of dextral Pacific-Australian plate convergence.

Energetics of normal earthquakes on dip-slip faults

That earthquakes release vast quantities of energy is widely accepted; however, the most commonly experienced component, radiated seismic energy, is a minor contribution to the total energy budget. The elastic rebound model for earthquakes recognizes that elastic strain energy does work displacing, deforming, and accelerating the crust, as well as causing frictional heating. In this paper we present an energy budget for dip-slip fault rupture in an extensional tectonic regime. A computational model of an elastic-plastic-viscous crust hosting a single fault, modeled as two surfaces in frictional contact, demonstrates contrasting energy flows between the hanging-wall and footwall fault blocks. Our analysis suggests that in the period leading up to an earthquake, the total strain energy contained within the crust decreases, although a local increase within the footwall at mid-crustal depths is observed. During an earthquake, the footwall is subject to an elastic rebound, whereupon uplift of the fault scarp is caused by a mid-crustal stress drop and elastic expansion that releases strain energy. In contrast, gravitational potential energy released from a subsiding hanging wall does work compressing the wider crust, particularly in the mid-crust at the fault tip. This has the unusual consequence of increasing strain energy throughout much of the upper crust during an earthquake. These counterintuitive energy flows suggest that extensional deformation is caused by stored gravitational potential and elastic strain energy, and not by the external tectonic forcing.

Emplacement of magma at shallow depth: insights from field relationships at Allan Hills, south Victoria Land, East Antarctica

Abstract Allan Hills nunatak, south Victoria Land, Antarctica, exposes an exceptional example of a shallow depth (< 500 m) intrusive complex formed during the evolution of the Ferrar large igneous province (LIP). Dyke distribution, geometries and relationships allow reconstruction of its history and mechanics of intrusion. Sills interconnect across host sedimentary layers, and a swarm of parallel inclined dolerite sheets is intersected by a radiating dyke-array associated with remnants of a phreatomagmatic vent, where the dolerite is locally quenched and mixed to form peperite. Intrusion geometries, and lack of dominant rift-related structures in the country rock indicate that magma overpressure, local stresses between mutually interacting dykes and vertical variations of host rock mechanical properties controlled the intrusive process throughout the thick and otherwise undeformed pile of sedimentary rocks (Victoria Group). Dolerite sills connected to one another by inclined sheets are inferred to record the preferred mode of propagation for magma-carrying cracks that represent the shallow portions of the Ferrar LIP plumbing system.

Rapid priming, accumulation, and recharge of magma driving recent eruptions at a hyperactive caldera volcano

A major challenge in volcanology is determining the factors that control the frequency and magnitude of eruptions at hazardous caldera volcanoes. Understanding the critical sequence of events that may lead to future eruptions is vital for volcanic monitoring and risk assessment. Here we use magma chemistry and mineral diffusion modeling to interpret the magmatic processes and time scales involved in the youngest three eruptions (2.15–1.7 ka) from Taupo volcano (New Zealand), which peaked with the voluminous A.D. 232 eruption. Of the rhyolites erupted since ca. 12 ka, the 2.75 ka) rhyolites that were tapped from the same magma reservoir. Orthopyroxene Fe-Mg diffusion time scales indicate that the onset of rapid heating and priming of the host silicic mush occurred <120 yr prior to the <2.15 ka eruptions, with subsequent melt accumulation occurring in only decades. Elevated mafic magma supply to the silicic mush pile, rapid melt accumulation, and high differential tectonic stress built up and culminated in the ∼105 km3 A.D. 232 eruption, one of the largest and most violent Holocene eruptions globally. These youngest eruptions demonstrate how Taupo’s magmatic system can rapidly change behavior to generate large eruptible melt bodies on time scales of direct relevance to humans and monitoring initiatives.

Age and eruptive center of the Paeroa Subgroup ignimbrites (Whakamaru Group) within the Taupo Volcanic Zone of New Zealand

We here explore the temporal and spatial relationships between the contrasting sources for two eruptive episodes that collectively represent the Whakamaru Group, the largest ignimbrite-forming sequence in the ∼2 m.y. history of the Taupo Volcanic Zone in New Zealand. At 349 ± 4 ka (weighted mean at 2σ), the >1500 km 3 widespread Whakamaru Group ignimbrites and ∼700 km 3 Rangitawa Tephra fallout were erupted in association with collapse of the 40 km long by 25 km wide rectilinear Whakamaru caldera. New 40 Ar/ 39 Ar age data presented here show that the co-magmatic >110 km 3 Paeroa Subgroup ignimbrites, previously included as part of the Whakamaru Group, are slightly younger and were erupted at 339 ± 5 ka (weighted mean at 2σ). New field evidence also presented here demonstrates that the Paeroa Subgroup ignimbrites came from a source geographically separated from vents for the widespread Whakamaru Group ignimbrites. The presence of co-ignimbrite lag breccias, sizes of vent-derived lithic clasts, thicknesses of exposed and subsurface deposits, and morphologies of deposits imply that eruptions of the Paeroa Subgroup occurred from a linear source (the Paeroa linear vent zone), coinciding with the present-day northeast-striking Paeroa fault, and outside (northeast) of the earlier Whakamaru caldera collapse area. No separate caldera has been recognized, although three nearby areas may have undergone eruption-related subsidence. Residual magma from the Whakamaru or adjacent Kapenga caldera areas may have migrated toward the Paeroa linear vent zone during eruptive episodes, resulting in subsidence in either, or both, of these areas. Shallow plutons are also inferred to lie beneath near source fault blocks (Paeroa and Te Weta) on each side of the fault, and eruption-related subsidence may have been expressed as movement across the Paeroa fault and localized subsidence in the southern Paeroa fault block. Subsequent secular, rift-related displacement along the Paeroa fault has obscured the Paeroa linear vent zone.

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.