Aquifer Structure Research Articles

Geophysical and geological contribution to determining the neritic limestone aquifer structure of Hammam Bradaa – El Fedjoudj (Seybouse medium), Northeastern Algeria

Purpose. The neritic aquifer that extends between Hammam Bradaa and El Fedjoudj, despite its faulted and in places karstified structure, due toits capacity and lateral extension represent a strategic resource for the region. Its waters are used for drinking water supply in the neighboring towns: Heliopolis, Guelaat Bousbaa, Nechmaya and El Fedjoudj, as well as in part of the Annaba wilaya. These neritic limestones outcrop at Djebel Debagh and Bouzitoune. Methods. In order to better study the potential of this aquifer, it is essential to identify its geometry and structure by analyzing the geological data, mechanical drilling data and geophysical data analysis through an electrical survey campaign. Findings. The main results indicate that the study area has identified two important aquifer formations: a formation in the Mio-Plio-Quaternary alluvium consisting of clays, marls, gravel and sand; the second is a deep carbonate aquifer composed of fissured and karstified neritic limestone of the Cretaceous age of variable depth ranging within 50 and 350 m. Originality. The originality of the study is in the fact that the studied area is characterized by the presence of thick, fractured and karstified carbonate formations, which are intensively tectonized and have significant aquifer potential. Practical implications. The correlation of geological data with geophysical data made it possible to conclude that the stu-died area is a sedimentary basin bounded by faults predominantly oriented to the south-west and north-east, forming a highly fractured unit consisting of horsts and grabens. The significant water potential reservoir is formed essentially by carbonate geological formations, highly fractured with the presence of karst forms represented by resistant horizons.

Development of a Distributed Control System for the Hydrodynamic Processes of Aquifers, Taking into Account Stochastic Disturbing Factors

Uncontrolled, intensive extraction of water from mineral water deposits can lead to negative consequences, the penetration of sewage beyond the water-resistant layers, the destruction of the structure of aquifers and the further loss of the mineral water source. This paper deals with the issues of control of hydrodynamic processes under random impacts in the aquifers of mineral water deposits using mathematical models. The analysis of the technological process of mineral water extraction and statistical analysis of retrospective monitoring data of the deposit were carried out in order to confirm the hypothesis of the stationarity of the processes under consideration. A mathematical model of hydrogeological processes occurring under random disturbances has been constructed. The influence of such disturbances is expressed in deviations of the head in the aquifer from the desired value. Controlling consists of maintaining the head of aquifers at a given level in order to preserve the hydro-mineral base of the region. Lumped and distributed controllers are proposed as corrective blocks, the parameters of which are calculated by frequency methods of controller synthesis. The use of a closed loop control system minimizes the influence of random effects.

Constrained 3D gravity modelling of a multilayered aquifer system in an arid region: Case of Sebkhat El Bhira basin (Central Tunisia)

In arid regions, the significant and increasing demand for water has caused intensive exploitation of groundwater resources. The Sebkhat El Bhira basin, shared between Kairouan and Sidi Bouzid governorates in Central Tunisia, presents a multi-layer aquifer system composed of four main reservoirs: the sandy Mio-Plio-Quaternary (MPQ), the Paleogene, the Upper Cretaceous limestone and the Lower Cretaceous limestone, sand and dolomite reservoirs. The intensive exploitation of the shallower aquifer has led to a drawdown of the water table and a deterioration of water quality. The purpose of this work is to characterize the geometry of these potential reservoirs for sustainable water management, by applying geophysical methods.We used land gravity data covering 1387 Km2, twelve seismic reflection profiles with a total length of approximately 308 Km, two petroleum wells and 39 water wells.The processing of the gravity data shows the subdivision of Sebkhat El Bhira basin into four sub-basins and the identification of principal and minor lineaments oriented N-S, NE-SW, NNW-SSE, and E-W. The interpretation of reflection seismic lines allowed the identification of five main horizons: the bottom of MPQ, the bottom of Paleogene, the bottom of Upper Cretaceous, the bottom of Lower Cretaceous and the bottom of Jurassic (Top Triassic). For a better understanding of the subsurface structuring, we built twelve two-dimensional (2D) forward models of density. These latter are constrained in the area inside the plain by the converted time-to-depth seismic reflection data and wells. The interpolation of the 2D models allowed us to visualize the lateral and vertical thickness variations of each of the chrono-stratigraphic units and to produce a three-dimensional (3D) density model. The geophysical modelling allowed us to highlight the post-Jurassic tectonic and the Triassic halokinesis activities responsible for the actual subsurface geometry.The Triassic evaporites, in the study area, flow and rise through the fault's subvertical branches. It outcrops locally at the North-South Axis mountains range (NOSA). Inside the plain, the evaporitic uprising shows different schemes like diapir, wall and pillow. This Triassic halokinesis, with the tectonics, has shaped the deep morphology of the basin and managed the available space for the deposition of the overlying stratigraphic formations.The general structure of the potential aquifers is characterized by a depth which increases from the East to the West and towards the central part of the plain. Inside the basin, the structure of these aquifers is complex.

Extension of CO2 storage life in the Sleipner CCS project by reservoir pressure management

In this study, we performed reservoir simulations to investigate CO2 storage in the saline aquifer in Sleipner by reservoir pressure management. Results show that by producing water at the lowest aquifer structure far away from the CO2-water contact in the vertical direction, an additional 73% (24 Mt) CO2 can be stored compared to the case without water production. This extra CO2 stored can generate a revenue of $800 millions at the prevailing carbon tax of $69 per ton in Norway. There is 5.31 Mtpa of natural CO2 production from 124 gas and oil fields in Norway in 2020. Of these, the Marulk, Dvalin, Skarv, Morvin, Åsgard, Kristin, Tyrihans and Mikkel fields in the Norwegian Sea, and the Kvitebjørn, Valemon, Visund, Gudrun fields in the North Sea produce a total of 4.55 Mtpa CO2 potentially supplying CO2 to Sleipner for project life extension.

Delineation of the aquifer structure and estimation of hydraulic properties on the flanks of Mount Meru, Northern Tanzania

Understanding of the aquifer structure and its hydraulic properties provides comprehensive knowledge for proper groundwater utilisation and management. This study delineated the aquifer structure using litho-hydrostratigraphical cross-sections, and estimated the hydraulic parameters using single well pumping tests for various locations on and around Mount Meru. Results show that, the aquifer system on the flanks of Mount Meru is a sloping aquifer with sloping beds. On the far east of the eastern flank, the aquifer is composed of debris avalanche deposits, while on the north-eastern and west flanks the aquifer is composed of weathered fractured lava, whereas on the south-western flank, the aquifer is composed of different layers: pyroclastics on the top, weathered fractured lava, weathered pyroclastics, and weathered fractured lava at the bottom. The aquifer is semi-confined on the north-eastern and western flanks; on the north-eastern flank, the overlying debris avalanche deposits acting as an aquitard, while on the western flank, the overlying layers: pyroclastics and unweathered lava acting as an aquitard and aquiclude, respectively. The aquifer is unconfined on the far east of the eastern flank and south-western flank. The transmissivity of the aquifer on the north-eastern flank is substantially increasing with increasing depth, while variable, both at the shallow depth of hand-dug wells and the larger depth of boreholes, on the south-western flank; indicating aquifer heterogeneity. On the north-eastern flank, the topmost part of the aquifer, exploited by hand-dug wells, has a low transmissivity (T = 1.3 m2/d) and potential for smaller withdrawals for local water supply with limited consumption, while the upper part of the aquifer, captured by boreholes, has an intermediate transmissivity (T = 35 m2/d) and potential for local water supply, whereas the deeper part of the aquifer has a high transmissivity (T = 788 m2/d) with potential of somewhat regional importance. On the western flank, the aquifer has a very low transmissivity (T = 0.4 m2/d) and potential for local water supply with limited consumption. On the south-western flank, on average, the topmost part of the aquifer, exploited by hand-dug wells, has very low to intermediate transmissivity (range of T: 0.3–21 m2/d), leading to variable potential for smaller withdrawals for local water supply (private consumption), whereas the deeper part of the aquifer, captured by boreholes, has low to intermediate transmissivity (range of T: 9–43 m2/d) and potential for local water supply.

Identification of hydrogeological features using remote sensing and electromagnetic methods in the hard- rock formations of the Cameroon coastal plain (Central Africa): implications for water borehole location

Sustainable access to groundwater in rural areas such as the coastal plain of Cameroon (Sub-Saharan Africa) is an ongoing challenge due to the lack of knowledge on hydrogeological features favorable to the identification of potential groundwater targets. In this study, remote sensing imagery and electromagnetic surveys were used to assess the hydrogeological potential of the Ngwei locality hard-rock. Automatic extraction of lineaments from Landsat-8 images indicated 200 lineaments with a total length of 293 km corresponding to 0.6 km of lineament per square km, suggesting a median groundwater potential. The main directions are NNW-SSE, N-S and NNE-SSW while the secondary directions are NW-SE; NE-SW and E-W. Surface geophysical investigations showed an electromagnetic field of 0 to 0.57 mV. These surface and borehole investigations confirmed basement fracturing and spatial variations in weathering profiles whose aquifer structure and yields (4.4 m3/h on average) are similar to what is commonly observed in the hard rock aquifer. The correlation between borehole productivity and the nearest lineament suggests that these lineaments are not (at least not all) surface traces of regional discontinuities that may act as major groundwater flow paths. This implies that the productivity of the Ngwei aquifers appears to be dependent on the fissured/fractured horizon, so exploration should be based on a contextual rather than general approach.Analysis of the consistency of the results with the field data allowed to propose a preliminary electromagnetic log model and also shows the value of using this approach to assess the availability of potential groundwater zones in hard rock areas. However, the methods used should be supported by further investigations to allow a better understanding of the hydrogeological characteristics of this part of the Central African basement so that no one is left behind in the achievement of the SDGs.

Characterization of Shallow Groundwater Circulation Based on Chemical Kinetics: A Case Study of Xiong’an New Area, China

Xiong’an New Area, located in the middle of the North China Plain, will have been built as a “city of the future.” Urban planning and construction need to comprehensively consider the constraints of hydrogeological conditions such as aquifer structure and parameters. As the main aquifer in this area, the paleo-channel is heterogeneous and anisotropic, and the two-dimensional hydraulic conductivity in each horizontal direction cannot be obtained from aquifer tests. Therefore, this study adopts a chemical kinetics method to calculate the ionic activity and mineral saturation indices of shallow groundwater, determine the groundwater chemical potential field, and construct a horizontal two-dimensional groundwater chemical kinetics model. This model is used to calculate the hydraulic conductivity, flow rate and retention time of groundwater in areas of different chemical kinetics, as well as evaluate horizontal heterogeneity of the Quaternary paleo-channel aquifer. The results indicate that the groundwater chemical potential field can reflect the characteristics of the groundwater seepage field in each horizontal direction. The paleo-channel is the main channel of groundwater circulation, which shows the statistical difference of its permeability. Alluvial and lacustrine strata affect groundwater circulation due to their different hydrogeological structures and permeability. The groundwater chemical kinetics results of hydraulic conductivity along the paleo-channel are approximately consistent with traditional hydrogeological calculation results derived from aquifer test data. Hydraulic conductivity is higher in the extension direction of the paleo-channel, and lower if the path crosscuts multiple paleo-channels. This feature can be used to build a hydrogeological structure model combined with drilling data. Furthermore, excessive groundwater exploitation will change the actual flow rate and retention time of groundwater, thereby affecting the groundwater circulation conditions. This study of groundwater circulation in Xiong’an New Area by means of chemical kinetics makes up for the deficiency in the study of the unconsolidated sedimentary aquifer anisotropy within the paleo-channel.

Application of hydrochemical and isotopic data to determine the origin and circulation conditions of karst groundwater in an alpine and gorge region in the Qinghai–Xizang Plateau: a case study of Genie Mountain

Complex structural systems, lithological differences, and extreme heterogeneity of aquifers combine to create a complex karst aquifer structure in alpine and gorge areas; however, because of the topography, direct investigation of aquifer structure is difficult. In this study, field survey, hydrochemical, and isotopic data are analyzed to reveal the development of karst groundwater and to describe the karst water cycle in Genie Mountain, Qinghai–Xizang Plateau. The results show that groundwater circulation is mainly controlled by active fracture. Atmospheric precipitation and melting ice and snow are the groundwater recharge sources, and recharge areas are mostly located in high mountains above 4500 m a.s.l. The direction of groundwater flow is mostly controlled by the Jinshajiang active Fault, with drainage areas at the intersection of multiple faults. There are two regional karst water runoff conduits. One is along the Edexi-Hongjunshan Fault, where groundwater runs from south to north; the other is along Gangtonglong Fault, where groundwater runs from north to south, and is discharged at Gangtonglong gully. The groundwater cycle can be divided into three levels: epikarst water circulation; mid to deep karst water circulation; and deep geothermal water circulation. The karst springs located in the outlet of the Huolong gully contain markedly higher levels of Na+ and SO42− than other karst springs because of the leaching effect of groundwater on mirabilite. The presence of evaporites in Huolong gully indicates that the plateau planation formed by the structural uplift will change the local climatic conditions, and then affect the groundwater circulation process and the water–rock reaction in the process.

Groundwater flow evaluation using a groundwater budget model and updated aquifer structures at an alluvial fan of Echi-gawa, Japan

Groundwater flow evaluation using a groundwater budget model and updated aquifer structures at an alluvial fan of Echi-gawa, Japan

An integrated inversion framework for heterogeneous aquifer structure identification with single-sample generative adversarial network

Generating reasonable heterogeneous aquifer structures is essential for understanding the physicochemical processes controlling groundwater flow and solute transport better. The inversion process of aquifer structure identification is usually time-consuming. This study develops an integrated inversion framework, which combines the geological single-sample generative adversarial network (GeoSinGAN), the deep octave convolution dense residual network (DOCRN), and the iterative local updating ensemble smoother (ILUES), named GeoSinGAN-DOCRN-ILUES, for more efficiently generating heterogeneous aquifer structures. The performance of the integrated framework is illustrated by two synthetic contaminant experiments. We show that GeoSinGAN can generate heterogeneous aquifer structures with geostatistical characteristics similar to those of the training sample, while its training time is at least 10 times faster than that of typical approaches (e.g., multi-sample-based GAN). The octave convolution layer and multi-residual connection enable the DOCRN to map the heterogeneity structures to the state variable fields (e.g., hydraulic head, concentration distributions) while reducing the computational cost. The results show that the integrated inversion framework of GeoSinGAN and DOCRN can effectively and reasonably generate the heterogeneous aquifer structures.

An Improved Groundwater Model Framework for Aquifer Structures of the Quaternary-Formed Sediment Body in the Southernmost Parts of the Mekong Delta, Vietnam

The Ca Mau peninsula (CMP) is a key economic region in southern Vietnam. In recent decades, the high demand for water has increased the exploitation of groundwater, thus lowering the groundwater level and leading to risks of degradation, depletion, and land subsidence, as well as salinity intrusion in the groundwater of the whole Mekong Delta region. By using a finite element groundwater model with boundary expansion to the sea, we updated the latest data on hydrogeological profiles, groundwater levels, and exploitation. The basic model setup covers seven aquifers and seven aquitards. It is determined that the inflow along the coastline to the mainland is 39% of the total inflow. The exploitation of the study area in 2019 was 567,364 m3/day. The most exploited aquifers are the upper-middle Pleistocene (qp2–3) and the middle Pliocene (n22), accounting for 63.7% and 24.6%, respectively; the least exploited aquifers are the upper Pleistocene and the upper Miocene, accounting for 0.35% and 0.02%, respectively. In the deeper aquifers, qp2–3 and n22, the change in storage is negative due to the high exploitation rate, leading to a decline in the reserves of these aquifers. These groundwater model results are the calculations of groundwater reserves from the coast to the mainland in the entire system of aquifers in the CMP. This makes groundwater decision managers, stakeholders, and others more efficient in sustainable water resources planning in the CMP and Mekong Delta (MKD).

Migration mechanism of pollutants in karst groundwater system of tailings impoundment and management control effect analysis: Gold mine tailing impoundment case

Tailing impoundments are sources of high-risk pollution due to the accumulation of mining tailings and wastewater. As a result of thin karst overburden as well as actions of precipitation leaching and scouring, pollutants in the tailing impoundment can easily enter the karst aquifer through sinkholes, corrosion fissures and other channels, causing pollution to groundwater and is not conducive to the utilization of groundwater resources. Through the analysis of hydrogeological conditions, the boundary condition, aquifer structure, source-sink term, and hydraulic characteristics of groundwater are generalized, the groundwater flow model and solute migration model were established, and the migration mechanism of the characteristic pollutant (arsenic ions) was analyzed. The research results show that with the passage of time, the area of high-concentration polluted areas showed a decreasing trend, and the area of low-concentration polluted areas showed a decreasing trend; Under the influence of groundwater flow and the pollutant concentration difference, the pollution plume presents an elliptical shape with the long axis along the direction of water flow; The migration law of pollutants is controlled by the direction of groundwater flow, the concentration difference of pollutants and the permeability of the aqueous medium; The direction of groundwater flow determines the migration direction of pollutants, the concentration difference of pollutants determines the diffusion range of pollutants, and the permeability of aqueous media determines the migration speed of pollutants; Through the prediction of the pollutant control effect, it can be found that a reasonable number of pumping wells that do not interfere with each other are arranged around the groundwater flow channel, which can effectively control the migration of pollutants to the target area. The above research reveals the migration mechanism of pollutants in karst groundwater and main control factors, and also provides a scientific basis for the protection of groundwater resources in karst areas.

Estimating Effective Porosity in Bedrock Aquifers.

Flow in many bedrock aquifers is through fracture networks. Point to point tracer tests using applied tracers provide a direct measure of time of travel and are most useful for determining effective porosity. Calculated values from these tests are typically between 10-4 and 10-2 (0.01% to 1%), with these low values indicating preferential flow through fracture and channel networks. Tracer tests are not commonly used in site investigations, and specific yield is often used as a proxy for effective porosity. The most popular methods have used centrifuge measurements, water table fluctuations, pumping tests, and packer tests. Specific yield varies substantially with the testing method. No method is as reliable as tracer testing for providing estimates of effective porosity, but all methods provide complementary insights on aquifer structure. Temporal and spatial scaling effects suggest that bedrock aquifers have hierarchical structures, with a network of more permeable fractures and channels, which are connected to less permeable fractures and to the matrix. Consequences of the low effective porosities include groundwater velocities that often exceed 100 m/d and more frequent microbial contamination than in aquifers in unconsolidated sediments. The large uncertainty over the magnitude of effective porosity in bedrock aquifers makes it an important parameter to determine in studies where time of travel is of interest.

Improving the Methods for Processing Hard Rock Aquifers Boreholes’ Databases. Application to the Hydrodynamic Characterization of Metamorphic Aquifers from Western Côte d’Ivoire

Statistical analysis of a borehole database, linear discharges, and water strikes processing enabled an understanding of the structure, geometry and hydrodynamic properties of the metamorphic hard rock aquifers from the Montagnes District, Western Côte d’Ivoire. The database comprises 1654 boreholes among which 445 only were usable for this research work after its pre-processing. Analysis shows that the structure of the aquifer is similar to that observed in several other areas in the world: it developed due to weathering processes, comprises the capacitive saprolite, 10–20 m thick on average, and an underlying transmissive fractured layer, overlying the unweathered impermeable hard rock. The fractured layer is 80 m thick, the first 40 to 45 metres being its most productive zone, with a 11.3 m3/h median productivity. This research shows that metamorphic aquifers exhibit similar aquifer properties (thickness, hydrodynamic parameters) as plutonic ones and that there is interest in using such databases for research and other purposes. However, a rigorous pre-treatment of the data is mandatory, and geological data from published maps must be used instead of the geological data from the database. A previous methodology aiming at processing the boreholes’ linear discharges was improved. It notably appears that the slope method must be preferred to the percentile method.

Identification of Preferential Recharge Zones in Karst Systems Based on the Correlation between the Spring Level and Precipitation: A Case Study from Jinan Spring Basin

The Jinan spring basin is located in the karst area of northern China, where springs serve as important sources of water supply. Several studies on spring protection and water supply have been carried out, and scholars have developed some laws on local groundwater flow dynamic and characteristics of aquifer structures. Unfortunately, there is a lack of detailed research on preferential recharge zones, which are the main recharge pathways of springs. Therefore, this research focuses on identifying preferential recharge zones based on the correlation between the spring level and precipitation. The results show that when precipitation is more intense or lasts longer, there is a stronger correlation between spring level and precipitation. It has been established that the precipitation at Donghongmiao station has the closest relationship with the dynamic of Baotu spring, which is found to be the most significant contribution to spring preservation. Two potential preferential recharge zones in the Jinan spring basin are detected through correlation analysis and geological exploration data. These findings support spring protection and water supply projects in karst regions.

Source, Partition and Ecological Risk of Polycyclic Aromatic Hydrocarbons in Karst Underground River Environment, Southern China

As the main source of drinking water in the karst areas of southern China, underground rivers play an important supporting role in local economic and social development. However, due to the special aquifer structure in the karst areas of southern China, polycyclic aromatic hydrocarbons (PAHs) can easily enter the water environment of underground rivers and cause serious pollution, which will affect the water quality safety. In this study, Qing-shui Spring underground river in Nanning city was selected as a representative of the typical underground river in southern China, and the pollution characteristics, source analysis, ecological risk assessment and diffusion law of PAHs were studied. The results showed that the total concentration of PAHs in groundwater and surface sediments was in the range of 282.42–464.88 ng/L and 400.14–1194.69 ng/g, respectively, and the overall concentrations of PAHs were below the moderate pollution level. Pollutant discharge and physical and chemical properties of PAHs caused the concentrations of PAHs to decrease gradually from upstream to downstream, but the proportion of high-cyclic PAHs increased gradually. According to the source analysis results of the ratio method, PAHs sources in groundwater were mainly biomass combustion sources in the upstream, oil sources in the middle and mixed sources in the downstream. Sediment and groundwater had some differences. The middle and upper reaches were the source of biomass combustion, while the middle and lower reaches were the mixed source. The ecological risk assessment results showed that the ecological risk of groundwater was at a moderate level, while the ecological risk of sediment was at a low level. Benzo[a]anthracene (BaA) and benzo[b]fluoranthene (BbF) in groundwater and acenaphthylene (Acy) and fluoranthene (Flu) in sediment should be monitored intensively. Appropriate control measures should be put in place to prevent further pollution. Partition coefficients varied from 0.73 to 3.85 L/g in the sediment–water interface, increasing with the rise of PAH compounds. All of the organic carbon partition coefficients (Koc) in the sediment–water interface were higher than predicate values based on the typical model of equilibrium distribution; this indicated that PAHs were strongly adsorbed in the sediment. The linear free energy relationship coefficient between the Koc value and the n-octanol-water partition coefficient (Kow) was 0.85 but the slope was lower than 1, indicating that sediment in the Qing-shui Spring underground river had weak lipophilic characteristics and adsorption ability for PAHs. The results can provide a reference for the risk assessment and pollution control of PAHs in the karst underground rivers in South China.

Comparative Assessment of Methods for Coupling Regional and Local Groundwater Flow Models: A Case Study in the Beijing Plain, China

A coupled regional and local model is required when groundwater flow and solute transport are to be simulated in local areas of interest with a finer grid while regional aquifer boundary and major stresses should be retained with a coarser grid. The coupled model should also maintain interactions between the regional and local flow systems. In the Beijing Plain (China), assessment of managed aquifer recharge (MAR), groundwater pollution caused by rivers, capture zone of well fields, and land subsidence at the cone of depression requires a coupled regional and local model. This study evaluates three methods for coupling regional and local flow models for simulating MAR in the Chaobai River catchment in the Beijing Plain. These methods are the conventional grid refinement (CGR) method, the local grid refinement (LGR) method and the unstructured grid (USG) method. The assessment included the comparison of the complexity of the coupled model construction, the goodness of fit of the computed and observed groundwater heads, the consistency of regional and local groundwater budgets, and the capture zone of a well filed influenced by the MAR site. The results indicated that the CGR method based on MODFLOW-2005 is the easiest to implement the coupled model, capable of reproducing regional and local groundwater heads and budget, and already coupled with density and viscosity dependent model codes for transport simulation. However, the CGR method inherits shortcomings of finite difference grids to create multiple local models with inefficient computing efforts. The USG method based on MODFLOW-USG has the advantage of creating multi-scale models and is flexible to simulate rivers, wells, irregular boundaries, heterogeneities and the MAR site. However, it is more difficult to construct the coupled models with the unstructured grids, therefore, a good graphic user interface is necessary for efficient model construction. The LGR method based on MODFLOW-LGR can be used to create multiple local models in uniform aquifer systems. So far, little effort has been devoted to upgrade the LGR method for complex aquifer structures and develop the coupled transport models.

Simulation of solute transport behaviors in saturated karst aquifer system

The karst development makes aquifer have strong anisotropy and heterogeneity. In order to reveal the characteristics of solute transport in the karst fissure–conduit aquifer system, this study presents a physical model of fissure–conduit in laboratory experiments to carry out the solute transport simulation. In this paper, the tracer tests of fissure–conduit combination, fissure, and conduit solute transport process in saturated flow are designed. We found that different aquifer structures and tracer injection points have an influence on the shape of the breakthrough curve. Besides, the two-dimensional dispersion model of tracer injection of the instantaneous point was used to calculate the dispersion parameters of each group of experiments. Then, the dynamic responses of the linear distance (x) between the injection point and the receiving point, initial time (t0), peak time (tm), peak concentration (cm), average tracer transport velocity (V), and porosity (p) of aqueous media to the longitudinal dispersion coefficient are discussed. In addition, according to the measured data, Gaussian multi-peak fitting can be used to reflect the overall shape and change trend of the multi-peak BTC. These results demonstrate the solute transport behaviors in the saturated karst aquifer system, which have important reference significance for solving the engineering environmental problems in the karst area.

Including Vertical Fault Structures in Layered Groundwater Flow Models.

Accurate representation of groundwater flow and solute transport requires a sound representation of the underlying geometry of aquifers. Faults can have a significant influence on the structure and connectivity of aquifers, which may allow permeable units to connect, and aquifers to seal when juxtaposed against lower permeability units. Robust representation of groundwater flow around faults remains challenging despite the significance of faults for flow and transport. We present a methodology for the inclusion of faults utilizing the unstructured grid features of MODFLOW-USG and MODFLOW 6. The method focuses on the representation of fault geometries using non-neighbor connections between juxtaposed layers. We present an illustration of the method for a synthetic fluvial aquifer. The combined impact of the heterogeneous aquifer and fault offset is clearly visible where channel features at different depths in the aquifer were connected at the fault. These results highlight the importance of representing fault features in groundwater flow models.

Characterization of a multi-layer karst aquifer through analysis of tidal fluctuation

Interaction between ocean tides and coastal aquifers is an important process influencing the quality and quantity of groundwater resources. In this study, high frequency groundwater fluctuation along two transects from the Gulf of Mexico to the Upper Floridan Aquifer (UFA) were analyzed to characterize the hydrogeologic structure of the multi-layer karst aquifer. Tidal efficiency and time lag of harmonic constants were calculated from groundwater and seawater level and aquifer diffusivity were estimated using the Jacob-Ferris equation. Measurable tidal signals were observed at a maximum inland distance of approximately 7–8 km due to the high conductivity of the unconfined karst aquifer. The largest tidal fluctuation was observed in the deeper unit of the aquifer (146–155 m below land surface), which indicate the hydraulic contact between seawater and the high permeability zone of Avon Park Formation. The estimated values of hydraulic diffusivities agreed well with the regional groundwater flow model despite the highly heterogeneous nature of karst aquifer. Our analysis suggested that tidal springs are important sea water/groundwater interfaces and the conduit networks of tidal creeks, caves, and fractures are preferential pathways of tidal propagation.