Review Article: Permafrost Trapped Natural Gas in Svalbard, Norway

Permafrost has become an increasingly important subject in the High Arctic archipelago of Svalbard. However, whilst the uppermost permafrost intervals have been well studied, the processes at its base and the impacts of the underlying geology have been largely overlooked. More than a century of coal, hydrocarbon and scientific drilling through the permafrost interval shows that accumulations of natural gas trapped at the base permafrost is common. They exist throughout Svalbard in several stratigraphic intervals and show both thermogenic and biogenic origins. These accumulations combined with the relatively young permafrost age indicate gas migration, driven by isostatic rebound, is presently ongoing throughout Svalbard. The accumulation sizes are uncertain, but one case demonstrably produced several million cubic metres of gas over eight years. Gas encountered in two boreholes on the island of Hopen appears to be situated in the gas hydrate stability zone and thusly extremely voluminous. While permafrost is demonstrably ice-saturated and acting as seal to gas in lowland areas, in the highlands it appears to be more complex, and often dry and permeable. Svalbard shares a similar geological and glacial history with much of the Circum-Arctic meaning that sub-permafrost gas accumulations are regionally common. With permafrost thawing in arctic regions, there is a risk that the impacts of releasing of sub-permafrost trapped methane is largely overlooked when assessing positive climatic feedback effects.

Numerical solutions of the flexure equation

The 4th order differential equation describing elastic flexure of the lithosphere is one of the cornerstones of geodynamics, key to understanding topography, gravity, glacial isostatic rebound, foreland basin evolution and a host of other phenomena. Despite being fully formulated in the 1940’s, a number of significant issues concerning the basic equation have remained overlooked to this day. We first explain the different fundamental forms the equation can take and their difference in meaning and solution procedures. We then show how numerical solutions to flexure problems in general as they are currently formulated, are potentially unreliable in an unpredictable manner for cases where the coefficient of rigidity varies in space due to variations of the elastic thickness parameter. This is due to fundamental issues related to the numerical discretisation scheme employed. We demonstrate an alternative discretisation that is stable and accurate across the broadest conceivable range of conditions and variations of elastic thickness, and show how such a scheme can simulate conditions up to and including a completely broken lithosphere more usually modelled as an end loaded, single, continuous plate. Importantly, our scheme will allow breaks in plate interiors, allowing for instance, the creation of separate blocks of lithosphere which can also share the support of loads. The scheme we use has been known for many years, but remains rarely applied or discussed. We show that it is generally the most suitable finite difference discretisation of fourth order, elliptic equations of the kind describing many phenomena in elasticity, including the problem of bending of elastic beams. We compare the earlier discretisation scheme to the new one in 1 dimensional form, and also give the 2 dimensional discretisation based on the new scheme.We also describe a general issue concerning the numerical stability of any second order finite difference discretisation of a fourth order differential equation like that describing flexure where contrasting magnitudes of coefficients of different summed terms lead to round off problems which in turn destroy matrix positivity. We explain the use of 128 bit, floating point storage for variables to mitigate this issue.

Hammerfest Basin Composite Tectono-Sedimentary Element

The Hammerfest Basin is an E -W trending graben located between the Loppa High and the Finnmark Platform in the southern part of the Norwegian Barents Sea. Mainly siliciclastic strata of Carboniferous to Cenozoic age cover the Caledonian basement and have a total estimated thickness of 5-8 km. The basin evolved through several tectonic phases: the Carboniferous rifting, Late Jurassic rifting, the opening of the Atlantic Ocean, Oligocene reorganisation of plate movements and postglacial isostatic rebound. An E-W trending dome in the centre of the basin developed during the main extensional tectonic event in Late Jurassic. Horst structures represent the main hydrocarbon traps. Erosional channels on the flanks of the basin represent entry points for Lower Cretaceous sands. For the rest of the Cretaceous and Cenozoic intervals no significant reservoir sands are expected.The first exploration well in the Barents Sea in 1980 was located in the Hammerfest basin, and by 2019 a total of 45 wells had been drilled in the basin where 34 are classified as exploration wells. The result is 18 oil and gas discoveries, which gives a discovery rate of 53%. Two fields are now in production: the Snøhvit gas-condensate fields and the Goliat oil field.A total of 340 Msm3 (2140 Mbbl) recoverable oil equivalents have been discovered. For the middle Jurassic Play, the yet-to-find potential may be around 50 Msm3, distributed in several small structures in the basin. Following the oil discovery in the Middle Triassic interval in the Goliat structure, and because several of the previously drilled structures only penetrated the Jurassic and the uppermost Triassic section, considerable exploration potential may exist in the deeper Triassic interval in structures with the best reservoir facies. Stratigraphic traps of Cretaceous age may have a moderate petroleum potential, with excellent reservoirs encountered along the flank of the basin. Exploration potential may also exist in Upper Permian sandstones along the southern and eastern flanks of the basin. However, in large parts of the basin, the remaining potential is in the deep structures and hence is gas prone.

Escarpment retreat on a viscoelastic lithosphere

<p>Rift escarpments have long been the subject of coupled geodynamic/landscape evolution studies.  Many of these studies have shown that the flexural unloading response of the lithosphere plays a significant role in the rate of divide migration and the longevity of the escarpment topography, with lower elastic thickness values allowing for more localized isostatic rebound of the lithosphere in response to erosional unloading. However, rift escarpments are thought to last hundreds of millions of years, and therefore the lithosphere may exhibit viscoelastic behavior on this timescale.  Here we present a simplified model of a viscoelastic response to erosional unloading during escarpment evolution, and show that this drastically alters the behavior of the escarpment system.  Specifically, the escarpment retreat rate is significantly reduced, and topography maintained, when compared to a purely flexural model. Additionally, the area in front of the retreating the scarp (i.e., seaward of the scarp) experience delayed uplift response and topographic rejuvenation many millions of years after the divide passes.      </p>

Glacial isostatic adjustment near the center of the former Patagonian Ice Sheet (48°S) during the last 16.5 kyr

<p>Our understanding of glacial isostatic rebound across Patagonia is highly limited, despite its importance to constrain past ice volume estimates and better comprehend relative sea-level variations. With this in mind, our research objective is to reconstruct the magnitude and rate of Late Glacial to Holocene glacial isostatic adjustment near the center of the former Patagonian Ice Sheet. We focus on Larenas Bay (48°S; 1.26 km<sup>2</sup>), which is connected to Baker Channel through a shallow (<em>ca.</em> 7.4 m) and narrow (<em>ca.</em> 150 m across) inlet, and hence has the potential to record periods of basin isolation and marine ingression. The paleoenvironmental evolution of the bay was investigated through a sedimentological analysis of a 9.2 m long, radiocarbon-dated, sediment core covering the last 16.8 cal. kyr BP. Salinity indicators, including diatom paleoecology, alkenone concentrations and CaCO<sub>3</sub> content, were used to reconstruct the bay’s connectivity to the fjord. Results indicate that Larenas Bay was a marine environment before 16.5 cal. kyr BP and after 9.1 cal. kyr BP, but that it was disconnected from Baker Channel in-between. We infer that glacial isostatic adjustment outpaced global sea-level rise between 16.5 – 9.1 cal. kyr BP. During the Late Glacial - Holocene transition, the center of the former Patagonian Ice Sheet rose <em>ca.</em> 96 m, at an average rate of 1.30 cm/year. During the remainder of the Holocene, glacial isostatic adjustment continued (<em>ca.</em> 19.5 m), but at a slower average pace of 0.21 cm/year. Comparisons between multi-centennial variations in the salinity indicators and existing records of global sea-level rise suggest that the glacial isostatic adjustment rate fluctuated during these time intervals, in agreement with local glacier dynamics. More specifically, most of the glacial isostatic adjustment registered between 16.5 – 9.1 cal. kyr BP seems to have occurred before meltwater pulse 1A (14.5 – 14.0 kyr BP). Likewise, it appears that the highest Holocene glacial isostatic rebound rates occurred during the last 1.4 kyr, most likely in response to glacier recession from Neoglacial maxima. This implies a relatively rapid response of the local solid earth to ice unloading, which agrees with independent modelling studies investigating contemporary uplift. We conclude that the center of the former Patagonian Ice Sheet experienced a glacial isostatic adjustment of <em>ca.</em> 115 m over the last 16.5 kyr, and that >80% occurred during the Late Glacial and early Holocene.</p>

The Late-Quaternary morphostratigraphic record of glaciation, paleoseismicity and postglacial environmental changes in fjord-lakes of Québec-Labrador

<p>Hydroacoustic surveys were conducted in eight fjord-lakes of Québec-Labrador in order to analyse their Late-Quaternary geomorphological and stratigraphic record of glaciation, paleoseismicity and postglacial environmental changes. This large morphostratigraphic dataset provided a unique opportunity to establish a conceptual model of the evolution for fjord-lakes in relation to deglaciation, glacio-isostatic rebound, sediment fluxes and paleoseismicity. The analysis of the morphology and distribution of many morainic deposits into the fjord-lakes (hummocky moraines, morainic sills and morainic complexes) allows relating their formation to the glacial erosion potential, as well as to climatic and topographic controls. During past glaciations, a topographic sill was left uneroded at the opening of valleys due to the decrease in the glacial erosion potential associated with the lateral extension of the glacier down-ice; this bedrock sill created in turn an anchoring point to the ice during deglaciation. Hummocky moraines were documented at the outlet of five fjord-lakes that are located within the deepest and narrowest valleys of the studied systems. Based on our analysis of these sublacustrine landform-sediment assemblages, fjord-lakes constitute distinct sedimentary systems that should be differentiated from typical fjord system (i.e., in marine waters). The large-scale landforms contained in the fjord-lakes of Québec-Labrador (i.e., esker, moraines, gullies, lateral banks, turbidity channels and circular cavities) are inherited from their past subglacial, glaciomarine and paraglacial conditions, while only small deltaic bedforms (i.e., sediment waves and crescent-shaped bedforms) were formed in postglacial times. The present-day hydrological regime of fjord-lakes of Québec-Labrador is considered river-driven, except for the lakes located near active seismic zones where widespread postglacial mass-movements are documented.</p>

Cryptic terrestrial fungus-like fossils of the early Ediacaran Period

The colonization of land by fungi had a significant impact on the terrestrial ecosystem and biogeochemical cycles on Earth surface systems. Although fungi may have diverged ~1500–900 million years ago (Ma) or even as early as 2400 Ma, it is uncertain when fungi first colonized the land. Here we report pyritized fungus-like microfossils preserved in the basal Ediacaran Doushantuo Formation (~635 Ma) in South China. These micro-organisms colonized and were preserved in cryptic karstic cavities formed via meteoric water dissolution related to deglacial isostatic rebound after the terminal Cryogenian snowball Earth event. They are interpreted as eukaryotes and probable fungi, thus providing direct fossil evidence for the colonization of land by fungi and offering a key constraint on fungal terrestrialization.

The “Passive” Margin of Eastern North America: Rifting and the Influence of Prerift Orogenic Activity on Postrift Development

We have analyzed and synthesized geologic and geophysical data from the onshore Newark rift basin and adjacent onshore and offshore basins to better understand the Mesozoic development of the eastern North American rift system and passive margin. Our work indicates that rifting had three phases: (1) an initial, prolonged phase of extension and subsidence; (2) a short-lived phase with higher rates of extension and subsidence, intrabasin faulting, and intense magmatism; and (3) a final phase with limited subsidence and deposition. Additionally, our work shows that anomalous uplift and erosion, associated with crustal-scale arching/warping subparallel to the prerift and syn-rift crustal fabric not the continent-ocean boundary, affected a region landward of the basement hinge zone. Uplift and erosion began during the final rifting phase and continued into early drifting with erosion locally exceeding 6 km. Subsequent subsidence was minimal. We propose that denudation unloading related to relic, prerift orogenic crustal thickness and elevated topography triggered the anomalous uplift and erosion. After the Paleozoic orogenies, postorogenic denudation unloading (cyclic erosion and isostatic rebound/uplift) significantly thinned the thickened crust and reduced topographic elevation. During rifting, extension stretched and tectonically thinned the crust, promoting widespread subsidence and deposition that dampened the postorogenic cycle of erosion and isostatic rebound/uplift. During the rift-drift transition, with extension focused near the breakup site, denudation unloading resumed landward of the basement hinge zone, producing significant erosion and uplift (related to isostatic rebound), crustal thinning, and topographic decay that left behind only eroded remnants of the once massive rift basins.

Fluvial incised networks on top of late Ordovician interglacial valley buried hills as the result of post glacial isostatic rebound; 3D seismic input

<p>It is classically assume that prior to deep glacial valleys incision below large scale ice cap, often interpreted as the results of ice flow melting during tidewater period, the initial glacial topography was flat or very low angle and created during a major phase of cold glaciers advance as suggested by quaternary studies. Therefore up to now we have assume that the top of late Ordovician buried hills separating major glacial valleys was the remains of this flat surface truncating the pre-glacial Ordovician Hawaz series, later on flooded by the Lower Silurian. Surprisingly by reinterpreting 3D seismic cubes using spectral decomposition technics on the Murzuk basin in SE Libya, it appears that the top of buried hills are not at all characterized by a flat erosional surface, but it is strongly irregular and shows the development of narrow valley networks displaying the classical dendritic erosional pattern diagnostic of fluvial erosion. These small valleys are organized into a tributary network and don’t flow toward the ice margin, i.e. toward the N-NW but most of the time flow at right angle toward the adjacent main glacial valleys which are pointing toward the NW. These narrow valley networks in this context could be either glacial tunnel valleys located at the periphery of the ice cap in close relationships with glacial fronts (their common settings) or could correspond to fluvial valleys developed later on, in a subaerial setting at some distance from glacial fronts; we retain this second interpretation because in addition to the geomorphic features: (1) they flow parallel to the fronts that we have already recognized, Moreau et al. (2005), Rubino et al.  2007 and (2) they are suspended in the sense that these lateral networks do not reach the bottom of the main glacial valley but, they appear to be connected within the upper part of the glacial infill, immediately below the early Silurian post glacial flooding characterized by the well-known Rhuddanian hot shales. As a result, the incision of the valley network appears quite late in the ice cap melting history. It is why we tend to interpret these valleys erosion as the result of post glacial melting during ice retreat at some distance from the ice front and strongly enhanced by isostatic rebound. Some possible modern analogs of such valley fringing highs may exist in Artic Canadian islands.</p>