Skip to main content
Log in

Physics-based fine-scale numerical model of a karst system (Milandre Cave, Switzerland)

Modèle numérique à échelle détaillée et base physique d’un système karstique (Grotte de Milandre, Suisse)

Modelo numérico de base física en escala de detalle de un sistema cárstico (Milandre Cave, Suiza)

基于物理机制的岩溶系统 (瑞士米兰德尔洞穴)精细尺度数值模型

Modelo numérico de escala fina baseado em física de um sistema cárstico (Caverna Milandre, Suíça)

Hydrogeology Journal Aims and scope Submit manuscript

Abstract

In karst aquifers, groundwater flow is often turbulent and occurs in variably saturated conduits. This implies a nonlinear response to recharge events which cannot be reproduced with the tools commonly used for groundwater flow modeling. Recent studies have shown the usefulness of using conduit flow models to simulate flow in karst systems. However existing models are built on abstract representations of the actual conduit network or at a rather coarse resolution. Such models cannot provide information on local flow conditions in the conduits or be used to simulate mass transport. In the catchment of the Milandre Cave, in the Jura Mountains of Switzerland, a good portion of the active drainage network is accessible by caving, which opens the possibility for a distributed fine-scale numerical model. This report presents the development of a hydraulic model of the downstream part of this system (1.5 km × 0.5 km) using EPA SWMM 5. The network geometry is based on a detailed cave survey. Calibration is achieved by fitting the observed hydraulic head vs. flow-rate curves. The model performs well when compared against the hydraulic heads that were monitored throughout the system. Simulated transit times are also in line with tracer test results. This confirms that the model reproduces well the physics of flow in this karst system, while being spatially distributed at a fine scale (median element size of 6 m). It can thus provide information on local flow conditions in the conduits at various water stages and be used to address mass transport problems.

Résumé

Dans les aquifères karstiques, l’écoulement des eaux souterraines est généralement turbulent et prend place dans des conduits à saturation en eau variable. Ceci conduit à une réponse non linéaire aux événements de recharge, qui ne peut pas être reproduite avec les outils utilisés de manière classique pour la modélisation des écoulements d’eaux souterraines. De récentes études ont montré l’utilité d’avoir recours à des modèles d’écoulement en conduits pour simuler l’écoulement dans les systèmes karstiques. Cependant les modèles existants sont construits sur des représentations abstraites du réseau de conduits réels ou à une résolution des plus approximative. De tels modèles ne peuvent pas informer sur les conditions locales d’écoulement dans les conduits ou être utilisés pour simuler le transport de masse. Dans le bassin versant de la grotte de Milandre, dans les montagnes du Jura en Suisse, une bonne portion du réseau de drainage actif est accessible par les spéléologues, ce qui offre la possibilité d’établir un modèle numérique distribué à une échelle détaillée. Cet article présente le développement d’un modèle hydraulique de la partie aval du système (1.5 km × 0.5 km) en utilisant le code EPA SWMM 5. La géométrie du réseau est définie à partir d’une cartographie détaillée de la grotte. La calibration est réalisée en calant les charges hydrauliques observées par rapport aux hydrogrammes de débit. Le modèle reproduit correctement lorsqu’on compare les charges hydrauliques calculées avec celles observées tout au long du système. Les temps de transit simulés sont aussi cohérents avec les résultats des essais de traçage artificiels. Cela confirme que le modèle qui est distribué spatialement à une échelle fine (taille d’un élément médian de 6 m), reproduit bien la physique de l’écoulement dans ce système karstique. Il peut également fournir une information sur les conditions locales d’écoulement dans les conduits pour différents niveaux d’eau et peut être utilisés pour aborder les problèmes de transport de masse.

Resumen

En los acuíferos cársticos, el flujo de agua subterránea es con frecuencia turbulento y ocurre en conductos variablemente saturados. Esto implica una respuesta no lineal a los eventos de recarga que no puede ser reproducida con las herramientas comúnmente usadas para el modelado de flujo de agua subterránea. Estudios recientes han demostrado la utilidad de utilizar modelos de flujo de conducto para simular el flujo en sistemas cársticos. Sin embargo, los modelos existentes se basan en representaciones abstractas de la red real de conductos o en una resolución bastante gruesa. Estos modelos no pueden informar sobre las condiciones de flujo local en los conductos ni utilizarse para simular el transporte de masas. En la cuenca de Milandre Cave, en las Jura Mountains, en Suiza, una buena parte de la red de drenaje activa es accesible por la espeleología, lo que abre la posibilidad de un modelo numérico distribuido en escala de detalle. Este artículo presenta el desarrollo de un modelo hidráulico del sector aguas abajo de este sistema (1.5 km × 0.5 km) utilizando el SWMM 5 de la EPA. La geometría de la red se basa en un estudio detallado de las cuevas. La calibración se realiza ajustando la carga hidráulica observada frente a las curvas de caudal. El modelo funciona bien cuando se compara con las cargas hidráulicas que fueron monitoreados en todo el sistema. Los tiempos de tránsito simulados también están en línea con los resultados de los ensayos de monitoreo. Esto confirma que el modelo reproduce bien la física del flujo en este sistema cárstico, a la vez que se distribuye espacialmente a escala de detalle (tamaño del elemento medio de 6 m). Por lo tanto, puede informar sobre las condiciones de flujo local en los conductos en las distintas etapas del agua y utilizarse para abordar los problemas de transporte de masas.

摘要

岩溶含水层中地下水流通常是紊流,并且发生在变饱和管道中。这表明对补给的非线性响应不能用地下水流建模的传统工具再现。最近的研究表明使用管道流动模型来模拟岩溶系统中的流动是有效的。然而,现有模型基于实际管道网络的抽象概化或以较粗略的分辨率构建。这些模型不能了解管道中局部流动的条件或用于模拟质量传输。在瑞士汝拉山脉的米兰德尔洞穴区,很大一部分现行排水网络可以通道进入,这为分布式精细数值模型提供了可能性。本文介绍了使用EPA SWMM 5开发该系统下游部分(1.5 km × 0.5 km)的水力模型。网络几何尺寸是基于详细的洞穴调查。通过拟合观察到的水头与流量曲线进行了模型识别。与系统中监测的水头相比,该模型模拟效果良好。模拟的运输时间也与示踪剂测试结果一致。这证实了该模型很好地再现了该岩溶系统中的流动特性,同时在空间上以精细尺度分布(中等单元大小为6 m)。因此该模型可以分析不同水位管道中局部流动状况,并应用于质量输移问题。

Resumo

Em aquíferos cársticos, o fluxo das águas subterrâneas é frequentemente turbulento e ocorre em condutos saturados variáveis. Isto implica numa resposta não linear de eventos de recarga que não podem ser reproduzidos com ferramentas comumente usadas para modelagem de fluxo das águas subterrâneas. Estudos recentes mostraram a utilidade do uso de modelos de fluxo de conduto para simular o fluxo em sistemas cársticos. No entanto, os modelos existentes são construídos em representações abstratas de rede de conduto real ou em uma resolução bastante grosseira. Tais modelos não podem informar sobre condições de fluxo locais nos condutos ou ser usados para simular o transporte de massa. Na bacia da Caverna Milandre, nas Montanhas Jura da Suíça, uma boa parte da rede de drenagem ativa é acessível por espeleologia, o que abre a possibilidade de um modelo numérico distribuído de escala fina. Este artigo apresenta o desenvolvimento de um modelo hidráulico da parte a jusante deste sistema (1.5 km × 0.5 km) usando o EPA SWMM 5. A geometria da rede é baseada em um levantamento detalhado da caverna. A calibração é obtida ajustando a carga hidráulica observada versus as curvas de taxa de vazão. O modelo tem um bom desempenho quando comparado com as cargas hidráulicas monitoradas em todo o sistema. Os tempos de trânsito simulados também estão de acordo com resultados de ensaios com traçadores. Isto confirma que o modelo reproduz bem a física do fluxo neste sistema cárstico, ao mesmo tempo que está espacialmente distribuído a uma escala fina (tamanho médio do elemento de 6 m). Isso pode, portanto, informar sobre as condições de fluxo local nos condutos em vários estágios de água e ser usado para resolver problemas de transporte de massa.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8
Fig. 9
Fig. 10

References

  • Campbell C, Sullivan SM (2002) Simulating time-varying cave flow and water levels using the storm water management model. Eng Geol 65(2–3):133–139. https://doi.org/10.1016/S0013-7952(01)00120-X

    Article  Google Scholar 

  • Chen Z, Goldscheider N (2014) Modeling spatially and temporally varied hydraulic behavior of a folded karst system with dominant conduit drainage at catchment scale, Hochifen–Gottesacker, Alps. J Hydrol 514:41–52. https://doi.org/10.1016/j.jhydrol.2014.04.005

    Article  Google Scholar 

  • Chen Z, Hartmann A, Goldscheider N (2017) A new approach to evaluate spatiotemporal dynamics of controlling parameters in distributed environmental models. Environ Model Softw 87:1–16. https://doi.org/10.1016/j.envsoft.2016.10.005

    Article  Google Scholar 

  • Favre I (2001) Base de données des essais de traçage du plateau karstique de Bure (JU), SIG, interprétations statistiques [Database of Bure karst plateau (JU) tracer tests, GIS, statistical interpretations]. MSc Thesis, University of Neuchâtel, Switzerland

  • Swisstopo (Federal Office of Topography) (2016) Digital elevation model swissALTI3D, Swisstopo, Wabern bei Bern, Switzerland

  • Gabrovšek F, Peric B (2006) Monitoring the flood pulses in the epiphreatic zone of karst aquifers: the case of Reka River system, karst plateau, SW Slovenia. Acta Carsolog 35(1). https://doi.org/10.3986/ac.v35i1.241

  • Gabrovšek F, Peric B, Kaufmann G (2018) Hydraulics of epiphreatic flow of a karst aquifer. J Hydrol 560:56–74. https://doi.org/10.1016/j.jhydrol.2018.03.019

    Article  Google Scholar 

  • Gigon R, Wenger R (1986) Inventaire spéléologique de la Suisse, tome II: Canton du Jura [Speleological inventory of Switzerland, vol II: Canton of Jura]. Commission de Spéléologie de la Société helvétique des Sciences naturelles, La Chaux-de-Fonds, Switzerland

  • Gill LW, Naughton O, Johnston PM (2013) Modelling a network of turloughs in lowland karst. Water Resour Res 49(6):3487–3503

    Article  Google Scholar 

  • Grasso DA, Jeannin P-Y (1994) Estimation des pertes dans la partie aval du réseau karstique de la Milandrine: bilan hydrique au sein d’un aquifère karstique [Estimation of flow losses in the downstream part of the karst network of the Milandrine: water balance in a karst aquifer]. Bull Hydrogéol 13:115–128

    Google Scholar 

  • Hauns M, Jeannin PY, Atteia O (2001) Dispersion, retardation and scale effect in tracer breakthrough curves in karst conduits. J Hydrol 241(3–4):177–193. https://doi.org/10.1016/S0022-1694(00)00366-8

    Article  Google Scholar 

  • James W, Rossman L, James WRC (2010) User’s guide to SWMM5, 13th edn. Computational Hydraulics Int, Guelph, ON, Canada

  • Jeannin P-Y (2001) Modeling flow in phreatic and epiphreatic karst conduits in the Hölloch Cave (Muotatal, Switzerland). Water Resour Res 37(2):191–200. https://doi.org/10.1029/2000WR900257

    Article  Google Scholar 

  • Jeannin P-Y, Maréchal J-C (1995) Lois de pertes de charges dans les conduits karstiques: base théorique et observations [Head-loss laws in karst conduits: theoretical basis and observations]. Bull Hydrogéol 14:149–176

    Google Scholar 

  • Jeannin P-Y, Malard A, Rickerl D, Weber E (2015) Assessing karst-hydraulic hazards in tunneling: the Brunnmühle spring system—Bernese Jura, Switzerland. Environ Earth Sci 74(12):7655–7670. https://doi.org/10.1007/s12665-015-4655-5

    Article  Google Scholar 

  • Kaufmann G, Gabrovšek F, Turk J (2016) Modelling cave flow hydraulics in Postojnska Jama, Slovenia. Acta Carsolog 45:57–70. https://doi.org/10.3986/ac.v45i1.3059

    Article  Google Scholar 

  • Kiraly L (1995) Effect of the epikarst on the hydrograph of karst springs: a numerical approach. Bull Hydrogéol 14:199–220

    Google Scholar 

  • Kovács A, Jeannin P-Y (2003) Hydrogeological overview of the Bure plateau, Ajoie, Switzerland. Eclogae Geol Helv 96(3):367–380. https://doi.org/10.1007/s00015-003-1099-9

    Article  Google Scholar 

  • Kresic N (2010) Modeling. In: Kresic N, Stevanovic Z (eds) Groundwater hydrology of springs, vol 2010. Butterworth-Heinemann, Burlington, MA, pp 165–230

  • Maloszewski P, Stichler W, Zuber A, Rank D (2002) Identifying the flow systems in a karstic-fissured-porous aquifer, the Schneealpe, Austria, by modelling of environmental 18O and 3H isotopes. J Hydrol 256(1):48–59. https://doi.org/10.1016/S0022-1694(01)00526-1

    Article  Google Scholar 

  • Maréchal J-C (1994) Étude et modélisation de l’hydraulique et du transport dans les drains karstiques [Study and modelling of hydraulics and transport in karst drains]. MSc Thesis, Université de Neuchâtel, Switzerland

  • Perrin J (2003) A conceptual model of flow and transport in a karst aquifer based on spatial and temporal variations of natural tracers. PhD Thesis, Université de Neuchâtel, Switzerland

  • Peterson EW, Wicks CM (2006) Assessing the importance of conduit geometry and physical parameters in karst systems using the storm water management model (SWMM). J Hydrol 329(1–2):294–305. https://doi.org/10.1016/j.jhydrol.2006.02.017

    Article  Google Scholar 

  • White WB (1969) Conceptual models for carbonate aquifers. Ground Water 7(3):15–21. https://doi.org/10.1111/j.1745-6584.1969.tb01279.x

  • Worthington SRH (2003) A comprehensive strategy for understanding flow in carbonate aquifer. Speleogenesis Evol Karst Aquif 1(1):1–8

    Google Scholar 

  • Worthington SRH (2009) Diagnostic hydrogeologic characteristics of a karst aquifer (Kentucky, USA). Hydrogeol J 17(7):1665–1678. https://doi.org/10.1007/s10040-009-0489-0

    Article  Google Scholar 

  • Wu Y, Jiang Y, Yuan D, Li L (2008) Modeling hydrological responses of karst spring to storm events: example of the Shuifang spring (Jinfo Mt., Chongqing, China). Environ Geol 55(7):1545–1553. https://doi.org/10.1007/s00254-007-1105-z

    Article  Google Scholar 

Download references

Acknowledgements

Some of the measurements analyzed in this study were supported by the A16 highway project, the Canton of Jura and the Swiss Federal Roads Office (FEDRO) within a contract with MFR Géologie-Géotechnique SA in Delémont. Special thanks is dedicated to Marc Hessenauer who took care of the data and shared them with us. The Spéléo-Club Jura and its president Pierre Xavier Meury have also been supportive of our research and made the access to the cave possible under good conditions. We thank R. de Rooij, F. Gabrovšek and two anonymous reviewers for their helpful comments.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Cécile Vuilleumier.

Electronic supplementary material

ESM 1

(PDF 3005 kb)

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Vuilleumier, C., Jeannin, PY. & Perrochet, P. Physics-based fine-scale numerical model of a karst system (Milandre Cave, Switzerland). Hydrogeol J 27, 2347–2363 (2019). https://doi.org/10.1007/s10040-019-02006-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s10040-019-02006-y

Keywords

Navigation