Search

Gravity surveys using relative gravity meters are often tied to accessible, accurate and stable gravity reference stations. Since 1947 gravity reference networks have included pendulum gravity stat...

The Taupō Volcanic Zone (TVZ) is characterised by a negative residual gravity anomaly that correlates with a zone of normal faulting, subsidence and voluminous silicic volcanism from c. 2 Ma. A ste...

Downhole data and cores collected during International Ocean Discovery Program (IODP) Expedition 376 at Brothers volcano, Kermadec arc, provide unprecedented, in situ views of volcanic facies and fluid pathways in an actively forming volcanogenic massive sulfide (VMS) ore deposit. Brothers volcano is a submarine caldera with extensive sea floor hydrothermal alteration. Downhole data were collected in two holes: Hole U1530A at the NW Caldera and Hole U1528D at the Upper Cone. Textural analysis of microresistivity images in Hole U1530A provides a continuous image facies record that greatly improves findings based upon sporadic and partial (18%) core recovery. Between 90 and 214 meters below sea floor (mbsf), the heterogeneous image facies with local pattern variations is consistent with the volcaniclastic facies interpreted from cores. Between 232 and 445 mbsf, a volcanic facies was not recognizable in cores because of overprinting alteration, apart from five intervals of coherent lava flows that were less altered. Based on the fairly constant petrophysical data, Vp-porosity relationship, and presence of five to six coherent image facies intervals on the microresistivity image, we propose that the apparent volcaniclastic textures observed on cores and microresistivity images beneath 232 mbsf are dominantly lava flows. The change from volcaniclastic to dominant lava flow facies occurs over a transition zone (214–232 mbsf) where all petrophysical properties gradually change. In Hole U1528D, cores and petrophysical data show a similar transition from deep coherent lava flows to shallower, largely volcaniclastic sequences at ~270 mbsf. Down to 232 mbsf in Hole U1530A and 360 mbsf in Hole U1528D, the overall first-order downward decrease in porosity is interpreted to be caused by compaction and increased alteration intensity. Volcanic facies and fractures exert a second-order local control on petrophysical properties. Beneath 232 mbsf in Hole U1530A, the prolonged hydrothermal activity is inferred to have diminished local petrophysical property variations within the proposed lava flow-dominated rock package. High downhole fluid temperatures in Hole U1528D contrast with the moderate temperatures in Hole U1530A. Permeable zones show a mix of structural (inferred fault in Hole U1530A) and lithological controls in both holes. Some low-permeability layers and/or lithological interfaces possibly focus fluids laterally in higher-permeability layers, which may act as a trap for metal-rich fluids to form stratabound massive sulfides and deposits. Matrix is likely too low in permeability to conduct fluids but provides perfect conditions for the storage of super saline brines. In Hole U1530A, located near active vents at the sea floor, the relatively low fluid temperature and the alteration overprint of moderate temperature demonstrate the high spatial and temporal variations at Brothers volcano. The implications of the new stratigraphy and controls on permeability proposed here for Brothers volcano include a better understanding of the following: 1) submarine volcanic eruption sequences, 2) permeability in active submarine volcanoes, and 3) the formation of volcanogenic massive sulfide deposits on (and near) the sea floor.

We show a simple and fast method for calculating geometric parameters of magnetic contacts from spatial gradients of magnetic field data. The method is based on well-established properties of the tangent of the tilt-angle of reduced-to-the-pole magnetic data, and extends the performance of existing methods by allowing direct estimation of depths, locations and dips of magnetic contacts. It uses a semi-automatic approach where the user interactively specifies points on magnetic maps where the calculation is to be performed. Some prior geologic knowledge and visual interpretation of magnetic anomalies is required to choose proper calculation points. We successfully tested the method on synthetic models of contacts at different depths and with different dip angles. We offer an example of the method applied to airborne magnetic data from Taranaki Basin located offshore the North Island of New Zealand.

The Hinepuia volcanic center is made up of two distinct edifices aligned northwest to southeast, with an active cone complex in the SE. Hinepuia is one of several active volcanoes in the northern segment of the Kermadec arc. Regional magnetic data shows no evidence for large-scale hydrothermal alteration at Hinepuia, yet plume data confirm present-day hydrothermal discharge, suggesting that the hydrothermal system may be too young to have altered the host rocks with respect to measurable changes in magnetic signal. Gravity data are consistent with crustal thinning and shallow mantle under the volcanic center. Following the discovery of hydrothermal plumes over Hinepuia, the submersible Shinkai 6500 was used to explore the SE cone and sample hydrothermal fluids. The chemistry of hydrothermal fluids from submarine arc and backarc volcanoes are typically dominated by water-rock interactions and/or magmatic degassing. Chemical analyses of vent fluids show that Hinepuia does not quite fit either traditional model. Moreover, the Hinepuia samples fall between those typically ascribed to both end-member fluid types when plotted on a K-Mg-SO4 ternary diagram. Due to evidence of strong degassing, abundant native sulfur deposition, and H2S presence, the vent sampled at Hinepuia is ultimately classified as a magmatic-hydrothermal system with a water-rock influence. This vent is releasing water vapor and magmatic volatiles with a notable lack of salinity due to subcritical boiling and phase separation. Magmatic-hydrothermal fluid chemistry appears to be controlled by a combination of gas flux, phase separation processes, and volcano evolution and/or distance from the magma source.

Little is known about the effects that subducting an oceanic large igneous province (LIP) has on the petrogenesis of submarine arc volcanoes and their geochemical composition. The southern Kermadec arc represents a rare example where an LIP—the Hikurangi Plateau—is currently subducting and where its effect on mantle composition, element recycling and arc volcanism can be studied. We present mineral chemistry and whole-rock major and trace element, and Sr–Nd–Pb isotope data from samples recovered from the southern Kermadec arc volcanoes Rumble II East and Rumble II West, together with shipboard gravity and magnetic measurements. The Rumble II volcanoes (including a volcanic cone ∼10 km further west) form an ∼23 km long arc–backarc transect located ∼250 km north of New Zealand above the subducting Hikurangi Plateau. Although only a short distance apart, rocks from the two volcanoes have different mineral and whole-rock geochemical compositions. Lavas from Rumble II East are predominantly basaltic and contain primitive olivine phenocrysts (≤Fo91), high-Mg# clinopyroxene (≤96) and anorthitic plagioclase (≤An97). Geochemically these lavas are very diverse and cover a spectrum from low Th/Yb ( 1014) to higher Th/Yb (>0·15) at lower Ba/Th (

High-resolution magnetic surveys have been acquired over the partially sedimented Palinuro massive sulfide deposits in the Aeolian volcanic arc, Tyrrhenian Sea. Surveys flown close to the seafloor using an autonomous underwater vehicle (AUV) show that the volcanic-arc-related basalt-hosted hydrothermal site is associated with zones of lower magnetization. This observation reflects the alteration of basalt affected by hydrothermal circulation and/or the progressive accumulation of a nonmagnetic deposit made of hydrothermal and volcaniclastic material and/or a thermal demagnetization of titanomagnetite due to the upwelling of hot fluids. To discriminate among these inferences, estimate the shape of the nonmagnetic deposit and the characteristics of the underlying altered area—the stockwork—we use high-resolution vector magnetic data acquired by the AUV Abyss (GEOMAR) above a crater-shaped depression hosting a weakly active hydrothermal site. Our study unveils a relatively small nonmagnetic deposit accumulated at the bottom of the depression and locked between the surrounding volcanic cones. Thermal demagnetization is unlikely but the stockwork extends beyond the limits of the nonmagnetic deposit, forming lobe-shaped zones believed to be a consequence of older volcanic episodes having contributed in generating the cones.

A one arc-minute resolution gravimetric quasigeoid model has been computed for New Zealand, covering the region $$25^{\circ }\hbox {S}$$ – $$60^{\circ }\hbox {S}$$ and $$160^{\circ }\hbox {E}$$ – $$170^{\circ }\hbox {W}$$ . It was calculated by Wong and Gore modified Stokes integration using the remove–compute–restore technique with the EIGEN-6C4 global gravity model as the reference field. The gridded gravity data used for the computation consisted of 40,677 land gravity observations, satellite altimetry-derived marine gravity anomalies, historical shipborne marine gravity observations and, importantly, approximately one million new airborne gravity observations. The airborne data were collected with the specific intention of reinforcing the shortcomings of the existing data in areas of rough topography inaccessible to land gravimetry and in coastal areas where shipborne gravimetry cannot be collected and altimeter-derived gravity anomalies are generally poor. The new quasigeoid has a nominal precision of $$\pm \,48\,\hbox {mm}$$ on comparison with GPS-levelling data, which is approximately $$14\,\hbox {mm}$$ less than its predecessor NZGeoid09.

Sea-floor imagery, volcanic rock, massive sulfide, and hydrothermal plume samples (δ 3 He, pH, dissolved Fe and Mn, and particulate chemistry) have been collected from the Rumble II West volcano, southern Kermadec arc, New Zealand. Rumble II West is a caldera volcano with an ∼3-km-diameter summit depression bounded by ring faults with a resurgent central cone. Rocks recovered to date are predominantly mafic in composition (i.e., basalt to basaltic andesite) with volumetrically lesser intermediate rocks (i.e., andesite). On the basis of its size, geometry, volcanic products, and composition, Rumble II West can be classified as a mafic caldera volcano. Rumble II West has a weak hydrothermal plume signature characterized by a small but detectable δ 3 He anomaly (25%). Time-series light scattering data though, obtained from vertical casts and tow-yos, do show that hydrothermal activity has increased in intensity between 1999 and 2011. Massive sulfides recovered from the eastern caldera wall and eastern flank of the central cone are primarily comprised of barite and chalcopyrite, with lesser sphalerite, pyrite, and traces of galena. The weak hydrothermal plume signal indicates that the volcano is in a volcanic-hydrothermal quiescent stage compared to other volcanoes along the southern Kermadec arc, although the preponderance of barite with massive sulfide mineralization indicates higher temperature venting in the past. Of the volcanoes along the Kermadec-Tonga arc known to host massive sulfides (i.e., Clark, Rumble II West, Brothers, Monowai, Volcano 19, and Volcano 1), the majority (five out of six) are dominantly mafic in composition and all but one of these mafic volcanoes form moderate-size to large calderas. To date, mafic calderas have been largely ignored as hosts to sea-floor massive sulfide deposits. That 75% of the presently known massive sulfide-bearing calderas along the arc are mafic in composition (the dacitic Brothers volcano is the exception) has important implications for sea-floor massive sulfide mineral exploration in the modern oceans and ancient rock record on land.

SUMMARY New high-resolution bathymetric and magnetic data from the western Aeolian sector, southern Tyrrhenian Sea, provide insights into structural and volcanic development of the area, suggesting a strong interaction between volcanism and tectonics. The analysis of these data combined with relocated earthquake distribution, focal plane solutions and strain rate evaluation indicates that the dextral strike-slip Sisifo-Alicudi shear zone is a complex and wide area of active deformation, representing the superficial expression of the deep seated lithospheric tear fault separating the subduction slab below Sicily and Calabria. Most of the observed volcanic features are aligned along a NW–SE trend, such as the Filicudi island-Alicudi North Seamount and Eolo-Enarete alignments, and are dissected by hundred-metre-high scarps along conjugate NNE–SSW trending fault systems. The magnetic field pattern matches the main trends of volcanic features. Spectral analysis and Euler deconvolution of magnetic anomalies show the existence of both deep and shallow sources. High-amplitude, high-frequency anomalies due to shallow sources are dominant close to the volcanic edifices of Alicudi and Filicudi, while the main contribution on the surrounding Eolo, Enarete, Alicudi North and Filicudi North seamounts is given by low-amplitude anomalies and/or deeper magnetic sources. This is probably related to different ages of the volcanic rocks, although hydrothermal processes may have played an important role in blanketing magnetic anomalies, in particular at Enarete and Eolo seamounts. Relative chronology of the eruptive centres and the inferred deformation pattern outline the Quaternary evolution of the western Aeolian Arc: Sisifo, Alicudi North and Filicudi North seamounts might have developed in an early stage, following the Late Pliocene–Early Pleistocene SE-ward migration of arc-related volcanism due to the Ionian subduction hinge retreat; Eolo, Enarete and Filicudi represent later manifestations that led volcanoes to develop during Mid-Late Pleistocene, when the stress regime in the area changed, due to the SSE-ward propagation of the subduction slab tear fault and the consequent reorientation and decrease of trench migration velocity. Finally, volcanic activity occurred in a very short time span at Alicudi, where an almost conical volcanic edifice emerged, suggesting negligible interactions with regional fault systems.

I thank Pedersen for raising some intriguing points that deserve a detailed discussion as well as finding a mistake due to a residual minus sign. Below I correct the spectrum that replaces equation 19 of Tontini (2005) using the same convention as Pedersen for the Fourier operator.

SUMMARY The final aim of a magnetic survey is to obtain quantitative information about the geological source generating the observed field. In this paper, I show how the employment of a suitable non-iterative operator in the Fourier domain permits to recover the horizontal magnetization distribution of a plausible equivalent source by a simple deconvolution. The results obtained in the synthetic case show a good agreement with the generating model, and moreover the recovered magnetization map appears more stable concerning high-frequency noise amplification with respect to other operators previously used in literature. The analytic justification of this stability will be described in detail by comparison with the traditional operators. The extension of this operator to the sector of field transformations, such as Upward Continuation or Reduction to the Pole, is straightforward because it is sufficient to re-calculate the field due to the apparent magnetization distribution. The results obtained by synthetic and real data tests demonstrate that the equivalent-source operator shown in this paper is able to give meaningful answers with a minimum amount of information by the user concerning depth of the source. These parameters as known, when unavailable, can be estimated in a stable way from the power spectrum of the magnetic anomaly, so that in the final analysis only a Fourier transformation of the data is needed to evaluate a realistic magnetization model.

I obtain a closed-form expression for the total-field magnetic anomaly generated by a 3D Gaussian magnetization distribution. The solution is evaluated by a novel contour integration in the complex-wavenumber domain. The spatial-domain solution is then computed by a numerical inverse Fourier transform. This calculated field can be used to interpret magnetic data because a 3D Gaussian function can represent a compact and smooth source. The method, moreover, allows for calculation of the Fourier components as a function of the distance between the depth of the source and the height of the observation level. The solution is thus useful for the inverse problem, spectral analysis, or data continuation.

[1] We obtained areal variations of crustal thickness, magnetic intensity, and degree of melting of the sub-axial upwelling mantle at Thetis and Nereus Deeps, the two northernmost axial segments of initial oceanic crustal accretion in the Red Sea, where Arabia is separating from Africa. The initial emplacement of oceanic crust occurred at South Thetis and Central Nereus roughly ∼2.2 and ∼2 Ma, respectively, and is taking place today in the northern Thetis and southern Nereus tips. Basaltic glasses major and trace element composition suggests a rift-to-drift transition marked by magmatic activity with typical MORB signature, with no contamination by continental lithosphere, but with slight differences in mantle source composition and/or potential temperature between Thetis and Nereus. Eruption rate, spreading rate, magnetic intensity, crustal thickness and degree of mantle melting were highest at both Thetis and Nereus in the very initial phases of oceanic crust accretion, immediately after continental breakup, probably due to fast mantle upwelling enhanced by an initially strong horizontal thermal gradient. This is consistent with a rift model where the lower continental lithosphere has been replaced by upwelling asthenosphere before continental rupturing, implying depth-dependent extension due to decoupling between the upper and lower lithosphere with mantle-lithosphere-necking breakup before crustal-necking breakup. Independent along-axis centers of upwelling form at the rifting stage just before oceanic crust accretion, with buoyancy-driven convection within a hot, low viscosity asthenosphere. Each initial axial cell taps a different asthenospheric source and serves as nucleus for axial propagation of oceanic accretion, resulting in linear segments of spreading.

The 500 m.y. cycle whereby continents assemble in a single supercontinent and then fragment and disperse again involves the rupturing of a continent and the birth of a new ocean, with the formation of passive plate margins. This process is well displayed today in the Red Sea, where Arabia is separating from Africa. We carried out geophysical surveys and bottom rock sampling in the two Red Sea northernmost axial segments of initial oceanic crust accretion, Thetis and Nereus. Areal variations of crustal thickness, magnetic intensity, and degree of melting of the subaxial upwelling mantle reveal an initial burst of active oceanic crust generation and rapid seafloor spreading below each cell, occurring as soon as the lid of continental lithosphere breaks. This initial pulse may be caused by edge-driven subrift mantle convection, triggered by a strong horizontal thermal gradient between the cold continental lithosphere and the hot ascending asthenosphere. The thermal gradient weakens as the oceanic rift widens; therefore the initial active pulse fades into steady, more passive crustal accretion, with slower spreading and along axis rift propagation.

We describe an inversion method for 3D geomagnetic data based on approximation of the source distribution by means of positive constrained Gaussian functions. In this way, smoothness and positivity are automatically imposed on the source without any subjective input from the user apart from selecting the number of functions to use. The algorithm has been tested with synthetic data in order to resolve sources at very different depths, using data from one measurement plane only.The forward modeling is based on prismatic cell parameterization, but the algebraic nonuniqueness is reduced because a relationship among the cells, expressed by the Gaussian envelope, is assumed to describe the spatial variation of the source distribution.We assume that there is no remanent magnetization and that the magnetic data are produced by induced magnetization only, neglecting any demagnetization effects. The algorithm proceeds by minimization of a χ2misfit function between real and predicted data using a nonlinear Levenberg‐Marquardt iteration scheme, easily implemented on a desktop PC, without any additional regularization. We demonstrate the robustness and utility of the method using synthetic data corrupted by pseudorandom generated noise and a real field data set.