Proglacial erosion rates and processes in a glacierized catchment in the Swiss Alps

In the Swiss Alps, climatic changes have not only caused glacier retreat, but also likely increased sedimentation downstream of glaciers. This material either originates from below the glacier or from periglacial environments, which are exposed as glaciers retreat, and often consist of easily erodible sediment. Griesgletscher's catchment in the Swiss Alps was examined to quantify erosion in the proglacial area, possible hydrological drivers and contributions of the sub‐ and periglacial sources. Digital elevation models, created from annual aerial photographs, were subtracted to determine annual volume changes in the proglacial area from 1986 to 2014. These data show a strong increase in proglacial erosion in the decade prior to 2012, coincident with increasing proglacial area size. However, examination of the gradient between discharge and sediment evacuation, and modeled sediment transport, could suggest that the proglacial area began to stabilize and sediment supply is limited. The large influx of sediment into the proglacial reservoir, which is roughly 2.5 times greater than the amount of sediment eroded from the proglacial area, demonstrates the importance of subglacial erosion to the catchment's sediment budget. Although far more sediment originates subglacially, erosion rates in the proglacial area are over 50 times greater than the rest of the catchment. In turn, both sub‐ and periglacial processes, in addition to constraining sediment supply, must be considered for assessing future sediment dynamics as glacier area shrinks and proglacial areas grow. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Interrogating the time and processes of development of the Las Liebres rock glacier,central Chilean Andes,using a numerical flow model

The Las Liebres rock glacier is a large (~2.2 km long) Andean rock glacier whose internal composition and kinematics are known from previous studies. We investigate its development by posing and testing the following null hypothesis: the rock glacier has developed from a constant supply of debris and ground ice in periglacial conditions and resulting creep of the ice‐rock mixture. A rheological model was formulated based on recent advances in the study of ice‐rock mixture rheology, and calibrated on the known surface velocities and internal composition of the rock glacier. We show that the rock glacier viscosity is inversely related to both water and debris fractions, in agreement with recent field and theoretical studies of ice‐rock mixture rheology. Taking into account the possible variations in water fraction, the model was used to estimate the time spans of development (0.91–7.11 ka), rates of rock wall retreat (0.44–4.18 mm/a), and rates of ground ice formation (0.004–0.026 m/a) for the rock glacier. These results support the null hypothesis of a periglacial origin of the Las Liebres rock glacier. Copyright © 2016 John Wiley & Sons, Ltd.     

Mountain permafrost dynamics within a recently exposed glacier forefield inferred by a combined geomorphological,geophysical and photogrammetrical approach

Geomorphological observations, geoelectrical soundings and photogrammetric measurements of surface movement on the Muragl glacier forefield were used to obtain an integrative analysis of a highly complex glacial and periglacial landform consisting of a push moraine, creeping permafrost and permafrost‐free glacial till in close proximity. Electrical resistivity tomography is considered as an important multifunctional geophysical method for research in periglacial permafrost related environments. Joint application with measurements of surface displacements offers a promising tool for investigating periglacial landforms related to ice‐rich permafrost for a more comprehensive characterization of permafrost characteristics and geomorphological interpretation of periglacial morphodynamics. The patchy permafrost distribution pattern described in this paper is determined by several factors, including the sediment characteristics, the snow cover distribution and duration, the aspect and the former glacier distribution and thermal regime. Recent and modern permafrost dynamics within the glacier forefield comprise aggradation, degradation and permafrost creep. Copyright © 2007 John Wiley & Sons, Ltd.     

A multi‐method field experiment to determine local groundwater flow in a glacier forefield

We implemented multiple independent field techniques to determine the direction and velocity of groundwater flow at a specific stream reach in a glacier forefield. Time‐lapse experiments were conducted using two electrical resistivity tomography (ERT) lines installed in a cross pattern. A circular array of groundwater tubes was also installed to monitor groundwater flow via discrete salt injections. Both inter‐borehole and ERT results confirmed this stream section as a losing reach and enabled quantification of the flow direction. Both techniques yielded advection velocities varying between 5.7 and 21.8 m/day. Estimates of groundwater flow direction and velocity indicated that groundwater infiltrates from the stream nearby and not from the adjacent lateral moraine. Groundwater age estimated from radon concentration measurements supported this hypothesis. Despite uncertainties inherent to each of the methods deployed, the combination of multiple field techniques allowed drawing consistent conclusions about local groundwater flow. We thus regard our multi‐method approach as a reliable way to characterize the two‐dimensional groundwater flow at sites where more invasive groundwater investigation techniques are difficult to carry out and local heterogeneities can make single measurements unreliable. Copyright © 2014 John Wiley & Sons, Ltd.     

Hydraulic subsurface measurements and hydrodynamic modelling as indicators for groundwater flow systems in the Rotondo granite,Central Alps (Switzerland)

Regional groundwater flow in high mountainous terrain is governed by a multitude of factors such as geology, topography, recharge conditions, structural elements such as fracturation and regional fault zones as well as man‐made underground structures. By means of a numerical groundwater flow model, we consider the impact of deep underground tunnels and of an idealized major fault zone on the groundwater flow systems within the fractured Rotondo granite. The position of the free groundwater table as response to the above subsurface structures and, in particular, with regard to the influence of spatial distributed groundwater recharge rates is addressed. The model results show significant unsaturated zones below the mountain ridges in the study area with a thickness of up to several hundred metres. The subsurface galleries are shown to have a strong effect on the head distribution in the model domain, causing locally a reversal of natural head gradients. With respect to the position of the catchment areas to the tunnel and the corresponding type of recharge source for the tunnel inflows (i.e. glaciers or recent precipitation), as well as water table elevation, the influence of spatial distributed recharge rates is compared to uniform recharge rates. Water table elevations below the well exposed high‐relief mountain ridges are observed to be more sensitive to changes in groundwater recharge rates and permeability than below ridges with less topographic relief. In the conceptual framework of the numerical simulations, the model fault zone has less influence on the groundwater table position, but more importantly acts as fast flow path for recharge from glaciated areas towards the subsurface galleries. This is in agreement with a previous study, where the imprint of glacial recharge was observed in the environmental isotope composition of groundwater sampled in the subsurface galleries. Copyright © 2012 John Wiley & Sons, Ltd.     

Analysis and interpretation of patterns within and between hydroclimatological time series in an alpine glacier basin

This paper explores patterns within and between climatological and hydrological time series from an alpine glacier basin. Time series recorded in the basin of the Haut Glacier d'Arolla over the 1989 ablation season are subdivided into five subperiods. Box-Jenkins ARIMA (AutoRegressive Integrated Moving Average) and TFN (Transfer Function-Noise) models are estimated for each of the five subperiods and differences between the models are interpreted in the context of changing glacier hydrology, particularly the changing nature and extent of the glacier drainage network.     

Influence of a rock glacier spring on the stream energy budget and cold‐water refuge in an alpine stream

The thermal regimes of alpine streams remain understudied and have important implications for cold‐water fish habitat, which is expected to decline due to climatic warming. Previous research has focused on the effects of distributed energy fluxes and meltwater from snowpacks and glaciers on the temperature of mountain streams. This study presents the effects of the groundwater spring discharge from an inactive rock glacier containing little ground ice on the temperature of an alpine stream. Rock glaciers are coarse blocky landforms that are ubiquitous in alpine environments and typically exhibit low groundwater discharge temperatures and resilience to climatic warming. Water temperature data indicate that the rock glacier spring cools the stream by an average of 3 °C during July and August and reduces maximum daily temperatures by an average of 5 °C during the peak temperature period of the first two weeks in August, producing a cold‐water refuge downstream of the spring. The distributed stream surface and streambed energy fluxes are calculated for the reach along the toe of the rock glacier, and solar radiation dominates the distributed stream energy budget. The lateral advective heat flux generated by the rock glacier spring is compared to the distributed energy fluxes over the study reach, and the spring advective heat flux is the dominant control on stream temperature at the reach scale. This study highlights the potential for coarse blocky landforms to generate climatically resilient cold‐water refuges in alpine streams.     

Groundwater dynamics in a till hillslope: flow directions,gradients and delay

Knowledge of groundwater dynamics is important for the understanding of hydrological controls on chemical processes along the water flow pathways. To increase our knowledge of groundwater dynamics in areas with shallow groundwater, the groundwater dynamics along a hillslope were studied in a boreal catchment in Southern Sweden. The forested hillslope had a 1‐ to 2‐m deep layer of sandy till above bedrock. The groundwater flow direction and slope were calculated under the assumption that the flow followed the slope of the groundwater table, which was computed for different triangles, each defined by three groundwater wells. The flow direction showed considerable variations over time, with a maximum variation of 75°. During periods of high groundwater levels the flow was almost perpendicular to the stream, but as the groundwater level fell, the flow direction became gradually more parallel to the stream, directed in the downstream direction. These findings are of importance for the interpretation of results from hillslope transects, where the flow direction usually is assumed to be invariable and always in the direction of the hillslope. The variations in the groundwater flow direction may also cause an apparent dispersion for groundwater‐based transport. In contrast to findings in several other studies, the groundwater level was most responsive to rainfall and snowmelt in the upper part of the hillslope, while the lower parts of the slope reached their highest groundwater level up to 40 h after the upper parts. This can be explained by the topography with a wetter hollow area in the upper part. Copyright © 2011 John Wiley & Sons, Ltd.     

Improved understanding of spring and stream water responses in headwaters of the Indian Lesser Himalaya using stable isotopes,conductivity and temperature as tracers

Stable isotope data are presented for precipitation, spring and stream water in a headwater catchments in the Indian Lesser Himalaya. Isotopic contents of phreatic groundwater followed the local meteoric water line and showed minimal alteration by evaporation, suggesting fast recharge. Mean isotopic values for springs and the stream were close to the weighted annual mean for precipitation, indicating recharge was in synchrony with seasonal rainfall distribution. Precipitation exhibited isotopic declines of ?0.6‰ and ?0.2‰ δ¹⁸O per 100 m rise in elevation in July and August (monsoon), respectively. The time lag of one month between rainfall and spring discharge, combined with the isotopic lapse rate indicated a recharge elevation of 70–165 m above the spring outflow point, implying the water originated within the catchment. Time series of electrical conductivity and temperature of spring, seepage and stream waters confirmed the rapid recharge and limited storage capacity of the shallow aquifers.     

Seasonal evolution of meltwater generation,storage and discharge at a non‐temperate glacier in Svalbard

In glacierized catchments, meteorological inputs driving surface melting are translated into runoff outputs mediated by the glacier hydrological system: analysis of the relationship between meteorology and diurnal and seasonal patterns of runoff should reflect the functioning of that system, with the role of meltwater storage likely to be of particular importance. Daily meltwater storage is determined for a glacier at 78 °N in the Svalbard archipelago, by comparing inputs calculated from a surface energy balance model with measured outputs (proglacial discharge). Solar radiation, air temperature, wind speed and proglacial discharge are then analysed by regression and time‐series methods, in order to assess the meteorology–discharge relationship and its variation at diurnal and seasonal time‐scales. The recorded discharge time‐series can be divided into two contrasting intervals: up to early August, proglacial discharge was high and variable, mean hydrographs showed little indication of diurnal cycling, ARIMA models of discharge indicated a non‐seasonal, moving‐average generating process, and there was a net loss of meltwater from storage; from early August, proglacial discharge was low and relatively invariable, but with clearer diurnal cycles, regression models of discharge showed substantially improved correlations with air temperature and solar radiation, ARIMA models indicated a non‐seasonal, autoregressive generating process, and eventually a seasonal component, and there was a net gain in meltwater storage. The transition between the two periods is brief compared with the duration of the melt season. The runoff response to meteorology therefore lacks the strongly progressive element previously identified in mid‐latitude glacierized catchments. In particular, the glacier hydrological system only appears responsive to diurnal forcing following the depletion of the seasonal snowpack meltwater store. Copyright © 2001 John Wiley & Sons, Ltd.     

Temporal and spatial variability of cation and silica export in an alpine watershed,Emerald Lake,California

A reaction set of possible mineral weathering reactions is proposed to explain observed cation and silica export for the Emerald Lake watershed, a small Sierra Nevada, California catchment. The reaction set was calculated through a stoichiometric mole‐balance method, using a multiyear record of stream flow and snowpack chemical analyses and site‐specific mineral compositions. Reaction‐set calculations were intended to explore how the processes controlling stream cation and silica export depend on differing bedrock mineralogy across the catchment as snowmelt and runoff patterns change over the year. Different regions within the watershed can be differentiated by lake inflow subdrainages, each exhibiting different stream‐flow chemistry and calculated weathering stoichiometry, indicating that different silica and cation generation processes are dominant in wet steep portions of the catchment. Short‐term differences in stream concentrations were assumed to reflect ion exchange equilibria and rapid biological processes, whereas long‐term persistent stream concentration differences in different areas of the catchment were assumed to reflect spatial variability in mineral weathering stoichiometry. Mineralogical analyses of rock samples from the watershed provided site‐specific chemical compositions of major mineral species for reaction calculations. Reaction sets were evaluated by linear regression of calculated versus observed differences between snowmelt and stream‐flow chemistry and by a combined measure. Initially, single weathering reactions were balanced and evaluated to determine the reactions that best explained observed stream chemical export. Next, reactions were combined, using mineral compositions from different rock types to estimate the dependence of ion fluxes on lithology. The seasonal variability of major solute calculated fluxes is low, approximately one order of magnitude, relative to the observed three orders of magnitude variability in basin discharge. Reaction sets using basin‐averaged lithology and Aplite lithologies gave superior explanations of stream chemical composition. Copyright © 2004 John Wiley & Sons, Ltd.     

Tracer‐based evidence of heterogeneity in subsurface flow and storage within a boreal hillslope

Runoff from boreal hillslopes is often affected by distinct soil boundaries, including the frozen boundary and the organic‐mineral boundary (OMB), where highly porous and hydraulically conductive organic material overlies fine‐grained mineral soils. Viewed from the surface, ground cover appears as a patchwork on sub‐meter scales, with thick, moss mats interspersed with lichen‐covered, silty soils with gravel inclusions. We conducted a decameter‐scale subsurface tracer test on a boreal forest hillslope in interior Alaska to quantify locations and mechanisms of transport and storage in these soils, focusing on the OMB. A sodium bromide tracer was added as a slug addition to a pit and sampled at 40 down‐gradient wells, screened primarily at the OMB and within a 7 × 12 m well field. We maintained an elevated head in the injection pit for 8.5 hr to simulate a storm. Tracer breakthrough velocities ranged from <0.12 to 0.93 m hr^?1, with the highest velocities in lichen‐covered soils. After 12 hr and cessation of the elevated head, the tracer coalesced and was only detected in thick mosses at a trough in the OMB. By 24 hr, approximately 17% of the tracer mass could be accounted for. The majority of the mass loss occurred between 4 and 12 hr, while the tracer was in contact with lichen‐covered soils, which is consistent with tracer transport into deeper flow paths via preferential flow through discrete gravelly areas. Slow breakthroughs suggest that storage and exchange also occurred in shallow soils, likely related to saturation and drainage in fine‐grained mineral soils caused by the elevated hydraulic head. These findings highlight the complex nature of storage and transmission of water and solutes from boreal hillslopes to streams and are particularly relevant given rapid changes to boreal environments related to climate change, thawing permafrost and increasing fire severity.     

Contemporary rates of chemical denudation and atmospheric CO2 sequestration in glacier basins: an Arctic perspective

This paper presents new estimates of solute fluxes from five high Arctic glacier basins in Svalbard. These estimates are combined with data from two other glacier basins to assess the effectiveness of chemical denudation on Svalbard and to estimate rates of temporary (or transient) CO₂ drawdown. We use a solute provenance model to partition solutes into marine, aerosol, atmospheric and crustal components and to estimate their annual fluxes. Crustally derived solute fluxes are equivalent to a mean chemical denudation rate of 350 Σmeq⁺ m⁻² a⁻¹ for Svalbard (range: 160–560 Σmeq⁺ m⁻² a⁻¹), which lies within the global range of 94–4200 Σmeq⁺ m⁻² a⁻¹ for 21 glacier basins in the northern hemisphere, and is close to the continental average of 390 Σmeq⁺ m⁻² a⁻¹. Specific annual discharge is the most significant control upon chemical denudation in the glacierized basins, and basin lithology is an important secondary control, with carbonate‐rich and basaltic lithologies currently showing the greatest chemical denudation rates. Estimates of transient CO₂ drawdown are also directly associated with specific annual discharge and rock type. On Svalbard transient CO₂ drawdown lies in the range 110–3000 kg C km⁻² a⁻¹, whilst the range is 110–13000 kg C km⁻² a⁻¹ for the northern hemisphere glacial data set. Transient CO₂ drawdown is therefore usually low in the Svalbard basins unless carbonate or basalt rocks are abundant. The analysis shows that a large area of uncertainty in the transient CO₂ drawdown estimates exists due to the non‐stoichiometric release of solute during silicate hydrolysis. Silicate hydrolysis is particularly non‐stoichiometric in basins where the extent of glacierization is high, which is most probably an artefact of high flushing rates through ice‐marginal and subglacial environments where K‐feldspars are undergoing mechanical comminution. Copyright © 2000 John Wiley & Sons, Ltd.     

River flow regime and snow cover of the Pamir Alay (Central Asia) in a changing climate

Abstract

The Vakhsh and Pyandj rivers, main tributaries of the Amu Darya River in the mountainous region of the Pamir Alay, play an important role in the water resources of the Aral Sea basin (Central Asia). In this region, the glaciers and snow cover significantly influence the water cycle and flow regime, which could be strongly modified by climate change. The present study, part of a project funded by the European Commission, analyses the hydrological situation in six benchmark basins covering areas of between 1800 and 8400 km², essentially located in Tajikistan, with a variety of topographical situations, precipitation amounts and glacierized areas. Four types of parameter are discussed: temperature, glaciation, snow cover and river flows. The study is based mainly on a long-time series that ended in the 1990s (with the collapse of the Soviet Union) and on field observations and data collection. In addition, a short, more recent period (May 2000 to May 2002) was examined to better understand the role of snow cover, using scarce monitored data and satellite information. The results confirm the overall homogeneous trend of temperature increase in the mountain range and its impacts on the surface water regime. Concerning the snow cover, significant differences are noted in the location, elevation, orientation and morphology of snow cover in the respective basins. The changes in the river flow regime are regulated by the combination of the snow cover dynamics and the increasing trend of the air temperature.

Editor Z.W. Kundzewicz     

Synoptic sampling reveals spatial stability in flow and chemistry patterns despite large seasonal variability in a subarctic mountain catchment

Seasonality plays a critical role in cold mountain regions as variation in air temperature, ground thermal status, and precipitation phase alter the rate, timing and magnitude of hydrological and chemical transport. Additionally, cold mountain catchments can have highly variable topography, geology, permafrost, and landcover, which intrinsically add to this irregularity. Understanding how external and internal variability act to control mass fluxes requires sampling at a high spatial resolution over time, which rarely occurs in cold remote regions. In this work, we conduct five snapshot sampling surveys across 34 subcatchments during the ice-free period in Wolf Creek Research Basin (a mesoscale montane subarctic catchment) and two additional winter surveys across a subset of sites to assess the drivers of variability in stream chemistry and discharge. We sampled for specific conductance (SpC), major ions, and dissolved organic carbon (DOC) and used statistical metrics and Bayesian mixing analysis to quantify patterns of flow and chemistry across space and time. Our results indicate patterns in both flow and chemistry remain largely consistent across seasons for all solutes. However, there was weaker correlation of chemistry between sites, suggesting asynchronous behaviour within the catchment. There was evidence of increasing production of ions and DOC along the stream network during high spring flows but not during low flows. Although concentrations and flows exhibit high seasonality in subarctic mountains, this seasonal variability does not alter spatial patterns that arise from highly variable catchment characteristics.     

Characterization of contrasting flow and thermal regimes in two adjacent subarctic alpine headwaters in Northwest Canada

Alpine headwaters in subarctic regions are particularly sensitive to climate change, yet there is little information on stream thermal regimes in these areas and how they might respond to global warming. In this paper, we characterize and compare the hydrological and thermal regimes of two subarctic headwater alpine streams within an empirical framework. The streams investigated are located within two adjacent catchments with similar geology, size, elevation and landscape, Granger Creek (GC) and Buckbrush Creek (BB), which are part of the Wolf Creek Research Basin in the Yukon Territory, Canada. Hydrometeorological and high-resolution stream temperature data were collected throughout summer 2016. Both sites exhibited a flow regime typical of cold alpine headwater catchments influenced by frozen ground and permafrost. Comparatively, GC was characterized by a flashier response with more extreme flows, than BB. In both sites, stream temperature was highly variable and very responsive to short-term changes in climatic conditions. On average, stream temperature in BB was slightly higher than in GC (respectively 5.8 and 5.7°C), but less variable (average difference between 75th and 25th quantiles of 1.6 and 2.0°C). Regression analysis between mean daily air and stream temperature suggested that a greater relative (to stream flow) groundwater contribution in BB could more effectively buffer atmospheric fluctuations. Heat fluxes were derived and utilized to assess their relative contribution to the energy balance. Overall, non-advective fluxes followed a daily pattern highly correlated to short-wave radiation. G1enerally, solar radiation and latent heat were respectively the most important heat source and sink, while air–water interface processes were major factors driving nighttime stream temperature fluctuations.     

Bacteria‐sediment associations in an alpine,freshwater environment and their effects on particle size,density and settling velocity

This study measures the presence of bacteria‐sediment associations (BSAs) in an alpine, glacier‐fed watershed in the Southern Coast Mountains of British Columbia, Canada. The impact of BSAs on the creation of flocculated particles and their settling velocity are quantified using a laser transmissometer. Results from the study indicate that BSAs are present in the watershed and vary over both space and time. The percentage of bacteria associated with sediment particles was found to range from < 1% to 40%. Major sources of planktonic bacteria such as agricultural land and wastewater treatment outflow co‐occur with large decreases in the BSA ratio. Laboratory analysis demonstrates that an increase in the concentration of bacteria was associated with a decrease in the volume concentration of small particles, and a decrease in both estimated density and measured settling velocity for particles in larger size classes; consistent with flocculated particles of increasing complexity arising from combinations of primary particles and/or BSAs. Paleoenvironmental reconstructions using laminated lake sediments in alpine, glacier‐fed systems benefit from a fuller understanding of the geomorphologic processes by which they formed. While bacteria are noted to enhance formation of flocculated particles in laboratory systems, their impact upon geomorphic processes in natural systems have yet to be fully explored. © 2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Constraining water age dynamics in a south-eastern Australian catchment using an age-ranked storage and stable isotope approach

Improving our knowledge of the travel times of water through catchments is critical for the management and protection of water resources and to improve our understanding of fundamental catchment behaviour. In this study we use the age-ranked storage framework StorAge Selection (SAS) to investigate travel times in the Corin catchment, a headwater catchment in the south-east of Australia covered by native Eucalyptus species. Few studies have applied the SAS framework globally and in energy-intensive areas where catchment losses are heavily in favour of evapotranspiration relative to streamflow. A combination of observed and modelled values of oxygen-18 (δ¹⁸O ), the stable isotope in water, are used to constrain storage selection preferences of streamflow and evapotranspiration and the size of the catchment active storage. The highest performing parameter combinations that could reproduce δ¹⁸O in streamflow were dependent on a strong preference for young water in evapotranspiration, and a mixture of weak young and old water preference in streamflow. The mean travel time of streamflow over the study period 2007–2019, weighted by the flow rate, is limited to within a probable range of 2.81–9.77 years. The size of the active storage, a key parameter in the SAS framework, was poorly identified, and in combination with the isotopic inputs into the model, contributed to the uncertainty of the results. We discuss the implications of the results with respect to the study area, as well as within the context of SAS research globally and identify ways to improve the modelling process.     

Hydrogeological characterization of groundwater storage and drainage in an alpine karst aquifer (the Kanin massif,Julian Alps)

The Kanin massif is an important trans‐boundary aquifer, which stretches between Slovenia and Italy. The groundwater is only partially exploited, mainly for water supply, but the aquifer exhibits great potential for future exploitation. Since no consistent regional overview of the hydrogeological functioning of the Kanin massif was available, the decision was made to perform a study of this area, using a pragmatic approach based on 3D geological and hydrogeological modelling. The so‐called KARSYS approach was applied, with the aim of characterizing the groundwater reserves within this karst massif and of locating the main drainage axes that carry groundwater from the recharge areas to the respective springs. Delineation of the catchment areas of the corresponding springs was carried out, and some new explanations were obtained, especially with regard to the Mo?nica spring, which is located in Slovenia and forms a potential source of drinking water. It was found that this spring's catchment area extends as far as the Italian ski resort of Sella Nevea. The conceptual model also provides a possible explanation about the underground drainage towards the Boka spring and waterfall, which has been a challenge for decades. This new explanation is based on the existence of a perched groundwater body that feeds the Boka spring via a system of conduits. Despite some limitations, the results, which consist of a visualization of the underground drainage and groundwater storage within the Kanin massif, can be used as a basis for planning the sustainable management of karst waters in the studied area. Copyright © 2014 John Wiley & Sons, Ltd.