the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 14 Jun 2024
Age, thinning and spatial origin of the Beyond EPICA ice from a 2.5D ice flow model
Abstract. The European Beyond EPICA – Oldest Ice consortium is currently conducting an ice core drilling project at Little Dome C (LDC) in Antarctica with the aim of retrieving a continuous ice core up to 1.5 Ma. In order to determine the age of the ice at a given depth, 1D numerical models are often employed. However, they do not take into account any effects due to horizontal flow. We present a 2.5D inverse model that determines the age–depth profile along a flow line from Dome C (DC) to LDC that is assumed to be stable in time. The model is constrained by dated radar internal reflecting horizons. Surface velocity measurements are used to determine the flow line and ascertain the flow tube width, which also allows the model to consider lateral divergence. This new model therefore improves on the results produced by 1D models previously applied to the DC area. By inferring a mechanical ice thickness, the model predicts either the thickness of a basal layer of stagnant ice or a basal melt rate.
Results show that the deepest ice at Beyond EPICA Little Dome C (BELDC) originates from around 15 km upstream. The oldest ice with useful age resolution, i.e. with an age density of 20 kyr m-1, is predicted to be 1.12 Ma at BELDC. Over the LDC area, the 2.5D model predicts a basal layer 200–250 m thick at the base of the ice sheet. Modelled ice particle trajectories suggest that this layer could be composed of stagnant ice, accreted ice or even disturbed ice containing debris. We explore the possibilities, though this is an open question that may only be answered by analysis the Beyond EPICA ice core once it has been drilled. Finally, we discuss in detail a thinning in the basal layer which is less than predicted by the model, as observed in other ice cores. This could mean that modelled ages are significantly over-estimated in the deepest part of the ice column. Given that the age estimate from the 2.5D model is younger than previous estimates, we suggest that horizontal flow is an important factor in this region. However, our model assumes that the flow line features such as flow direction and dome location have not change over the time period considered, which might not be the case.
- Preprint
(5179 KB) - Metadata XML
-
Supplement
(1158 KB) - BibTeX
- EndNote
Status: final response (author comments only)
-
RC1: 'Comment on egusphere-2024-1650', Anonymous Referee #1, 06 Aug 2024
This study present a so called 2.5D inverse model for simulating the age–depth relationships along a flow line from Dome C to Little Dome C in Antarctica. This is an important topic in glaciology, e.g., finding a location of an old ice core, which matchs the scope of The Cryosphere. However, I find it is very difficult for me to understand this manuscript, and at some places I feel not much senses to me.
The main reason is that the authors do not even put some basic explanations of their model. I do not undertand why “a full description of the forward model will be available in a separate subsequent article” is possible if they do not make it clear in this paper? For example, Eqn 1 gives the vertical velocity profile, but what does the horizontal ice velocity (flux) look like? Both horizontal and vertical ice velocity are critical in determing the shapes of age-depth profiles. The authors also do not explain basically the inverse model, e.g., what is p_prior and why they need to use a term “mechanical ice thickess” since they have already the observed ice thickness data - this do not make much sense not using the observed data!
Actually, the age-depth modeling is not new. If the authors can dig a bit in past literatures, they can easily find some nice and important papers, like Greve et al. (2002) and Rybak and Huybrechts (2003). Based on ice flow models and the Eulerian or Lagrangian methods, we can determine the age-depth relationships in a more physical way. I do not see the authors even mention these previous work in the Introduction section. I am not sure what the “forward model” is about if they do not use a physical ice flow model. If I am not convinced the velocity field is correct, it is also hard for me to believe the inversed age-depth profiles are correct either. The inverse model is then just an optimization approach to find the numbers that match the radar chronology record, but without much physics inside.
The writings also make me confused often. For example,
L13, “the 2.5D model predicts a basal layer 200–250 m thick at the base of the ice sheet”, what is “basal layer 200-250 thick”?
L68, “There is no direct thermal representation in the forward model as we use an inferred mechanical ice thickness to determine a basal melt rate”, this sentence comes from no where, and has no references and no explanations.
L80: what is the form of “horizontal flux shape function”?
L97: no definitions for p and Hw
L110: What is delta q and delta x?
Table 1: what is the spatial locations for these 19 IRHs?
These kind of major and minor issues make me feel very difficult to read and understand this manuscript. Thus, I suggest the authors take another careful round of modifications and re-submit the manuscript after they add the necessary inputs.
References:
Ralf Greve, Yongqi Wang, and Bernd Mügge. Comparison of numerical schemes for the solution of the advective age equation in ice sheets. Annals of Glaciology, 35:487–494, 2002.
Oleg Rybak and Philippe Huybrechts. A comparison of eulerian and lagrangian methods for dating in numerical ice-sheet models. Annals of Glaciology, 37:150–158, 2003.
Citation: https://doi.org/10.5194/egusphere-2024-1650-RC1 -
RC2: 'Comment on egusphere-2024-1650', Anonymous Referee #2, 10 Sep 2024
Chung et al. use a 2.5D flowband model and dated internal layers to estimate the age of basal ice along a flowline from Dome C to the Beyond EPICA drill site of Little Dome C. The model finds the best fit spatial pattern of accumulation rate, velocity shape function, and mechanical ice thickness, which can then be translated to either stagnant ice thickness or basal melt rate. The primary result is that the age of basal ice (for Little Dome C this is just above the layer of stagnant ice) is younger than previously suggested, which has implications for the Beyond EPICA ice core.
The modeling is novel and well described. The manuscript is clearly written with a clear and well supported conclusion. There is acknowledgement of the unknown basal process which adds to the excitement of what Beyond EPICA will find when the drilling is completed this season (assuming good fortune). The discussion of ice cores having younger ages than the model predicts is interesting and useful component. The paper is ready for publication, but I hope the authors will consider the including the points below.
The one area that I suggest more of is a discussion of the age results with Lilien et al., 2021. While Lilien et al., 2021 is mentioned in multiple places, it is not clear how this age scale differs. Lilien et al. suggested that a 1.5Ma age, rather than a 1.1Ma age, is likely to be reached, but I think this is not so much a difference in the depth-age relationship as in the definition of "interpretable ice", with this paper using a value of 20 ka/m while Lilien et al. find 14 ka/m. It is a bit difficult for people outside of Beyond EPICA to keep the differences straight, so providing discussions and summaries of the differences more clearly is quite helpful.
A few other minor comments:
- The discussion of ice fabric is appreciated and not including ice fabric in your model is understandable. However, it seems like the effects of fabric could be parameterized through constraints on the shape of the velocity, i.e. the p value. This is obviously future work, but I think considerable progress could be made without trying to model fabric evolution.
- The last sentence of the abstract is frustrating. Why end on such a sour note? Can you finish instead with a concluding note about the progress you have made?
Citation: https://doi.org/10.5194/egusphere-2024-1650-RC2
Viewed
| HTML | XML | Total | Supplement | BibTeX | EndNote | |
|---|---|---|---|---|---|---|
| 333 | 116 | 32 | 481 | 43 | 19 | 20 |
- HTML: 333
- PDF: 116
- XML: 32
- Total: 481
- Supplement: 43
- BibTeX: 19
- EndNote: 20
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1