the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Regularization and L-curves in ice sheet inverse models, a case study in the Filchner-Ronne catchment
Abstract. Over the past three 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 show that the optimal regularization level converges towards a finite non-zero limit in the continuous problem, associated with a best knowable basal drag field. We find that, when inversion results are judged by a metric that accounts for both the variance of the result and the quality of the fit, then they support nonlinear as opposed to linear sliding laws. We also find that geometry-based approximations for effective pressure degrade inversion performance, but that an actual hydrology model may marginally improve performance in some cases. Our results with 3D inversions suggest that the additional model complexity is not justified by the 2D nature of the velocity data, but we are hopeful that inversions of 3D models may be better situated to take advantage of new constraints in the future. We conclude with recommendations for best practices in glaciological inversions moving forward.
-
Notice on discussion status
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
-
Preprint
(39068 KB)
-
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
- Preprint
(39068 KB) - Metadata XML
- BibTeX
- EndNote
- Final revised paper
Journal article(s) based on this preprint
Interactive discussion
Status: closed
-
AC1: 'Comment on egusphere-2023-741', Michael Wolovick, 05 May 2023
There is an error in our equations 1 and 2, defining the SSA stress balance equations and the effective viscosity computation.
Equation 1 should be,
∇·(Hµ (∇u + ∇Tu + 2∇·u I ) ) + τd - τb = 0,
where I is the identity matrix and the other symbols are as defined in the manuscript. This version differs from the version in the manuscript by a removal of a factor of 1/2, by the inclusion of the divergence of horizontal velocity, and by moving the driving and basal stress terms to the left hand side. Meanwhile, equation 2 should be,
µ = 1/2 B ε0(1-n)/n,
where all symbols are defined as in the manuscript. This version differs from the version in the manuscript by a factor of 1/2.
We will correct both of these equations in the next revision.
Much thanks to Christian Schoof for pointing out our mistake in an email.
Citation: https://doi.org/10.5194/egusphere-2023-741-AC1 -
RC1: 'Comment on egusphere-2023-741', Anonymous Referee #1, 07 Jun 2023
- AC2: 'Reply on RC1', Michael Wolovick, 25 Jul 2023
-
RC2: 'Comment on egusphere-2023-741', Anonymous Referee #2, 14 Jun 2023
- AC3: 'Reply on RC2', Michael Wolovick, 25 Jul 2023
Interactive discussion
Status: closed
-
AC1: 'Comment on egusphere-2023-741', Michael Wolovick, 05 May 2023
There is an error in our equations 1 and 2, defining the SSA stress balance equations and the effective viscosity computation.
Equation 1 should be,
∇·(Hµ (∇u + ∇Tu + 2∇·u I ) ) + τd - τb = 0,
where I is the identity matrix and the other symbols are as defined in the manuscript. This version differs from the version in the manuscript by a removal of a factor of 1/2, by the inclusion of the divergence of horizontal velocity, and by moving the driving and basal stress terms to the left hand side. Meanwhile, equation 2 should be,
µ = 1/2 B ε0(1-n)/n,
where all symbols are defined as in the manuscript. This version differs from the version in the manuscript by a factor of 1/2.
We will correct both of these equations in the next revision.
Much thanks to Christian Schoof for pointing out our mistake in an email.
Citation: https://doi.org/10.5194/egusphere-2023-741-AC1 -
RC1: 'Comment on egusphere-2023-741', Anonymous Referee #1, 07 Jun 2023
- AC2: 'Reply on RC1', Michael Wolovick, 25 Jul 2023
-
RC2: 'Comment on egusphere-2023-741', Anonymous Referee #2, 14 Jun 2023
- AC3: 'Reply on RC2', Michael Wolovick, 25 Jul 2023
Peer review completion
Journal article(s) based on this preprint
Data sets
Inverse Model Results for Filchner-Ronne Catchment Michael Wolovick, Angelika Humbert, Thomas Kleiner, and Martin Rueckamp https://zenodo.org/record/7798650#.ZDkxqZFBxH4
Model code and software
Inverse Model Results for Filchner-Ronne Catchment Michael Wolovick, Angelika Humbert, Thomas Kleiner, and Martin Rueckamp https://zenodo.org/record/7798650#.ZDkxqZFBxH4
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
342 | 148 | 25 | 515 | 17 | 12 |
- HTML: 342
- PDF: 148
- XML: 25
- Total: 515
- BibTeX: 17
- EndNote: 12
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1
Cited
1 citations as recorded by crossref.
Michael Wolovick
Angelika Humbert
Thomas Kleiner
Martin Rückamp
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
- Preprint
(39068 KB) - Metadata XML