in silico
Phase 1 : In Silico Validation
- Objectif
- Determine whether hierarchical clustering on historical open meteorological data yields stable, physically interpretable regimes with statistically distinct ventilation coefficients, and whether a simulated control rule would have produced a measurable pollution reduction.
- Coût estimé
- €500-2000
- Durée estimée
- 4-6 weeks
- Critères de réussite
- Regime stability (centroid reproducibility across years) · Cosine similarity > 0.85 between year 1 and year 2 centroids for all regimes · (Compute cosine similarity between centroid vectors of corresponding regimes (matched by hierarchical clustering on pooled data))
- Ventilation coefficient separation between accumulation-prone and dispersion-prone regimes · Mean ventilation coefficient of accumulation-prone regime < 60% of dispersion-prone regime (as per Prediction 2) · (Bootstrap ratio test (10,000 resamples) on mean ventilation coefficients from ERA5 + station data)
- Simulated PM2.5 reduction under control rule · Statistically significant reduction (p < 0.05, one-tailed Welch’s t-test) of at least 10% in mean PM2.5 during accumulation-prone hours · (Two-sample t-test comparing actual vs simulated hourly PM2.5 during accumulation-prone regimes in year 3)
- On continue si
- All three success criteria are met (regime stability, ventilation separation, simulated reduction ≥10% with p<0.05)
- On arrête si
- Any criterion fails: e.g., silhouette score < 0.25, or simulated reduction <5%, or p > 0.10
- On réoriente si
- If regime stability is good but simulated reduction is 5-10% (p<0.10), consider increasing simulated scrubber efficiency to 60% or testing a different clustering algorithm (e.g., Gaussian Mixture Model) before proceeding
- Risques
- Open meteorological data does not have sufficient spatial coverage (e.g., only one station available)Probabilité : mediumParade : Use ERA5 grid point (0.25° resolution) as a pseudo-station; or use airport ASOS data which is widely available. If only one station, simulate a 10-station network by adding spatially correlated noise (spatial correlation decay length 500 m).
- Mixing layer height from ERA5 is too coarse or inaccurate for urban scaleProbabilité : mediumParade : Cross-validate with radiosonde data (IGRA) if available; or use a simpler proxy: ventilation coefficient = WS * (T_surface - T_850hPa) as a surrogate for stability.
- No clear elbow or silhouette peak for KProbabilité : lowParade : Use gap statistic or Davies-Bouldin index; if still ambiguous, fix K=4 based on physical expectations (stable, convective, advective, frontal).