Introduction

Declining groundwater levels are a global concern, as a large part of the human population depends on fresh water from aquifers (UNESCO 2009). Van Loon et al. (2016) showed that humans are an important factor in the causes of droughts and the depletion of reservoirs. The depletion of reservoirs is mainly a result of groundwater pumping for drinking-water supply (Willis and Garrod 1998), irrigation (Amelung et al. 1999; Foster et al. 2004; Konikow and Kendy 2005; Hoque et al. 2007; Qureshi et al. 2010; Wada et al. 2012) and industrial use (Hayashi et al. 2009). Globally, groundwater withdrawal amounts to 750–800 km3 year−1 and exceeds annual recharge in many places of the world (Shah et al. 2000).

Urgent action is needed to slow down the depletion of groundwater reservoirs. Currently, managed aquifer recharge (MAR) methods are available to recover depleted aquifers (Dillon 2005; Maliva et al. 2014). Examples are known from many aquifers in the world—e.g. Burdekin Delta and Angas-Bremer area in Australia (Gerges et al. 2002; Dillon 2009); Satlasana and Kodangipalayam in India (Gale et al. 2006); and the Columbia River Basalt Group aquifer in USA (Eaton et al. 2009). However, in order to improve the success of MAR projects (Gale et al. 2006), advanced aquifer characterization is needed to identify suitable locations for MAR practices (Dillon 2005; Rahman et al. 2012; Maliva et al. 2014). Basaltic aquifers which are widely used for water abstraction face water conflicts between irrigation demand, drinking-water supply and sustainable aquifer management (Macdonald et al. 1995; El-Naqa et al. 2007; Pavelic et al. 2012)—for example on the basaltic Deccan Traps in India, irrigation by groundwater has developed from shallow dug wells to mechanical abstraction since the 1950s. In the Maharashtra state in India alone, the total abstraction volume from more than 1.5 million wells has increased by over 700 % since then (Macdonald et al. 1995).

One of the areas where severe groundwater depletion exists is the Moscow-Pullman region on the border of Idaho and Washington State in the USA (Fig. 1). Here, groundwater decline started with exploitation of upper aquifers in the early 20th century (Laney et al. 1923). The Moscow sub-basin is part of the Palouse basin, which is located on the eastern part of the Columbia River Basalt Group (CRBG; Hooper 1982). The CRBG aquifer is a main source of fresh water in the area (Tolan et al. 1989). Despite curtailed groundwater pumping since 1992, aquifer levels have continued to decline at 40–60 cm year−1 (Robischon 2006, 2007). The Moscow sub-basin boundary is characterized by eolian loess deposits overlying eastward-dipping basalts to the west with granitic highlands and Miocene/Pliocene sediments to the east. These unique geologic zones interfinger in the subsurface (Fig. 2). Groundwater levels in the shallow basaltic aquifer (Wanapum Formation) recover in response to reduced pumping, indicating that recharge occurs to the shallow aquifer (Leek 2006). The presence of tritium in water samples recently taken from the Wanapum Formation and in some wells in the upper portion of the deeper basaltic aquifer (Grande Ronde), suggest that both these aquifers received recharge to some extent over the last 50 years (Carey 2011).

Fig. 1
figure 1

a The Moscow sub-basin (study area), located in the Palouse Basin, on the border of Idaho/Washington state, USA. b The Palouse Basin, subdivided in the North Fork and South Fork of the Palouse River. The red line shows the Moscow sub-basin Figures 1 and 3 contains poor quality of text. Otherwise, please provide replacement figure file.correct

Fig. 2
figure 2

Geological cross-section of the Moscow sub-basin (after Bush and Garwood 2005). Letters AC and dashed arrows correspond to the recharge pathways mentioned in the text. Cross-section is based on approximate W–E dashed line A–A′ in Fig. 4

Identifying groundwater-recharge-source areas is critical to local managers and municipalities for targeting management which preserves and potentially enhances aquifer recharge. Previous work suggested that no vertical recharge occurs in the Moscow sub-basin (arrow A in Fig. 2). This was concluded based on the presence of thick impermeable clay deposits (up to 30 m) that were identified through deep well logs (Lum et al. 1990; Bush 2005; Fairley et al. 2006). Rather than vertical percolation, lateral saturated flow occurs on top of these impermeable clays (Brooks et al. 2004). Using a spatially distributed hydrologic analysis approach, Dijksma et al. (2011) showed that it is possible that recharge occurs only in certain regions. They, and Bush (2005), suggested that the greatest potential for recharge might be via spatially discrete areas of high permeability (e.g. paleo-channels) draining from forested regions on the Moscow Mountain range front, recharging upper and lower aquifers at the granite/basalt interface (arrow B and C in Fig. 2). In addition, Dijksma et al. (2011) and Fairley et al. (2006) found evidence in boreholes of the existence of thick layers of coarse material, of which some end up at the surface as perennial springs (De Graaf 2011). Using a numeric groundwater model, De Graaf (2011) confirmed that aquifer recharge through paleo-channels could provide a relatively large part of the aquifer recharge that occurs in the Paradise Creek watershed. It is very likely that paleo-channels are present and provide lateral conduits of groundwater recharge to the deeper aquifers (Brooks and Grader 2011). One of the remaining challenges is the identification of these groundwater recharge source areas in order to use these lateral conduits in MAR applications.

The aim of this study is to identify these groundwater recharge pathways in the Moscow sub-basin. A soil moisture routing (SMR) model (Dijksma et al. 2011) was used first, to provide a spatial representation of percolation below the root zone across the entire basin. Secondly, stable isotopes δ18O and δ2H as well as deuterium excess (d-excess) were used to identify linkages between the surface water and groundwater systems. δ18O and δ2H have widely been recognized as useful tracers in providing insights into water movements in watersheds (McDonnell and Kendall 1992; Kendall and MacDonnell 1998). Seasonal and event-based variability in δ18O and δ2H in rain and snowmelt provide a unique composition which can be traced in streamflow and observed in groundwater wells. δ18O and δ2H tracers have been used in the Moscow sub-basin to describe variability in groundwater sources; however, since only one sample was taken from each well at different moments in time (Larson et al. 2000; Carey 2011; Moxley 2012) these studies were not able to investigate temporal variability in δ18O and δ2H. If a hydrologic connection exists between two different sources of water (e.g. groundwater, precipitation or surface water), then it is expected that the down-gradient water source would have a similar but damped δ18O and δ2H that is lagged in time relative to up-gradient sources (Changnon 1987; Soulsby and Tetzlaff 2008; Katsuyama et al. 2010; Speed et al. 2010; Wassenaar et al. 2011). If a hydrologic connection between different sources of water exists, temporal fluctuations in δ18O and δ2H will be hydrologically similar (Scanlon et al. 2002).

In addition to δ18O fluctuations, the δ18O can be compared to the Global Meteoric Water Line (GMWL: δ2H = 8 × δ18O + 10), that serves as a reference to determine deviations of the relation between δ18O and δ2H (Craig 1961). The offset of 10 is determined by kinetic isotope fractionation that occurs during non-equilibrium processes such as evaporation (Craig 1961; Cappa et al. 2003). This offset can be calculated by the d-excess (d = δ2H − 8 × δ18O), and is a measure for the relative proportions of δ2H and δ18O that are affected by evaporation (Dansgaard 1964). The correlation with precipitation is higher when the d-excess of a water source is close to 10 ‰. A d-excess close to 10 ‰ also indicates that precipitation has infiltrated and percolated below the root zone quickly, and has not been subject to evaporation (Froehlich et al. 2001). Sánchez-Murillo et al. (2015) measured δ18O and δ2H in precipitation (N samples = 203) and streamflow (N = 244 and N = 195, in Crumarine Creek and the South Fork of the Palouse River, respectively) from June 2011 to January 2014 (Fig. 3). Their record provides the isotopic baseline necessary to investigate hydrologic connections to groundwater.

Fig. 3
figure 3

Local meteoric water lines for various water sources collected in this study in comparison to the global meteoric water line (GMWL). Crumarine Creek and precipitation is according to Sánchez-Murillo et al. 2015

Area description

Figure 4a shows the digital elevation model for the Moscow sub-basin, including the drainages from Moscow Mountain (part of the east–west trending Palouse Range). These drainages include Paradise Creek and the South Fork of the Palouse River. Discharge measurements by the US Geological Survey (USGS) of Paradise Creek near the University of Idaho Campus from 1979 to 2012 indicate the mean annual discharge of 0.12 mm/year. The South Fork of the Palouse River at Pullman (Washington) has a mean annual discharge of 0.40 mm/year, which includes the Paradise Creek discharge (US Geological Survey 2012). Non-irrigated agriculture is the dominant land use in the Moscow sub-basin (71 %); forest (15 %), grass fields (6 %), and urban areas (8 %) are the other land use types present (Fig. 4b; Bara and Shaw 1995).

Fig. 4
figure 4

a Digital elevation model (DEM) of the Moscow sub-basin. Black line indicates the extent of the CRBG basalt flows, defined by Grader (2011). The dashed line A–A′ represents the geological cross-section as shown in Fig. 2, b Land use map of the Moscow sub-basin

Pullman-Moscow groundwater systems are fairly complex with multiple basalt flows underlain by and laterally resting against protruding basement metamorphic/granitic ridges. Basin margin paleo-valleys occur at different scales, and are stacked and altered through geologic time by successive basalt dam effects and are interwoven with clay and sand/gravel lacustrine and stream/ swamp inter-beds (Grader 2011). High-porosity interbeds (associated with subordinate paleo-channels) are known to occur at lower stratigraphic levels (Grader 2011).

The extent of the upper CRBG basalt flows (Wanapum Basalt) is indicated in Fig. 4. Underlying older, thicker Grande Ronde flows occur westward of this transitional line (shown in cross-section, Fig. 2). Associated aquifers use the same names “upper Wanapum” or “lower Grande Ronde” formations, but involve heterogeneous geology. The CRBG basalt flows are inter-fingered with Miocene and younger sediments (Latah Formation including Sediments of Moscow, Vantage Member, and Sediments of Bovill), and are a result of alluvial deposition when basalt flows blocked former drainage patterns (Fig. 2). A detailed Quaternary surficial geological map is dominated by eolian deposited loess silts (Palouse Formation) as shown in Fig. 5. Arable stoneless soils show that no ice age tillites were deposited at this latitude; however, major volcanic airfall ashes derived from the Cascade Mountains were admixed into both Miocene and Quaternary sediments. Lower through upper Miocene sediments are a mixture of all textures from clay to sand and gravel (Lum et al. 1990; Bush 2005; Fairley et al. 2006). This research focuses on recharge in this part of the basin, where paleo-channels of different orders of size should be contained by and will laterally migrate within the larger-scaled subsurface paleo-valleys.

Fig. 5
figure 5

Surficial geological map of the Moscow sub-basin including the well locations, well sample IDs (160) and well depth (Othberg et al. 2001)

The Wanapum and Grande Ronde aquifers are deformed, and tectonic features or feeder dikes may be present resulting in isolated groundwater reservoirs with low connectivity (McVay 2007). Locally open horizontal and vertical fractures are present and are, at present, clay- or even sand-filled. These fractures have not been extensively studied and are likely to vary spatially (Fairley et al. 2006). Unpredictable rock/sediment relationships are the result of invasive basalts that flowed over wet sediments as observed in ancient lakes or fluvial streams along paleo-valley boundaries. Previous attempts to predict groundwater flow in the Wanapum and Grande Ronde aquifers by numerical modelling (Barker 1979; Lum et al. 1990; Hansen et al. 1994; Whiteman et al. 1994; Vaccaro 1999; Reidel et al. 2003) were often inconsistent due to the combination of both aquifer’s hydrogeological heterogeneity and the complex groundwater extractions from different locations and intervals (Leek 2006; McVay 2007; Bennett 2009).

Material and methods

In this study, the hydrologic connections between precipitation, streamflow, and groundwater are investigated by testing groundwater in 22 wells and 2 springs throughout the Moscow sub-basin (Fig. 5; Table 1). These locations have been selected based on a preliminary sampling of 49 wells and 3 springs (Candel 2014). Well information was derived from the online database of the Idaho Department of Water Resources (IDWR; Idaho Department of Water Resources 2013). When well information was not present in the online database, it was derived from the well owners. Each of the wells and springs were sampled once every 2 weeks. In total, 352 samples have been analysed for δ18O and δ2H. Wells were selected based on geology, well depth (Fig. 4b), and proximity to streams, with the intention to have distributed samples throughout the Moscow sub-basin (Idaho Department of Water Resources 2013). Water samples are taken directly from outdoor spigots after flushing the system of old water. Water lines were flushed until the water temperature was constant. Water samples were collected in 30-ml glass E-C borosilicate bottles with TFE-lined caps (Wheaton Science Products, USA) and stored upside down at 5 °C. Bottles were filled with no head space and covered with parafilm (Thermo Scientific, USA) to avoid exchange with atmospheric moisture.

Table 1 Characteristics of all wells and springs that were sampled

Stable isotope analyses were conducted at the Idaho Stable Isotope Laboratory, University of Idaho (Idaho, USA) using a cavity ring-down spectroscopy (CRDS) water isotope analyser L1120-i (Picarro, USA). Laboratory standards, previously calibrated to the reference waters Vienna Standard Mean Ocean Water and Standard Light Antarctic Precipitation (VSMOW2-SLAP2), were glacier water (δ2H = −255.0 ‰, δ18O = −30.8 ‰), commercial bottled water (δ2H = −64.2 ‰, δ18O = −8.3 ‰), and Moscow tap water (δ2H = −125.5 ‰, δ18O = −15.4 ‰). Glacier water and commercial bottled water were used to normalize the results to the VSMOW2-SLAP2 scale, while the Moscow tap water was a quality control standard. The laboratory instrument precision on average was ±0.5 ‰ (1σ) for δ2H and ±0.1 ‰ (1σ) for δ18O. The estimated d-excess analytical uncertainty was ±0.6 ‰ (1σ).

A check to the suitability of using outdoor spigots was completed by quantifying changes in δ18O with time at five spigots. The temporal variation in δ18O after 2, 5 and 10 min of flushing was always smaller than the standard error of the Picarro Instrument, indicating that differences found in isotopic composition cannot be explained by different water being sampled within the pipelines of the wells.

Soil-moisture routing model

In this study, a SMR model was used to provide spatial predictions of the mean annual percolation rates below the root zone throughout the Moscow sub-basin. The SMR model used is a spatially distributed, grid-based hydrologic model which operates within a geographic information system (GIS) environment, originally developed at Cornell University (Frankenberger et al. 1999). The model simulated interception, evapotranspiration, subsurface lateral flow, deep vertical percolation and saturation-excess surface runoff through a multi-layer soil profile. Snow accumulation and melt is simulated using an energy balance approach (Brooks and Boll 2005; Brooks et al. 2007). The SMR model has been developed as a spatially explicit management tool and therefore has been developed to rely primarily on publicly available data with minimal calibration (Frankenberger et al. 1999; Brooks et al. 2007). It uses spatial explicit topographic, land cover, and soil maps to represent the hydrologic mass balance within the root zone of a watershed. It is particularly well suited in landscapes having restrictive soil horizons (e.g. argillic and fragipan) where saturation excess runoff is generated by the subsurface lateral redistribution of water following variable source area hydrology concepts (Brooks et al. 2004, 2007; McDaniel et al. 2008). Despite the modest input requirements, modelling results are at least as good as more complex hydrology models (Johnson et al. 2003; Mehta et al. 2004; Dijksma et al. 2011). Brooks et al. (2007) demonstrated the ability of the model to accurately represent the development of shallow, dynamic perched water tables over restrictive fragipan soil horizons, including snow accumulation and melt, and surface runoff, from a small grassland catchment on the eastern edge of the Moscow sub-basin.

Dijksma et al. (2011) applied and assessed the ability of the SMR model to represent the hydrology of the Paradise Creek watershed. The Paradise Creek watershed is also located in the Moscow sub-basin (Fig. 2), and covers approximately one third of the total (150 km2) area. The Nash-Sutcliffe efficiency (NSE; Nash and Suthcliffe 1970) between observed and simulated streamflow from Paradise creek watershed was 0.57, which is a “good” result according to the qualitative assessment of NSE by Foglia et al. (2009). More details can be found in Dijksma et al. (2011).

Since the focus of this study was on the entire Moscow sub-basin, the spatial extent of the Dijksma et al. (2011) model extended to cover the entire region. The model was extended using the same topographic, land cover, soil maps and parameters as had been used as by Dijksma et al. (2011). The daily average percolation was calculated by SMR for each 30 × 30 m cell in the Moscow sub-basin, over the period 2001–2008.

Results

Model

As seen in Fig. 6, the SMR model indicates that percolation rates vary spatially with the greatest average percolation rate of 2 mm d−1, whereas there are other areas that have negligible percolation. The overall average daily percolation over the Moscow sub-basin was 0.50 mm d−1.

Fig. 6
figure 6

Daily average percolation in the Moscow sub-basin in m d−1. A result from the SMR model

Isotopes

Table 2 shows the results of the δ18O analysis for all wells and springs. The δ18O values for the wells varied from −16.7 ‰ in well 23 to −14.1 ‰ in well 39. The temporal fluctuations in δ18O varied in magnitude at each of the sampled wells and springs. The variability in δ18O was significant and consistent across many of the locations. Figure 7a shows highly dynamic δ18O fluctuations at four wells that had the highest standard deviation in δ18O (wells 9, 14, 24 and 27). In contrast, δ18O in some other wells is relatively stable. Figure 7b shows relatively stable δ18O fluctuations at four wells that had the lowest standard deviation in δ18O (wells 11, 35, 36, 37). The standard deviation correlates with the maximum range in δ18O values (R 2 = 0.91). Figure 8b shows the spatial distribution of the maximum range in δ18O in all wells and springs. Similarly, Fig. 8c shows the spatial distribution of the average d-excess in all wells and springs. Interestingly the fluctuations in δ18O indicated by Fig. 8b, c do not necessarily correlate with well depth (Fig. 8a). For example, deep wells (e.g. wells 35, 29, 9, 14) can have high fluctuations in δ18O and a high d-excess and shallow wells can have stable fluctuations and a low d-excess (e.g. wells 1, 41, 36).

Table 2 Results of water samples from all wells and springs, from 17-04-2013 until 02-01-2014. SD standard deviation
Fig. 7
figure 7

a δ18O readings for the four most dynamic wells in the study, including the precipitation δ18O, b δ18O readings for the four most stable wells in the study, including the precipitation δ18O

Fig. 8
figure 8

Well locations plotted on the DEM. Numbers indicate the well sample IDs. a the depth of the wells, b the maximum magnitude of the δ18O composition, c the average d-excess for all samples, d indications of recharge

Discussion

Recharge in the Moscow sub-basin

According to the SMR model, most of the percolation occurs at higher elevations near Moscow Mountain. The relatively large percolation at the higher elevations can be partially attributed to the fact that Moscow Mountain receives nearly twice the precipitation (1,270 mm) as the city of Moscow located at relatively low elevation in the study area (WRCC 2013). The model predicts zero percolation in areas having strong hydraulically restrictive argillic soil characteristics. Instead of water moving vertically through these horizons, perched water tables develop in these soils and the water then will run laterally downslope, producing runoff at toe-slope positions. Soils having low percolation exist throughout a large part of the eastern Palouse region between Moscow and the Moscow Mountain, where thick clay layers are present (Lum et al. 1990; Fairley et al. 2006; Dijksma et al. 2011). Percolation is also limited within the city of Moscow due to impervious surfaces which route water through storm drain networks.

The δ18O data indicate that connections exist between precipitation, surface water and groundwater. It is suggested here that there are two characteristics in the δ18O data that indicate recharge. Firstly, the magnitude of δ18O fluctuations, whereby if fluctuations are large, then the groundwater is hydrologically connected to a source of water causing these variations in δ18O. This different source can either be groundwater recharge from precipitation or surface water. As seen in Fig. 7, there is considerable overlap in the range of δ18O values in wells exhibiting both relatively large and small fluctuations, which suggests that no distinction in recharge can be made based on the average annual δ18O; hence, the time record is needed for this distinction. Interestingly, wells that have small δ18O fluctuations still have a dampened δ18O, suggesting a slower connection or smaller amount of recharge may be occurring in these wells. The second characteristic that suggests that recharge is occurring is if the d-excess is close to 10 ‰. If water percolation occurs rapidly without excessive evaporation, then the groundwater will have a d-excess close to 10 ‰.

In order to examine the spatial patterns in recharge, three ranges were defined for each of these recharge indicators varying from no indications of recharge to strong indications of recharge. These categories are shown in Fig. 8b, c with different colors. Yellow, orange and red represent that the indications for recharge are absent, present and strong, respectively. Each of these categories was given a value of 0 to 2, respectively. Figure 8d shows the sum of these values for both recharge indicators in similar categories ranging from yellow (0–1), orange (2–3) and red (4). This method, which assumes that both recharge indicators are equally important, provides a clear overview of potential groundwater recharge source areas in the Moscow sub-basin.

In order to provide a perspective of the importance on well depth and location relative to specific geologic features and groundwater recharge, a geological cross-section is provided of the well locations across the Moscow sub-basin (see Fig. 9). The well locations are simplified and schematically presented in the 2D cross-section. The color of each well represents the same recharge potential as used in Fig. 8d. In order to properly interpret these data, it is important to recognize that well logs indicate most of these wells pull groundwater from coarse interbeds below poorly permeable (e.g. thick clay) layers (Idaho Department of Water Resources 2013). The wells showing the strongest indications of recharge (35, 11, 29 and 14) are all located along the granite/basalt interface near Moscow Mountain and vary largely in depth (Figs. 8d and 9). The location of these wells agrees well with the simulated percolation maps from the SMR model. Further west towards Moscow, the indications of recharge decline, although both shallow springs exhibit strong indications of recharge (6, 19). These springs are likely supplied by upstream recharge. Some wells show no indications of recharge (1, 22, 23, 36). These wells are all relatively shallow and get water from the Sediments of Bovill, which lay above the Wanapum. As most of the water is recharged at the granite/basalt interface (arrow B and C in Fig. 2), a relatively small amount of water recharges the top of the Wanapum where these wells are located. The lack of recharge indications in these wells, combined with the SMR output, also underlines that little to no vertical percolation and recharge pathways occur in this region (arrow A in Fig. 2). Moreover, wells located relatively far from the recharging area do not exhibit any indications of recharge (see Wanapum and Vantage Member wells 20, 41 and 60); however, some Wanapum wells (39, 5) have indications of recharge, explained by the heterogeneity of the Wanapum aquifer (Leek 2006; McVay 2007; Bennett 2009). Overall, the isotope and modelling data both suggest that groundwater recharge is occurring at the basalt/granite interface near the forested uplands in the eastern fringe of the Moscow sub-basin. The well data suggest that this groundwater recharge flows laterally below the poorly permeable argillic horizons and thick clay deposits typically found in the lower agricultural production regions in the Moscow Sub-basin. This lateral flow likely occurs through complex coarse interbed layers or conduits at the granite/basalt interface where these geological units interfinger. Some of these lateral flow paths end up as surficial springs (Dijksma et al. 2011), whereas others likely recharge both aquifers. Other than the forested uplands, the modelling and isotope data indicate that little vertical percolation occurs over much of the Moscow sub-basin. This agrees well with the Fairley et al. (2006) groundwater characterization study which concluded little vertical recharge occurs over most of the Moscow sub-basin due to the presence of widespread, thick clay deposits; however, this study provides added insight that the groundwater recharge is likely occurring as lateral flow beneath these clay deposits along the granite-basalt interface.

Fig. 9
figure 9

Geological cross-section of the Moscow sub-basin (after Bush and Garwood 2005). All wells and springs have been included. The well locations are simplified and schematic. The numbers corresponds to the well sample ID. The well colours represent the same as in Fig. 8d. Red is strong indications of recharge, orange is indications of recharge, yellow is no indications of recharge. Cross-section is based on approximate W–E dashed line A–A’ in Fig. 4

Temporally the δ18O peaks in groundwater are lagged and damped relative to δ18O peaks in precipitation. Although the dataset is relatively short, Fig. 7a indicates there is a 2–5-week lag between a δ18O peak in precipitation and a δ18O peak in groundwater. This lagged δ18O suggests that the recharge response time may be as short as 2–5 weeks. A longer more complete dataset would be required to determine a characteristic recharge response time for specific regions in the basin.

Wider implications

This study shows that permeable linear elements on the interface between impermeable (granitic) bedrock and permeable (basalt) formations can act as recharge pathways. Inducing recharge by artificial infiltration ponds may help to reduce the water level decline or even cause groundwater level rise in such basaltic aquifers. This principle also can be applied to other basaltic regions such as Deccan Traps. Isotope fractioning can help finding suitable locations for infiltration ponds.

Conclusion

Biweekly isotope measurements of 22 wells and two springs combined with a soil moisture routing (SMR) model showed strong indications of recharge in the proximity of the Moscow Mountain. The SMR model results suggest that precipitation percolates mainly at the Moscow Mountain, where maximum daily average percolation rates are reached (up to 2.0 mm d−1). The Moscow Mountain provides the source of aquifer recharge, from where water flows through lateral conduits below the poorly permeable layers recharging the basalt aquifers at the granite/basalt interface where these geological units interfinger. This study shows the importance of granite/basalt interface areas for recharge to the basalt aquifers. These locations may be the main sources from where these types of aquifers receive recharge. As the local economy and many individuals depend on groundwater from basalt aquifers, these interface areas should be further explored towards the possibilities of using managed aquifer recharge and aquifer storage methods.