Formation of phreatic caves in an eogenetic karst aquifer by CO2 enrichment at lower water tables and subsequent flooding by sea level rise

Formation of extensive phreatic caves in eogenetic karst aquifers is widely believed to require mixing of fresh and saltwater. Extensive phreatic caves also occur, however, in eogenetic karst aquifers where fresh and saltwater do not mix, for example in the upper Floridan aquifer. These caves are thought to have formed in their modern settings by dissolution from sinking streams or by convergence of groundwater flow paths on springs. Alternatively, these caves have been hypothesized to have formed at lower water tables during sea level low‐stands. These hypotheses have not previously been tested against one another. Analyzing morphological data and water chemistry from caves in the Suwannee River Basin in north‐central Florida and water chemistry from wells in the central Florida carbonate platform indicates that phreatic caves within the Suwannee River Basin most likely formed at lower water tables during lower sea levels. Consideration of the hydrological and geochemical constraints posed by the upper Floridan aquifer leads to the conclusion that cave formation was most likely driven by dissolution of vadose CO₂ gas into the groundwater. Sea level rise and a wetter climate during the mid‐Holocene lifted the water table above the elevation of the caves and placed the caves tens of meters below the modern water table. When rising water tables reached the land surface, surface streams formed. Incision of surface streams breached the pre‐existing caves to form modern springs, which provide access to the phreatic caves. Phreatic caves in the Suwannee River Basin are thus relict and have no causal relationship with modern surficial drainage systems. Neither mixing dissolution nor sinking streams are necessary to form laterally extensive phreatic caves in eogenetic karst aquifers. Dissolution at water tables, potentially driven by vadose CO₂ gas, offers an underappreciated mechanism to form cavernous porosity in eogenetic carbonate rocks. Copyright © 2012 John Wiley & Sons, Ltd.     

Impact of grout curtains on karst groundwater behaviour: an example from the Dinaric karst

Grout curtains are vertical grout walls installed in the ground. In karst terrains, their construction is primarily connected with dams and reservoirs. Their main role is to increase water tightness and to prevent progressive erosion, blocking possible seepage paths along karst fissures and conduits. In this article, changes in the behaviour of the groundwater level (GWL) and the water temperature in nine deep piezometers, which were caused by the construction of a grout curtain at the ?ale Reservoir on the Cetina River (Croatia), were analysed. The total length of the grout curtain is 3966 m. It spreads 120 m below the dam. The most analysed data are from the period after the dam had been built. Only few data and figures concern the comparison between pre‐ and post‐dam periods. The hourly data of the GWL and the water temperature were analysed for the period between 1 September 2008 at 02:00 h to 31 December 2009 at 23:00 h (11 687 h total) in six deep piezometers (marked in the text and figures as 1, 2, 3, 4, 5 and 6). For three piezometers (marked in the text and figures as A, B and C), some discontinuous measurements of the GWL and the water temperature were available for analysis. The construction of the grout curtain made strong, sudden and possibly dangerous changes to the characteristics of the aquifer and the circulation of groundwater in the local area. Special attention is paid to analyses of the behaviour of the hourly GWL data measured in the piezometers pairs (two neighbouring piezometers, one inside and the other outside of the grout curtain). During more than 80% of the analysed period, the GWL was higher in the piezometer inside the grout curtain than the one outside of it. The intensity and range of the dynamics of GWL was higher in piezometer outside the grout curtain than the inside ones. After the construction of the grout curtain, the maximum measured hydrostatic pressure on some parts of the grout curtain was approximately 40 m. It changes quickly in both time and direction. The water temperature was found to be similar in all of the measured piezometers, and it varies between 10.2 and 15.7 °C with an average value of 12.7 °C. Copyright © 2011 John Wiley & Sons, Ltd.     

Characterizing the heterogeneity of karst critical zone and its hydrological function: An integrated approach

Spatial heterogeneity in the subsurface of karst environments is high, as evidenced by the multiphase porosity of carbonate rocks and complex landform features that result in marked variability of hydrological processes in space and time. This includes complex exchange of various flows (e.g., fast conduit flows and slow fracture flows) in different locations. Here, we integrate various “state‐of‐the‐art” methods to understand the structure and function of this poorly constrained critical zone environment. Geophysical, hydrometric, and tracer tools are used to characterize the hydrological functions of the cockpit karst critical zone in the small catchment of Chenqi, Guizhou Province, China. Geophysical surveys, using electrical resistivity tomography (ERT), inferred the spatial heterogeneity of permeability in the epikarst and underlying aquifer. Water tables in depression wells in valley bottom areas, as well as discharge from springs on steeper hillslopes and at the catchment outlet, showed different hydrodynamic responses to storm event rainwater recharge and hillslope flows. Tracer studies using water temperatures and stable water isotopes (δD and δ¹⁸O) could be used alongside insights into aquifer permeability from ERT surveys to explain site‐ and depth‐dependent variability in the groundwater response in terms of the degree to which “new” water from storm rainfall recharges and mixes with “old” pre‐event water in karst aquifers. This integrated approach reveals spatial structure in the karst critical zone and provides a conceptual framework of hydrological functions across spatial and temporal scales.     

Simulation of spring flows from a karst aquifer with an artificial neural network

In China, 9·5% of the landmass is karst terrain and of that 47,000 km² is located in semiarid regions. In these regions the karst aquifers feed many large karst springs within basins of thousands of square kilometres. Spring discharges reflect the fluctuation of ground water level and variability of ground water storage in the basins. However, karst aquifers are highly heterogeneous and monitoring data are sparse in these regions. Therefore, for sustainable utilization and conservation of karst ground water it is necessary to simulate the spring flows to acquire better understanding of karst hydrological processes. The purpose of this study is to develop a parsimonious model that accurately simulates spring discharges using an artificial neural network (ANN) model. The karst spring aquifer was treated as a non‐linear input/output system to simulate the response of karst spring flow to precipitation and applied the model to the Niangziguan Springs, located in the east of Shanxi Province, China and a representative of karst springs in a semiarid area. Moreover, the ANN model was compared with a previous time‐lag linear model and it was found that the ANN model performed better. Copyright © 2007 John Wiley & Sons, Ltd.     

Locating the zone of saline intrusion in a coastal karst aquifer using springflow data

Coastal fresh water aquifers are an increasingly desirable resource. In a karstic aquifer, sea water intrusion occurs as a salt water wedge, like in porous media. However, preferential flow conduits may alter the spatial and temporal distribution of the salt water. This is typically the case when the outlet of the aquifer is a brackish spring. This paper shows that salinity and flow rate variations at a spring, where salinity is inversely proportional to discharge, can help to understand the hydrodynamic functioning of the aquifer and to locate the fresh water-sea water mixing zone deep inside the aquifer. The volume of water-filled conduit between the sea water intrusion zone and the spring outlet is calculated by the integral over time of the flow rate during the time lag between the flow rate increase and the salinity decrease as measured at the spring. In the example of the spring at Almyros of Heraklio (Crete, Greece), this time lag is variable, depending on the discharge, but the volume of water-filled conduit appears to be constant, which shows that the processes of salt water intrusion and mixing in the conduit are constant throughout the year. The distance between the spring and the zone where sea water enters the conduit is estimated and provides an indication of the position where only fresh water is present in the conduit.     

Hydrogeological investigation of karst system properties by common use of diverse methods: a case study of Lička Jesenica springs in Dinaric karst of Croatia

Li?ka Jesenica is a sinking river situated in the Dinaric karst environment of the Lika region of Croatia. The two main karst springs, Veliko Vrelo and Malo Vrelo, provide the majority of Li?ka Jesenica's water. Because of the quality and abundance of the water they contain, these springs are strategically important to the public water supply. Previous knowledge regarding the springs` karst system has been negligible. Therefore, a bulk hydrogeological research program was conducted with the purposes of establishing protection zones and assessing the springs' water‐supply potential. The research work comprised the following methods: hydrogeological mapping, setting up of monitoring stations for evaluating discharge and water quality, time series analyses (hydrograph–chemograph), hydrochemical analyses, and geophysical surveys. The main results and conclusions of the research include the following: (i) confirmation that both springs drain a common aquifer; (ii) determination of the highly vulnerable nature of the karst systems; (iii) identification of the Veliko Vrelo Spring's more direct connection to the aquifer‐drainage system in comparison to the Malo Vrelo; (iv) estimation of the positions and depths of the springs` main feeding conduits; (v) confirmation of the extent of the geologically presumed catchment area; and (vi) finally, establishment of sanitary protection zones. Copyright © 2012 John Wiley & Sons, Ltd.     

Changes in hydrological regime caused by human intervention in karst: the case of the Rumin Springs

ABSTRACT

Any human intervention or action in karst terrains can unexpectedly, suddenly, strongly and, generally, dangerously change a local and/or a regional hydrological regime. A characteristic example of the Dinaric karst is given in this paper. The operation of two reservoirs in Livanjsko Polje at an altitude of about 702 m above sea level (m a.s.l.) and hydro-electric development of the Cetina River system started in 1973. This year marked a drastic and instantaneous change in the regional hydrological regime. A significant drop in the minimum, mean and maximum annual discharges of two neighbouring karst springs, Rumin Mali and Rumin Veliki, was caused by this anthropogenic construction. The exits of the two analysed karst springs, Rumin Mali and Veliki, are located at altitudes of about 326.8 m a.s.l. and 307.6 m a.s.l. respectively. The areal distance between them is 640 m. Their hydrological regimes have been altered in different ways. The drop in discharges is stronger for the Rumin Veliki than for the Rumin Mali. The analysis shows that the system development influenced the redistribution of the regional karst aquifer and, in this manner, contributed to the reduction of both karst springs’ recharge areas.

Editor D. Koutsoyiannis Associate editor E. Rozos     

Using Alternate Pumping and Cooperative Game Theory to Reduce Sea Water Intrusion

In this article, alternate pumping is studied as a means used to reduce the salinity concentration in coastal aquifers, pumped using a system of wells. This approach has been applied to a hypothetical confined coastal aquifer. Flow has been modeled, using SEAWAT. Two strategies are proposed based on cooperative game theory, to promote alternate pumping. In both strategies an external player will compensate the users that will pump during an unpopular pumping period. In the first strategy it is supposed that this external player aims at protecting a critical well, e.g. a municipal well, reducing its maximum salinity concentration by pumping alternately. In the second strategy proposed, the target is to reduce the overall salinity of the water pumped by the wells. In applying the cooperative game theory, the Shapley value is used to distribute the benefits of cooperation between the players (well users), according to their marginal contribution. Overall, well users can reduce sea water intrusion by cooperatively changing their pumping time schedules. The game theoretical model developed is a useful tool to promote cooperation toward this direction. The methods applied in the hypothetical aquifer, can be tested in actual aquifers to reduce sea water intrusion.     

Dissolution in a variably confined carbonate platform: effects of allogenic runoff,hydraulic damming of groundwater inputs,and surface–groundwater exchange at the basin scale

In variably confined carbonate platforms, impermeable confining units collect rainfall over large areas and deliver runoff to rivers or conduits in unconfined portions of platforms. Runoff can increase river stage or conduit heads in unconfined portions of platforms faster than local infiltration of rainfall can increase groundwater heads, causing hydraulic gradients between rivers, conduits and the aquifer to reverse. Gradient reversals cause flood waters to flow from rivers and conduits into the aquifer where they can dissolve limestone. Previous work on impacts of gradient reversals on dissolution has primarily emphasized individual caves and little research has been conducted at basin scales. To address this gap in knowledge, we used legacy data to assess how a gradient of aquifer confinement across the Suwannee River Basin, north‐central Florida affected locations, magnitudes and processes of dissolution during 2005–2007, a period with extreme ranges of discharge. During intense rain events, runoff from the confining unit increased river stage above groundwater heads in unconfined portions of the platform, hydraulically damming inputs of groundwater along a 200 km reach of river. Hydraulic damming allowed allogenic runoff with SI_CAL < ?4 to fill the entire river channel and flow into the aquifer via reversing springs. Storage of runoff in the aquifer decreased peak river discharges downstream and contributed to dissolution within the aquifer. Temporary storage of allogenic runoff in karst aquifers represents hyporheic exchange at a scale that is larger than found in streams flowing over non‐karst aquifers because conduits in karst aquifers extend the area available for exchange beyond river beds deep into aquifers. Post‐depositional porosity in variably confined carbonate platforms should thus be enhanced along rivers that originate on confining units. This distribution should be considered in models of porosity distribution used to manage water and hydrocarbon resources in carbonate rocks. Copyright © 2013 John Wiley & Sons, Ltd.     

A high‐altitude temporary spring in a compartmentalized carbonate aquifer: the role of low‐permeability faults and karst conduits

The aim of this research was to refine the actual conceptual model related to the activation of high‐altitude temporary springs within the carbonate Apennines in southern Italy. The research was carried out through geophysical, hydrogeological, hydrochemical and isotopic investigations at the Acqua dei Faggi experimental site during five hydrologic years. The research demonstrated that, in carbonate aquifers where low‐permeability faults cause the aquifer system to be compartmentalized, high‐altitude temporary springs may be recharged by groundwater. In such settings, neither surface water infiltration in karst systems nor perched temporary aquifers play a role of utmost importance. The rare (once or a few time a year) activation of such springs is due to the fact that groundwater unusually reach the threshold head that allows the spring to flow. The activation of the studied high‐altitude temporary spring also depended on relationships between a low‐permeability fault core and a karst system that locally interrupts the low‐permeability barrier. In fact, when the hydraulic head did not reach the karst system, the concentrated head loss within the fault core did not allow the spring to flow, because the groundwater entirely flowed through the fault towards the downgradient compartment. Copyright © 2009 John Wiley & Sons, Ltd.     

Application of an Analytical Solution as a Screening Tool for Sea Water Intrusion

Sea water intrusion into aquifers is problematic in many coastal areas. The physics and chemistry of this issue are complex, and sea water intrusion remains challenging to quantify. Simple assessment tools like analytical models offer advantages of rapid application, but their applicability to field situations is unclear. This study examines the reliability of a popular sharp‐interface analytical approach for estimating the extent of sea water in a homogeneous coastal aquifer subjected to pumping and regional flow effects and under steady‐state conditions. The analytical model is tested against observations from Canada, the United States, and Australia to assess its utility as an initial approximation of sea water extent for the purposes of rapid groundwater management decision making. The occurrence of sea water intrusion resulting in increased salinity at pumping wells was correctly predicted in approximately 60% of cases. Application of a correction to account for dispersion did not markedly improve the results. Failure of the analytical model to provide correct predictions can be attributed to mismatches between its simplifying assumptions and more complex field settings. The best results occurred where the toe of the salt water wedge is expected to be the closest to the coast under predevelopment conditions. Predictions were the poorest for aquifers where the salt water wedge was expected to extend further inland under predevelopment conditions and was therefore more dispersive prior to pumping. Sharp‐interface solutions remain useful tools to screen for the vulnerability of coastal aquifers to sea water intrusion, although the significant sources of uncertainty identified in this study require careful consideration to avoid misinterpreting sharp‐interface results.     

Effects of Recharge Wells and Flow Barriers on Seawater Intrusion

The installation of recharge wells and subsurface flow barriers are among several strategies proposed to control seawater intrusion on coastal groundwater systems. In this study, we performed laboratory‐scale experiments and numerical simulations to determine the effects of the location and application of recharge wells, and of the location and penetration depth of flow barriers, on controlling seawater intrusion in unconfined coastal aquifers. We also compared the experimental results with existing analytical solutions. Our results showed that more effective saltwater repulsion is achieved when the recharge water is injected at the toe of the saltwater wedge. Point injection yields about the same repulsion compared with line injection from a screened well for the same recharge rate. Results for flow barriers showed that more effective saltwater repulsion is achieved with deeper barrier penetration and with barriers located closer to the coast. When the flow barrier is installed inland from the original toe position however, saltwater intrusion increases with deeper barrier penetration. Saltwater repulsion due to flow barrier installation was found to be linearly related to horizontal barrier location and a polynomial function of the barrier penetration depth.     

Hydrogeochemical characteristics of surface water and groundwater in the karst basin,southwest China

Groundwater is a very significant water source used for irrigation and drinking purposes in the karst region, and therefore understanding the hydrogeochemistry of karst water is extremely important. Surface water and groundwater were collected, and major chemical compositions and environmental isotopes in the water were measured in order to reveal the geochemical processes affecting water quality in the Gaoping karst basin, southwest China. Dominated by Ca²⁺, Mg²⁺, HCO₃^? and SO₄^2?, the groundwater is typically characterized by Ca? Mg? HCO₃ type in a shallow aquifer, and Ca? Mg? SO₄ type in a deeper aquifer. Dissolution of dolomite aquifer with gypsiferous rocks and dedolomitization in karst aquifers are important processes for chemical compositions of water in the study basin, and produce water with increased Mg²⁺, Ca²⁺ and SO₄^2? concentrations, and also increased TDS in surface water and groundwater. Mg²⁺/Ca²⁺ molar ratios in groundwater decrease slightly due to dedolomitization, while the mixing of discharge of groundwater with high Mg²⁺/Ca²⁺ ratios may be responsible for Mg²⁺/Ca²⁺ ratios obviously increasing in surface water, and Mg²⁺/Ca²⁺ ratios in both surface water and groundwater finally tending to a constant. In combination with environmental isotopic analyses, the major mechanism responsible for the water chemistry and its geochemical evolution in the study basin can be revealed as being mainly from the water–rock interaction in karst aquifers, the agricultural irrigation and its infiltration, the mixing of surface water and groundwater and the water movement along faults and joints in the karst basin. Copyright © 2009 John Wiley & Sons, Ltd.     

The role of antecedent groundwater heads in controlling transient aquifer storage and flood peak attenuation in karst watersheds

Transient storage of floodwaters in aquifers is known to attenuate peak flows in rivers and drive subsurface dissolution. Transient aquifer storage could be enhanced in watersheds overlying karst aquifers where caves facilitate surface and groundwater exchange. Few studies, however, have examined controls on, or magnitudes of, transient aquifer storage or flood peak attenuation in karstic watersheds. Here we evaluate flood peak attenuation with multiple linear regression analyses of 10 years of river and groundwater data from the Suwannee River, which flows over the karstic upper Floridan aquifer in north-central Florida and experiences frequent flooding. Regressions show antecedent river stage exerts the dominant control on magnitudes of transient aquifer storage, with recharge and time to peak having secondary controls. Specifically, low antecedent stages result in larger magnitudes of transient aquifer storage and thus greater flood attenuation than conditions of elevated antecedent stage. These findings suggest subsurface weathering, including cave formation and enlargement, caused by transient aquifer storage could occur on a more frequent basis in aquifers where groundwater table elevation is lowered due to anthropogenic or climatic influences. Our work also shows that measures of groundwater table elevation prior to an event could be used to improve predictive flood models. © 2018 John Wiley & Sons, Ltd.     

Characterizing moldic and vuggy pore space in karst aquifers using borehole-wall,slabbed-core and thin-section images

Carbonate aquifers are prolific and important sources of potable water in many parts of the world owing to enlarged dissolution features that enhance porosity and interconnectivity. To better understand the variations of pore space in different karst aquifers, image and geospatial analyses are used to analyze pore attributes (i.e., pore area and perimeter) in images of vuggy aquifers. Pore geometry and 2D porosity derived from images of the moldic Castle Hayne and vuggy Biscayne aquifers are analyzed at three scales of observation: borehole televiewer, core and thin-section. The Castle Hayne and Biscayne aquifers are the foci of this study because the pore spaces that control the hydrologic properties in each of these aquifers are markedly different even though both of these carbonate reservoirs are prolific aquifers. Assessments of pore area, perimeter and shape index (a measure of shape complexity) indicate that pore geometries and pore complexities vary as a function pore type and scale of observation. For each aquifer type, the areas, perimeters and complexities of pores are higher at the larger scale of observation (e.g., borehole) than the smaller scale of observation (e.g., thin section). When the complexity of the moldic pores is compared to the complexity of vuggy pores, the results indicate that moldic pores are generally more complex than vuggy pores at the same scale of observation. Whereas estimates of 2D porosity from the borehole televiewer image of the vuggy aquifer are higher than those derived from the moldic aquifer, the range of 2D porosities is larger in core and thin section images for the vuggy aquifer than the moldic aquifer. A model for the development of pores is presented that suggests that the coalescence of small pores with simple shapes leads to the growth of larger pores with more complex shapes. The model suggests that the younger Biscayne aquifer is a more mature karst than the Castle Hayne aquifer.     

Improved understanding of dynamic water and mass budgets of high-alpine karst systems obtained from studying a well-defined catchment area

Large areas of Europe, especially in the Alps, are covered by carbonate rocks and in many alpine regions, karst springs are important sources for drinking water supply. Because of their high variability and heterogeneity, the understanding of the hydrogeological functioning of karst aquifers is of particular importance for their protection and utilisation. Climate change and heavy rainfall events are major challenges in managing alpine karst aquifers which possess an enormous potential for future drinking water supply. In this study, we present research from a high-alpine karst system in the UNESCO Biosphere Reserve Großes Walsertal in Austria, which has a clearly defined catchment and is drained by only one spring system. Results show that (a) the investigated system is a highly dynamic karst aquifer with distinct reactions to rainfall events in discharge and electrical conductivity; (b) the estimated transient atmospheric CO₂ sink is about 270 t/a; (c) the calculated carbonate rock denudation rate is between 23 and 47 mm/1000a and (d) the rainfall-discharge behaviour and the internal flow dynamics can be successfully simulated using the modelling package KarstMod. The modelling results indicate the relevance of matrix storage in determining the discharge behaviour of the spring, particularly during low-flow periods. This research and the consequent results can contribute and initiate a better understanding and management of alpine karst aquifers considering climate change with more heavy rainfall events and also longer dry periods.     

Mass transport triggered by heavy rainfall: the role of endorheic basins and epikarst in a regional karst aquifer

This study aims at recognizing the mechanisms of mass transport between the karst surface and the saturated zone in a morphostructural relief of the Mesozoic karst carbonate platform of Murgia (Puglia, Southern Italy). The large dimension of the karst aquifer, the regional scale of the flow system, the boundary condition constituted by the sea, and the lack of freshwater springs constrain to the use of wells as monitoring points and limit the study area to the recharge area comprising 986 endorheic basins. The concentrations of non‐reactive tracers (nitrates) in the waters of autogenic recharge (from endorheic basins) have been modeled through the evaluation of effective infiltration, land use, and nitrogen surplus, with reference to a time window, which includes a low precipitation period followed by significant rainfall events. The comparison between the modeled nitrate concentrations and the nitrate concentrations measured in ground waters, coupled with the analysis of groundwater chemograms and records of hydraulic heads (all referred to the same time window), allows inferring the mechanism of mass transport between the karst surface and the groundwater table. The mass transport conceptual model requires the presence of the epikarst. The infiltration of significant rainfall in the endorheic basins after a low precipitation period displaces waters stored in the epikarst toward the saturated zone. Ground waters in the post‐event period show higher concentrations of nitrates, lower concentrations of total organic carbon, and higher Mg/Ca ratios than both those of the pre‐event period and the autumn‐winter recharge period. The post‐event recharge from epikarst storage determines a transient hazard of groundwater pollution with a time lag from the occurrence of the heavy rainfall.     

The impact of aquifer stratification on saltwater intrusion characteristics. Comprehensive laboratory and numerical study

Laboratory experiments and numerical simulations were utilized in this study to assess the impact of aquifer stratification on saltwater intrusion. Three homogeneous and six layered aquifers were investigated. Image processing algorithms facilitated the precise calculation of saltwater wedge toe length, width of the mixing zone, and angle of intrusion. It was concluded that the length of intrusion in stratified aquifers is predominantly a function of permeability contrast, total aquifer transmissivity and the number of heterogeneous layers, being positively correlated to all three. When a lower permeability layer overlays or underlays more permeable zones its mixing zone widens, while it becomes thinner for the higher permeability strata. The change in the width of the mixing zone (WMZ) is positively correlated to permeability contrast, while it applies to all strata irrespectively of their relative vertical position in the aquifer. Variations in the applied hydraulic head causes the transient widening of WMZ. These peak WMZ values are larger during saltwater retreat and are negatively correlated to the layer's permeability and distance from the aquifer's bottom. Moreover, steeper angles of intrusion are observed in cases where low permeability layers overlay more permeable strata, and milder ones in the inverse aquifer setups. The presence of a low permeability upper layer results in the confinement of the saltwater wedge in the lower part of the stratified aquifer. This occurs until a critical hydraulic head difference is applied to the system. This hydraulic gradient value was found to be a function of layer width and permeability contrast alike.     

Identifying Aquifer Type in Fractured Rock Aquifers using Harmonic Analysis

Determining aquifer type, unconfined, semi‐confined, or confined, by drilling or performing pumping tests has inherent problems (i.e., cost and complex field issues) while sometimes yielding inconclusive results. An improved method to cost‐effectively determine aquifer type would be beneficial for hydraulic mapping of complex aquifer systems like fractured rock aquifers. Earth tides are known to influence water levels in wells penetrating confined aquifers or unconfined thick, low‐porosity aquifers. Water‐level fluctuations in wells tapping confined and unconfined aquifers are also influenced by changes in barometric pressure. Harmonic analyses of water‐level fluctuations of a thick (～1000 m) carbonate aquifer located in south‐central Oklahoma (Arbuckle‐Simpson aquifer) were utilized in nine wells to identify aquifer type by evaluating the influence of earth tides and barometric‐pressure variations using signal identification. On the basis of the results, portions of the aquifer responded hydraulically as each type of aquifer even though there was no significant variation in lithostratigraphy. The aquifer type was depth dependent with confined conditions becoming more prevalent with depth. The results demonstrate that harmonic analysis is an accurate and low‐cost method to determine aquifer type.