Deformation of salt structures by ice-sheet loading: insights into the controlling parameters from numerical modelling

Abstract

Subsurface salt flow is driven by differential loading, which is typically caused by tectonics or sedimentation. During glaciations, the weight of an ice sheet represents another source of differential loading. In salt-bearing basins affected by Pleistocene glaciations, such as the Central European Basin System, ice loading has been postulated as a trigger of young deformation at salt structures. Here, we present finite-element simulations (ABAQUS) with models based on a simplified 50-km long and 10-km-deep two-dimensional geological cross-section of a salt diapir subject to the load of a 300-m-thick ice sheet. The focus of our study is to evaluate the sensitivity of the model to material parameters, including linear and non-linear viscosity of the salt rocks and different elasticities. A spatially and temporarily variable pressure was applied to simulate ice loading. An ice advance towards the diapir causes lateral salt flow into the diapir and diapiric rise. Complete ice coverage leads to downward displacement of the diapir. After unloading, displacements are largely restored. The modelled displacements do not exceed few metres and are always larger in models with linear viscosity than in those with non-linear viscosity. Considering the low stresses caused by ice-sheet loading and the long time-scale, the application of linear viscosity seems appropriate. The elastic parameters also have a strong impact, with lower Young's moduli leading to larger deformation. The impact of both the viscosity and the elasticity highlights the importance of a careful parameter choice in numerical modelling, especially when aiming to replicate any real-world observations.