The artificial lighting of caves adapted for touristic visits, leads to the appearance and propagation of a complex community of phototrophic organisms known as "lampenflora". Formed mainly by algae and cyanobacteria, they produce the degradation of the colonized substrates and decrease the show value of the caves. This phenomenon became famous worldwide in the 1960s due to the damage caused to the paintings in the Lascaux Cave (France). Since then it has become an issue of serious concern to both managers of show caves and to the international scientific community. Over time, the problem has been approached following two complementary strategies: preventing colonization by the invading organisms or eliminating them once they have become established through the use of chemical products, mainly biocides and strong oxidants. This kind of treatment generates pollutant effluents that can move the problem from the walls of the caves to the groundwater. This paper presents a critical literature review of the problem and the proposed solutions, and emphasizes the need for further study of the optimal doses of treatment chemicals and to develop quantitative methods to determine their effectiveness. . The lampenflora in show caves and its treatment: an emerging ecological problem.

Simulatingflow through multilayer aquifer systems is essential for groundwater management applications such evaluating natural recharge, assessing managed aquifer recharge (MAR), the characterization of surface water–groundwater interactions, and groundwater pollution risk analysis. However, dealing with flow through variably saturated porous media with perched water tables is a difficult task. Usual codes for groundwater flow modeling such as MODFLOW, TOUGH2, and FEFLOW require tedious procedures based on simplified assumptions and approximations. This paper presents a simple approach to simulate two‐dimensional flow through perched aquifers and aquitards in deep vadose zones by combining the Dirichlet and Neumann boundary conditions that may apply to any code simulating unsaturated flow. The proposed approach is illustrated with the unsaturated flow code VS2DTI. The main strengths of the proposed methodology are (a) variably saturated flow in the multilayer system is solved using Richards’ equation, taking into account the lateral flow that sustains observed water levels in perched aquifers and variations in recharge; (b) a Dirichlet‐type boundary condition at the top perched water table is substituted by a Neumann‐type boundary condition, allowing for the representation of any disturbance (e.g., infiltration from ponds, pumping from wells, etc.), regardless of duration and intensity; and (c) the impact of the disturbance is evaluated by comparing the responses of the undisturbed and the disturbed systems. The versatility of this methodology is applied to a MAR case study of a deep aquifer in a sedimentary basin where aquitards limit its feasibility.