in silico
Phase 1: In Silico Validation
- Objective
- Test whether MFS-informed hydrostratigraphic models explain significantly more variance in recharge blockage and flow compartmentalization than topographic models, using existing public datasets, before any field campaign.
- Estimated cost
- €500-2000
- Estimated duration
- 4-8 weeks
- Success criteria
- ΔR² (MFS-informed vs topographic model) · ΔR² ≥ 0.15 in at least 2 of 3-5 basins · (Cross-validated R² on held-out wells; Bayesian posterior P(ΔR²≥0.20) > 0.5)
- Sensitivity of K contrast on simulated residence time · Sobol first-order index S1 ≥ 0.20 · (SALib Sobol analysis on 10^4 Monte Carlo runs)
- MFS continuity required for compartmentalization · ≥1 km lateral continuity yields ≥30% head difference · (MODFLOW scenario runs with variable MFS patch size)
- Go if
- ΔR² ≥ 0.15 in ≥2 basins AND Sobol S1(K contrast) ≥ 0.20 AND at least 2 candidate sites identified with adequate well control
- No-go if
- ΔR² < 0.05 in all basins OR MFS layer has negligible effect on heads/residence times (S1 < 0.05)
- Pivot if
- ΔR² 0.05-0.15 in one basin: pivot to testing TS/LFS surfaces or to a single-basin focused study with higher-resolution stratigraphic data
- Risks
- Public well logs lack sufficient stratigraphic resolution to confidently place MFSProbability: highMitigation: Use biostratigraphic/isotope markers from published literature; fall back to TS/LFS as proxy surfaces; restrict to basins with published sequence-stratigraphic frameworks
- MODFLOW cannot resolve thin (1-5 m) MFS aquitard at regional scaleProbability: mediumMitigation: Use local grid refinement (MODFLOW 6 DISV) or FEFLOW unstructured mesh; test equivalent K via harmonic averaging
- Confounding by anthropogenic pumping in public head dataProbability: highMitigation: Filter wells by distance from pumping centers; use pre-1970 head data where available; include pumping as covariate
- Publication bias in A/S ratio compilationsProbability: mediumMitigation: Use only studies with quantitative accommodation curves; triangulate with Wheeler diagrams