PyGEM-scripts — Frontal Ablation Calibration and Simulation branch (Zenodo Release v1.0.0)
Latest
🧊 PyGEM-scripts — Frontal Ablation Calibration and Simulation Branch (Zenodo Release v1.0.0)
This release corresponds to the calibration and simulation scripts (FA_AdaPBS_RT branch) used for the frontal ablation calibration and simulation experiments in the associated publication.
📄 Paper
Joint Bayesian Calibration of Frontal Ablation and Surface Mass Balance in Global Glacier Models
(GMD Preprint, 2026)
📦 Zenodo Archive
🔗 Related Components
This release is part of a three-component framework. The companion repositories and their tagged versions are:
| Component | Repository | Branch | Tag |
|---|---|---|---|
| OGGM core | Ruitangtang/oggm | SERMeQ_RT |
v1.0.0-zenodo-frontal-ablation |
| PyGEM | Ruitangtang/PyGEM | PyGEM_RT |
v1.0.0-zenodo-pygem-rt |
| Calibration scripts (this release) | Ruitangtang/PyGEM-scripts | FA_AdaPBS_RT |
v1.0.0-zenodo-fa-adapbs-rt |
🔧 Key Modifications
This branch introduces the following enhancements to the calibration and simulation workflow:
- Bayesian calibration framework (AdaPBS) — Integrated the Adaptive Particle Batch Smoother (Aalstad et al.,2026 ) for joint calibration for frontal ablation and climatic mass balance parameters estimations
- Frontal ablation observation integration — Added data loaders for annual terminus position change time series (Li et al., 2024)
- Parallel processing support — Enabled regional-scale calibration and simulation with Bash scripts for high-performance computing
- Custom output handlers — Extended outputs to include frontal ablation flux alongside standard mass balance variables
- Multi-GCM support — Configured for 12 CMIP6 models under SSP126 and SSP585 scenarios
📂 Contents
This repository provides the execution layer of the coupled modeling framework, with the following key components:
Calibration scripts:
run_calibration_AMIS_SERMeQ.py— Bayesian calibration of frontal ablation parameters (2000–2010) using the AdaPBS (Adaptive Particle Batch Smoother) framework
Simulation script:
run_simulation_AMIS_SERMeQ.py— Future projections (2000–2100) under CMIP6 scenarios (12 GCMs, SSP126/SSP585)
Parallel execution scripts (Bash):
run_calibration_AMIS_Region.sh— Regional calibration in parallelrun_simulation_AMIS_Region.sh— Regional simulations in parallelrun_simulation_AMIS_SingleGlac.sh— Single glacier simulation
Jupyter notebooks for paper figures:
Analysis_Hindcast_Sabinebreen_Fig1_2.ipynb— Hindcast analysis and model validation (Fig. 1, 2)Postprocess_calibration_20002020_Fig3.ipynb— Calibration results and parameter posteriors (Fig. 3)Postprocess_simulation_20002100_Fig4_5.ipynb— Future projections (Fig. 4, 5)
Supporting files:
- Configuration files for the Sabinebreen glacier (RGI60-07.00036)
- Utility functions for parallel processing and data handling
- Custom output handlers to export frontal ablation flux alongside standard mass balance outputs
🚀 Getting Started
This branch is the source-code reference for the modified calibration and simulation backend used by the frontal-ablation glacier workflow.
For a lightweight, review-friendly interface with YAML configs, dry-run quickstart commands, tests, method/results visuals, and reproducibility links, visit:
👉 Ruitangtang/frontal-ablation-glacier-demo
The full paper-related research package, including the single-glacier test data and execution materials, is archived at:
📦 Zenodo concept DOI 10.5281/zenodo.18761729
📌 Citation
If you use this code in your research, please cite both the paper and the Zenodo archive:
Joint Bayesian Calibration of Frontal Ablation and Surface Mass Balance in Global Glacier Models. GMD Preprint, 2026. DOI: 10.5194/egusphere-2026-1081
Frontal Ablation Module for OGGM-PyGEM (v1.0.0). Zenodo, 2026. DOI: 10.5281/zenodo.18761729