Regularization and L-curves in ice sheet inverse models: a case study in the Filchner–Ronne catchment
<p>Over the past 3 decades, inversions for ice sheet basal drag have become commonplace in glaciological modeling. Such inversions require regularization to prevent over-fitting and ensure that the structure they recover is a robust inference from the observations, confidence which is require...
| Main Authors: | , , , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
Copernicus Publications
2023-11-01
|
| Series: | The Cryosphere |
| Online Access: | https://tc.copernicus.org/articles/17/5027/2023/tc-17-5027-2023.pdf |
| _version_ | 1827630743604428800 |
|---|---|
| author | M. Wolovick A. Humbert A. Humbert T. Kleiner M. Rückamp |
| author_facet | M. Wolovick A. Humbert A. Humbert T. Kleiner M. Rückamp |
| author_sort | M. Wolovick |
| collection | DOAJ |
| description | <p>Over the past 3 decades, inversions for ice sheet basal drag have become commonplace in glaciological modeling. Such inversions require regularization to prevent over-fitting and ensure that the structure they recover is a robust inference from the observations, confidence which is required if they are to be used to draw conclusions about processes and properties of the ice base. While L-curve analysis can be used to select the optimal regularization level, the treatment of L-curve analysis in glaciological inverse modeling has been highly variable. Building on the history of glaciological inverse modeling, we demonstrate general best practices for regularizing glaciological inverse problems, using a domain in the Filchner–Ronne catchment of Antarctica as our test bed. We show a step-by-step approach to cost function normalization and L-curve analysis. We explore the spatial and spectral characteristics of the solution as a function of regularization, and we test the sensitivity of L-curve analysis and regularization to model resolution, effective pressure, sliding nonlinearity, and the flow equation. We find that the optimal regularization level converges towards a finite non-zero limit in the continuous problem, associated with a best knowable basal drag field. Nonlinear sliding laws outperform linear sliding in our analysis, with both a lower total variance and a more sharply cornered L-curve. By contrast, geometry-based approximations for effective pressure degrade inversion performance when added to a sliding law, but an actual hydrology model may marginally improve performance in some cases. Our results with 3D inversions suggest that the additional model complexity may not be justified by the 2D nature of the surface velocity data. We conclude with recommendations for best practices in future glaciological inversions.</p> |
| first_indexed | 2024-03-09T14:11:03Z |
| format | Article |
| id | doaj.art-a22d9b38fe1c4ff998360c21495f35e0 |
| institution | Directory Open Access Journal |
| issn | 1994-0416 1994-0424 |
| language | English |
| last_indexed | 2024-03-09T14:11:03Z |
| publishDate | 2023-11-01 |
| publisher | Copernicus Publications |
| record_format | Article |
| series | The Cryosphere |
| spelling | doaj.art-a22d9b38fe1c4ff998360c21495f35e02023-11-29T10:02:13ZengCopernicus PublicationsThe Cryosphere1994-04161994-04242023-11-01175027506010.5194/tc-17-5027-2023Regularization and L-curves in ice sheet inverse models: a case study in the Filchner–Ronne catchmentM. Wolovick0A. Humbert1A. Humbert2T. Kleiner3M. Rückamp4Glaciology Section, Alfred-Wegener-Institut Helmholtz-Zentrum für Polar- und Meeresforschung, Bremerhaven, GermanyGlaciology Section, Alfred-Wegener-Institut Helmholtz-Zentrum für Polar- und Meeresforschung, Bremerhaven, GermanyDepartment of Geosciences, University of Bremen, Bremen, GermanyGlaciology Section, Alfred-Wegener-Institut Helmholtz-Zentrum für Polar- und Meeresforschung, Bremerhaven, GermanyDepartment of Geodesy and Glaciology, Bavarian Academy of Sciences and Humanities, Munich, Germany<p>Over the past 3 decades, inversions for ice sheet basal drag have become commonplace in glaciological modeling. Such inversions require regularization to prevent over-fitting and ensure that the structure they recover is a robust inference from the observations, confidence which is required if they are to be used to draw conclusions about processes and properties of the ice base. While L-curve analysis can be used to select the optimal regularization level, the treatment of L-curve analysis in glaciological inverse modeling has been highly variable. Building on the history of glaciological inverse modeling, we demonstrate general best practices for regularizing glaciological inverse problems, using a domain in the Filchner–Ronne catchment of Antarctica as our test bed. We show a step-by-step approach to cost function normalization and L-curve analysis. We explore the spatial and spectral characteristics of the solution as a function of regularization, and we test the sensitivity of L-curve analysis and regularization to model resolution, effective pressure, sliding nonlinearity, and the flow equation. We find that the optimal regularization level converges towards a finite non-zero limit in the continuous problem, associated with a best knowable basal drag field. Nonlinear sliding laws outperform linear sliding in our analysis, with both a lower total variance and a more sharply cornered L-curve. By contrast, geometry-based approximations for effective pressure degrade inversion performance when added to a sliding law, but an actual hydrology model may marginally improve performance in some cases. Our results with 3D inversions suggest that the additional model complexity may not be justified by the 2D nature of the surface velocity data. We conclude with recommendations for best practices in future glaciological inversions.</p>https://tc.copernicus.org/articles/17/5027/2023/tc-17-5027-2023.pdf |
| spellingShingle | M. Wolovick A. Humbert A. Humbert T. Kleiner M. Rückamp Regularization and L-curves in ice sheet inverse models: a case study in the Filchner–Ronne catchment The Cryosphere |
| title | Regularization and L-curves in ice sheet inverse models: a case study in the Filchner–Ronne catchment |
| title_full | Regularization and L-curves in ice sheet inverse models: a case study in the Filchner–Ronne catchment |
| title_fullStr | Regularization and L-curves in ice sheet inverse models: a case study in the Filchner–Ronne catchment |
| title_full_unstemmed | Regularization and L-curves in ice sheet inverse models: a case study in the Filchner–Ronne catchment |
| title_short | Regularization and L-curves in ice sheet inverse models: a case study in the Filchner–Ronne catchment |
| title_sort | regularization and l curves in ice sheet inverse models a case study in the filchner ronne catchment |
| url | https://tc.copernicus.org/articles/17/5027/2023/tc-17-5027-2023.pdf |
| work_keys_str_mv | AT mwolovick regularizationandlcurvesinicesheetinversemodelsacasestudyinthefilchnerronnecatchment AT ahumbert regularizationandlcurvesinicesheetinversemodelsacasestudyinthefilchnerronnecatchment AT ahumbert regularizationandlcurvesinicesheetinversemodelsacasestudyinthefilchnerronnecatchment AT tkleiner regularizationandlcurvesinicesheetinversemodelsacasestudyinthefilchnerronnecatchment AT mruckamp regularizationandlcurvesinicesheetinversemodelsacasestudyinthefilchnerronnecatchment |