Seismic imaging of mid-Pleistocene tunnel-valleys in the North Sea Basin—high resolution from low frequencies

Seismic reflection methods provide continuous access both to stratigraphy (vertical) and to subsurface morphology (horizontal), for which the scales of interest may differ by orders of magnitude. Seismic surveys of Quaternary successions have generally sought to optimise vertical resolution, through the use of higher source frequency content. Here, I show that low-frequency bandwidth is not necessarily a limiting factor for the seismic resolution of glacigenic morpho-sedimentary features. Observations are presented from a buried network of large mid-Pleistocene (Elsterian) tunnel-valleys in the southern North Sea Basin, across a 100×130 km study area with water depths less than 30 m. Low-frequency 2D and 3D seismic multi-channel data, acquired for deeper hydrocarbon exploration, are compared with previously available high-frequency single- and multi-channel profiles (5–15 km grid spacing). The low-frequency data contribute to a new understanding of the basal morphology and fill stratigraphy of the tunnel-valleys, in part due to higher data densities (≥1 km grid spacing), but also to improved imaging of reflectors at depth. The tunnel-valleys are seen to be overdeepened troughs, shallow (≤0.5 km) relative to their widths (≤6 km). The basal unconformity defines a series of arborescent elements, convergent to the south; erosional overlap by younger elements to the north has resulted in anastomosing patterns in places. The fill is dominated by axially downlapping clinoforms, descending to the north, onlapped and overlain by subhorizontal reflectors. Well data show that sand-dominated glaciofluvial sediments are overlain by glaciolacustrine to marine muds. Better definition of the clinoforms on low- versus high-frequency multi-channel data is suggested to reflect the coarse spatial scale of the backset glaciofluvial strata. The results support a simple interpretation of time-transgressive tunnel-valley formation by coeval glaciofluvial erosion and backfill beneath the outer tens of kilometres of the northward receding Elsterian ice sheet margin. Comparable submarginal interpretations have been proposed for drainage features (tunnel-valleys and eskers) of the last deglaciation of both northern Europe and North America using integrated geomorphologic and stratigraphic methods. Commercial 2D and 3D seismic data are widely available from exploration areas such as the North Sea and are argued to constitute an underexploited resource for Quaternary research.     

A predictive numerical model for potential mapping of the gas hydrate stability zone in the Gulf of Cadiz

This paper presents a computational model for mapping the regional 3D distribution in which seafloor gas hydrates would be stable, that is carried out in a Geographical Information System (GIS) environment. The construction of the model is comprised of three primary steps, namely: (1) the construction of surfaces for the various variables based on available 3D data (seafloor temperature, geothermal gradient and depth-pressure); (2) the calculation of the gas function equilibrium functions for the various hydrocarbon compositions reported from hydrate and sediment samples; and (3) the calculation of the thickness of the hydrate stability zone. The solution is based on a transcendental function, which is solved iteratively in a GIS environment.The model has been applied in the northernmost continental slope of the Gulf of Cadiz, an area where an abundant supply for hydrate formation, such as extensive hydrocarbon seeps, diapirs and fault structures, is combined with deep undercurrents and a complex seafloor morphology. In the Gulf of Cadiz, the model depicts the distribution of the base of the gas hydrate stability zone for both biogenic and thermogenic gas compositions, and explains the geometry and distribution of geological structures derived from gas venting in the Tasyo Field (Gulf of Cadiz) and the generation of BSR levels on the upper continental slope.     

Tectonically-driven mud volcanism since the late Pliocene on the Calabrian accretionary prism,central Mediterranean Sea

Mud volcanoes recently discovered on the offshore Calabrian Arc are investigated at two sites 60 km apart, in water depths of 1650--2300 m, using swath bathymetry, 2D&3D multichannel seismic and cores. The seabed and subsurface data provide information on their formation and functioning in relation to tectonic activity during the rapid Plio-Quaternary advance of the accretionary prism. Fore-arc extension and thrust-belt compression are seen to have involved two main phases of activity, separated by a regional unconformity recording a mid-Pliocene (3.5–3.0 Ma) tectonic reorganization. The two sites of mud volcanism lie in contrasting tectonic settings (inner fore-arc basin vs central fold-and-thrust belt) and record differing forms of seabed extrusive activity (twin mud cones and a caldera vs a broad mud pie). At both sites, subsurface data show that mud volcanism took place throughout the second tectonic phase, since the late Pliocene; differing forms of mud extrusion were accompanied by subsidence to form depressions beneath and within extrusive edifices up to 1.5 km thick. The basal subsidence depressions point to sources within the succession of thrusts underlying the inner to central Arc, consistent with microfossils within cored mud breccias from both sites that are derived from strata as old as Late Cretaceous.     

Seafloor distribution and last glacial to postglacial activity of mud volcanoes on the Calabrian accretionary prism,Ionian Sea

Mud volcanoes (MVs) are abundant along the eastern Mediterranean subduction zones, recording mud breccia extrusion over long timescales (10⁶ years), but to date relatively few have been recognised in the northern Ionian Sea on the Calabrian accretionary prism (CAP). In the present study, the seafloor distribution and recent activity of MVs is investigated across a 35,600 km² sector of the CAP using a regional acoustic dataset (multibeam bathymetric and backscatter imagery, integrated with subbottom profiles) locally ground-truthed by sediment cores. A total of 54 MVs are identified across water depths of 150–2,750 m using up to four geophysical criteria: distinctive morphology, high backscatter, unstratified subbottom facies and, in one case, a hydroacoustic flare. Fourteen MVs are identified from 3–4 criteria, of which five have been previously proven by cores containing mud breccia beneath up to 1.6 m of hemipelagic sediments (Madonna dello Ionio MVs 1–3, Pythagoras MV and the newly named Sartori MV), while nine others are identified for the first time (Athena, Catanzaro, Cerere, Diana, Giunone, Minerva, ‘right foot’, Venere 1 and 2). Forty other as yet unnamed MVs are inferred from 1–2 geophysical criteria (three from distinctive morphology alone). All but one possible MV lie on the inner plateau of the CAP, landwards of the Calabrian Escarpment in a zone up to 120 km wide that includes the inner pre-Messinian wedge and the fore-arc basins, where they are interpreted to record the ascent from depth of overpressured fluids that interacted with tectonic structures and with evaporitic or shale seals within the fore-arc basins. The rise of fluids may have been triggered by post-Messinian out-of-sequence tectonism that affected the entire pre-Messinian prism, but Plio-Quaternary sedimentation rates and depositional styles support the inference that significant mud volcanism has taken place only on the inner plateau. Sedimentation rates across the CAP applied to a 12 khz sonar detection depth of 225 cm imply that all MVs with backscatter signatures (50 of 54) have erupted mud breccias within the last 56 ka, and within the last 12.5 ka in the fore-arc basins. Ages of eruption estimated from the depth of cored mud breccias at five MVs, and a seismo-stratigraphic relationship at a sixth, indicate episodes at the last glacial maximum ca. 20 ka BP and during the postglacial period. Eruptive episodes within the Calabrian MV province constitute recurrent geohazards, separated by longer periods of quiescent (subdued) fluid seepage that are likely to support gas hydrate formation and chemosynthetic ecosystems.     

Episodic Cenozoic tectonism and the development of the NW European ‘passive’ continental margin

The North Atlantic margins are archetypally passive, yet they have experienced post-rift vertical movements of up to kilometre scale. The Cenozoic history of such movements along the NW European margin, from Ireland to mid-Norway, is examined by integrating published analyses of uplift and subsidence with higher resolution tectono-stratigraphic indicators of relative movements (including results from the STRATAGEM project). Three episodes of epeirogenic movement are identified, in the early, mid- and late Cenozoic, distinct from at least one phase of compressive tectonism. Two forms of epeirogenic movement are recognised, referred to as tilting (coeval subsidence and uplift, rotations <1° over distances of 100s of Kilometres) and sagging (strongly differential subsidence, rotations up to 4° over distances <100 km). Each epeirogenic episode involved relatively rapid (<10 Ma) km-scale tectonic movements that drove major changes in patterns of sedimentation to find expression in regional unconformity-bounded stratigraphic units. Early Cenozoic tilting (late Paleocene to early Eocene, c. 60–50 Ma) caused the basinward progradation of shelf-slope wedges from elongate uplifts along the inner continental margin and from offshore highs. Mid-Cenozoic sagging (late Eocene to early Oligocene, c. 35–25 Ma) ended wedge progradation and caused the onset of contourite deposition in deep-water basins. Late Cenozoic tilting (early Pliocene to present, <4±0.5 Ma) again caused the basinward progradation of shelf-slope wedges, from uplifts along the inner margin (including broad dome-like features) and from offshore highs. The early, mid- and late Cenozoic epeirogenic episodes coincided with Atlantic plate reorganisations, but the observed km-scale tectonic movements are too large to be accounted for as flexural deflections due to intra-plate stress variations. Mantle–lithosphere interactions are implied, but the succession of epeirogenic episodes, of differing form, are difficult to reconcile with the various syn-to post-rift mechanisms of permanent and/or transient movements proposed in the hypothetical context of a plume beneath Iceland. The epeirogenic movements can be explained as dynamic topographic responses to changing forms of small-scale convective flow in the upper mantle: tilting as coeval upwelling and downwelling above an edge-driven convection cell, sagging as a loss of dynamic support above a former upwelling. The inferred Cenozoic succession of epeirogenic tilting, sagging and tilting is proposed to record the episodic evolution of upper mantle convection during ocean opening, a process that may also be the underlying cause of plate reorganisations. The postulated episodes of flow reorganisation in the NE Atlantic region have testable implications for epeirogenic movements along the adjacent oceanic spreading ridge and conjugate continental margin, as well as on other Atlantic-type ‘passive’ margins.     

Neogene stratigraphy and the sedimentary and oceanographic development of the NW European Atlantic margin

A regional correlation of Neogene stratigraphy has been attempted along and across the NW European Atlantic continental margin, between Mid-Norway and SW Ireland. Two unconformity-bounded successions are recognised. These are referred to as the lower and upper Neogene successions, and have been dated as Miocene–early Pliocene and early Pliocene–Holocene, respectively, in age. Their development is interpreted to reflect plate-wide, tectonically driven changes in the sedimentary, oceanographic and latterly climatic evolution of the NE Atlantic region. The lower Neogene succession mainly preserves a record of deep-water sedimentation that indicates an expansion of contourite sediment drifts above submarine unconformities, within this succession, on both sides of the eastern Greenland–Scotland Ridge from the mid-Miocene. This is interpreted to record enhanced deep-water exchange through the Faroe Conduit (deepest part of the Southern Gateway), and can be linked to compressive inversion of the Wyville–Thomson Ridge Complex. Thus, a pervasive, interconnected Arctic–North Atlantic deep-water circulation system is a Neogene phenomenon. The upper Neogene succession records a regional change, at about 4 Ma, in the patterns of contourite sedimentation (submarine erosion, new depocentres) coeval with the onset of rapid seaward-progradation of the continental margin by up to 100 km. This build-out of the shelf and slope is inferred to record a marked increase in sediment supply in response to uplift and tilting of the continental margin. Associated changes in deep-water circulation may be part of an Atlantic-wide reorganisation of ocean bottom currents. Glacial sediments form a major component of the prograding shelf margin (shelf-slope) sediment wedges, but stratigraphic data indicate that the onset of progradation pre-dates significant high-latitude glaciation by at least 1 Ma, and expansive Northern Hemisphere glaciation by at least 3 Ma.     

Distribution and geological control of mud volcanoes and other fluid/free gas seepage features in the Mediterranean Sea and nearby Gulf of Cadiz

Existing knowledge on the distribution of mud volcanoes (MVs) and other significant fluid/free gas-venting features (mud cones, mud pies, mud-brine pools, mud carbonate cones, gas chimneys and, in some cases, pockmark fields) discovered on the seafloor of the Mediterranean Sea and in the nearby Gulf of Cadiz has been compiled using regional geophysical information (including multibeam coverage of most deepwater areas). The resulting dataset comprises both features proven from geological sampling, or in situ observations, and many previously unrecognized MVs inferred from geophysical evidence. The synthesis reveals that MVs clearly have non-random distributions that correspond to two main geodynamic settings: (1) the vast majority occur along the various tectono-sedimentary accretionary wedges of the Africa-Eurasia subduction zone, particularly in the central and eastern Mediterranean basins (external Calabrian Arc, Mediterranean Ridge, Florence Rise) but also along its westernmost boundary in the Gulf of Cadiz; (2) other MVs characterize thick depocentres along parts of the Mesozoic passive continental margins that border Africa from eastern Tunisia to the Levantine coasts, particularly off Egypt and, locally, within some areas of the western Mediterranean back-arc basins. Meaningfully accounting for MV distribution necessitates evidence of overpressured fluids and mud-rich layers. In addition, cross-correlations between MVs and other GIS-based data, such as maps of the Messinian evaporite basins and/or active (or recently active) tectonic trends, stress the importance of assessing geological control in terms of the presence, or not, of thick seals and potential conduits. It is contended that new MV discoveries may be expected in the study region, particularly along the southern Ionian Sea continental margins.     

A comparison of the NW European glaciated margin with other glaciated margins

The present paper provides an overview of glacial related seabed features and sedimentary sequences found along the formerly glaciated NW European margin and compare it with those found on contemporary glaciated margins from both the Southern and Northern Hemispheres. A brief review of the seabed physiography and strata architecture of the margins under consideration is followed by comparison of the most relevant similarities and differences. Comparison of the present-day bathymetric setting of both former and contemporary glaciated margins reveals no clear link to the effect of neither ice sheet or sediment load. Three different types of glacially eroded shelf transverse troughs have been identified, while marginal troughs seem connected to similar geological settings everywhere. Beyond the shelf edge interaction between downslope and alongslope processes has occurred resulting, amongst others, in the formation of large sedimentary mounds on the rise. More frequent large-scale mass wasting occurs on the former glaciated NW European margin than the Greenland and Antarctic margins in the latest Neogene to recent times. A two-stage evolution of the shelf prograding wedges is observed on all margins under consideration, which may reflect a general development of an ice cover from an initial phase of non- to restricted glaciation, evolving to a mature stage of expansive glaciation.     

Anomalous Cenozoic subsidence along the ‘passive’ continental margin from Ireland to mid-Norway

Two-dimensional flexural backstripping and thermal modelling (assuming uniform stretching and cooling) is applied to four interpreted, depth-converted seismic profiles across the Rockall, Faroe–Shetland and Vøring basins, along 1600 km of the Atlantic continental margin of NW Europe. The results reveal a significant discrepancy between the modelled palaeo-depths for the base of the Cenozoic succession and those proven by geological evidence at control points (subaerial conditions or depositional depth ranges in wells). The discrepancy is of Rm-scale, much larger than the possible range of parameter error determined by sensitivity tests (up to 0.5 km). Assuming a Cretaceous rift episode (100 Ma), the discrepancy is at least 1.7 km in the Rockall Basin, up to 2.1 km in the Faroe–Shetland Basin and at least 1 km in the Vøring Basin (which also contains evidence of kilometre-scale uplift of the inner margin). Assuming (unproven) a second rift in the early Cenozoic (60 Ma), the discrepancy remains of kilometre-scale in the Rockall and Faroe–Shetland basins. The restorations also provide evidence of uplift, both above compressive structures and across the modelled profiles as seaward rotations of palaeo-bathymetric records. The palaeo-bathymetric discrepancy corresponds to an anomaly in subsidence that is the cumulative product of all the tectonic episodes that have affected the NW European margin, and may incorporate both permanent effects of the last episode of lithospheric extension and transient responses to the interaction of the margin with mantle convective flow. Any explanation must accommodate both the large magnitude of anomalous subsidence along the margin and evidence of its episodic character.     

Maximum extent and readvance dynamics of the Irish Sea Ice Stream and Irish Sea Glacier since the Last Glacial Maximum

The BRITICE-CHRONO Project has generated a suite of recently published radiocarbon ages from deglacial sequences offshore in the Celtic and Irish seas and terrestrial cosmogenic nuclide and optically stimulated luminescence ages from adjacent onshore sites. All published data are integrated here with new geochronological data from Wales in a revised Bayesian analysis that enables reconstruction of ice retreat dynamics across the basin. Patterns and changes in the pace of deglaciation are conditioned more by topographic constraints and internal ice dynamics than by external controls. The data indicate a major but rapid and very short-lived extensive thin ice advance of the Irish Sea Ice Stream (ISIS) more than 300 km south of St George's Channel to a marine calving margin at the shelf break at 25.5 ka; this may have been preceded by extensive ice accumulation plugging the constriction of St George's Channel. The release event between 25 and 26 ka is interpreted to have stimulated fast ice streaming and diverted ice to the west in the northern Irish Sea into the main axis of the marine ISIS away from terrestrial ice terminating in the English Midlands, a process initiating ice stagnation and the formation of an extensive dead ice landscape in the Midlands.