Search

Subsurface fluid escape structures are geological features which are commonly observed in sedimentary basins worldwide. Their identification and description have implications for various subsurface fluid flow applications, such as assuring integrity of overburden rocks to geological CO2 storage sites. In this study, we applied 3-D first-arrival traveltime tomography to a densely sampled wide-azimuth and wide-angle ocean bottom seismometer (OBS) data set collected over the Scanner Pockmark complex, a site of active gas venting in the North Sea. Seismic reflection data show a chimney structure underlying the Scanner Pockmark. The objective of this study was to characterize this chimney as a representative fluid escape structure in the North Sea. An area of 6× 6 km2 down to a depth of 2 km below sea level was investigated using a regularized tomography algorithm. In total, 182 069 manually picked traveltimes from 24 OBS were used. Our final velocity model contains compressional wave velocity perturbations ranging from −125 to +110 ms−1 relative to its average 1-D model and compares favourably with a coincident seismic reflection data set. The tomographic velocity model reveals that the chimney as observed in seismic reflection data is part of a larger complex fluid escape structure, and discriminates the genuine chimney from seismic artefacts. We find that part of the seeping gas migrates from a deep source, accumulates beneath the Crenulate Reflector unconformity at ∼250 m below seafloor (mbsf) before reaching the porous sediments of the Ling Bank and Coal Pit formation at <100 mbsf. In addition, the model shows that the venting gas at Scanner Pockmark is also being fed laterally through a narrow NW–SE shallow channel. Quantitative velocity analysis suggests a patchy gas saturation within the gas-charged sediments of the Ling Bank and the Coal Pit formations. Confined to the well-resolved regions, we estimate a base case average gas saturation of ∼9 per cent and in-si

Serpentinized mantle peridotites form prominent mountains, including the highest elevations of the Troodos ophiolite in Cyprus (Mount Olympus, 1,952 m), but to date, only qualitative mechanisms have been proposed to explain the uplift of mantle rocks to high altitudes. Serpentinization reactions between mantle rocks and water result in profound changes to the rheology and physical properties of peridotites including significant density reduction (∼900 kg/m3). Field observations, density measurements, and isostatic uplift and erosional modeling provide new constraints on the contribution of serpentinization to the uplift of the Troodos Mantle Sequence. Different serpentinization styles have resulted in two distinct serpentinite domains with contrasting densities. Our modeling shows that the Troodos Mountains can form within the geologically constrained uplift time frame (∼5.5 Myr) exclusively through partial serpentinization reactions. We interpret the serpentinite domains as two nested diapirs that formed due to different extents of serpentinization and density reduction. Differential uplift and exhumation have decoupled the two serpentinite diapirs from the originally overlying ocean crustal rocks. Once at high altitudes the incursion of meteoric water reinforced coupled deformation-alteration-recrystallization processes in the shallow subsurface producing a localized low density completely serpentinized diapir. A second decoupling between the contrasting serpentinite diapirs results in localized differential uplift and exhumation, extruding deep materials to the east of Mount Olympus. Application of our modeling to other serpentinite massifs (e.g., St. Peter and St. Paul Rocks, New Idria, California) highlights the contribution of isostasy to the uplift of serpentinized massifs.

The equatorial eastern Indian Ocean hosts a diffuse plate boundary, where widespread deformation accommodates the relative motion between the Indian, Australian and Capricorn sub-plates. We integrate IODP Expedition 362 borehole data, which for the first time provides an accurate, ground-truthed chronostratigraphy of the sedimentary sequence east of the Ninety East Ridge (NER), with 2D seismic reflection profiles and multibeam bathymetry to assess the styles of faulting between the NER and the Sunda subduction zone, timing of activity and comparison with physical and rheological properties. We identify four distinct fault sets east of the NER in the northern Wharton Basin. N-S (350-010°) orientated faults, associated with the N-S fracture zones formed at the now extinct Wharton spreading centre, are still active and have been continuously active since at least 10 Ma. NNE- and WNW-trending fault fabrics develop between the fracture zones. The orientations and likely sense of displacement on these three sets of faults defines a Riedel shear system responding to ∼NNE-SSW left-lateral strike-slip activity at depth, demonstrated by the recent 2012 great intraplate earthquakes. We also find evidence of ∼NE-SW reverse faults, similar in style to E-W reverse faults observed west of the NER, where reverse faulting is more dominant. We show that the activity of this strike-slip system increased ca. 7-9 Ma, contemporaneous with reverse faulting and intraplate deformation west of the NER.

Drill sites in the southern Bay of Bengal at 3 degrees N 91 degrees E (International Ocean Discovery Program Expedition 362) have sampled for the first time a complete section of the Nicobar Fan and below to the oceanic crust. This generally overlooked part of the Bengal-Nicobar Fan System may provide new insights into uplift and denudation rates of the Himalayas and Tibetan Plateau. The Nicobar Fan comprises sediment gravity-flow deposits, mostly turbidites, that alternate with hemipelagite drapes and pelagite intervals of varying thicknesses. The decimetre-thick to metre-thick oldest pre-fan sediments (limestones/chalks) dated at 69 Ma are overlain by volcanic material and slowly accumulated pelagites (0.5 g cm(-2) kyr(-1)). At Expedition 362 Site U1480, terrigenous input began in the early Miocene at ca 22.5 Ma as muds, overlain by very thin-bedded and thin-bedded muddy turbidites at ca 19.5 Ma. From 9.5 Ma, sand content and sediment supply sharply increase (from 1-5 to 10-50 g cm(-2) kyr(-1)). Despite the abundant normal faulting in the Nicobar Fan compared with the Bengal Fan, it offers a better-preserved and more homogeneous sedimentary record with fewer unconformities. The persistent connection between the two fans ceased at 0.28 Ma when the Nicobar Fan became inactive. The Nicobar Fan is a major sink for Himalaya-derived material. This study presents integrated results of International Ocean Discovery Program Expedition 362 with older Deep Sea Drilling Project/Ocean Drilling Program/International Ocean Discovery Program sites that show that the Bengal-Nicobar Fan System experienced successive large-scale avulsion processes that switched sediment supply between the Bengal Fan (middle Miocene and late Pleistocene) and the Nicobar Fan (late Miocene to early Pleistocene). A quantitative analysis of the submarine channels of the Nicobar Fan is also presented, including their stratigraphic frequency, showing that channel size/area and abundance peaked at ca 2 to 3

Author Posting. © American Geophysical Union, 2017. This article is posted here by permission of [publisher] for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 122 (2017): 8675–8708, doi:10.1002/2017JB014225., In the North Atlantic Ocean, the geometry of diachronous V-shaped features that straddle the Reykjanes Ridge is often attributed to thermal pulses which advect away from the center of the Iceland plume. Recently, two alternative hypotheses have been proposed: rift propagation and buoyant mantle upwelling. Here we evaluate these different proposals using basin-wide geophysical and geochemical observations. The centerpiece of our analysis is a pair of seismic reflection profiles oriented parallel to flow lines that span the North Atlantic Ocean. V-shaped ridges and troughs are mapped on both Neogene and Paleogene oceanic crust, enabling a detailed chronology of activity to be established for the last 50 million years. Estimates of the cumulative horizontal displacement across normal faults help to discriminate between brittle and magmatic modes of plate separation, suggesting that crustal architecture is sensitive to the changing planform of the plume. Water-loaded residual depth measurements are used to estimate crustal thickness and to infer mantle potential temperature which varies by ±25°C on timescales of 3–8 Ma. This variation is consistent with the range of temperatures inferred from geochemical modeling of dredged basaltic rocks along the ridge axis itself, from changes in Neogene deep-water circulation, and from the regional record of episodic Cenozoic magmatism. We conclude that radial propagation of transient thermal anomalies within an asthenospheric channel that is 150 ± 50 km thick best accounts for the available geophysical and geochemical observations., NERC Grant Number: NE/G007632/1; Girdler Fund, University of Cambridge; BP Exploration, 2018-05-14

Characterization of the top 10–50 m of the subseabed is key for landslide hazard assessment, offshore structure engineering design and underground gas-storage monitoring. In this paper, we present a methodology for the stochastic inversion of ultra-high-frequency (UHF, 0.2–4.0 kHz) pre-stack seismic reflection waveforms, designed to obtain a decimetric-resolution remote elastic characterization of the shallow sediments with minimal pre-processing and little a priori information. We use a genetic algorithm in which the space of possible solutions is sampled by explicitly decoupling the short and long wavelengths of the P-wave velocity model. This approach, combined with an objective function robust to cycle skipping, outperforms a conventional model parametrization when the ground-truth is offset from the centre of the search domain. The robust P-wave velocity model is used to precondition the width of the search range of the multiparameter elastic inversion, thereby improving the efficiency in high-dimensional parametrizations. Multiple independent runs provide a set of independent results from which the reproducibility of the solution can be estimated. In a real data set acquired in Finneidfjord, Norway, we also demonstrate the sensitivity of UHF seismic inversion to shallow subseabed anomalies that play a role in submarine slope stability. Thus, the methodology has the potential to become an important practical tool for marine ground model building in spatially heterogeneous areas, reducing the reliance on expensive and time-consuming coring campaigns for geohazard mitigation in marine areas.

A holistic view of the Bengal-Nicobar Fan system requires sampling the full sedimentary section of the Nicobar Fan, which was achieved for the first time by International Ocean Discovery Program (IODP) Expedition 362 west of North Sumatra. We identified a distinct rise in sediment accumulation rate (SAR) beginning similar to 9.5 Ma and reaching 250-350 m/Myr in the 9.5-2 Ma interval, which equal or far exceed rates on the Bengal Fan at similar latitudes. This marked rise in SAR and a constant Himalayan-derived provenance necessitates a major restructuring of sediment routing in the Bengal-Nicobar submarine fan. This coincides with the inversion of the Eastern Himalayan Shillong Plateau and encroachment of the west-propagating Indo-Burmese wedge, which reduced continental accommodation space and increased sediment supply directly to the fan. Our results challenge a commonly held view that changes in sediment flux seen in the Bengal-Nicobar submarine fan were caused by discrete tectonic or climatic events acting on the Himalayan-Tibetan Plateau. Instead, an interplay of tectonic and climatic processes caused the fan system to develop by punctuated changes rather than gradual progradation.

In the North Atlantic Ocean, the geometry of diachronous V-shaped features that straddle the Reykjanes Ridge is often attributed to thermal pulses which advect away from the center of the Iceland plume. Recently, two alternative hypotheses have been proposed: rift propagation and buoyant mantle upwelling. Here, we evaluate these different proposals using basin-wide geophysical and geochemical observations. The centerpiece of our analysis is a pair of seismic reflection profiles oriented parallel to flowlines that span the North Atlantic Ocean. V-shaped ridges and troughs are mapped on both Neogene and Paleogene oceanic crust, enabling a detailed chronology of activity to be established for the last 50 million years. Estimates of the cumulative horizontal displacement across normal faults help to discriminate between brittle and magmatic modes of plate separation, suggesting that crustal architecture is sensitive to the changing planform of the plume. Water-loaded residual depth measurements are used to estimate crustal thickness and to infer mantle potential temperature which varies by 25◦C on timescales of 3–8 Ma. This variation is consistent with the range of temperatures inferred from geochemical modeling of dredged basaltic rocks along the ridge axis itself, from changes in Neogene deep-water circulation, and from the regional record of episodic Cenozoic magmatism. We conclude that radial propagation of transient thermal anomalies within an asthenospheric channel that is 150 50 km thick best accounts for the available geophysical and geochemical observations.

The full waveform inversion (FWI) of seismic reflection data aims to reconstruct a detailed physical properties model of the subsurface, fitting both the amplitude and the traveltime of the reflections generated at physical discontinuities in the propagation medium. Unlike reservoir-scale seismic exploration, where seismic inversion is a widely adopted remote characterization tool, ultrahigh-frequency (UHF, 0.2–4.0 kHz) multichannel marine reflection seismology is still most often limited to a qualitative interpretation of the reflections’ architecture. Here we propose an elastic FWI methodology, custom-tailored for pre-stack UHF marine data in vertically heterogeneous media to obtain a decimetric-scale distribution of P-impedance, density and Poisson’s ratio within the shallow subseabed sediments. We address the deterministic multiparameter inversion in a sequential fashion. The complex trace instantaneous phase is first inverted for the P-wave velocity to make up for the lack of low frequency in the data and reduce the nonlinearity of the problem. This is followed by a short-offset P-impedance optimization and a further step of full offset range Poisson’s ratio inversion. Provided that the seismogram contains wide reflection angles (>40°), we show that it is possible to invert for density and decompose a posteriori the relative contribution of P-wave velocity and density to the P-impedance. A broad range of synthetic tests is used to prove the potential of the methodology and highlights sensitivity issues specific to UHF seismic. An example application to real data is also presented. In the real case, trace normalization is applied to minimize the systematic error deriving from an inaccurate source wavelet estimation. The inverted model for the top 15 m of the subseabed agrees with the local lithological information and core-log data. Thus, we can obtain a detailed remote characterization of the shallow sediments using a multichannel sub-bottom profiler within a r