the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Impact of ocean vertical mixing parameterization on Arctic sea ice and upper ocean properties using the NEMO-SI3 model
Abstract. We evaluate the vertical turbulent kinetic energy (TKE) mixing scheme of the NEMO-SI3 ocean–sea ice model in sea ice-covered regions of the Arctic Ocean. Specifically, we assess the parameters involved in the TKE mixed layer penetration (MLP) parameterization. This ad-hoc parameterization aims to capture processes like near-inertial oscillations, ocean swells, and waves that impact the ocean surface boundary layer, often not well-represented in the default TKE scheme. We evaluate this parameterization for the first time in three regions of the Arctic Ocean: the Makarov, Eurasian, and Canada Basins.
We demonstrate the strong effect of the scaling parameter that accounts for the presence of sea ice. Our results confirm that the TKE MLP must be scaled down below sea ice to avoid unrealistic deep mixed layers. The other parameters evaluated are the percentage of energy penetrating below the mixed layer and the length scale of its decay with depth. All these parameters affect the mixed layer depth and its seasonal cycle, the surface temperature and salinity, as well as the underlying stratification. Shallow mixed layers are associated with stronger stratification and fresh surface anomalies, and deeper mixed layers correspond to weaker stratification and salty surface anomalies.
Notably, we observe significant impacts on sea ice thickness across the Arctic Ocean in two scenarios: when the scaling parameter due to sea ice is absent and when the TKE mixed layer penetration process vanishes. In the former case, we observe an increase of several meters in the mixed layer depth together with a reduction in sea ice thickness ranging from 30 to 40 centimeters, reflecting the impact of stronger mixing. Conversely, in the latter case, we notice that a smaller mixed layer depth is accompanied by a moderate increase in sea ice thickness, ranging from 10 to 20 centimeters, as expected from a weaker mixing. Additionally, inter-annual variability suggests that experiments incorporating a scaling parameter based on sea ice concentration display an increased mixed layer depth during periods of reduced sea ice, which is consistent with observed trends. These findings underscore the influence, through specific parameterizations, of enhanced ocean mixing on the physical properties of the upper ocean and sea ice.
- Preprint
(11268 KB) - Metadata XML
- BibTeX
- EndNote
Status: open (until 10 Jul 2024)
-
RC1: 'Comment on gmd-2024-49', Anonymous Referee #1, 11 May 2024
reply
This is a useful modelling study with a finding that differences in ocean stratification due to TKE MLP parameters have significant implications for the Arctic freshwater content. For global NEMO modellers would be nice to know how these various simulations perform in other sea-ice covered seas, not just in the Arctic. Also, is there a specific configuration that is recommended and differs from default settings?
The text could be improved significantly. Its expressions are often confusing and the writing style could be clearer, in particular scientific terms and simulation names could be used in a more uniform manner, which would help reading a lot. I give some detailed comments to assist on this:
- line 58: Instead of writing 'several regions', be specific and write '3 regions'
- line 58: 'by varying the mixing scheme' is more clear, if that is meant here. 'varying the mixing just below the mixed layer TKE
parameters' is hard to understand. - line 61: Change to 'The vertical turbulent kinetic energy (TKE) closure scheme'
- line 63: Change to 'integrated by Blanke and Delecluse (1993) into the OPA model'
- line 75: 'ocean's surface boundary layer'
- line 84: Can f_r be > 0.1? If it can, would it make sense to try higher values?
- line 84: Do you mean that 5% of TKE is redistributed below MLD or in the MLD and below it?
- line 94: Explain what ORCA1 is.
- line 97: 'Rathore et al. (2024)' [year is missing]
- Table 1 and others: Table captions are usually above tables.
- line 114: Is pss Practical Salinity Scale? Does it equal psu which is more commonly used?
- line 120: 'de Boyer Montegut (2024)' [typos in the reference]
- line 125: 0.03 kg/m$^3$ [there should be space between number and unit, change also line 127]
- line 127: '(e.g. de Boyer Montegut et al. 2004)' [remove double parenthesis, change also line 130]
- Table 2: It is difficult to follow in text what simulation is what. Perhaps number simulations in this table from 1 to 9 and use
those numbers consistently in text when refererring to particular simulations. - line 137: 'We computed'
- line 140: 'The ITP data are from'
- Figure 2 caption: 'Basins following Peralta-Ferriz and Woodgate (2015).'
- lines 155-157: Seems that latitude-longitude boxes can not be used to define the regions drawn in Figure 2.
- lines 157-162: This text does not belong here but fits better in Introduction.
- line 172: 'deeper than modelled ML'
- lines 176-178: It is hard to see these salinity differences from Figure 4 which shows absolute salinities. A solution is to add
Control-LOPS salinity difference panels to Figure 4. - line 181: 'data are averaged spatially and temporally'
- line 189: Here is a place where indicating the parameter would be a reader: 'Increasing the fraction of surface TKE (rn_efr) that'
- line 190: Do you mean 0.75 rather than 0.08?
- lines 191-192: The effect of nn_tau is not shown in Figure 5, so this sentence in this context is a bit misleading.
- line 194: 'could arguably be caused'
- line 195: 'However, we recomputed'
- line 198: Is the MLD std modelled and/or regional, what are the summer months used? This information should be added.
- line 201: 'and is only 8 m'
- line 204: 'with MDL at least 20 m'
- line 207: '20 m' [add space between number and unit]
- line 207: 'in the Canada Basin and the Eurasian Basin in March'
- line 209: The spatial distribution looks similar in September too, not only in March.
- lines 215-216: The effect of deeper ML on surface salinity is mentioned. What about temperature?
- line 226: Could continue that the TKE MLP parameters do not improve the temperature maximum below 200 m because the temperature maximum is due to the heat advection at that depth.
- line 266: Note that there is also negative sea-ice thickenss bias in March.
- line 266: Is the East coast the Siberian coast?
- line 267: 'negative biases between'
- line 268: Discrepancies between OSI-SAF and model are not relatively minor in September but look quite substantial. What is the OSI-SAF uncertainty? Would make sense to use other satellite sea-ice concentration data for comparison. It is known that different data sets provide somewhat different sea-ice concentration and the ice edge.
- lines 295-296: The link between stratification and sea-ice melt is interesting. What would be the role of ice-albedo feedback strengthening this process?
- line 299: Would be reader friendly to remind what these three sensitivity experiments are, not just list their parameter values. Here Tables 1 and 2 could be used/ cited.
- Figure 11. Differences seem often not very clear and there can be biases both ways. Perhaps calculating RMS or mean absolute differences per basin would give clearer and more quantitative measure.
- Figure 12: Caption: 'MLD, sea-ice concentration and sea-ice thickness'. How much these timeseries correlate? Seems that year-to-ear variability is captured rather well due to the same atmospheric forcing and free drift in summer.
- line 317: 'is close to'
- line 325: What is short-term?
- line 332: 'exhibit closer resemblances'
- line 346: Can you link the MLD increase to the erosion of the Arctic halocline?
- line 370: 'in these regions than the 1-f_i option'
- line 371: 'the LOPS climatology'
- line 372: 'upper ocean vertical properties'
- line 398: 'their coverage'
- line 401: 'to entirely remove TKE MLP'
- line 405: 'this parameterization lacks a'
- lines 406-407: Do you have a reference for the claim of time dependence of the NEMO TKE MLP scheme?
Citation: https://doi.org/10.5194/gmd-2024-49-RC1 -
RC2: 'Comment on gmd-2024-49', Anonymous Referee #2, 05 Jun 2024
reply
The comment was uploaded in the form of a supplement: https://gmd.copernicus.org/preprints/gmd-2024-49/gmd-2024-49-RC2-supplement.pdf
Data sets
Sensitivity experiments of the parameters involved in the turbulent kinetic energy mixed layer penetration scheme of the NEMO ocean model S. Allende https://doi.org/10.14428/DVN/NZSKTU
Model code and software
nemo: v4.2.1 S. Allende https://doi.org/10.5281/zenodo.10732752
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
445 | 84 | 27 | 556 | 24 | 22 |
- HTML: 445
- PDF: 84
- XML: 27
- Total: 556
- BibTeX: 24
- EndNote: 22
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1