Review of revised manuscript "A 2700-year annual timescale and accumulation history for an ice core from Roosevelt Island, West Antarctica" by Mai Winstrup et al.
In their revised manuscript the authors made a serious effort to address the issues raised in the first round of reviews. The manuscript has been restructured for clarity and several potential sources of misunderstanding have been eliminated. In particular the independence of the RICE17 age scale w.r.t the WAIS Divide ice core chronology is now made clear. Likewise the manuscript benefits from several related publications now being published or accessible in their discussion stage. This especially concerns the companion paper of Lee et al.
However, there are still two fundamental issues to be fully resolved before the manuscript can be considered for final publication in CP.
1. The automated annual layer counting by the StratiCounter algorithm is the backbone of the RICE17 chronology, but the authors provide no evidence of its performance in identifying annual layers, especially for the deeper and more thinned core sections. At present, the manuscript only contains one figure (Fig. 4) showing the manual assignment of annual layers in a very shallow section of the core, which is not enough. A separate figure would be needed showing exemplarily the automated layer identification. Of special interest here would be the depth section that has been manually counted to initialize the algorithm, as well as part of the deep sections with thinned layers.
As an added benefit, such a figure would also provide evidence of how clear an annual signal can be counted in the various proxies - here referring to the discussion in the first round of reviews about using BC vs. H+ for counting.
Although the main text would be adequate, I suppose the figure could also be placed in the supplement part S2 if the authors would not like to further increase the volume of the main manuscript.
2. The volcanic identification and matching against WAIS remains elusive to me and needs to be clarified further. The detection of a presumably volcanic signal in the various proxies seems to rest on manually selecting even small peaks in the signals. The manual approach appears somewhat arbitrary by comparison to the automated annual layer counting procedure. This is especially so since the authors are aware of the difficulties associated with the high background of non-volcanic sulphur/sulfate and acidity and state to have removed a number of peaks previously considered as volcanic markers. At present, it appears questionable and unclear from the manuscript how reliably volcanic signals are in fact recorded at RICE. Hence, selecting volcanic signals would benefit from a more quantitative detection technique, or at least from employing clearly stated criteria.
The development of the non-sea-salt conductivity (nss-cond) may provide remedy in this context, but is not convincingly shown to record volcanic signals. For instance, Figure 9 shows three distinct peaks in nss-cond (around 162, 164 and 168 m RICE depth). The other match points, however, have little or no sign of a peak, including the Pleiades tephra horizon. There is an additional peak around 169 m not considered. Moreover, judging from Figure 9 the signals of H+ and ECM appear to perform equal or even better (less noise) in indicating peaks at the presumed volcanic markers - which would argue against the need to use nss-cond. Based on Figure 9, I am not convinced that the acidity signal provides more information than the nss-S signal, other than its higher depth resolution.
The authors are aware of potential pitfalls associated with nss-cond calculated as a secondary quantity, but do not quantify how much (or how little) the signal is affected by these problems. For example, this concerns potential false peaks produced by uncertainty in depth-alignment between Cond and Ca++, as well as peaks through very low Ca++.
Consequently, I believe the volcanic identification and matching needs to be expanded and clarified, ideally including: i) a quantitative approach to peak detection, e.g. by employing a local "peak-over-threshold" criterion, ii) a more convincing demonstration of the new volcanic proxy nss-cond, including addressing the pitfalls of a secondary quantity, and iii) showing a more extended version of Figure 9 as a side-by-side comparison between RICE and WAIS volcanic signals, also at depths not including the distinct tephra marker of Pleiades.
Detailed comments:
Introduction: The companionship to the Lee et al. paper should be mentioned explicitly when stating the scope of this work.
P6 L37-38: " Some uncertain layers were counted as a year in the timescale, while others were not."
This procedure remains unclear. As I am sure the authors are aware, a typical approach here is to count uncertain layers as 0.5+/-0.5 years. Why was this procedure not adopted in this case? It is also not clear what uncertainty estimate was obtained and why this should be a 95% confidence interval.
P8 L5-6: "The calcium and conductivity records frequently displayed multiple peaks per year, limiting their contribution to the annual layer interpretations."
This is another reason why it would be very important to see the performance of StratiCounter in such a case.
P8 L20-23: "Using ... the modelled density profile fits well the observed values".
How were the parameters used here obtained, i.e. why these values selected?
Also, from Figure S2 it can be seen that this statement is true only for the top 50 m. Below that depth there is a substantial misfit in the lower core - including depth of the firn-ice transition. This should at least be mentioned - I assume consequences of the misfit would be negligible if only an extrapolation of near-surface values was concerned.
P9 L15: "estimated using an ice-flow model to match the dated architecture of the Raymond stack"
This is rather vague. What ice-flow model was used? Is this work done by the authors or the approach of Kingslake et al. (2014)? Please clarify.
P9 L19-20: "a linear combination with divide-type velocities weighted by 0.7 and flank-type velocities weighted by 0.3."
Please state why these values were used. Was this an arbitrary selection? If so, would uncertainty in these values contribute to the uncertainty in thinning function? This approach remains unclear.
P9 L29-30: "Except for the uppermost part of the record, uncertainty in the thinning function dominates the total uncertainty, and only this factor will be considered here."
This is an important statement but is not supported by any evidence. Please explain.
P9 L31 ff: This paragraph has some vague statements making it hard to fully comprehend. E.g. Was the near-surface uncertainty in density (it is only extrapolated) considered here? And: "(b) variation of the vertical velocity profile over time in ways not accounted for." E.g. I assume you are referring here to your change in weights of the linear combination?
P10 L3-4: "We suggest to interpret this uncertainty as a 95% (2σ) confidence interval."
It is not clear why this needs to be exactly a 95% confidence level, or why this confinement to 95% is even needed.
P10 L13-14: "The CFA acidity record is driven primarily by the influx of non-sea-salt sulfur-containing compounds, as evident by its high resemblance to the IC non-sea-salt sulfate and ICP-MS non-sea-salt sulfur records in the top part of the core (Fig. 4)."
This may be true, but how can it be reconciled with the later statements about using acidity as a more reliable tracer of volcanic events, e.g. the statements about halide acids from regional volcanism? See comment for P19 L23 below.
P10 L28 and L33: Should be now Figure 9, not 8.
P13 L24: Please give a definition of how the conductivity-to-calcium excess was calculated.
P13 L27-28: "Nevertheless, peaks in the conductivity-to-calcium excess showed high consistency with peaks in total acidity, and it proved to be a reliable tracer for volcanic activity."
Referring to main comment 2., this needs to be shown in more convincing detail. For instance, how does the excess signal look for the recent historical eruptions used in the paper, such as Krakatau?
P13 L36: The Kalteyer, 2015 reference is not accessible without a login from the University of Maine.
P14 L44: "Fig 9a" should be Figure 10.
P15 L40: "It has been suggested that this shift is due to other factors than temperature (Bertler et al., 2018),..."
Vague statement, please clarify what "other factors" are.
P16 L35-36: "Present accumulation rates show a distinct decrease on the downwind (western) side of the ice divide with a gradient of ~5 10-3 m w.e./km yr-1, although muted around the summit area."
Where does this data come from? Is there a missing citation?
P19, L23 ff. "The halide acids are highly soluble, and will be removed from the atmosphere relatively quickly during transport. Hence, they will contribute to increased ice acidity in ice cores located close to the eruption site, whereas only sulfate is deposited from distant volcanic eruptions. By focusing on acidity as volcanic tracer instead of sulfur, the RICE volcanic proxies may thus be more sensitive to regional volcanism than to larger far-field eruptions."
As pointed out above, if acidity is dominated by non-sea-salt sulfur as claimed earlier in the manuscript, how would additional, presumably volcanic acids, leave a detectable imprint? In other words, if the H+ data shown in Figure 9 would be down-sampled to the depth-resolution of nss-S shown in the same Figure, would there be in fact a discernible difference between the two datasets, especially at the locations of the assumed eruptions? The reader is unable to evaluate this based on the presented data in Figure 9. |