the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Satellite-based Analysis of Ocean-Surface Stress across the Ice-free and Ice-covered Polar Oceans
Abstract. Ocean-surface stress is a critical driver of polar sea ice dynamics, air-sea interactions, and ocean circulation. This work provides a daily analysis of ocean-surface stress on 25-km Equal-Area Scalable Earth (EASE) Grids across the ice-free and ice-covered regions of the polar oceans (2011–2018 for Arctic, 2013–2018 for Antarctic), covering latitudes north of 60° N in the Arctic and south of 50° S in the Antarctic and Southern Ocean. Ocean-surface stress is calculated using a bulk parameterization approach that combines ocean-surface winds, ice motion vectors, and sea surface height (SSH) data from multiple satellite platforms. The analysis captures significant spatial and temporal variability in ocean-surface wind stress and the resultant wind-driven Ekman transport, while providing enhanced spatiotemporal resolution. Two sensitivity analyses are conducted to address key sources of uncertainty. The first addresses the fine-scale variability in SSH fields, which was mitigated using a 150-km Gaussian filter to smooth three-day SSH datasets and enhance compatibility with the other monthly product, followed by linear interpolation to achieve daily resolution. The second investigates uncertainty in the ice-water drag coefficient, which revealed that variations in the coefficient have a proportional influence on the computed ocean-surface stress under the tested conditions. These uncertainties are most pronounced during winter, with median values reaching 20 % in the Arctic and 40 % in the Southern Ocean. Validation efforts utilized Ice-Tethered Profiler velocity records, revealing moderate correlations (r = 0.6–0.8) at monthly timescales, effectively capturing low-frequency signals but with small northward biases. Satellite-derived velocity fields, including both Ekman and geostrophic components, explain 40–50 % of the total variance. The unexplained variance reflects unresolved processes, such as mesoscale dynamics and other unparameterized factors. This dataset is publicly available at https://doi.org/10.5281/zenodo.14750492 (Liu & Yu, 2024).
- Preprint
(11789 KB) - Metadata XML
- BibTeX
- EndNote
Status: final response (author comments only)
-
RC1: 'Comment on essd-2025-14', Anonymous Referee #1, 19 Apr 2025
Revised Review of Satellite-based Analysis of Ocean-Surface Stress across the Ice-free and Ice-covered Polar Oceans by Liu and Yu (2025)
The authors present an analysis of satellite-derived ocean-surface stress across the ice-free and ice-covered polar oceans. They utilize a bulk parameterization approach, combining multiple satellite datasets, to estimate ocean-surface stress and Ekman transport. The study provides valuable insights into the spatial and temporal variability of these key drivers of polar climate dynamics. However, I have several concerns regarding the methodology, data processing, and presentation, which must be adequately addressed before the manuscript is suitable for publication.
General Comments:
- The authors need to update their references to include some relevant papers especially when it relates to relevant dataset (CPOM extended dataset was recently used in Lin et al (2023)), parameterisations (Brenner et al 2021, Lupkes et al 2015), modelling work that can server as motivation (Tsamados et al 2014, Sterlin et al 2023).
- The validation with ITP observations is a valuable component, but its limitations need to be discussed in more detail. In particular more discussion is needed and this section needs rewriting to better account for the impact of time and space sampling in the input datasets on the derived ocean surface stress.
- The authors perform a sensitivity analysis to value of drag and gaussian spatial filters but do not discuss errors due to model parameterisation (i.e. impact of ue vs ug vs nothing), impact of uncertainty of reanalysis, impact of ocean (and atmospheric) stratification (see papers by Lupkes and Brenner for example), impact of time and length scale over which stresses are calculated.
- The authors could consider extending the analysis period to include more recent available data (CPOM or CLS) and by removing their current constraint due to the reanalysis they use. Please justify why this reanalysis is really so good to justify loosing 5+ years of analysis.
- The authors use open access satellite products but fail to share their code! Share code / Github
References:
Lin, Peigen, Robert S. Pickart, Harry Heorton, Michel Tsamados, Motoyo Itoh, and Takashi Kikuchi. "Recent state transition of the Arctic Ocean’s Beaufort Gyre." Nature Geoscience 16, no. 6 (2023): 485-491.
Brenner, Samuel, Luc Rainville, Jim Thomson, Sylvia Cole, and Craig Lee. "Comparing observations and parameterizations of ice‐ocean drag through an annual cycle across the Beaufort Sea." Journal of Geophysical Research: Oceans 126, no. 4 (2021): e2020JC016977.
Lüpkes, Christof, and Vladimir M. Gryanik. "A stability‐dependent parametrization of transfer coefficients for momentum and heat over polar sea ice to be used in climate models." Journal of Geophysical Research: Atmospheres 120, no. 2 (2015): 552-581.
Tsamados, Michel, Daniel L. Feltham, David Schroeder, Daniela Flocco, Sinead L. Farrell, Nathan Kurtz, Seymour W. Laxon, and Sheldon Bacon. "Impact of variable atmospheric and oceanic form drag on simulations of Arctic sea ice." Journal of Physical Oceanography 44, no. 5 (2014): 1329-1353.
Sterlin, Jean, Michel Tsamados, Thierry Fichefet, François Massonnet, and Gaia Barbic. "Effects of sea ice form drag on the polar oceans in the NEMO-LIM3 global ocean–sea ice model." Ocean Modelling 184 (2023): 102227.
Specific Comments:
- L6 EASE2? clarify if so
- L15 clarify why you need both CPOM and CLS?
- L19 clarify / correlation of what
- L21 not showed in this paper?
- L40 ice-ocean governor is Menenghello citation?
- L47 add references here i.e Lin et al 2022 but also more references in introduction for motivation / modelling work
- L51 Seems a waste to limit yourselves to 2018 due to reanalysis
- L77 cite recent papers (including modelling)
- L92 discuss impact on ocean surface stress of time sampling i.e. daily vs 3d vs monthly
- L92 discuss why this reanalysis so important
- L117 add references to paper Lin that used more extended CPOM extension reference
- L126 is it worth limiting to that reanalisis then?
- Figure 2 why gaussian filtering to reconcile with CPOM?
- L144 why not just use CLS if better? You criticise CPOM product but you fail to discuss all issues in these CLS and CPOM product when it comes to continuity of SSH at the ice-ocean boundary.
- L186 clarify
- L190 not so clear on the map. Shouldn’t the maps be monthly?
- Figure 4 mean is noisy so what about monthly / daily products. Is that mean over all winter and all months?
- Figure 5 what is the total region used? Not clear what is aaw
- L277 cite recent papers and modelling ones
- L312 what about error from time sampling, model physics parameterization
- Figure 9 is based on two errors only not sure I buy this. Also not clear how the range of tau_w chose 10^-3 to 10^-2.
- L340 test impact on derived stress
- L366 shouldn't you compare at the time averageds used in stress calculations? i.e. daily
- Figure 11 is u really daily?
- Figure 12 representation issue here clearly seen (space or time or both?)
- Figure 14 obvious representativity issue (space, time or both?). I would redo all this section on monthly products and discuss better
- Figure 15 what about plot for daily timescales
- L433 share code!
Citation: https://doi.org/10.5194/essd-2025-14-RC1 -
AC1: 'Reply on RC1', Chao Liu, 17 May 2025
We thank the reviewer for their time and for providing positive feedback and constructive comments on our manuscript. We have carefully revised the manuscript accordingly and incorporated all the proposed changes. The attached PDF file provides a detailed response addressing each of the reviewer’s concerns.
-
RC2: 'Comment on essd-2025-14', Anonymous Referee #2, 21 Apr 2025
This is a well written manuscript and the analysis seems legitimate and thoughtful. The first thing that came to me after reading the manuscript is why the new dataset of ocean surface stress is needed? Since there are many assumptions have been made in order to calculate the ocean-surface stress. It will be great if the authors can explain the motivation and reasoning that this data is needed for the community. Second, the value of Ekman depth used in the calculation is not explained, eqn. 4 & 5. The Ekman depth is needed in order to calculate the vertical Ekman velocity. What are the values of Ekman depth that you used? Third, for the scatterplots of ITP vs satellite (Fig. 12 & 13), how the ITP data is processed to be comparable to the daily satellite product? Did you do time average on ITP data to daily data or interpolate the satellite daily product to the ITP profile time? From the number of samples, it seems the author interpolated the daily satellite data to the ITP data (profile) time. If so, the scatter plots are meaningless. The comparison is wrong, since the satellite temporal resolution is under sampled. You can not draw conclusions from these scatterplots. Also, the 30-day low pass resulted in too few data points (Fig. 14, 15 & Table 3). I suggest to do a time-average ITP data to daily or weekly and , and redo the analysis.
Other small comments
Line 40: "ice-ocean governor," → "ice-ocean governor",
Line 60: Figure 1(a) caption: there is no dashed magenta has shown. Is the September ice extent boundaries out of the scope of the plot? The Solid line near the Novaya Zemlya has a portion of dashed line. Is this supposed to be solid?
Line 84 It is unclear to me how the De: Ekman layer depth is determined in the calculation. Can you explain?
Acronym stands for?
Line 98 acronym OAFLux2?
Line 106 acronym SSM/I, AMSR-E, AVHRR, MODIS, QuikSCAT?
Line 107 acronym IABP?
Line 111 acronym CLS/PML?
Line 120 acronym SMMR DMSP, NSIDC?
Line 139 "Noting the 25 km resolution could introduce uncertainties near the 15% sea ice concentration boundary." --> How did you arrive this conclusion?
Line 152 "where W represents the local vertical Ekman velocity W_e" → W and W_e were used interchangeably. Can you decide which one to use and keep it consistently?
Line 227 "in the Indian Ocean sector and the southeast Pacific" → Can you label it on Fig 6b?
Line 232 of 50°S (Figure 6d), → Can you label the latitude and longitude on Fig 6?
Line 293 Baffin Bay, the Chukchi Sea, and north of Fram Strait → Can you label these places in Figure 8?.
Line 309 What do you mean “ full eight-year/six-year period” period?
Lien 324: the Mendeleev Ridge. → please label it on the Figures
Line 327-328 Weddell and Ross Seas, Antarctic Peninsula and western Ross Sea, Enderby Land and the Amundsen Sea → please label these places on the Figures.
Line 381 What do you mean “paired” observations?
Line 385 “ITP-80 pair shows better agreement in the zonal component after 200 days.” → I can't tell if this is true from the figure 12. How did you arrive this conclusion?
Citation: https://doi.org/10.5194/essd-2025-14-RC2 -
AC2: 'Reply on RC2', Chao Liu, 17 May 2025
We thank the reviewer for their time and for providing positive feedback and constructive comments on our manuscript. We have carefully revised the manuscript accordingly and incorporated all the proposed changes. The attached PDF file provides a detailed response addressing each of the reviewer’s concerns.
-
AC2: 'Reply on RC2', Chao Liu, 17 May 2025
Data sets
Arctic/Antarctic Ocean-Surface Stress Analysis, 2011-2018/2013-2018 Chao Liu and Lisan Yu https://doi.org/10.5281/zenodo.14750492
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
257 | 32 | 14 | 303 | 14 | 13 |
- HTML: 257
- PDF: 32
- XML: 14
- Total: 303
- BibTeX: 14
- EndNote: 13
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1