Dissolution-induced Subsidence Research Articles

Closed depressions in Kotido crater, Arabia Terra, Mars. Possible evidence of evaporite dissolution-induced subsidence

The identification of karst sinkholes in Mars may provide evidence of dissolution processes caused by liquid water and information on paleoclimatic and paleohydrological conditions. This work presents a comprehensive cartographic inventory of 513 closed depressions developed on evaporite-bearing Equatorial Layered Deposits (ELDs) within Kotido crater, Arabia Terra. Detailed mapping, morphometric analyses and spatial distribution relationships reveal a number of features supporting that the depressions correspond to collapse sinkholes related to evaporite dissolution: (1) suitable topographic and litho-structural conditions for the development of a fracture-controlled epigene evaporite karst; (2) presence of open fissures at the foot of the scarped margins; (3) dimensions and frequency-size distributions comparable with those reported on Earth; (4) spatial association with high-permeability zones (i.e., fractures). Some characteristics of the depressions indicate that they have been re-shaped and enlarged by wind erosion: (1) dominant orientation consistent with the prevalent one-directional winds; (2) differing morphological characteristics on the downwind- and upwind-sides; and (3) nested depressions associated with the upwind sector. The relatively fresh appearance of the depressions and the lack of impact craters suggest a poorly constrained Amazonian karstification phase in the region.

Gypsum scarps and asymmetric fluvial valleys in evaporitic terrains. The role of river migration, landslides, karstification and lithology (Ebro River, NE Spain)

Most of the Spanish fluvial systems excavated in Tertiary evaporitic gypsum formations show asymmetric valleys characterized by a stepped sequence of fluvial terraces on one valley flank and kilometric-long and >100-m high prominent river scarp on the opposite side of the valley. Scarp undermining by the continuous preferential lateral migration of the river channel toward the valley margin leads to vertical to overhanging unstable slopes affected by a large number of slope failures that become the main geological hazard for villages located at the toe of the scarps. Detailed mapping of the gypsum scarps along the Ebro and Huerva Rivers gypsum scarps demonstrates that landslides and lateral spreading processes are predominant when claystones crop out at the base of the scarp, while rockfalls and topples become the dominant movement in those reaches where the rock mass is mainly constituted by evaporites. The dissolution of gypsum nodules, seasonal swelling and shrinking, and dispersion processes contribute to a decrease in the mechanical strength of claystones. The existence of dissolution-enlarged joints, sinkholes, and severely damaged buildings at the toe of the scarp from karstic subsidence demonstrates that the interstratal karstification of evaporites becomes a triggering factor in the instability of the rock mass. The genesis of asymmetric valleys and river gypsum scarps in the study area seem to be caused by the random migration of the river channel in the absence of lateral tilting related to tectonics or dissolution-induced subsidence. Once the scarp is developed, its preservation depends on the physicochemical properties of the substratum, the ratio between bedrock erosion and river incision rates, and climatic conditions that favour runoff erosion versus dissolution.

Railway deformation detected by DInSAR over active sinkholes in the Ebro Valley evaporite karst, Spain

Abstract. Subsidence was measured for the first time on railway tracks in the central sector of Ebro Valley (NE Spain) using Differential Synthetic Aperture Radar Interferometry (DInSAR) techniques. This area is affected by evaporite karst and the analysed railway corridors traverse active sinkholes that produce deformations in these infrastructures. One of the railway tracks affected by slight settlements is the Madrid–Barcelona high-speed line, a form of transport infrastructure highly vulnerable to ground deformation processes. Our analysis based on DInSAR measurements and geomorphological surveys indicates that this line shows dissolution-induced subsidence and compaction of anthropogenic deposits (infills and embankments). Significant sinkhole-related subsidence was also measured by DInSAR techniques on the Castejón–Zaragoza conventional railway line. This study demonstrates that DInSAR velocity maps, coupled with detailed geomorphological surveys, may help in the identification of the railway track sections that are affected by active subsidence.

Assessing sinkhole activity in the Ebro Valley mantled evaporite karst using advanced DInSAR

Sinkholes in karst areas may cause subsidence damage in transportation infrastructures, demolition of buildings and even the loss of human lives when they occur in a catastrophic way. Differential Interferometry (DInSAR) is a promising technology for detecting and characterizing sinkholes, as well as for reducing the associated risk when combined with other sources of data such as a sinkhole inventory. In this work, the usefulness of InSAR techniques and data for sinkhole risk management has been analyzed through the comparison of three DInSAR-derived velocity maps with a comprehensive sinkhole inventory in the Ebro Valley, NE Spain. The DInSAR maps have contributed to improve the sinkhole inventory in different ways: (1) detection of non-inventoried sinkholes; (2) revision of sinkhole areas previously classified as inactive as active; and (3) refinement of underestimated sinkhole boundaries. The obtained results suggest that DInSAR products are suitable for analyzing active dissolution-induced subsidence. The application of these techniques may help in recognizing and better characterizing previously unknown karst subsidence problems and in preventing personal and property losses. However, the analysis reveals that the available DInSAR maps combined overlook about 70% of the previously mapped active sinkholes mainly due to decorrelation.

Stratigraphy of the Arriaga Palaeolithic sites. Implications for the geomorphological evolution recorded by thickened fluvial sequences within the Manzanares River valley (Madrid Neogene Basin, Central Spain)

The Arriaga Palaeolithic sites, located within the Middle–Late Pleistocene thickened terrace (TCMZ: +18–22m) of the Manzanares River valley (Madrid, Central Spain), were subject to intensive archaeological and palaeontological prospecting during the 1980s. Compilation of documents from these old excavations, together with new geoarchaeological, sedimentological, pedological and geophysical data, allow us to locate the morpho-stratigraphic position of the analysed sites within the overall stratigraphy of the TCMZ. This thickened terrace comprises two main fluvial sequences (Lower and Upper) topped by a thick (2.5–5m) alluvial–colluvial formation. The fluvial sequences are stacked in the study site located in the lowermost reach of the valley, but display complex inset relationships upstream, where they are individualized in two different terrace levels at +18–22 and +12–15m. Terrace thickening was primarily controlled by synsedimentary subsidence caused by dissolution of the evaporitic substratum and locally influenced and backfed by tectonic activity. The regional analysis of the dated (TL and OSL) fluvial sequences containing Palaeolithic sites within the TCMZ, together with new TL dates provided in this study, indicate that the three sedimentary sequences in the TCMZ are time-transgressive valley-fill bodies. Terrace thickening started before the Last Interglacial Period (MIS 6 or older) and continued during whole MIS 5 (lower fluvial sequence) and MIS 4 (upper fluvial sequence) reaching the MIS 3 (top alluvial formation), the latter characterized by the accumulation of alluvial–colluvial sequences derived from the main tributaries and valley slopes. The TCMZ records the Middle–Late Pleistocene boundary, but also the transition between the Lower and Middle Palaeolithic periods during the Late MIS 5 (ca. 96 to 74ka). The studied Arriaga sites contain evolved Lower Palaeolithic industry (evolved Acheulean techno-complexes) and warm faunal assemblages located within the Lower fluvial sequence, but apparently well constrained Middle Palaeolithic sites are placed within the Upper fluvial sequence at other upstream locations. Deposition of the thickened alluvium was mainly controlled by the upstream advance of dissolution-induced subsidence phenomena, blurring the impact of Late Pleistocene climatic cycles and producing time-transgressive longitudinal valley-fill bodies (i.e. sedimentary sequences). Late Quaternary climatic changes only seem to control the incision/aggradation cycles after the termination of the TCMZ from the Late MIS 3. Dates related to the development of younger inset terraces indicate that they are apparently linked with cold Heinrich events H4 to H1. These younger inset terraces yield cold faunal assemblages and abundant Middle Palaeolithic “Mousterian” assemblages.

Geomorphic and stratigraphic evidence of incision-induced halokinetic uplift and dissolution subsidence in transverse drainages crossing the evaporite-cored Barbastro–Balaguer Anticline (Ebro Basin, NE Spain)

The evaporite-cored Barbastro–Balaguer Anticline, located in the northeastern sector of the Ebro Tertiary Basin, NE Spain, is traversed perpendicularly by several drainages coming from the Pyrenees, including the Cinca, Noguera-Ribagorzana and Segre rivers. The terraces associated with these discordant fluvial systems, deposited upon evaporitic and detrital formations across the anticline, allow the analysis of geomorphic and stratigraphic anomalies related to the soluble and ductile nature of the halite-bearing rocks in the core of the anticline. Geomorphological mapping reveals that the gravels discontinuously capping the evaporitic core of the anticline, together with some alluvial mantles disconnected from the present-day river valleys, constitute the oldest sediments deposited under exorheic conditions in this sector of the Ebro Tertiary Basin. These alluvial deposits inset into the basin fill record an early post-capture palaeogeographic stage during which the drainage was dominated by unconfined distributary channels developed on extensive alluvial fans fed by perennial rivers coming from the Pyrenees. The terraces associated with the present-day fluvial valleys correspond to a later palaeogeographic stage characterized by a well-integrated drainage confined to fluvial valleys entrenched in the basin fill. These terraces show evidence of both dissolution-induced subsidence and halokinesis restricted to areas underlain by the evaporites of the anticline core. The deposits of the highest terrace levels of the Noguera-Ribagorzana and Segre rivers and its tributaries, Lo Reguer Creek and Farfanya River, are locally thickened filling basins generated by dissolution-induced synsedimentary subsidence up to several kilometers long and more than 100m deep. Subsidence caused by the karstification of the evaporites reached higher magnitude and extent during the early stages of development of the new external drainage network, when the more soluble halite-bearing units were closer to the new incisional base level. Most likely, the development of a thick caprock, rich in less soluble sediments, has hampered dissolution and subsidence processes in subsequent stages. Incision-induced halokinetic deformation is locally recorded by differentially uplifted terraces dipping away from the valley. The increasing differential loading caused by the entrenchment of the fluvial systems has driven the lateral and upward flow of salt-rich evaporites from the valley margins towards the valley bottom. A minimum long-term uplift rate of 0.3mm/yr has been calculated for a titled OSL-dated Upper Pleistocene terrace of the Cinca River valley. The capture of the Ebro Basin and its change to exorheic conditions have played an instrumental role in the development of the Quaternary gravitational deformations associated with the evaporitic bedrock. The new incisional drainage network propitiated the karstification and evacuation of large amounts of evaporites, and its entrenchment has created adequate differential loading conditions for the salt to flow towards the valley floors.

Environmental problems and geological implications derived from evaporite dissolution in the Barbastro salt anticline (NE Spain)

The halite-bearing Barbastro Formation crops out in the core of the Barbastro Anticline (Ebro Tertiary Basin). This anticline is traversed perpendicularly by some of the most important Pyrenean drainages such as the Cinca and Noguera-Ribagorzana Rivers. The terrace sequences of these fluvial systems have been used as markers to identify and assess dissolution-induced subsidence and salt tectonics. In the limbs of the anticline, terrace deposits underlain by detrital bedrock do not show any evidence of deformation and have a consistent thickness of less than 10 m. The deposits of certain terrace levels of the Noguera-Ribagorzana River and its tributary, the Lo Reguer Creek, are locally thickened filling basins generated by dissolution-induced synsedimentary subsidence up to several kilometers long and more than 100 m deep. Conversely, terraces of the Cinca River do not show anomalously high thicknesses, but local uplifts related to differential upward flow of the halite-bearing bedrock. Locally, a minimum uplift rate of 0.3 mm/year has been estimated from a 64-ka terrace tilted away from the valley. The subsidence hazards occur chiefly in areas where the ground receives artificial water recharge. Serviceability of some canals has been notoriously affected by evaporite karstification. The problem has been mitigated to acceptable levels by grouting. Numerous buildings of Ivars de Noguera are severely damaged by dissolution subsidence, and possibly, by hydrocompaction of gypsiferous silts. The pipe network has been replaced to ameliorate the subsidence risk. In the Cinca River valley, cavities with a total volume of about 180,500 m3 have been created by solution mining at depths greater than 500 m. No investigation methods are applied in the brine field to monitor the distribution and evolution of artificial voids. Substantial increase in salinity of the Cinca River is another evidence of subjacent evaporite dissolution.

Origin of the salt valleys in the Canyonlands section of the Colorado Plateau: Evaporite-dissolution collapse versus tectonic subsidence

The salt valleys over the axis of the salt-cored anticlines in the Paradox fold and fault belt (Canyonlands, Utah and Colorado) are created by subsidence of the anticline crests. Traditionally, the collapse of the anticlinal crests was attributed to dissolution of the salt walls (diapirs) forming the anticline cores. Recent studies based on scaled physical models and field observations propose that the salt valleys are a result of regional extension and that salt dissolution had only a minor influence in the development of the axial depressions. This paper presents several arguments and lines of evidence that refute the tectonic model and support the salt dissolution subsidence interpretation. The development of contractional structures in salt dissolution experiments led the advocates of the tectonic interpretation to reject the dissolution-induced subsidence explanation. However, these salt dissolution models do not reproduce the karstification of salt walls in a realistic way, since their analog involves removal of salt from the base of the diapirs during the experiments. Additionally, numerous field examples and laboratory models conducted by other authors indicate that brittle subsidence in karst settings is commonly controlled by subvertical gravity faults. Field evidence against the regional extension model includes (1) a thick cap rock at the top of the salt walls, (2) the concentration of subsidence deformation structures along the crest of the anticlines (salt walls), (3) deformational structures not consistent with the proposed NNE extension, like crestal synforms and NE–SW grabens, (4) dissolution-induced subsidence structures controlled by ring faulting, revealing deep-seated dissolution, (5) large blocks foundered several hundred meters into the salt wall, (6) evidence of recent and active dissolution subsidence, and (7) the aseismic nature of the recently active collapse faults. Although underground salt dissolution seems to be the main cause for the generation of the salt valleys, this phenomenon may have been favored by regional extension tectonics that enhance the circulation of groundwater and salt dissolution.

Paleosubsidence and active subsidence due to evaporite dissolution in Spain

Evaporite formations crop out or are at shallow depth present in an extensive area of Spain. These soluble sediments occur in diverse geological domains and were deposited over a long time span, from the Triassic up to the present day. Broadly, the Mesozoic and Paleogene formations (Alpine cycle) are affected by compressional structures, wheras the Neogene (post-orogenic) sediments remain underformed. Subsidence caused by subsurface dissolution of evaporites (subjacent karst) takes place in three main types of stratigraphic settings: a) subsidence affecting evaporite-bearing Mesozoic and Tertiary successions (interstratal karst); b) subsidence in Quaternary alluvial deposits related to the exorheic evolution of present-day fluvial systems (alluvial or mantled karst); and c) subsidence in exposed evaporites (uncovered karst). These types may be represented by paleosubsidence phenomena (synsedimentary and/or postsedimentary) recognizable in the stratigraphic record, or by equivalent, currently active or modern examples which have a surface expression. Interstratal karstification of Mesozoic marine evaporites, and the consequent subsidence of overlying strata, is revealed by stratiform collapse breccias and wedge outs of the evaporites grading into unsoluble residues. In several Tertiary basins, the sediments overlying evaporites locally show synsedimentary and/or postsedimentary subsidence structures. Dissolution-induced subsidence coeval with sedimentation is accompanied by local thicknening of strata in basin-like structures with convergent dips and cumulative wedge-out systems. This sinking process controls the generation of depositional environments and lithofacies distribution. Postsedimentary subsidence produces a great variety of gravitational deformations in Tertiary supra-evaporitic units, including both ductile and brittle structures (flexures, synforms, fractures, collapse, and brecciation). Quaternary fluvial terrace deposits overlying evaporites show anomalous thickenings (>150m) caused by a dissolution-induced subsidence process in the alluvial plain, which is balanced by alluvial aggradation. The complex evolution (in time and space) of paleosubsidence leads to intricate and chaotic structures in the alluvium, which may be erroneously interpreted as pure tectonic deformations. The current subsidence and generation of sinkholes due to suballuvial karstification constitutes a geohazard which affects large, densely populated areas, and thus endangers human safety and poses limitations on development. An outstanding example can be seen in Calatayud, an important historical city where subsidence has severely damaged highly valuable monuments. Subsidence resulting from the underground karstification of evaporites has caused or influenced the generation of some important modern lacustrine basins, such as Gallocanta, Fuente de Piedra, and Banyoles Lakes. The sudden formation of sinkholes due to collapse of cave roofs is fairly frequent in some evaporite outcrops. Very harmful and spectacular subsidence activity is currently occurring in the Cardona salt diapir, where subsidence has been dramatically exacerbated by mining practices.

Geomorphological and sedimentological features in Quaternary fluvial systems affected by solution-induced subsidence (Ebro Basin, NE-Spain)

The Quaternary evolution and the morpho-sedimentary features of some of the most important rivers in Spain (Ebro and Tagus rivers among others) have been controlled by subsidence due to alluvial karstification of the evaporitic bedrock. The subsidence mechanism may range from catastrophic collapse to slow sagging of the alluvium by passive bending. In the Ebro Basin, the mechanisms and processes involved in karstic subsidence were studied through the analysis of present-day closed depressions as well as through old subsidence depressions (palaeocollapses and solution-induced basins) and associated deformations recorded in the Quaternary alluvial sediments. The Gállego–Ebro river system is presented as a case study of channel adjustments and geomorphic and sedimentary evolution of fluvial systems in dissolution-induced subsidence areas. In this fluvial system, evaporite dissolution during particular Quaternary time intervals (namely early and middle Pleistocene) have lead to the development of a solution-induced basin, approximately 30 km-long by 8 km-wide, filled by Quaternary deposits with a total thickness in excess of 190 m. The main river response to balance the subsidence in the alluvial plain was aggradation in the central reach of the subsiding area, and degradation both in the upstream reach and in the valley sides where alluvial fans and covered pediments may prograde over the fluvial sediments. The main sinking areas are recognized in the sedimentary record by anomalous thickenings in the alluvial deposits and fine-grained sediments deposited in backswamp and ponded areas.

Subsidence beneath a playa basin on the Southern High Plains, U.S.A.: Evidence from shallow seismic data

Shallow seismic data from Sevenmile Basin, a large ephemeral lake (playa) basin in the Texas Panhandle, reveal that subsidence has been an important agent in basin formation. Several hypotheses exist for the origin of thousands of playa basins on the Southern High Plains of Texas and New Mexico, including eolian deflation, evaporite or carbonate dissolution and subsidence, piping, and animal activity. Sevenmile Basin is 5.5 x 3.6 km across and 14 m deep and contains 20 m of lacustrine and eolian sediments that interfinger with the Quaternary Blackwater Draw Formation. Below these sediments is the upper Tertiary Ogallala Formation, which overlies Permian or Triassic bedrock. Seismic reflection and refraction data were collected from the unlithified and variably saturated clastic sequence beneath Sevenmile Basin to investigate the geological history and hydrogeological framework of playa basins, which recharge the regionally important Ogallala aquifer. Three-layer velocity models provide good solutions for reversed refraction data. Near-surface p-wave velocities (layer 1) range from 349 to 505 m/s, layer 2 velocities range from 806 to 851 m/s, and layer 3 velocities range from 2,037 to 2,161 m/s. Shallow test holes and drillers' logs suggest that layer 1 is composed of playa and upper Blackwater Draw Formation deposits, layer 2 consists of lower Blackwater Draw Formation and upper Ogallala Formation deposits, and layer 3 represents a competent and partly saturated zone near the top of the Ogallala aquifer. Reflection sections show a middle Ogallala reflector, a reflector at the top of Permian or Triassic bedrock, and internal bedrock reflectors that indicate a structural low beneath Sevenmile Basin. Increasing relief with age, from 14 m at the surface to 70 m on the middle Ogallala reflector to 110 m at the base of the Ogallala, is interpreted as evidence of subsidence of underlying Permian evaporite-bearing strata before or during Ogallala deposition. Ogallala and Blackwater Draw Formation thicknesses greater than bedrock Ogallala deposition and may continue to the present. Virtually all playa basins on the Southern High Plains are underlain by Permian evaporite-bearing strata; some basins have been affected by dissolution-induced subsidence. Shallow seismic methods are an ideal approach to determine the relative importance of subsidence in basin formation.