Résumé
The Barrois region corresponds to a Tithonian-age limestone plateau located to the east of the Paris Basin, and constituting the most karstified area in the basin. This limestone is primarily exposed, except to the north of the plateau, where it is found under a semi-permeable Lower Cretaceous-age discordant cover where karstification is currently active. Karstification likely began during a phase of emersion of the summit of the carbonate reservoir (Jaillet, 2000) and has continued to the present day in various forms (ghost-rock karstification, hypogenic alteration, crypto-alteration, and epigenic alteration under cover) depending on the regional base level drop and the exposure of the Tithonian limestones through erosion of the Cretaceous cover.In the past, an attempt to model this network was made in the context of surface pollution issues (Jaquet, 2004). Although this modeling was applied across the entire plateau, it was limited by a 2D modeling of the phreatic network along different horizons compartmentalizing the Barrois limestones, a purely epigenic karstification mechanism, and coarse drain resolution.With the goal of better understanding the geometry of the karst network in the Barrois region to enhance the protection of the associated water resource, we propose to use the KarstNSim method (Gouy et al., 2024) to stochastically simulate various plausible 3D network skeletons within the Barrois limestones. This method allows constraining the network's morphology based on simple input parameters but also enables the replication of polygenic karsts through sequential simulations where each iteration influences subsequent ones.