Description: Highlight • It is important to develop systems able to detect and classify mineralized zones from waste materials while drilling deep-water; • Seismic P-wave velocities (Vp) were measured on 40 SMS and unmineralized mini-cores; • The porosity was back-calculated from Vp; • The results were compared with electrical resistivity measurements; • Using Archie’s Law, it is possible to observe that metallic conduction exists. Abstract Deep-sea mining exploration for seafloor massive sulfide (SMS) deposits is currently increasing. At present, most exploration activities are surficial and use indirect methods to identify potential sites and perform first assessments. For a proper resource estimate, however, drilling is inevitable. By using seabed drill rigs, exploration costs can be reduced considerably. SMS deposits are normally found at depths between 1000 and 4000 m and in order for deep sea mining to be implemented, reliable technologies are needed. Additionally, the development of geophysical systems that can detect and classify mineralized zones from waste materials while drilling could decrease costs and speed up offshore operations by limiting the amount of drilling of unmineralized materials. This paper shows how the physical properties of SMS can be used to discriminate between host rocks and mineralization. Seismic P-wave velocities (Vp) were measured on 40 SMS and unmineralized mini-cores. By back-calculating the porosity from Vp, comparing the results with electrical resistivity measurements, and using Archie's Law, it is possible to observe that metallic conduction exists. For deep-sea mineral exploration, the combination of seismic tests, electrical resistivity and magnetic susceptibility could support the preliminary discrimination of mineralized samples in the cores while drilling at the seafloor.

Description: One of the main problem of 3D time domain controlled source electromagnetic (CSEM) inversion is the high runtimes of forward modeling codes. We reduced the runtime of the 3D time domain finite difference CSEM code TEMDDD by the GPU-parallelization of expensive algorithms. The code solves the electromagnetic diffusion equation by discretization of spatial operators and subsequent calculation of eigenpairs. These eigenpairs are found by approximation of the eigenspace in a Kiylov subspace using the spectral Lanczos decomposition Method. This algorithm was in its original form not parallelizable due to implementation of the upper boundary condition at the air-water interface. We show for the marine case that replacing the original boundary condition at the air-water surface by a discretized air layer allows GPU parallelization of every time consuming algorithm of the code in the marine case. Speedups between 20 and 60 have been achieved compared to the original code for a larger 3-D model. In this model the bathymetry from a survey area offshore Egypt is used as an example demonstrating that the parallelized version of the code is applicable to real survey scenarios.

Description: Several known gas seep sites along the Hikurangi Margin off the east coast of New Zealand were surveyed by marine controlled source electromagnetic (CSEM) experiments. A bottom-towed electric dipole–dipole system was used to reveal the occurrence of gas hydrate and methane related to the seeps. The experiments were part of the international multidisciplinary research program “New Vents” carried out on German R/V Sonne in 2007 (cruise SO191) to study key parameters controlling the release and transformation of methane from marine cold vents and shallow gas hydrate deposits. Two CSEM lines have been surveyed over known seep sites on Opouawe Bank in the Wairarapa region off the SE corner of the North Island. The data have been inverted to sub-seafloor apparent resistivity profiles and one-dimensional layered models. Clearly anomalous resistivities are coincident with the location of two gas seep sites, North Tower and South Tower on Opouawe Bank. A layer of concentrated gas hydrate within the uppermost 100 m below the seafloor is likely to cause the anomalous resistivities, but free gas and thick carbonate crusts may also play a role. Seismic data show evidence of fault related venting which may also indicate the distribution of gas hydrates and/or authigenic carbonate. Geochemical profiles indicate an increase of methane flux and the formation of gas hydrate in the shallow sediment section around the seep sites. Takahe is another seep site in the area where active venting, higher heat flow, shallow gas hydrate recovered from cores, and seismic fault planes, but only moderately elevated resistivities have been observed. The reasons could be a) the gas hydrate concentration is too low, even though methane venting is evident, b) strong temporal or spatial variation of the seep activity, and c) the thermal anomaly indicates rather temperature driven fluid expulsion that hampers the formation of gas hydrate beneath the vent.

Description: Highlights • We report on marine 3D Magnetotelluric study on Walvis Ridge • Derived 3D electrical resistivity model shows a large scale resistive zone, which we link to crustal extension due to local uplift. It might indicate the location where the hot-spot impinged on the crust prior to rifting • Smaller scale resistive region is attributed to magma ascent during rifting • Rift basin is identified by low resistivity region The Namibian continental margin marks the starting point of the Tristan da Cunha hotspot trail, the Walvis Ridge. This section of the volcanic southwestern African margin is therefore ideal to study the interaction of hotspot volcanism and rifting, which occurred in the late Jurassic/early Cretaceous. Offshore magnetotelluric data image electromagnetically the landfall of Walvis Ridge. Two large-scale high resistivity anomalies in the 3-D resistivity model indicate old magmatic intrusions related to hot-spot volcanism and rifting. The large-scale resistivity anomalies correlate with seismically identified lower crustal high velocity anomalies attributed to magmatic underplating along 2-D offshore seismic profiles. One of the high resistivity anomalies (above 500 Ωm) has three arms of approximately 100 km width and 300 km to 400 km length at 120 degree angles in the lower crust. One of the arms stretches underneath Walvis Ridge. The shape is suggestive of crustal extension due to local uplift. It might indicate the location where the hot-spot impinged on the crust prior to rifting. A second, smaller anomaly of 50 km width underneath the continent ocean boundary may be attributed to magma ascent during rifting. We attribute a low resistivity anomaly east of the continent ocean boundary and south of Walvis Ridge to the presence of a rift basin that formed prior to the rifting.

Description: The active volcanic island Tristan da Cunha, located at the southwestern and youngest end of the Walvis Ridge - Tristan/Gough hotspot track, is believed to be the surface expression of a huge thermal mantle anomaly. While several criteria for the diagnosis of a classical hotspot track are met, the Tristan region also shows some peculiarities. Consequently it is vigorously debated if the active volcanism in this region is the expression of a deep mantle plume, or if it is caused by shallow plate tectonics and the interaction with the nearby Mid-Atlantic Ridge. Because of a lack of geophysical data in the study area, no model or assumption has been completely confirmed. We present the first amphibian P-wave finite-frequency travel time tomography of the Tristan da Cunha region, based on cross-correlated travel time residuals of teleseismic earthquakes recorded by 24 ocean-bottom seismometers. The data can be used to image a low velocity structure southwest of the island. The feature is cylindrical with a radius of ~ 100 km down to a depth of 250 km. We relate this structure to the origin of Tristan da Cunha and name it the Tristan conduit. Below 250 km the low velocity structure ramifies into narrow veins, each with a radius of ~ 50 km. Furthermore, we imaged a linkage between young seamounts southeast of Tristan da Cunha and the Tristan conduit.

Description: Highlights • The electrical structure beneath the Tristan da Cunha (TDC) hotspot was investigated. • Plume-like structure was not imaged beneath TDC by 3-D inversion analysis. • The plume may be small and/or weak or take place elsewhere outside of the study area. • Conductivity and bathymetry anomalies show a contrast across the TDC fracture zone. • Mantle temperature and melting process at ridge may cause the conductivity anomaly. Abstract The Tristan da Cunha (TDC) is a volcanic island located above a prominent hotspot in the Atlantic Ocean. Many geological and geochemical evidences support a deep origin of the mantle material feeding the hotspot. However, the existence of a plume has not been confirmed as an anomalous structure in the mantle resolved by geophysical data because of lack of the observations in the area. Marine magnetotelluric and seismological observations were conducted in 2012–2013 to examine the upper mantle structure adjacent to TDC. The electrical conductivity structure of the upper mantle beneath the area was investigated in this study. Three-dimensional inversion analysis depicted a high conductive layer at ~ 120 km depth but no distinct plume-like vertical structure. The conductive layer is mostly flat independently on seafloor age and bulges upward beneath the lithospheric segment where the TDC islands are located compared to younger segment south of the TDC Fracture Zone, while the bathymetry is rather deeper than prediction for the northern segment. The apparent inconsistency between the absence of vertical structure in this study and geochemical evidences on deep origin materials suggests that either the upwelling is too small and/or weak to be resolved by the current data set or that the upwelling takes place elsewhere outside of the study area. Other observations suggest that 1) the conductivity of the upper mantle can be explained by the fact that the mantle above the high conductivity layer is depleted in volatiles as the result of partial melting beneath the spreading ridge, 2) the potential temperature of the segments north of the TDC Fracture Zone is lower than that of the southern segment at least during the past ~ 30 Myr.

Description: Highlights • Petrophysical joint inversion of MT, seismic and gravity for sub-basalt imaging • Use of adaptive coupling contraints result in a robust convergence behavior. • Adaptive coupling allows deviations from used parameter relationships in models. • Scheme is applied successfully on synthetic and real data examples. Abstract Joint inversion strategies for geophysical data have become increasingly popular as they allow for the efficient combination of complementary information from different data sets. The algorithm used for the joint inversion needs to be flexible in its description of the subsurface so as to be able to handle the diverse nature of the data. Hence, joint inversion schemes are needed that 1) adequately balance data from the different methods, 2) have stable convergence behavior, 3) consider the different resolution power of the methods used and 4) link the parameter models in a way that they are suited for a wide range of applications. Here, we combine active source seismic P-wave tomography, gravity and magnetotelluric (MT) data in a petrophysical joint inversion that accounts for these issues. Data from the different methods are inverted separately but are linked through constraints accounting for parameter relationships. An advantage of performing the inversions separately is that no relative weighting between the data sets is required. To avoid perturbing the convergence behavior of the inversions by the coupling, the strengths of the constraints are readjusted at each iteration. The criterion we use to control the adaption of the coupling strengths is based on variations in the objective functions of the individual inversions from one to the next iteration. Adaption of the coupling strengths makes the joint inversion scheme also applicable to subsurface conditions, where assumed relationships are not valid everywhere, because the individual inversions decouple if it is not possible to reach adequately low data misfits for the made assumptions. The coupling constraints depend on the relative resolutions of the methods, which leads to an improved convergence behavior of the joint inversion. Another benefit of the proposed scheme is that structural information can easily be incorporated in the petrophysical joint inversion (no additional terms are added in the objective functions) by using mutually controlled structural weights for the smoothing constraints. We test our scheme using data generated from a synthetic 2-D sub-basalt model. We observe that the adaption of the coupling strengths makes the convergence of the inversions very robust (data misfits of all methods are close to the target misfits) and that final results are always close to the true models independent of the parameter choices. Finally, the scheme is applied on real data sets from the Faroe-Shetland Basin to image a basaltic sequence and underlying structures. The presence of a borehole and a 3-D reflection seismic survey in this region allows direct comparison and, hence, evaluate the quality of the joint inversion results. The results from joint inversion are more consistent with results from other studies than the ones from the corresponding individual inversions and the shape of the basaltic sequence is better resolved. However, due to the limited resolution of the individual methods used it was not possible to resolve structures underneath the basalt in detail, indicating that additional geophysical information (e.g. CSEM, reflection onsets) needs to be included.

Description: Highlights • Receiver functions from ocean-bottom seismometer stations reveal no significant crustal thickening in the surrounding of the Tristan da Cunha hot spot. • The mantle transition zone to the NW of Tristan da Cunha is thickened and cool. • The mantle transition zone is potentially thinned to the south/southwest of Tristan da Cunha. • A thickness of 60 to 75 km beneath Tristan da Cunha argues for a compositional control on the seismological lithosphere in the South Atlantic. Abstract The most prominent hotspot in the South Atlantic is Tristan da Cunha, which is widely considered to be underlain by a mantle plume. But the existence, location and size of this mantle plume have not been established due to the lack of regional geophysical observations. A passive seismic experiment using ocean bottom seismometers aims to investigate the lithosphere and upper mantle structure beneath the hotspot. Using the Ps receiver function method we calculate a thickness of 5 to 8 km for the oceanic crust at 17 ocean-bottom stations deployed around the islands. Within the errors of the method the thickness of the oceanic crust is very close to the global mean. The Tristan hotspot seems to have contributed little additional magmatic material or heat to the melting zone at the mid-oceanic ridge, which could be detected as thickened oceanic crust. Magmatic activity on the archipelago and surrounding seamounts seems to have only effected the crustal thickness locally. Furthermore, we imaged the mantle transition zone discontinuities by analysing receiver functions at the permanent seismological station TRIS and surrounding OBS stations. Our observations provide evidence for a thickened (cold) mantle transition zone west and northwest of the islands, which excludes the presence of a deep-reaching mantle plume. We have some indications of a thinned, hot mantle transition zone south of Tristan da Cunha inferred from sparse and noisy observations, which might indicate the location of a Tristan mantle plume at mid-mantle depths. Sp receiver functions image the base of lithosphere at about 60 to 75 km beneath the islands, which argues for a compositionally controlled seismological lithosphere-asthenosphere boundary beneath the study area.

Description: Porangahau Ridge, located offshore the Wairarapa on the Hikurangi Margin, is an active ocean-continent collision region in northeastern New Zealand coastal waters. Bottom simulating reflections (BSRs) in seismic data indicate the potential for significant gas hydrate deposits across this part of the margin. Beneath Porangahau Ridge a prominent high-amplitude reflection band has been observed to extend from a deep BSR towards the seafloor. Review of the seismic data suggest that this high-amplitude band is caused by local shoaling of the base of gas hydrate stability due to advective heat flow and it may constitute the location of elevated gas hydrate concentrations. During R/V Tangaroa cruise TAN0607 in 2006 heat flow probing for measurements of vertical fluid migration, sediment coring for methane concentrations, and additional seismic profiles were obtained across the ridge. In a subsequent 2007 expedition, on R/V Sonne cruise SO191, a controlled source electromagnetic (CSEM) experiment was conducted along the same seismic, geochemical, and heat flow transect to reveal the electrical resistivity distribution. CSEM data highlight a remarkable coincidence of anomalously high resistivity along the western, landward flank of the ridge which point to locally higher gas hydrate concentration above the high amplitude reflection band. Measured sediment temperature profiles, also along the western flank, consistently show non-linear and concave geothermal gradients typical of advective heat flow. Geochemical data reveal elevated methane concentrations in surface sediments concomitant with a rapid decline in sulfate concentrations indicating elevated methane flux and oxidation of methane in conjunction with sulfate reduction at the landward ridge base. Together, these data sets suggest that the western rim of Porangahau Ridge is a tectonically driven zone of rising fluids that transport methane and cause an upward inflection of the base of gas hydrate stability and the formation of locally enriched gas hydrate above the reflective zone.

Description: Highlights • Joint seismic, CSEM and geochemical estimate of gas and hydrate saturation. • Converted shear wave analyses for gas hydrate distribution. • Wide angle full waveform inversion. • Paleo-BSRs and gas hydrate system in disequilibrium. • Numerical framework for joint inversion of multiple geophysical methods. Abstract Gas hydrates are naturally-occurring solid compounds of gas and water within almost all sediment-rich continental margins. Due to the large amounts of methane stored in submarine gas hydrates, they might serve as future reservoirs for offshore marine gas production. Assessing the reservoir characteristics requires reliable estimates of both the gas and gas hydrate concentration, which can be best addressed using geophysical and geological investigations. Here, we demonstrate the power of joint interpretation of interdisciplinary geophysical techniques and geological laboratory experiments. Regional 2D multichannel seismic data provide the broad overview of a hydrate-bearing area. High-resolution 2D and 3D seismic reflection data provide detailed images of two working areas, the buried S1 channel-levee system at 1500 m water depth (well within the gas hydrate stability zone) and a slope failure location, located at 665 m water depth (top limit of the hydrate formation) next to the S2 channel. Detailed compressional and shear wave (Vs) velocity-depth models were derived from four component ocean-bottom seismic data, the latter from P- to S-conversion upon reflection. Due to their steep reflection angles, shear wave events result in less resolved Vs models. Nevertheless, in case of a change in elasticity of the sediment matrix due to gas hydrate cementation, shear wave events can be used as an indicator. As such, Vs can give insight into the nature of hydrate formation throughout the GHSZ. We present new developments in the application of common reflection surface, normal-incidence-point tomography and full waveform inversion techniques to enhance model resolution for the seismic data sets. 2D and 3D controlled-source electromagnetic measurements provide volume information of the resistivity-depth distribution models. Electrical resistivity of the sediment formation depends on its porosity and the resistivity of the pore fluid. Gas hydrate and free gas generally have much higher electrical resistivities than saline pore fluid, and can be assessed using empirical relationships if the porosity and pore fluid salinity are known. Calibration with logging data, laboratory experiments on hydrate- or ice-bearing sediments, and resulting velocity and resistivity values, guide the joint interpretation into more accurate saturation estimations. Beyond that, a joint inversion framework supporting forward calculation of specialized geophysical methods at distributed locations is under development. In this paper, we summarize these individual components of a multi-parameter study, and their joint application to investigate gas hydrate systems, their equilibrium conditions and preservation of bottom-simulating-reflectors. We analyze data from two working areas at different locations and depth levels along the slope of the Danube Fan, which are both characterized by multiple bottom simulating reflectors indicating the presence of gas hydrate. In the first working area we located two depth windows with indications for moderate 16%–24% gas hydrate formation, but no vertical gas migration. In the second working area we observed fluid migration pathways and active gas seepage, limiting gas hydrate formation to less than 10% at the BSR. Some discrepancies remain between seismic-based and electromagnetic-based models of gas and gas hydrate distribution and saturation estimates, indicating that further in-situ investigations are likely required to better understand the gas hydrate systems at our study areas and to calibrate the inversion processes, which will be required for a joint inversion framework as well.