the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 19 Jan 2026
Improved recovery of sub ice shelf bathymetry from gravity data using an isostatic correction: A case study from the Dotson and Crosson ice shelves, West Antarctica
Abstract. Bathymetry beneath ice shelves is challenging to observe yet is vitally important for modelling how ice sheets will evolve into the future. An alternative to direct observation of bathymetry is to invert airborne gravity data for the bathymetric signal. Appropriate gravity data can be collected via remote sensing above the ice shelf and be used to provide an initial estimate of sub-ice-shelf bathymetry, typically at wavelengths of ~5 km and above. However, lateral variations in density associated with the underlying geology can distort the gravity field biassing the results. We show that techniques which tie inversion results to known bathymetry and topography, although solving some of these issues, may be insufficient in the case of large and deep basins lacking centrally located tie points. Using new direct observations of the Dotson and Crosson ice shelves as a case study, we show that gravity inversion for bathymetry can be improved by considering and removing a model of the gravity field due to crustal isostatic compensation prior to inversion. We finally present our updated and improved bathymetric model for the Dotson-Crosson and Thwaites Glacier Ice Shelf system and discuss where our method can be best applied in future.
- Preprint
(2176 KB) - Metadata XML
- BibTeX
- EndNote
Status: final response (author comments only)
-
RC1: 'Comment on egusphere-2025-6001', Matthew Tankersley & Jörg Ebbing (co-review team), 19 Feb 2026
-
AC1: 'Reply on RC1', Tom Jordan, 26 Mar 2026
We thank the reviewers for their time and comments on our work – they are useful and highlight sections/assumptions which have not been adequately described. Our response is below and our manuscript will be updated in-line with our responses. Reviewers comment in italics.
General comments:
This work provides a novel addition to the field of sub-ice shelf bathymetry modelling. It addresses the major limitation of this technique, the effective removal of the often long-wavelength signal associated with crustal gravity effects. I think the paper can be improved with a better explanation or justification of some of the methods, and some more discussion added about the assumptions made and their limitations.
Calculation of the isostatic Moho: I have a few reservations about the methods used to calculate the isostatic Moho model. I'm not entirely sure why the free-air gravity was used with a Bouguer slab correction to define the topography used in the isostatic calculation. Why not just use an initial (interpolated) topography model like BedMap2? I assume the differences would be small, but you would be removing what seems to be the most tenuous step in this process.
We agree with the reviewer that the initial isostatic calculation is the critical and novel step. However, we do feel that using the free air anomaly as an input is both important and justified, especially when considering the alternatives.
The basic assumption is that the free air anomaly contains the signature of the bathymetry, although distorted by an isostatic signal (i.e. relatively negative values indicate the presence of a basin). Therefore, using the free air as an initial proxy for bathymetry/topography in the isostatic calculation at-least includes the ‘real’ bathymetric pattern.
In contrast an initial interpolated topography model like Bedmap2 includes no information about the bathymetry, therefore can-not be used to predict the pattern of associated isostatic compensation. In-fact BEDMAP2 is a case in point in the Dotson-Crosson area. As there was no knowledge of bathymetry the floor of the basin was simply assumed to be ~200 m below the base of the floating ice shelf. This is totally wrong in this location and made the bathymetry beneath Dotson/Crosson even shallower than the surrounding ocean (and even shallower that what we now know is an incorrect gravity inversion). From an isostatic perspective using this as an initial isostatic model would have been worse than using nothing at all.
If you retain the Bouguer slab correction, it would be good to have some more clarity on how you calculated the topography from the Bouguer slab, as I don't follow how you accounted for locations where the topography is a density contrast between rock and water, rock and ice, or rock and air.
For simplicity we assumed a slab of rock displacing water when calculating the initial bathymetry from the free air anomaly using the Bouguer slab formula. This assumption is correct for the ice shelf regions which are the focus of this study. Away from marine areas topographic amplitude will be over-estimated, by ~6% in ice covered areas and ~40% in ice free areas. As bathymetry/topography away from the ice shelf regions is not a focus of the study we chose to ignore this error. We acknowledge that ignoring the correct factor for calculating the amplitude of the topography in surrounding regions may distort the amplitude of the isostatic correction in the adjacent marine areas. However, there should not be a significant shift in the pattern of the isostatic model from our assumption. As the amplitude of the isostatic correction is further adjusted by subsequent assumptions to give a model which fitted our new observational data we believe that overall our assumptions are valid.
As the isostatic gravity correction is the main novel concept in this paper, I think a slightly more thorough discussion of the assumptions is required so future users of the method understand the limitations. For example, how does the assumption of no lateral strength affect the results? With some lateral strength (as expected for a small study region), the Moho model would likely be more subdued, and therefore less of the shorter-wavelength gravity signal would be accounted for by the Moho model, meaning the inverted bathymetry model may contain more short-to-mid wavelength features than if no lateral strength is assumed.
We agree with the reviewer on the point that the isostatic correction is the key novel part and will update the main text to clarify the method. Specifically, regarding the amount of lateral strength of the lithosphere when considering the isostatic model. The reviewer is correct that a more rigid lithosphere would smooth the response at the Moho. However, in the application of our correction the isostatic correction effectively adds signal (negative amplitude) which increases the depth of the recovered basin. As rigidity increases the isostatic signal is spread over a broader region, removing shorter wavelengths from the correction. The Airy model is therefore the one which provides the greatest possible short wavelength adjustment, and hence the deepest basin.
If your ice shelf basin contained some sediments, how much would this affect the isostatic Moho and the resulting bathymetry?
Sediments within basins will always be an issue for gravity inversion of bathymetry. Generally, the low-density sediment will give an extra negative gravity signal, compared to the water filled bathymetry with no sediment. The resulting bathymetry would then tend to be over-estimated. For many glacially over deepened basins close to the current ice margin this may be less of an issue than might be expected. The glacial erosion which formed them likely removed much of the low-density sediment, and the few thousand years (typically <10 ka) since the ice sheet retreated from its LGM position has not provided enough time for significant deposition, as evidenced for example by the bedrock exposed in many of the in-board parts of the Amundsen Sea Embayment. There are settings where tectonically controlled basins with several km of sediment could be present, potentially with complex changes in isostatic compensation with time, as suggested in the Ross Sea (e.g. Karner et al 2005). Such settings are by definition complex and may cause simple gravity inversion for bathymetry to fail. Such regions could be identified from complex structurally aligned magnetic anomaly patterns, or prior knowledge of regional tectonics.
Gravity processing: It would be nice to have some more details on the gravity reduction. Did you remove the effects of the ice and water loads from the free-air gravity before the inversion, or before the preliminary topography estimate? It would be helpful to understand how the method works to include more figures of the intermediate steps, mainly the isostatic gravity anomaly, but maybe also the preliminary topography and the initial gravity-derived topography. For the various Bouguer calculations, did you use a density of 2670, and if so, why?
Regarding gravity processing: We work directly with the free-air gravity anomaly, without consideration of ice or water loads prior to the inversion. These are in effect the unknown parameters we are trying to recover. Inverting directly from the free air anomaly therefore minimises any corrections which may have been done with inaccurate or incomplete topographic/bathymetric data. We note many authors do use topography/bathymetry as a prior observation, and then refine the bathymetric surface. However, our approach is to create a ‘pure’ gravity derived surface, then adjust it to fit available topographic data.
For density we use 2670 kgm-3 for rock and 1028 kgm-3 for water. The density used is the ‘standard’ value for rock used in the Bouguer correction based on global averages for typical topography. Choice of the density value will directly impact the amplitude of the recovered bathymetry, and this remains an uncertainty in the method. However, use of non-standard density values would also require justification. In the case of the Dotson-Crosson region exposed geology on Bear Peninsula reveals metamorphic and igneous rocks which would be expected to have typical densities of 2670 kgm-3.
Specific comments:
Line 39: bathymetry only dominates for the free-air anomaly, so I would either change gravity to free-air gravity or change dominates to contains.Will amend
Line 64 and 68: with "initial" clarify if you are talking about past published inversion (and which ones, Jordan et al 2020?) or the first attempt of your inversion in this paper.We are referring here to the Jordan et al 2020 and will clarify this point.
Line 70: possible rewording: which is changing the observed gravity -> which is a component of the observed gravityWill amend this line.
Line 71: reword for clarity: The underestimate of the bathymetry suggests that the expected negative gravity anomaly associated with deep bathymetry is being offset by something causing a positive gravity anomaly.Will re-word as recommended.
Fig 3: can you include flight lines over the gravity data? For the gravity, since it is a free air anomaly, the sign is important. Could you use a linear diverging colormap (i.e., blue to red) that is centered on 0 and has the same positive and negative limits (-70 to +70 mGal) just like the colormap in Fig 2a. I'd also suggest a linear colormap for the topography in b, the current colormap exaggerates the difference between 0 and 100 m, and 0 m is relatively meaningless in Antarctica since neither groundingline or coastline is at 0 m.We will add the flight lines of the used gravity data to Fig. 3a.
We have adjusted the colour scale for the free air gravity anomaly as suggested by the reviewers.
We have retained the colour scale on the topography image. This colour bar is derived from a standard GMT colour bar optimised for display of bathymetry and topography datasets when split linearly between its end points. While the scale it is not strictly linear it has some distinct advantages. First it conveniently visually breaks the bathymetry into three broad regions, deep (<~-1000, intermediate, and shallow >~-500), without using red/green distinction. Secondly it places a distinct colour switch at zero. Although disputed by the reviewers, we feel this is a useful distinction as it separates areas which cannot un-ground from those vulnerable to marine incursion. We agree the current coast/grounding line will often ignore this specific depth contour, but we feel it is an important distinction when considering the general pattern of past/future ice sheet evolution and vulnerability.
Line 103: It would be good to have additional details on the gravity processing. What was the general flight line spacing, speed, altitude etc. Cross over errors? Levelling procedure. Was the free-air anomaly corrected for the gravity effects of the ice sheet or the water thickness? In typical inversions, the water layer is included in the initial model and therefore accounted for. But with your technique, I don't see where the gravity effect of the water column is accounted for.
We will update our description of the free-air data collection.
Regarding subsequent correction for ice sheet or water/ice thickness. We chose not to conduct any correction for water or ice thickness prior to inversion. These values (bathymetric depth) are the key unknown values in this region and the aim of the inversion is to constrain these. In our inversion we consider the full free air anomaly is due to an interface between rock and water. The full water column thickness is therefore implicitly included in the inversion.
Line 147: Maybe replace reference elevation model with starting elevation model for clarity.
We consider “reference elevation” is the correct term. We do not adjust these elevations, as they are derived from direct observations. The gravity derived surface is ‘twisted’ to fit this reference data.
Line 165: For such a small region it is harder to justify assuming 0 lateral strength in the lithosphere. Given some lateral strength is expected, this would essentially smooth out your isostatic Moho. A smooth Moho would account for less short wavelength gravity variations, leaving those shortwavelength anomalies within your signal used in the inversion. Therefore a smoother Moho will result in more variation in your topography and vice versa. Due to this, your assumption about lateral strength has a direct effect on the amplitude of your bathymetry variation. By assuming no strength, you are putting much more of the short-wavelength variation into Moho, and therefore less into the bathymetry. While I agree it is a fine assumption, I think this aspect of it is important for you to state.As the reviewers note increased crustal rigidity would flatten the predicted Moho. However, we disagree that increasing rigidity (flattening the Moho) will increase the recovered bathymetry. Because the isostatic gravity model is subtracted from the free air anomaly the higher the amplitude of the isostatic effect beneath the basin the more negative the corrected gravity anomaly will become, and hence the deeper the predicted depth of the basin will be. At the most extreme for an infinitely rigid plate there would be no isostatic correction and the result would match that of Jordan et al 2020, which underestimates the real bathymetry by ~400 m.
Line 173: Please describe how you did the Bouguer slab transformation. Did you take free air gravity (not corrected for ice or water) and divided it by (2 pi G (2670 (crust) - 1 (air)). Wouldn't you need to use a spatially variable density contrast, where it is 2670 - 1 on land, and 2670-1040 for regions of topography<0. By using density of air instead of water, you are increasing the density contrast, which should result in more subdued topography, not amplified as I think you stated. Where did you get 1.6 from?For the Bouguer slab correction we used the effective contrast of rock against water, as this is the density across the boundary we are trying to recover. In a revision of the paper we will better acknowledge in the text that this distorts the predicted topographic amplitude into onshore areas, but as these are not the focus of the study, we consider it a reasonable simplification.
At this point there is a complexity in our method/logic which we acknowledge was poorly explained.
The preliminary topographic estimate is calculated using a density of rock against water. This maximises the amplitude of the topography/bathymetry. However, we know this topographic amplitude is an underestimate (from observations). Calculating the isostatic correction from topography with underestimated amplitude will give an incorrect (too small) correction. To amplify the isostatic effect of the topography we arbitrarily assume the topography recovered above is displacing air rather than water for the isostatic calculation. This amplifies the recovered Moho and associated correction by ~1.6 times. We chose this adjustment as it provided recovered bathymetry which fitted the observed data well. However, we acknowledge that the scaling of the isostatic correction is an assumption which may well be location dependent.
Line 174: The max underestimate is ~400 m, and only in the ice shelf region, as other regions are much better constrained. So does a maximum of 400 m really justify creating a new topography from the Bouguer slab equation?The fact that the bathymetry shows such a significant and one-sided error supports re-calculation of the bathymetry. An alternative would be to use the topography/bathymetry from Jordan et al 2020 as a starting point. However, we chose to re-start the inversion procedure from scratch, rather than try and adjust a partially successful inversion.
Fig. 4: replace isostatic moho gravity anomaly with isostatic moho gravity effect.We will update the test in line with this suggestion.
Fig. 7: for b, you should use the same magnitude of upper and lower limits for the colorbar, and currently it distorts it to make it seem like the grid is dominated by negative values.The negative part of the colour bar is mislabelled – the colour-bar is symmetrical about zero between -400 and 400. We note the errors are dominated by negative values because isostatic compensation was not considered in the original Jordan et al 2020 model and the pattern of isostatic compensation in the region is complex. It is interesting to note that where isostatic pattern is simpler adjusting the bathymetry based on observations alone works well.
Line 284: Change from Moho gravity anomaly to Moho gravity effect, or isostatic gravity effect.Will update revised paper.
Technical corrections:
We will update the revised paper in-line with these technical suggestions.
Ensure consistent use of either sub ice shelf or sub-ice-shelf in title and text
Line 25: add word "bathymetry": Using new direct bathymetry observation ... to clarify.
Line 75: by basin do you mean ice shelf region?
Line 94: add 'past' before gravity-derived to clarify these are not the inversion results from this paper.
Line 170: Mantle should be lowercaseCitation: https://doi.org/10.5194/egusphere-2025-6001-AC1
-
AC1: 'Reply on RC1', Tom Jordan, 26 Mar 2026
-
RC2: 'Comment on egusphere-2025-6001', Mirko Scheinert, 30 Mar 2026
The paper presents a case study on some of the currently most-discussed issues in Antarctic sciences – how to determine sub-glacial / ocean bottom topography especially in the transition zone from the Antarctic continent to the ocean. Reliable and precise knowledge of this topography has utmost importance for a number of questions, such as for the control of ocean water inflow beneath ice shelves (as the authors stated) but also, for example, for the definition of fluxgates for the determination of the ice-mass loss across the grounding line (in the input-output method).
The authors present a method to improve estimates of ocean-bottom topography beneath ice shelves described in an earlier paper (Jordan et al. 2020). They show that the inversion results can be improved considering isostatic compensation which accounts for larger amplitudes as well as for larger variations (gradients) in the working area. In this way, long-wave errors can be mitigated. As such, the paper provides a valuable method to re-assess and improve ocean-bottom topography resulting from gravity inversion (where no direct measurements are available).
However, there are a number of substantial issues such that I suggest a major revision of the manuscript.
First of all, the paper should be fully self-explanatory without the need to read previous paper(s) of the same author(s), which it isn’t. Of course, it should not be necessary to repeat entire sections of previous publications, but a short and concise description would be advisable where needed. For instance, it is not clear, how the free-air gravity anomalies were formed and which surface they refer to.
Second, in my view the paper is not that well written. It lacks thoroughness in many formulations, and consistency in using special terms. For instance, the usage of the terms “error”, “accuracy” and “uncertainty” seems to be arbitrarily.
Third, to me it is surprising that there isn’t any formula given at all. I think it would be much easier to explain and clarify relations when respective formulas are used.
Fourth, I miss an overall discussion of the state of the art, how the presented investigation fits to current research and results obtained by other authors. The resulting product could be compared to available products on sub-glacial and ocean bottom topography to assess similarities and discrepancies, e.g. Pritchard et al. 2025 (BedMap3), Morlighem et al. 2025 (BedMachine Antarctica, version 4), Ockenden et al. 2026 (subglacial topography inferred from ice flow perturbation analysis) or Charrassin et al. 2025.
In this way, also the structure of the paper needs to be improved. This refers especially to Section 1 “Introduction”. There, some general background is given only in the first short paragraph, whereas the largest part deals with the discussion of the results inferred so far by the same author (Jordan et al. 2020) and the discrepancies detected when new direct observational data have become available. This discussion should not be part of an introduction.
A further issue I would like to discuss is the usage of the term “bathymetry”. I know that it is widely used as synonym for “ocean-bottom topography”. In my opinion, however, “bathymetry” refers to the measurements, i.e. to the applied methods/techniques to measure ocean-bottom topography – such as “gravimetry” refers to measurement techniques to observe gravity or gravity differences. Terms such as “ocean-bottom topography” or “seafloor topography” should be preferred.
Detailed remarks
- Please check orthography to be consistent throughout the paper – e.g. “sub ice shelf bathymetry” in the title, “sub-ice-shelf bathymetry” at line 22; “free air gravity anomaly” vs. “free-air gravity anomaly”, etc.
- Line 24: “bathymetry and topography” – see my remark above; clarify if you mean both seafloor and/or land surface topography and/or subglacial topography.
- Line 73: Which “gravity anomaly”? Be consistent to use the correct term(s) throughout the manuscript.
- Line 80/81: “mostly likely” – please correct.
- Line 94 (Fig. 2): Subfig. 2a presents a difference of the previous and the new ocean-bottom topography. Please state, how the difference was formed (w.r.t. A minus B). Moreover, I wouldn’t call it “error” – if you use the term “error” it implies that one of the datasets represents the truth which I doubt is the case. The same holds true for Subfigs. 2b and 2c. (difference vs. error). Please name the abscissa axes in both subfigures (e.g. “profile length”).
- Line 116: What does “positional accuracy” mean? Is it 3D or 2D (horizontal position)? Again, see remark on the usage of the terms “accuracy” and “error”. Maybe you mean “precision” in this case? “Relative errors” (or precision) should not have a unit (it can be given as percentage). Where are the values (“<40 m”) coming from?
- Lines 120-122: Term “accuracy”? What do you like to state with “the horizontal accuracy is relatively high (<1 km) compared with the resolution of the gravity data (> 5 km)”?
- Line 128: “Uncertainty”? (see remark above)
- Line 146: Where is this value (37.2 m) coming from? How was it calculated?
- Line 147: Here you introduce the yet new term “elevation”? Please clarify – are “surface elevation” and “ocean-bottom topography” compiled into one dataset? If done so, why does it make sense? How was the conversion from raster to point data done?
- Line 150: How was the “weighted mean” calculated?
- Lines 165…: Please increase clarity when discussing the application of the isostatic correction. Do you imply isostatic compensation in this region, such that the measured gravity does not contain sufficient signal generated by topography? How do you come to the assumption to use the density of air (hence ~0) instead of water – is this just “try and error”? What term is really calculated in the end – a correction due to the Airy isostatic model, an “isostatic Moho gravity effect”, …? (Correction and effect should normally differ by the sign, a correction being the negative of an effect.)
- Line 173-174: What do you mean with “…a simple Bouguer slab transformation of the observed free air gravity anomaly to elevation”? Please clarify.
- Line 171: correct “O’Donnell et al., 2019” (capital D).
- Line 177…: Why can you use Fourier Transform to calculate the isostatic correction?
- Line 180: How was the inversion done? A short paragraph should be sufficient, where you may refer to the previous paper for details.
- Line 210 / Fig. 5: In the caption, it is stated that “sub-shelf bathymetry” is shown. However, what about the land-based ice regions?
- Line 230 / Fig. 6: Again, instead of “error” I would use the term “difference” (and which quantity has been subtracted from which?). Explain what the lines and yellow triangles denote.
- Line 233: correct “in the and Crosson”.
- Line 242: What do you understand by “model error” and “standard deviation” in this context? Where is the value “90 m” coming from?
- Line 243: Correct “un-resolved” --> “unresolved”, and “short wavelength bathymetric features” and “… topography” --> “short-wavelength” (with hyphen).
- Line 244: What do you mean by “intrinsic errors in the gravity observations”? What uncertainty measure are they given?
- Line 245: Correct “may-not” --> “may not”
- Line 246: Correct “Tankersley et al. 2025” --> “Tankersley et al. (2025)”, same in line 273
- Line 249: Correct “un-known” --> “unknown”
- Line 280 / Fig. 7: Again, a new naming “Moho gravity anomaly”??
- Line 289: Correct “Free air gravity anomaly” --> “free-air gravity anomaly”.
Citation: https://doi.org/10.5194/egusphere-2025-6001-RC2
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 172 | 86 | 26 | 284 | 24 | 33 |
- HTML: 172
- PDF: 86
- XML: 26
- Total: 284
- BibTeX: 24
- EndNote: 33
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
General comments:
This work provides a novel addition to the field of sub-ice shelf bathymetry modelling. It addresses the major limitation of this technique, the effective removal of the often long-wavelength signal associated with crustal gravity effects. I think the paper can be improved with a better explanation or justification of some of the methods, and some more discussion added about the assumptions made and their limitations.
Calculation of the isostatic Moho: I have a few reservations about the methods used to calculate the isostatic Moho model. I'm not entirely sure why the free-air gravity was used with a Bouguer slab correction to define the topography used in the isostatic calculation. Why not just use an initial (interpolated) topography model like BedMap2? I assume the differences would be small, but you would be removing what seems to be the most tenuous step in this process. If you retain the Bouguer slab correction, it would be good to have some more clarity on how you calculated the topography from the Bouguer slab, as I don't follow how you accounted for locations where the topography is a density contrast between rock and water, rock and ice, or rock and air. As the isostatic gravity correction is the main novel concept in this paper, I think a slightly more thorough discussion of the assumptions is required so future users of the method understand the limitations. For example, how does the assumption of no lateral strength affect the results? With some lateral strength (as expected for a small study region), the Moho model would likely be more subdued, and therefore less of the shorter-wavelength gravity signal would be accounted for by the Moho model, meaning the inverted bathymetry model may contain more short-to-mid wavelength features than if no lateral strength is assumed. If your ice shelf basin contained some sediments, how much would this affect the isostatic Moho and the resulting bathymetry?
Gravity processing: It would be nice to have some more details on the gravity reduction. Did you remove the effects of the ice and water loads from the free-air gravity before the inversion, or before the preliminary topography estimate? It would be helpful to understand how the method works to include more figures of the intermediate steps, mainly the isostatic gravity anomaly, but maybe also the preliminary topography and the initial gravity-derived topography. For the various Bouguer calculations, did you use a density of 2670, and if so, why?
Specific comments:
Line 39: bathymetry only dominates for the free-air anomaly, so I would either change gravity to free-air gravity or change dominates to contains.
Line 64 and 68: with "initial" clarify if you are talking about past published inversion (and which ones, Jordan et al 2020?) or the first attempt of your inversion in this paper.
Line 70: possible rewording: which is changing the observed gravity -> which is a component of the observed gravity
Line 71: reword for clarity: The underestimate of the bathymetry suggests that the expected negative gravity anomaly associated with deep bathymetry is being offset by something causing a positive gravity anomaly.
Fig 3: can you include flight lines over the gravity data? For the gravity, since it is a free air anomaly, the sign is important. Could you use a linear diverging colormap (i.e., blue to red) that is centered on 0 and has the same positive and negative limits (-70 to +70 mGal) just like the colormap in Fig 2a. I'd also suggest a linear colormap for the topography in b, the current colormap exaggerates the difference between 0 and 100 m, and 0 m is relatively meaningless in Antarctica since neither groundingline or coastline is at 0 m.
Line 103: It would be good to have additional details on the gravity processing. What was the general flight line spacing, speed, altitude etc. Cross over errors? Levelling procedure. Was the free-air anomaly corrected for the gravity effects of the ice sheet or the water thickness? In typical inversions, the water layer is included in the initial model and therefore accounted for. But with your technique, I don't see where the gravity effect of the water column is accounted for.
Line 147: Maybe replace reference elevation model with starting elevation model for clarity.
Line 165: For such a small region it is harder to justify assuming 0 lateral strength in the lithosphere. Given some lateral strength is expected, this would essentially smooth out your isostatic Moho. A smooth Moho would account for less short wavelength gravity variations, leaving those shortwavelength anomalies within your signal used in the inversion. Therefore a smoother Moho will result in more variation in your topography and vice versa. Due to this, your assumption about lateral strength has a direct effect on the amplitude of your bathymetry variation. By assuming no strength, you are putting much more of the short-wavelength variation into Moho, and therefore less into the bathymetry. While I agree it is a fine assumption, I think this aspect of it is important for you to state.
Line 173: Please describe how you did the Bouguer slab transformation. Did you take free air gravity (not corrected for ice or water) and divided it by (2 pi G (2670 (crust) - 1 (air)). Wouldn't you need to use a spatially variable density contrast, where it is 2670 - 1 on land, and 2670-1040 for regions of topography<0. By using density of air instead of water, you are increasing the density contrast, which should result in more subdued topography, not amplified as I think you stated. Where did you get 1.6 from?
Line 174: The max underestimate is ~400 m, and only in the ice shelf region, as other regions are much better constrained. So does a maximum of 400 m really justify creating a new topography from the Bouguer slab equation?
Fig. 4: replace isostatic moho gravity anomaly with isostatic moho gravity effect.
Fig. 7: for b, you should use the same magnitude of upper and lower limits for the colorbar, and currently it distorts it to make it seem like the grid is dominated by negative values.
Line 284: Change from Moho gravity anomaly to Moho gravity effect, or isostatic gravity effect.
Technical corrections:
Ensure consistent use of either sub ice shelf or sub-ice-shelf in title and text
Line 25: add word "bathymetry": Using new direct bathymetry observation ... to clarify.
Line 75: by basin do you mean ice shelf region?
Line 94: add 'past' before gravity-derived to clarify these are not the inversion results from this paper.
Line 170: Mantle should be lowercase