Abstract

Abstract. The 2-D distinct element method (DEM) code (PFC2D_V5) is used here to simulate the evolution of subsidence-related karst landforms, such as single and clustered sinkholes, and associated larger-scale depressions. Subsurface material in the DEM model is removed progressively to produce an array of cavities; this simulates a network of subsurface groundwater conduits growing by chemical/mechanical erosion. The growth of the cavity array is coupled mechanically to the gravitationally loaded surroundings, such that cavities can grow also in part by material failure at their margins, which in the limit can produce individual collapse sinkholes. Two end-member growth scenarios of the cavity array and their impact on surface subsidence were examined in the models: (1) cavity growth at the same depth level and growth rate; (2) cavity growth at progressively deepening levels with varying growth rates. These growth scenarios are characterised by differing stress patterns across the cavity array and its overburden, which are in turn an important factor for the formation of sinkholes and uvala-like depressions. For growth scenario (1), a stable compression arch is established around the entire cavity array, hindering sinkhole collapse into individual cavities and favouring block-wise, relatively even subsidence across the whole cavity array. In contrast, for growth scenario (2), the stress system is more heterogeneous, such that local stress concentrations exist around individual cavities, leading to stress interactions and local wall/overburden fractures. Consequently, sinkhole collapses occur in individual cavities, which results in uneven, differential subsidence within a larger-scale depression. Depending on material properties of the cavity-hosting material and the overburden, the larger-scale depression forms either by sinkhole coalescence or by widespread subsidence linked geometrically to the entire cavity array. The results from models with growth scenario (2) are in close agreement with surface morphological and subsurface geophysical observations from an evaporite karst area on the eastern shore of the Dead Sea.

Highlights

  • Karstification occurs worldwide in rocks like limestone, dolomite, gypsum, anhydrite and salt primarily by chemical dissolution (BGR et al, 2017)

  • While subsurface-solutionbased drainage networks and connected void spaces resulting from karstification are hydrologically important for groundwater provision (Chen et al, 2017), such features reduce the mechanical stability of the geologic material and so may Published by Copernicus Publications on behalf of the European Geosciences Union

  • We examine two end-member growth scenarios of model cavity arrays, and we look at the surface morphologies, subsurface structure and stress patterns developed by subsidence of the overburden as those cavity arrays grow

Read more

Summary

Introduction

Karstification occurs worldwide in rocks like limestone, dolomite, gypsum, anhydrite and salt primarily by chemical dissolution (BGR et al, 2017). Systems develop into agglomerations of closely spaced or coalesced dolines and elongated valley-like depressions, potentially revealing linear patterns of drainage (Waltham et al, 2005) Such sinkhole cluster development can be highly dynamic and partly accelerating, and may affect large areas in short times Such clusters commonly lie within gentler, larger-scale (uvala-like) depressions of up to several hundreds of metres in diameter, as depicted in Fig. 1b and c These karst landforms develop as small localised subsidence zones, with single sinkholes that form in heterogeneous material made of Dead Sea mud, alluvial fan sediments and salt (Watson et al, 2019). We show that our more complex end-member modelling scenario is able to explain complementary observations from surface morphology to subsurface hydrology and subsurface geophysics

Distinct element method numerical modelling
Cavity growth in a DEM model
Material parameters by
Geophysical parameter tracking
Modelling results
Stresses and strains in a multiple void space system
Generic geophysical parameters
Discussion
Subsurface patterns of sinkhole clusters and subrosion
Summary and conclusions
Cavity growth implementation
Optimal model development
Principal stresses in a multiple void space system
Findings
Incremental shear strain evolution
Maximum shear stress evolution
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call