Next Article in Journal
Asymmetric Craters on the Dwarf Planet Ceres—Results of Second Extended Mission Data Analysis
Next Article in Special Issue
The Influence of Magmatic Intrusions on Diagenetic Processes and Stress Accumulation
Previous Article in Journal
The Role of Shear Wave Velocity and Non-Linearity of Soil in the Seismic Response of a Coupled Tunnel-Soil-Above Ground Building System
Previous Article in Special Issue
An Assessment of Stress States in Passive Margin Sediments: Iterative Hydro-Mechanical Simulations on Basin Models and Implications for Rock Failure Predictions
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Geosciences
Volume 9
Issue 11
10.3390/geosciences9110474
geosciences-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tilting and Flexural Stresses in Basins Due to Glaciations—An Example from the Barents Sea

by
Ingrid F. Løtveit
*,
Willy Fjeldskaar
and
Magnhild Sydnes
Tectonor AS, P.O. Box 8034, N-4068 Stavanger, Norway
*
Author to whom correspondence should be addressed.
Geosciences 2019, 9(11), 474; https://doi.org/10.3390/geosciences9110474
Submission received: 21 September 2019 / Revised: 6 November 2019 / Accepted: 7 November 2019 / Published: 11 November 2019
(This article belongs to the Special Issue Future Advances in Basin Modeling: Suggestions from Current Observations, Analyses, and Simulations)
Browse Figures
Review Reports Versions Notes

Abstract

:
Many of the Earth’s sedimentary basins are affected by glaciations. Repeated glaciations over millions of years may have had a significant effect on the physical conditions in sedimentary basins and on basin structuring. This paper presents some of the major effects that ice sheets might have on sedimentary basins, and includes examples of quantifications of their significance. Among the most important effects are movements of the solid Earth caused by glacial loading and unloading, and the related flexural stresses. The driving factor of these movements is isostasy. Most of the production licenses on the Norwegian Continental Shelf are located inside the margin of the former Last Glacial Maximum (LGM) ice sheet. Isostatic modeling shows that sedimentary basins near the former ice margin can be tilted as much as 3 m/km which might significantly alter pathways of hydrocarbon migration. In an example from the SW Barents Sea we show that flexural stresses related to the isostatic uplift after LGM deglaciation can produce stress changes large enough to result in increased fracture-related permeability in the sedimentary basin, and lead to potential spillage of hydrocarbons out of potential reservoirs. The results demonstrate that future basin modeling should consider including the loading effect of glaciations when dealing with petroleum potential in former glaciated areas.

1. Introduction

It is commonly accepted that major changes in Earth’s climate started in Gelasian (~2.6 million years ago), and initiated the growth of ice sheets in the Northern Hemisphere (e.g., [1]). Over the last 2.6 million years, there have probably been more than 30 glaciations with cycles of ice sheets advancing and retreating on time scales of 40,000 to 100,000 years. A glaciation cycle of 40,000 years dominated in the early period, while the 100,000 year cycle has dominated in the past million years or so [2]. The most extensive glaciations are traditionally referenced in the last million years (Late Quaternary; cf. [3]) and the covers on North America and Europe have been estimated to be around 2–4 km thick near the centers of maximum ice accumulation (cf. e.g., [4]), and tapered towards the ice sheet margins. Repeated episodes of growth and withdrawal of huge ice sheets, glacial erosion, and associated sediment deposition lead to geologically dramatic changes in surface loading. The lithosphere responds to these changes by subsiding under loading, and by uplifting when loading is removed. The driving mechanism for this is isostasy.
Differential vertical movement of the lithosphere related to glacial isostasy lead to repeated tilting of sedimentary formations and potential petroleum reservoirs therein, which may have greatly affected hydrocarbon migration pathways in the former glaciated areas [5]. In addition, the upward and downward bending of the lithosphere due to glaciation leads to flexural stresses [6] likely to affect faults and their permeability, which could add to the changes in hydrocarbon migration pathways.
In parts of the Barents Sea, there are clear indications of hydrocarbon spillage (e.g., [7]). Some of the hydrocarbon spillage is linked to isostatic movements due to cycles of Pliocene-Pleistocene ice sheet loading/unloading and glacial sediment redistribution [8,9,10,11,12]. It is suggested that the recent discovery of the giant Johan Sverdrup oil field in the Norwegian North Sea was oil charged during Quaternary, and that the area more than likely underwent tilting and possible leakage several times over the last one million years [13]. Other research, related to the North Sea, showed that influence of ice sheet loading does not result in a significant tilting [14], which may be the result of smaller gradients of the ice thickness. Anyway, detailed control on the glacial isostasy is important and so far an insufficiently utilized factor for identification of the remaining hydrocarbon resources in sedimentary basins formerly covered by ice sheets.
Faulting related to the last glaciation is reported from northern Scandinavia, where several fault with up to 30 m high fault scarps [15]. Examples include the 80 km long Stuoragurra fault in northern Norway [16] and the 150 km long Pärvie fault in northern Sweden [17]. There is also evidence of such faults in Denmark [18], Germany, Poland, and the United Kingdom ([19,20,21,22,23]). Thus, such faults are found in the center of former ice sheets, as well as at their margins and beyond [24,25].
For the offshore areas, the situation is more complicated. However, the number of seismic events recorded in the Barents Sea in the last 20 years is significant ([26,27]). For the Barents Sea, the impact of the glacial isostasy on triggering earthquakes has not yet been investigated. A better understanding of the origin of the seismicity will help to recognize which faults and related areas exhibit a risk of leakage of hydrocarbons from traps.
Stress changes induced by glaciations are related to two main effects: 1) direct load introduced by the weight of the ice on the surface and 2) flexural load caused by the subsidence and uplift of the lithosphere due to glacial isostasy [10,28]. Horizontal stress due to stress migration, which is depending on loading time and mantle viscosity, is an additional effect that is modeled by Steffen et al. [29]. Grollimund and Zoback [30] investigated, by a 2D approach, lithospheric flexure as an alternative mechanism to explain the local stress perturbations observed in the northern North Sea. They found that flexural stresses due to glaciations/deglaciation could explain the observed lateral stress variations.
Johnston [31] used the Mohr–Coulomb theory for the study of faulting and explained the lack of seismic activity in areas of large continental ice sheets. Turpeinen [32] use finite-element models of rheologically layered lithosphere with a thrust fault to investigate its slip evolution during the glacial/ postglacial period. Their modelling results suggest that the rate of thrusting decreases during ice sheet loading and strongly increases during deglaciation. Steffen et al. [33] introduced a new approach that is based on a Glacial Isostatic Adjustment (GIA) model. They combined the GIA model with a fault surface to investigate the fault slip and fault activation time during a glacial cycle. Their model allows estimation of the fault throw for former glaciated areas. One of their conclusions is that faults start to move close to the end of deglaciation.
Our approach is flexural isostasy modelling, and we are dealing with two effects of the glaciations which we believe future basin modelling has to consider:
  • Possible glacially induced tilting of reservoirs. Here, we calculate the glacially induced tilting of the Norwegian Continental Shelf due to the last ice age, with particular focus on the SW Barents Sea.
  • Assessment of the flexural stresses induced by the lithospheric flexure related to glacial isostasy and tilt, and their effect on faults during glaciations and interglacials, here illustrated by an example from the SW Barents Sea where the stress modelling is limited to the flexural stress effects of deglaciation after LGM. The stress effects of vertical ice load and related stress migration are not calculated here.

2. Glaciations

The last glacial period began about 110,000 years ago and ended 10,000 years ago. A model of the deglaciation history for our area is given in [34]. The glaciations that occurred during this glacial period covered many areas of the Northern Hemisphere. The extent of the Eurasian Ice Sheet during the Last Glacial Maximum (LGM) has been long debated in the literature. Researchers [35] argue for a highly asynchronous ice sheet with regional differences in the age of maximum ice-sheet positions. However, in this paper we use the LGM (~20,000 BP) reconstruction of [36], which suggests the SW Barents Sea was glaciated from 24,000 to 18,000 BP. In the following calculations, we assume the Earth reached isostatic equilibrium during this period.

2.1. Ice Thickness

Glacial isostasy is calculated based on the palaeo ice extent (from dated marginal moraines) and models of the ice thicknesses. Observational data increasingly constrain the extent of the Quaternary ice sheets. The thickness and volume of these ice sheets are much harder to reconstruct and generally need to be inferred from indirect evidence and modeling.
Several models of ice thickness of LGM and the late-glacial deglaciation are published; the loading scenarios are produced using different methods and sets of constraints. Some of them consider thermomechanical coupling (e.g., [37]). A number of models of the LGM ice thickness are presented in [38]. We follow the work of [36]; their LGM extent model is shown in Figure 1.

2.2. Glacial Isostasy and Tilting of Reservoirs

The Earth’s response to deglaciation shows that the lithosphere acts as an elastic shell overlying a viscous mantle. If a load is applied to the elastic lithosphere, part of the applied load will be supported by the elastic lithosphere, and part by buoyant forces of the mantle underneath, acting through the lithosphere.
In the literature, there is disagreement on the elastic thickness of the lithosphere; the elastic thickness of the lithosphere of Fennoscandia in GIA (Glacial Isostatic Adjustment) modeling varies from 30 to 160 km [39]. Our modeling of the post-glacial isostatic response, calibrated with the observations onshore Norway, gave best fit with an elastic thickness of 30 km; see [34]. Calibration with the Barents Sea gave similar results; see [40]. The method used in the isostatic calculations is described in [41] and in the Appendix. The spatial resolution is 10 km.
There will be significant subsidence of potential sedimentary basins formerly covered by ice sheets. Sedimentary basins that are located in the periphery of the former ice sheet will be notably affected by tilting of the Earth´s surface as the subsidence (or uplift after the melting) is gradually decreasing towards the periphery of the former ice sheet. Therefore, the largest tilts of the Earth´s surface are found in the more peripheral areas of former ice sheets. This is clearly seen along the entire Norwegian coast and in the Barents Sea (Figure 2). Sedimentary basins located in the peripheral areas in the SW Barents Sea could be tilted by up to 2.7 m/km due to the LGM ice sheet (Figure 2), and even more in the previous greater glaciations [3]. The ice sheets will act as seesaw during the total glaciation period (Quaternary); the Earth will subside under the ice sheets and rise again in interglacial periods.
This process is, of course, important for migration modelling. Figure 3 shows how a simple structure behaves during one glacial cycle. It is assumed, for illustration purposes, that the structure is filled with hydrocarbons to the spill point at the onset of glaciation. When the lithosphere and the structure subside due to the glacial load, the reservoir is tilted hydrocarbons start to spill out of the structure. This goes on until the maximum tilting is reached (Figure 3, upper right). The closure volume will be reduced.
When the lithosphere and structure are uplifted due to deglaciation, a new hydrocarbon-water contact (HCWC) will be established (Figure 3, bottom). Paleo-HCWCs can theoretically be found at the spill point level for the positions before and after deglaciation. Residual oil can be found in the water zone down to the original spill point level. In a gas-filled reservoir with an oil leg originally, residual oil can also be found in the gas-filled part, limited by the uppermost HCWC.
The effect of the glacially induced hydrocarbon migration is determined by several factors:
  • The number of glacial events. Tilting will result from each isostatic event. A number of glacial events will create chaotic patterns of paleo contacts and residual oil which will be very difficult to interpret.
  • The geometry of the receiving area. If a saddle area is opened, the fluids can migrate into higher lying traps or be lost to the surface. A surrounding flat area, on the other hand, can become part of the trap during the subsidence, and fluid can migrate back into the trap under uplift.
  • Additional potential factors that determine the migration effect are: i) the geometry of the deformation, ii) the geometry of the hydrocarbon filled structures, iii) the initial fill of the trap, and iv) the spill point location and orientation of the structure.
The isostatic response on the Pleistocene sediment redistribution and ice loading has been evaluated for the Bjørnøyrenna Fault Complex in the Barents Sea [11]. The isostatic impact on hydrocarbon trap capacity changes and hydrocarbon maximum spillage was assessed in that study and the conclusion was that tilting together with gas volume expansion might have been responsible for some part of the hydrocarbon loss during the Cenozoic.

3. Flexure-Related Stress Effects

As mentioned above, the flexural stresses that occur during glaciations will be related to the isostatic response of the lithosphere, which is a time-dependent process. The lithosphere is subsiding under the weight of ice, and uplifted back to original state during unloading. During this process, the stresses in the elastic part of the lithosphere will change.
The estimation of stress regime in the lithosphere due to flexural stress during glaciations depends on the timing of crustal stress equilibrium. Two elastic models have been proposed for the estimation of stress changes by flexural stresses during glaciations. Stephansson [42] assumes that the lithosphere is in equilibrium before the ice load is introduced (Figure 4, left). The resulting stress state during ice loading then becomes compressive underneath the ice sheet, and extensional outside the margin of the ice sheet (forebulge area). During ice withdrawal, the horizontal stresses slowly diminish until the uplift ceases and a new equilibrium is reached. Stein et al. [43] present a model where the lithosphere is assumed to have reached equilibrium after some time of ice loading (Figure 4, right). Thus, during melting and uplift the area originally covered by ice experiences extensional stresses, and the area outside the former margin of the ice sheet becomes subject to compressive stresses. These stress changes decrease with depth and are inverted in the lower half of the lithosphere.
If we assume that the duration of ice loading and unloading is long enough for the lithosphere to reach equilibrium during both glaciation and deglaciation, then both models are applicable during one glacial cycle. The average duration of glaciations and interglacials during Quaternary has been 40,000 to 100,000 years. The last ice age ended around 10,000 years ago, and at present the crust is nearly in isostatic equilibrium. The assumption of crustal equilibrium during both glaciations and interglacials is therefore not unrealistic, as the relaxation time for the Fennoscandian uplift is less than 4000 years [34,44].

3.1. Identification of Glacially Flexured Areas

To identify if a certain study area has been affected by glacially-induced flexure, we need information of the glaciation history. The lithospheric flexure is calculated based on the glacial isostasy (Figure 5), which is derived from the ice extent and thickness, as well as the properties and thickness of the lithosphere and mantle.
Here, we will provide numerical models of flexure-related stress effects in the Hoop Fault Complex, located in the SW Barents Sea (Figure 6). The modelling will be limited to the stress effects related to the isostatic response during the deglaciation after LGM. For identifying the areas of SW Barents Sea that are affected by flexural stresses, we use the glacial isostasy of the LGM ice sheet calculated in Section 2 as input.
Based on the resulting glacial isostasy, we can identify the areas of maximum bending of the crust, and the hinge zones with minimum flexure. The derivation results of a function are illustrated in Figure 5. The first function illustrates the glacial isostasy on a vertical section through the lithosphere near the ice margin. The gradient (or tilt) can be derived from the first derivative of the glacial isostasy, the curvature or flexure of the lithosphere can be found by the second derivative.
The second derivative of the glacial isostasy of the LGM ice for the SW Barents Sea is shown in Figure 6. According to this, it is clear that the SW Barents Sea was greatly affected by lithospheric flexure during the deglaciation after LGM. There is a pronounced hinge zone (white area) extending from west of Bjørnøya over the Stappen High and into Loppa High. From here it extends eastwards over the Samson Dome onto the Finnmark Platform. Areas in the hinge zone are less likely to be affected by the flexural stress after the LGM deglaciation.

3.2. Model Setup

The stress modelling is done in 2D using the state of plane stress, and carried out using the commercial software Comsol Multiphysics 5.2 and the add-on structural mechanics module. Comsol Multiphysics is a FEA (Finite Element Analysis) software, where the problem domain is divided by a mesh of triangular elements in which the calculations are done. The set-up of a stress model in Comsol Multiphysics includes defining model geometry, then adding material properties, loads and boundary conditions.
The model geometry includes faults geometries based on a published continental shelf map of Norway [45] with scale 1:2 mill, interpreted at the base of Top Jurassic. The model includes an area of approximately 300 × 300 km in the SW Barents Sea that was covered by ice during LGM. All fault geometries are manually traced in a CAD tool, and then imported into Comsol Multiphysics.
Every geometry unit of the model is assigned material properties, that is, Young´s modulus and Poisson´s ratio for the lithologies included. In the models presented here, the host rock is modelled as soft shale with a Young´s modulus of 5 GPa and a Poisson´s ratio of 0.25, which is within the typical values for shales [46]. Fault zones normally consist of numerous fractures and smaller faults which lower the general stiffness, or Young´s modulus, of the host rock [47]. The faults are therefore given a Young´s modulus of 0.1 GPa and Poisson´s ratio of 0.25.
Next, boundary conditions are assigned to the model. As the model represents the base of a geological unit at depth, the sides of the model cannot deform freely. Two opposite boundaries of the model are therefore fixed. The other two boundaries are assigned a load representing the flexure-related extension originating from the isostatic response during the deglaciation after LGM. The orientation of the flexure-related extension is perpendicular to the curvature trends, with an orientation ~SW–NE. (Figure 6). In a study of the post-glacial lithospheric flexure in the North Sea, analytical and numerical calculation of glacial unloading (assuming isostatic equilibrium before the onset of unloading) of a 1 km thick ice sheet and a lithosphere with elastic thickness of 30 km resulted in bending stresses in the order of 20 MPa [30]. We have used 5 MPa, which is a more conservative magnitude for this effect. The model does not take into consideration possible tectonic background stress at the time of glacial unloading. This modelling may therefore be considered as a study of the effect of stress changes due to the glacially induced flexure.

3.3. Results

In the resulting models (Figure 7), the stress is calculated based on the assumption of a linear elastic material. Rocks normally behave elastic up to 1%–2% strain at low temperature and pressure [46,48,49]. The elastic behavior of rocks is limited by their strength, that is, the maximum applied compressive or tensile stress the rock can withstand before failure occurs. The tensile strength of most rocks is in the range 0.5 to 6 MPa, most commonly 2–3 MPa [46,50]. This means that crustal tensile stress exceeding these values is not very common, as it is likely to result in tensile failure of the rock. The shear strength of rocks is twice their tensile strength, as follows from the Griffith failure criterion [46,51]. Thus, in the models, the tensile stress is truncated at 6 MPa, and the shear stress at 12 MPa. Areas concentrating stress above these values are likely to fail in tension or shear. However, due to uncertainties related to the material properties for the involved sedimentary units and the flexure-related load, the stress magnitudes in the models below should be read in terms of relative values. The stress is likely to be released by failure in the high magnitude areas (marked in yellow and red), and less likely to lead to failure in low magnitude areas (marked in blue).
The first model (Figure 7, top) shows the resulting maximum principal tensile stress (σ3) when the area is subject to extension during LGM deglaciation. Tensile stress can be described as stretching stress, and can lead to opening up of new extension fractures (when the tensile strength of the rock is exceeded) or opening up of preexisting fractures. The SW–NE extension causes large areas between the faults to concentrate high magnitude tensile stress, whereas other parts of the host rock are located in a “stress shadow”. Here, the host rock surrounding SW–NE trending faults accumulate high magnitude tensile stresses, whereas the fault zones themselves show minimal effect. This indicates that due to the LGM deglaciation, the host rock surrounding SW–NE trending faults is likely to experience increased local fracture-related permeability, which may lead potential hydrocarbons towards and along the faults. Focusing on the exploration wells, the result suggests that wells 7324/2-1 and 7324/8-2 are placed within areas of high fracture-related permeability during the LGM deglaciation. Potential hydrocarbon may have migrated along faults out of a potential reservoir during the deglaciation. According to the results, the other wells are located in areas of low tensile stress, and therefore less likely to have been affected by flexural stresses during and after LGM.
The von Mises shear stress results (Figure 7, bottom) show similar tendency. Most of the wells are located in areas of low magnitude shear stress, whereas wells 7324/2-1 and 7324/8-2 are located in areas of higher shear stress magnitudes. Well data for exploration wells drilled in the area show that all wells included here are oil or gas discoveries, except for well 7324/2-1 and 7324/8-2 which are dry with shows, indicating former presence of hydrocarbons.

4. Discussion

4.1. Glaciation Model

Our glaciation model is not intended to include time-varying ice sheet geometries. However, LGM was the ice sheet configuration that was latest influencing the SW Barents Sea. The SW Barents Sea was glaciated around 24,000 BP and starting to melt away at 18,000 BP [36]. The LGM ice sheet is thus a reasonable ice sheet configuration for modelling the flexural stress effects in SW Barents Sea.

4.2. Earth Model

In the isostatic calculations we have assumed an elastic lithosphere thickness of 30 km (flexural rigidity 5 × 1023 Nm). As mentioned above this Earth model has been shown to match other Fennoscandian post-glacial data, but most of the published GIA Earth models of Fennoscandia have a thick elastic lithosphere, between 80 and 160 km [34,39]. However, a thick lithosphere is not a viable option for data on tilted palaeo shorelines of coastal southwestern Norway and not for the Barents Sea [34,52]. However, a recent study of the Tapes transgression along the Norwegian coast indicates a thicker effective elastic lithosphere thickness of 70 km (flexural rigidity 3 × 1024 Nm) in Northern Norway [53].
The glacial isostatic tilt caused by a 70 km thick elastic lithosphere is shown in Figure 8. Sedimentary basins located in SW Barents Sea could be tilted by up to 2.0 m/km due to the LGM ice sheet, which is about 70% of the tilt with lithospheric thickness of 30 km (cf. Figure 2). The flexural stress due to a 70 km thick lithosphere is thus somewhat lower than for a thinner lithosphere, but the stresses will be located basically in the same area. We believe, however, that a model with thin elastic lithosphere of 30 km better describes the Earth rheology of the Barents Sea than a model with thicker elastic lithosphere (see [52,54]).

4.3. Stress Models

The way the stress distributes in a material is depending on existing fractures and weaknesses and their orientation to the added load. In the stress model examples we have shown here, a published regional interpreted fault map was used as basis. A more detailed fault interpretation will be needed for a more local study of the stress effects of glacially induced flexure.

4.4. Implications for Basin Modeling

Repeated ice loading during the Quaternary could have significant impact on sedimentary basins and add a degree of difficulty to petroleum exploration. The isostatic effects of ice sheets, (in addition to possible glacial erosion/deposition) could give significant tilting of sedimentary basins in the peripheral areas of former glaciated area. The tilting of the reservoirs at the Norwegian Continental Shelf could be up to almost 3 m/km and could alter hydrocarbon migration pathways. The ice sheets will act as seesaw during glacial/interglacial period, while the isostatic response of erosion and sedimentation will grow iteratively during the glaciations, depending on the stability of the depocenters.
Differential vertical movement of the lithosphere related to glacial isostasy leads to repeated tilting of sedimentary formations and potential petroleum reservoirs therein, which may have greatly affected hydrocarbon migration pathways in the former glaciated area [3]. In addition, the upward and downward bending of the lithosphere lead to flexural stresses likely to affect faults and their permeability, which could add to changes in hydrocarbon migration pathways. The flexural stresses are adding to background stresses in the area.
The map of the flexural effect based on glacial isostasy shows that the SW Barents Sea is affected by flexural stresses due to the last ice age. Stress models suggest that this effect may have caused high magnitude stresses to accumulate and therefore lead to an increased fracture-related permeability and escape of hydrocarbons in certain areas.
Most of the production licenses on the Norwegian Continental Shelf are located inside the former Last Glacial Maximum ice sheet (Figure 1). The results of our study demonstrate the importance of including the effect of glaciations in basin modeling studies of areas of petroleum potential. Adding to this are the temperature effects of glaciations and the effects of glacial erosion on temperature and stress discussed in more detail in [3].

Author Contributions

Conceptualization, I.F.L., W.F.; methodology, I.F.L., W.F., M.S.; software, I.F.L., W.F.; investigation, I.F.L., W.F., M.S.; writing—original draft preparation, I.F.L., W.F.; writing—review and editing, I.F.L., W.F.; M.S.; project administration, I.F.L.

Funding

Part of this study was supported through project no. 200657 “Neogene Uplift of the Barents Sea” by Norwegian Research Council and the oil companies ConocoPhillips, Det Norske, E.ON E&P, GDF Suez E&P Norge AS, NORECO, OMV (Norge), Statoil, Total E&P and Wintershall.

Acknowledgments

The authors would like to thank four anonymous reviewers for constructive and fruitful feedback, which improved the paper.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Glacial Isostatic Model

The ultimate isostatic compensation achieved by any harmonic component is given by [44,55]:
h ( k )   =   F ( k ) α − 1 ρ g
where F(k) = transformed ice load, ρ = density of the upper mantle, and g = gravity.
The ’lithosphere filter’ is α = 2 μ k ρ g ( S 2 −   k 2 H 2 )   + ( CS +   kH ) ] / ( S +   kHC ) , where k = wavenumber, H = elastic thickness of the lithosphere, μ = Lame’s parameter, S = sinh kH and C = cosh kH. The relation between the flexural rigidity (D) and the elastic thickness of the lithosphere (H) is given by the equation D = E H 3 12 ( 1 − ν 2 ) , where E = Young’s modulus and ν = Poisson’s ratio.

References

  1. Gibbard, P.L.; Head, M.J.; Walker, M.J.C. Formal ratification of the quaternary system/period and the pleistocene series/Epoch with a base at 2.58 Ma. J. Quat. Sci. 2009, 25, 96–102. [Google Scholar] [CrossRef]
  2. Mangerud, J.; Gyllencreutz, R.; Lohne, Ø.; Svendsen, J.I. Glacial history of Norway. In Developments in Quaternary Science; Ehlers, J., Gibbard, P.L., Hughes, P.D., Eds.; Elsevier: Amsterdam, The Netherlands, 2011; Volume 15, pp. 279–298. ISBN 978-0-444-53447-7. [Google Scholar]
  3. Fjeldskaar, W.; Amantov, A. Effects of glaciations on sedimentary basins. J. Geodyn. 2018, 118, 66–81. [Google Scholar] [CrossRef]
  4. Argus, D.F.; Drummond, R.; Peltier, W.R. Space geodesy constrains ice age terminal deglaciation: The global ICE-6G_C (VM5a) model. J. Geophys. Res. Solid Earth 2015, 120, 450–487. [Google Scholar] [CrossRef]
  5. Kjemperud, A.; Fjeldskaar, W. Pleistocene glacial isostasy-implications for petroleum geology. In Structural and Tectonic Modelling and its Application to Petroleum Geology; Larsen, R.M., Brekke, H., Larsen, B.T., Talleraas, E., Eds.; Elsevier: Amsterdam, The Netherlands, 1992; pp. 187–195. [Google Scholar]
  6. Grunnaleite, I.; Fjeldskaar, W.; Wilson, J.; Faleide, J.; Zweigel, J. Effect of local variations of vertical and horizontal stresses on the Cenozoic structuring of the mid-Norwegian shelf. Tectonophysics 2009, 470, 267–283. [Google Scholar] [CrossRef]
  7. Vadakkepuliyambatta, S.; Bünz, S.; Mienert, J.; Chand, S. Distribution of subsurface fluid-flow systems in the SW Barents Sea. Mar. Pet. Geol. 2013, 43, 208–221. [Google Scholar] [CrossRef] [Green Version]
  8. Cavanagh, A.J.; Di Primio, R.; Scheck-Wenderoth, M.; Horsfield, B. Severity and timing of Cenozoic exhumation in the southwestern Barents Sea. J. Geol. Soc. 2006, 163, 761–774. [Google Scholar] [CrossRef]
  9. Doré, A.G.; Jensen, L.N. The impact of late Cenozoic uplift and erosion on hydrocarbon exploration: offshore Norway and some other uplifted basins. Glob. Planet. Change 1996, 12, 415–436. [Google Scholar] [CrossRef]
  10. Duran, E.R.; Di Primio, R.; Anka, Z.; Stoddart, D.; Horsfield, B. 3D-basin modelling of the Hammerfest Basin (southwestern Barents Sea): A quantitative assessment of petroleum generation, migration and leakage. Mar. Pet. Geol. 2013, 45, 281–303. [Google Scholar] [CrossRef]
  11. Zieba, K.J.; Grøver, A. Isostatic response to glacial erosion, deposition and ice loading. Impact on hydrocarbon traps of the southwestern Barents Sea. Mar. Pet. Geol. 2016, 78, 168–183. [Google Scholar] [CrossRef]
  12. Lerche, I.; Yu, Z.; Tørudbakken, B.; Thomsen, R. Ice loading effects in sedimentary basins with reference to the Barents sea. Mar. Pet. Geol. 1997, 14, 277–338. [Google Scholar] [CrossRef]
  13. Stoddart, D.; Jørstad, A.; Rønnevik, H.-C.; Fjeldskaar, W.; Løtveit, I.F. Recent glacial events in the Norwegian North Sea-implications towards a better understanding of charging/leakage of oil fields and its impact oil exploration. In Proceedings of the IMOG 2015—27th International Meeting on Organic Geochemistry, Prague, Czech Republic, 13–18 September 2015. [Google Scholar] [CrossRef]
  14. Medvedev, S.; Hartz, E.H.; Schmid, D.W.; Zakariassen, E.; Varhaug, P. Influence of Glaciations on North Sea Petroleum Systems; Geological Society: London, UK, Special Publications; 2019; Volume 494, p. SP494-2018-2183. [Google Scholar] [CrossRef]
  15. Mikko, H.; Smith, C.A.; Lund, B.; Ask, M.V.; Munier, R. LiDAR-derived inventory of post-glacial fault scarps in Sweden. GFF 2015, 137, 334–338. [Google Scholar] [CrossRef]
  16. Dehls, J.F.; Olesen, O.; Olsen, L.; Blikra, L.H. Neotectonic faulting in northern Norway; the Stuoragurra and Nordmannvikdalen postglacial faults. Quat. Sci. Rev. 2000, 19, 1447–1460. [Google Scholar] [CrossRef]
  17. Lagerbäck, R. Neotectonic structures in northern Sweden. Geol. Foren. Stockh. Förh. 1979, 100, 271–278. [Google Scholar] [CrossRef]
  18. Brandes, C.; Winsemann, J.; Roskosch, J.; Meinsen, J.; Tanner, D.C.; Frechen, M.; Steffen, H.; Wu, P. Activity along the Osning Thrust in Central Europe during the Lateglacial: Ice-sheet and lithosphere interactions. Quat. Sci. Rev. 2012, 38, 49–62. [Google Scholar] [CrossRef]
  19. Sandersen, P.B.E.; Jørgensen, F. Neotectonic deformation of a Late Weichselian outwash plain by deglaciation-induced fault reactivation of a deep-seated graben structure. Boreas 2015, 44, 413–431. [Google Scholar] [CrossRef]
  20. Grube, A. Palaeoseismic structures in Quaternary sediments, related to an assumed fault zone north of the Permian Peissen-Gnutz salt structure (NW Germany)—Neotectonic activity and earthquakes from the Saalian to the Holocene. Geomorphology 2019, 328, 15–27. [Google Scholar] [CrossRef]
  21. Hoffmann, G.; Reicherter, K. Soft-sediment deformation of late Pleistocene sediments along the southwestern coast of the Baltic Sea. Intern. J. Earth Sci. 2012, 101, 351–363. [Google Scholar] [CrossRef]
  22. Van Loon, A. (Tom); Pisarska-Jamroży, M. Sedimentological evidence of Pleistocene earthquakes in NW Poland induced by glacio-isostatic rebound. Sediment. Geol. 2014, 300, 1–10. [Google Scholar] [CrossRef]
  23. Stewart, I.S.; Firth, C.R.; Rust, D.J.; Collins, P.E.; Firth, J.A. Postglacial fault movement and palaeoseismicity in western Scotland: A reappraisal of the Kinloch Hourn fault, Kintail. J. Seism. 2001, 5, 307–328. [Google Scholar] [CrossRef]
  24. Fejerskov, M.; Lindholm, C. Crustal stress in and around Norway: an evaluation of stress-generating mechanisms. Geol. Soc. London Spéc. Publ. 2000, 167, 451–467. [Google Scholar] [CrossRef]
  25. Druzhinina, O.; Bitinas, A.; Molodkov, A.; Kolesnik, T. Palaeoseismic deformations in the Eastern Baltic region (Kaliningrad District of Russia). Estonian J. Earth Sci. 2017, 66, 119–129. [Google Scholar] [CrossRef]
  26. Antonovskaya, G.; Konechnaya, Y.; Kremenetskaya, E.O.; Asming, V.; Kvaerna, T.; Schweitzer, J.; Ringdal, F. Enhanced earthquake monitoring in the European Arctic. Polar Sci. 2015, 9, 158–167. [Google Scholar] [CrossRef]
  27. Gibbons, S.J.; Antonovskaya, G.; Asming, V.; Konechnaya, Y.V.; Kremenetskaya, E.; Kværna, T.; Schweitzer, J.; Vaganova, N.V. The 11 October 2010 Novaya Zemlya Earthquake: Implications for Velocity Models and Regional Event Location. Bull. Seismol. Soc. Am. 2016, 106, 1470–1481. [Google Scholar] [CrossRef] [Green Version]
  28. Neuzil, C.E. Hydromechanical effects of continental glaciation on groundwater systems. Geofluids 2012, 12, 22–37. [Google Scholar] [CrossRef]
  29. Steffen, R.; Steffen, H.; Wu, P.; Eaton, D.W. Reply to comment by Hampel et al. on “Stress and fault parameters affecting fault slip magnitude and activation time during a glacial cycle.”. Tectonics 2015, 34, 2359–2366. [Google Scholar]
  30. Grollimund, B.; Zoback, M.D. Post glacial lithospheric flexure and induced stresses and pore pressure changes in the northern North Sea. Tectonophysics 2000, 327, 61–81. [Google Scholar] [CrossRef]
  31. Johnston, A.C. Suppression of earthquakes by large continental ice sheets. Nature 1987, 330, 467–469. [Google Scholar] [CrossRef]
  32. Turpeinen, H.; Hampel, A.; Karow, T.; Maniatis, G. Effect of ice sheet growth and melting on the slip evolution of thrust faults. Earth Planet. Sci. Lett. 2008, 269, 230–241. [Google Scholar] [CrossRef]
  33. Steffen, R.; Steffen, H.; Wu, P.; Eaton, D.W. Stress and fault parameters affecting fault slip magnitude and activation time during a glacial cycle. Tectonics 2014, 33, 1461–1476. [Google Scholar] [CrossRef] [Green Version]
  34. Fjeldskaar, W.; Amantov, A. Tilted Norwegian post-glacial shorelines require a low viscosity asthenosphere and a weak lithosphere. Reg. Geol. Metallogeny 2017, 70, 48–59. [Google Scholar]
  35. Larsen, E.; Fredin, O.; Lyså, A.; Amantov, A.; Fjeldskaar, W.; Ottesen, D. Causes of time-transgressive glacial maxima positions of the last Scandinavian Ice Sheet. Nor. J. Geol. 2016, 96. [Google Scholar]
  36. Hughes, A.L.C.; Gyllencreutz, R.; Lohne, Ø.S.; Mangerud, J.; Svendsen, J.I. The last Eurasian ice sheets—A chronological database and time-slice reconstruction, DATED-1. Boreas 2016, 45, 1–45. [Google Scholar] [CrossRef]
  37. Zweck, C.; Huybrechts, P. Modeling of the northern hemisphere ice sheets during the last glacial cycle and glaciological sensitivity. J. Geophys. Res. Space Phys. 2005, 110, 1984–2012. [Google Scholar] [CrossRef]
  38. Auriac, A.; Whitehouse, P.L.; Bentley, M.J. Glacial isostatic adjustment associated with the Barents Sea ice sheet: A modeling inter-comparison. Quat. Sci. Rev. 2016, 147, 122–135. [Google Scholar] [CrossRef]
  39. Steffen, H.; Wu, P. Glacial isostatic adjustment in Fennoscandia—A review of data and modeling. J. Geodyn. 2011, 52, 169–204. [Google Scholar] [CrossRef]
  40. Amantov, A.; Fjeldskaar, W. Meso-Cenozoic exhumation and relevant isostatic process: The Barents and Kara shelves. J. Geodyn. 2018, 118, 118–139. [Google Scholar] [CrossRef]
  41. Fjeldskaar, W. The flexural rigidity of Fennoscandia inferred from the post-glacial uplift. Tectonics 1997, 16, 596–608. [Google Scholar] [CrossRef]
  42. Stephansson, O. Stress Measurements And Modelling of Crustal Rock Mechanics in Fennoscandia; In Earthquakes at North-Atlantic Passive Margins: Neotectonics and Postglacial, Rebound, Gregersen, S., Basham, P.W., Eds.; Springer: Dordrecht, The Netherlands, 1989; pp. 213–229. [Google Scholar]
  43. Stein, S.; Cloetingh, S.; Sleep, N.H.; Wortel, R.; Gregersen, S.; Basham, P.W. Passive Margin Earthquakes, Stresses and Rheology. In Earthquakes at North-Atlantic Passive Margins: Neotectonics and Postglacial Rebound; Springer: Dordrecht, The Netherlands, 1989; pp. 231–259. [Google Scholar] [Green Version]
  44. Cathles, L.M. The Viscosity of the Earth´s Mantle; Princeton University Press: Princeton, NJ, USA, 1975; p. 386. [Google Scholar]
  45. Brekke, H.; Kalheim, J.E.; Riis, F.; Egeland, B.; Blystad, P.; Johnsen, S.; Ragnhildstveit, J. Two-Way Time-Map of the Unconformity at the Base of Top Jurassic (North of 69° N) and the Unconformity at the Base of Cretaceous (South of 69°N), Offshore Norway, Including the Main Geological Trends Onshore. NPD Continental Shelf Map No. 1; The Norwegian Petroleum Directorate/The Geological Survey of Norway: Stavanger/Trondheim, Norway, 1992. [Google Scholar]
  46. Gudmundsson, A. Rock Fractures in Geological Processes, 1st ed.; Cambridge University Press: New York, NY, USA, 2011. [Google Scholar]
  47. Gudmundsson, A.; Simmenes, T.H.; Larsen, B.; Philipp, S.L. Effects of internal structure and local stresses on fracture propagation, deflection, and arrest in fault zones. J. Struct. Geol. 2009, 32, 1643–1655. [Google Scholar] [CrossRef]
  48. Paterson, M.S.; Paterson, D.M.S. Experimental Rock Deformation—The Brittle Field; Springer: Berlin, Germany, 1978; Volume 13. [Google Scholar]
  49. Farmer, I. Engineering Behaviour of Rocks; Chapman and Hal: London, UK, 1983. [Google Scholar]
  50. Amadei, B.; Stephansson, O. Rock Stress and its Measurement; Chapman and Hall: London, UK, 1997; p. 490. [Google Scholar]
  51. Jaeger, J.C.; Cook, N.G.W.; Zimmerman, R.W. Fundamentals of Rock Mechanics, 2nd ed.; Chapman and Hall: London, UK, 2007. [Google Scholar]
  52. Fjeldskaar, W.; Bondevik, S.; Amantov, A. Glaciers on Svalbard survived the Holocene thermal optimum. Quat. Sci. Rev. 2018, 199, 18–29. [Google Scholar] [CrossRef]
  53. Fjeldskaar, W.; Bondevik, S. The Early-Mid-Holocene transgression (Tapes) at the Norwegian coast–comparing observations with numerical modelling. Quat. Sci. Rev. (submitted).
  54. Medvedev, S.; Hartz, E.H.; Faleide, J.I. Erosion-driven vertical motions of the circum Arctic: Comparative analysis of modern topography. J. Geodyn. 2018, 119, 62–81. [Google Scholar] [CrossRef]
  55. Fjeldskaar, W.; Cathles, L. Rheology of Mantle and Lithosphere Inferred from Post-Glacial Uplift in Fennoscandia. In Glacial Isostasy, Sea-Level and Mantle Rheology; Springer: Berlin, Germany, 1991; pp. 1–19. [Google Scholar]
Geosciences 09 00474 g001 550
Figure 1. Extent (based on [36]) and maximum thickness model for the Last Glacial Maximum (LGM), from [3]. The production licenses (from www.npd.no) in the Norwegian continental shelf are shown by grey polygons.
Figure 1. Extent (based on [36]) and maximum thickness model for the Last Glacial Maximum (LGM), from [3]. The production licenses (from www.npd.no) in the Norwegian continental shelf are shown by grey polygons.
Geosciences 09 00474 g001
Geosciences 09 00474 g002a 550Geosciences 09 00474 g002b 550
Figure 2. Calculated tilts (m/km) due to glacial isostasy for the LGM ice sheet (based on [3]). The white box in a) marks the Norwegian part of the Barents Sea that is open for petroleum exploration (production licenses are marked in grey). The tilts of this area are shown in b).
Figure 2. Calculated tilts (m/km) due to glacial isostasy for the LGM ice sheet (based on [3]). The white box in a) marks the Norwegian part of the Barents Sea that is open for petroleum exploration (production licenses are marked in grey). The tilts of this area are shown in b).
Geosciences 09 00474 g002aGeosciences 09 00474 g002b
Geosciences 09 00474 g003 550
Figure 3. Illustration of tilting during one glacial cycle (redrawn from [5]). Upper left shows the situation prior to tilting; upper right during tilting. Lower panel shows the situation after deglaciation. “Spill point” is the structurally lowest point in a hydrocarbon trap that can retain hydrocarbons. Once a trap has been filled to its spill point, further storage or retention of hydrocarbons will not occur. The illustration assumes the trap filled to the spill point. HCWC is the hydrocarbon-water contact.
Figure 3. Illustration of tilting during one glacial cycle (redrawn from [5]). Upper left shows the situation prior to tilting; upper right during tilting. Lower panel shows the situation after deglaciation. “Spill point” is the structurally lowest point in a hydrocarbon trap that can retain hydrocarbons. Once a trap has been filled to its spill point, further storage or retention of hydrocarbons will not occur. The illustration assumes the trap filled to the spill point. HCWC is the hydrocarbon-water contact.
Geosciences 09 00474 g003
Geosciences 09 00474 g004 550
Figure 4. Left: Flexural stress model ([42]) assuming crustal stress equilibrium prior to the ice load. The resulting stresses are large during loading and diminish during unloading. The resulting stress state is compressional beneath the ice sheet and extensional in the hinge zones. Right: Flexural stress model ([43]) assuming crustal stress equilibrium after a period of ice loading. From [28].
Figure 4. Left: Flexural stress model ([42]) assuming crustal stress equilibrium prior to the ice load. The resulting stresses are large during loading and diminish during unloading. The resulting stress state is compressional beneath the ice sheet and extensional in the hinge zones. Right: Flexural stress model ([43]) assuming crustal stress equilibrium after a period of ice loading. From [28].
Geosciences 09 00474 g004
Geosciences 09 00474 g005 550
Figure 5. Derivative of the glacial isostasy used to identify areas of high flexural stresses.
Figure 5. Derivative of the glacial isostasy used to identify areas of high flexural stresses.
Geosciences 09 00474 g005
Geosciences 09 00474 g006 550
Figure 6. Second derivative, or curvature (m/km2), of the isostatic response of the LGM ice sheet (Figure 1) in SW Barents Sea license area (white square of Figure 2). Yellow, orange, and green areas have high curvature and are likely to be influenced by stress effects during and after the glaciation. (cf. Figure 5). The black square marks the area included in the stress modeling, and the blue arrows show the orientation of maximum flexure-related extension after LGM deglaciation within that area.
Figure 6. Second derivative, or curvature (m/km2), of the isostatic response of the LGM ice sheet (Figure 1) in SW Barents Sea license area (white square of Figure 2). Yellow, orange, and green areas have high curvature and are likely to be influenced by stress effects during and after the glaciation. (cf. Figure 5). The black square marks the area included in the stress modeling, and the blue arrows show the orientation of maximum flexure-related extension after LGM deglaciation within that area.
Geosciences 09 00474 g006
Geosciences 09 00474 g007 550
Figure 7. Stress models covering an area of 120 × 90 km in the Hoop Fault Complex, SW Barents Sea. The results include related exploration wells with data made public prior to 2017. Top: Resulting distribution of maximum principal tensile stress (σ3) during SW–NE extension due to deglaciation of LGM. Bottom: Resulting distribution of von Mises shear stress during flexure-related SW–NE extension. Scale in megapascals.
Figure 7. Stress models covering an area of 120 × 90 km in the Hoop Fault Complex, SW Barents Sea. The results include related exploration wells with data made public prior to 2017. Top: Resulting distribution of maximum principal tensile stress (σ3) during SW–NE extension due to deglaciation of LGM. Bottom: Resulting distribution of von Mises shear stress during flexure-related SW–NE extension. Scale in megapascals.
Geosciences 09 00474 g007
Geosciences 09 00474 g008 550
Figure 8. The isostatic tilts for lithosphere thickness of 70 km.
Figure 8. The isostatic tilts for lithosphere thickness of 70 km.
Geosciences 09 00474 g008

Share and Cite

MDPI and ACS Style

Løtveit, I.F.; Fjeldskaar, W.; Sydnes, M. Tilting and Flexural Stresses in Basins Due to Glaciations—An Example from the Barents Sea. Geosciences 2019, 9, 474. https://doi.org/10.3390/geosciences9110474

AMA Style

Løtveit IF, Fjeldskaar W, Sydnes M. Tilting and Flexural Stresses in Basins Due to Glaciations—An Example from the Barents Sea. Geosciences. 2019; 9(11):474. https://doi.org/10.3390/geosciences9110474

Chicago/Turabian Style

Løtveit, Ingrid F., Willy Fjeldskaar, and Magnhild Sydnes. 2019. "Tilting and Flexural Stresses in Basins Due to Glaciations—An Example from the Barents Sea" Geosciences 9, no. 11: 474. https://doi.org/10.3390/geosciences9110474

APA Style

Løtveit, I. F., Fjeldskaar, W., & Sydnes, M. (2019). Tilting and Flexural Stresses in Basins Due to Glaciations—An Example from the Barents Sea. Geosciences, 9(11), 474. https://doi.org/10.3390/geosciences9110474

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop