Ultrasonic and seismic constraints on crystallographic preferred orientations of the Priestley Glacier shear margin, Antarctica
- 1Department of Geology, University of Otago, Dunedin, New Zealand
- 2School of Surveying, University of Otago, Dunedin, New Zealand
- 3Antarctica New Zealand, Christchurch, New Zealand
- 4Institute for Marine and Antarctic Sciences, University of Tasmania, Hobart, TAS, Australia
- 5Division of Earth Sciences, Korea Polar Research Institute, Incheon, Republic of Korea
- 6British Antarctic Survey, Natural Environment Research Council, Cambridge, United Kingdom
- 1Department of Geology, University of Otago, Dunedin, New Zealand
- 2School of Surveying, University of Otago, Dunedin, New Zealand
- 3Antarctica New Zealand, Christchurch, New Zealand
- 4Institute for Marine and Antarctic Sciences, University of Tasmania, Hobart, TAS, Australia
- 5Division of Earth Sciences, Korea Polar Research Institute, Incheon, Republic of Korea
- 6British Antarctic Survey, Natural Environment Research Council, Cambridge, United Kingdom
Abstract. Crystallographic preferred orientations (CPOs) are particularly important in controlling the mechanical properties of glacial shear margins. Logistical and safety considerations often make direct sampling of shear margins difficult and geophysical measurements are commonly used to constrain the CPOs. We present here the first direct comparison of seismic and ultrasonic data with measured CPOs in a polar shear margin. The measured CPO from ice samples from a 58 m deep borehole in the left lateral shear margin of the Priestley Glacier, Antarctica, is dominated by horizontal c-axes aligned sub-perpendicular to flow. A vertical seismic profile experiment with hammer shots up to 50 m away from the borehole, in four different azimuthal directions, shows velocity anisotropy of both P-waves and S-waves. Matching P-wave data to the anisotropy corresponding to CPO models defined by horizontally aligned c-axes gives two possible solutions for c-axis azimuth, one of which matches the c-axis measurements. If both P-wave and S-wave data are used, there is one best fit for azimuth and intensity of c-axis alignment that matches well the measurements. Azimuthal P-wave and S-wave ultrasonic data recorded in the laboratory on the ice core show clear anisotropy that matches that predicted from the CPO of the samples. With good quality data, azimuthal increments of 30° or less will constrain well the orientation and intensity of c-axis alignment. Our experiments provide a good framework for planning seismic surveys aimed at constraining the anisotropy of shear margins.
Franz Lutz et al.
Status: final response (author comments only)
-
CC1: 'Comment on tc-2021-382', Nicholas Rathmann, 23 Feb 2022
Dear authors,
Thank you for the opportunity to read this preprint of your work. I’m very intrigued by your methodology, which has made me think a lot about future practical applications. Please find below a few comments that you may (or may not) consider useful for the final version of your manuscript.
1) L23-25: I would suggest writing “can affect”, since basal friction might also exert important resistance (unlike over e.g. ice shelfs). At least from a modeling perspective, it is not yet clear (to me) how little basal resistance is needed before shear margins significantly affect the character of ice-stream flow.
2) Page 7 between eqn. (1) and (2): Maybe it would be useful to explicitly state that you, therefore, are considering a homogeneous, constant-fabric slab/layer (unless I misunderstood)?
3) Discussion section: I would be very interested in a discussion that also treats (possibly briefly) how your VSP method might be applied to regions where a vertical preferred direction CPO is expected. Is it correctly understood that in order to get the largest P and S phase-velocity variations, you would have to sample in a plane perpendicular to the CPO plane-of-isotropy (horizontal plane), and hence consider very wide angle experiments at the surface?
4) Discussion section: I wonder whether you have thought about (or possibly modelled?) what kind of ultrasound sampling scheme could be adopted if the (approximate) CPO is not known a priori? That is, what set of sampling directions would I need to consider to make an unambiguous estimate of the CPO using your method without knowing the plane that the single maximum is in?
Kind regards, Nicholas Rathmann
- AC1: 'Reply on CC1', Franz Lutz, 08 May 2022
-
RC1: 'Comment on tc-2021-382', Thomas S Husdon, 23 Mar 2022
Review for the paper:
Ultrasonic and seismic constraints on crystallographic preferred orientations of the Priestley Glacier shear margin, AntarcticaGeneral comments:
The paper presents a comparison of seismic Vertical Seismic Profile (VSP) data to ultrasonic lab measurements and a model of Crystal Preferred Orientations (CPO). The work is undoubtably a useful and novel contribution to the field. The introduction, methods and results are well described, albeit in a perhaps unconventional structure (methods and results introduced together for each technique). The discussion lacks the detail I expected given the comprehensiveness of the earlier sections of the paper. Indeed, it reads more like an extension of the results rather than a discussion. I’d recommend that the authors dive into a little more detail for why their results show what they show. Similarly, for the conclusions, I’d really like to see a few sentences on how these results have changed the field. They certainly are interesting and rigorous results, so improving the discussion and conclusions shouldn’t require too much work. No more analysis is required, just some additional thought on why the observations and models agree/don’t agree. As such, I am on the fence between suggesting minor or major revisions. Regardless of whether these changes constitute minor or major revisions, it shouldn’t take much time/work to improve this work. With the suggested changes, it will make a useful contribution to the field and so I would then be happy to recommend it for swift publication.
More general comments/further details on the above:Would be nice somewhere to state why horizontal cluster CPO is assumed in this study (valid assumption I think, but needs justifying).
Discussion is a little lacking in “discussion”. Much of it reads more like results. I’d like to see more reasoning for why the results show what they show than is currently included. The results are sound, so the discussion shouldn’t be hard to develop further.
The conclusions would benefit from a final sentence or two summarising the impact on the field and the new insights that their analysis could provide at other glaciers and ice streams. The work is a valuable contribution to the field, so I feel the authors should really sell that briefly in the conclusions.
Specific comments:
L12: “good” – subjective. Consider removing.L17: Large strain compared to what? Comparative language should be avoided unless referring to something else.
L15-29: Very nice opening paragraph!
Figure 1d: Make the labelling of P and S waves clearer (maybe use different colours?). It is obvious to me as a seismologist, but would not necessarily be obvious to other glaciologists.
L53: Is the borehole seismometer three-component? I think it is from what is written, but can you state this explicitly?
L65: Is the sampling interval 1 / the sampling rate? Maybe better to state that you sampled with 8000 Hz rather than the rather odd unit of 0.125 ms.
L66: Is the hard ice surface an ice lense from surface melting or the actual ice column surface? Important to clarify whether you have a firn layer at the site.
L74-75: Does P wave energy in the S wave explicitly preclude shear-wave splitting? If one assumed the P wave energy is linearised (which could be easily tested) then I don’t think it should affect the shear-wave splitting measurement? I would also expect a resonance of the instrument to be linear and so not affect shear-wave splitting?
L87: The underlying cause of the low SNR for zero-offset S waves will be the radiation pattern of the source. For a hammer shot directly above an instrument, you will be in a null for the S wave radiation pattern. I think it would be worth stating this rather than just the SNR effect of this.
L87-89: This sentence doesn’t make grammatical sense. Needs restructuring.
Figure 2: The data show some interesting things. I like the indication of increasing travel time with depth, but wonder whether it might be clearer to present the difference in travel times, relative to a reference profile (e.g. in-flow) as it is difficult to tease out velocity perturbations in this figure.
L97-98: Good that you state that there is no firn layer here. However, I think it would benefit from stating explicitly earlier in the text too (but don’t remove from here as it is key to your assumption used to calculate velocity).
L105: I think the uncertainty in tp and x are probably underestimated, but maybe its not too important.
Equation 3: Factor of 2 shouldn’t be in first term of velocity uncertainty (check the differentiation of Eq. 1). Otherwise equation seems sound.
L108: Not quite sure what Table 2 adds? I don’t quite understand why in one direction one can observe such a high range of velocities? There is no physical reason for this I don’t think, yet your uncertainties suggest that it is real. I think this again points to the uncertainties mentioned above being underestimated.
L109: How does this temperature compare to the ice column? Anisotropy can be very sensitive to temperature. Just worth stating that this is a similar temperature to the ice column in the field, assuming is approximately is.
Figure 3c-e: Nice results! Very clear anisotropic signal.
Figure 4: It would be good to see the corresponding seismic observations also plotted on Figure 4. Unless I am mistaken, you have seismic measurements corresponding approximately to each of your ice core results? Why not plot them together?
Figure 4: Ah, vp does appear to vary azimuthally. This is an interesting result and should probably be highlighted in the abstract. It is potentially expected, but not really been observed at typical seismic wavelengths. This is maybe a reason why shear wave splitting might not work (nullifying one of the assumptions I made in a previous comment).
L151: Need units for the velocity measurements.
L157: Are any bubbles or cracks observed in the samples?
L161: There hasn’t been clear evidence yet of the high degree of seismic anisotropy in the VSP. Could you create a figure similar to figure 4 but for the seismic data? I can see from Figure 2 that there are velocity differences, but it is unclear whether they are coherent/correspond to the various orientations of the ray paths with respect to flow. In summary, To make the statement in L161, I think you need to display your results more clearly.
L165-167: This sentence seems a little random. Is it necessary?
L186: I was really getting convinced by the merits of the CPO modelling, until you mentioned the issues with comparing GHz vs. Hz frequency regimes. I am not well versed in the validity of comparing such disparate frequency regimes, so could you elaborate a little more in the text on why the CPO modelling can be compared to the observations.
Table 4: The misfit uncertainties are larger than the actual values. What does this mean?!
Section 4.3: I have to confess that I got a bit lost here. Lots of data without much text explaining what it all means. It feels tough to read as so many symbols in the prose. I think it could do with some reworking and replacing all the “chi”s with “phase-misfit” or something. Could also consolidate figures as there is a lot of data without much perceived interpretation. The rest of the paper to here reads so well, so worth brushing up on this section.
Section 4.4: Much better than 4.3. Could you restructure Section 4.3 to be similar to 4.4?
Figure 10: Consolidate this figure by putting the measured and modelled data together (as different coloured scatter points). Should make the results look more compelling too.
L284-285: Make this sentence clearer, as the “however” in the middle makes it unclear.
L287: What’s a side minimum? Be more precise.
L290: This sentence should be reworded. You are basically stating that when you remove some of your observations, the results become clearer. You are actually stating, I think, that the data is noisy and so removing so many azimuths reduces the apparent “noise” in the results. Also, you are not really “sparse sampling” as this is a technical term for something completely different (randomly sampling then reconstructing the signal based on assumptions about the frequency content), so try to avoid this terminology.
L296-304: This is more like a results section than discussion. If this remains in the discussion, you need to suggest why the 90 degree results are still very scattered and “unrealistic”.
Figure 13: What do the size of the scatter points represent? It should be mentioned in the caption.
L312-314: Again, more of a results style text than discussion. Why is the VSP CPO modelling ambiguous if only P waves are considered? Is it because P waves shouldn’t exhibit significant anisotropy in this case? I think more detail, i.e. some discussion is needed here.
L331: Do you mean “measured … (CPOs)” or “modelled… (CPOs)”? I think the latter. If so, make the change.
L333-337: Change the word “matches” to “agrees” or similar throughout this paragraph. The data does not exactly match, but does broadly agree.
L342: Change the word “degraded” to something else. Doesn’t do what you’ve done justice and sounds negative rather than positive.
Technical comments:
L10: “that matches well the measurements” -> “that matches the measurements well”L16: “As result” -> “As a result”
L41: comma required before “which”
L92: “times” -> “time” otherwise it is a confusing sentence.
L94: “registered later” -> “slower”?
Figure 2: The figure labels need to be clearer (i.e. a,b,c etc at top of figure and bigger font).
L95: Better not to refer to Equation 1 before it appears in the text. Instead, just state that this is the equation used to calculate seismic velocities. Same for other equations, especially Eq. 3.
L125: Small point, but sampling rates should be in units of Hz. You are stating a sampling interval.
L163: Remove the word “very” – subjective.
Generally, lots of examples of “…. however…” in sentences. Makes it difficult to read. Need a comma or two when using, or ideally only use at beginning of sentences.
Table 4: Could the uncertainties be displaced consistently with the rest of the paper (i.e. + and – rather than []).
L325-326: Sentence could be improved to communicate the idea better. Basically: greater azimuthal sampling of VSP required to improve CPO model constraint.
- AC2: 'Reply on RC1', Franz Lutz, 08 May 2022
- AC4: 'Reply on RC1', Franz Lutz, 08 May 2022
-
RC2: 'Comment on tc-2021-382', Chao Qi, 30 Mar 2022
These authors conducted a combination of field, laboratory and modeling studies on the shear margin of Priestley glacier. In the field, they shot seismic waves in different directions from the borehole and receive the waves in the borehole. In the laboratory, they shot and received ultrasonic waves around the cored sample. Using a forward model, they could find a CPO that fits the observed seismic velocities. This agreement suggests that the CPO and its resulted anisotropy can be constrained by seismic survey. Understanding the CPOs of fast-moving glaciers is crucial to the predictions of the velocities of these glaciers and therefore the balance of ice masses. These authors are establishing a method that could constrain the CPOs beneath the surface without drilling the glaicers. This manuscript continues the work published by the authors in 2020 on JGR Earth Surface. I am happy to see that these authors are making improvements and progresses on this project. I highly recommend publishing this manuscript but would ask the authors to consider the comments listed below, which in my humble opinion could improve the manuscript.
- Writing. I can basically understand the wording of this manuscript, but in some cases there could be misunderstanding, which can be easily clarified by rewriting the sentence. Besides, a large portion of this manuscript was not written in proper scientific written language. I hope the authors could make some improvements. In the manuscript file, I have commented some places I think would cause confusions. Some general issues are included but not limited to the followings. (1) When a multi-word compound noun (or adjective plus noun) is being used as an adjective to describe another noun, then one need to hyphenate it. For example, “vertical seismic profile experiment” in line 43, “Horizontal Cluster CPO type” in line 48, etc. (2) Please try to avoid indefinite antecedent, that is, the ambiguous “this” and “that”. For example, “This precludes …” in line 73. (3) Some verbs often infer logical relationships, so extra cautions are needed. I have commented some places in the manuscript file.
- Some figures lack sufficient descriptions in the main text. Figures often contain a lot of information, useful and useless. It would be better if the authors could lead the readers through a figure, so that we know what the important features are. (1) Figure 4. Most of the descriptions of Figure 4 are in Section 3.2. When I look at Figure 4, the first thing that caught my eye is the velocity difference between samples, the second thing is that all three sample share the same trend of velocity with degree, and the third thing is that the peaks of red and yellow data are mostly at the same angle, while there seems a “phase lag” in angle for the blue data. However, these features were not mentioned in Section 3.2. So I don’t know if I get the correct messages from the figure. I see some of what I concerned were discussed in the Discussion section, but it is necessary to describe these in Section 3.2, which is a “results” section. (2) Figures 7 and 8 are very complicated but are only mentioned in two paragraphs from 208 to 222. The figures are very pretty with cold to warm colors indicating the depth. However, I could not get how to relate the depths with the three model curves. The authors did not mention the use of depth in these figures (apology if the authors actually used the depth info but I missed it). If the depth info is useless, then there is no need to show it in the figure. If the depth is useful, then please describe the results in the text. (3) Figures 11 and 12 should belong to discussion section but were only briefly mentioned in lines 243 to 248. I think it would be interesting to discuss possible causes for the differences between the three.
- Because three aspects of work were combined in this study, the manuscript is structured a bit differently. Based on my understanding, section 2, 3 and 4 are methods and results for field work, laboratory work and modelling, respectively. In section 2, the authors talked about data collection method, travel-time differences and velocity differences. The section is entitled “Seismic anisotropy informed by a vertical seismic profile shooting”. When we talk about seismic anisotropy, two aspects should be included: the radial anisotropy, i.e., velocity difference, and the azimuthal anisotropy, i.e., fast direction. If the authors would like to keep the section title as it, I would ask them to add some descriptions of the fast direction in Section 2.3, so that this subsection is a complete description for anisotropy. If it is not easy to reconstruct the fast direction, I would suggest changing the title of the section and the subsection. The same issue applies to Section 3, in which only wave velocities were described.
- OK, Section 4 is my biggest concern of this manuscript. Here is my understanding for section 4. (1) A CPO pattern was decided: horizontal clustering of c-axes. (2) Calculate the seismic anisotropy induced by this CPO, given different cluster orientation angle and opening angle. (3) Find the best fit for the actual anisotropy. (4) The CPO that gives the best fit is the same the CPO analyzed from the core. The logic is clear. However, I cannot agree with step (1). I think one purpose of this manuscript is to demonstrate that CPOs in a glacier can be constrained by seismic experiments in the field without analyzing the actual core. So I think a reasonable approach is to start with different CPO patterns, like what these authors did in their 2020 JGR paper. In sheared ice, two CPO patterns could form, a single cluster of c-axes and double clusters of c-axes. Given that in the real case the glacial valley may narrow or widen, compressional or tensional deformation may also occur, leading to elongated clusters. So there could be different patterns to start with. I understand that considering multiple clusters may significantly increase the amount of calculations, but it is at least necessary to let readers know why this pattern was decided (Please do not start by mentioning the actual CPO, which is the final answer to check the models). Also, if no other CPO patterns were considered, some discussions are necessary to explain the current limitations of this modeling approach.
- One thing could be helpful for readers to understand this manuscript better is to include some introductions of the wave velocities along different crystallographic axes of ice. In many papers on seismic anisotropy of the upper mantle for example, the wave velocities along different crystal axes of olivine were always introduced, although these people know olivine very well. Ice is not studied as much as olivine for seismic anisotropy. So I think it is necessary to briefly introduce the anisotropy of ice single crystal. Moreover, during the introduction of anisotropy of single crystal, it is also necessary to explain why only c-axes were considered in the modeling, but not a-axes or others.
- Currently, the seismometers were placed vertically in a borehole. The authors also found that no significant bending of the wave pathway (line 100). So the application is a bit different from the seismic survey people do for the crust. I wonder will it be possible to apply this method without a borehole, that is, seismometers on the surface receiving reflected waves? I know this question is steps ahead of the current stage of the method. But the future of this method would be worth discussing.
Line by line comments are marked on the manuscript file.
- AC3: 'Reply on RC2', Franz Lutz, 08 May 2022
- AC5: 'Reply on RC2', Franz Lutz, 08 May 2022
Franz Lutz et al.
Franz Lutz et al.
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
748 | 149 | 21 | 918 | 11 | 11 |
- HTML: 748
- PDF: 149
- XML: 21
- Total: 918
- BibTeX: 11
- EndNote: 11
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1