Large grain-size-dependent rheology contrasts of halite at low differential stress: evidence from microstructural study of naturally deformed gneissic Zechstein 2 rock salt (Kristallbrockensalz) from the northern Netherlands

Jessica Barabasch,Joyce Schmatz,Jop Klaver,Alexander Schwedt,and Janos L. Urai


Constitutive laws to predict long-term deformation of solution-mined caverns and radioactive-waste repositories in rock salt play an important role in the energy transition. Much of this deformation is at differential stresses of a few megapascals, while the vast majority of laboratory measurements are at much higher differential stress and require extrapolation. This can be much improved by including microstructural data of samples deformed in natural laboratories. Deformation of rock salt can occur by dislocation creep and grain-size-dependent dissolution–precipitation creep processes (pressure solution); this mechanism is not commonly included in current engineering predictions.

Here we show evidence for large grain-size-dependent differences in rock salt rheology based on microstructural observations from Zechstein rock salt cores of the northern Netherlands that experienced different degrees of tectonic deformation. We studied the relatively undeformed horizontally layered Zechstein 2 (Z2) salt (Stassfurt Formation) from Barradeel and compared it with a much more strongly deformed equivalent in diapiric salt from Winschoten, Zuidwending, and Pieterburen. We used optical microscopy of thin gamma-irradiated sections for microtectonic analysis, recrystallized grain-size measurements and subgrain-size piezometry, electron microscopy with energy-dispersive X-ray spectroscopy, and X-ray diffraction analysis for second-phase mineralogy. Subgrain-size piezometry shows that this deformation took place at differential stresses between 0.5 and 2 MPa.

In the undeformed, layered salt from Barradeel we find centimetre-thick layers of single crystalline halite (Kristalllagen or megacrystals) alternating with fine-grained halite and thin anhydrite layers. The domal salt samples are typical of the well-known “Kristallbrockensalz” and consist of centimetre-size tectonically disrupted megacrystals surrounded by fine-grained halite with a grain size of a few millimetres. We infer high strains in the fine-grained halite as shown by folding and boudinage of thin anhydrite layers, as compared to the megacrystals, which are internally much less deformed and develop subgrains during dislocation creep. Subgrain size shows comparable differential stresses in Kristallbrocken as in matrix salt. The fine-grained matrix salt is dynamically recrystallized to some extent and has few subgrains and microstructures, indicating deformation by solution–precipitation processes. We infer that the finer-grained halite deformed dominantly via pressure solution and the megacrystals dominantly by dislocation creep.

The samples show that the fine-grained matrix salt is much weaker than Kristallbrocken because of different dominant deformation mechanisms. This is in agreement with microphysical models of pressure solution creep in which grain size has a significant effect on strain rate at low differential stress. Our results point to the importance of pressure solution creep in rock salt at low differential stresses around engineered structures but also in most salt tectonic settings. We suggest that including results of microstructural analysis can strongly improve engineering models of rock salt deformation.

We recommend that this mechanism of grain-size-dependent rheology is included more consistently in the constitutive laws describing the deformation of rock salt.

Please read the corrigendum first before continuing.

How to cite. 

Barabasch, J., Schmatz, J., Klaver, J., Schwedt, A., and Urai, J. L.: Large grain-size-dependent rheology contrasts of halite at low differential stress: evidence from microstructural study of naturally deformed gneissic Zechstein 2 rock salt (Kristallbrockensalz) from the northern Netherlands, Solid Earth, 14, 271–291,, 2023.Received: 15 Jul 2022

 – Discussion started: 21 Jul 2022

 – Revised: 29 Dec 2022

 – Accepted: 17 Jan 2023

 – Published: 09 Mar 2023

1 Introduction

1.1 Salt rheology, deformation processes, and associated microstructures – relevance of natural laboratories

Salt – unlike most sedimentary rocks – is very weak and deforms like a viscous fluid, even at low temperatures and shallow depths, and its presence in sedimentary basins fundamentally changes their evolution by salt tectonic processes (Jackson et al., 1994; Urai et al., 2008; Jackson and Hudec, 2017). This has long been recognized for salt diapirs (Trusheim, 1957; Schultz-Ela et al., 1993; Kukla et al., 2019), salt glaciers (Wenkert, 1979; Talbot, 1979; Talbot and Aftabi, 2004; Závada et al., 2012), salt decoupling (Peel, 2014; Rowan and Krzywiec, 2014; Tămaṡ et al., 2021), and at a smaller scale for internal salt flow creating folding or boudinage (van Gent et al., 2011; Fiduk and Rowan, 2012; Strozyk et al., 2014; Dooley et al., 2015; Rowan et al., 2019; Adamuszek et al., 2021).

Rock salt is the host of many large engineering constructions like salt caverns and radioactive-waste repositories, and the long-term operation and abandonment of these structures requires predictions of deformation up to thousands of years (Brouard et al., 2004; Chemia et al., 2009; Bräuer et al., 2011; Baumann et al., 2018). This deformation takes place at differential stresses of below a few megapascals; however the measurement of salt rheology in the laboratory is problematic because of the long duration of the required experiments (Brouard et al., 2004; Bérest et al., 2010). Because laboratory experiments are limited in time, they are commonly conducted at relatively high deviatoric stresses and strain rates at which dislocation creep processes are expected and inferred to be the dominant deformation mechanism (Wawersik and Zeuch, 1986; Carter et al., 1993; Hunsche and Hampel, 1999; Hampel et al., 2007; Urai and Spiers, 2007; Urai et al., 2008; Hampel, 2016). Therefore, most of the long-term predictions of salt deformation have been made based on extrapolation of such experimentally derived power law creep rheology to low differential stress (Albrecht et al., 1993; Bräuer et al., 2011), although some predictions include both mechanisms (Buijze et al., 2022). However, it has been documented in the geological and materials science literature that below a certain differential stress in rock salt there can be a dramatic grain-size-dependent change in deformation mechanism from dislocation creep to pressure solution creep, and this can change the creep rate by many orders of magnitude (Urai et al., 1986b; Spiers et al., 1990; van Keken et al., 1993; Spiers and Carter, 1998).

In recent years, the validity of dislocation creep-based extrapolations was extensively discussed in the engineering community, and very slow creep tests were carried out that indeed show much faster deformation than power law creep rates would predict (Popp and Hansen, 2018; Bérest et al., 2019).

An important contribution to this discussion comes from materials science, where microphysical processes are related to constitutive laws by state variables such as water content (at grain boundaries), solid solution, and second-phase impurities and grain size (Carter and Hansen, 1983; Heard and Ryerson, 1986; Urai et al., 1986b; Hunsche et al., 2003; Bérest et al., 2019; Spiers et al., 1990). Recent experiments show that the presence of secondary minerals such as phyllosilicates adds further complexity in the dissolution–precipitation process with potential influence on pressure solution creep rates (Macente et al., 2018; Schwichtenberg et al., 2022).

For steady-state and non-dilatant deformation the microphysics-based deformation rate of rock salt is given by Eq. (1)

for dislocation creep and Eq. (2)

for solution–precipitation creep, with the sum of both being the total strain rate (Eq. 3):

(Wawersik and Zeuch, 1986; Spiers et al., 1990; Carter et al., 1993; van Keken et al., 1993; Hunsche and Hampel, 1999; Urai et al., 2008). Thereby A and B are material constants, QDC and QPS are apparent activation energies for dislocation creep and pressure solution creep, R is the gas constant, T is absolute temperature, σ1σ3 is the differential stress, DG is the grain size, and n and m are stress and grain size exponents. For a compilation of parameters for the two mechanisms, we refer to Urai et al. (2008) and references therein (e.g. Spiers et al., 1990; van Keken et al., 1993; Brouard and Bérest, 1998; Hunsche and Hampel, 1999; Hunsche et al., 2003).

Because of the large differences in the stress and grain size exponents of the two mechanisms, pressure solution creep will dominate deformation at low differential stress and small grain size, and dislocation creep will dominate deformation at high differential stress and large grain size. Key parameters in the evolution of grain size are the structure and mobility of brine-filled grain boundaries, which can change as a function of differential stress (Drury and Urai, 1990; van Noort et al., 2008), and the migration of grain boundaries, which can change the grain size (Urai et al., 1986a; Peach et al., 2001; Schenk and Urai, 2004; Schenk et al., 2006; Schmatz et al., 2011).

Because of the extremely long experiments (several years) required to measure pressure solution creep rates in rock salt at low differential stress (Bérest et al., 2019) and the very small strains reached in these experiments, it is useful to compare these results with naturally deformed samples. The differential stresses during natural deformation can be measured using subgrain-size piezometry and are commonly found to be between 0.5 and 5 MPa (Carter et al., 1982; Schléder and Urai, 2005; Leitner et al., 2011; Rowan et al., 2019), but at much higher strains than can be achieved in low-stress laboratory experiments. In addition, relative differences in rheology can be measured in layered salt (Talbot, 1979; Carter et al., 1982; Schmalholz and Podladchikov, 2001; Schléder and Urai, 2005; Hudleston and Treagus, 2010; Leitner et al., 2011; Komoróczi et al., 2013; Zulauf et al., 2014; Závada et al., 2015; Schmalholz and Mancktelow, 2016; Adamuszek et al., 2021; Schléder and Urai, 2007). Microstructural observations of naturally deformed samples show that dislocation creep, dynamic recrystallization, and pressure solution creep can be readily distinguished, indicated by plastic crystal deformation forming slip bands and subgrains (dislocation creep), grain boundary migration recrystallization (dynamic recrystallization), and growth bands and crystal elongation (pressure solution creep) (Desbois et al., 2010, 2012; Závada et al., 2012, 2015).

1.2 Kristallbrockensalz

The Kristallbrockensalz (z2HS3) in the Zechstein 2 (Stassfurt Formation in Germany) is a well-known tectonite in the Permian Basin (e.g. Czapowski, 1987; Popp et al., 1999; Urai et al., 2019). Richter-Bernburg (1953) described Trümmersalze (German Trümmersalze: debris salt) with Kristallbrocken (German Kristallbrocken: crystal fragments) as large, sharp-edged crystal porphyroclasts inside a fine recrystallized matrix (Küster et al., 2008). Vintage mine photographs from northern Germany show layer-parallel boudins of Kristallbrocken, large rotations of fragments, and regular boudinage in multiple 20 cm thick layers (Richter-Bernburg, 1953; Simon, 1972). It is interpreted to have formed by sedimentation of fine-grained matrix salt, layers of single crystalline halite (German Kristalllagen: crystal layer; Seidl, 1914), and thin anhydrite layers. Although the genesis of extraordinarily large Kristalllagen is not clear, different models have been proposed that suggest post-sedimentary diagenetic grain growth or coalescence of fine-grained halite (Küster et al., 2011). Deformation of the formation resulted in an augen gneiss structure, with the Kristallbrocken decreasing in size with tectonic strain (Küster et al., 2008, 2009; Urai et al., 2019). Pape et al. (2002) presented samples from the Gorleben salt dome and showed that the depositional setting of Kristalllagen halite was shallow marine based on the occurrence of abundant fluid inclusion trails in chevron patterns. Another argument for Kristallbrocken being fragments from a primary layering is based on observations of internal laminations of fluid and solid inclusions (Richter-Bernburg, 1953; Simon, 1972; Pape et al., 2002). Following earlier interpretations (Richter-Bernburg, 1953; Simon, 1972), Pape et al. (2002) hypothesized that rupturing of halite crystal layers is attributable to early diagenetic destruction of layers from subaquatic gliding. This is in contrast to other interpretations (Löffler, 1962) where the so-called Augensalz (another name for Kristallbrockensalz) is interpreted to have formed from rupturing of halite crystal layers during tectonic flow, with clear halite crystallizing in the boudin necks. Extensive studies of distribution of bromide content, Kristallbrocken texture and microstructure, abundant sulfate and fluid inclusions in Kristallbrocken, and interpretation of their formation and deformation mechanisms were presented by Küster et al. (2008, 2009, 2011; Küster, 2011). The structural characteristics of the Kristallbrockensalz were attributed to brittle deformation, dislocation glide, and a strong competence contrast between porphyroclastic Kristallbrocken and fine-grained mylonitic matrix halite (Küster et al., 2008). Küster et al. (2009) described the matrix halite and inferred a secondary and dynamically recrystallized microstructure of sub-structured and sub-structure-free grains. Lobate grain boundaries were interpreted to have formed by grain boundary migration recrystallization, and it was later stated that due to the lack of crystallographic preferred orientation and high water contents, fluid-assisted solution–precipitation creep or grain boundary migration might have been controlling factors for the dominant deformation mechanisms (p. 139 in Küster, 2011). It was hypothesized that second-phase inclusions, large grain size, and the monocrystallinity of Kristallbrocken contribute to the observed rheology contrasts (Küster et al., 2008). In a more recent report from the KEM-17 study (Urai et al., 2019) boudinaged and folded anhydrite layers as well as boudinaged Kristallbrocken surrounded by recrystallized halite and original fine-grained halite were described in Zuidwending and Winschoten samples.

Figure 1(a) Zechstein salt structures in the Permian basin and Zechstein 2 (Z2) carbonate facies distribution after Laier et al. (1998) and Geluk (2000). (b) Study area (indicated in a) after Strozyk et al. (2014) showing top salt depth from seismic interpretation and well locations. BAS: Barradeel; PBN: Pieterburen; WSN: Winschoten; ZWD: Zuidwending. (c) Stratigraphy of Zechstein salt in the Netherlands after Geluk et al. (2007) with stratigraphic position of studied samples.

1.3 Aim of the study

In this study we build on the observations briefly reviewed above; add our new observations of samples with Kristalllagen and Kristallbrocken from the Zechstein of the Netherlands (Fig. 1); and hypothesize that the large competence contrast and different deformation styles between Kristallbrocken and fine-grained matrix salt are caused by the different microphysical deformation mechanisms, with dislocation creep and ductile rupturing (boudinage) being dominant in the Kristallbrocken and grain-size-dependent pressure solution creep being dominant in the fine-grained matrix. If this hypothesis is supported by our data, it provides further evidence for the operation of grain-size-dependent creep in rock salt at the conditions relevant for the operation of engineered structures in rock salt.

2 Materials and methods

We studied 21 samples from drill cores from the areas of Barradeel (BAS), Pieterburen (PBN), Zuidwending (ZWD), and Winschoten (WSN, near Heiligerlee) Zechstein salt structures (Fig. 1). We use the term domal salt to refer to the subsurface diapirs in the Zechstein basin. BAS, ZWD and WSN samples were collected and prepared for the KEM-17 study without gamma-irradiating them (Urai et al., 2019). The samples were selected in the TNO (Netherlands Organisation for Applied Scientific Research) central core storage facility in Zeist and core shed of the mining company (BAS). They were cut perpendicular to the bedding in a dry laboratory with a diamond saw cooled by a small amount of slightly undersaturated salt brine to reduce mechanical damage (Schléder and Urai, 2005).

To decorate crystal defect structures in NaCl, samples were irradiated in the research reactor FRM-II at the TU Munich in Garching with varying dose rates between 6 and 11 kGy h−1 to a total dose of 4 MGy at a constant temperature of 100 C (Urai et al., 1986b; Garcia Celma et al., 1988; Schléder and Urai, 2005, 2007). Thin sections of un-irradiated samples were dry-polished to a thickness of approximately 1 mm, and gamma-irradiated thin sections were dry-polished to a thickness of approximately 50 µm to reach optical transparency. To decorate grain boundaries and subgrain boundaries, the samples were chemically etched with slightly undersaturated brine, which was removed with a jet of n-hexane using the technique described in Spiers et al. (1986) and Urai et al. (1987). The thin sections were imaged in reflected and transmitted light using a Zeiss optical microscope (Axioscope) with the stitching panorama function of the ZEN imaging software. Reflected-light panoramas with 25× magnification were used for grain boundary digitization (Fig. 3) after making sure that the grain size distribution is sufficiently well captured.

Halite grain and subgrain boundaries were manually traced with a touch pen and tablet for statistical analysis. Grain and subgrain sizes for piezometry were analysed with Fiji (Schindelin et al., 2012) and calculated as equivalent circular diameter (Schléder and Urai, 2005; Lopez-Sanchez and Llana-Fúnez, 2015), and differential stresses σ were calculated as follows according to Carter et al. (1993) and Schléder and Urai (2005) with subgrain size =DSG:

Electron microscopy with energy-dispersive X-ray spectroscopy (SEM-EDS) mapping has been used to identify the chemical composition and distribution of second-phase impurities in halite. We used an Oxford Instruments X-Max 150 EDS system at the Institute of Structural Geology, Tectonics and Geomechanics (RWTH Aachen University) at 15 kV acceleration voltage. The samples were sputter-coated with approximately 7 nm of tungsten for conductivity.

Crystallographic orientations were measured by means of electron backscatter diffraction (EBSD) on one thin section using a Symmetry EBSD camera by Oxford Instruments attached to a GeminiSEM 300 SEM by Carl Zeiss Microscopy. The measurement was performed at an acceleration voltage of 20 kV and a probe current of approximately 10 nA. To avoid charging of the uncoated specimen, the analysis was performed in variable-pressure mode using nitrogen at a pressure of 30 Pa in the sample chamber. In a first measurement series, approximately half of the entire thin section was scanned by 30 small measurement areas with a step size of 20 µm, which afterwards were manually merged to a large dataset, finally covering an area of approximately 1.8 cm × 3 cm. The central Kristallbrocken region was subsequently measured again with higher spatial resolution by four measurements using a step size of 5 µm, which after merging covered an area of approximately 0.7 cm × 1.4 cm. All EBSD data were collected and indexed in “refined accuracy” mode with AZtec V 5.1 by Oxford Instruments. The final indexed datasets were further evaluated with OIM Analysis V 8.0 by Ametek–EDAX. Prior to the electron microscopical investigations, the samples were freshly polished and etched again using the method described above and coated with tungsten (in the case of the EDS measurements) or left uncoated (in the case of the EBSD measurements).

Qualitative X-ray diffraction analysis (XRD) to identify impurity content was performed on a Bruker D8 equipped with a graphite monochromator and a scintillation counter. Scans were measured with Cu–Kα radiation.

Figure 2(a) Reflected-light photograph of slabs (KS-07 from Barradeel) used for microstructural analysis showing layered white- and honey-coloured salt with variable grain sizes and milky Kristalllagen with internal lamination. (b) Slab specimen KS-15 of Barradeel used for microstructural analysis showing layered transparent and white salt with Kristalllagen up to 5 cm thickness with cracks. Small fold in 5 mm thick anhydrite band is bend following the displaced crystal layer. (c) Overview of thin gamma-decorated section of sample KS-07 in transmitted light (location indicated in a) showing Kristalllagen with abundant inclusions and decorated slip bands, dark-blue decorated bands parallel to crack, and layers of fine- and very fine-grained halite with white cores and blue rims. (d) overview of gamma-decorated thin section of sample KS-15 in transmitted light (location indicated in b) showing inclusion-poor crystal layers that are slightly displaced and layers of fine-grained halite, with some grains showing a white overgrowth. Anhydrite–halite band is slightly bent and contains abundant halite mineral inclusions of up to 2 mm. (e) Micrograph (location indicated in c) showing gamma-decorated slip bands and blue cellular structures at two different scales, abundant solid inclusions up to 200 µm, and fluid inclusions with gas bubbles. (f) Reflected-light image of Kristallbrocken and fine-grained halite with fluid inclusions and anhydrite at grain boundaries. (g) Transmitted-light micrograph showing white cores in halite with growth bands and slip bands. Finely dispersed anhydrite bands are next to very fine-grained halite. (h) Photomontage of transmitted-light image (with λ plate) and reflected-light image showing fluid inclusions at grain boundaries of fine-grained halite, polyhalite, anhydrite, and anhydrite band enriched with small anhydrite minerals.


Figure 3Overview of gamma-decorated samples from Barradeel, Pieterburen, Zuidwending, and Winschoten, together with maps of the interpreted microstructures. Reflected-light images that were used for interpretation of grain boundaries are presented in supplement 2 in Barabasch et al. (2022). Anhydrite, polyhalite, and gypsum layers were all mapped, but dispersed particles were mapped only partly as practicable.


3 Results

3.1 Relatively undeformed Barradeel samples

The Zechstein salt at Barradeel is sub-horizontally layered as seen in seismic data (Strozyk et al., 2014; Barabasch et al., 2019) and has a total thickness of 800 m with a 580 m thick Zechstein 2 (Z2) salt at a depth of up to 3 km based on well data (BAS-01; Kaart boringen | NLOG, 2022). Salt cores have mostly sub-horizontal layering, except for sample KS-15, where layering is vertical, indicating local folding, which is common in these settings (Fig. 2a and b).

The salt is layered with centimetre-scale milky to honey-coloured halite layers interbedded with millimetre-thick anhydrite bands (Fig. 2a). There are large variations in grain size of halite, with 1 to 5 cm thick single megacrystal layers (Kristalllagen) alternating with layers of fine-grained (about 3 mm) and very fine-grained (about 0.2 mm) halite layers (Fig. 2a and c).

Kristalllagen are clear or milky with internal laminae (Fig. 2a and b) due to variable content of fluid and solid inclusions of sulfate minerals (Fig. 2e) as previously described (Simon, 1972; Küster et al., 2011). The Kristalllagen have bedding perpendicular or inclined cleavage cracks without displacement of layering, the cracks presumably caused by drilling (Fig. 2a and c). However, in sample KS-15 the Kristalllage is displaced by a few millimetres along cracks next to small folds in anhydrite bands (Fig. 2b and d), which points to tectonic origin. Gamma decoration in Kristalllagen is mostly homogenous blue (Fig. 2c and d), except for a few dark laminae, which are interpreted as healed cracks in sample KS-07 that show the same orientation as an open crack in the same crystal. A Kristalllage in sample KS-07 has abundant slip bands oriented at 45 to the cleavage crack, consistent with slip on the [110] system (Fig. 2c). Kristalllagen are subgrain-free as seen in reflected light on etched surfaces (Fig. 2c). In sample KS-16 (Fig. 3) microstructures are very similar, with thicker anhydrite laminae. Cellular structures formed by dark-blue gamma decoration (Fig. 2e) locally coincide with barely visible subgrain boundaries on etched surfaces (supplement 1 in Barabasch et al., 2022) but do not correspond to the abundant inclusions. Some megacrystals have single 1 mm large subgrains with small angle misorientation of slip bands and hopper crystal cores indicated by arrays of cubic fluid inclusions with fluid and gas bubbles (supplement 1 in Barabasch et al., 2022).

The fine- and very fine-grained halite grains are equigranular, in layers of two distinct grain size classes of fine-grained (about 3 mm) and very fine-grained (about 0.2 mm) (Fig. 2c, g, h and supplement 1 in Barabasch et al., 2022). Grain boundaries are slightly curved, forming 120 angles at triple junctions (Fig. 2f and h). In reflected light, etched surfaces show that the only porosity in the salt is that of isolated fluid inclusions in grain boundaries (supplement 1 in Barabasch et al., 2022). Fine-grained halite grains have white cores and blue growth bands visible through gamma decoration (Fig. 2c and g). Most grains have bright crystal cores with a characteristic arrangement of cubic and chevron-shaped fluid inclusion trails and abundant impurities (Fig. 2g).

The layers with very fine grain sizes are rich in dispersed anhydrite and polyhalite (0.1 to 1 mm) located usually at grain boundaries (Fig. 2c, f, g, h).

Anhydrite layers of up to 5 mm thickness consist of micrometre- to millimetre-sized anhydrite and polyhalite grains and occasionally include 0.1 to 1 mm sized halite crystals (Fig. 2c, d, g, h). These layers are straight or locally folded and continuous.

3.2 Diapiric salt samples from Pieterburen, Winschoten, and Zuidwending

The studied Z2 rock salt comes from cores of three different salt diapirs (Winschoten, WSN; Zuidwending, ZWD; and Pieterburen, PBN) (Juez-Larré et al., 2019). Base salt reaches a depth of 3000 m, and the salt pierces the overburden up to depths of 100 m as in the case of Zuidwending diapir (Geluk et al., 2007). The bedding as seen in cores is mostly vertical and strongly deformed, with the older Z2 salt in the centre of the sampled diapirs, except for the Pieterburen structure, where the younger salt is in the centre (Geluk et al., 2007). The samples come from 260 to 1800 m depth.

The deformed diapiric salt consists of Kristallbrocken, surrounded by fine-grained and very fine-grained matrix halite and folded or boudinaged anhydrite bands (Fig. 3) together with dispersed anhydrite and polyhalite inclusions, also confirmed by XRD measurements (supplement 5 in Barabasch et al., 2022). These inclusions locally have a lighter-blue-coloured halite rim, but in the majority of cases the inclusions are not associated with colour changes in the surrounding halite crystal. The rock salt has a typical gneissic appearance consistent with earlier descriptions of Kristallbrockensalz (Küster et al., 2008, 2011, and references therein).

Figure 4(a) Transmitted-light salt core slab (Zuidwending) showing Kristallbrocken (KB), folded 3 mm thick anhydrite band in a fine-grained halite matrix, and unfolded but boudinaged 3 mm thick anhydrite band next to Kristallbrocken. (b) Gamma-irradiated transmitted-light scan (samples KS-13A and B) location indicated in (b) showing white cores and blue growth rims in fine-grained halite. (c) Salt core slab photograph of Pieterburen slab and photo-montaged transmitted-light image of thick section showing Kristallbrocken (KB) with clear halite on the grain edges of Kristallbrocken grains (modified after Sadler, 2012). (d) Transmitted-light scan of gamma-irradiated sample BS-A showing lighter-blue cores and dark-blue growth rims in a fine-grained halite matrix.

The Kristallbrocken megacrystals are between 5 mm and 3 cm thick (comparable to the thickness of Kristalllagen in Barradeel) with a strongly elongated shape up to decimetres long (possibly longer, but this could not be confirmed due to limitations of core dimensions; Fig. 4a–d). Kristallbrocken grains can be correlated to adjacent ones by thickness and inclusion content (Figs. 3, 4). In a number of cases, Kristallbrocken can be seen to contain a fracture along which the two megacrystals are displaced, and in other cases the two megacrystals are boudinaged with the boudin necks filled with inclusion-free overgrowth on the megacrystals or with fine-grained matrix halite in the boudin neck (Figs. 3, 4).

A good example is presented in Fig. 4d (BS-A), where two 1 cm thick Kristallbrocken parts can be correlated by thickness, enrichment of inclusions in the upper part, and adjacent anhydrite bands. The 2 cm boudin neck contains epitaxially grown (see EBSD below), inclusion-free halite from both crystals and fine-grained matrix halite (Fig. 4d). Both inclusion-free overgrowths also have a lighter-blue gamma irradiation colour. Two thin adjacent anhydrite layers are also discontinuous, in coherence with the boudinaged Kristallbrocken. A 3 cm thick Kristallbrocken below is also displaced and has a small inclusion-free overgrowth. Both inclusion-free parts of the crystal also have lighter-blue gamma irradiation colour.

Table 1Differential stresses from subgrain-size piezometry calculated with  (supplement 2 in Barabasch et al., 2022) after Carter et al. (1993) and Schléder and Urai (2005). KB: Kristallbrocken.

Download Print Version | Download XLSX

Figure 5(a) Transmitted-light micrograph of gamma-irradiated KS-08 (Winschoten) sample showing abundant dark inclusions in Kristallbrocken, slip bands, and white gamma-decorated subgrain boundaries. (b) Transmitted-light micrograph of gamma-irradiated KS-05 (Winschoten) sample showing abundant dark inclusions, slip bands, and white gamma-decorated subgrain boundaries in Kristallbrocken next to anhydrite band and inclusion-poor light halite in between. (c) Reflected-light image of sample KS-08 (Winschoten) Kristallbrocken with abundant subgrains and mineral inclusions used for piezometry. Exact location indicated in Fig. 3 (KS-08, s2). (d) Digitized subgrain boundaries are an example of the data used for piezometry. (e) Transmitted-light micrograph of gamma-irradiated sample KS-13B (Zuidwending) showing Kristallbrocken (KB) grain with gamma-decorated subgrains and abundant primary dark solid inclusions (white arrow). Clear halite grain shows array of reworked elongated clear anhydrite minerals parallel to grain boundary (white arrow) as well as abundant fluid inclusions at grain boundaries. (f) Transmitted-light micrograph of gamma-irradiated Zuidwending sample KS-02 showing subgrain structures, abundant dark inclusions, and clear elongated halite grains next to it. (g) Transmitted-light micrograph of gamma-decorated halite grains showing light cores and directional overgrowth of elongated fine-grained halite in sample BS-B (Pieterburen). (h) Transmitted-light micrograph of thick section showing abundant fluid inclusions at fine-grained halite grain boundary of Pieterburen sample (Sadler, 2012).

Gamma-decorated slip bands and subgrains are present in all Kristallbrocken, visible as white or dark-blue subgrain boundaries and in reflected light on etched surfaces (Fig. 5a, supplement 2 in Barabasch et al., 2022). Subgrain sizes are on average 100 to 200 µm (Table 1), and subgrains can be either rounded or elongated with preferred orientation of subgrain boundaries forming a subgrain boundary network. Figure 5d shows one of the interpreted subgrain networks, and Table 1 shows the subgrain-size statistics of all samples. Images of all other samples are included in supplement 2 in Barabasch et al. (2022), together with their interpretation and subgrain size distribution. Planar arrays of fluid inclusions are present in some Kristallbrocken (Fig. 5a). These are interpreted to be healed microcracks (supplement 1 in Barabasch et al., 2022). All Kristallbrocken contain dispersed inclusions of anhydrite and/or polyhalite and/or brine (Figs. 5a–f, 6), very similar to those observed and discussed by Küster et al. (2011).

The fine-grained matrix halite has a grain size of approximately 1 mm (Fig. 7). Halite grains are mostly elongated (KS-06, Kbx1, BS-B; Fig. 3). The fine-grained matrix halite grains commonly have white cores and blue rims (BS-A, KS-13A, KS-05, kbx2; Figs. 3, 4b, and d and 5g). An interesting observation is that in some of these grains (Fig. 5g), which are elongated, the white core has the blue overgrowth only in the direction of elongation (e.g. p. 36 in supplement 1 in Barabasch et al., 2022). Very fine-grained matrix halite grains are equiaxed with a similar white core and blue mantle to the fine-grained matrix halite and are associated with a high fraction of impurities, which have grain sizes of approximately 0.2 mm (KS-13A and B; Fig. 4b).

Figure 6(a) EDS map of KS-05 (Winschoten) showing anhydrite–polyhalite band, as well as Kristallbrocken with smaller anhydrite and polyhalite inclusions. (b) SEM micrograph of polyhalite and anhydrite inclusions in fine-grained halite. (c) SEM micrograph of Kristallbrocken broken along grain boundary showing fluid inclusion imprints and abundant magnesium and potassium sulfate inclusions.

Figure 7(a) Grain-size histograms of fine-grained matrix halite for each sample presented in Fig. 3 including statistical parameters. Grains smaller than 0.117 mm were not mapped due to image resolution. (b) Fractions of matrix halite, Kristallbrocken, and anhydrite/polyhalite and gypsum for each sample.


Figure 8(a) Inverse pole figure maps for halite grains with large Kristallbrocken grains labelled 1–5. (b) Pole figures of fine-grained matrix halite excluding five Kristallbrocken grains show no significant crystallographic preferred orientation (CPO). (c) Kernel average misorientation (KAM) map overlaid with white low- and high-angle grain boundaries (misorientation >5). KAM map was calculated over a distance of 40 µm (second neighbour) with a threshold of 3 in order to enhance the small-angle subgrain boundaries. KAM shows subgrain-free matrix halite with few exceptions in large matrix halite grains, Kristallbrocken 1 and 5 with subgrains, and subgrains in boudin necks of Kristallbrocken 3 and 4. (d) Cumulative reference orientation deviation map over the areas of Kristallbrocken 3 and 4, based on more highly resolved EBSD measurements. Reference points for each of the two grains are indicated. Panels (a) and (b) consist of 30 individual measurements, which due to image distortion under 70 tilt cannot be stitched perfectly. Therefore, in some cases an artificial separation of areas belonging to the same grain is visible.

While in the Barradeel samples the boundary between the Kristalllagen and the fine-grained halite is sharp, in the domal salt this boundary is less sharp, and the fine grains (without subgrains) are locally also present in the Kristallbrocken around their edges (KS-02, KS-11, KS-08; Fig. 3). These are interpreted to have formed by fluid-assisted grain boundary migration recrystallization, following the interpretation of Küster et al. (2008, 2011). The recrystallizing halite grains are free of fluid inclusions, with occasional second-phase inclusions at grain boundaries (BS-A, KS-13A, KS-13B, KS-05, and KS-02; Figs. 4b and d and 5b, e, f). Such recrystallized parts and the presumed primary fine-grained matrix halite cannot clearly be distinguished, and grain-size measurements regarding fine-grained matrix halite in Fig. 7 comprise both classes, referring to the halite-labelled grains in Fig. 3. The grain boundaries in halite in all samples are rich in fluid inclusions and occasional second-phase inclusions. In contrast to the observations of Küster et al. (2008), subgrains are infrequent in fine-grained salt in the samples studied, making up less than 1 % of the total matrix halite volume. When possible, subgrains were digitized (Fig. 3, supplement 1 in Barabasch et al., 2022), and piezometry results are presented in Table 1.

Thin anhydrite layers, when enclosed in fine-grained or very fine-grained matrix halite, are strongly deformed by folding and boudinage. The concentric folding in sample KS-13A (Fig. 4b) indicates that the viscosity of the anhydrite is much higher than that of the surrounding halite (Adamuszek et al., 2011, 2021). These structures indicate that the matrix halite was clearly more strongly deformed than the Kristallbrocken. Anhydrite layers directly adjacent to Kristallbrocken are much less deformed but are ruptured together with the Kristallbrocken to form coherent boudins.

3.3 EBSD results

EBSD data from sample BS-A are presented in Fig. 8. Figure 8a shows inverse pole figure maps for the points indexed as halite with respect to the horizontal direction and the direction perpendicular to the image plane, respectively. The large Kristallbrocken grains are labelled 1–5. The two parts of the halite boudin show different crystallographic orientation with a deviation of approximately 60 in the horizontal plane (supplement 3 in Barabasch et al., 2022). Further each Kristallbrocken (1–5) has different crystallographic orientation (Fig. 8a). Figure 8b shows pole figures of the matrix halite after excluding these five Kristallbrocken grains from the calculation. The small orientation distribution densities (<1.8 times random) show that the fine-grained matrix halite has no significant crystallographic preferred orientation (CPO), especially considering the fact that with approximately 1400 grains covered, the statistical base for calculation is comparatively small. The kernel average misorientation (KAM) map in Fig. 8c shows that the matrix halite is more or less subgrain-free, with very few exceptions in large matrix halite grains, whereas subgrains are present in Kristallbrocken 1 and 5 and in boudin necks of Kristallbrocken 3 and 4. An interesting aspect of the subgrain boundary network is their locally fibrous morphology in the overgrowths in boudin necks; these are interpreted as grown-in subgrains (Fig. 8) (Means and Ree, 1988). Figure 8d shows the cumulative reference orientation deviation over the areas of Kristallbrocken 3 and 4, based on the more highly resolved EBSD measurements of the central Kristallbrocken region. The corresponding reference points are marked with white crosses for each of the two grains. Colour jumps at the boundary of the individual measurement areas appear due to the movement of the specimen and subsequent recalibration between the individual measurements and should be neglected.

4 Discussion

4.1 Comparison to previous studies

Our observations of the Kristallbrockensalz samples correspond closely with those of earlier studies, the most extensive of which is the work of Küster (2011) on samples from the same stratigraphic unit, but several hundred kilometres further towards the basin interior. The key shared observations of both studies are internally deformed megacrystals with solid (polyhalite and anhydrite) and fluid inclusions, subgrains between 0.05 and 0.5 mm, and ruptured in extension to form boudins with dissolution–precipitation in the boudin necks. These megacrystals are surrounded by finer-grained halite, with a grain size around 1 mm and some grains containing subgrains and others being subgrain-free. Grain boundaries contain fluid inclusion arrays, and the matrix contains boudinaged and strongly folded thin anhydrite layers indicating high strain in the matrix. Following the interpretation of Löffler (1962) we interpret rupturing of boudins to be tectonic and in contrast to other hypotheses (Pape et al., 2002; Richter-Bernburg, 1953) that rupturing of halite crystal layers is not formed diagenetically from subaquatic gliding; although we cannot completely exclude fractures during diagenesis, the presence of abundant subgrains in Kristallbrocken requires differential stresses that cannot be generated close to the surface, and this deformation is consistent with formation of boudins.

The interpretation of these observations by Küster et al. (2008, 2011) was that Kristallbrockensalz resembles porphyroclasts in gneissic or mylonitic rocks with a strong rheology contrast between Kristallbrocken and matrix. Kristallbrocken formed from Kristalllagen (of diagenetic origin) by plastic deformation with active [110] slip accompanied by subgrain formation and by rupturing in extension, while the matrix salt was weaker and more strongly deformed than the Kristallbrocken and deformed by dislocation creep and dynamic recrystallization. We measured subgrain size in the studied micrographs and found comparable average subgrain sizes in Kristallbrocken (Fig. 2.10 in Küster, 2011) and in matrix halite (Fig. 7f in Küster et al., 2008), both from Teutschenthal, that indicate similar stress values of a few megapascals when compared to our samples (Table 1, Fig. 9). The same studies also suggest that matrix halite was secondary, replacing Kristallbrocken by recrystallization. The inferred strong rheology contrast was suggested to be related to “monocrystallinity” (large grain size) and the abundant inclusions in Kristallbrocken by restricting dislocation mobility (Küster et al., 2008).

Figure 9Differential stress vs. strain rate diagram plotting selected flow laws at 60 C. For dislocation creep BGRa KK, Kriechklassen 5, 4, and 3 calculated with Eq. (1),  s−1QDC=54 kJ mol−1, and n=5 (Liu et al., 2017) as well as Avery Island samples calculated with Eq. (1),  s−1QDC=68 kJ mol−1, and n=5.3 from Carter et al. (1993) are included. Dotted lines show combined pressure solution and dislocation creep flow laws for different halite grain sizes calculated with Eq. (3); previous values for dislocation creep (Eq. 1); and  s−1QPS=24.53 kJ mol−1, and m=3 (Eq. 2) (Spiers et al., 1990). Results from this study are plotted based on measured median fine-grained halite grain sizes (Fig. 7) and differential stresses from subgrain-size piezometry of Kristallbrocken for each sample (supplement 2 in Barabasch et al., 2022, Table 1).

Figure 10Comparison of differential stresses measured in Kristallbrocken and matrix halite for samples where both were available (), D= subgrain size; Carter et al., 1993; Schléder and Urai, 2005) with 95 % confidence intervals based on all measured subgrains per sample. Measurements indicate comparable differential stresses for both halite types, but slightly lower values for differential stress in Kristallbrocken of this study. Differential stresses from Teutschenthal were measured based on micrographs presented in Küster et al. (2008) for matrix halite and Küster (2011) for Kristallbrocken and show comparable, slightly lower differential stresses.

Our observations have added to the data of Küster et al. (2008, 2011) by providing high-resolution microstructures of gamma-irradiated samples and crystallographic orientation data by EBSD. In reflected-light micrographs we measured subgrain- and grain-size distributions (supplement 2 in Barabasch et al., 2022, Table 1, Fig. 7). The measured grain sizes were comparable, but slightly smaller than grain sizes measured on similar samples from the KEM-17 report (Urai et al., 2019), and stresses from piezometry between 0.8 and 1.5 MPa (Urai et al., 2019) also fit our observations. Our data suggest that subgrains in fine-grained halite are slightly smaller than subgrains in the associated Kristallbrocken (Fig. 9).

Based on results of both studies we also conclude that Kristallbrocken deformed by dislocation creep and that there is a large rheology contrast between Kristallbrocken and matrix, a widely occurring phenomenon in the Zechstein basin. However, our results of the Barradeel samples (which we interpret as an undeformed equivalent of the diapiric salt samples) clearly show that much of the fine- and very fine-grained halite was already present before the onset of salt tectonic deformation and not only formed by recrystallization of the Kristallbrocken. Because Küster et al. (2008, 2011) did not investigate undeformed equivalents, they did not have access to this information.

The original hypothesis by Küster et al. (2008, 2011) implies that the matrix salt is a dynamically recrystallized equivalent of the Kristallbrocken and deforms by equal contributions of dislocation creep and pressure solution creep (De Bresser et al., 1998, 2001; Ter Heege et al., 2005a, b). If this is the case, the predicted weakening would be relatively small (a factor of 2 in strain rate at the same differential stress) and could not explain the large differences in strain between Kristallbrocken and matrix. Additionally, if one would take the mean grain sizes from the matrix (Fig. 7), assume that these are recrystallized grain sizes, and then apply the recrystallized grain-size piezometer from Ter Heege et al. (2005a), one would get unrealistically high stresses between 5 and 10 MPa for the deformation. Large differences in strain could also be enhanced by the presence of solid inclusions in the Kristallbrocken as proposed by Küster et al. (2008, 2011); however this would lead to larger density of dislocation around the inclusions with corresponding darker colour after gamma irradiation (Garcia Celma et al., 1988), which was not observed. Finally, solid solution impurities could make the Kristallbrocken more competent (Heard and Ryerson, 1986). The most common solid solution impurity in halite is bromide, but the Kristallbrocken in a previous study have a comparable structure in samples where bromide concentration in Kristallbrocken is the same as in the matrix, and in samples where the bromide concentrations are different (Küster, 2011), again not supporting this hypothesis. However, we note that a more extensive analysis of solid solution impurity in these samples would help to further test this hypothesis. Another effect that may make the Kristallbrocken stronger is the presence of thin adjacent anhydrite layers which can form a stronger sandwich; however microstructures with or without sandwiched anhydrite layers are very similar, so this effect cannot be major either (samples BS-A and KS-13A in Fig. 3).

4.2 Dominant grain-size-dependent deformation mechanism and rheology contrast

Here we present an alternative hypothesis to explain the rheology contrast: we propose that the matrix deformed by dominant pressure solution creep, while the Kristallbrocken deformed by dislocation creep. Because pressure solution creep is strongly grain-size-dependent (see Eq. 2), this difference is caused by the large difference in grain size between Kristallbrocken and matrix, which we interpret to be sedimentary or early diagenetic, somewhat modified by dynamic recrystallization. Evidence for dominant pressure solution creep in the fine-grained and very fine-grained matrix is provided by the oriented overgrowth bands (Poirier, 1985) (Fig. 5g), the epitaxial overgrowths in boudin necks (Fig. 4d), and the absence of crystallographic preferred orientation as shown by EBSD (Fig. 8). Examples of rock salt microstructures indicating pressure solution creep were presented by Desbois et al. (2010) and Závada et al. (2012), and our microstructures are very similar to these. The presence of abundant fluid inclusions at grain boundaries and the evidence for fluid-assisted grain boundary migration recrystallization in both studies show that the required fluid is widely available in this stratigraphic unit. The studied diapiric samples have comparable mean grain sizes in the fine-grained salt around 1 or 2 mm (Fig. 7), while folding of anhydrite layers was best observed in Zuidwending samples, and overgrowth in boudin necks is especially pronounced in Pieterburen samples (Figs. 3, 4).

To quantify the stress–strain rate conditions corresponding to this hypothesis, in Fig. 10 we plot the data from this study in a differential stress vs. strain rate diagram for a reference temperature of 60 C. For dislocation creep, we used the BGR (German Federal Institute for Geosciences and Natural Resources)-recommended values (Liu et al., 2017; Eickemeier et al., 2021; Popp, 2022) for Kristallbrockensalz, together with a small selection of characteristic values for rock salts from the Kriechklassen (German Kriechklasse: creep classes) and plotted our samples for the values of differential stress obtained from subgrain-size piezometry, using the grain sizes of Kristallbrocken and matrix to estimate the expected strain rates. It can be seen that the corresponding strain rates differ by 4 orders of magnitude, in agreement with the inferred large differences in strain between Kristallbrocken and matrix (Fig. 10).

4.3 Dynamic recrystallization and grain boundary mobility

The dynamic recrystallization of Kristallbrocken is also recognized in grains directly next to Kristallbrocken (Figs. 4b and d and 5b, e, and f), with a very interesting change in the solid inclusions present in the Kristallbrocken; these are completely reworked by grain boundary migration (Fig. 5e) and not present in the same configuration in the new grains. In our interpretation the fraction of newly recrystallized grains in the matrix is relatively minor, and grain size is similar to the ones already present in the undeformed samples, providing the small grains for dominant pressure solution creep. The small grain size in the matrix is interpreted to be pre-deformational as in relatively undeformed Barradeel samples; hardly influenced by dynamic recrystallization as seen by only minor subgrains and a small grain size, which would indicate unrealistically high differential stresses when plotted in the Ter Heege et al. (2005a) recrystallized grain-size piezometer (Sect. 4.1). Recrystallization reduces the size of large single Kristallbrocken crystals and increases the relative content of finer-grained salt. Overgrowth of large Kristallbrocken crystals (increasing their size) has been observed in boudin necks. Hopper crystals indicating primary grains (Pape et al., 2002) were preserved in some BAS samples but not found in highly deformed domal salt. The relatively minor newly recrystallized grain fraction in the strongly deformed diapiric salt samples is interesting and suggests that grain boundary mobility in nature is lower than in the models of Peach et al. (2001) and Schléder and Urai (2005) and that recrystallization and grain growth in salt is more sluggish than previously thought.

4.4 Kristallbrockensalz bulk rheology

It is interesting to speculate on the deformation of initial layered salt with the very long Kristalllagen. During shearing, this rock must have had an extremely anisotropic rheology due to the weak fine-grained and very fine-grained salt depending on the orientation of the original layering with respect to the shortening direction. In coaxial deformation the Kristalllagen could have carried most of the stress in the salt rock and deformed accordingly by dislocation creep (Bons and Urai, 1994). Once the Kristalllagen were sufficiently fragmented, the rheology of the Kristallbrockensalz becomes a mixture rheology (Bons and Urai, 1994), such that when the Kristallbrocken can form a load-bearing framework, the rock salt deforms by power law creep dominated by dislocation creep in the Kristallbrocken, but when there is a sufficient matrix present, the rock salt deforms by Newtonian viscous creep dominated by pressure solution (Jessell et al., 2009). Based on subgrain-size analysis we inferred similar differential palaeostresses in Kristallbrocken and matrix halite (Fig. 9; even though we had a limited number of measurements in matrix halite), which suggest that in our samples the second model was dominant. The different crystallographic orientation of individual Kristallbrocken with respect to the bedding (Fig. 8a) has not been studied previously and might contribute to the understanding of Kristalllagen formation (Küster, 2011). The absence of a crystallographic preferred orientation of individual Kristallbrocken neither promotes nor prevents the activation of different slip systems and hence would not significantly influence the bulk rheology (Linckens et al., 2016). However, to substantiate this observation and hypothesis, more studies of Kristallbrocken crystallographic orientation with respect to bedding on different samples are required.

4.5 Implication of our results on engineering predictions

In current constitutive models used in salt engineering (Albrecht et al., 1993; Hunsche et al., 2003; Bräuer et al., 2011; Kukla et al., 2011; Liu et al., 2017; Popp, 2022) it was long recognized that the creep strain rate of rock salt (as measured at relatively high differential stress) can show differences of several orders of magnitude at the same differential stress and temperature. Based on an extensive dataset, such observations provided the basis for the definition of Kriechklassen, which were used to model the evolution of engineered structures (Liu et al., 2017) and are still used to describe a non-Newtonian rheology in salt tectonics numerical modelling (Granado et al., 2021). However, as is reviewed above, extrapolation of these data to low differential stress predicts creep rates that are orders of magnitude lower than measured in experiments (Brouard and Bérest, 1998; Bérest et al., 2019) and predicted by microphysical models of pressure solution creep (Spiers et al., 1986; Urai et al., 1986b; Spiers et al., 1990). Integration of these mechanisms into engineering constitutive equations is sometimes conducted (e.g. Zill et al., 2022; Buijze et al., 2022). However, because data on the grain size of the rock salt to be modelled are usually not available, these could still be improved considerably. Together with microstructural characterization, we recommend including grain sizes in engineering predictions based on grain-size-dependent pressure solution creep classes.

5 Conclusion

In this study, we present examples of flat-lying and diapiric Zechstein salt from the same formation (Z2), which was naturally deformed under low differential stresses between 1 and 2 MPa, providing a natural laboratory to study salt rheology under conditions which are difficult to study in the laboratory but relevant for predicting the evolution of engineered structures over long timescales. The gneissic Kristallbrockensalz deforms by dislocation creep and pressure solution processes depending on the grain size. The fine-grained matrix halite is weaker and deforms with a higher strain rate by pressure solution and dynamic recrystallization, while the Kristallbrocken mega grains tectonically boudinage and deform by dislocation creep. We infer that this large rheology contrast in halite deformation at low differential stresses is caused by grain-size-dependent dissolution–precipitation creep.

Further studies to better define solid solution impurity contents, the role of mineral impurities, and the contribution of dynamic recrystallization with grain boundary migration to solution precipitation processes will help to test the operating deformation mechanisms in more detail.

At present, grain-size measurements are not available for most of the Zechstein salts used in salt engineering, and we recommend the creation of a salt microstructure knowledge base which will help predict creep rates at low differential stress.

Data availability

The data set has been uploaded to the Zenodo repository and is available at (Barabasch et al., 2022). It is organized as follows:

  1. additional micrographs;
  2. high-resolution sample overview scans, piezometry images, and results;
  3. EBSD material;
  4. BIB-SEM material;
  5. XRD analysis results.

Author contributions

JB carried out the study and designed figures. JLU and JS designed the study. JK and JLU selected samples. AS did EBSD measurements and Fig. 8. JB and JLU prepared the manuscript with contributions from all co-authors, who were all involved in scientific discussions.

Competing interests

The contact author has declared that none of the authors has any competing interests.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


We would like to acknowledge the support of Nedmag Industries and RWTH Aachen University funding for completion of this study. Sample preparation by Werner Kraus was much appreciated. We further thank Marc Sadler for use of images in Figs. 4c and 5h from his MSc thesis (Sadler, 2012). We thank Prokop Závada and Hans de Bresser for reviews, which helped improve the manuscript.

Review statement

This paper was edited by Florian Fusseis and reviewed by Hans de Bresser and Prokop Závada.


Adamuszek, M., Schmid, D. W., and Dabrowski, M.: Fold geometry toolbox – Automated determination of fold shape, shortening, and material properties, J. Struct. Geol., 33, 1406–1416,, 2011. 

Adamuszek, M., Tămaş, D. M., Barabasch, J., and Urai, J. L.: Rheological stratification in impure rock salt during long-term creep: morphology, microstructure, and numerical models of multilayer folds in the Ocnele Mari salt mine, Romania, Solid Earth, 12, 2041–2065,, 2021. 

Albrecht, H., Hunsche, U., and Schulze, O.: Results from the application of the laboratory test program for mapping homogeneous parts in the Gorleben salt dome, 10th national rock mechanics symposium in Aachen, November 1992; Geotechnik-Sonderheft 1993, 155–158, Glückauf, Essen, 1993. 

Barabasch, J., Urai, J. L., Raith, A. F., and de Jager, J.: The early life of a salt giant: 3D seismic study on syntectonic Zechstein salt and stringer deposition on the Friesland Platform, Netherlands, Z. Dtsch. Geol. Ges., 170, 273–288,, 2019. 

Barabasch, J., Schmatz, J., Klaver, J., Schwedt, A., and Urai, J. L.: Supplementary material of Kristallbrockensalt microstructures of the Zechstein-2 rock salt from the Northern Netherlands, Zenodo [data set],, 2022. 

Baumann, T. S., Kaus, B. J. P., and Popov, A. A.: Deformation and stresses related to the Gorleben salt structure: Insights from 3D numerical models, in: Mechanical behavior of salt IX, edited by: Fahland, S., Hammer, J., Hansen, F., Heusermann, S., Lux, K.-H., and Minkley, W., BGR, Hannover, Germany, 597–609, 2018. 

Bérest, P., Béraud, J. F., Brouard, B., Blum, P. A., Charpentier, J. P., de Greef, V., Gharbi, H., and Valès, F.: Very slow creep tests on salt samples, EPJ Web of Conferences, 6, 22002,, 2010. 

Bérest, P., Gharbi, H., Brouard, B., Brückner, D., DeVries, K., Hévin, G., Hofer, G., Spiers, C., and Urai, J. L.: Very Slow Creep Tests on Salt Samples, Rock. Mech. Rock. Eng., 52, 2917–2934,, 2019. 

Bons, P. D. and Urai, J. L.: Experimental deformation of two-phase rock analogues, Mater. Sci. Eng, 175, 221–229,, 1994. 

Bräuer, V., Eickenmeier, R., Eisenburger, D., Grissemann, C., Hesser, J., Heusermann, S., Kaiser, D., Nipp, H.-K., Nowak, T., Plischke, I., Schnier, H., Schulze, O., Sönnke, J., and Weber, J. R.: Description of the Gorleben site part 4: Geotechnical exploration of the Gorleben salt dome, Bundesanstalt für Geowissenschaften und Rohstoffe (BGR), Hannover, Germany, ISBN 978-3-9814108-0-8, 2011. 

Brouard, B. and Bérest, P.: A tentative classification of salts according to their creep properties, SMRI Spring 1998 Meeting, 19–22 April 1998, New Orleans, Louisiana, USA, 18–38, 1998. 

Brouard, B., Bérest, P., Héas, J.-Y., Fourmaintraux, D., de Laguérie, P., and You, T.: An in situ creep test in advance of abandoning a salt cavern, SMRI Fall 2004 Technical Meeting, 3–5 October 2004, Berlin, Germany, 45–69, 2004. 

Buijze, L., Heege, J. T., and Wassing, B.: Finite Element modeling of natural sealing of wellbores in salt using advanced, laboratory-based salt creep laws, in: The Mechanical Behavior of Salt X, edited by: de Bresser, J. H. P., Drury, M. R., Fokker, P. A., Gazzani, M., Hangx, S. J. T., Niemeijer, A. R., and Spiers, C. J., CRC Press, ISBN 978-1-032-28220-6, 2022. 

Carter, N. L. and Hansen, F. D.: Creep of rocksalt, Tectonophysics, 92, 275–333,, 1983. 

Carter, N. L., Hansen, F. D., and Senseny, P. E.: Stress magnitudes in natural rock salt, J. Geophys. Res., 87, 9289,, 1982. 

Carter, N. L., Horseman, S. T., Russell, J. E., and Handin, J.: Rheology of rocksalt, J. Struct. Geol., 15, 1257–1271,, 1993. 

Chemia, Z., Schmeling, H., and Koyi, H.: The effect of the salt viscosity on future evolution of the Gorleben salt diapir, Germany, Tectonophysics, 473, 446–456,, 2009. 

Czapowski, G.: Sedimentary facies in the oldest rock salt (Na1) of the Leba elevation (Northern Poland), in: The Zechstein Facies in Europe: Lecture Notes in Earth Sciences, edited by: Peryt, T. M., Springer, Berlin, 10, 207–224, 1987. 

De Bresser, J. H. P., Peach, C. J., Reijs, J. P. J., and Spiers, C. J.: On dynamic recrystallization during solid state flow: Effects of stress and temperature, Geophys. Res. Lett., 25, 3457–3460,, 1998. 

De Bresser, J. H. P., Ter Heege, J. H., and Spiers, C. J.: Grain size reduction by dynamic recrystallization: can it result in major rheological weakening?, Int. J. Earth Sciences (Geol. Rundsch.), 90, 28–45,, 2001. 

Desbois, G., Závada, P., Schléder, Z., and Urai, J. L.: Deformation and recrystallization mechanisms in actively extruding salt fountain: Microstructural evidence for a switch in deformation mechanisms with increased availability of meteoric water and decreased grain size (Qum Kuh, central Iran), J. Struct. Geol., 32, 580–594,, 2010. 

Desbois, G., Urai, J. L., Schmatz, J., Závada, P., and de Bresser, J. H. P.: The distribution of fluids in natural rock salt to understand deformation mechanism, in: Mechanical behaviour of Salt VII, edited by: Bérest, P., Ghoreychi, M., Hadj-Hassen, F., and Tijani, M., Taylor & Francis Group, London, United Kingdom, 3–12, 2012. 

Dooley, T. P., Jackson, M. P. A., Jackson, C. A.-L., Hudec, M. R., and Rodriguez, C. R.: Enigmatic structures within salt walls of the Santos Basin – Part 2: Mechanical explanation from physical modelling, J. Struct. Geol., 75, 163–187,, 2015. 

Drury, M. R. and Urai, J. L.: Deformation-related recrystallization processes, Tectonophysics, 172, 235–253,, 1990. 

Eickemeier, R., Beese, S., Maniatis, G., and Fahland, S.: Sensitivitätsstudien zur gebirgsmechanischen Beurteilung der Integrität der Salzstockbarriere im Südfeld Bartensleben – Teil 5, Bundesanstalt für Geowissenschaften und Rohstoffe (BGR), Hannover, Germany, 2021. 

Fiduk, J. C. and Rowan, M. G.: Analysis of folding and deformation within layered evaporites in Blocks BM-S-8 & -9, Santos Basin, Brazil, in: Salt tectonics, sediments and prospectivity, edited by: Alsop, G. I., Archer, S. G., Hartley, A. J., Grant, N. T., and Hodgkinson, R., The Geological Society, London, United Kingdom, 471–487,, 2012. 

Garcia Celma, A., Urai, J. L., and Spiers, C. J.: A laboratory investigation into the interaction of recrystallization and radiation damage effects in polycrystalline salt rocks, Office for Official Publications of the European Communities, Luxembourg, 125 pp., 1988. 

Geluk, M. C.: Late Permian (Zechstein) carbonate-facies maps, the Netherlands, Geol. Mijnbouw, 79, 17–27,, 2000. 

Geluk, M. C., Paar, W. A., and Fokker, P. A.: Salt, in: Geology of the Netherlands, edited by: Wong, T. E., Batjes, D. A. J., and de Jager, J., Royal Netherlands Academy of Arts and Sciences, Amsterdam, Netherlands, 283–294, 2007. 

Granado, P., Ruh, J. B., Santolaria, P., Strauss, P., and Muñoz, J. A.: Stretching and Contraction of Extensional Basins With Pre-Rift Salt: A Numerical Modeling Approach, Front. Earth Sci., 9, 1–24,, 2021. 

Hampel, A.: Verbundprojekt: Vergleich aktueller Stoffgesetze und Vorgehensweisen anhand von Modellberechnungen zum thermo-mechanischen Verhalten und zur Verheilung von Steinsalz – Ergebnisbericht zum Teilvorhaben 1, Projektträger Karlsruhe, Wassertechnologie und Entsorgung (PTKA-WTE), Karlsruher Institut für Technologie (KIT), Mainz, Germany, 2016. 

Hampel, A., Schulze, O., Heemann, U., Zetsche, F., Günther, R.-M., Salzer, K., Minkley, W., Hou, M. Z., Wolters, R., Düsterloh, U., Zapf, D., Rokahr, R. B., and Pudewills, A.: BMBF-Verbundvorhaben: Die Modellierung des mechanischen Verhaltens von Steinsalz: Vergleich aktueller Stoffgesetze und Vorgehensweisen. Synthesebericht, Projektträger Forschungszentrum Karlsruhe (PTKA), Bereich Wassertechnologie und Entsorgung (WTE), Karlsruhe, Germany, 2007. 

Heard, H. C. and Ryerson, F. J.: Effect of cation impurities on steady-state flow of salt, in: Mineral and rock deformation: Laboratory studies, edited by: Hobbs, B. E. and Heard, H. C., American Geophysical Union, Washington, D.C., USA, 99–115,, 1986. 

Hudleston, P. J. and Treagus, S. H.: Information from folds: A review, J. Struct. Geol., 32, 2042–2071,, 2010. 

Hunsche, U. and Hampel, A.: Rock salt – the mechanical properties of the host rock material for a radioactive waste repository, Eng. Geol., 52, 271–291,, 1999. 

Hunsche, U., Schulze, O., Walter, F., and Plischke, I.: Thermomechanisches Verhalten von Salzgestein, report, Bundesanstalt für Geowissenschaften und Rohstoffe (BGR), Hannover, Germany, 2003. 

Jackson, M. P. A. and Hudec, M. R.: Salt Tectonics: Principles and Practice, Cambridge University Press,, 2017. 

Jackson, M. P. A., Vendeville, B. C., and Schultz-Ela, D. D.: Structural dynamics of salt systems, Annu. Rev. Earth Planet. Sci., 22, 93–117,, 1994. 

Jessell, M. W., Bons, P. D., Griera, A., Evans, L. A., and Wilson, C. J. L.: A tale of two viscosities, J. Struct. Geol., 31, 719–736,, 2009. 

Juez-Larré, J., van Gessel, S., Dalman, R., Remmelts, G., and Groenenberg, R.: Assessment of underground energy storage potential to support the energy transition in the Netherlands, First Break, 37, 57–66,, 2019. 

Kaart boringen | NLOG:, last access: 15 February 2022. 

Komoróczi, A., Abe, S., and Urai, J. L.: Meshless numerical modeling of brittle–viscous deformation: first results on boudinage and hydrofracturing using a coupling of discrete element method (DEM) and smoothed particle hydrodynamics (SPH), Comput. Geosci., 17, 373–390,, 2013. 

Kukla, P. A., Pechnig, R., and Urai, J. L.: Sichtung und Bewertung der Standortdaten Gorleben: Bericht zum Arbeitspaket 2: Vorläufige Sicherheitsanalyse für den Standort Gorleben, Gesellschaft für Anlagen- und Reaktorsicherheit (GRS) mbH, Cologne, Germany, ISBN 978-3-939355-52-6, 2011. 

Kukla, P. A., Urai, J. L., Raith, A., Li, S., Barabasch, J., and Strozyk, F.: The European Zechstein Salt Giant – Trusheim and Beyond, 37th Annual GCSSEPM Foundation Perkins-Rosen Research Conference, 3–6 December 2019, Huston, Texas, edited by: Fiduk, J. C. and Rosen, N. C., (last access: 7 February 2023), 2019. 

Küster, Y.: Bromide characteristics and deformation mechanisms of naturally deformed rock salt of the German Zechstein Basin, Dissertation, Georg-August-Universität, Göttingen, Germany, 221 pp.,, 2011. 

Küster, Y., Leiss, B., and Schramm, M.: Structural characteristics of the halite fabric type `Kristallbrocken' from the Zechstein Basin with regard to its development, Int. J. Earth Sci., 99, 505–526,, 2008. 

Küster, Y., Schramm, M., Bornemann, O., and Leiss, B.: Bromide distribution characteristics of different Zechstein 2 rock salt sequences of the Southern Permian Basin: a comparison between bedded and domal salts, Sedimentology, 56, 1368–1391,, 2009. 

Küster, Y., Schramm, M., and Leiss, B.: Compositional and microstructural characterisation of solid inclusions in the laminated halite type “Kristallbrocken” with regard to its formation in the Central European Zechstein Basin, Z. Dt. Ges. Geowiss., 162, 277–294,, 2011. 

Laier, T., Kockel, F., Geluk, M. C., Pokorsky, J., and Lott, G.K.: Section A (Geology), in: NW European Gas Atlas, edited by: Lokhorst, A., NITG-TNO, Haarlem, ISBN 90-72869-60-5, 1998. 

Leitner, C., Neubauer, F., Urai, J. L., and Schoenherr, J.: Structure and evolution of a rocksalt-mudrock-tectonite: The haselgebirge in the Northern Calcareous Alps, J. Struct. Geol., 33, 970–984,, 2011. 

Linckens, J., Zulauf, G., and Hammer, J.: Experimental deformation of coarse-grained rock salt to high strain, J. Geophys. Res.-Sol. Ea., 121, 6150–6171,, 2016. 

Liu, W., Völkner, E., Minkley, W., and Popp, T.: Zusammenstellung der Materialparameter für THM-Modellberechnungen: Ergebnisse aus dem Vorhaben KOSINA, report, Bundesanstalt für Geowissenschaften und Rohstoffe (BGR), Hannover, Germany, 2017. 

Löffler, J.: Zur Genese der Augensalze im Zechstein der Deutschen Demokratischen Republik, Z. Angew. Geol., 11, 583–589,, 1962. 

Lopez-Sanchez, M. A. and Llana-Fúnez, S.: An evaluation of different measures of dynamically recrystallized grain size for paleopiezometry or paleowattometry studies, Solid Earth, 6, 475–495,, 2015. 

Macente, A., Fusseis, F., Butler, I. B., Tudisco, E., Hall, S. A., and Andò, E.: 4D porosity evolution during pressure-solution of NaCl in the presence of phyllosilicates, Earth Planet. Sci. Lett., 502, 115–125,, 2018. 

Means, W. D. and Ree, J. H.: Seven types of subgrain boundaries in octachloropropane, J. Struct. Geol., 10, 765–770,, 1988. 

Pape, T., Michalzik, D., and Bornemann, O.: Chevronkristalle im Kristallbrockensalz (Zechstein 2) des Salzstocks Gorleben – Primärgefüge salinarer Flachwassersedimentation im Zechsteinbecken, Z. Dt. Ges. Geowiss., 115–129,, 2002. 

Peach, C. J., Spiers, C. J., and Trimby, P. W.: Effect of confining pressure on dilatation, recrystallization, and flow of rock salt at 150 C, J. Geophys. Res., 106, 13315–13328,, 2001. 

Peel, F. J.: The engines of gravity-driven movement on passive margins: Quantifying the relative contribution of spreading vs. gravity sliding mechanisms, Tectonophysics, 633, 126–142,, 2014. 

Poirier, J.-P.: Creep of Crystals: High-Temperature Deformation Processes in Metals, Ceramics and Minerals, 1st edn., Cambridge University Press,, 1985. 

Popp, T.: Eigenschaften und Potential stratiformer Salz-Formationen für die Endlagerung hochradioaktiver Abfälle, report, Institut für Gebirgsmechanik IfG, Leipzig, Germany, 2022. 

Popp, T. and Hansen, F. D.: Creep at low deviatoric stress, in: Proceedings of the 8th US/German workshop on salt repository research, design, and operation, edited by: Hansen, F. D., Steininger, W., Bollingerfehr, W., Kuhlman, K., and Dunagan, S., Sandia National Laboratories, Albuquerque, New Mexico, USA, 21–24, 2018. 

Popp, T., Kern, H., and Schulze, O.: Lithologische Variabilität der petrophysikalischen und mineralogisch-gefügekundlichen Eigenschaften des älteren Steinsalzes (z2) aus dem Salzstock Gorleben, Meyniana, 51, 55–75, 1999. 

Richter-Bernburg, G.: Über salinare Sedimentation, Z. Dt. Ges. Geowiss., 105, 593–645, 1953. 

Rowan, M. and Krzywiec, P.: The Szamotuły salt diapir and Mid-Polish Trough: Decoupling during both Triassic-Jurassic rifting and Alpine inversion, Interpretation, 2, SM1–SM18,, 2014. 

Rowan, M. G., Urai, J. L., Fiduk, J. C., and Kukla, P. A.: Deformation of intrasalt competent layers in different modes of salt tectonics, Solid Earth, 10, 987–1013,, 2019. 

Sadler, M.: The distribution and morphology of grain boundary fluids in natural rock salt, MSc, RWTH Aachen University, 2012. 

Schenk, O. and Urai, J. L.: Microstructural evolution and grain boundary structure during static recrystallization in synthetic polycrystals of sodium chloride containing saturated brine, Contrib. Mineral. Petrol., 146, 671–682,, 2004. 

Schenk, O., Urai, J. L., and Piazolo, S.: Structure of grain boundaries in wet, synthetic polycrystalline, statically recrystallizing halite – evidence from cryo-SEM observations, Geofluids, 6, 93–104,, 2006. 

Schindelin, J., Arganda-Carreras, I., Frise, E., Kaynig, V., Longair, M., Pietzsch, T., Preibisch, S., Rueden, C., Saalfeld, S., Schmid, B., Tinevez, J.-Y., White, D. J., Hartenstein, V., Eliceiri, K., Tomancak, P., and Cardona, A.: Fiji: an open-source platform for biological-image analysis, Nat. Methods, 9, 676–682,, 2012. 

Schléder, Z. and Urai, J. L.: Microstructural evolution of deformation-modified primary halite from the Middle Triassic Röt Formation at Hengelo, The Netherlands, Int. J. Earth Sci. (Geol. Rundsch.), 94, 941–955,, 2005. 

Schléder, Z. and Urai, J. L.: Deformation and recrystallization mechanisms in mylonitic shear zones in naturally deformed extrusive Eocene–Oligocene rocksalt from Eyvanekey plateau and Garmsar hills (central Iran), J. Struct. Geol., 29, 241–255,, 2007. 

Schmalholz, S. M. and Mancktelow, N. S.: Folding and necking across the scales: a review of theoretical and experimental results and their applications, Solid Earth, 7, 1417–1465,, 2016. 

Schmalholz, S. M. and Podladchikov, Y. Yu.: Strain and competence contrast estimation from fold shape, Tectonophysics, 340, 195–213,, 2001. 

Schmatz, J., Schenk, O., and Urai, J. L.: The interaction of migrating grain boundaries with fluid inclusions in rock analogues: the effect of wetting angle and fluid inclusion velocity, Contrib. Mineral. Petrol., 162, 193–208,, 2011. 

Schultz-Ela, D. D., Jackson, M. P. A., and Vendeville, B. C.: Mechanics of active salt diapirism, Tectonophysics, 228, 275–312,, 1993. 

Schwichtenberg, B., Fusseis, F., Butler, I. B., and Andò, E.: Biotite supports long-range diffusive transport in dissolution–precipitation creep in halite through small porosity fluctuations, Solid Earth, 13, 41–64,, 2022. 

Seidl, E.: Die permische salzlagerstätte im Graf Moltke schalcht und in der umgebung von Schönebeck a. d. Elbe: Beziehung zwischen mechanismus der gebirgsbildung und innerer umformung der salzlagerstätte, Königlich Preußische Geologische Landesanstalt, Berlin, Germany, 104 pp., 1914. 

Simon, P.: Stratigraphie und Bromgehalt des Staßfurt-Steinsalzes (Zechstein 2) im hannoverschen Kalisalzbergbaugebiet, Geol. Jhrbch, 90, 67–126, 1972. 

Spiers, C. J. and Carter, N. L.: Microphysics of rocksalt flow in nature, in: 4th Conference on the Mechanical Behavior of Salt, edited by: Aubertin, M. and Hardy, Jr., H. R., Trans Tech Publications, Clausthal-Zellerfeld, Germany, 115–128, 1998. 

Spiers, C. J., Urai, J. L., Lister, G. S., Boland, J. N., and Zwart, H. J.: The influence of fluid-rock interaction on the rheology of salt rock, Commission of the European Communities, Luxembourg, EUR 10399 EN, ISBN 92-825-6214-X, 1986. 

Spiers, C. J., Schutjens, P. M. T. M., Brzesowsky, R. H., Peach, C. J., Liezenberg, J. L., and Zwart, H. J.: Experimental determination of constitutive parameters governing creep of rocksalt by pressure solution, in: Deformation mechanisms, rheology and tectonics, edited by: Knipe, R. J. and Rutter, E. H., The Geological Society, London, United Kingdom, 215–227,, 1990. 

Strozyk, F., Urai, J. L., van Gent, H., de Keijzer, M., and Kukla, P. A.: Regional variations in the structure of the Permian Zechstein 3 intrasalt stringer in the northern Netherlands: 3D seismic interpretation and implications for salt tectonic evolution, Interpretation, 2, SM101–SM117,, 2014. 

Talbot, C. J.: Fold trains in a glacier of salt in southern Iran, J. Struct. Geol., 1, 5–18,, 1979. 

Talbot, C. J. and Aftabi, P.: Geology and models of salt extrusion at Qum Kuh, central Iran, J. Geol. Soc., London, 161, 321–334,, 2004. 

Tămaṡ, D. M., Tăma?, A., Barabasch, J., Rowan, M. G., Schléder, Z., Krézsek, C., and Urai, J. L.: Low-Angle Shear Within the Exposed Mânzăleṡti Diapir, Romania: Salt Decapitation in the Eastern Carpathians Fold-and-Thrust Belt, Tectonics, 40, e2021TC006850,, 2021. 

Ter Heege, J. H., De Bresser, J. H. P., and Spiers, C. J.: Dynamic recrystallization of wet synthetic polycrystalline halite: dependence of grain size distribution on flow stress, temperature and strain, Tectonophysics, 396, 35–57,, 2005a. 

Ter Heege, J. H., de Bresser, J. H. P., and Spiers, C. J.: Rheological behaviour of synthetic rocksalt: the interplay between water, dynamic recrystallization and deformation mechanisms, J. Struct. Geol., 27, 948–963,, 2005b. 

Trusheim, F.: Über Halokinese und ihre Bedeutung für die strukturelle Entwicklung Norddeutschlands, Z. Dt. Ges. Geowiss., 109, 111–158, 1957. 

Urai, J. L. and Spiers, C. J.: The effect of grain boundary water on deformation mechanisms and rheology of rocksalt during long-term deformation, in: The mechanical behavior of salt – understanding of THMC processes in salt, edited by: Wallner, M., Lux, K.-H., Minkley, W., and Hardy Jr., H. R., Taylor & Francis, Hannover, Germany, 149–158,, 2007. 

Urai, J. L., Means, W. D., and Lister, G. S.: Dynamic Recrystallization of Minerals, in: Mineral and Rock Deformation, American Geophysical Union (AGU), 161–199,, 1986a. 

Urai, J. L., Spiers, C. J., Zwart, H. J., and Lister, G. S.: Weakening of rock salt by water during long-term creep, Nature, 324, 554–557,, 1986b. 

Urai, J. L., Spiers, C. J., Peach, C. J., Franssen, R. C. M. W., and Liezenberg, J. L.: Deformation mechanisms operating in naturally deformed halite rocks as deduced from microstructural investigations, Geol. Mijnbouw, 66, 165–176, 1987. 

Urai, J. L., Schléder, Z., Spiers, C. J., and Kukla, P. A.: Flow and transport properties of salt rocks, in: Dynamics of complex intracontinental basins: The central European basin system, edited by: Littke, R., Bayer, U., Gajewski, D., and Nelskamp, S., Springer, Berlin, Germany, 277–290, 2008. 

Urai, J. L., Schmatz, J., and Klaver, J.: Report, Project KEM-17 Over-pressured salt solution mining caverns and leakage mechanisms, report, Ministry of Economic Affairs and Climate, the Netherlands, 2019. 

van Gent, H. W., Urai, J. L., and de Keijzer, M.: The internal geometry of salt structures – A first look using 3D seismic data from the Zechstein of the Netherlands, J. Struct. Geol., 33, 292–311,, 2011. 

van Keken, P. E., Spiers, C. J., van den Berg, A. P., and Muyzert, E. J.: The effective viscosity of rocksalt: Implementation of steady-state creep laws in numerical models of salt diapirism, Tectonophysics, 225, 457–476,, 1993. 

van Noort, R., Visser, H. J. M., and Spiers, C. J.: Influence of grain boundary structure on dissolution controlled pressure solution and retarding effects of grain boundary healing, J. Geophys. Res., 113, B03201,, 2008. 

Wawersik, W. R. and Zeuch, D. H.: Modeling and mechanistic interpretation of creep of rock salt below 200 C, Tectonophysica, 121, 125–152,, 1986.  

Wenkert, D.: The flow of salt glaciers, Geophys. Res. Lett., 6, 523–526, 1979. 

Závada, P., Desbois, G., Schwedt, A., Lexa, O., and Urai, J. L.: Extreme ductile deformation of fine-grained salt by coupled solution-precipitation creep and microcracking: Microstructural evidence from perennial Zechstein sequence (Neuhof salt mine, Germany), J. Struct. Geol., 37, 89–104,, 2012. 

Závada, P., Desbois, G., Urai, J. L., Schulmann, K., Rahmati, M., and Lexa, O.: Impact of solid second phases on deformation mechanisms of naturally deformed salt rocks (Kuh-e-Namak, Dashti, Iran) and rheological stratification of the Hormuz Salt Formation, J. Struct. Geol., 74, 117–144, 2015. 

Zill, F., Wang, W., and Nagel, T.: Influence of THM process coupling and constitutive models on the simulated evolution of deep salt formations during glaciation, in: The Mechanical Behavior of Salt X, edited by: de Bresser, J. H. P., Drury, M. R., Fokker, P. A., Gazzani, M., Hangx, S. J. T., Niemeijer, A. R., and Spiers, C. J., CRC Press, London, 353–362,, 2022. 

Zulauf, J., Zulauf, G., Göttlich, J., and Peinl, M.: Formation of chocolate-tablet boudins: Results from scaled analogue models, J. Struct. Geol., 68, 97–111,, 2014.