Geological and environmental implications of the evaporite karst in Spain

Disentangling the impact of climate change, human activities, vegetation dynamics and atmospheric CO2 concentration on soil water use efficiency in global karst landscapes

Soil Water Use Efficiency (SWUE), which quantifies the carbon gain against each unit of soil moisture depletion, represents an essential ecological parameter that delineates the carbon-water coupling within terrestrial ecosystems. However, the spatiotemporal dynamics of SWUE, its sensitivity to environmental variables, and the underlying driving mechanisms across various temporal scales in the global karst region are largely uncharted. This study utilized the sensitivity algorithm of partial least squares regression, partial differential equations, and elasticity coefficients to investigate the characteristics of SWUE variations across different climatic zones in the global karst region and their responsiveness to environmental variables. Moreover, the study quantified the individual contributions of climate variability, atmospheric carbon dioxide concentration, human activities, and vegetation changes to SWUE variations. The results indicated that SWUE across different climatic zones in the global karst region demonstrated an increasing trend from 2000 to 2018, with the most notable improvement observed in the humid zone. SWUE presented regular distribution and variation characteristics across different latitudinal zones at a monthly scale. The sensitivity of SWUE to precipitation was significantly higher compared to its responsiveness to other environmental factors. Additionally, the trend in SWUE's sensitivity to precipitation demonstrated the most significant change. The sensitivity of SWUE to various environmental factors and the trend of this sensitivity in the arid zone revealed significant variation compared to other climatic zones. Gross primary productivity and soil moisture were identified as the intrinsic factors influencing SWUE changes, contributing 16 % and - 84 %, respectively. Climate variability and human activities were identified as the primary exogenous factors contributing to the increase in SWUE, accounting for 76 % and 16 %, respectively. This study advances the understanding of carbon-water coupling in karst regions, providing significant insights into the ecological management of global karst environments amidst climate variations.

Reactive Transport Model of Gypsum Karstification in Physically and Chemically Heterogeneous Fractured Media

Gypsum dissolution leads to the development of karstic features within much shorter timescales than in other sedimentary rocks, potentially leading to rapid deterioration of groundwater quality and increasing the risk of catastrophes caused by subsidence. Here, we present a 2-D reactive transport model to evaluate gypsum karstification in physically and chemically heterogeneous systems. The model considers a low-permeability rock matrix composed mainly of gypsum and a discontinuity (fracture), which acts as a preferential water pathway. Several scenarios are analyzed and simulated to investigate the relevance for gypsum karstification of: (1) the dynamic update of flow and transport parameters due to porosity changes; (2) the spatial distribution of minerals in the rock matrix; (3) the time evolution of water inflows through the boundaries of the model; (4) the functions relating permeability, k, to porosity, ϕ. The average porosity of the matrix after 1000 years of simulation increases from 0.045 to 0.29 when flow, transport, and chemical parameters and the water inflows through the boundary are dynamically updated according to the porosity changes. On the contrary, the porosity of the matrix hardly changes when the porosity feedback effect is not considered, while its average increases to 0.13 if the water inflow occurs through the discontinuity. Moreover, the dissolution of small amounts of highly soluble sulfate minerals plays a major role in the development of additional fractures. The increase in hydraulic conductivity is largest for the power law with an exponent of n = 5, as well as the Kozeny-Carman and the modified Fair-atch k-ϕ relationships. The gypsum dissolution front propagates into the matrix faster when the power law with n = 2 and 3 and the Verma–Pruess k-ϕ relationships are used.

Non-Invasive Methodological Approach to Detect and Characterize High-Risk Sinkholes in Urban Cover Evaporite Karst: Integrated Reflection Seismics, PS-InSAR, Leveling, 3D-GPR and Ancillary Data. A NE Italian Case Study

Sinkholes linked to cover evaporite karst in urban environments still represent a challenge in terms of their clear identification and mapping considering the rehash and man-made structures. In the present research, we have proposed and tested a methodology to identify the subsiding features through an integrated and non-invasive multi-scale approach combining seismic reflection, PS-InSAR (PSI), leveling and full 3D Ground Penetrating Radar (GPR), and thus overpassing the limits of each method. The analysis was conducted in a small village in the Alta Val Tagliamento Valley (Friuli Venezia Giulia region, NE Italy). Here, sinkholes have been reported for a long time as well as the hazards linked to their presence. Within past years, several houses have been demolished and at present many of them are damaged. The PSI investigation allowed the identification of an area with higher vertical velocities; seismic reflection imagined the covered karst bedrock, identifying three depocenters; leveling data presented a downward displacement comparable with PSI results; 3D GPR, applied here for the first time in the study and characterization of sinkholes, defined shallow sinking features. Combining all the obtained results with accurate field observations, we identified and mapped the highest vulnerable zone.

Assessing the effects of rainfall, groundwater downward leakage, and groundwater head differences on the development of cover-collapse and cover-suffosion sinkholes in central Florida (USA)

Cover-collapse and cover-suffosion sinkholes are widely distributed in central Florida (USA) karst terrains and have been recognized as the primary geo-hazard threatening human lives and destroying infrastructure. Previous studies indicated that the development of cover-collapse and cover-suffosion sinkholes in central Florida might be related to hydrologic/hydrogeologic conditions such as rainfall, groundwater downward leakage and groundwater hydraulic head differences (groundwater level differences between the water tables in unconfined aquifer and the potentiometric levels in confined aquifer). Here, a case study in central Florida urban areas is conducted to quantify the effects of rainfall, groundwater downward leakage and groundwater head differences on the development of cover-collapse and cover-suffosion sinkholes in central Florida with a focus on the timing of their occurrences. Results indicate that heavy rainfall/storm(s) and rapid increase of head differences within a relatively short period of time are major factors affecting the timing of sinkhole occurrences, and the spatial variation of groundwater downward leakage rate can be used to generate sinkhole susceptibility zonation maps for serving as a useful indicator of the likelihood of sinkhole development at certain areas. Results caution that the groundwater pumping and mining dewatering rate should be setup properly and the starting time of groundwater pumping and/or mining dewatering should be selected carefully in central Florida, i.e., the activities should be put into abeyance after a heavy rainfall/storm(s) to reduce the probability of sinkhole occurrences.

Stratigraphic Relations of the Ipubi Formation: Siliciclastic-Evaporitic Succession of the Araripe Basin

The Ipubi Formation represents the Aptian-Albian siliciclastic-evaporitic succession of Araripe Basin, NE Brazil. This succession comprises siliciclastic rocks (bituminous shales and claystones) and evaporites (gypsum and secondary anhydrite) and represents part of the lacustrine-shallow marine post-rift phase I. This study used sequence stratigraphy concepts to define the relations between changes in the relative lake level and the formation of Ipubi deposits. Results show that the organic-rich shales of the Ipubi Formation formed during a transgressive pulse that covered large areas of the proximal domains. These deposits overlie a regional unconformity that marks the end of the deposition of the underlying Crato Formation. A High Stand stage that followed the transgression influenced the formation of evaporitic deposits. Climate conditions played a major role in influencing the triggering and stopping of evaporite deposition. Thus, a new relative lake level fall event caused the exposure of the Ipubi Formation deposits, and created another regional subaerial unconformity accompanied by widespread karstification of evaporite beds. A posterior transgression caused the deposition of siliciclastic rocks of the Romualdo Formation over the Ipubi Formation strata, and also promoted a new event of karstification of the Ipubi upper evaporite beds.

Karst tiankengs as refugia for indigenous tree flora amidst a degraded landscape in southwestern China

We conducted floristic and community analyses to compare the floristic composition, forest structure, taxonomic richness, and species diversity between two tiankeng (large doline, or sinkhole) habitats and two outside-tiankeng habitats of forest fragments in a degraded karst area in southwestern China. We found remarkably higher taxonomic richness in the tiankeng habitats than in the outside-tiankeng habitats at the species, generic, and familial levels. The inside-tiankeng habitats had higher floristic diversity but lower dominance. The remarkably higher uniqueness at all taxonomic levels and the much larger tree size in the two tiankeng habitats than in the outside-tiankeng habitats demonstrated the old-growth and isolated nature of the tiankeng flora. Plot-scale species richness, Shannon-Wiener index, Pielou’s evenness, and Berger-Parker dominance significantly differed across habitats. Heterogeneity in floristic composition at the species, generic, and familial levels was extremely significant across habitats. In pairwise comparisons, except for the Chuandong Tiankeng-Shenmu Tiankeng pair, all the pairs showed significant between-habitat heterogeneity in floristic composition. Our results suggest that as oases amidst the degraded karst landscape, tiankengs serve as modern refugia that preserve old-growth forest communities with their rich floristic diversity, and can provide a model for habitat conservation and forest restoration in that area.

Salt lakes of La Mancha (Central Spain): A hot spot for tiger beetle (Carabidae, Cicindelinae) species diversity

The tiger beetle assemblage of the wetlands of La Mancha (central Spain) comprises nine species: Calomera littoralis littoralis, Cephalota maura maura, Cephalota circumdata imperialis, Cephalota dulcinea, Cicindela campestris campestris, Cicindela maroccana, Cylindera paludosa, Lophyra flexuosa flexuosa, and Myriochila melancholica melancholica. This assemblage represents the largest concentration of tiger beetles in a single 1º latitude / longitude square in Europe. General patterns of spatial and temporal segregation among species are discussed based on observations of 1462 specimens registered during an observation period of one year, from April to August. The different species of Cicindelini appear to be distributed over space and time, with little overlapping among them. Three sets of species replace each other phenologically as the season goes on. Most of the species occupy drying or dried salt lakes and salt marshes, with sparse vegetation cover. Spatial segregation is marked in terms of substrate and vegetation use. Calomera littoralis and Myriochila melancholica have been observed mainly on wet soils; Cephalota circumdata on dry open saline flats; Cephalota dulcinea and Cylindera paludosa in granulated substrates with typical halophytic vegetation; Cephalota maura is often present in man-modified areas. Cephalota circumdata and Cephalota dulcinea are included as species of special interest in the list of protected species in Castilla-La Mancha. Conservation problems for the Cicindelini assemblage arise from agricultural activities and inadequate use of sport vehicles. Attempts at restoring the original habitat, supressing old semi-industrial structures, may affect the spatial heterogeneity of the lakes, and have an effect on Cicindelinae diversity.

Assessment of groundwater quality of the Tatlicay aquifer and relation to the adjacent evaporitic formations (Cankiri, Turkey)

One of the most important hydrogeologic problems in and adjacent areas of evaporitic formations is severe quality degradation of groundwaters. These kinds of groundwaters contain high content of dissolved solids and generally have some limitations for use. Tatlicay basin (north-central Turkey) is an example to effects of the evaporites on groundwater quality in the adjacent alluvium aquifer. Gypsum and anhydrites in the two evaporite formations (Bayindir and Bozkir) effect of the groundwater quality in the alluvium adversely, by dissolution of the evaporites by surface drainage and infiltration into the alluvium aquifer (widespread effect) and by infiltration of low quality gypsum springs (local effect) into the aquifer. Evaporitic formations significantly increased EC, TDS, Ca and SO(4) parameters in the alluvium aquifer in the central and downstream regions. EC has increased roughly from 500-800 to 1,700-2,000 μS/cm, Ca has roughly increased from 3-4 to 10 meq/l, SO(4) has increased 0.5-1 to 11-12 meq/l. Consequently, three clusters were distinguished in the basin; (1) nonevaporitic waters in low TDS, Na, Ca, Mg, Cl and SO(4), (2) diluted waters in high TDS and relatively high Cl, moderate-relatively high Na, Ca, Mg, SO(4), (3) gypsum springs in highest TDS, Ca, SO(4), but moderate Mg and low Na, Cl.