The science
Depth-averaged steady groundwater flow: div(T grad h) + R - Q = 0, solved on a 49 x 49 grid by Gauss-Seidel with fixed heads on the outer boundary; T = K x 5 m, scaled locally by connected fissures; R is rainfall x recharge fraction (12 %).
Carbonate dissolution is a first-order rate law that grows with acidity, rainfall, saturation and flow, saturating as the void fraction approaches 97 %. Roof sag follows a smoothed 5 x 5 cavity-support proxy; once it exceeds 0.35 m a collapse event is counted.
The conduit drawn in the slice is a schematic at fixed depth whose height scales with the model void fraction; surface, head and void profiles are real model fields at the chosen z. Simplifications: 2.5-D aquifer, accelerated time, no real geochemistry. Educational, not a geotechnical risk assessment.