Forecasting changes of water resources in large rivers following land drainage / Méthode de prévision des modifications des ressources en eau des grands fleuves sous l'effet du drainage

ABSTRACT

Techniques are described for annual forecasts of the water balance after drainage of large river basins. In the development of these techniques precipitation was assumed to be constant and unaffected by drainage. It is shown that the effect of drainage upon the annual runoff of the improved basins is to decrease the groundwater and swamp water resources which leads to evaporation changes. According to experimental data on the hydrophysical properties of peats, mineral soils and subsoils and how they change after drainage, the decrease in the groundwater resources was estimated for each per cent of the basin drained. This allowed account to betaken of this effect while making forecasts of runoff changes. Evaporation changes are computed as the difference between the maximum possible evaporation (potential evaporation) from cultivated areas and that from undisturbed swamps.     

Hydrological characteristics of groundwater in a subtropical coastal plain with large variations in salinity: Pimpama,Queensland, Australia

Abstract

The Pimpama coastal plain is situated in southern Moreton Bay, in subtropical eastern Australia. The plain is low lying and tidal and is situated behind a large sand barrier island. Largely due to recent (30 years) drainage networks within the flood plain, surface water quality has declined. Groundwater hydrographs have enabled the determination of different flow systems: a deeper system responding to seasonal weather patterns and a shallower flow system more responsive to individual rainfall events. Elevated potentiometric heads in semi-confined aquifers reflect upward movement of saline to hypersaline groundwaters. However, interaction of this deeper groundwater with shallower groundwater and the surface drains is yet to be determined. Recharge to the shallower system is by direct infiltration while recharge to the deeper system includes a component from landward ranges or bedrock outcrops within the plain. Discrimination between groundwater bodies is possible using salinity, ionic ratios and stable isotopes. Features of groundwater hydrology, the distribution of salinity and variations in water chemistry all suggest that under current conditions infiltration has increased, plus there is a greater landward migration of groundwaters of marine origin.     

Recharge processes during snowmelt: An isotopic and hydrometric investigation

Results from hydrometric and isotopic investigations of unsaturated flow during snowmelt are presented for a hillslope underlain by well-sorted sands. Passage of melt and rainwater through the vadose zone was detected from temporal changes in soil water ²H concentrations obtained from sequential soil cores. Bypassing flow was indicated during the initial snowmelt phase, but was confined to the near-surface zone. Recharge below this zone was via translatory flow, as meltwater inputs displaced premelt soil water. Estimates of premelt water fluxes indicate that up to 19 per cent of the premelt soil water may have been immobile. Average water particle velocities during snowmelt ranged from 6.2 × 10^?7 to 1.1 × 10^?6 ms^?1, suggesting that direct groundwater recharge by meltwater during snowmelt was confined to areas where the premelt water table was within 1 m of the ground surface. Soil water ²H signatures showed a rapid response to isotopically-heavy rain-on-snow inputs late in the melt. In addition, spatial variations in soil moisture content at a given depth induced a pronounced lateral component to the predominantly vertical transport of water. Both factors may complicate isotopic profiles in the vadose zone, and should be considered when employing environmental isotopes to infer recharge processes during snowmelt.     

SOME ASPECTS OF THE HYDRIC BALANCE OF LAKES IN THE EASTERN ROMANIAN PLAIN IN THE YEARS 1956–1970

ABSTRACT

Because of the late withdrawal of the Levantine lake waters and because of low relief the Eastern Romanian Plain was fragmented only by big alochthonous rivers (Ialomi?a, C?lm??ui and Buz?u).

The tabular-like, 40–50-km-wide interfluve areas covered by loessoid deposits and eolian sands on the periphery are deprived of surface drainage which accounts for their present evolution.

The major relief forms in these interfluves are depressions called in Romanian ‘crov’ (sink-holes) in the central areas and short valleys formed initially by erosion processes and now modelled by mechanical and chemical weathering at their periphery: in these depressions (sink-holes) and in the secondary valleys, peripheral to the interfluve areas, lakes had started to be formed.

Because of the semiarid climate sink-hole lakes have an intermittent hydrological regime, whereas those located in the small fluviatile liman-type valleys, enjoy a permanent regime. By the absence of surface drainage, by the loss of significant amounts of water through evaporation and the degree of mineralization, these lakes fall within the group of salt lakes.

In the past few years (since 1966 and especially since 1969) the level of these lakes has continually risen and the depressions formerly lacking water started being flooded by the rising of the piezometric level.

An analysis was made of the water balance of the Amara-Ialomi?a lake to investigate this phenomenon.

Level and evaporation recordings were made in the period 1956–1970. The findings revealed that the supply of underground water to the lake amounts to 47·3 per cent exceeding the water supply produced by the rains that fell on the surface of the lake (46·7 per cent). A close relationship was established (with a lapse of 8–12 months) between the surface supply of the basin (through rainfalls) and the flow of underground water to the lake.

Extending the precipitation-induced level changes over a longer period (1896–1915 and 1921–1970) it was found that level increases are cyclic, as a direct consequence of the corresponding precipitation regime.     

Interactions between groundwater and surface water in a virginia coastal plain watershed. 1. Hydrological flowpaths

A field study of surface water and groundwater interactions during baseflow and stormflow conditions was performed at the Reedy Creek watershed in the Virginia Coastal Plain. Three estimates of the average saturated hydraulic conductivity (K_s) of the unconfined aquifer were in reasonable agreement (ranging from 0.0033 to 0.010 cm/s), indicating that baseflow in the creek is entirely from the drainage of shallow groundwater from the relatively thin (1–6 m thick) unconfined aquifer. This relatively permeable surficial aquifer was found to be underlain by dark, olive grey, clay-silt and diatomaceous Miocene deposits of low permeability known as the Calvert Formation, which is believed to function as a confining bed in the area. A chemical hydrograph separation technique was used to resolve the contributions of [old] (pre-event) and [new] (event) water to stormflow. Results from a major rainstorm indicated that old water dominated the stormflow response of the watershed, although the new water contribution approached 40% at the hydrograph peak. Stormflow at Reedy Creek appears to result from saturation overland flow from variable source areas which include the stream channels and a significant part of the riparian wetland area. This response appears to be attributable to the transient dynamics of the shallow groundwater flow system and to the formation of localized groundwater mounds which raise the water-table to the wetland surface.     

Geomorphology and hydrology of upper sinking cove,cumberland plateau,Tennessee

Upper Sinking Cove, dissecting the eastern escarpment of the Cumberland Plateau, is characterized by a multiple aquifer, predominantly vadose hydrologic system with minor surface components. There is a central trunk channel along the axis of the cove and a network of independent tributaries. Aquitards within the limestones, particularly Hartselle Formation shales, have influenced both cave and surface landform development by perching ground waters and slowing the vertical growth of closed depressions. Long-term solutional denudation in the portion of the cove underlain by limestones (40 per cent) is an estimated 56 mm per 1000 years, suggesting that karst development began 15–16 million years ago. Despite lower soil CO₂ and spring water hardness, 61 per cent of annual denudation occurs in the six winter months when 76 per cent of yearly runoff occurs. Landform development in Upper Sinking Cove appears to have begun as stream erosion carved a valley first in the sandstone caprock of the escarpment and later in the underlying Pennington Formation limestones containing numerous shale layers which promoted surface stream flow. Eventually stream erosion exposed the massive Bangor limestones which allowed deep ground water flow. Surface streams were pirated underground with the eventual formation of the chain of three closed depressions which constitute Upper Sinking Cove.     

A hydrogeochemical model for stream chemistry,Cascade range,Washington, U.S.A.

A simple mixing model demonstrates that chemical variations in Cascade surface waters reflect flow from three general zones: alpine areas, forested colluvial slopes, and seasonally saturated areas. The chemistry of weathering solutions in alpine portions of the Williamson Creek catchment (North Cascade Range) results from alteration of plagioclase, hornblende, and biotite to kaolinitic material and vermiculite. Surface and shallow groundwater in forested portions of the catchment reflect these reactions, dissolution of small quantities of carbonate, and biologic activity. Both at-a-point and downstream chemical variations are explained quantitatively by the volume of water that originates in each of the hydrogeochemical source areas. Water from the forested colluvial slopes is most significant on an annual basis. However, summer low-flow is a mixture of colluvial waters and dilute solutions from the alpine zone, whereas 10 to 30 per cent of peak flow in snowmelt and rainstorms is produced from seasonally saturated areas. Poor concentration/discharge (C/Q) correlations, typical of Cascade rivers, result from mixing of significant C/Q relations for water leaving each source area. Model predictions could be substantially improved by better data for the effects of temperature, water-contact time, and biologic cycling on the chemistry of soil water from forested zones.     

Effects of Groundwater Development on Uranium: Central Valley,California, USA

Uranium (U) concentrations in groundwater in several parts of the eastern San Joaquin Valley, California, have exceeded federal and state drinking water standards during the last 20 years. The San Joaquin Valley is located within the Central Valley of California and is one of the most productive agricultural areas in the world. Increased irrigation and pumping associated with agricultural and urban development during the last 100 years have changed the chemistry and magnitude of groundwater recharge, and increased the rate of downward groundwater movement. Strong correlations between U and bicarbonate suggest that U is leached from shallow sediments by high bicarbonate water, consistent with findings of previous work in Modesto, California. Summer irrigation of crops in agricultural areas and, to lesser extent, of landscape plants and grasses in urban areas, has increased Pco₂ concentrations in the soil zone and caused higher temperature and salinity of groundwater recharge. Coupled with groundwater pumping, this process, as evidenced by increasing bicarbonate concentrations in groundwater over the last 100 years, has caused shallow, young groundwater with high U concentrations to migrate to deeper parts of the groundwater system that are tapped by public-supply wells. Continued downward migration of U-affected groundwater and expansion of urban centers into agricultural areas will likely be associated with increased U concentrations in public-supply wells. The results from this study illustrate the potential long-term effects of groundwater development and irrigation-supported agriculture on water quality in arid and semiarid regions around the world.     

Characteristics of preferential flow and groundwater discharge to Shingobee Lake,Minnesota, USA

Small‐scale heterogeneities and large changes in hydraulic gradient over short distances can create preferential groundwater flow paths that discharge to lakes. A 170 m² grid within an area of springs and seeps along the shore of Shingobee Lake, Minnesota, was intensively instrumented to characterize groundwater‐lake interaction within underlying organic‐rich soil and sandy glacial sediments. Seepage meters in the lake and piezometer nests, installed at depths of 0·5 and 1·0 m below the ground surface and lakebed, were used to estimate groundwater flow. Statistical analysis of hydraulic conductivity estimated from slug tests indicated a range from 21 to 4·8 × 10^?3 m day^?1 and small spatial correlation. Although hydraulic gradients are overall upward and toward the lake, surface water that flows onto an area about 2 m onshore results in downward flow and localized recharge. Most flow occurred within 3 m of the shore through more permeable pathways. Seepage meter and Darcy law estimates of groundwater discharge agreed well within error limits. In the small area examined, discharge decreases irregularly with distance into the lake, indicating that sediment heterogeneity plays an important role in the distribution of groundwater discharge. Temperature gradients showed some relationship to discharge, but neither temperature profiles nor specific electrical conductance could provide a more convenient method to map groundwater–lake interaction. These results suggest that site‐specific data may be needed to evaluate local water budget and to protect the water quality and quantity of discharge‐dominated lakes. Copyright © 2002 John Wiley & Sons, Ltd.     

Estimation of vertical water and heat fluxes in the semi-confined aquifers in tokyo metropolitan area,japan

Groundwater circulation is known to be one of the agents responsible for the redistribution of geothermal energy by acting as a source or sink in the course of its movement through porous media. Heat transport in groundwater systems is considered to be a coupled process and the theory based on this was used to analyse temperature profiles of 30 thermally stable observation wells in a deep, semi-confined aquifer system in the Tokyo Metropolitan area. Vertical water fluxes in the semi-confined aquifers and the associated upward heat fluxes were estimated from a heat flux equation that describes convection and conduction processes of heat transport in one dimension. The vertical downward water fluxes in Shitamachi lowland, Musashino and Tachikawa terraces were 0.69.26.91 × 10^?9, 1.46-70.92 × 10^?9 and 2.61.2204 × 10^?9 m/s, respectively. A vertical upward water flux of 1.80-33.60 × 10^?9 m/s was estimated in Shitamachi lowland. The water flux generally decreased with increasing depth for observation wells which intercepted more than one semi-confining layer. The estimated upward heat fluxes for Shitamachi lowland, Musashino and Tachikawa terraces were 0.32-1.12, 0.49-1.21 and 1.00-11.62 W/m², respectively. The heat flux was highest in Tachikawa terrace where a major fault, the Tachikawa fault, is located. Generally, the estimated heat flux was higher in the semi-confining layers than in the aquifers. Areas with heat sources and sinks as well as groundwater flow patterns in the semi-confined aquifers were revealed by heat flux and temperature distributions in the study area.     

Vertical soil water movement in a tropical rainforest catchment in northeast queensland

In a tropical rainforest catchment, shallow piezometers respond almost instantaneously to rainfall, but the dominant ground water recharge mechanisms are not well understood. To improve understanding, the downward movement of soil water on a runoff plot was traced using tritiated water injected at 0·20 m below the surface which marks the lower boundary of active subsurface storm flow. The tritium pulse was translated slowly down the profile, apparently dominated by interstitial piston flow on the lines described by Zimmermann's theoretical model. This recharge mechanism accounted for about 35 per cent of rainfall or 50 per cent of throughfall. The pulse's advance may have also been delayed by the upward movement of soil water indicated by the distribution of hydraulic potential under different hydrological conditions. The result was an increase in soil water transit time particularly below 1·0 m. There was also evidence in the tracer profiles for rapid by-pass flow but the volumes concerned could not be quantified in this experiment.     

Groundwater recharge and discharge in a hyperarid alluvial plain (Akesu,Taklimakan Desert,China)

Groundwater recharge and discharge in the Akesu alluvial plain were estimated using a water balance method. The Akesu alluvial plain (4842 km²) is an oasis located in the hyperarid Tarim River basin of central Asia. The land along the Akesu River has a long history of agricultural development and the irrigation area is highly dependent on water withdrawals from the river. We present a water balance methodology to describe (a) surface water and groundwater interaction and (b) groundwater interaction between irrigated and non‐irrigated areas. Groundwater is recharged from the irrigation system and discharged in the non‐irrigated area. Uncultivated vegetation and wetlands are supplied from groundwater in the hyperarid environment. Results show that about 90% of groundwater recharge came from canal loss and field infiltration. The groundwater flow from irrigated to non‐irrigated areas was about 70% of non‐irrigated area recharge and acted as subsurface drainage for the irrigation area. This desalinated the irrigation area and supplied water to the non‐irrigated area. Salt moved to the non‐irrigation area following subsurface drainage. We conclude that the flooding of the Akesu River is a supplemental groundwater replenishment mechanism: the river desalinates the alluvial plain by recharging fresh water in summer and draining saline regeneration water in winter. Copyright © 2007 John Wiley & Sons, Ltd.     

Preferential groundwater flow through a sorted net landscape,Arctic Canada

Preferential suprapermafrost groundwater flow was observed in deepened channels lying between raised frost mounds. Here, saturated hydraulic conductivity, k, ranged from 90 to 1000 m/day but was only 0·1–1·0 m/day in the mound centres. A high proportion of fines occurs in the frost mound centre due to particle sorting, while channels contain gravels. Three approaches of areal weighting of k and groundwater flow, Q _s, across a wetland–upland boundary were explored. When percentage area covered by channels, mounds or gravel was considered, estimates of water flow on a daily and seasonal basis fell by 30 to 50 per cent. This study is of relevance to northern scientists who require reliable estimates of groundwater flow across patterned ground landscapes. Copyright © 2001 John Wiley & Sons, Ltd.     

Imaging Hydrological Processes in Headwater Riparian Seeps with Time‐Lapse Electrical Resistivity

Delineating hydrologic and pedogenic factors influencing groundwater flow in riparian zones is central in understanding pathways of water and nutrient transport. In this study, we combined two‐dimensional time‐lapse electrical resistivity imaging (ERI) (depth of investigation approximately 2 m) with hydrometric monitoring to examine hydrological processes in the riparian area of FD‐36, a small (0.4 km²) agricultural headwater basin in the Valley and Ridge region of east‐central Pennsylvania. We selected two contrasting study sites, including a seep with groundwater discharge and an adjacent area lacking such seepage. Both sites were underlain by a fragipan at 0.6 m. We then monitored changes in electrical resistivity, shallow groundwater, and nitrate‐N concentrations as a series of storms transitioned the landscape from dry to wet conditions. Time‐lapse ERI revealed different resistivity patterns between seep and non‐seep areas during the study period. Notably, the seep displayed strong resistivity reductions (～60%) along a vertically aligned region of the soil profile, which coincided with strong upward hydraulic gradients recorded in a grid of nested piezometers (0.2‐ and 0.6‐m depth). These patterns suggested a hydraulic connection between the seep and the nitrate‐rich shallow groundwater system below the fragipan, which enabled groundwater and associated nitrate‐N to discharge through the fragipan to the surface. In contrast, time‐lapse ERI indicated no such connections in the non‐seep area, with infiltrated rainwater presumably perched above the fragipan. Results highlight the value of pairing time‐lapse ERI with hydrometric and water quality monitoring to illuminate possible groundwater and nutrient flow pathways to seeps in headwater riparian areas.     

F) Personalia

Abstract

Knowledge of the natural resources of a country is essential for the assessment of its development possibilities. This is particularly true for the evaluation of groundwater resources in a country like Uruguay, where less than 5 per cent of the potential is exploited and this in an inefficient way, and where no drastic future increase in demand is predicted.

In the present paper the generalized hydrogeological model of the country was reconstructed, and recharge values, obtained by an approach adapted from nearby areas, were introduced.     

Defining hydrochemical evolution of streamflow through flowpath dynamics in Kawakami headwater catchment,Central Japan

The hydrochemical behaviour of catchments is often investigated by inferring stream chemistry through identification of source areas involved in hydrograph separation analysis, yet its dynamic evolution of hydrologic pathways has received little attention. Intensive hydrometric and hydrochemical measurements were performed during two different storms on March 29, 2001 and August 21–22, 2001 to define hydrochemical evolution under the dynamic of flow pathways in a 5·2 ha first‐order drainage of the Kawakami experimental basin (KEB), Central Japan, a forested headwater catchment with various soil depths (1·8 to 5 m) overlying late Neogene of volcanic bedrocks. The hydraulic potential distribution and flow lines data showed that the change in flow direction, which was controlled by rainfall amount and antecedent wetness of the soil profile, agreed well with the hydrochemical change across the slope segment during the storm. Hydrograph separation predicted by end‐member mixing analysis (EMMA) using Ca²⁺ and SiO₂ showed that near surface riparian, hillslope soil water and deep riparian groundwater were important in stream flow generation. The evidence of decrease in solutes concentration at a depth of 1 m in the hillslope and 0·6 m in the near surface riparian during peak storm suggested a flushing of high solutes concentration. Most of the solutes accumulated in the deep riparian groundwater zone, which was due to prominent downward flow and agreed well with the residence time. The distinct flow pathways and chemistry between the near surface riparian and deep riparian groundwater zones and the linkage hillslope aquifer and near surface riparian reservoir, which controls rapid flow and solutes flushing during the storm event, are in conflict with the typical assumption that the whole riparian zone resets flow pathways and chemical signature of hillslope soil water, as has been reported in a previous study. Copyright © 2005 John Wiley & Sons, Ltd.     

Effect of bedrock flow on catchment rainfall‐runoff characteristics and the water balance in forested catchments in Tanzawa Mountains,Japan

In this paper, we examined the role of bedrock groundwater discharge and recharge on the water balance and runoff characteristics in forested headwater catchments. Using rigorous observations of catchment precipitation, discharge and streamwater chemistry, we quantified net bedrock flow rates and contributions to streamwater runoff and the water balance in three forested catchments (second‐order to third‐order catchments) underlain by uniform bedrock in Japan. We found that annual rainfall in 2010 was 3130 mm. In the same period, annual discharge in the three catchments varied from 1800 to 3900 mm/year. Annual net bedrock flow rates estimated by the chloride mass balance method at each catchment ranged from ?1600 to 700 mm/year. The net bedrock flow rates were substantially different in the second‐order and third‐order catchments. During baseflow, discharge from the three catchments was significantly different; conversely, peak flows during large storm events and direct runoff ratios were not significantly different. These results suggest that differences in baseflow discharge rates, which are affected by bedrock flow and intercatchment groundwater transfer, result in the differences in water balance among the catchments. This study also suggests that in these second‐order to third‐order catchments, the drainage area during baseflow varies because of differences between the bedrock drainage area and surface drainage area, but that the effective drainage area during storm flow approaches the surface drainage area. Copyright © 2012 John Wiley & Sons, Ltd.     

Groundwater age rejuvenation caused by excessive urban pumping in Jakarta area,Indonesia

Rapid urbanization in the Jakarta area has become a severe subsurface environmental issue as it entails groundwater level decline and land subsidence caused by excessive groundwater pumping. In this study, apparent groundwater age rejuvenation in the deep aquifer under DKI Jakarta was found by comparing ¹⁴C activities between 1985 and 2008. We discussed the use of a numerical groundwater flow model to evaluate the rejuvenation process in this urbanized area. When considering the deep aquifer in the DKI Jakarta area, we can assume six direction fluxes toward the aquifer: two vertical fluxes (downward and upward flux) and four horizontal fluxes (northern, southern, western, and eastern flux). Results of model calculations show that the greatest groundwater flux among six flux directions became ‘vertical downward flux’, which means that shallower groundwater intrudes into the deep one because of excessive groundwater pumping from the mid‐1980s. This flux grows about 50% during the 2000s. This result is consistent with the detection of CFC‐12 and SF₆, which functions as an indicator of young groundwater even in the deep groundwater. The rejuvenation ratio ‘R’ was determined using ¹⁴C activity in the groundwater; R increases with the CFC‐12 concentration and both show good correlation. Furthermore, we estimated the ‘vertical downward flux’ at each well's screen depth using model estimation. Results show that this flux is greater in the urban groundwater depression area and especially at shallower parts of the deep aquifer, and that it affects the magnitude of the shallow groundwater intrusion. Copyright © 2012 John Wiley & Sons, Ltd.     

Groundwater potential of the Egbe-Mopa basement area,central Nigeria

Abstract

An investigation on the groundwater potentials of the Egbe-Mopa area in central Nigeria, underlain by the Basement Complex, is presented. The investigation involved mapping of the subsurface by use of vertical electrical soundings; measurement of depth to groundwater; and evaluation of hydraulic conductivity, transmissivity and yield by means of pumping test interpretation. The results indicate subsurface units that range from three to five resistivity layers; depth to groundwater of 0–10 m; overburden thickness of 3–16 m; hydraulic conductivity of 6.2?×?10^?6 to 3.4?×?10^?4 m/s; transmissivity of 4.3?×?10^?7 to 2?×?10^?3 m²/s; and groundwater yield of 0.2–2.5 L/s. The hydraulic head assessments revealed a general northward groundwater flow direction. The study identified three aquifer potential types, of high, medium and low productivity, respectively. Based on the longitudinal conductance of the overburden units, four distinct Aquifer Protective Capacity zones were delineated, namely, poor, weak, moderate and good.

Citation Okogbue, C.O. and Omonona, O.V., 2013. Groundwater potential of Egbe-Mopa basement area, central Nigeria. Hydrological Sciences Journal, 58 (4), 826–840.     

Toward a new conceptual model for groundwater flow in merokarst systems: Insights from multiple geophysical approaches

Merokarst aquifers — relatively thin (<1–2 m) karstified carbonate units interbedded between mudstone, shale, or sandstone — constitute a significant proportion of carbonate terrain and underlie a large portion of the west- and south-central USA, yet few advances have been made in our understanding of porosity development and flow-path generation in these complex systems in decades. Toward this end, we used a multi-geophysical approach at the well-studied Konza Prairie Biological Station (KPBS), a part of the larger Flint Hills (25,734 km²), underlain by thin limestone units (1–2 m thick) interbedded with mudstone/shale units (2–4 m thick), to elucidate hydrologic connectivity and potential controls on known groundwater flow directions. We combined electrical resistivity tomography (ERT), surface and borehole nuclear magnetic resonance (NMR), and ground penetrating radar (GPR) measurements across a low order catchment where over 25 boreholes and groundwater wells sampling perched aquifers could be used to constrain interpretation of lithology, potential flow paths, and permeability. Data revealed that groundwater export may be an unappreciated component of lateral-flow-dominated models used to represent merokarst in that: (a) potentiometric surfaces from two limestone units showed groundwater flows toward a hydrologic depression, opposite the direction of stream flow, in the upstream portion of the catchment, (b) long term measures of groundwater levels revealed a greater variance and overall water storage in this same upstream area compared to wells near the outlet, and (c) ERT and NMR results indicate pronounced lateral heterogeneity within a given unit, suggestive of a greater degree of vertical hydrological connectivity than usually considered for horizontally-layered merokarst. Our data suggest vertical connectivity can shunt water to depth in these “sandwiched” merokarst aquifers, yielding atypical groundwater flow directions and unrealized deep export of weathering solutes and carbon.