the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Subglacial valleys preserved in the highlands of south and east Greenland record restricted ice extent during past warmer climates
Abstract. The Greenland Ice Sheet is a key contributor to contemporary global sea level rise, but its long-term history and response to episodes of warming in Earth’s geological past remain uncertain. The terrain covered by the ice sheet comprises ~79 % of Greenland and ~1.1 % of the Earth’s land surface and contains geomorphological records that may provide valuable insights into past ice-sheet behaviour. Here we use ice surface morphology and radio-echo sounding data to identify ice-covered valleys within the highlands of southern and eastern Greenland and use numerical ice-sheet modelling to constrain the climatological and glaciological conditions responsible for valley incision. Our mapping reveals intricate subglacial valley networks with morphologies that are indicative of substantial glacial modification of an inherited fluvial landscape, yet many of these valleys are presently situated beneath cold-based, slow-moving (i.e., non-erosive) ice. We use the morphology of the valleys and our simple ice-sheet model experiments to infer that incision likely occurred beneath erosive mountain valley glaciers during one or more phases of Greenland’s glacial history when ice was restricted to the southern and eastern highlands, and Greenland’s contribution to barystatic sea level was up to +7 metres relative to today. We infer that this valley incision primarily occurred prior to the growth of a continental-scale ice sheet, most likely during the late Miocene (ca. 7–5 Ma) and/or late Pliocene (ca. 3.6–2.6 Ma). Our findings therefore provide new data-based constraints on early Greenland Ice Sheet extent and dynamics that can serve as valuable boundary conditions in models of regional and global palaeoclimate during past warm periods that are important analogues for the 21st Century and beyond.
-
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
(2910 KB)
-
Supplement
(764 KB)
-
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
- Preprint
(2910 KB) - Metadata XML
-
Supplement
(764 KB) - BibTeX
- EndNote
- Final revised paper
Journal article(s) based on this preprint
Interactive discussion
Status: closed
-
RC1: 'Comment on egusphere-2023-2502', Anonymous Referee #1, 20 Nov 2023
This study aims to reconstruct the broad climatic -- ice-sheet conditions which enabled the development (and then preservation) of subglacial landscapes along the south-eastern margins of the present-day Greenland Ice Sheet. The study has two primary foci: first, the identification and mapping of 'relict' subglacial valleys using a combination of optical remote sensing and airborne ground penetrating radar, then second, using the configuration of these valleys to explore through ice sheet modelling how they have been formed and preserved - and to explore constraints on the timing of these processes. In doing so, the authors conclude that these features were formed when the ice sheet had a configuration restricted to only a few percent of its contemporary size and was geographically restricted to the mountainous south-eastern area. They conclude that the features were formed around sometime in the range 2.6-7 Ma, so likely well before the onset of Quaternary glaciation.
My research experience primarily concerns contemporary Greenland surface processes and ice dynamics by fieldwork and remote sensing. I am not active in paleo-glaciology nor in ice-sheet-wide simulations, so I write this review from a more 'generalist' perspective. I want to stress that I welcome rebuttals if the authors feel I've misunderstood their point(s).
Overall I find the manuscript to be in excellent shape. The title is suitable and the paper addresses relevant scientific questions, particularly for this special issue. Substantial contributions on the timing of previous ice sheet extent are reached. The Methods are well described and extensive. I believe that the key areas of potential sensitivity in the modelling have been investigated. I find the results sufficient to support the interpretation. The paper is well-written and overall easy to follow. The figures are of good quality. It was a fascinating read!
I have some specific minor comments as follows:
L119: The authors state that MODIS MoG imagery records the intensity of the reflection of a satellite-emitted radar signal. This is untrue. MODIS is a passive sensor which records visible/near-visible solar radiation reflections.
L125: 'Conversely' is confusing here. Both the preceding sentence and this one focus on shortcomings of the MoG approach, whereas conversely would suggest that we're about to be told something about what it is good for. 'Furthermore' instead?
L136: Please add a further reference that cites the OIB project.
L139-140: Either expand this methods explanation or remove incompletely - there isn't enough information to judge whether the fact that an ML approach was used is important to the present study or not. If it is, expand why, otherwise I suggest removing this methodological detail.
L277: add a reference to Fig. 4g for concerning the thermal state analysis.
L279: (subjective) - I suggest starting a new paragraph here for clarity
L332: is 'mapped mountain valley networks' missing the term 'glacial'? (to distinguish from fluvial)
L480: The first part of this sentence is really the conclusion of the previous paragraph, so perhaps would be better off added there. Then the paragraphs will match points (a) and (b) introduced in L466-469.
L494: would 'also' result in higher rates of mass accumulation and turnover 'there', or similar (i.e. is the intended meaning that the conditions which enable higher rates of mass accumulation in turnover in the EH also cause the ~same conditions in the SH?)
L572-580: I struggled to understand this paragraph. This might be my shortcomings in being able to 'imagine' isostacy, but I think nonetheless that some rephrasing would be beneficial.
Citation: https://doi.org/10.5194/egusphere-2023-2502-RC1 - AC1: 'Reply on RC1', Guy Paxman, 08 Feb 2024
-
RC2: 'Comment on egusphere-2023-2502', Henry Patton, 15 Jan 2024
In this study, Paxman et al. identify subglacial valleys beneath the highlands of southern and eastern Greenland from MODIS imagery and RES data, and use numerical modelling to constrain the likely age of their incision during mountain-scale glaciation.
The majority of the valleys are presently preserved beneath cold-based/non-erosive ice, indicating their formation is linked to glaciological conditions prior to the growth of a continental-scale ice sheet. Likely timeframes proposed include the late Miocene (7-5 Ma) and/or late Pliocene (3.6-2.6 Ma).
Overall I found the manuscript and figures to be of a very high standard, with the arguments and analyses following in a logical manner. Furthermore, the new data-based constraints provided will serve as useful boundary condition data for future Neogene modelling studies of the GrIS. I have only some minor suggestions/questions below that could help clarify some aspects of the manuscript.Â
194: This process for producing a preQuaternary topopgraphy is unclear to me - why are the refilled subglacial valleys included in the isostatic correction if they are subsequently left open? Are any eroded sediments from troughs/marine sectors on the adjacent shelf included in the isostatic response?
344: I think some further clarification on how you calculate these MAD values would be useful e.g., by showing the equation used. I presume its something like (Σ | (Actual) – (Forecast)|)n, but are you keeping the domain size constant when comparing the two binary fields? What is the region being compared? For example, the MAD values in 7u for b&c look identical but the model misfit in terms of extent in reality is very different.
443: Rather state the fraction than force the reader to lookup a figure.
475: Given that ~30% of the glacial valleys are not cold based today and are seemingly being eroded by the present day ice configuration (Figure 5) does this depth contrast between east and south still hold if comparing just the current cold-based valleys? From Fig 6e it seems these coast-facing valleys of the eastern highlands add some skew.Â
478: An interesting observation on consistent maximum valley depths - possibly related to negative feedbacks related to sediment evacuation? (cf. Fig 13 Patton et al., 2016).Â
517: This asymmetric long-term development - coast v inland - is well-recognised on the Norwegian glaciated passive margin too, and a similar hypothesis was put forward by Kleman et al., 2008 (section 4.1), and Hall et al (2013), during mountain-scale glaciation, and could be useful context here.
544: Arguably I would not have included the PDD-driven models in this Fig 9a given the conclusions of the Plach 2018 paper, which also strengthens your argument here.ÂCitation: https://doi.org/10.5194/egusphere-2023-2502-RC2 - AC2: 'Reply on RC2', Guy Paxman, 08 Feb 2024
Interactive discussion
Status: closed
-
RC1: 'Comment on egusphere-2023-2502', Anonymous Referee #1, 20 Nov 2023
This study aims to reconstruct the broad climatic -- ice-sheet conditions which enabled the development (and then preservation) of subglacial landscapes along the south-eastern margins of the present-day Greenland Ice Sheet. The study has two primary foci: first, the identification and mapping of 'relict' subglacial valleys using a combination of optical remote sensing and airborne ground penetrating radar, then second, using the configuration of these valleys to explore through ice sheet modelling how they have been formed and preserved - and to explore constraints on the timing of these processes. In doing so, the authors conclude that these features were formed when the ice sheet had a configuration restricted to only a few percent of its contemporary size and was geographically restricted to the mountainous south-eastern area. They conclude that the features were formed around sometime in the range 2.6-7 Ma, so likely well before the onset of Quaternary glaciation.
My research experience primarily concerns contemporary Greenland surface processes and ice dynamics by fieldwork and remote sensing. I am not active in paleo-glaciology nor in ice-sheet-wide simulations, so I write this review from a more 'generalist' perspective. I want to stress that I welcome rebuttals if the authors feel I've misunderstood their point(s).
Overall I find the manuscript to be in excellent shape. The title is suitable and the paper addresses relevant scientific questions, particularly for this special issue. Substantial contributions on the timing of previous ice sheet extent are reached. The Methods are well described and extensive. I believe that the key areas of potential sensitivity in the modelling have been investigated. I find the results sufficient to support the interpretation. The paper is well-written and overall easy to follow. The figures are of good quality. It was a fascinating read!
I have some specific minor comments as follows:
L119: The authors state that MODIS MoG imagery records the intensity of the reflection of a satellite-emitted radar signal. This is untrue. MODIS is a passive sensor which records visible/near-visible solar radiation reflections.
L125: 'Conversely' is confusing here. Both the preceding sentence and this one focus on shortcomings of the MoG approach, whereas conversely would suggest that we're about to be told something about what it is good for. 'Furthermore' instead?
L136: Please add a further reference that cites the OIB project.
L139-140: Either expand this methods explanation or remove incompletely - there isn't enough information to judge whether the fact that an ML approach was used is important to the present study or not. If it is, expand why, otherwise I suggest removing this methodological detail.
L277: add a reference to Fig. 4g for concerning the thermal state analysis.
L279: (subjective) - I suggest starting a new paragraph here for clarity
L332: is 'mapped mountain valley networks' missing the term 'glacial'? (to distinguish from fluvial)
L480: The first part of this sentence is really the conclusion of the previous paragraph, so perhaps would be better off added there. Then the paragraphs will match points (a) and (b) introduced in L466-469.
L494: would 'also' result in higher rates of mass accumulation and turnover 'there', or similar (i.e. is the intended meaning that the conditions which enable higher rates of mass accumulation in turnover in the EH also cause the ~same conditions in the SH?)
L572-580: I struggled to understand this paragraph. This might be my shortcomings in being able to 'imagine' isostacy, but I think nonetheless that some rephrasing would be beneficial.
Citation: https://doi.org/10.5194/egusphere-2023-2502-RC1 - AC1: 'Reply on RC1', Guy Paxman, 08 Feb 2024
-
RC2: 'Comment on egusphere-2023-2502', Henry Patton, 15 Jan 2024
In this study, Paxman et al. identify subglacial valleys beneath the highlands of southern and eastern Greenland from MODIS imagery and RES data, and use numerical modelling to constrain the likely age of their incision during mountain-scale glaciation.
The majority of the valleys are presently preserved beneath cold-based/non-erosive ice, indicating their formation is linked to glaciological conditions prior to the growth of a continental-scale ice sheet. Likely timeframes proposed include the late Miocene (7-5 Ma) and/or late Pliocene (3.6-2.6 Ma).
Overall I found the manuscript and figures to be of a very high standard, with the arguments and analyses following in a logical manner. Furthermore, the new data-based constraints provided will serve as useful boundary condition data for future Neogene modelling studies of the GrIS. I have only some minor suggestions/questions below that could help clarify some aspects of the manuscript.Â
194: This process for producing a preQuaternary topopgraphy is unclear to me - why are the refilled subglacial valleys included in the isostatic correction if they are subsequently left open? Are any eroded sediments from troughs/marine sectors on the adjacent shelf included in the isostatic response?
344: I think some further clarification on how you calculate these MAD values would be useful e.g., by showing the equation used. I presume its something like (Σ | (Actual) – (Forecast)|)n, but are you keeping the domain size constant when comparing the two binary fields? What is the region being compared? For example, the MAD values in 7u for b&c look identical but the model misfit in terms of extent in reality is very different.
443: Rather state the fraction than force the reader to lookup a figure.
475: Given that ~30% of the glacial valleys are not cold based today and are seemingly being eroded by the present day ice configuration (Figure 5) does this depth contrast between east and south still hold if comparing just the current cold-based valleys? From Fig 6e it seems these coast-facing valleys of the eastern highlands add some skew.Â
478: An interesting observation on consistent maximum valley depths - possibly related to negative feedbacks related to sediment evacuation? (cf. Fig 13 Patton et al., 2016).Â
517: This asymmetric long-term development - coast v inland - is well-recognised on the Norwegian glaciated passive margin too, and a similar hypothesis was put forward by Kleman et al., 2008 (section 4.1), and Hall et al (2013), during mountain-scale glaciation, and could be useful context here.
544: Arguably I would not have included the PDD-driven models in this Fig 9a given the conclusions of the Plach 2018 paper, which also strengthens your argument here.ÂCitation: https://doi.org/10.5194/egusphere-2023-2502-RC2 - AC2: 'Reply on RC2', Guy Paxman, 08 Feb 2024
Peer review completion
Journal article(s) based on this preprint
Data sets
Subglacial valleys preserved in the highlands of south and east Greenland record restricted ice extent during past warmer climates: Datasets Guy Paxman, Stewart Jamieson, Aisling Dolan, Mike Bentley https://doi.org/10.5281/zenodo.10013440
Viewed
HTML | XML | Total | Supplement | BibTeX | EndNote | |
---|---|---|---|---|---|---|
312 | 108 | 20 | 440 | 31 | 10 | 11 |
- HTML: 312
- PDF: 108
- XML: 20
- Total: 440
- Supplement: 31
- BibTeX: 10
- EndNote: 11
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1
Cited
Guy J. G. Paxman
Stewart S. R. Jamieson
Aisling M. Dolan
Michael J. Bentley
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
- Preprint
(2910 KB) - Metadata XML
-
Supplement
(764 KB) - BibTeX
- EndNote
- Final revised paper