the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 19 Aug 2024
Internal-wave-induced dissipation rates in the Weddell Sea Bottom Water gravity current
Abstract. This study investigates the role of wave-induced turbulence in the dynamics of the Weddell Sea Bottom Water gravity current. The current transports dense water from its formation sites on the shelf to the deep sea and is a crucial component of the Southern Ocean overturning circulation. The analysis is based on velocity records from a mooring array deployed across the continental slope between January 2017 and January 2019 and salinity and temperature (CTD) profiles measured by various ship expeditions. To quantify the importance of internal waves for entrainment into the gravity current along the continental slope, we employ three independent methods for estimating turbulence. First, we use a Thorpe scale approach to compute turbulence from density inversions in density profiles in order to calculate total, process-independent dissipation rate. Second, we apply the finestructure parameterization to estimate wave-induced mixing from vertical profiles. Third, we estimate wave energy levels from moored velocity time series and deduce turbulent kinetic energy dissipation rates by applying a formulation that is at the heart of the finestructure parameterization. On this transect, turbulence is highest on the shelf break and decreases towards the deep sea, in line with a decreasing strength of wave-induced turbulence. We observe a 2-layer structure of the gravity current, a strongly turbulent about 60–80 m thick bottom layer and an upper, more quiescent interfacial layer. In the interfacial layer, internal waves induce an important part of the dissipation rate and therefore to entrainment of warmer upper water into the gravity current. A literature comparison with turbulence measurements up- and downstream of our study site suggests that the question of which turbulent process is dominant may be dependent on the location along the Weddell Sea Bottom Water gravity current. On the shelf, trapped waves are most important, on the slope, we see the effect of breaking internal waves and in the basin, symmetric instabilities are identified as the main driver of turbulence.
- Preprint
(3647 KB) - Metadata XML
- BibTeX
- EndNote
Status: final response (author comments only)
-
RC1: 'Comment on egusphere-2024-2444', Anonymous Referee #1, 22 Sep 2024
Overall, the manuscript is well-written. Their science is solid, and the results based on field observations in the Weddel Sea are pretty interesting. I’m almost satisfied with the current manuscript, but there are some sentences hard to follow. I hope the following specific comments will help improve the manuscript.
Specific comments:
L4-5: The sentence “On the continental shelves…” is hard to follow.
L65-66, Table 1, and L83: Which months were the observations conducted? The observation months should be specified in the main text for readers who want know what the "background mean" stands for. (I found the sentence “all CTD measurements were collected in the same season of austral summer” in L325. I think this information should be mentioned earlier.)
L178-182: It’s hard to follow the methodology because there are vague directives “this” and “it”.
L194: The mooring time series are 2-years long. Does the averaged CTD profiles represent the time-mean density structures?
L202-203: I could not understand the sentence “Buoyancy frequency…”. Could you rephrase it?
L212: I could not follow the sentence “A second fit …”.
L284-285: “This results in …”: I could not follow this sentence because the authors described that the CTD profiles were depth-binned at 1 or 2 dbar resolution (L68-69). They did not mention the vertical resolution of LADCP profiles.
L319-320: The finestructure parameterization can calculate TKE dissipation rates for all profiles. How did the authors use the background dissipation rates?
L466-467: This is not clear form Fig.4.
L620-621: I could not follow the sentence.
Typos:
L147: (Fig. 2a) -> (Fig. 2b)
L585, RHS of the equation (A1): The constant coefficient should be ½, not 2.
Citation: https://doi.org/10.5194/egusphere-2024-2444-RC1 -
AC1: 'Reply on RC1', Ole Pinner, 30 Oct 2024
The authors thank the reviewer for the overall very positive assessment as well as for highlighting of inconsistencies, ambiguities and the occasional lack of clarity.
Overall, we are convinced that all wording criticized for lacking clarity can be rephrased and/or extended to improve the readability of the text, following the guidance of the review.
Following the review, we looked into a possible seasonal bias of the CTD profiles and compared them to vertically interpolated profiles from the mooring records, which provide year-round data. We find that CTD profiles are able to represent the variability of temperature and salinity profiles revealed by the moored instruments. The seasonal bias in availability of CTD profiles (mostly data measured in austral summer) therefore does not translate to a seasonal bias in our dissipation rate estimates. We will point this out in the revised manuscript. However, during this analysis we found 5 (previously undetected) CTD profiles, which deviate suspiciously far from an average profile, that they were deemed physically implausible, and were therefore removed from the analysis. We find that this correction to play a minor role for the derived dissipation.We agree with the reviewer’s statement on the finestructure parameterization. The latter indeed provides complete dissipation rate profiles, any assumption of a background value is not necessary. We assured ourselves, that the actual averaging into bins is done correctly without any false background assumptions. The mistake is therefore solely in the text and not in the calculation, and we are grateful to the reviewer for spotting this. We will remove the erroneous statement.
Following the guidance, we will add the information about the vertical resolution of the LADCP to the text. Additionally, will add a clearer explanation of what the given length scales mean, namely the integration limits in the finestructure method. This will be extended by an added paragraph in the discussion, where we compare our choice of integration limits with the choices in comparable literature.
During this peer-review, we found our integration limits for the shear spectra to be relatively narrow compared to other studies. Therefore, we tested the effect of an increase of the range of wave numbers in the integration from (0,1,2,3) to (0,1,2,…,7,8) in order to be more consistent with previous applications.
The new integration limits of the shear-based finestructure parameterization lead to a slightly improved agreement of shear- and strain-based results. The corresponding figure D1 will be updated accordingly in the revised manuscript. The average shear-to-strain variance ratio R_w value changes from 7.9±10.3 to 7.8±9.2. This therefore does not impact our decision to use the literature value of R_w = 7 for the strain-based formulation. In summary, because the change to the shear integration limits does not affect the strain-based calculations, our interpretation, and discussion of the results remains unchanged.Citation: https://doi.org/10.5194/egusphere-2024-2444-AC1
-
AC1: 'Reply on RC1', Ole Pinner, 30 Oct 2024
-
RC2: 'Comment on egusphere-2024-2444', Anonymous Referee #2, 14 Nov 2024
This manuscript presents vertical mixing estimates near and inside a bottom water gravity current in the Weddell Sea based on three different methods to estimate the mixing, to attempt to identify the contribution of internal gravity waves to the total ocean mixing.
This work is challenging because the mixing parameterizations that are used as not straight forward and involve a lot of tuning and careful choice of parameters. In addition, the errors associated with these methods are difficult to estimate and the authors did not have ‘direct’ microstructure mixing observations to benchmark their mixing estimates.
Nevertheless, the authors have taken great care to explain the three methods well and provided an excellent discussion of assumptions and choices of parameters. While the findings about the relative role of internal waves to the mixing in the gravity current and surrounding waters are limited, the novelty of using the wave energy method to estimate the wave-induced dissipation rate from a velocity time series is one to highlight, and worth further investigating in future work.
The scientific quality of the work is excellent and so is the presentation quality. The significance of the work needs to be better presented.
Main comments:
To compare the three methods, the vertical mean profiles in Fig6 are useful but they only show the data for part of the transect and misses half of the gravity current it seems (text is a bit unclear as to what is included in the data for Fig 6; see other individual comments). It would be useful to have a similar figure showing the mean vertical mixing profiles of all the data, including closer to the shelf past 51.5oW . If you want to separate mixing estimates within and outside the gravity current (because the finestructure param doesn’t work inside the gravity current), then show two mean vertical profile figures.
If you know that in the homogenously mixed BL, the assumptions for the finescale para are violated, why still present the mixing estimates from that method for that layer in some of the figures?
Limitations and resolution of the Thorpe scale method: There are some ways to use your own dataset to work out if it’s the sampling resolution (vertical sampling of CTD) or the instrument accuracy and noise level (CTD resolution) that limit the resolvable density inversions. The parameter R can represent this (Stansfield et al., 2001 and Johnson and Garrett, 2004). Comparing your LT data with au Gaussian fit can also help you estimate how much of the Thorpe scales you have resolved in your data set. This is because the distribution of LT is expected to be lognormal since it is theoretically the result of a multiplicative series of independent events (Stansfield et al., 2001).
Great to see you didn’t estimate diffusivity and thanks for adding the small discussion you provide on this Line 471-474.
Section called ‘Connection to larger scales’: You start this section by saying that you ‘want to set our results in a greater context’. By greater context, I think of the Southern Ocean or global ocean. Restricting that discussion section to the Weddell Sea is not really connecting to ‘larger scales’ in my mind. I would suggest renaming that section to better reflect the content.
The last paragraph in your Discussion brings in the topic of climate change and discusses changes in stratification potentially leading to increased vertical mixing. Currently you seem to summarise the findings from Zhou et al. (2023). Are you able to relate better these statements to your findings? Have you tried to see if the mean mixing along the transect has increased between 1989 and 2022? Or are the uncertainties in mixing estimates too high to be able to do that?
-------------------------------------------------
Specific comments:
Line 34: You might want to remove the word ‘vertical’ since all types of mixing, not just vertical, will entrain ambient waters.
Line 104-105: Add relevant references such as Dillon 1982, Crawford 1986 and Ferron et al., 1998.
Line 118: How did you estimate your density noise level of 3x10-4 kgm-3? Have you considered applying a minimum thickness test based on the Galbraith and Kelley (1996) definition which puts a limit on the resolution of the data set? This minimum height of a density overturn is defined partly on the density accuracy of the instrument.
Line 129-130: Please add a few references here of other people having successfully applied this technique.
Line 136: Is this P=10 value similar to what is usually applied?
Line 154-155: add reference to Fig 2b
Line 315-316: What do you base this statement on? ‘Luckily, these observed higher modes contain the energy that is dissipated locally through turbulence.’ How do you know this?
Line 323-324: Maybe a little more discussion around that choice of neutral density =28.40 for the gravity current definition is needed: is this a common definition used by more than Naveira et al 2002b?
Line 325: if you don’t show the approx. 100 m variation in a table or figure, please add ‘not shown’.
Line 338-342: If this is not shown in a table or figure, please add ‘not shown’
Line 345-346: That single outlier profile looks dubious. Have you got any other CTD data to check the buoyancy frequency profile from another source? Would there be any reason for such large values of dissipation at that place and that time, like increased wind forcing (storm) or something else?
Line 358: ‘around 52oW’ . The elevated mixing in the whole water column is at 53oW on Figure 4. Please either change the value in the text of fix the figure.
Line 375-376: How did you estimate variables like the dissipation rate, inside the gravity current? Did you do it qualitatively ‘by eye’ on the figures or did you quantitatively average values within the core of the current based on a core definition? I suggest you try doing the quantitative approach.
Line 388-389: Did you only use data between 48.5 and 51.5oW for Figure 6? I think that is what you mean by this sentence. If so, please add that info in Fig 6 caption.
Line 400: add a depth range for what you mean by intermittent layer in brackets please.
Line 404: same as above but for ‘interfacial layer’; pls add a depth range.
Line 464: This has been observed before. Add ref and discussion, with here or in your 5.3 ‘relation to other studies’ section, based on existing literature on this topic such as Waterman, S., K. L. Polzin, A. C. Naveira Garabato, K. L. Sheen, and A. Forryan, 2014: Suppression of Internal Wave Breaking in the Antarctic Circumpolar Current near Topography. J. Phys. Oceanogr., 44, 1466–1492, https://doi.org/10.1175/JPO-D-12-0154.1.
Line 469: Add a depth range for what you call the ‘inner water column’ please.
Line 485-494: This section would benefit from being tidied up. Currently not very convincing and unclear what you can actually demonstrate based on your data.
Line 525-530: This paragraph is maybe a bit oversimplified? There are likely some appropriate models that resolve the gravity current and in which the wave propagation would be simulated. Instead of saying it is not possible, maybe say this could be part of future studies when the right tools are identified.
Line 550-558: in this paragraph please better separate your own statements from Zhou et al (2023) findings. Can you better relate what you say here to your own results?
Line 552-553: ‘The parameterization yields results comparable in value to the long-tested method of finestructure analysis.’ This is mostly true but not completely. In the bottom layer, the IDEMIX epsi estimates and the finescale epsi estimates differ significantly in my opinion (compared to the rest of the water column), and in a way that is currently unexplained. I would suggest to temper that statement.
Line 570: add ‘… is complicated by large uncertainties in the mixing estimates, …’
Acronyms throughout: Please define acronyms the first time they are used.
-------------------------------------------------
Technical corrections:
Line 28: there is something missing in that sentence, like a word and it does not make sense. Please fix that sentence.
Line 51-52: Consider rephrasing the beginning of that sentence, which is currently awkward ‘Due to its remote and difficult to access location at high latitudes, …’.
Line 719: Here and elsewhere in the references, the hyperlinks to the datasets on Pangaea currently include a comma (‘,’) at the end of the link, which makes the link invalid when you click on it. Please remove the comma from within the link so the links can be used to access the data.
-------------------------------------------------
Figures:
Figure 1: Very nice. Subplot a and b would benefit from being bigger. Currently it is difficult to look at features on figure 1a as it is too small.
Figure 2b insert: What is the second most prominent frequency that is not labelled, after M2?
Figure 3: Can you remind the reader in the caption what the measurement period is? Is it January 2017 to January 2019?
Figure 4: Nice figure! It would be useful to see a contour of the core of the gravity current based on the mean velocity field shown on Fig 3. Maybe a contour of 0.30 m/s or 0.25 m/s? In the caption, add info about the grey rectangles which probably mean no data available .
Figure 5: Same as above: add a mean velocity contour to show the location of the core of the gravity current. In the caption, add info about the grey rectangles which probably mean no data available .
Figure 6: Nice figure! See main comments for more feedback.
Citation: https://doi.org/10.5194/egusphere-2024-2444-RC2
Model code and software
Analysis code to: Internal-wave-induced dissipation rates in the Weddell Sea Bottom Water gravity current, (Pinner et al., 2024) Ole Pinner https://doi.org/10.5281/zenodo.13134608
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 189 | 52 | 16 | 257 | 6 | 6 |
- HTML: 189
- PDF: 52
- XML: 16
- Total: 257
- BibTeX: 6
- EndNote: 6
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1