Karst Aquifer Research Articles

Hydrotectonics of Grand Canyon Groundwater

The Grand Canyon provides a deeply dissected view of the aquifers of the Colorado Plateau and its public and tribal lands. Stacked sandstone and karst aquifers are vertically connected by a network of faults and breccia pipes creating a complex groundwater network. Hydrochemical variations define structurally controlled groundwater sub-basins, each with main discharging springs. North Rim (N-Rim), South Rim (S-Rim), and far-west springs have different stable isotope fingerprints, reflecting different mean recharge elevations. Variation within each region reflects proportions of fast/slow aquifer pathways. Often considered perched, the upper Coconino (C) aquifer has a similar compositional range as the regional Redwall-Muav (R-M) karst aquifer, indicating connectivity. Natural and anthropogenic tracers show that recharge can travel 2 km vertically and tens of kilometers laterally in days to months via fracture conduits to mix with older karst baseflow. Six decades of piping N-Rim water to S-Rim Village and infiltration of effluent along the Bright Angel fault have sustained S-Rim groundwaters and likely induced S-Rim microseismicity. Sustainable groundwater management and uranium mining threats require better monitoring and application of hydrotectonic concepts. ▪Hydrotectonic concepts include distinct structural sub-basins, fault fast conduits, confined aquifers, karst aquifers, upwelling geothermal fluids, and induced seismicity.▪N-Rim, S-Rim, and far-west springs have different stable isotope fingerprints reflecting different mean recharge elevations and residence times.▪The upper C and lower R-M aquifers have overlapping stable isotope fingerprints in a given region, indicating vertical connectivity between aquifers.▪S-Rim springs and groundwater wells are being sustained by ∼60 years of piping of N-Rim water to S-Rim, possibly inducing seismicity.

Assessing the impact of climate change on a coastal karst aquifer in a semi-arid area

This study investigates the potential future state of a coastal karst aquifer located in the Salento peninsula of southern Italy, examining the impact of climate change. After a historical and future overview of the climate in the study area, an analysis of meteorological indices (SPIs and SPEIs) and groundwater levels (GWLs) during the historical period revealed a strong correlation between them. Then, straightforward linear regression models were applied to establish relationships between SPIs/SPEIs and GWLs for the analysed wells. To estimate future SPIs and SPEIs, an ensemble of regional climate models was used, and estimates of GWLs until the end of the century were obtained by incorporating these projections into the regression models. Two emission pathways, the RCP4.5 and the RCP8.5, were considered in the analysis. The results based on SPIs suggest no significant changes in future GWLs across all wells. However, when considering also temperature through the use of SPEIs, a general decline in groundwater levels is expected, with the decrease being more pronounced under the RCP8.5.

Large-scale spatial reconstitution of groundwater fluxes in a karst aquifer on the basis of hydrodynamic and hydrodispersive measurements

This article presents a multiphysics-based approach for a large-scale spatial reconstitution of head and tracer tests responses measured in a regional scale karst aquifer in Southern France. The main dataset consists in hydraulic heads measured in 11 wells during 6 weeks and 4 tracer tests. A pump station withdraws groundwater from this aquifer at varying flow rates during days and nights to supply water to Montpellier agglomeration, generating a periodic signal which propagates in the whole aquifer. Some of the measurements show a periodic response, some others not. As a first step, this specificity is exploited to localize spatially the preferential flow paths network with a structural deterministic inversion approach. Then, inversions of the weekly averaged head variations are performed on the first and the last weeks of the dataset (showing different drawdown trends) to assess the hydraulic diffusivity field. Finally, the tracer restitutions at the spring are inverted in order to optimize the conduits diameters associated to the fastest flow paths.The proposed approach permits a satisfying reproduction of the hydraulic heads and tracer measurements. Simulations with the optimized model on a week of the dataset not used in the inversion allows to validate the inversion result. The inverted model suggests that daily periodic responses represent efficient constraints for the localization of the preferential flows, while the integration of weekly head variations provides a characterization of the hydraulic properties between each well and the closest preferential flow paths, and tracer restitutions permit to characterize the velocities in the preferential paths.

The connection between Submarine Groundwater Discharge and seawater quality: The threat of treated wastewater injected into coastal aquifers

Submarine Groundwater Discharge (SGD) delivers nutrients to the coastal sea triggering phytoplankton blooms, eutrophication, and can also serve as a pathway for contaminants. Wastewater treatment plants (WWTP) including injection wells in coastal areas influence coastal aquifers and might impact the composition and magnitude of SGD fluxes. In tourist areas, wastewater treatment may be less efficient and larger in volume during high seasons, potentially impacting nutrient fluxes from SGD and exacerbating environmental impacts. This study analyzes the nutrient transfer from treated wastewater injection in karstic aquifers to the coastal sea via SGD, considering the impacts of tourism seasonality. This study is conducted in Cala Deià, a small cove in the Balearic Islands, a Mediterranean tourist destination. The findings suggest that the seasonality of tourism, leading to variations in the volume of wastewater treated in the WWTP, influences the dynamics of the coastal aquifer. This leads to increased SGD water and nutrient fluxes to the sea in summer, i.e. the peak tourist season. The measured DIN, DIP, and DSi inventories in the cove are much larger in August than in April (3, 10, and 1.5 times higher, respectively) due to higher input of nutrients in summer due to SGD impacted by the WWTP. These elevated nutrient flows can support algal blooms in the cove, compromising water quality for local swimmers and tourists. Indeed, in August, shoreline stations exhibited eutrophic Chl-a concentrations, with peaks reaching approximately 4 mg Chl-a L−1. These elevated levels suggest the presence of an algal bloom during the survey. The anthropogenic origin of SGD-driven nutrients is traced in seawater and seagrass meadows, as evidenced by high ∂15N signatures indicative of polluted areas. Thus, the high pressure exerted on coastal areas by tourism activities increased the magnitude of SGD nutrient fluxes, thereby threatening coastal ecosystems and the services they provide.

Recharge Assessment in Greek Karst Systems: Methodological Considerations and Implications

The recharge of karst aquifers is significant for the effective management of groundwater resources, and its estimation should be tailored to accommodate their specific hydrogeological characteristics. This study provides a two-step methodological approach for the determination of karst aquifer recharge. Initially, mean annual recharge rates were estimated in the karst system of Ziria (Southern Greece) utilizing the APLIS and modified APLIS methods in order to decipher which was the most suitable version for recharge assessments. The results indicated similar mean recharge rate values at 42.7% and 41.4%, respectively, but significant differences in the spatial distribution. The modified methodology emerged as a more accurate and realistic approach, mainly due to the incorporation of permeability assessments. The final phase of the methodological approach involved the application of modified APLIS in two additional karst hydrosystems, Planitero and Xiromero, while a quantitative cross-comparison of the recharge rates was obtained for a deeper understanding of the factors controlling the groundwater recharge process. In Ziria, recharge rates exhibit a relatively uniform distribution throughout the area, with a median value of 46.7%. Conversely, in Planitero, High recharge rates (60–80%) occupy 56.8% of the surface, while in Xiromero, Moderate recharge rates (40–60%) dominate, representing 53.4% of the land coverage. These variations underscore the spatial heterogeneity of recharge within the karst systems, highlighting the importance of considering local geological and hydrological conditions in its assessments. The methodological approach of this study is flexible and can be adapted to different karst sites for the determination of recharge regimes, contributing to the alleviation of the groundwater depletion issue.

Comparing the Fate and Transport of MS2 Bacteriophage and Sodium Fluorescein in a Karstic Chalk Aquifer.

Groundwater flow and contaminant migration tracing is a vital method of identifying and characterising pollutant source-pathway-receptor linkages in karst aquifers. Bacteriophages are an attractive alternative tracer to non-reactive fluorescent dye tracers, as high titres (>1012 pfu mL-1) can be safely released into the aquifer, offering improved tracer detectability. However, the interpretation of bacteriophage tracer breakthrough curves is complicated as their fate and transport are impacted by aquifer physicochemical conditions. A comparative tracer migration experiment was conducted in a peri-urban catchment in southeast England to characterise the behaviour of MS2 bacteriophage relative to sodium fluorescein dye in a karstic chalk aquifer. Tracers were released into a stream sink and detected at two abstraction boreholes located 3 km and 10 km away. At both sites, the loss of MS2 phage greatly exceeded that of the solute tracer. In contrast, the qualitative shape of the dye and phage breakthrough curves were visually very similar, suggesting that the bacteriophage arriving at each site was governed by comparable transport parameters to the non-reactive dye tracer. The colloid filtration theory was applied to explain the apparent contradiction of comparable tracer breakthrough patterns despite massive phage losses in the subsurface. One-dimensional transport models were also fitted to each breakthrough curve to facilitate a quantitative comparison of the transport parameter values. The model results suggest that the bacteriophage migrates through the conduit system slightly faster than the fluorescent dye, but that the former is significantly less dispersed. These results suggest that whilst the bacteriophage tracer cannot be used to predict receptor concentrations from transport via karstic flow paths, it can provide estimates for groundwater flow and solute contaminant transit times. This study also provides insight into the attenuation and transport of pathogenic viruses in karstic chalk aquifers.

KarstNSim: A graph-based method for 3D geologically-driven simulation of karst networks

In karst aquifers, groundwater flow is highly influenced by the interconnected underground cavities and conduits that form the karst network. Modeling karst flows requires the use of spatially distributed approaches accounting for these networks. Their exploration is, however, often complex, and mapping them using indirect methods such as geophysical ones has proven challenging. To overcome these limitations, stochastically simulating discrete karst networks should account for the uncertainties on conduit position and geometry. Only a few existing methods can reproduce realistic and diverse karst morphologies. In this work, we propose a new algorithm, KarstNSim, for simulating discrete karst networks, that incorporates field data to generate a range of possible karst network geometries. It relies on the computation of the shortest path between the inlets and outlets of the network with the use of an anisotropic cost function defined on an n-nearest neighbor graph conformal to geological and structural heterogeneities. This cost function represents the physico-chemical processes that govern speleogenesis – such as erosion and chemical weathering – providing simplified control over the geometry of the generated networks. Our approach reproduces the vadose-phreatic partition visible in the karst networks, by generating sub-vertical conduits in the unsaturated zone and sub-horizontal ones in the saturated part. It encompasses geological parameters such as inception surfaces, fractures, permeability, and solubility of layers, along with considering the hydrological context of recharge by assigning relative weights to the inlets. We simulate various synthetic models to demonstrate the influence of different input parameters on the spatial organization of the past and present karst flows. We also apply KarstNSim to a real case study – the Ribeaucourt system (Meuse, France) – and use available data to generate realistic karst networks. The lack of knowledge on the Ribeaucourt network geometry is alleviated by the use of multiple scenarios corresponding to different sets of weights in the cost function, resulting in a variety of network morphologies which represent the uncertainty on the real geometry. The geometrical and topological metrics computed on the simulated networks generally overlap with those of real ones, suggesting the reliability of the method.

A New Index to Assess the Effect of Climate Change on Karst Spring Flow Rate

Karstic aquifers, because of their conduit system, are susceptible to climate change. Ten karst springs in the Zagros region were selected to investigate the impact of climate change under three CMIP6 scenarios: SSP1-1.9, SSP2-4.5, and SSP5-8.5. This study was conducted in three steps: downscaling climate projection, analyzing spring discharge time series, and introducing a new index to assess the impact of climate change on spring flow rate. Applying LARS-WG6, precipitation was downscaled at 14 stations in the study area. Moreover, time series and trend analysis showed that the selected springs have experienced a decrease in their flow rate. Assuming the covariance function between precipitation and spring discharge is constant, new indices (i.e., IQd, IdQd, and Icc) were introduced to highlight the effect of climate change according to the three scenarios. dQd is the variability of spring discharge from past to future, IdQd is spring discharge variability over the historical data, and Icc is the effect of precipitation and spring discharge change together. Icc has a range from −0.25 to 0.25 below and above, which is indicative that two extreme conditions including the spring dryness and overflow are in effect, respectively. The main results revealed that the degree of impact at each spring is a function of climate change scenarios and hydrogeological characteristics of the karstic systems. A more noticeable negative trend in spring flow rate is observed for the karst springs characterized by a dominant conduit flow regime and low matrix storage, located in the areas with low cumulative rainfall, and has a stronger relationship with precipitation. Based on the results, decisions on the management of karst water resources should be made considering where the springs bear free surface and pressurized flow conditions.

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.

Using multivariate statistics to link major ion chemistry changes at karst springs to agriculture

Agricultural contamination is a pervasive problem in karst landscapes that negatively impacts groundwater quality, and subsequently harms hydraulically connected springs and associated ecosystems. While most research on agricultural contamination in karst aquifers has focused on nitrate, agricultural operations also alter major ion chemistry in groundwater, which can have significant implications on water rock interactions, geochemical interpretability, and enrichment of harmful contamination. Additionally, recent evidence has questioned the singular role of nutrient enrichment on aquatic ecosystem degradation, thus prompting further investigation into the broader extent of agricultural contamination for improved water quality restoration. In karst aquifers, quantifying the disruption of groundwater quality from agricultural contamination at springs can be challenging as major ions associated with agricultural contamination may be conflated with natural processes and interpretation further compounded by complex flow dynamics. These factors can render traditional statistical methods inadequate for identifying more comprehensive groundwater quality disruptions. To overcome these challenges, we use principal component analysis to classify the major ion covariance signature of agriculturally impaired and unimpaired groundwater and develop a new impairment metric to characterize the relative major ion impairment from agriculture to karst spring water quality. We apply our method to a primarily rural region in northern Florida whereby major ion chemistry has been changing over the last several decades. Our results identify a significant difference in major ion covariance structure between agriculturally impaired and unimpaired groundwater which is exploited to directly link major ion changes at karst springs to agriculture via relating a derived impairment metric to measured nitrate. Our newly developed method is easy to apply and can be used to identify more comprehensive contamination from pollutant sources and prioritize relative water quality restoration.

Spatiotemporal variation in cave percolation waters: A functional approach

Understanding the mechanisms controlling spatial heterogeneity of drip water percolation into caves is essential for interpreting karst aquifer recharge and speleothem isotopic and geochemical records for paleoclimate analyses. Here we present the first analysis of drip rate variability using a novel time-varying Functional Principal Component Analysis (FPCA), validated against drip water stable isotope composition. Twenty-six drip sites were monitored across Harrie Wood Cave, south-east Australia, over a 2.5 year period. A positive relationship of cave drip water hydrology with precipitation and soil moisture was identified, with soil moisture recording the strongest relationship. FPCA was used to classify drip-water flow (percolation) pathways based on temporal shifts in the drip rate time series. Our results reveal that three percolation classes can be used to explain water movement within the cave: storage baseflow, fracture baseflow and overflow. The successful application of FPCA in this study suggests that this statistical technique will be useful for the analysis and interpretation of other large, discontinuous hydrological datasets.

Identification of high-permeability and water-rich zones in a fractured karst water source area based on the hydraulic tomography method

Determining the spatial hydraulic parameters of aquifers is of paramount importance for the development and utilization of groundwater in water source areas. However, most methods are unable to accurately characterize the heterogeneity of aquifers because of the influences of complex geology, hydrogeology, structure and insufficient information. Hydraulic tomography (HT) is an effective and robust method for hydraulic parameter inversion to characterize heterogeneity. In this study, a two-dimensional HT model was established to simulate hydraulic parameter fields for a fractured karst aquifer in the Zhangji water source area (Xuzhou city, China). Steady-state hydraulic tomography (SSHT) and transient hydraulic tomography (THT) with different prior information conditions were applied to determine the high-permeability area and water-rich zones, and to infer possible locations of karst conduits and strong runoff paths in the water source area. The results indicate that: (1) THT can better reveal larger-scale heterogeneous characteristics of hydraulic properties compared to SSHT. The hydraulic conductivity distribution estimated by THT correlates well with the geological conditions of the area, especially in identifying water-rich zones formed by faults and other formations. (2) When combined with geological conditions, HT possesses the ability to identify strong runoff paths in fractured karst aquifers. (3) Accurate prior information, which is more relevant to the head response information, is conducive to obtaining more acceptable HT results. (4) THT can simultaneously identify high-permeability and high-storage zones, and the analysis of the overlap of these two regions has the potential to predict underground water reservoirs, which is of the utmost importance for the development and utilization of groundwater resources in water source areas.

Characteristics and Influencing Factors of Dissolved Organic Matter in Cave Drip Water—A Case Study from Furong Cave, Southwest China

Understanding the hydrological processes of dissolved organic matter (DOM) in the surface karst zone is crucial for the utilization and safety of groundwater resources. However, research on DOM in drip water from karst caves is limited. In this study, continuous monitoring was conducted for four years at four drip water monitoring sites (MP1, MP2, MP3, and MP9) in Furong Cave, Southwest China. The three-dimensional fluorescence excitation–emission spectroscopy (3D-EEM) and fluorescence region integration (FRI) methods were employed, along with correlation analysis, to investigate the spectral characteristics, composition, sources, and influencing factors of the DOM in the drip water of Furong Cave. The results indicated that (1) the three-dimensional fluorescence peaks in the drip water were classified into six categories: A, B, C, T, Ti, and M. (2) The dominant source of the drip water DOM is endogenous organic matter. (3) The FRI analysis indicates a relatively high proportion of Type II substances in the drip water, predominantly composed of tryptophan-like substances. (4) The DOM in the drip water of Furong Cave was influenced by various factors, including the mixing effect of “new” and “old” water, water residence time in karst systems, and water–rock interactions (WRI), resulting in the complex responses of drip rates and DOM to surface precipitation and temperature at different drip water sites. This study provides a reference for comparative research on DOM in cave drip water in karst regions, which contributes to a better understanding of the migration mechanism of DOM in karst aquifers under different climate and karst ecological conditions.

Yeast Biodiversity of Karst Waters: Interest of Four Culture Media and an Improved MALDI-TOF MS Database

Karst aquifers are a significant source of drinking water and highly vulnerable to pollution and microbial contamination. Microbiological regulations for the quality of drinking water mostly focus on bacterial levels and lack guidance concerning fungal contamination. Moreover, there is no standardised microbial analysis methodology for identifying fungi in water. Our main objective was to establish the most effective culture and identification methodology to examine yeast diversity in karst waters. We assessed the comparative efficacy of four culture media (CHROMagar Candida, dichloran glycerol 18% [DG18], dichloran rose Bengal chloramphenicol [DRBC], and SYMPHONY agar) for yeast isolation from karst water samples. Furthermore, we investigated the comprehensiveness of databases used in MALDI-TOF mass spectrometry (MALDI-TOF MS) for identifying environmental yeast species. In total, we analysed 162 water samples, allowing the identification of 2479 yeast isolates. We demonstrate that a combination of four culture media, each with distinct specifications, more efficiently covers a wide range of yeast species in karst water than a combination of only two or three. Supplementation of a MALDI-TOF MS database is also critical for analysing environmental microbial samples and improved the identification of yeast biodiversity. This study is an initial step towards standardising the analysis of fungal biodiversity in karst waters, enabling a better understanding of the significance of this environmental reservoir in relation to public health.

Comparative assessment of groundwater quality challenges confronting Mediterranean coastal aquifers: A synoptic review

Coastal aquifers, essential for agriculture, domestic, and industrial water supply, face significant sustainability challenges, particularly in Mediterranean regions, due to the adverse effects of human activities, climatic threats, and natural geological conditions causing salinization. This paper synthesizes findings from studies on three Mediterranean coastal aquifers— the Pleistocene aquifer in the Gaza Strip, the karstic aquifer in Bokkoya, and the MioPlio-Quaternary aquifer in Djeffara of Medenine—highlighting the escalating deterioration due to over-exploitation and arid conditions. It compares the methodologies used in these studies, such as statistical analysis, geostatistical modeling, and groundwater quality indices, to assess groundwater quality and address salinization challenges. This review aims to identify methodological limitations and propose strategies to fill research gaps, particularly in combating seawater intrusion, to enhance the sustainable management of these critical water resources.

Karstic Aquifers and Climate Refugia: A Preliminary Outline History of Water-Management Strategy in Bronze and Iron Age Southeast Arabia

Karstic Aquifers and Climate Refugia: A Preliminary Outline History of Water-Management Strategy in Bronze and Iron Age Southeast Arabia

Study on Dynamic Evolution Law of Fracture and Seepage in Overlying Strata during Mining under Karst Aquifer

Study on Dynamic Evolution Law of Fracture and Seepage in Overlying Strata during Mining under Karst Aquifer

Tracking flowpaths in a complex karst system through tracer test and hydrogeochemical monitoring: Implications for groundwater protection (Gran Sasso, Italy)

Groundwater in karst aquifers is frequently tapped for drinking purposes, due to frequent huge volumes of resources. Unfortunately, vulnerability of these aquifers can be high, due to possible fast transfer of recharge water on springs by the karst network. On Gran Sasso Mountain regional aquifer, several springs are subjected to drinking withdrawal and an updated evaluation of their potential is now a fundamental issue to be considered, facing climate change effects, which reflect on variation of discharge regimen and values. To distinguish between different contribution of spring recharge, a tracer test has been carried out on the Vitella d’Oro spring, fed both by the regional aquifer and by a local system exposed to karst features developed in the Rigopiano Conglomerates formation. Thanks to hydrogeological, hydrogeochemical and isotopic data, a conceptual model of spring recharge has been proposed and subsequently validated by the tracer test results. All information confirms the superimposition on the regional base flow, by a relevant contribution of the karst network, influencing the spring discharge in recharge periods. In detail, a fast flow component is responsible for discharge peaks and frequently of turbidity events, having a mean velocity ranging from 30 to 70 m/h in the aquifer. Besides of this fast flow, an additional aliquot of the recharge is due to the same local aquifer, but slower flow clearly identifiable by hydrochemistry and isotopic data. Thanks to these findings, a renewed management of the spring has been suggested, considering the different degrees of aquifer vulnerability (turbidity occurrence) directly related to the discharge regimen.

Exploratory study of groundwater infiltration in lough arrow using thermal imagery and geochemical tracers

abstract: Groundwater dynamics were investigated in Lough Arrow in the north-west of Ireland. The lake is designated under the EU Habitats Directive (92/43/EEC) for its representative Annex I habitat from which charophytes are the main feature of interest. Groundwater interactions are considered a significant component of the lough's hydrology with its Annex I archetype deemed vulnerable to the pressures of nutrient pollution. This study utilised thermal imagery and geochemical tracers to detect and quantify groundwater inputs. Landsat 7 Thermal Infrared Imagery from two dates in both summer and winter was processed to investigate the potential occurrence of localised seepage points of groundwater captured as anomalous spatial plumes. In situ surveys of radon were undertaken in summer to confirm the presence of groundwater in the lake. Results were subsequently modelled to estimate the total groundwater discharge into the lake during the survey period. Additionally, a conceptual site model was created to characterise potential pressures to the lake from groundwater interactions within the catchment. Thermal imaging failed to conclusively identify thermal plumes consistent with localised groundwater inputs in Lough Arrow. However, the radon inventory applied to the mass balance estimated the lake was subject to 6,722m3 of groundwater influx per day at the time of the survey. The conceptual site model determined the lake was subject to potential pollution pressures from groundwater due to the catchment's physical attributes and land use activities. This is primarily due to the presence of a karst aquifer in the majority of the catchment coupled with significant agricultural and forestry practices. The vulnerability of the aquifer and the influx of groundwater to the lake highlights the potential threat groundwater contamination may pose to the favourable conservation condition of Lough Arrow's Annex I habitat.

An approach to the global well-posedness of a coupled 3-dimensional Navier-Stokes-Darcy model with Beavers-Joseph-Saffman-Jones interface boundary condition

<abstract><p>This study focused on investigating the global well-posedness of a coupled Navier-Stokes-Darcy model with the Beavers-Joseph-Saffman-Jones interface boundary condition in the three-dimensional Euclidean space. By utilizing this approach, we successfully obtained the global strong solution of the system in the three-dimensional space. Furthermore, we demonstrated the exponential stability of this strong solution. The significance of such coupled systems lies in their pivotal role in the analysis of subsurface flow problems, particularly in the context of karst aquifers.</p></abstract>