Preprints
https://doi.org/10.5194/egusphere-2025-3556
https://doi.org/10.5194/egusphere-2025-3556
08 Oct 2025
 | 08 Oct 2025
Status: this preprint is open for discussion and under review for The Cryosphere (TC).

Data-driven equation discovery of a sea ice albedo parametrisation

Diajeng W. Atmojo, Katja Weigel, Arthur Grundner, Marika M. Holland, Dmitry Sidorenko, and Veronika Eyring

Abstract. In the Finite-Element Sea Ice Model (FESIM), a part of the Finite-Element Sea ice Ocean Model (FESOM), sea ice albedo is treated as a tuning parameter defined by four constant values depending on snow cover and surface temperature. This parametrisation is too simple to capture the spatiotemporal variability of observed sea ice albedo. Here, we aim for an improved parametrisation by discovering an interpretable, physically consistent equation for sea ice albedo using symbolic regression, an interpretable machine learning technique, combined with physical constraints. Leveraging daily pan-Arctic satellite and reanalyses data from 2013 to 2020, we apply sequential feature selection which identifies snow depth, surface temperature, sea ice thickness and 2 m air temperature as the most informative features for sea ice albedo. As a function of these features, our data-driven equation identifies two critical mechanisms for determining sea ice albedo: the high sensitivity of sea ice albedo to small changes in thin snow and a weighted difference of the sea ice surface and 2 m air temperature, serving as a seasonal proxy that indicates the transition between melting and freezing conditions. To understand how additional model complexity reduces errors, we evaluate our discovered equation against baseline models with different complexities, such as multilayer perceptron neural networks (NNs) and polynomials on an error-complexity plane, showing that the equation excels in balancing error and complexity and reduces the mean squared error by about 51 % compared to the current FESIM parametrisation. Unlike NNs, our discovered equation allows for further regional and seasonal analyses due to its inherent interpretability. By fine-tuning its coefficients we uncover differences in physical conditions that drive sea ice albedo. This study demonstrates that learning an equation from observational data can deepen the process-level understanding of the Arctic Ocean's surface radiative budget and improve climate projections.

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.
Share
Diajeng W. Atmojo, Katja Weigel, Arthur Grundner, Marika M. Holland, Dmitry Sidorenko, and Veronika Eyring

Status: open (until 19 Nov 2025)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
Diajeng W. Atmojo, Katja Weigel, Arthur Grundner, Marika M. Holland, Dmitry Sidorenko, and Veronika Eyring
Diajeng W. Atmojo, Katja Weigel, Arthur Grundner, Marika M. Holland, Dmitry Sidorenko, and Veronika Eyring
Metrics will be available soon.
Latest update: 08 Oct 2025
Download
Short summary
This study improves the sea ice albedo parametrisation in the Finite-Element Sea Ice Model by discovering an equation using symbolic regression, an interpretable machine learning method. Leveraging satellite and reanalyses data, our discovered equation identifies high sensitivity to thin snow and the weighted temperature difference between sea ice surface and 2 m air as critical to determine sea ice albedo. Our findings contribute to improving Arctic climate projections and understanding.
Share