Interpretation of gravity and magnetic anomalies at Lake Rotomahana: geological and hydrothermal implications

Author Posting. © The Author(s), 2015. This is the author's version of the work. It is posted here for personal use, not for redistribution. The definitive version was published in Journal of Volcanology and Geothermal Research 314 (2016): 84-94, doi:10.1016/j.jvolgeores.2015.07.002.We investigate the geological and hydrothermal setting at Lake Rotomahana, using recently collected potential-field data, integrated with pre-existing regional gravity and aeromagnetic compilations. The lake is located on the southwest margin of the Okataina Volcanic Center (Haroharo caldera) and had well-known, pre-1886 Tarawera eruption hydrothermal manifestations (the famous Pink and White Terraces). Its present physiography was set by the caldera collapse during the 1886 eruption, together with the appearance of surface activities at the Waimangu Valley. Gravity models suggest subsidence associated with the Haroharo caldera is wider than the previously mapped extent of the caldera margins. Magnetic anomalies closely correlate with heat-flux data and surface hydrothermal manifestations and indicate that the west and northwestern shore of Lake Rotomahana are characterized by a large, well-developed hydrothermal field. The field extends beyond the lake area with deep connections to the Waimangu area to the south. On the south, the contact between hydrothermally demagnetized and magnetized rocks strikes along a structural lineament with high heat-flux and bubble plumes which suggest hydrothermal activity occurring west of Patiti Island. The absence of a well-defined demagnetization anomaly at this location suggests a very young age for the underlying geothermal system which was likely generated by the 1886 Tarawera eruption. Locally confined intense magnetic anomalies on the north shore of Lake Rotomahana are interpreted as basalts dikes with high magnetization. Some appear to have been emplaced before the 1886 Tarawera eruption. A dike located in proximity of the southwest lake shore may be related to the structural lineament controlling the development of the Patiti geothermal system, and could have been originated from the 1886 Tarawera eruption.Science funding provided by GNS Science Strategic Development Fund

Reactivation of tectonics, crustal underplating, and uplift after 60 Myr of passive subsidence, Raukumara Basin, Hikurangi-Kermadec fore arc, New Zealand: implications for global growth and recycling of continents

We use seismic reflection and refraction data to determine crustal structure, to map a fore-arc basin containing 12 km of sediment, and to image the subduction thrust at 35 km depth. Seismic reflection megasequences within the basin are correlated with onshore geology: megasequence X, Late Cretaceous and Paleogene marine passive margin sediments; megasequence Y, a similar to 10,000 km(3) submarine landslide emplaced during subduction initiation at 22 Ma; and megasequence Z, a Neogene subduction margin megasequence. The Moho lies at 17 km beneath the basin center and at 35 km at the southern margin. Beneath the western basin margin, we interpret reflective units as deformed Gondwana fore-arc sediment that was thrust in Cretaceous time over oceanic crust 7 km thick. Raukumara Basin has normal faults at its western margin and is uplifted along its eastern and southern margins. Raukumara Basin represents a rigid fore-arc block > 150 km long, which contrasts with widespread faulting and large Neogene vertical axis rotations farther south. Taper of the western edge of allochthonous unit Y and westward thickening and downlap of immediately overlying strata suggest westward or northwestward paleoslope and emplacement direction rather than southwestward, as proposed for the correlative onshore allochthon. Spatial correlation between rock uplift of the eastern and southern basin margins with the intersection between Moho and subduction thrust leads us to suggest that crustal underplating is modulated by fore-arc crustal thickness. The trench slope has many small extensional faults and lacks coherent internal reflections, suggesting collapse of indurated rock, rather than accretion of > 1 km of sediment from the downgoing plate. The lack of volcanic intrusion east of the active arc, and stratigraphic evidence for the broadening of East Cape Ridge with time, suggests net fore-arc accretion since 22 Ma. We propose a cyclical fore-arc kinematic: rock moves down a subduction channel to near the base of the crust, where underplating drives rock uplift, oversteepens the trench slope, and causes collapse toward the trench and subduction channel. Cyclical rock particle paths led to persistent trench slope subsidence during net accretion. Existing global estimates of fore-arc loss are systematically too high because they assume vertical particle paths. Citation: Sutherland, R., et al. (2009), Reactivation of tectonics, crustal underplating, and uplift after 60 Myr of passive subsidence, Raukumara Basin, Hikurangi-Kermadec fore arc, New Zealand: Implications for global growth and recycling of continents, Tectonics, 28, TC5017, doi: 10.1029/2008TC002356

Gsolve, a Python computer program with a graphical user interface to transform relative gravity survey measurements to absolute gravity values and gravity anomalies

A Python program (Gsolve) with a graphical user interface has been developed to assist with routine data processing of relative gravity measurements. Gsolve calculates the gravity at each measurement site of a relative gravity survey, which is referenced to at least one known gravity value. The tidal effects of the sun and moon, gravimeter drift and tares in the data are all accounted for during the processing of the survey measurements.The calculation is based on a least squares formulation where the difference between the absolute gravity at each surveyed location and parameters relating to the dynamics of the gravimeter are minimized with respect to the relative gravity observations, and some supplied gravity reference site values. The program additionally allows the user to compute free air gravity anomalies, with respect to the GRS80 and GRS67 reference ellipsoids, from the determined gravity values and calculate terrain corrections at each of the surveyed sites using a prism formula and a user supplied digital elevation model.This paper reviews the mathematical framework used to reduce relative gravimeter survey observations to gravity values. It then goes on to detail how the processing steps can be implemented using the software. Keywords: Relative gravity, Python, Least squares, Free air anomalies, Terrain correction

Zealandia: Earth’s Hidden Continent

A 4.9 Mkm2 region of the southwest Pacific Ocean is made up of continental crust. The region has elevated bathymetry relative to surrounding oceanic crust, diverse and silica-rich rocks, and relatively thick and low-velocity crustal structure. Its isolation from Australia and large area support its definition as a continent— Zealandia. Zealandia was formerly part of Gondwana. Today it is 94% submerged, mainly as a result of widespread Late Cretaceous crustal thinning preceding supercontinent breakup and consequent isostatic balance. The identification of Zealandia as a geological continent, rather than a collection of continental islands, fragments, and slices, more correctly represents the geology of this part of Earth. Zealandia provides a fresh context in which to investigate processes of continental rifting, thinning, and breakup