Pith. sign in

REVIEW 3 major objections 5 minor 90 references

Glacier data assimilation on an Arctic glacier: Learning from large ensemble twin experiments

T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

Pith's one-line read Joint assimilation of synthetic albedo and snow depth observations into a full energy-balance glacier model improves annual surface mass balance skill by up to 86% in CRPS terms, with either of two ensemble smoothers.

desk verdict Solid twin-experiment evaluation of joint albedo and snow depth assimilation for glacier SMB, but the abstract overclaims universal improvement when the authors' own results show degradation in the accumulation area under low snowfall. read the letter →

arxiv 2502.09314 v1 pith:OBVJDSLB submitted 2025-02-13 physics.geo-ph physics.data-an

classification physics.geo-phphysics.data-an
keywords dataassimilationglaciermassbalanceparticlebatchsmootherensemblealbedosnowdepthtwinexperimentsCryoGrid
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper tries to establish that feeding noisy, satellite-like observations of albedo and snow depth into a full energy-balance glacier model, through either of two ensemble data-assimilation schemes, can substantially reduce the uncertainty of simulated surface mass balance on an Arctic glacier. In identical-twin experiments on Kongsvegen, Svalbard, jointly assimilating both observation types lowered the continuous ranked probability score (CRPS) of annual mass balance by up to 86% relative to the prior ensemble in the ablation zone, with average joint improvements of roughly 63% across all glacier zones. The authors argue that this makes particle batch smoothing and ensemble smoothing practical tools, requiring no differentiable model, for constraining the two parameters that drive much of CryoGrid's mass-balance uncertainty: the snow albedo decay rate and the snowfall bias factor. Joint assimilation is presented as the most robust choice across the tested climatic scenarios, since it automatically balances the complementary information in the two observation types.

What carries the argument

The load-bearing mechanism is a two-parameter Bayesian inversion carried out with ensemble smoothers: the parameter vector $\theta=(\tau_a,\beta_s)$ pairs the snow albedo decay rate $\tau_a$ with the multiplicative snowfall bias factor $\beta_s$, both treated as constant within a mass-balance year. A generalized logit-normal prior on each parameter is combined with a Gaussian observation-error likelihood, and the two schemes differ in how they approximate the posterior: the particle batch smoother resamples prior particles by their likelihood weights, while the ensemble smoother performs an ensemble Kalman update in logit-transformed parameter space and reruns the model with the updated parameters. The data-generating map $G$ is the CryoGrid glacier configuration with the CROCUS albedo scheme and full surface energy balance, and the continuous ranked probability score against synthetic truth is the evaluation metric.

What would settle it

Run the same assimilation on Kongsvegen using real MODIS albedo and ICESat-2 snow depth retrievals instead of synthetic ones, and compare the posterior annual mass balance against the glacier's independent stake measurements; a posterior CRPS or RMSE that fails to beat the prior by a margin comparable to the twin-experiment gains would show that the perfect-model assumption is broken.

Watch

Extended reading notes

Core claim

The central claim is that, when the model is assumed perfect and all uncertainty is confined to two time-invariant parameters, jointly assimilating synthetic albedo and snow depth observations into CryoGrid improves annual surface mass balance skill by up to 86% in CRPS terms, and that this improvement holds across glacier zones with either the particle batch smoother or the ensemble smoother. The particle batch smoother is better at representing albedo dynamics, while the ensemble smoother is better for snow depth under low snowfall conditions where the prior ensemble sometimes fails to bracket the truth. The paper also claims that joint assimilation yields the best or tied-best results in the majority of experiments, and that CRPS gains level off once the ensemble size reaches about 100 members, so the method is computationally achievable for a model as expensive as CryoGrid.

Load-bearing premise

The whole result rests on the identical-twin assumption that CryoGrid is a perfect model of the glacier and that all uncertainty lives in two fixed parameters plus Gaussian observation noise; if the real system behaves differently, the reported 86% improvement may not appear with real data.

Editorial extensions

If this is right

  • Joint assimilation of albedo and snow depth is the robust default: it gave the best or tied-best CRPS in 10 of 12 particle-batch-smoother experiments and the best overall performance for both schemes.
  • Albedo assimilation alone is the cheaper option when only one variable is available, with the particle batch smoother requiring just one CryoGrid run per ensemble member; snow depth assimilation may require the ensemble smoother and its second set of model runs to avoid overconfident posteriors.
  • A prior ensemble of about 100 members captures most of the possible error reduction; enlarging it to 1000 gives diminishing returns in mean CRPS and in Monte Carlo resampling variance.
  • Because the workflow only needs an ensemble of parameter draws and likelihood evaluations, it is transferable to other glaciers and other energy-balance models without requiring a differentiable model.
  • The twin experiments suggest that the choice of observation type should depend on snowfall regime: albedo assimilation helps more under low snowfall, while snow depth assimilation helps more under high snowfall.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the twin-experiment gains survive contact with real data, the same framework could turn MODIS and ICESat-2 retrievals into a routine mass-balance constraint for unmonitored Arctic glaciers, where in situ stake networks are sparse; the paper does not test this.
  • The complementary strengths of the two smoothers suggest an untested hybrid: use an ensemble Kalman update when snow depth is informative and the prior does not bracket the truth, and particle weighting when albedo dominates, switching by snowfall regime.
  • Relaxing the strong-constraint assumption to include forcing bias or initial snow-state uncertainty would make the method robust to model structural error; the twin setup deliberately avoids this, so the real-world performance ceiling is unknown.
  • The bootstrap analysis implies a practical design rule: sample a large prior ensemble once, then use resampling to choose a cost-effective ensemble size before committing to expensive full model runs.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper presents a set of identical-twin (Observing System Simulation) experiments in which synthetic albedo and snow depth observations, generated by the CryoGrid glacier mass balance model, are assimilated into the same model using two ensemble-based batch smoothers: the Particle Batch Smoother (PBS) and the Ensemble Smoother (ES). The experiments cover three grid cells representing ablation, equilibrium-line, and accumulation zones of Kongsvegen glacier, under four combinations of true albedo evolution rate and snowfall factor. The authors report large reductions in the Continuous Ranked Probability Score (CRPS) of posterior surface mass balance relative to the prior, including an up-to-86% improvement for joint assimilation in the ablation area, and they compare the two schemes in terms of accuracy, precision, and computational cost. A sensitivity analysis of the PBS to ensemble size is performed by bootstrapping from a 1000-member ensemble.

Significance. If the reported results hold, the paper provides a useful proof-of-concept that ensemble-based batch smoothers can constrain glacier surface mass balance parameters from albedo and snow depth observations, with potential relevance to future assimilation of MODIS and ICESat-2 products. The use of a full energy-balance model (CryoGrid), a 12-year forcing period, a 1000-member ensemble, and a comparison of two contrasting smoothers are strengths. The paper also includes machine-checkable algorithmic descriptions and a sensitivity analysis based on explicit bootstrap resampling. However, the central claim as stated in the abstract and comparison section is broader than what the paper's own figures and text support, and the identical-twin design limits the transferability of the quantitative improvements to real-world applications. These issues are fixable by qualification, but they currently affect the headline message.

major comments (3)
  1. [Abstract; Results (PBS section, Fig. 5); Conclusions] The abstract states that joint assimilation 'achieves improved mass balance simulations across different glacier zones using either assimilation scheme,' and the Comparison section states that 'the posterior always improved over the prior in terms of mass balance CRPS.' These claims are contradicted by the paper's own results: the PBS Results section explicitly says that in the accumulation area, snow depth assimilation under low snowfall scenarios increases CRPS relative to the prior, and that joint assimilation behaves similarly to snow depth assimilation alone in that area. The Conclusions also carve out 'all experiments except those given by low snowfall level in the accumulation area.' The central claim must be qualified to exclude the low-snowfall accumulation-area case, or the abstract, Results, Comparison, and Conclusions must be made consistent with Fig. 5.
  2. [Table 2; Comparison of two data assimilation schemes] Table 2 reports average CRPS improvements over four scenarios without per-scenario breakdowns, uncertainties, or confidence intervals. Given that at least one scenario (low snowfall in the accumulation area) shows degradation relative to the prior for snow depth and joint assimilation under the PBS, the reported averages (e.g., 25.5% for PBS ACC joint assimilation) can hide systematically negative cases. The authors should report per-scenario values or provide error bars/confidence intervals around the averages, and avoid claiming universal improvement based on averaged numbers.
  3. [Data assimilation (identical twin, strong constraint); Conclusions] The paper acknowledges the identical-twin and perfect-model assumptions in the Data assimilation section, but the Conclusions go on to state that the approach is 'potentially transferable for estimating mass balance of all glaciers on Svalbard.' This extrapolation is not supported by the experimental design, since the synthetic truth is generated with the same model whose parameters are being inferred, and structural model error is excluded by construction. The authors should either soften the transferability claim or add a prominent caveat that real-world structural errors may change both the magnitude of improvements and the PBS-versus-ES ranking.
minor comments (5)
  1. [Comparison of two data assimilation schemes] The text refers to 'Fig. ??' when discussing posterior annual mass balance results in the ELA region; the correct figure reference (presumably Fig. 6) should be inserted.
  2. [Data assimilation, Eq. (6)] The notation 'by' for the predicted observations is unconventional and could be confused with the observation vector y; consider using ŷ or G(θ) throughout for clarity.
  3. [Table 1] The prior scale σ0 is set to 1 for both τa and βs, but no justification or sensitivity analysis is provided for this choice; a brief comment on how this scale was selected would help.
  4. [Sensitivity of data assimilation performance to the ensemble size] The description of the bootstrap procedure is clear, but the sentence 'the resampling variation remains non-zero even when the ensemble size is 1000' could be sharpened by noting that bootstrapping with replacement never exactly reproduces the original pool; the current wording is slightly redundant.
  5. [Results, PBS section] The statement that 'joint assimilation tends to yield the best (including ties) results across the majority of experiments (10 out of 12)' is useful, but it would be even more informative to state explicitly how many of those 10 cases are in the ablation/ELA regions versus the accumulation area, given the known exception.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the posterior SMB improvements are obtained from noisy synthetic observations through standard DA, not from the target SMB itself.

full rationale

The paper's central claim is that joint assimilation of albedo and snow depth improves CryoGrid SMB skill relative to the prior in twin experiments. The synthetic truth is generated with the same CryoGrid model used for assimilation, so the evaluation is internal rather than external; the authors state this explicitly ('we avoid model structural uncertainty by confining ourselves to so-called identical twin experiments'). However, this does not reduce the result to its inputs by construction: the noisy observations are of albedo and snow depth, the target quantity is annual SMB, and the posterior SMB emerges from the PBS/ES update, not from a direct fit of the reported CRPS improvements. The improvement percentages are diagnostics of the posterior, not fitted parameters renamed as predictions. Citations to prior work by the authors (Aalstad et al. 2018; Alonso-González et al. 2022; Schmidt et al. 2023) support methodology and model configuration but are not invoked as external proof of the present result. The internal inconsistency between the abstract's unqualified 'across different glacier zones' and the paper's own Fig. 5/Conclusions exception for low-snowfall accumulation-area cases is a claims-consistency concern, not a circularity. The identical-twin limitation is real for real-world transferability but belongs under validity/correctness risk rather than circular reasoning.

Assumptions & free parameters 8 free parameters · 7 assumptions · 0 invented entities

The central claim rests on a long list of chosen prior hyperparameters, observation error scales, and the identical-twin assumption that CryoGrid is the truth. None of these are derived in the paper; they are imported from prior literature or the authors' earlier work. No new physical entities are introduced. The experiment's quantitative conclusions are therefore conditional on these inputs.

free parameters (8)
  • Prior median for tau_a (albedo evolution rate) = 0.005 day^-1
    Sets the center of the logit-normal prior in Table 1; chosen from prior modeling literature, not from Kongsvegen-specific calibration.
  • Prior median for beta_s (snowfall factor) = 1
    Assumes no snowfall bias correction by default; center of the prior in Table 1.
  • Prior scale sigma0 for tau_a and beta_s = 1
    Controls the spread of the prior ensemble; a wider or narrower prior would change the likelihood weights and the reported CRPS improvements.
  • Prior bounds for tau_a = [0.0001, 0.05] day^-1
    Double-bounded support of the prior; if the true value is outside these bounds, PBS cannot represent it, which is implicated in the accumulation-area collapse under low snowfall.
  • Prior bounds for beta_s = [0.5, 2]
    Support of the prior for the snowfall factor; determines whether the synthetic truth is bracketed by the prior ensemble.
  • Albedo observation error standard deviation = 0.1
    Set from the upper limit of MODIS albedo RMSE (Stroeve and others, 2005); directly determines the albedo likelihood sharpness.
  • Snow depth observation error standard deviation = 0.5 m
    Set from ICESat-2 snow depth retrieval errors on low slopes (Deschamps-Berger and others, 2023); determines the snow depth likelihood.
  • True scenario parameter values (four theta* combinations) = Not stated numerically
    The four scenarios are described qualitatively (rapid/slow albedo evolution, high/low snowfall), but the exact numeric truth values are not reported; these define the synthetic truth and therefore underlie every reported CRPS improvement and method ranking.
assumptions (7)
  • domain assumption CryoGrid is a perfect model of the glacier system (identical-twin assumption).
    Invoked in the Data assimilation section: the authors 'avoid model structural uncertainty by confining ourselves to so-called identical twin experiments' and 'assumes a perfect data generating model'. This is the main premise that makes the CRPS improvements interpretable as an upper bound.
  • domain assumption All model uncertainty is captured by the two time-invariant parameters tau_a and beta_s.
    The parameter vector is Np=2; forcing, initial state, and structural errors are treated as negligible. Set out in the Prior and likelihood and Twin experiments sections.
  • domain assumption CARRA reanalysis forcing is unbiased except for a multiplicative snowfall correction.
    The model applies a relative bias correction beta_s; no other meteorological bias is considered. Described in the Mass balance model section.
  • domain assumption Observation errors are independent, zero-mean Gaussian with known variances.
    Used to define the likelihood in Eq. (6) and the diagonal covariance matrix R in the Prior and likelihood section; a standard but unverified assumption for the synthetic observations.
  • domain assumption The generalized logit-normal prior with the hyperparameters in Table 1 is appropriate.
    Prior defined in Eqs. (3)-(5); hyperparameters are chosen 'building on several related studies' (citations) rather than direct Kongsvegen measurements.
  • standard math CRPS is a valid skill score for evaluating probabilistic ensemble predictions.
    Standard definition and properties of CRPS from Gneiting and others (2005) and Hersbach (2000); used throughout for evaluation.
  • domain assumption Three CARRA grid cells are representative of Kongsvegen's glacier zones and other Svalbard glaciers.
    The experiments are run for one cell per zone (ablation, ELA, accumulation); the conclusion that the method is 'potentially transferable for estimating mass balance of all glaciers on Svalbard' extrapolates from these three cells.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Glacier data assimilation on an Arctic glacier: Learning from large ensemble twin experiments." pith.science (2026). https://pith.science/paper/OBVJDSLB

@misc{pith2026250209314,
  author       = {Pith},
  title        = {Pith review of: Glacier data assimilation on an Arctic glacier: Learning from large ensemble twin experiments},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/OBVJDSLB}},
  note         = {Machine review of arXiv:2502.09314}
}
read the original abstract

Glacier modeling is crucial for quantifying the evolution of cryospheric processes. At the same time, uncertainties hamper process understanding and predictive accuracy. Here, we suggest improving glacier mass balance simulations for the Kongsvegen glacier in Svalbard through the application of Bayesian data assimilation techniques in a set of large ensemble twin experiments. Noisy synthetic observations of albedo and snow depth, generated using the multilayer CryoGrid community model with a full energy balance, are assimilated using two ensemble-based data assimilation schemes: the particle batch smoother and the ensemble smoother. A comprehensive evaluation exercise demonstrates that the joint assimilation of albedo and snow depth improves the simulation skill by up to 86% relative to the prior in specific glacier regions. The particle batch smoother excels in representing albedo dynamics, while the ensemble smoother is particularly effective for snow depth under low snowfall conditions. By combining the strengths of both observations, the joint assimilation achieves improved mass balance simulations across different glacier zones using either assimilation scheme. This work underscores the potential of ensemble-based data assimilation methods for refining glacier models by offering a robust framework to enhance predictive accuracy and reduce uncertainties in cryospheric simulations. Further advances in glacier data assimilation will be critical to better understanding the fate and role of Arctic glaciers in a changing climate.

Figures

Figures reproduced from arXiv: 2502.09314 by the authors.

Figure 1
Figure 1. Atmospherically corrected shortwave infrared false color image over the area surrounding Kongsvegen glacier near Ny-Ålesund in the Svalbard archipelago captured by the Sentinel-2B satellite at 13:07 UTC on the 25th of August 2020. The image shows the locations of Ny-Ålesund (yellow star) and the Kongsvegen glacier outline from RGI (white) as well as the locations of 2.5 by 2.5 km grid cells that were extracted from … view at source ↗
Figure 2
Figure 2. Workflow in the twin experiments involving the sequential generation of: synthetic truth runs (orange), noisy synthetic observa￾tions (green), and data assimilation experiments (blue) followed by the evaluation of each experiment (purple) [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Structure of the large ensemble twin experiments based on permutations of four parameter scenarios, three types of assimilated observation vectors, and three experimental areas generating a total of 36 twin experiments. The scenarios combine either a rapid or slow albedo evolution rate with either a high or low snowfall factor. The assimilated observation vectors are either albedo only, snow depth only, or joint ass… view at source ↗
Figures from the paper (5 more)
Figure 4
Figure 4. Figure 4: Comparison of prior, posterior, and true surface mass balance in the ablation area when using the PBS to assimilate albedo only, snow depth only, and both observations jointly. The figure presents four scenarios based on the snow albedo evolution rates and snowfall fac…
Figure 5
Figure 5. Figure 5: Continuous ranked probability score (CRPS) for the prior and posterior mass balance after assimilating albedo, snow depth, and both observations jointly using the PBS under all scenarios, compared to synthetic true mass balance [PITH_FULL_IMAGE:figures/full_fig_p014_5.png]
Figure 6
Figure 6. Figure 6: Comparison of the performance of two assimilation schemes applied to the ELA area in terms of RMSE (top row) and ensemble standard deviation (bottom row) for the Particle Batch Smoother (left panels a and c) and the Ensemble Smoother (right panels b and d) [PITH_FULL_…
Figure 7
Figure 7. Figure 7: Sensitivity of the posterior surface mass balance CRPS to ensemble size following joint assimilation of albedo and snow depth using the PBS scheme under the scenario of R&H in the ablation area. For each ensemble size (Ne ) the CRPS statistics were estimated by resampl…
Figure 8
Figure 8. Figure 8: Comparison of the prior, posterior, and true annual surface mass balance in the accumulation area when assimilating different types of observations with the PBS (a) and the ES (b) under a low snowfall factor and slow albedo evolution rate with low snowfall where the po…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

90 extracted references · 60 canonical work pages

  1. [1]

    Cryosphere, 12(1), 247--270, ISSN 19940424 ( 10.5194/tc-12-247-2018 )

    Aalstad K, Westermann S, Schuler TV, Boike J and Bertino L (2018) Ensemble-based assimilation of fractional snow-covered area satellite retrievals to estimate the snow distribution at Arctic sites . Cryosphere, 12(1), 247--270, ISSN 19940424 ( 10.5194/tc-12-247-2018 )

  2. [2]

    Alonso-Gonz \' a lez E, Aalstad K, Baba MW, Revuelto J, L \' o pez-Moreno JI, Fiddes J, Essery R and Gascoin S (2022) The Multiple Snow Data Assimilation System (MuSA v1.0) . Geosci. Model Dev., 15(24), 9127--9155, ISSN 19919603 ( 10.5194/GMD-15-9127-2022 )

  3. [3]

    Arnold CP and Dey CH (1986) Observing-systems simulation experiments: past, present and future. Bull. - Am. Meteorol. Soc., 67(6), 687--695, ISSN 0003-0007 ( 10.1175/1520-0477(1986)067<0687:OSSEPP>2.0.CO;2 )

  4. [4]

    Bengtsson L, Andrae U, Aspelien T, Batrak Y, Calvo J, Rooy Wd, Gleeson E, Hansen-Sass B, Homleid M, Hortal M, Ivarsson KI, Lenderink G, Niemel \" a S, Nielsen KP, Onvlee J, Rontu L, Samuelsson P, Mu \ n oz DS, Subias A, Tijm S, Toll V, Yang X and Køltzow M (2017) The HARMONIE–AROME Model Configuration in the ALADIN–HIRLAM NWP System . Mon. Weather. Rev., ...

  5. [5]

    Bertino L, Evensen G and Wackernagel H (2003) Sequential data assimilation techniques in oceanography. Int. Stat. Rev., 71(2), 223--241 ( https://doi.org/10.1111/j.1751-5823.2003.tb00194.x )

  6. [6]

    Authorea Prepr

    Bertoncini A, Pomeroy J and Pomeroy JW (2024) Assimilation of Satellite Albedo to Improve Simulations of Glacier Hydrology . Authorea Prepr. ( 10.22541/AU.171128344.49725407/V1 )

  7. [7]

    Brinkerhoff DJ (2022) Variational inference at glacier scale . J. Comput. Phys., 459, 111095, ISSN 0021-9991 ( 10.1016/J.JCP.2022.111095 )

  8. [8]

    Tellus, 21(5), 611--619 ( https://doi.org/10.1111/j.2153-3490.1969.tb00466.x )

    Budyko MI (1969) The effect of solar radiation variations on the climate of the Earth . Tellus, 21(5), 611--619 ( https://doi.org/10.1111/j.2153-3490.1969.tb00466.x )

Show all 90 references
  1. [9]

    WIREs Clim

    Carrassi A, Bocquet M, Bertino L and Evensen G (2018) Data assimilation in the geosciences: An overview of methods, issues, and perspectives . WIREs Clim. Chang., 9, e535 ( 10.1002/wcc.535 )

  2. [10]

    ISSN 01672789 ( 10.1016/j.physd.2024.134436 )

    Choi B and Lee Y (2024) Sampling error mitigation through spectrum smoothing: first experiments with ensemble transform Kalman filters and Lorenz models . ISSN 01672789 ( 10.1016/j.physd.2024.134436 )

  3. [11]

    Cryosphere, 17(12), 5499--5517, ISSN 19940424 ( 10.5194/TC-17-5499-2023 )

    Choi Y, Seroussi H, Morlighem M, Schlegel NJ and Gardner A (2023) Impact of time-dependent data assimilation on ice flow model initialization and projections: a case study of Kjer Glacier, Greenland . Cryosphere, 17(12), 5499--5517, ISSN 19940424 ( 10.5194/TC-17-5499-2023 )

  4. [12]

    Springer ( 10.1007/978-3-030-47845-2 )

    Chopin N and Papaspiliopoulos O (2020) An Introduction to Sequential Monte Carlo . Springer ( 10.1007/978-3-030-47845-2 )

  5. [13]

    Copernicus Climate Change Service (C3S) Climate Data Store (CDS) (2024) Arctic regional reanal- ysis on single levels from 1991 to present, Copernicus Climate Change Service (C3S) Climate Data Store (CDS) ( 10.24381/cds.713858f6 )

  6. [14]

    Bootstrap Method Their Appl

    Davison AC and Hinkley DV (1997) Bootstrap Methods and their Application . Bootstrap Method Their Appl. ( 10.1017/CBO9780511802843 )

  7. [15]

    Cryosphere, 17(7), 2779--2792, ISSN 19940424 ( 10.5194/TC-17-2779-2023 )

    Deschamps-Berger C, Gascoin S, Shean D, Besso H, Guiot A and L \' o pez-Moreno JI (2023) Evaluation of snow depth retrievals from ICESat-2 using airborne laser-scanning data . Cryosphere, 17(7), 2779--2792, ISSN 19940424 ( 10.5194/TC-17-2779-2023 )

  8. [16]

    Dumont M, Durand Y, Arnaud Y and Six D (2012) Variational assimilation of albedo in a snowpack model and reconstruction of the spatial mass-balance distribution of an alpine glacier . J. Glaciol., 58(207), 151--164, ISSN 0022-1430 ( 10.3189/2012JOG11J163 )

  9. [17]

    Evensen G, Vossepoel FC and van Leeuwen PJ (2022) Data Assimilation Fundamentals ( 10.1007/978-3-030-96709-3 )

  10. [18]

    Ezat MM, Fahl K and Rasmussen TL (2024) Arctic freshwater outflow suppressed Nordic Seas overturning and oceanic heat transport during the Last Interglacial . Nat. Commun. 2024 15:1, 15(1), 1--9, ISSN 2041-1723 ( 10.1038/s41467-024-53401-3 )

  11. [19]

    Forbes R, Tompkins AM and Untch A (2011) A new prognostic bulk microphysics scheme for the IFS . Tech. Mem., (649), 22 ( 10.21957/bf6vjvxk )

  12. [20]

    Gabarr \' o C, Hughes N, Wilkinson J, Bertino L, Bracher A, Diehl T, Dierking W, Gonzalez-Gambau V, Lavergne T, Madurell T, Malnes E and Wagner PM (2023) Improving satellite-based monitoring of the polar regions: Identification of research and capacity gaps . Front. In Remote....

  13. [21]

    Nature, 601(7893), 374--379, ISSN 0028-0836 ( 10.1038/s41586-021-04314-4 )

    Geyman EC, J J van Pelt W, Maloof AC, Aas HF and Kohler J (2022) Historical glacier change on Svalbard predicts doubling of mass loss by 2100 . Nature, 601(7893), 374--379, ISSN 0028-0836 ( 10.1038/s41586-021-04314-4 )

  14. [22]

    Cryosphere, 14(3), 811--832, ISSN 19940424 ( 10.5194/TC-14-811-2020 )

    Gillet-Chaulet F (2020) Assimilation of surface observations in a transient marine ice sheet model using an ensemble Kalman filter . Cryosphere, 14(3), 811--832, ISSN 19940424 ( 10.5194/TC-14-811-2020 )

  15. [23]

    Gneiting T, Raftery AE, Westveld AH and Goldman T (2005) Calibrated probabilistic forecasting using ensemble model output statistics and minimum crps estimation. Mon. Weather. Rev., 133(5), 1098 -- 1118 ( 10.1175/MWR2904.1 )

  16. [24]

    Cryosphere, 17(8), 3505--3533, ISSN 19940424 ( 10.5194/TC-17-3505-2023 )

    Groenke B, Langer M, Nitzbon J, Westermann S, Gallego G and Boike J (2023) Investigating the thermal state of permafrost with Bayesian inverse modeling of heat transfer . Cryosphere, 17(8), 3505--3533, ISSN 19940424 ( 10.5194/TC-17-3505-2023 )

  17. [25]

    Guidicelli M, Aalstad K, Treichler D and Salzmann N (2024) A combined data assimilation and deep learning approach for continuous spatio-temporal swe reconstruction from sparse ground tracks. J. Hydrol. X, 25, 100190 ( https://doi.org/10.1016/j.hydroa.2024.100190 )

  18. [26]

    Cryosphere, 17(9), 3955--3986, ISSN 19940424 ( 10.5194/TC-17-3955-2023 )

    Gunnarsson A, Gardarsson SM and P \' a lsson F (2023) Modeling of surface energy balance for Icelandic glaciers using remote-sensing albedo . Cryosphere, 17(9), 3955--3986, ISSN 19940424 ( 10.5194/TC-17-3955-2023 )

  19. [27]

    Hagen JO, Melvold K, Eiken T, Isaksson E and Lefauconnier B (1999) Mass balance methods on Kongsvegen, Svalbard . Geogr. Ann. A: Phys. Geogr., 81(4), 593--601, ISSN 04353676 ( 10.1111/j.0435-3676.1999.00087.x )

  20. [28]

    Hersbach H (2000) Decomposition of the continuous ranked probability score for ensemble prediction systems. Weather. Forecast., 15(5), 559 -- 570 ( 10.1175/1520-0434(2000)015<0559:DOTCRP>2.0.CO;2 )

  21. [29]

    Hock R (2003) Temperature index melt modelling in mountain areas . J. Hydrol., 282(1-4), 104--115, ISSN 0022-1694 ( 10.1016/S0022-1694(03)00257-9 )

  22. [30]

    Hock R and Holmgren B (2005) A distributed surface energy-balance model for complex topography and its application to Storglaci \" a ren, Sweden . J. Glaciol., 51(172), 25--36, ISSN 0022-1430 ( 10.3189/172756505781829566 )

  23. [31]

    Hopwood MJ, Carroll D, Dunse T, Hodson A, Holding JM, Iriarte JL, Ribeiro S, Achterberg EP, Cantoni C, Carlson DF, Chierici M, Clarke JS, Cozzi S, Fransson A, Juul-Pedersen T, Winding MH and Meire L (2020) Review article: How does glacier discharge affect marine biogeochemistr...

  24. [32]

    Huss M and Hock R (2015) A new model for global glacier change and sea-level rise . Front. In Earth Sci., 3(September), 1--22, ISSN 22966463 ( 10.3389/feart.2015.00054 )

  25. [33]

    Immerzeel WW, Lutz AF, Andrade M, Bahl A, Biemans H, Bolch T, Hyde S, Brumby S, Davies BJ, Elmore AC, Emmer A, Feng M, Fern \' a ndez A, Haritashya U, Kargel JS, Koppes M, Kraaijenbrink PD, Kulkarni AV, Mayewski PA, Nepal S, Pacheco P, Painter TH, Pellicciotti F, Rajaram H, Ru...

  26. [34]

    IPCC (2022) IPCC Special Report on the Ocean and Cryosphere in a Changing Climate — ( https://doi.org/10.1017/9781009157964 )

  27. [35]

    Isaksen K, Nordli , F rland EJ, upikasza E, Eastwood S and Nied \' z wied \' z T (2016) Recent warming on Spitsbergen—Influence of atmospheric circulation and sea ice cover . J. Geophys. Res. Atmospheres, 121(20), 913--11, ISSN 2169-8996 ( 10.1002/2016JD025606 )

  28. [36]

    Cambridge University Press ( 10.1017/CBO9780511790423 )

    Jaynes E (2003) Probability Theory . Cambridge University Press ( 10.1017/CBO9780511790423 )

  29. [37]

    Karner F, Obleitner F, Krismer T, Kohler J and Greuell W (2013) A decade of energy and mass balance investigations on the glacier Kongsvegen, Svalbard . J. Geophys. Res. Atmospheres, 118(10), 3986--4000, ISSN 21698996 ( 10.1029/2012JD018342 )

  30. [38]

    ESS Open Arch

    Keetz L, Aalstad K, Fisher R, Poppe Teran C, Naz B, Pirk N, Yilmaz Y and Skarpaas O (2024) Inferring parameters in a complex land surface model by combining data assimilation and machine learning . ESS Open Arch. ( 10.22541/essoar.172070530.05098424/v1 )

  31. [39]

    Atmospheric Meas

    Kotarba AZ (2022) Impact of the revisit frequency on cloud climatology for CALIPSO, EarthCARE, Aeolus, and ICESat-2 satellite lidar missions . Atmospheric Meas. Tech., 15(14), 4307--4322, ISSN 18678548 ( 10.5194/AMT-15-4307-2022 )

  32. [40]

    Cryosphere, 15(11), 5017--5040, ISSN 19940424 ( 10.5194/tc-15-5017-2021 )

    Landmann JM, K \" u nsch HR, Huss M, Ogier C, Kalisch M and Farinotti D (2021) Assimilating near-real-time mass balance stake readings into a model ensemble using a particle filter . Cryosphere, 15(11), 5017--5040, ISSN 19940424 ( 10.5194/tc-15-5017-2021 )

  33. [41]

    Cryosphere, 14(7), ISSN 19940424 ( 10.5194/tc-14-2253-2020 )

    Lenaerts JT, Drew Camron M, Wyburn-Powell CR and Kay JE (2020) Present-day and future Greenland Ice Sheet precipitation frequency from CloudSat observations and the Community Earth System Model . Cryosphere, 14(7), ISSN 19940424 ( 10.5194/tc-14-2253-2020 )

  34. [42]

    Lenaerts JTM, Medley B, van den Broeke MR and Wouters B (2019) Observing and Modeling Ice Sheet Surface Mass Balance . Rev. Geophys., 57(2), 376--420, ISSN 8755-1209 ( 10.1029/2018RG000622 )

  35. [43]

    Lind S, Ingvaldsen RB and Furevik T (2018) Arctic warming hotspot in the northern Barents Sea linked to declining sea-ice import . Nat. Clim. Chang., 8(7), 634--639, ISSN 17586798 ( 10.1038/s41558-018-0205-y )

  36. [44]

    Cambridge University Press

    MacKay DJC (2003) Information Theory, Inference, and Learning Algorithms . Cambridge University Press

  37. [45]

    Water Resour

    Magnusson J, Winstral A, Stordal AS, Essery R and Jonas T (2017) Improving physically based snow simulations by assimilating snow depths using the particle filter . Water Resour. Res., 53(2), 1125--1143, ISSN 1944-7973 ( 10.1002/2016WR019092 )

  38. [46]

    Malone AG, Doughty AM and MacAyeal DR (2019) Interannual climate variability helps define the mean state of glaciers . J. Glaciol., 65(251), 508--517, ISSN 0022-1430 ( 10.1017/JOG.2019.28 )

  39. [47]

    Margulis SA, Girotto M, Cort\' e s G and Durand M (2015) A particle batch smoother approach to snow water equivalent estimation . J. Hydrometeorol., 16, 1752–1772 ( 10.1175/JHM-D-14-0177.1 )

  40. [48]

    Markus T, Neumann T, Martino A, Abdalati W, Brunt K, Csatho B, Farrell S, Fricker H, Gardner A, Harding D, Jasinski M, Kwok R, Magruder L, Lubin D, Luthcke S, Morison J, Nelson R, Neuenschwander A, Palm S, Popescu S, Shum CK, Schutz BE, Smith B, Yang Y and Zwally J (2017) The ...

  41. [49]

    Marshall G, Rees W and Dowdeswell J (1993) Limitations imposed by cloud cover on multitemporal visible band satellite data sets from polar regions . Ann. Glaciol., 17, 113--120, ISSN 0260-3055 ( 10.3189/S0260305500012696 )

  42. [50]

    Cryosphere, 6 ( 10.5194/tc-6-1295-2012 )

    Marzeion B, Jarosch AH and Hofer M (2012) The Cryosphere Past and future sea-level change from the surface mass balance of glaciers . Cryosphere, 6 ( 10.5194/tc-6-1295-2012 )

  43. [51]

    Earth's Future, 8(7), e2019EF001470, ISSN 2328-4277 ( 10.1029/2019EF001470 )

    Marzeion B, Hock R, Anderson B, Bliss A, Champollion N, Fujita K, Huss M, Immerzeel WW, Kraaijenbrink P, Malles JH, Maussion F, Radi \' c V, Rounce DR, Sakai A, Shannon S, van de Wal R and Zekollari H (2020) Partitioning the Uncertainty of Ensemble Projections of Global Glacie...

  44. [52]

    Springer ( 10.1007/978-3-540-74703-1_24 )

    Masutani M, Schlatter TW, Errico RM, Stoffelen A, Andersson E, Lahoz W, Woollen JS, Emmitt GD, Riish jgaard LP and Lord SJ (2010) Observing System Simulation Experiments, 647--679. Springer ( 10.1007/978-3-540-74703-1_24 )

  45. [53]

    Mazzolini M, Aalstad K, Alonso-Gonz \' a lez E, Westermann S and Treichler D (2024) Spatio-temporal snow data assimilation with the ICESat-2 laser altimeter ( 10.5194/EGUSPHERE-2024-1404 )

  46. [54]

    Dyn., 11(4), 885--901, ISSN 21904987 ( 10.5194/ESD-11-885-2020 )

    Milinski S, Maher N and Olonscheck D (2020) How large does a large ensemble need to be? Earth Syst. Dyn., 11(4), 885--901, ISSN 21904987 ( 10.5194/ESD-11-885-2020 )

  47. [55]

    Navari M, Margulis SA, Tedesco M, Fettweis X and van de Wal RS (2021) Reanalysis Surface Mass Balance of the Greenland Ice Sheet Along K-Transect (2000–2014) . Geophys. Res. Lett., 48(17), e2021GL094602, ISSN 1944-8007 ( 10.1029/2021GL094602 )

  48. [56]

    Neuenschwander AL and Magruder LA (2019) Canopy and Terrain Height Retrievals with ICESat-2: A First Look . Remote. Sens. 2019 Vol. 11 Page 1721, 11(14), 1721, ISSN 2072-4292 ( 10.3390/RS11141721 )

  49. [57]

    No \" e l B, Jakobs CL, van Pelt WJ, Lhermitte S, Wouters B, Kohler J, Hagen JO, Luks B, Reijmer CH, van de Berg WJ and van den Broeke MR (2020) Low elevation of Svalbard glaciers drives high mass loss variability . Nat. Commun. 2020 11:1, 11(1), 1--8, ISSN 2041-1723 ( 10.1038...

  50. [58]

    Polar Res., 33(1 SUPPL), 1--23, ISSN 17518369 ( 10.3402/polar.v33.21349 )

    Nordli , Przybylak R, Ogilvie AE and Isaksen K (2014) Long-term temperature trends and variability on spitsbergen: The extended svalbard airport temperature series, 1898-2012 . Polar Res., 33(1 SUPPL), 1--23, ISSN 17518369 ( 10.3402/polar.v33.21349 )

  51. [59]

    stby TI, Schuler TV and Westermann S (2014) Severe cloud contamination of MODIS Land Surface Temperatures over an Arctic ice cap, Svalbard . Remote. Sens. Environ., 142, 95--102, ISSN 0034-4257 ( 10.1016/J.RSE.2013.11.005 )

  52. [60]

    Cryosphere, 11(1), 191--215, ISSN 19940424 ( 10.5194/tc-11-191-2017 )

    stby TI, Vikhamar Schuler T, Ove Hagen J, Hock R, Kohler J and Reijmer CH (2017) Diagnosing the decline in climatic mass balance of glaciers in Svalbard over 1957-2014 . Cryosphere, 11(1), 191--215, ISSN 19940424 ( 10.5194/tc-11-191-2017 )

  53. [61]

    Atmospheric Meas

    Pirk N, Aalstad K, Westermann S, Vatne A, van Hove A, Tallaksen LM, Cassiani M and Katul G (2022) Inferring surface energy fluxes using drone data assimilation in large eddy simulations. Atmospheric Meas. Tech., 15(24), 7293--7314 ( 10.5194/amt-15-7293-2022 )

  54. [62]

    Pontes GM and Menviel L (2024) Weakening of the Atlantic Meridional Overturning Circulation driven by subarctic freshening since the mid-twentieth century . Nat. Geosci. 2024 17:12, 17(12), 1291--1298, ISSN 1752-0908 ( 10.1038/s41561-024-01568-1 )

  55. [63]

    Polar Res., 38, ISSN 1751-8369 ( 10.33265/POLAR.V38.3364 )

    Pramanik A, Kohler J, Schuler TV, van Pelt W and Cohen L (2019) Comparison of snow accumulation events on two High-Arctic glaciers to model-derived and observed precipitation . Polar Res., 38, ISSN 1751-8369 ( 10.33265/POLAR.V38.3364 )

  56. [64]

    Rantanen M, Karpechko AY, Lipponen A, Nordling K, Hyv \" a rinen O, Ruosteenoja K, Vihma T and Laaksonen A (2022) The Arctic has warmed nearly four times faster than the globe since 1979 . Commun. Earth & Environ. 2022 3:1 , 3(1), 1--10, ISSN 2662-4435 ( 10.1038/s43247-022-00498-3 )

  57. [65]

    Cryosphere, 17(7), 2705--2724 ( 10.5194/tc-17-2705-2023 )

    Raoult N, Charbit S, Dumas C, Maignan F, Ottl\'e C and Bastrikov V (2023) Improving modelled albedo over the greenland ice sheet through parameter optimisation and modis snow albedo retrievals. Cryosphere, 17(7), 2705--2724 ( 10.5194/tc-17-2705-2023 )

  58. [66]

    Robinson G, Grooms I and Kleiber W (2018) Improving Particle Filter Performance by Smoothing Observations . Mon. Weather. Rev., 146(8), 2433--2446, ISSN 1520-0493 ( 10.1175/MWR-D-17-0349.1 )

  59. [67]

    Rounce DR, Khurana T, Short MB, Hock R, Shean DE and Brinkerhoff DJ (2020) Quantifying parameter uncertainty in a large-scale glacier evolution model using Bayesian inference: application to High Mountain Asia . J. Glaciol., 66(256), 175--187, ISSN 0022-1430 ( 10.1017/JOG.2019.91 )

  60. [68]

    Science, 379(6627), 78--83, ISSN 10959203 ( 10.1126/SCIENCE.ABO1324/SUPPL \_ FILE/SCIENCE.ABO1324 \_ SM.PDF )

    Rounce DR, Hock R, Maussion F, Hugonnet R, Kochtitzky W, Huss M, Berthier E, Brinkerhoff D, Compagno L, Copland L, Farinotti D, Menounos B and McNabb RW (2023) Global glacier change in the 21st century: Every increase in temperature matters . Science, 379(6627), 78--83, ISSN 1...

  61. [69]

    Part I: Theory

    Sacher W and Bartello P (2008) Sampling Errors in Ensemble Kalman Filtering. Part I: Theory . Mon. Weather. Rev., 136(8), 3035--3049, ISSN 1520-0493 ( 10.1175/2007MWR2323.1 )

  62. [70]

    Sandven S, Spreen G, Heygster G, Girard-Ardhuin F, Farrell SL, Dierking W and Allard RA (2023) Sea Ice Remote Sensing—Recent Developments in Methods and Climate Data Sets . Surv. In Geophys., 44(5), 1653--1689, ISSN 15730956 ( 10.1007/S10712-023-09781-0/METRICS )

  63. [71]

    London Math

    Sanz-Alonso D, Stuart A and Taeb A (2023) Inverse Problems and Data Assimilation. London Math. Soc. Student Texts, Cambridge University Press ( 10.1017/9781009414319 )

  64. [72]

    Cambridge University Press, 2 edition ( 10.1017/9781108917407 )

    S \"a rkk \"a S and Svensson L (2023) Bayesian Filtering and Smoothing . Cambridge University Press, 2 edition ( 10.1017/9781108917407 )

  65. [73]

    o rnsson H (2017) The importance of accurate glacier albedo for estimates of surface mass balance on Vatnaj \

    Schmidt LS, Aoalgeirsdottir G, Gumundsson S, Langen PL, P \' a lsson F, Mottram R, Gascoin S and Bj \" o rnsson H (2017) The importance of accurate glacier albedo for estimates of surface mass balance on Vatnaj \" o kull: Evaluating the surface energy budget in a regional clim...

  66. [74]

    EGUsphere, 2023, 1--32

    Schmidt LS, Schuler TV, Thomas EE and Westermann S (2023) Meltwater runoff and glacier mass balance in the high Arctic: 1991-2022 simulations for Svalbard . EGUsphere, 2023, 1--32

  67. [75]

    Sjursen KH, Dunse T, Tambue A, Schuler TV and Andreassen LM (2023) Bayesian parameter estimation in glacier mass-balance modelling using observations with distinct temporal resolutions and uncertainties . J. Glaciol., 1--20, ISSN 0022-1430 ( 10.1017/JOG.2023.62 )

  68. [76]

    Stroeve J, Box JE, Gao F, Liang S, Nolin A and Schaaf C (2005) Accuracy assessment of the MODIS 16-day albedo product for snow: Comparisons with Greenland in situ measurements . Remote. Sens. Environ., 94(1), 46--60, ISSN 00344257 ( 10.1016/j.rse.2004.09.001 )

  69. [77]

    van Leeuwen P (2009) Particle Filtering in Geophysical Systems . Mon. Weather. Rev., 137, 4089--4114 ( 10.1175/2009MWR2835.1. )

  70. [78]

    van Leeuwen PJ (2020) A consistent interpretation of the stochastic version of the ensemble kalman filter. Q. J. R. Meteorol. Soc., 146(731), 2815--2825 ( https://doi.org/10.1002/qj.3819 )

  71. [79]

    van Leeuwen PJ and Evensen G (1996) Data Assimilation and Inverse Methods in Terms of a Probabilistic Formulation . Mon. Weather. Rev., 124(12), 2898 -- 2913 ( 10.1175/1520-0493(1996)124<2898:DAAIMI>2.0.CO;2 )

  72. [80]

    Cryosphere, 13(9), 2259--2280, ISSN 19940424 ( 10.5194/tc-13-2259-2019 )

    Van Pelt W, Pohjola V, Pettersson R, Marchenko S, Kohler J, Luks B, Ove Hagen J, Schuler TV, Dunse T, No \" e l B and Reijmer C (2019) A long-term dataset of climatic mass balance, snow conditions, and runoff in Svalbard (1957-2018) . Cryosphere, 13(9), 2259--2280, ISSN 199404...

  73. [81]

    Vionnet V, Brun E, Morin S, Boone A, Faroux S, Le Moigne P, Martin E and Willemet JM (2012) The detailed snowpack scheme Crocus and its implementation in SURFEX v7.2 . Geosci. Model Dev., 5(3), 773--791, ISSN 1991959X ( 10.5194/GMD-5-773-2012 )

  74. [82]

    Wang S, Lu X, Cheng X, Li X, Peichl M and Mammarella I (2018) Limitations and Challenges of MODIS-Derived Phenological Metrics Across Different Landscapes in Pan-Arctic Regions . Remote. Sens. 2018 Vol. 10 Page 1784, 10(11), 1784, ISSN 2072-4292 ( 10.3390/RS10111784 )

  75. [83]

    Wei M, Qiao F, Guo Y, Deng J, Song Z, Shu Q and Yang X (2019) Quantifying the importance of interannual, interdecadal and multidecadal climate natural variabilities in the modulation of global warming rates . Clim. Dyn., 53(11), 6715--6727, ISSN 14320894 ( 10.1007/S00382-019-0...

  76. [84]

    u ller B, Filhol S, K \

    Westermann S, Ingeman-Nielsen T, Scheer J, Aalstad K, Aga J, Chaudhary N, Etzelm \" u ller B, Filhol S, K \" a \" a b A, Renette C, Schmidt LS, Schuler TV, Zweigel RB, Martin L, Morard S, Ben-Asher M, Angelopoulos M, Boike J, Groenke B, Miesner F, Nitzbon J, Overduin P, Stuenz...

  77. [85]

    Yang C, Leonelli FE, Marullo S, Artale V, Beggs H, Nardelli BB, Chin TM, De Toma V, Good S, Huang B, Merchant CJ, Sakurai T, Santoleri R, Vazquez-Cuervo J, Zhang HM and Pisano A (2021) Sea Surface Temperature Intercomparison in the Framework of the Copernicus Climate Change Se...

  78. [86]

    Ye F, Cheng Q, Hao W, Hu A and Liang D (2024) Unveiling Glacier Mass Balance: Albedo Aggregation Insights for Austrian and Norwegian Glaciers . Remote. Sens. 2024 Vol. 16 Page 1914, 16(11), 1914, ISSN 2072-4292 ( 10.3390/RS16111914 )

  79. [87]

    Nature, 568(7752), 382--386, ISSN 14764687 ( 10.1038/S41586-019-1071-0 )

    Zemp M, Huss M, Thibert E, Eckert N, McNabb R, Huber J, Barandun M, Machguth H, Nussbaumer SU, G \" a rtner-Roer I, Thomson L, Paul F, Maussion F, Kutuzov S and Cogley JG (2019) Global glacier mass changes and their contributions to sea-level rise from 1961 to 2016 . Nature, 5...

  80. [88]

    Zhang Q, Shen Z, Pokhrel Y, Farinotti D, Singh VP, Xu CY, Wu W and Wang G (2023) Oceanic climate changes threaten the sustainability of Asia’s water tower . Nat. 2023 615:7950, 615(7950), 87--93, ISSN 1476-4687 ( 10.1038/s41586-022-05643-8 )

  81. [89]

    , " * write output.state after.block = add.period write newline

    ENTRY address author booktitle chapter doi edition editor eid howpublished institution isbn issn journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentenc...

  82. [90]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 gl...

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.