Multi-model analysis of the Adriatic dense-water dynamics

Pranić, Petra; Denamiel, Cléa; Janeković, Ivica; Vilibić, Ivica

doi:https://doi.org/10.5194/os-19-649-2023

Articles | Volume 19, issue 3

https://doi.org/10.5194/os-19-649-2023

© Author(s) 2023. This work is distributed under
the Creative Commons Attribution 4.0 License.

https://doi.org/10.5194/os-19-649-2023

© Author(s) 2023. This work is distributed under
the Creative Commons Attribution 4.0 License.

Articles | Volume 19, issue 3

Research article

|

16 May 2023

Research article |

| 16 May 2023

Multi-model analysis of the Adriatic dense-water dynamics

Petra Pranić, Cléa Denamiel, Ivica Janeković, and Ivica Vilibić

Download

Final revised paper (published on 16 May 2023)
Supplement to the final revised paper
Preprint (discussion started on 06 Dec 2022)
Supplement to the preprint

Interactive discussion

Status: closed

RC1:
'Comment on egusphere-2022-1274', Anonymous Referee #1, 05 Jan 2023

The MS presents a helpful analysis of the formation and collection of dense water in the Norther Adriatic during the bora (cold strong winds) events, and some insight into the transport of the newly formed dense bottom water to a greater depth. The authors post-processed previously obtained outputs form four models and calculated basic statistics in order to validate the models against archived observations in relation to the cascading process. They also used a combination of model outputs and observation to analyse the physics of bora-driven cascading.

The structure of the MS is slightly misleading. The Method section describes in much detail the methodology used in previous research, such as model set ups and parameters, which may give a false impression that the models were configured and run as part of this particular study. At the same time, the methods used in the present study are described only embryonically. The MS could be published after major revision to address the comments in this report.

Specific comments

The transport of dense water is the key element of dense water cascades. In the abstract, the authors claim ‘Finally, we analyse in detail the numerical reproduction of the dense water dynamics as seen by the four simulations.’ However, the dense water transport is described only qualitatively in the main text, and the assessment of the transport only appears in the appendix as ‘additional information’. The transport should be given much greater focus.

Lines 18-20 ‘This study confirms that … are prerequisites for appropriate modelling of the ocean circulation in the Adriatic basin’ . This is an overstatement. The study only reveals that one of the four models represents better some variables while other models are better in some other aspects. The paper does not proof that the named parameters are prerequisites for future research. Please re-word.

Introduction. The authors concentrate entirely on the Adriatic Sea. In order to make the results helpful for a wider oceanographic community the authors are advised to place the bora-driven cascading in a broader context. For example to compare briefly with cascades from other shelves, e.g. as discussed in (Ivanov et al, 2003, https://doi.org/10.1016/j.pocean.2003.12.002; GarciaâQuintana et al, 2021. https://doi.org/10.1029/2020JC016951)

Lines 69-70. ‘ ..the river climatology used in previous studies … has been replaced by a new climatology’. Please clarify, what is the difference between ‘old’ and ’new’ and give a reference. It seems from the text in Line 71 that the ‘new’ climatology was used in some ‘previous studies’.

Lines 78-79. ‘the most advanced variational scheme, the Four-Dimensional Variational…’. Which of the many versions and sub-versions of DA schemes is ‘the most advanced’ is a matter of discussion. Please re-word.

Line 81. ‘a 31-year evaluation simulation’ . Please define what is ‘evaluation simulation’. Was it run with or without DA?

Line 85 ‘the newest reanalysis product for the Mediterranean Sea’ . Please give a reference.

Line 95 and Line 98 ‘…between late autumn 2014 and summer 2015’ . ‘…between late November 2014 and mid-August 2015’ . Please exact dates as you do in Line 100.

Line 109. ‘..of various types of models’. This statement is too wide. Based on the results presented in the MS, the authors could only assess, contrast and compare the specific models they used, not the ‘types of models’.

Line 113.’ reanalysis product for the Mediterranean Sea’. Please give the exact product ID. Is it available from CMEMS catalogue? Is it MEDSEA_MULTIYEAR_PHY_006_004 ? If yes, please use the name for product given by the originators, namely ’Med MFC’ to avoid confusion. This product is generated not just by NEMO as stated in Lines 115-16 but ‘The Med MFC physical reanalysis product is generated by a numerical system composed of an hydrodynamic model, supplied by the Nucleous for European Modelling of the Ocean (NEMO) and a variational data assimilation scheme (OceanVAR)’ (https://data.marine.copernicus.eu/product/MEDSEA_MULTIYEAR_PHY_006_004/description)

Lines 113-200. The description of models from previous studies is too extensive, it should be reduced in size and moved from Material and methods to Introduction. Material and methods should present methodology used specifically for this study and in much greater detail than it is now.

Line 137. Please give a reference to the ALADIN/HR atmospheric model.

Lines 143-144. ‘In this model, the horizontal ROMS grid resolution is 2 km and there are 20 vertically spaced sigma levels controlled by the following parameters …’. Ocean model results are strongly dependent on the vertical resolution of a model. This is hinted by the authors in Lines 622-623 ’ In addition, the lack of vertical resolution in the ROMS-full model probably contributes to the improper representation of the dense water dynamics’. Has the sensitivity study been performed by the authors or other researchers? The authors attempt to compare ‘the types of the models’ however the skill of the same model can change significantly when model governing parameters or boundary conditions (e.g. river discharge) are changed. The authors should demonstrate that parameters of the models used for comparison provide the best results within the limitation of the specific ‘type’ of model in order to suggest which ‘type’ of the model is the best.

Lines 201-206. Have you noticed the ‘double penalty effect’ ( see e.g. https://doi.org/10.5194/os-16-831-2020 ) in the higher resolution models? If yes, how it impacts on the final results?

Line 202-203. ‘In order to compare different simulations, model results with grid resolution coarser than 1 km are interpolated to the AdriSCROMS 1 km grid…’ .Does it mean that all model outputs were also interpolated in the vertical to the AdriSCROMS sigma-coordinate grid with 35 levels? Please clarify.

Line 204-206’ For the ocean simulations, MEDSEA, ROMS-hind and ROMS-full results are regridded to 1 km resolution, while for the atmosphere, ERA5, ALADIN/HR-hind, ALADIN/HR-full and AdriSC-WRF results are all regridded to 1 km resolution.’ The use of the word ‘while’ is strange as all model outputs were regridded to the same scale. Please re-word.

Line 214-215. ‘ probability density functions of the biases (i.e., differences) between the results of the simulations and the in-situ temperature and salinity observations’ . This is the core component of methodology and it has to be described in much more detail. How the bias is calculated? If it is the average of all differences at all locations and all times then it would be impossible to calculate the PDF. Or is it an average of daily differences? Or something else?

Line 218. ‘The probability density functions are obtained with a kernel-smoothing method..’ Which kernel was used? What was the size of the smoothing window? The results may be sensitive to these factors. Please give more details here (including a reference) as otherwise your results cannot be replicated. An estimate of the sensitivity of results to the size of the smoothing window will be helpful.

Line 223. ‘minimum turbulent heat fluxes in the atmosphere’. Do you mean ‘minimum downward turbulent heat fluxes in the atmosphere at the sea surface’ Please clarify.

Line 230. ‘bottom PDA time series is presented without the seasonal signal which is removed from the series using the least-squares method.’ Please give more details of how PDAs are calculated otherwise the statement in Line 232 ‘the time evolution of the spatial distributions of the bottom PDAs’ is difficult to comprehend.

Lines 233-234. ‘An additional analysis (only presented and commented in Supplementary Material) quantifies the total daily volume transport of the outflowing dense waters..’ The near-bottom transport of dense waters is a major parameter quantifying the intensity of dense water cascades. Therefore it has to be included in the main text ( both results and discussion) in sufficient detail.

Lines 238-248. The validation of the four models against CTD casts is very helpful. In order to help a reader to interpret the figures given in this section, the methods of calculating biases given in the previous section have to be presented in much more detail. It is advisable to extend the basic stats ( mean and standard deviation) to include more advanced tools of model validation, e.g. Pearson correlation, Willmott skill parameter or Taylor diagram.

Lines 238-255. The four models have different resolutions and some of them may not resolve the processes of the scale of the baroclinic Rossby radius. Please provide a map of Rossby radius for your area. You may wish to use a simplified method presented in Chelton et al . 1998: Geographical variability of the first-baroclinic Rossby radius of deformation. J. Phys. Oceanogr., 28, 433-460.

Lines 301-302. ‘Overall, for all models, maximums of wind stresses are associated with bora events, while downward turbulent heat fluxes seem to be more influenced by the seasonal variations of the sea surface temperature’ . Please clarify the second part of this statement. From the qualitative point of view, stronger and colder winds (the bora) should have a greater influence on downward heat fluxes from the atmosphere to the ocean than smooth and therefore weaker seasonal variations.

Lines 345-416. The presentation of results is mostly concentrated on the atmosphere-ocean heat fluxes and the processes of formation of dense water. However, in contrast to CTD observations the models provide an opportunity to calculate dense water transport, which is a key component of cascading. This sections gives a good quantitative description of heat fluxes, while it describes the transports only qualitatively. This omission has to be rectified.

Citation: https://doi.org/10.5194/egusphere-2022-1274-RC1
- AC1: 'Reply on RC1', Petra Pranić, 01 Mar 2023
  
  Dear Referee #1,
  Thank you very much for your detailed review which will highly contribute to the improvement of our article.
  
  The transport of dense water is the key element of dense water cascades. In the abstract, the authors claim ‘Finally, we analyse in detail the numerical reproduction of the dense water dynamics as seen by the four simulations.’ However, the dense water transport is described only qualitatively in the main text, and the assessment of the transport only appears in the appendix as ‘additional information’. The transport should be given much greater focus.
  Response: The assessment of the transport will be placed in the main article as subsection 3.3.3.
  Lines 18-20 ‘This study confirms that … are prerequisites for appropriate modelling of the ocean circulation in the Adriatic basin’ . This is an overstatement. The study only reveals that one of the four models represents better some variables while other models are better in some other aspects. The paper does not proof that the named parameters are prerequisites for future research. Please re-word.
  Response: The sentences will be modified accordingly.
  Introduction. The authors concentrate entirely on the Adriatic Sea. In order to make the results helpful for a wider oceanographic community the authors are advised to place the bora-driven cascading in a broader context. For example to compare briefly with cascades from other shelves, e.g. as discussed in (Ivanov et al, 2003, https://doi.org/10.1016/j.pocean.2003.12.002; Garcia Quintana et al, 2021. https://doi.org/10.1029/2020JC016951)
  Response: A brief comparison will be added.
  Lines 69-70. ‘ ..the river climatology used in previous studies … has been replaced by a new climatology’. Please clarify, what is the difference between ‘old’ and ’new’ and give a reference. It seems from the text in Line 71 that the ‘new’ climatology was used in some ‘previous studies’.
  Response: A clarification will be added.
  Lines 78-79. ‘the most advanced variational scheme, the Four-Dimensional Variational…’. Which of the many versions and sub-versions of DA schemes is ‘the most advanced’ is a matter of discussion. Please re-word.
  Response: The sentence will be changed accordingly.
  Line 81. ‘a 31-year evaluation simulation’ . Please define what is ‘evaluation simulation’. Was it run with or without DA?
  Response: An explanation of the evaluation simulation will be added.
  Line 85 ‘the newest reanalysis product for the Mediterranean Sea’ . Please give a reference.
  Response: A reference will be added.
  Line 95 and Line 98 ‘…between late autumn 2014 and summer 2015’ . ‘…between late November 2014 and mid-August 2015’ . Please exact dates as you do in Line 100.
  Response: The Authors couldn't find the exact dates as they are not available in the original paper (Vilibić et al., 2018).
  Line 109. ‘..of various types of models’. This statement is too wide. Based on the results presented in the MS, the authors could only assess, contrast and compare the specific models they used, not the ‘types of models’.
  Response: The statement will be changed accordingly.
  Line 113.’ reanalysis product for the Mediterranean Sea’. Please give the exact product ID. Is it available from CMEMS catalogue? Is it MEDSEA_MULTIYEAR_PHY_006_004 ? If yes, please use the name for product given by the originators, namely ’Med MFC’ to avoid confusion. This product is generated not just by NEMO as stated in Lines 115-16 but ‘The Med MFC physical reanalysis product is generated by a numerical system composed of an hydrodynamic model, supplied by the Nucleous for European Modelling of the Ocean (NEMO) and a variational data assimilation scheme (OceanVAR)’ (https://data.marine.cop ernicus.eu/product/MEDSEA_MULTIYEAR_PHY_006_004/description)
  Response: Yes, the reanalysis product is MEDSEA_MULTIYEAR_PHY_006_004. The name of the product stated in Escudier et al. (2020) is „Med MFC“. The Authors agree that „Med MFC“ is the best name in order to avoid any confusion. However, Escudier et al. (2021) referred to the same reanalysis as „MEDREA24“ for the purpose of their study. Thus, the Authors chose the name „MEDSEA“ just for the purpose of this study and would like to keep the name for practical reasons. The product descripton will be changed accordingly.
  Lines 113-200. The description of models from previous studies is too extensive, it should be reduced in size and moved from Material and methods to Introduction. Material and methods should present methodology used specifically for this study and in much greater detail than it is now.
  Response: A summary of the main features of the models will be added in the form of a table in Material and methods, while a short paragraph about the models will be added to the Introduction. Also, the full description of the models will be modified and transfered to the Supplemet S2. The main focus will be placed on the methodology which will be described in greater detail.
  Line 137. Please give a reference to the ALADIN/HR atmospheric model.
  Response: A reference will be added.
  Lines 143-144. ‘In this model, the horizontal ROMS grid resolution is 2 km and there are 20 vertically spaced sigma levels controlled by the following parameters …’. Ocean model results are strongly dependent on the vertical resolution of a model. This is hinted by the authors in Lines 622-623 ’ In addition, the lack of vertical resolution in the ROMS-full model probably contributes to the improper representation of the dense water dynamics’. Has the sensitivity study been performed by the authors or other researchers? The authors attempt to compare ‘the types of the models’ however the skill of the same model can change significantly when model governing parameters or boundary conditions (e.g. river discharge) are changed. The authors should demonstrate that parameters of the models used for comparison provide the best results within the limitation of the specific ‘type’ of model in order to suggest which ‘type’ of the model is the best.
  Response: A sensitivity study has not been performed as we only use already well validated models – i.e., the skills of each model have previously been assessed (= the best parametrization possible are used) and this is not the aim of this study to redo this assessment. However, the forcing used in the different models (rivers, boundary conditions, etc.) are discussed at length in the paper.
  Lines 201-206. Have you noticed the ‘double penalty effect’ ( see e.g. https://doi.org/10.5194/os-16-831-2020 ) in the higher resolution models? If yes, how it impacts on the final results?
  Response: A comment about the „double-penalty effect“ will be added.
  Line 202-203. ‘In order to compare different simulations, model results with grid resolution coarser than 1 km are interpolated to the AdriSCROMS 1 km grid…’ .Does it mean that all model outputs were also interpolated in the vertical to the AdriSCROMS sigma-coordinate grid with 35 levels? Please clarify.
  Response: Model outputs were not interpolated in the vertical, thus the sentence will be modified.
  Line 204-206’ For the ocean simulations, MEDSEA, ROMS-hind and ROMS-full results are regridded to 1 km resolution, while for the atmosphere, ERA5, ALADIN/HR-hind, ALADIN/HR-full and AdriSC-WRF results are all regridded to 1 km resolution.’ The use of the word ‘while’ is strange as all model outputs were regridded to the same scale. Please re-word.
  Response: The sentence will be modified.
  Line 214-215. ‘ probability density functions of the biases (i.e., differences) between the results of the simulations and the in-situ temperature and salinity observations’ . This is the core component of methodology and it has to be described in much more detail. How the bias is calculated? If it is the average of all differences at all locations and all times then it would be impossible to calculate the PDF. Or is it an average of daily differences? Or something else?
  Response: An explanation of the bias calculation will be added to the Methods section.
  Line 218. ‘The probability density functions are obtained with a kernel-smoothing method..’ Which kernel was used? What was the size of the smoothing window? The results may be sensitive to these factors. Please give more details here (including a reference) as otherwise your results cannot be replicated. An estimate of the sensitivity of results to the size of the smoothing window will be helpful.
  Response: A more detailed description of PDFs will be added to the Methods section.
  Line 223. ‘minimum turbulent heat fluxes in the atmosphere’. Do you mean ‘minimum downward turbulent heat fluxes in the atmosphere at the sea surface’ Please clarify.
  Response: Yes, minimum downward turbulent heat fluxes in the atmosphere at the sea surface are presented. The sentence will be modified accordingly.
  Line 230. ‘bottom PDA time series is presented without the seasonal signal which is removed from the series using the least-squares method.’ Please give more details of how PDAs are calculated otherwise the statement in Line 232 ‘the time evolution of the spatial distributions of the bottom PDAs’ is difficult to comprehend.
  Response: The PDA calculation will be described in greater detail in the Methods section.
  Lines 233-234. ‘An additional analysis (only presented and commented in Supplementary Material) quantifies the total daily volume transport of the outflowing dense waters..’ The near-bottom transport of dense waters is a major parameter quantifying the intensity of dense water cascades. Therefore it has to be included in the main text ( both results and discussion) in sufficient detail.
  Response: The section about dense water transport will be included in the main text.
  Lines 238-248. The validation of the four models against CTD casts is very helpful. In order to help a reader to interpret the figures given in this section, the methods of calculating biases given in the previous section have to be presented in much more detail. It is advisable to extend the basic stats ( mean and standard deviation) to include more advanced tools of model validation, e.g. Pearson correlation, Willmott skill parameter or Taylor diagram.
  Response: The validation statistics will be extended with a Taylor diagram and a paragraph will be added to subsection 3.1.
  Lines 238-255. The four models have different resolutions and some of them may not resolve the processes of the scale of the baroclinic Rossby radius. Please provide a map of Rossby radius for your area. You may wish to use a simplified method presented in Chelton et al . 1998: Geographical variability of the first-baroclinic Rossby radius of deformation. J. Phys. Oceanogr., 28, 433-460.
  Response: A map and time series of the baroclinic Rossby radii for the Adriatic Sea will be added to the article.
  Lines 301-302. ‘Overall, for all models, maximums of wind stresses are associated with bora events, while downward turbulent heat fluxes seem to be more influenced by the seasonal variations of the sea surface temperature’ . Please clarify the second part of this statement. From the qualitative point of view, stronger and colder winds (the bora) should have a greater influence on downward heat fluxes from the atmosphere to the ocean than smooth and therefore weaker seasonal variations.
  Response: The statement refers to the fact that in December and January the days are shorter than the rest of the year. Consequently, there are less solar radiations during these months and the sea surface temperature is colder than usual. This seems to affect the values of the downward heat fluxes. The Authors however agree that the above formulation might not be so clear.
  Lines 345-416. The presentation of results is mostly concentrated on the atmosphereocean heat fluxes and the processes of formation of dense water. However, in contrast to CTD observations the models provide an opportunity to calculate dense water transport, which is a key component of cascading. This sections gives a good quantitative description of heat fluxes, while it describes the transports only qualitatively. This omission has to be rectified.
  Response: The Authors agree and a quantitative description of transports will be added to the main article.
  
  Citation: https://doi.org/10.5194/egusphere-2022-1274-AC1
RC2:
'Comment on egusphere-2022-1274', Anonymous Referee #2, 15 Feb 2023

The MS presents a nice analysis of the dynamics of the bora driven circulation in the Norther Adriatic during 2014-2015. The authors work with outputs form four simulations. The statistical study is interesting to validate the numerical results against the real observations.
Although the article is good overall, it needs some improvements especially in the methods used in this study which should be more detailed, but also the circulation of the dense water should be more emphasized.
The last line of the abstract is a bit strong in the sense that the article does not prove that 1-the spatial scale, 2-non hydrostatic, 3-fine resolution and 4-forcing at the open boundary are pre-requisites, it only shows us that one simulation, of the four taken into account, is better than the others and it is closest to the dynamic theory of the behavior of the northern Adriatic basin.
The calculations of the probability density functions should be given in more detail, perhaps a sensitivity study of the method should be mentioned.
Also, the seasonal signal removed from the PDA time series would be shown.
The numerical models used in the study have different resolutions, so a comment regarding the numerical resolution of processes with different Rossby's radius deformation would be important.
Vertical resolution can also play a large role in improperly representing the dynamics of dense water. Therefore, the same model can change significantly as the parameters that govern it vary as well as the boundary conditions. This should be described in more detail even if not verified.
Also, the 31-years simulation should be described in more detail.
The authors are very focused on heat fluxes, but a good model would also allow to estimate the transport of dense water

Citation: https://doi.org/10.5194/egusphere-2022-1274-RC2
- AC2: 'Reply on RC2', Petra Pranić, 01 Mar 2023
  
  Dear Referee #2,
  Thank you very much for your review which will contribute to the improvement of our article.
  
  The last line of the abstract is a bit strong in the sense that the article does not prove that 1-the spatial scale, 2-non hydrostatic, 3-fine resolution and 4-forcing at the open boundary are pre-requisites, it only shows us that one simulation, of the four taken into account, is better than the others and it is closest to the dynamic theory of the behavior of the northern Adriatic basin.
  Response: The abstract will be modified accordingly.
  The calculations of the probability density functions should be given in more detail, perhaps a sensitivity study of the method should be mentioned.
  Response: Calculations of the probability density functions will be described in more detail.
  Also, the seasonal signal removed from the PDA time series would be shown.
  Response: Time series with the seasonal signal are shown on Fig. 7.
  The numerical models used in the study have different resolutions, so a comment regarding the numerical resolution of processes with different Rossby's radius deformation would be important.
  Response: A comment about the resolution and the Rossby radii will be added.
  Vertical resolution can also play a large role in improperly representing the dynamics of dense water. Therefore, the same model can change significantly as the parameters that govern it vary as well as the boundary conditions. This should be described in more detail even if not verified.
  Response: The influence of the vertical resolution will be described.
  Also, the 31-years simulation should be described in more detail.
  Response: The 31-year simulation will be described in more detail.
  The authors are very focused on heat fluxes, but a good model would also allow to estimate the transport of dense water
  Response: More focus will be placed on the transport of dense water.
  
  Citation: https://doi.org/10.5194/egusphere-2022-1274-AC2

Peer review completion

AR: Author's response | RR: Referee report | ED: Editor decision | EF: Editorial file upload

AR by Petra Pranić on behalf of the Authors (08 Mar 2023) Author's response Author's tracked changes Manuscript

ED: Publish subject to minor revisions (review by editor) (13 Mar 2023) by John M. Huthnance

Dear Authors
Thank-you for your revised manuscript. In view of the considerable changes, I expect to send it back to the two reviewers. However, to clarify the manuscript and in some cases your response to them, I think it best if you first take account of my “Detailed comments” below. You should be able to deal with these quickly.
Please note that the main text will eventually be copy-edited (and formatted) but the Supplement will appear exactly as you upload it (no APCs for the Supplement means no Copernicus work on it).
Yours sincerely
John Huthnance (editor)

Detailed comments
Too frequently you use of “low”, “high”, “above”, below” etc. for quantities unrelated to elevation. Some specific cases are noted below but there are many other cases in the text.
Related to this, the negative signs for heat flux values makes difficulties in comparing them as large or small. Fluxes from sea to atmosphere would be positive and comparisons less confusing.
Line 13. Last “the” –> “a” and line 14, “the” –> “a”
Line 18. “The prerequisites” –> “Likely prerequisites”? [The reviewers might still complain that you have not proved that these aspects are necessary as you assume, although most would agree that they are.]
Line 21. “the 31-year long evaluation run of the Adriatic Sea climate model” –> “a 31-year long run of the fine-resolution Adriatic Sea climate model”? (to avoid “evaluation run”, undefined at this stage)

Line 33. “the water masses” –> “water”?
Line 37. Better “NAdDW is known . .”
Line 40. “. . mountain range (e.g.) . .”?
Line 87. “the 31-year” –> “a 31-year” (because not yet defined in this Introduction).
Line 90. Better “several-decade-long”
Line 92. “contrarily to” –> “contrary to” or “as distinct from”

Table 1. You could reduce the line spacing and perhaps the column width within the Table so that it fits better.
Line 140. To satisfy Referee 1 you might add explicitly “Model outputs were not interpolated in the vertical.”
Line 145. Better “. . in this study “dense waters” are defined as having PDAs equal or larger”?
Line 153. Better “Consequently” –> “I.e.” or “That is” or “Specifically” or “More precisely”?

Line 206. “an overestimation of the temperatures in MEDSEA” seems to contradict the MEDSEA median bias 0.00 (line 203).
Line 221. “a higher” –> “larger” and “lower” –> “smaller”. However, larger variability leading to smaller standardised deviations is a puzzle. How are the deviations “standardised / normalised”?
Line 242. “highly” –> “strongly”
Line 243. “return” –> “turn”
Line 250. “above” –> “less in magnitude than”?
Lines 255, 262. “below” –> “larger in magnitude than”?
Line 263. “lower” –> “smaller”.
Line 264. “below” –> “less than”
Line 267. “above” –> “less in magnitude than”?
Line 274. “highest” –> “strongest”. “lowest” –> “smallest” or “weakest”
Line 275. “. . smallest heat loss.”

Line 349. “Consequently” –> “That is”?
Line 355. “6b” –> “7b”
Line 396. Unclear what is “also (shown?) by ROMS-full”
Line 403. “of the” –> “in”?
Lines 462, 464. Omit “down to”?
Lines 466, 468, 473. Better “highest” –> “largest”. Similarly lines 478, 479 (“low” –> “small”).

Line 538. “highest” –> “largest”. Line 539; “lower” –> “less”. Etc.
Line 544. “. . varies between 0.5 and 1.0 km . .” or simply “. . is 0.5-1.0 km . .”
Line 546. “. . very small (less than 500 m) . .”
Lines 547-548. “. . 4 km to extremely small values . .”
Line 553. “below” –> “finer than”
Line 562. “. . (by 20-40 m). . .”
Line 611. “contrarily” –> “contrary”

Lines 647-648. Better “. . (1) kilometre- or finer-resolution atmosphere-ocean models and non-hydrostatic atmospheric models, . .”
Line 679. “Middle-range” –> “Medum-Range”. Line 680 “Forecasts . .”

Figure 6 caption. Better “. . series of wind stresses and daily turbulent fluxes averaged . .” (order as in figure)

Supplement
S2.1, paragraph “Nucleus . .” line 5. “Medium-Range Weather Forecasts . .”
S3.1, line 3 of text: “largest” –> “highest” – conventionally “low” and “high” are used with temperature!
S3.2, page 9 line 5. “form” –> “than in”

Hide

AR by Petra Pranić on behalf of the Authors (16 Mar 2023) Author's response Author's tracked changes Manuscript

ED: Publish subject to technical corrections (05 Apr 2023) by John M. Huthnance

Dear Authors
Thank-you for your re-revised manuscript. I now have comments back from the two referees. Unfortunately (owing to an accident) these are outside the Web system so I am copying the content here.
REFEREE 1. “My only minor comment is that the Conclusion section is too lengthy to my liking. Some of the text can be transferred to the Discussion section. The authors may follow the examples below from JGR-Oceans and JPO:
DOI: 10.1029/2021JC018161
DOI: 10.1175/JPO-D-19-0134.1
DOI: 10.1029/2022JC019055”
My (EDITOR’s) comment: I think all the present content of the “Conclusions” section is valid, but perhaps the content of the bullet points is quite detailed. I don’t think the bullet points are exactly “Discussion” but some of their content might appear usefully as a summary at the end of the “Results” section; content depending on the “Discussion” should remain in “Conclusions”. The “Conclusions” could then point to that summary. I leave this to you (the authors) to decide what to do (if anything).

REFEREE 2. “I really appreciated the considerable editing and rewriting of the article.
My overall comment is very good and the paper can be published.
I would like the authors to meditate on the role of the evaporation fluxes (E-P) in the process of producing dense water which instead seems to be mainly conducted by 'only' cooling.
A small discussion, with a couple of sentences, on the role of numerical conservation of the total volume of the Adriatic respect to the volume of the new and dense water produced, taking into account the different role of evaporation (E-P) compared to cooling, could be a good link also for future articles.”
My (EDITOR’s) comments: It should be clarified whether latent heat flux due to E-P is included in the cooling rates stated. Indeed E-P affects total volume but only a couple of sentences are suggested.

EDITOR’s “Technical corrections”.
Line 193. Please check the depth ranges; there seems to be a gap between 100 and 200 m.
Line 201. Better “lower” –> “poorer”.
Line 258. Better “high” –> “strong”.
Line 263. “-1100” –> “1100”?
Line 544. “3c” –> “14c”

Please consider all of the above noting that on publication all this discussion will be available to readers. The decision “Technical Corrections” means that you should upload your “corrected” manuscript directly to the editorial system for publication. I do not need to see it again. However, the main text will be copy-edited and you should check that this process retains your intended meaning. [The supplement will NOT be edited at all but will remain in the form that you upload it.]
Thank-you for publishing in Ocean Science.
Yours sincerely
John Huthnance (editor).

Hide

AR by Petra Pranić on behalf of the Authors (07 Apr 2023) Author's response Manuscript

Short summary

In this study, we analyse and compare the results of four different approaches in modelling bora-driven dense-water dynamics in the Adriatic. The study investigated the likely requirements for modelling the ocean circulation in the Adriatic and found that a 31-year run of a fine-resolution Adriatic climate model is able to outperform most aspects of the newest reanalysis product, a short-term hindcast and data-assimilated simulation, in reproducing the dense-water dynamics in the Adriatic Sea.