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.

Calderas are major volcanic features with large volcanic and seismic hazards. They also host diverse microbiota, provide heat, energy, mineral and economic benefits. Despite their scientific and socio-economic importance, we still do not completely understand calderas and the interactions between volcanism, tectonism, fluid circulation and the deep biosphere because in-situ and subsurface observations are sparse.The Okataina Volcanic Centre (OVC) in Aotearoa New Zealand, is one of two active giant calderas of the Taupō Volcanic Zone within the rapidly extending continental intra-arc Taupō Rift. This superb natural laboratory has: 1) numerous past eruptions of varied size and style, 2) documented co-eruptive earthquakes, 3) vigorous hydrothermal manifestations, 4) diverse microbial communities in hot springs but unknown in the subsurface.We propose to establish a scientific drilling programme at the OVC to address:What are the conditions leading to volcanic eruptions; and volcano-tectonic feedbacks in intra-rift calderas? What controls fluid circulations in active calderas/rift regions? Does subsurface microbial community composition vary with tectonic and/or volcanic activity? High temperatures complicate drillhole design, restrict data collection and prevent exploration of the biosphere. By targeting the cooler parts of the caldera, this project will use conventional engineering to maximise sampling (drill cores and fluids), downhole logging and establish long-term observatories.Two preliminary drill targets are suggested: (1) in the centre of the caldera; (2) through the caldera margin. Drill data will provide a comprehensive record of past activity, establishing eruption frequency-magnitude relationships and precursors. Combined with well-known fault rupture history, the relative timing of tectonic and magmatic activity will be untangled. Drill data will unravel the relationships between the groundwater and hydrothermal systems, magma, faults and stress, informing thermo-hydro-mechanical regional caldera models with findings applicable worldwide. Drill cores and a dedicated fluid sampler triggered by nearby earthquakes will reveal the composition, function and potential change of microbial activity in response to rock and fluid variations.The programme is informed by indigenous Māori, regulatory authorities and emergency managers to ensure scientific, cultural, regulatory and resilience outcomes. The programme will underpin 1) community resilience to volcanic and seismic hazards; 2) sustainable management of groundwater and geothermal resources, and 3) understanding of subsurface microbial diversity, function and geobiological interactions. At these early stages of planning, we invite the scientific community to contribute to the concept of this project in the exceptional OVC settings and strengthen linkages with other ongoing research and scientific drilling programmes.

Significance Much of Earth’s volcanism occurs in the deep sea, yet little is known about the microbial communities inhabiting such extreme and dynamic systems. Using a multidisciplinary approach to study distinct hydrothermal systems at Brothers submarine arc volcano, we provide insights into how microbial community composition and function reflect subtly different fluid chemistries resulting from subsurface fluid interactions with distinct alteration mineral assemblages. These variations can be traced to the subsurface hydrogeologic history beneath Brothers volcano. Further, we show that these systems represent oases of phylogenetically diverse Archaea and Bacteria. Our findings highlight the importance of geologic legacy in understanding drivers of microbial diversity, assembly, and evolution and may have insights into processes that drove early diversification of life on Earth., Hydrothermally active submarine volcanoes are mineral-rich biological oases contributing significantly to chemical fluxes in the deep sea, yet little is known about the microbial communities inhabiting these systems. Here we investigate the diversity of microbial life in hydrothermal deposits and their metagenomics-inferred physiology in light of the geological history and resulting hydrothermal fluid paths in the subsurface of Brothers submarine volcano north of New Zealand on the southern Kermadec arc. From metagenome-assembled genomes we identified over 90 putative bacterial and archaeal genomic families and nearly 300 previously unknown genera, many potentially endemic to this submarine volcanic environment. While magmatically influenced hydrothermal systems on the volcanic resurgent cones of Brothers volcano harbor communities of thermoacidophiles and diverse members of the superphylum “DPANN,” two distinct communities are associated with the caldera wall, likely shaped by two different types of hydrothermal circulation. The communities whose phylogenetic diversity primarily aligns with that of the cone sites and magmatically influenced hydrothermal systems elsewhere are characterized predominately by anaerobic metabolisms. These populations are probably maintained by fluids with greater magmatic inputs that have interacted with different (deeper) previously altered mineral assemblages. However, proximal (a few meters distant) communities with gene-inferred aerobic, microaerophilic, and anaerobic metabolisms are likely supported by shallower seawater-dominated circulation. Furthermore, mixing of fluids from these two distinct hydrothermal circulation systems may have an underlying imprint on the high microbial phylogenomic diversity. Collectively our results highlight the importance of considering geologic evolution and history of subsurface processes in studying microbial colonization and community dynamics in volcanic environments.

In 2018, Bunn and Nolden published a paper purporting to have established the true position of the Pink and White Terraces of Lake Rotomahana using forensic cartography, ‘reverse engineering’ an 18...

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.

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 corrections

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.

Hydrothermal systems hosted by submarine arc volcanoes commonly include a large component of magmatic fluid. The high Cu-Au contents and strongly acidic fluids in these systems are similar to those that formed in the shallow parts of some porphyry copper and epithermal gold deposits mined today on land. Two main types of hydrothermal systems occur along the submarine portion of the Kermadec arc (offshore New Zealand): magmatically influenced and seawater-dominated systems. Brothers volcano hosts both types. Here, we report results from a series of drill holes cored by the International Ocean Discovery Program into these two types of hydrothermal systems. We show that the extent of hydrothermal alteration of the host dacitic volcaniclastics and lavas reflects primary lithological porosity and contrasting spatial and temporal contributions of magmatic fluid, hydrothermal fluid, and seawater. We present a two-step model that links the changes in hydrothermal fluid regime to the evolution of the volcano caldera. Initial hydrothermal activity, prior to caldera formation, was dominated by magmatic gases and hypersaline brines. The former mixed with seawater as they ascended toward the seafloor, and the latter remained sequestered in the subsurface. Following caldera collapse, seawater infiltrated the volcano through fault-controlled permeability, interacted with wall rock and the segregated brines, and transported associated metals toward the seafloor and formed Cu-Zn-Au–rich chimneys on the caldera walls and rim, a process continuing to the present day. This two-step process may be common in submarine arc caldera volcanoes that host volcanogenic massive sulfide deposits, and it is particularly efficient at focusing mineralization at, or near, the seafloor.

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 (

Back-arc basins are found at convergent plate boundaries. Nevertheless, they are zones of significant crustal extension that show volcanic and hydrothermal processes somewhat similar to those of mid-ocean ridges. Accepted models imply the initial rifting and thinning of a pre-existing volcanic arc until seafloor spreading gradually develops over timescales of a few million years. The Havre Trough northeast of New Zealand is a unique place on Earth where the early stages of back-arc basin formation are well displayed in the recent geological record. Here we present evidence that, in this region, rifting of the original volcanic arc occurred in a very narrow area about 10–15 km wide, which could only accommodate minimal stretching for a very short time before mass balance required oceanic crustal accretion. An initial burst of seafloor spreading started around 5.5–5.0 million years ago and concluded abruptly about 3.0–2.5 million years ago, after which arc magmatism dominated the crustal accretion. The sudden transition between these different tectonomagmatic regimes is linked to trench rollback promoted by gradual sinking of the subducting lithosphere, which could have diverted the arc flux outside the region of seafloor spreading and induced the vertical realignment of surface volcanism with the source of arc melts at depth. The Havre Tough back-arc basin, New Zealand, formed rapidly and in two phases: initial, limited seafloor spreading was followed by a transition to arc magmatism, as shown by geophysical data and modelling.

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 submersible study of the products of a large submarine eruption demonstrates the influence of the ocean on eruption dynamics., The 2012 submarine eruption of Havre volcano in the Kermadec arc, New Zealand, is the largest deep-ocean eruption in history and one of very few recorded submarine eruptions involving rhyolite magma. It was recognized from a gigantic 400-km2 pumice raft seen in satellite imagery, but the complexity of this event was concealed beneath the sea surface. Mapping, observations, and sampling by submersibles have provided an exceptionally high fidelity record of the seafloor products, which included lava sourced from 14 vents at water depths of 900 to 1220 m, and fragmental deposits including giant pumice clasts up to 9 m in diameter. Most (>75%) of the total erupted volume was partitioned into the pumice raft and transported far from the volcano. The geological record on submarine volcanic edifices in volcanic arcs does not faithfully archive eruption size or magma production.

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.

Large igneous province subduction is a rare process on Earth. A modern example is the subduction of the oceanic Hikurangi Plateau beneath the southern Kermadec arc, offshore New Zealand. This segment of the arc has the largest total lava volume erupted and the highest volcano density of the entire Kermadec arc. Here we show that Kermadec arc lavas south of B32°S have elevated Pb and Sr and low Nd isotope ratios, which argues, together with increasing seafloor depth, forearc retreat and crustal thinning, for initial Hikurangi Plateau—Kermadec arc collision B250 km north of its present position. The combined data set indicates that a much larger portion of the Hikurangi Plateau (the missing Ontong Java Nui piece) than previously believed has already been subducted. Oblique plate convergence caused southward migration of the thickened and buoyant oceanic plateau crust, creating a buoyant ‘Hikurangi’ me´lange beneath the Moho that interacts with ascending arc melts.

[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.