Conduit Network Research Articles

Influence of structural properties and connectivity of initial fracture network on incipient karst genesis

A karst system is generally difficult to characterize especially regarding its conduits pattern and morphology. Even though speleology allows to map the geometry of cave systems, a large part of karst conduit networks remains unreachable underground. In this study, we use a reactive transport model that simulates dissolution processes in fracture networks to investigate the effect of discrete fracture network (DFN) properties (i.e., anisotropy, connectivity, fracture intersection types) on incipient karst conduit geometry under different boundary conditions (constant hydraulic head versus constant flow rate conditions, directional versus concentrated recharge conditions). We focused on single fracture (with/without cementation) and ladder like fracture networks which are common on limestone. Results show that fracture network properties, play a major role in the final geometry and topology of the incipient karst conduit patterns. Also, the presence of cement on a single fracture favors tortuous and curvilinear wormholes, while the anisotropy of the fracture network favors the formation of anisotropic incipient karst geometries. Finally, the type of hydraulic BC remains a major factor that governs dissolution processes and thus incipient karst conduit geometry.

Hydrochemical characterization of a Dinaric karst catchment in relation to emerging organic contaminants

The main findings of a hydrochemical investigation conducted within a typical Dinaric karst catchment located in Southern Croatia are outlined. The studied aquifer is drained by the Jadro and Žrnovnica springs, which are important for the regional and local water supplies, respectively. Presumably, there is intercatchment groundwater flow coming from the neighbouring Cetina River catchment. Various factors governing aquifer hydrochemistry and their interplay with emerging organic contaminants (EOCs) that were detected at different water resources in ng/L concentrations was assessed. A total of 26 sampling campaigns (October 2019 – October 2022) were conducted at two springs, in a river and at a deep borehole, all representative of this complex hydrogeological system. Assessment of major ion constituents and saturation indices calculated with PHREEQC revealed the sampled water resources are of a Ca-HCO3 type due to the predominant weathering of the carbonate mineral calcite. Sharp spikes observed in chemographs indicated a highly karstified system with an effective conduit network allowing rapid spring responses to precipitation events. Water resources are of good chemical status, as affirmed by anthropogenic contamination indicators, with nitrates, chlorides and sulphates all below maximum threshold values. Strong positive correlations were found between EOCs concentration, number of detected compounds, and nitrates in the Cetina River, indicating a common origin, most likely wastewater. Identification of persistent EOCs including widely used repellent N,N-diethyl- metatoluamide (DEET) during base flow conditions and its strong positive correlation with the Ca2+ content in both the Cetina and Jadro samples, suggests potential storage in the epikarst and aquifer matrix. This coupling of conventional hydrochemical indicators and novel markers of anthropogenic impacts, including EOCs, in vulnerable karst water resources is a crucial advancement in the assessment and management of emerging environmental and potential human health risks. Such an approach is pivotal for the sustainable protection of hydrogeologically intricate sites.

Connecting past, present, and future trends of hydraulic and phosphorus loading in the Bay of Quinte tributaries, Ontario, Canada

Study regionThe Bay of Quinte watershed, located on the northeastern shore of Lake Ontario, Canada. Special focus is placed on the Napanee River and Wilton Creek catchments, where the presence of a subterranean network of conduits that facilitate groundwater-surface water interactions is conducive to an increase in the delivery of nutrients during storm events and can exacerbate the risk of eutrophication. Study focusDevelopment of a rigorous modelling framework that can advance our understanding of past and present trends of streamflow and tributary phosphorus (P) concentration and loads. We also analyze future streamflow rates and P export trends in order to guide the long-term watershed management in the area. Our analysis is carried out through a comprehensive combination of data-driven and process-based models. Novel hydrological insights for the regionRetrospective analysis provides evidence of a recent increase in the impact of nonpoint-source P pollution. General Circulation Models suggest a future increase in both temperature minima and maxima relative to present conditions (2002–2018). Precipitation projections are indicative of increased frequency of occurrence of high extreme precipitation events during the summer and mid-fall, when the Bay of Quinte is more susceptible to undesirable ecological shifts. Our study predicts a seasonality shift in the streamflow, with the spring freshet occurring earlier in the year accompanied by higher flow rates in the local creeks for the winter and early spring. We also examined the relationship between precipitation and streamflow across varying levels of flow regulation. Our analysis indicates that our predictive capacity declines with increased flow regulation as well as when streamflow rates are predicted by contemporaneous precipitation values.

The Trans-Crustal Evolution of Miocene Magmatism Parental to Epithermal Mineralization in the Sierra Madre del Sur, Southern Mexico

Abstract Arc-related volcanism results from the interplay of magmatic processes occurring in trans-crustal systems that consolidate over time in a given location. Trans-crustal systems comprise extensive networks of magma reservoirs and conduits at different depths, wherein melts cyclically differentiate and segregate before reaching the surface. The study of these systems provides valuable insights into crustal-scale phenomena, such as the evolution of the continental crust and the formation of metallogenic regions. In this study, we address the trans-crustal magmatic evolution of lower Miocene magmatism parental to several intermediate-sulfidation epithermal deposits in the eastern Sierra Madre del Sur igneous province, southern Mexico. Using a multi-methodological approach, we document changes over ca. 1.1 Myr in the magmatic system that fed andesitic-to-felsic volcanism in this region. We employ whole-rock REE ratios and λ parameters—which are used to quantify the shape of a REE pattern—to track the involvement of pressure-sensitive minerals in the deep-crustal magmatic evolution. The andesitic rocks consist of lava flows, porphyries, and dikes that collectively show REE patterns suggestive of control by fractionated or residual (i.e. in crustal melting) amphibole and/or clinopyroxene. In contrast, the felsic rocks consist of rhyolitic–dacitic ignimbrites, domes, and dikes that show contrasting REE patterns suggestive of control by plagioclase, clinopyroxene, amphibole, and/or garnet. The distinct pressure-sensitive mineral assemblages in the andesitic and felsic rocks indicate that the locus of deep-crustal magma evolution varied within the middle–lower crust. These magmas were episodically injected into ephemeral shallow crustal reservoirs (shortly?) before being erupted, inducing a progressive thermomechanical maturation of the middle–upper crust. Meanwhile, low degrees of crustal assimilation occurred as recorded by Mesozoic inherited zircon ages and Sr–Nd–Pb radiogenic isotopes. An extensive middle–lower crustal magma evolution has been linked to the formation of porphyry Cu deposits (i.e. ‘fertile’ magmatism). Given that intermediate-sulfidation epithermal deposits may be genetically linked with porphyry Cu deposits, the documented processes could have contributed to the formation of epithermal deposits in the region. However, magmatic fertility proxies resemble those from infertile magmas worldwide. Since these proxies have been exclusively applied to porphyry-type deposits, our results highlight the importance of developing new geochemical exploration tools applicable to a wider range of ore deposits.

From Recharge to Cave to Spring: Transmission of a Flood Pulse through a Complex Karst Conduit Network, Castleton, Derbyshire (UK)

Storm Babet (18–21 October 2023) brought heavy and persistent rain (80–100 mm) to the English Peak District, causing widespread surface and underground flooding. The village of Castleton experienced groundwater flooding from springs that drain a complex mixed allogenic–autogenic karst catchment. Transmission of the flood pulse was monitored using high-resolution (2 and 4 min intervals) logging of (a) the hydraulic head at five underground locations in the karst conduits and (b) the water depth at three springs and in the surface river fed by the springs. Underground, there were large increases in the hydraulic head (9–35 m), which resulted in two types of flow switching. Firstly, the increased head at the input end of a phreatic (water-filled) conduit system removed an underwater permeability barrier in a relatively low-elevation conduit, resulting in a dramatic increase in flow out of the conduit and a corresponding decrease in flow from a linked higher-elevation conduit that had dominated before the storm. Secondly, the increased head upstream of two conduits with limited hydraulic conductivity allowed water to spill over into conduits that were inactive prior to the storm. As expected, the conduits fed by sinking streams from the allogenic catchment responded rapidly to the recharge, but there was also a rapid response from the autogenic catchment where there are no surface streams and only a small number of dolines. The complex signals measured underground are not apparent from the spring hydrographs.

Groundwater flow modelling in Neoproterozoic carbonate karst based on dye tracer pathways

The current study consisted of the development of a hydrogeological numerical model in the Lagoa Santa Karst Environmental Protection region, in order to represent the conduit-flow karst springs in the most well-studied karst system in Brazil. Once carbonate aquifers are characterized by a complex 3D pattern with a combination of a double porosity of matrix and conduits, two different modelling approaches were used: Equivalent Porous Medium (EPM), based on Darcy’s law, to represent laminar flow in the fractured massif and Conduit Network (CN), using Manning-Stricker Law and the discrete feature tool, to represent the karst conduits with turbulent flow. Modeling was calibrated in a steady-state and simulations of groundwater flow and an evaluation of water exchanges between the karstic aquifer and the other hydrogeological units, were conducted. The model was developed based on data available for many studies carried out in the area, including tracer-derived route data. The calibration results for spring discharges, water levels and mass balance were satisfactory and aligned with modeling good practices. In turn, the hydraulic conductivities and recharge rates obtained by back-analysis were in agreement with literature references. The regional groundwater flow was well represented and indicated how external water entrances can be relevant to explain the higher discharge rates observed in some karstic springs. In addition, the representation of conduits allowed for the simulation of tracer pathways, enabling the inference of potential flow rates to the karstic springs from the particle simulations.

Natural responses of Neoproterozoic dynamic karst springs to rainfall events, São Miguel Watershed, Minas Gerais, Brazil

AbstractKarst aquifers consist of complex networks of conduits in which groundwater flows and recharge/discharge processes are generally more dynamic than in other types of aquifers. Due to their intrinsic heterogeneity and anisotropy, monitoring, quantifying, and analysing natural responses of karst springs is an efficient tool. Unlike Cenozoic and Mesozoic rocks, in Neoproterozoic karst systems, groundwater circulates and stores generally in dissolution features known as tertiary porosity, as the rock's primary porosity is recrystallized, considered negligible. This article studies the hydrodynamics of a karst portion of the São Miguel River basin, southwest of the state of Minas Gerais, Brazil. The region is predominantly composed of Neoproterozoic carbonate rocks, dating from about 570 to 540 million years ago. During a hydrological year (2019–2020), three karst springs (S1, S2, and S3) were daily monitored through their natural responses (variations of electrical conductivity, EC, temperature, T, and discharge, Q) to rainfall episodes. The data were interpreted based on the analysis of spring hydrographs, time series, recession curves (seasonal and intra‐annual), and statistics of EC, T, and Q variations. The results show the three springs generally exhibit quick flow, typically karstic, in the case of hydrosystems with a well‐structured and functional underground drainage network. The time series indicate the hydrosystem drained by S1 presents slower circulation and a lower degree of linearity, resulting from the higher sinuosity of the system, while the hydrosystems of S2 and S3 have similar behaviours, of quick water circulations immediately after a rainy episode. The degrees of karstification classify S1 and S2 as complex and extensive karst systems consisting of several subsystems, and S3 as a system in which the conduit network is more developed at the upper epiphreatic zone than near the outlet.

Expanding Karst Groundwater Tracing Techniques: Incorporating Population Genetic and Isotopic Data to Enhance Flow-Path Characterization

Karst aquifers are unique among groundwater systems because of variable permeability and flow-path organization changes resulting from dissolution processes. Over time, changes in flow-path connectivity complicate interpretations of conduit network evolution in karst hydrogeology. Natural and artificial tracer techniques have long provided critical information for protecting karst aquifers and understanding the potential impacts on ecosystems and human populations. Conventional tracer methods are useful in karst hydrogeologic studies for delineating flow paths and defining recharge, storage, and discharge properties. However, these methods only provide snapshots of the current conditions and do not provide sufficient information to understand the changes in interconnection or larger-scale evolution of flow paths in the aquifer over time. With advances in population genetics, it is possible to assess population connectivity, which may provide greater insights into complex groundwater flow paths. To assess this potential, we combined the more traditional approaches collected in this and associated studies, including artificial (dye) and natural (geochemistry, isotopes, and discharge) tracers, with the population genetic data of a groundwater crustacean to determine whether these data can provide insights into seasonal or longer changes in connections between conduits. The data collected included dye trace, hydrographs, geochemistry, and asellid isopod (Caecidotea bicrenenta) population genetics in Fern Cave, AL, USA, a 25 km-long cave system. Combined, these data show the connections between two separate flow paths during flood events as the downstream populations of isopods belonging to the same subpopulation were measured in both systems. Additionally, the sub-populations found in higher elevations of the cave suggest a highly interconnected unsaturated zone that allows for genetic movement in the vadose zone. Although upstream populations show some similarities in genetics, hydrologic barriers, in the form of large waterfalls, likely separate populations within the same stream.

Evaluating physical controls on conduit flow contribution to spring discharge

The physical controls on the contribution of conduit to spring discharge, quantified as the ratio of conduit flow to spring discharge (Qc/Q), are evaluated through numerical simulations using MODFLOW-Conduit Flow Process. The controlling factors include the ratio of conduit surface roughness to conduit diameter (ε/D), the hydraulic conductivity ratio between the conduit wall and the porous matrix (Kc/K), the ratio of groundwater recharge to the hydraulic conductivity of the porous matrix, and conduit density (Cd) defined as the ratio of total conduit length to springshed area. For a single conduit, Qc/Q increases with Kc/K when flow in the conduit is laminar. When flow is turbulent, Qc/Q is stable and then decreases with ε/D, increases with Kc/K, and decreases with R/K due to the limited exchange capacity between matrix and conduit (i.e., small Kc/K). The results of hypothetical conduit networks show that Qc/Q increases with Cd but is not sensitive to the two sets of network structure evaluated. The sensitivity analysis for the Silver Springs in Florida shows that the spatial heterogeneity of conduit diameter causes higher Qc/Q, and Qc/Q is most sensitive to the conduit density, followed by R/K, and then Kc/K.

Hu-Qi-Zheng-Xiao Decoction Inhibits the Metastasis of Hepatocellular Carcinoma Cells by Suppressing the HIF-1α Signaling Pathway to Inhibit EMT, LCSC, and Angiogenic Process.

Hepatocellular carcinoma (HCC) is a common clinical malignant tumor of the digestive system. Hu-Qi-Zheng-Xiao (HQZX) decoction has been clinically found to prolong the survival of patients with hepatocellular carcinoma and improve the quality of patients' survival, but its antitumor biological mechanism is still unclear. A nude mouse hollow fiber hepatocellular carcinoma model was constructed to analyze the in vivo efficacy of HQZX decoction against 7 different hepatocellular carcinoma cells. The subcutaneous graft tumor model was again validated. In vitro, the effect of HQZX decoction on the growth and metastasis of the cell line with the highest growth inhibition was evaluated. The cell line with the best efficacy response screened was again used to construct a hollow fiber hepatocellular carcinoma model and hollow fiber conduit cells were extracted to detect the expression of HIF-1α, VEGF, EMT-related molecules, LCSCs-related molecules, and to observe the density of the subcutaneous vascular network of hollow fiber conduits. The liver metastasis model of splenic injection was constructed to observe the effect of HQZX decoction on tumor metastasis. The hollow fiber hepatocellular carcinoma model was evaluated for the efficacy of HQZX decoction, and it was found to have the highest growth inhibition of LM3-luc cells. In vitro, the CCK8 assay revealed that HQZX decoction could inhibit tumor migration and invasion and promote apoptosis. In addition, the mechanism study of extracting cells from hollow fiber tubes found that HQZX decoction could inhibit metastasis-associated HIF-1α, VEGF, EMT-related molecules, and LCSCs-related molecules expression. capillary network around subcutaneous fiber tubes was reduced in the HQZX decoction gavage group of mice. It inhibited tumor metastasis in nude mice. HQZX decoction inhibited the growth of a variety of hepatocellular carcinoma cells. HQZX decoction suppressed the expression of metastasis-associated VEGF, EMT-related molecules, and LCSCs-related molecules and inhibited tumor angiogenesis and growth and metastasis, which may be related to the inhibition of the HIF-1α signaling pathway. It reveals that HQZX decoction may be a promising herbal compound for anti-HCC therapy, and also reveals the accurate feasibility of the hollow fiber hepatocellular carcinoma model for in vivo pharmacodynamic evaluation and mechanism study.

Stochastic simulation and uncertainty analysis of karst conduit network using pyKasso

Owing to the intrinsic heterogeneity and anisotropy of karst systems, traditional hydrological exploration methods face significant challenges when investigating karst conduit networks. This study employs pyKasso for the stochastic simulation of karst conduit networks in Panzhou City, focusing on uncertainty analysis through local sensitivity analysis and Monte Carlo methods. The simulation process incorporates geological, topographic, and fracture data to create a realistic representation of the karst network. We found that the spatial configuration and characteristics of the karst network are significantly influenced by various input parameters such as fracture parameters, inlets, outlets, and cost ratios. We highlight the minimal influence of fracture densities and the substantial impact of the count of inlets and outlets on crucial network metrics. The results demonstrate the critical role of parameter sensitivity and variability in modeling the intricate karst systems, providing valuable insights for hydrogeological studies and the management of karst water.

Effects of global and climate change on the freshwater-seawater interface movement in a Mediterranean karst aquifer of Mallorca Island

Karst aquifers are globally prized freshwater sources, posing a significant preservation challenge. These aquifers typically exhibit dual or even triple porosities, encompassing matrix, fractures-fissures and conduits, rendering them highly responsive to variations in chemical characteristics and hydraulic head. In coastal regions, these aquifers often possess extensive subsurface conduit networks intricately linked to the rock matrix, facilitating groundwater discharge into the sea. Therefore, they display acute sensitivity to seawater intrusion, swiftly reacting to changes in precipitation and pumping regimes. This makes them exceptionally vulnerable to short-term meteorological fluctuations and long-term climate change. Their high heterogeneity leads to uneven penetration of the freshwater-seawater interface, causing rapid seawater intrusion inland over significant distances. The Mediterranean region, characterized by water deficit and water stress, faces strong impacts from climate change, featuring a warming atmospheric trend exceeding the global average, along with diminished rainfall exacerbating water scarcity. Increasing water demands for agriculture, urban development, and the growing tourism industry, because of global change, are worsening water stress. Our primary research objectives were analyzing the environmental consequences of global and climate change on seawater intrusion in Mediterranean coastal karst aquifers, with a focus on the role of the double-flow model, thus contributing to the understanding of the processes involved. To achieve this, we selected a study region on Mallorca Island in the western Mediterranean, where a karst aquifer system discharges into the sea. We employed various study methods, notably hydrochemical techniques and multi-isotopic analysis, encompassing the examination of 2H and 18O isotopes in water, 87Sr/86Sr ratio, Sr and B concentrations, and δ11B in water. A key finding is the rebound effect, wherein aquifers recontaminate due to solute molecular back-diffusion following cessation of extractions and the retreat of marine intrusion, providing insight into the impact of climate and global change on Mediterranean karst aquifers.

Influence of sediments burying the discharge area of a karst aquifer on the groundwater flow field—Numerical testing of conceptual models

AbstractKarst springs are a natural result of karst water discharging to the surface through unimpeded pathways where the water table meets the surface. This study investigates the impact of alluvial deposits of varying thicknesses and permeabilities burying the main outlet (karst spring) of a well‐developed conduit network on karst drainage, including the development of hydraulic heads, drainage patterns and conduit‐matrix interactions in response to a positive base‐level shift. Numerical testing using FEFLOW on a simplified conceptual model of a hypothetical karst aquifer with six different model configurations was used to examine various drainage structures (with and without flow through a conduit), spring conditions (free vs. partially/fully clogged), sediment cover thickness (20 and 50 m), and hydraulic conductivity of the sediments (low and high). The numerical testing model incorporated one‐dimensional discrete feature elements to simulate conduit flow and coupled conduit‐matrix interactions. Results indicate that even with a fully plugged outlet, the conduit network remains a significant contributor to the drainage system, collecting water from the matrix in the recharge zone. As the outlet becomes buried, the hydraulic head increases along the conduit, forcing water back up into the matrix. The elevated hydraulic head in the karst system will cause new conduits to form at the contact between limestone and sediments, creating new potential spring sites (or reactivating existing paleo‐phreatic levels). Artesian conditions will occur below the low permeability sediments. These findings provide valuable insights into the responses of natural karst systems.

Improving understanding of groundwater flow in an alpine karst system by reconstructing its geologic history using conduit network model ensembles

Abstract. Reconstructing the geologic history of a karst area can advance understanding of the system's present-day hydrogeologic functioning and help predict the location of unexplored conduits. This study tests competing hypotheses describing past conditions controlling cave formation in an alpine karst catchment, by comparing an ensemble of modeled networks to the observed network map. The catchment, the Gottesacker karst system (Germany and Austria), is drained by three major springs and a paleo-spring and includes the partially explored Hölloch cave, which consists of an active section whose formation is well-understood and an inactive section whose formation is the subject of debate. Two hypotheses for the formation of the inactive section are the following: (1) glaciation obscured the three present-day springs, leaving only the paleo-spring, or (2) the lowest of the three major springs (Sägebach) is comparatively young, so its subcatchment previously drained to the paleo-spring. These hypotheses were tested using the pyKasso Python library (built on anisotropic fast-marching methods) to generate two ensembles of networks, one representing each scenario. Each ensemble was then compared to the known cave map. The simulated networks generated under hypothesis 2 match the observed cave map more closely than those generated under hypothesis 1. This supports the conclusion that the Sägebach spring is young, and it suggests that the cave likely continues southwards. Finally, this study extends the applicability of model ensemble methods from situations where the geologic setting is known but the network is unknown to situations where the network is known but the geologic evolution is not.

Microbial biogeography of the eastern Yucatán carbonate aquifer.

The extensive Yucatán carbonate aquifer, located primarily in southeastern Mexico, is pockmarked by numerous sinkholes (cenotes) that lead to a complex web of underwater caves. The aquifer hosts a diverse yet understudied microbiome throughout its highly stratified water column, which is marked by a meteoric lens floating on intruding seawater owing to the coastal proximity and high permeability of the Yucatán carbonate platform. Here, we present a biogeographic survey of bacterial and archaeal communities from the eastern Yucatán aquifer. We apply a novel network analysis software that models ecological niche space from microbial taxonomic abundance data. Our analysis reveals that the aquifer community is composed of several distinct niches that follow broader regional and hydrological patterns. This work lays the groundwork for future investigations to characterize the biogeochemical potential of the entire aquifer with other systems biology approaches.

Mineral chemistry of the Água Limpa suite: Insights into petrological constraints and magma ascent rate of Mesoarchean sanukitoids from the Sapucaia terrane (Carajás Province, southeastern Amazonian Craton, Brazil)

The Água Limpa suite is comprised of two E‒W elongated sanukitoid-affinity plutons named the Água Azul and Água Limpa, exposed along with other intrusive granitic bodies and typifying Mesoarchean granitic activity in the Sapucaia terrane (Carajás Province). In the present study, the T-P-fO2-H2Omelt conditions, viscosity of magma and upward magma transport rate were retrieved on the basis of mineral chemistry and partial dissolution of magmatic epidote. The mafic mineral assemblage is dominated by amphibole (magnesio-hornblende, tremolite, pargasite and magnesiohastingsite compositions), in addition to magmatic epidote (Ps of 25–32%) and magnesio-biotite. Titanite, apatite, zircon, magnetite, ilmenite and chalcopyrite are the main accessory phases. The apatite saturation thermometer indicates at least initial crystallization temperatures between 903 and 850 °C. The Al-in-amphibole geobarometer provided temperatures varying from 1039 to 718 °C and 984 to 707 °C for the Água Limpa and Água Azul plutons, respectively. A variation in pressure from 625 to ∼100 MPa was obtained for both plutons, related to the deep crustal level of formation and shallow depths of emplacement. Experimental studies carried out on similar magma compositions indicate the parental Água Limpa suite magmas were water-rich (> 6 wt% dissolved H2O) and crystallized under high oxygen fugacity conditions (0.3 ≤ NNO ≤ 2.5). Partial corrosion of the magmatic epidote in these two plutons occurred within a timeframe of no more than 1–9 years, corresponding to upward magma transport rates of 1481–13,333 m/year. A structurally controlled network of conduits associated at the mesoscale pervasive migration developed during the active deformation stage (faults and shear zones), and small-scale diking probably acted as pathways for magma transfer through distinct pulses that underwent subsequent crystallization before reaching the hot continental crust at peak metamorphism, and permitted survival of the epidote that was out of equilibrium at the low pressure of final emplacement. The Água Limpa suite reveals strong mineralogical and geochemical affinities with the sanukitoids from the Rio Maria suite and the Ourilândia complex of the Carajás Province and differs from the sanukitoids from other cratons in terms of age and geological setting.

The Quaternary activity of the Estella diapir from the uplift record of fluvial terraces, pediments and cave sediments in the Western Pyrenees, Spain

This work analyses the Quaternary landforms, geomorphic changes and stratigraphic evidences of evaporite dissolution and halokinesis in the Estella salt diapir (Basque-Cantabrian Basin). The mapping of Quaternary deposits demonstrates that karstic subsidence and diapirism happen at the same time. The diapir geometry is characterised by a prominent 4.5 km-long fault scarp attributed to the uplift of the diapir and a fluvial gypsum escarpment along the northern and eastern diapir edges. The dissolution of the gypsiferous caprock is responsible for the development of a large number of collapse sinkholes, the local thickening of terrace deposits and a network of endokarstic conduits connected to saline springs that cause a notable increase in the salinity of the Ega river across the diapir (70% electrical conductivity increase). The 40 m vertical displacement of T4 terrace and P5 pediment deposits, the anomalous longitudinal profile of the old terraces, the deflection of infilled-valleys at the base of the fault scarp, the occurrence of wind gaps and the continuous lateral migration of the Ega river away from the diapir support salt upwelling since 644 ± 81 ka ago at a long-term, average uplift rate of 0.05–0.07 mm/yr. The 12 m offset of T8 terrace, the steepening of the Ega river floodplain and the tilting of the upper level and the antigravitative erosion genesis of the Longinos cave point to a faster diapir rising rate of 0.23–0.57 mm/yr in the last 36 ± 15 ka coinciding with a period of increasing river incision. The lack of tectonic activity in the region suggests that diapirism is enhanced by erosional unloading related to fluvial entrenchment. The thickening of pediment P5 supports the north-northeast flow of salt towards the river valley.

Frontiers of thermobarometry: GAIA, a novel Deep Learning-based tool for volcano plumbing systems

The anatomy of the plumbing system of active volcanoes is fundamental to understand how magma is stored and channeled to the surface. Reliable geothermobarometric estimates are, therefore, critical to assess the depths and temperatures of the complex system of magmatic reservoirs that form a volcano apparatus. Here, we developed a novel Machine Learning approach (named GAIA, Geo Artificial Intelligence thermobArometry) based upon Feedforward Neural Networks to estimate P-T conditions of magma (clinopyroxene) storage and migration within the crust. Our Feedforward Neural Network method applied to clinopyroxene compositions yields better uncertainties (Root-Mean-Square Error and R2 score) than previous Machine Learning methods and set the basis for a novel generation of reliable geothermobarometers, which extends beyond the paradigm associated to crystal-liquid equilibrium. Also, the bootstrap procedure, inherent to the Feedforward Neural Network architecture, permits to perform a rigorous assessment of the P-T uncertainty associated to each clinopyroxene composition, as opposed to the Root-Mean-Square Error representing the P-T uncertainty of whole set of clinopyroxene compositions.As a test, we applied GAIA to assess P-T conditions of five Italian volcanoes (Somma-Vesuvius, Campi Flegrei, Etna, Stromboli, Volcano), which are among the most dangerous volcanic centres in Europe. The results on the depths of the plumbing systems are in excellent agreement with those obtained with independent geophysical and geodetic surveys, and provide further evidence to current models of volcano plumbing systems consisting of physically-separated reservoirs interconnected by a network of conduits channelling magma en route to the surface. The results on the magma (clinopyroxene crystallization) temperatures are also in agreement with other estimates, albeit obtained considering - mainly but not only - thermodynamically-based clinopyroxene-liquid geothermometers.GAIA can set robust estimates of magma storage, segregation, and ascent conditions within the plumbing system of active volcanoes, helping to unravel P-T variations, if any, during their eruptive history and providing robust clues to volcanic hazard assessment.

Advances, challenges and perspective in modelling the functioning of karst systems: a review

We present a discussion of the state-of-the art on modelling geometrical characteristics, hydrogeological behavior and geochemical evolution of karst aquifers of meteoric origin. The considered key karst processes include: (1) the recharge processes, (2) the various hydrodynamic processes in the vadose and phreatic zones, (3) the related transport processes and (4) the speleogenesis processes. Different types of approaches for modelling geometrical characteristics of karst conduit networks are summarized. Integrated numerical studies on hydrogeological behavior of karst aquifers with functional and physically based models are then reviewed. Challenging issues in karst modelling are further discussed based on previous modelling progresses. The paper finally offers recommendations for advancing the modelling of hydrogeological behavior of karst systems and lists several open questions for future research.

Fire, frost, and drought constrain the structural diversity of wood within southern African Erica (Ericaceae)

Abstract Erica comprises ~860 species of evergreen shrubs and trees ranged from Europe to southern Africa and Madagascar. Wood structure of the around 20 European species is well studied, but despite its relevance to adaptation across the wider geographic range, it has not yet been explored across the much greater diversity, particularly of southern African lineages. In this study, we examine wood structure of 28 Erica species from southern Africa. In the African Erica clade, loss of scalariform perforation plates could be driven by increased aridity and seasonality in the mid-Miocene, and its re-gain can represent an adaptation to freezing in the high elevation species E. nubigena. As vessels in Erica are mostly solitary, imperforate tracheary elements probably form a subsidiary conduit network instead of vessel groups. Increase of ray frequency in habitats with a prominent dry and hot season probably facilitates refilling of vessels after embolism caused by water stress. Wider rays are ancestral for the lineage comprising African Erica and the Mediterranean E. australis. The negative correlation between ray width and expression of summer drought is consistent with Ojeda’s model explaining the diversification of seeders and resprouters among southern African Erica.