Regional Exploration and Characterisation of CO2 Storage Prospects in the Utsira-Skade Aquifer, North Viking Graben, North Sea

Subsurface CO2 storage is considered a key element of reducing anthropogenic emissions in virtually all scenarios compatible with limiting global warming to 1.5°C. The Utsira-Skade Aquifer (Utsira, Eir and Skade Formations), northern North Sea, has been identified as a suitable reservoir. Although the overall storage capacity of the full aquifer has been estimated based on regional data, it is lacking an integrated assessment of containment and internal heterogeneity, to identify optimal areas for injection and for calculation of site-specific storage capacities. A high-resolution, broadband 3D seismic reflection dataset, full waveform inverted velocity data and 102 exploration wells are utilised to provide a catalogue of CO2 storage prospects in the northern Utsira-Skade Aquifer. This is achieved through: 1) definition of the aquifer’s spatial limits; 2) calculation of porosity distribution; 3) assessment of the extent, geomorphology, thickness variability, and containment confidence (CC) of mudstones; and 4) mapping of closures through fill-to-spill simulations. CO2 storage capacity was calculated for the prospects using two approaches; using the full reservoir thickness (FRT) beneath the closures and using only the thickness from the closure top to the spill point (TSP), i.e., within structural traps. Porosity ranges from 29 to 39% across the aquifer and is higher in the Utsira and Eir Fms. relative to the underlying Skade Fm. The mudstone separating the Skade and Eir/Utsira Fm. has a thickness > 50 m, and is a potential barrier for CO2. Other intra-aquifer mudstones were mainly interpreted to act as baffles to flow. Structural traps at the top Utsira and Skade Fms. yield fifteen prospects, with criteria of > 700 m depth and FRT storage capacity of > 5 Mt CO2. They have a combined storage capacity of 330 Mt CO2 (FRT) or 196 Mt CO2 (TSP). Five prospects have a positive CC score (total capacity: 54 Mt CO2 FRT or 39 Mt CO2 TSP). Additional storage capacity could be achieved through more detailed analysis of the seal to upgrade the CC scores, or through use of a network of the mapped closures with a fill-to-spill approach, utilising more of the aquifer.

Paleozoic basin reactivation and inversion of the underexplored northern North Sea platforms: a cross-border approach

Well penetrations on the UK East Shetland Platform (ESP) prove 1-8 km thick Devonian post-orogenic extensional collapse-related successions. Conversely, extremely thick (1-6 km) Permo-Triassic basin fills without Devono-Carboniferous units were in the past interpreted west of the Utsira High, on the Norwegian Horda Platform and Stord Basin, albeit Pre-Triassic well penetrations are here very rare. In this work, the nature and age of Paleozoic-Triassic strata and structures in these underexplored platform regions are tentatively constrained by performing cross-border regional seismic interpretation east and west of the Viking Graben.We highlight cross-border analogies in structural style and seismic facies, with a similar evolution dominated by polyphase inversion tectonics and structural grain preservation. In the Norwegian study areas, much of the half-graben sedimentary fills may be interpreted as Devonian-?Carboniferous in age as in the ESP, rather than overly thick Permo-Triassic successions. Major graben-bounding extensional faults are low-angle (∼25-33°), approximately northerly-striking and likely rooting downwards into reactivated Caledonian shear zones. Rifting development occurred in multiple episodes, possibly creating different traps. Prior to Permian-Jurassic rifting, many low-angle Caledonian thrusts were subject to extensional inversion in the Devonian and then to Variscan compressional reactivation, causing vertical extrusion and deformation of Devonian syn-rift wedges.

Resolving the spatial distributions of Dipturus intermedius and Dipturus batis—the two taxa formerly known as the ‘common skate’

Batoid fishes are among the most endangered marine vertebrates, yet conservation efforts have been confounded by incomplete taxonomy. Evidence suggest that the critically endangered ‘common skate’ actually represents two species: the flapper skate (Dipturus intermedius) and the blue skate (Dipturus batis). However, knowledge of the geographic range of these two nominal species is limited. Here, DNA sequencing is used to distinguish these species, allowing their spatial distributions to be clarified. These records were also used as the basis for species distribution modelling, providing the first broad scale models for each species across the Northeast Atlantic. Samples were obtained from Iceland, the UK (specifically Shetland), the North Sea and the Azores. Results suggest that D. batis was commonly distributed in the Western Approaches and Celtic Sea, extending out to Rockall and Iceland. D. intermedius generally appears to be less abundant, but was most frequent around northern Scotland and Ireland, including the northern North Sea, and was also present in Portugal. Two individuals were also identified from seamounts in remote areas of the Atlantic around the Azores, the furthest south and west the species has been found. This supports reports that the flapper skate historically had a much wider distribution (which was also highlighted in the distribution model), emphasising the large scale over which fisheries may have led to extirpations. Furthermore, these Azorean samples shared a unique control region haplotype, highlighting the importance of seamounts in preserving genetic diversity.

Weakening and warming of the European Slope Current since the late 1990s attributed to basin-scale density changes

Oceanic influences on shelf seas are mediated by flow along and across continental slopes, with consequences for regional hydrography and ecosystems. Here we present evidence for the variable North Atlantic influence on European shelf seas over the last four decades, using ocean analysis and reanalysis products, and an eddy-resolving ocean model hindcast. To first order, flows oriented along isobaths at the continental slope are related to the poleward increase of density in the adjacent deep ocean that supports a geostrophic inflow towards the slope. In the North Atlantic, this density gradient and associated inflow has undergone substantial, sometimes abrupt, changes in recent decades. Inflow in the range 10–15 Sv is identified with eastward transport in temperature classes at 30° W, in the latitude range 45–60° N. Associated with major subpolar warming around 1997, a cool and fresh branch of the Atlantic inflow was substantially reduced, while a warm and more saline inflow branch strengthened, with respective changes of the order 5 Sv. Total inflow fell from ~ 15 Sv pre-1997 to ~ 10 Sv post-1997. In the model hindcast, particle tracking is used to trace the origins of poleward flows along the continental slope to the west of Ireland and Scotland, before and after 1997. Backtracking particles up to 4 years, a range of subtropical and subpolar pathways is identified from a statistical perspective. In broad terms, cold, fresh waters of subpolar provenance were replaced by warm, saline waters, of subtropical provenance. These changes have major implications for the downstream shelf regions that are strongly influenced by Atlantic inflow, the northern North Sea in particular, where “subtropicalization” of ecosystems has already been observed since the late 1990s.

Integrating CVMix into GOTM (v6.0): a consistent framework for testing, comparing, and applying ocean mixing schemes

The General Ocean Turbulence Model (GOTM) is a one-dimensional water column model, including a set of state-of-the-art turbulence closure models, and has widely been used in various applications in the ocean modeling community. Here, we extend GOTM to include a set of newly developed ocean surface vertical mixing parameterizations of Langmuir turbulence via coupling with the Community Vertical Mixing Project (CVMix). A Stokes drift module is also implemented in GOTM to provide the necessary ocean surface waves information to the Langmuir turbulence parameterizations, as well as to facilitate future development and evaluation of new Langmuir turbulence parameterizations. In addition, a streamlined workflow with Python and Jupyter notebooks is also described, enabled by the newly developed and more flexible configuration capability of GOTM. The newly implemented Langmuir turbulence parameterizations are evaluated against theoretical scalings and available observations in four test cases, including an idealized wind-driven entrainment case and three realistic cases at Ocean Station Papa, the northern North Sea, and the central Baltic Sea, and compared with the existing general length scale scheme in GOTM. The results are consistent with previous studies. This development extends the capability of GOTM towards including the effects of ocean surface waves and provides useful toolsets for the ocean modeling community to further study the effects of Langmuir turbulence in a broader scope.

Rift kinematics preserved in deep-time erosional landscape below the northern North Sea

Our understanding of continental rifting is largely derived from the stratigraphic record. This archive is, however, incomplete as it does not capture the geomorphic and erosional record of rifting. New 3D seismic reflection data reveals a Late Permian-Early Triassic landscape incised into the pre-rift basement of the northern North Sea. This landscape, which covers at least 542 km2, preserves a drainage system bound by two major tectonic faults. A quantitative geomorphic analysis of the drainage system reveals 68 catchments, with channel steepness and knickpoint analysis of catchment-hosted paleo-rivers showing that the landscape preserved a >2 Myrs long period of transient tectonics. We interpret that this landscape records punctuated uplift of the footwall of a major rift-related normal fault at the onset of rifting. The landscape was preserved by a combination of relatively rapid subsidence in the hangingwall of a younger fault and burial by post-incision sediments. We show how and why erosional landscapes are preserved in the stratigraphic record, and how they can help us understand the tectono-stratigraphic evolution of ancient continental rifts.

Ommastrephid squid spawning in the North Sea: oceanography, climate change and species range expansion

The lesser flying squid (Todaropsis eblanae) and the shortfin squid (Illex coindetii) are two abundant ommastrephids of the northeast Atlantic. Spawning ground existence was inferred from the captures of mature, mated females in summer 2016–2019 and their occurrences were compared with respective oceanographic data from international surveys to gain insight into environmental predictors of their presence throughout the North Sea. Spawning T. eblanae were found in relatively cooler and more saline waters (6–8°C, 34.2–35.1 psu) in the northern North Sea linked to the Fair Isle Current and East Shetland Atlantic Inflow, whilst spawning I. coindetii occurred across the entire North Sea (mostly at 9–10.5°C, 34.1–34.8 psu). We hypothesize that a combination of water salinity and water temperature are key factors in the spatiotemporal distribution of spawning ommastrephid squids as they define water density that is crucial for pelagic egg mass to attain neutral buoyancy.