scieee AI-readable full text Open interactive document viewer

A New Method to Statistically Combine Ice Sheet Mass Balance Estimates

Ozerov, Anthony; Pérez, Fernando; Hills, Benjamin

Abstract

We explore how to statistically combine multiple datasets estimating a common physical quantity. Three main methods are used to estimate the mass balance of the polar ice sheets: gravimetry, altimetry, and the input-output method. These yield different results, which leads to the question of how to combine them into one “best” mass balance estimate. Such measurement aggregation is a form of “meta-analysis”, which is a mature area of statistics. Using real-world data, we find that prominent meta-analysis methods, originally developed within medicine and social sciences, should not always be applied to problems in the physical sciences; that is, because their assumptions are too strong, they can yield intervals which are too narrow and understate our uncertainty. We propose new methods for meta-analysis with weaker assumptions that are more likely to be met in reality. The Ice sheet Mass Balance Inter-comparison Exercise (IMBIE) statistically aggregates mass balance estimates through a simple error-weighted mean and provides uncertainties which consider neither the heterogeneity between estimates nor the correlation in errors across time and between ice sheets. To sidestep the issue of error correlation, we estimate all final quantities (for example, mass trend in a given window for a given ice sheet) individually for the different mass balance time series datasets. We then perform a meta-analysis on the level of the final estimate using our proposed methods, and show illustrative results. Finally, we ask whether we should try to statistically reconcile different results in the first place, and on what basis we might trust the results.

Full text

Anthony Ozerov1 [email protected] Benjamin H. Hills2 Fernando Pérez1 1University of California, Berkeley (Statistics) 2Colorado School of Mines (Geophysics) A New Method to Statistically Combine Ice Sheet Mass Balance Estimates CURRENT METHOD TO COMBINE ESTIMATES: 1. Combine time series into one 2. Estimate trend [These are all of the publicly-available Antarctic mass-balance time series we could locate] ↓Assumes errors are independent across time and space. ↓Does not account for studies underestimating uncertainty. ↑Does not consider independence between different methods. This method also requires each underlying study to report uncertainties, and assumes these are calculated correctly and all mean the same thing. PROPOSED STATISTICAL METHOD FOR COMBINING TRENDS We can combine estimates without relying on their uncertainty quantification. Cumulative mass balance is the direct object of measurement for all methods except those using input-output / flux accounting. We can directly estimate the trend from the cumulative mass balance, or fit a linear model and estimate the slope. Maybe just making such a figure is enough and we can stop here. Really need a confidence interval! [e.g. for IPCC] Null hypothesis of the Sign Test. Confidence interval on median via inverting the test. Null hypothesis of the Signed-Rank Test. Confidence interval on median via inverting the test. Possible sets of assumptions for the combination ( θ : truth, : trend from study i.) PROPOSED METHOD TO COMBINE ESTIMATES: 1. Separately estimate trends 2. Combine trends ●No need to consider the different time resolutions ●Without even considering the errors attached to each method, we can visually see how different the results are. Monthly difference series (for each study) Cumulative mass balance series Cumulative mass balance series Cumulative mass balance series Trend Combine within each method (Grav., Alt., IOM) Combine between methods Simple average Addition Some parts of the calculation make uncertainties ↑larger than they should be, while others make them ↓smaller than they should be. The result is something that sometimes looks reasonable. Conversion from cumulative series to difference series and back reveals arbitrary ↑widening of uncertainties. Trend uncertainty × time does not coincide with the uncertainty on the cumulative series. Monthly difference series (for each method) Monthly difference series This is the method used by the Ice Sheet Mass Balance Intercomparison Exercise (IMBIE) team. [Otosaka, Inès & IMBIE Team, 2023] Antarctic Ice Sheet Mass Balance AGU 2025, #C13D-0819 Monday, 15 December, 14:15 - 17:45 Within each method (Grav., Alt., IOM), the results are clearly not independent. We can pick a “worst-case” set of three results: Altimetry B, Gravimetry B, and Input-Output. Then, we can construct a 75% confidence interval using one of the above sets of assumptions. ●When considering only three estimates, both sets of assumptions give the same result. ●Instead of a very conservative worst-case analysis, we could also take the mean or median within each of the groups. 75% confidence interval Scan for online copy of poster