Does the wind stress always damp an oceanic eddy?

Up to now, the literature has shown that the relative wind stress does negative work on ocean mesoscale eddies. In other words, the relative wind stress inhibits the development of the eddies. However, based on a newly derived simplified theoretical model, the present study finds that under the action of a steady and uniform wind field, eddies can rapidly obtain kinetic energy from the wind field following several hours of adaption and adjustment, in which the wind stress transitions from doing negative to positive work. The finding is supported by the fact that the relative wind stress work on oceanic eddies over the northeastern tropical Pacific ocean is positive with the nearly constant gap wind. This implies that energy input from the wind is sensitive to eddy velocity structure, and hence, wind stress is not always a killer of eddies.


Introduction
Sea surface winds are a major energy source for the oceanic motions. How the wind energy is transferred into the ocean is an important topic in the climatic studies. Since mesoscale oceanic eddies are ubiquitous phenomena in world oceans (Dong et al. 2014), understanding how oceanic eddies affects the wind energy input into the ocean has received significant attentions in oceanographic community, e.g. Zhai et al (2012), Xu et al (2016), Renault et al. (2016Renault et al. ( , 2018. The power that wind energy inputs into the ocean is estimated by the product of the wind stress and the sea surface velocity (Faller 1968;Oort et al. 1994;Wunsch 1998) and is as follows: where P is the power per unit area of work done on the ocean by the wind, ⇀ τ is the sea surface wind stress, ⇀ U is the sea surface velocity. The sea surface wind stress is usually calculated in the form of the second power of relative wind speed (Qin 1980;Dewar and Flierl 1987;Pacanowski 1987;Scott and Xu 2009): where ρ a is the air density above the sea surface ( 1.25kg · m −3 ), ⇀ U a indicates the atmospheric wind speed at a height of 10 m above the sea surface, C d represents the drag coefficient of sea surface, which performs the momentum transfer efficiency between atmosphere and ocean, and is dependent on wind speed, currents, waves, sea surface roughness, and stability amongst other factors (Charnock 1955;Garratt 1977;Wu 1980;Powell et al. 2003;Donelan et al. 2004;Makin 2005;Kara et al. 2007;Zhao et al. 2015). As ocean current velocities are usually less than 1m · s −1 , which is much smaller than the sea surface wind speed, the surface wind stress exerted by the atmosphere on the ocean can be approximately expressed as a function of the atmospheric wind (Bunker 1976;Bye 1985): (2) (3)

Open Access
*Correspondence: cmdong@nuist.edu.cn 1 Oceanic Modeling and Observation Laboratory, School of Marine Sciences, Nanjing University of Information Science and Technology, Nanjing, China Full list of author information is available at the end of the article The primary difference between Eqs. (3) and (2) is that Eq. (3) ignores the contribution of sea surface velocity and only considers the wind speed. In this study, Eq. (3) is termed the absolute wind stress, while Eq. (2) considering the sea surface current is termed the relative wind stress.
It has been suggested by the literature that the work of the relative wind stress on the ocean could be decreased significantly when oceanic mesoscale eddies are presented in oceanic models with the relative wind stress used as suggested by Duhaut and Straub (2006), Dawe and Thompson (2006), Zhai and Greatbatch (2007), Scott andXu (2009), Zhai et al (2012), Xu et al (2016), Renault et al. (2016Renault et al. ( , 2018, Yu et al (2018). These studies identify that if the relative wind stress is used rather than the absolute wind stress, a decrease of 17-35% in the input of the energy to the ocean could be observed and this is linked to the influence of the ubiquitous mesoscale ocean eddies. Particularly, Hughes and Wilson (2008), Zhai et al. (2012), Xu et al (2016), and Renault et al (2018) find that relative wind stress can damp the development of ocean mesoscale eddies. Xu et al. (2016) used satellite observational data to quantify the work of the global wind stress on mesoscale eddies, and their study shows that the work of relative wind stress on mesoscale eddies is negative in most oceanic regions. Considering the role of wind stress in damping the eddy intensity and activity, Renault et al. (2018) term the wind stress to be an 'eddy killer' .
However, so far, these studies have only considered the transient response of the eddies to the relative wind stress and ignored the dynamic response of ocean mesoscale eddies to wind stress over time. This study hypothesizes that the primary reason why the wind stress does negative work on eddies is that the direction of the wind stress is inconsistent with the direction of an eddy's velocity. However, once the eddy gradually evolves in response to the wind stress, the wind would be able to perform positive work on the eddy through an injection of energy. Consequently, relative wind stress, will no longer be an 'eddy killer' . To test this hypothesis, this study considers an idealized eddy under the influence of a continuous and steady wind field to calculate the power of the relative wind stress on an oceanic eddy. Through the time integration of the wind work on the eddy, the mechanism behind the adaptive adjustment of the eddy flow field due to wind energy input can be revealed.
The rest of this paper is structured as follows: "Assumptions" section describes the assumptions. "Response of an eddy to wind stress" section presents the main results.
"Discussion" section is the discussion, and "Conclusions" section provides the conclusions.

Assumptions
Suppose the radius ( R e0 ) of an ideal cyclonic eddy is 30 km, where r x, y r = x 2 + y 2 ≤ R e0 represents the distance from the center of the eddy, the velocity at the center of the eddy is 0, and the maximum value is located at the edge of the eddy, where the speed is U e0 = 1m · s −1 , the eddy velocity ⇀ U increases linearly with r , and the velocity component at any position the unit is m · s −1 , it tangentially rotates counterclockwise, and the velocity distribution can be given as shown in Fig. 1a.
Since the scale of the spatial variation of the atmospheric wind field is much larger than that of the ocean mesoscale eddy, it can be assumed that the background wind field is approximately uniformly distributed over the eddy, and that the background wind field is stationary, the intensity of wind field ⇀ U a is U a0 which blows from the north to the south, the magnitude of the wind is 10 m · s −1 , the components of the wind velocity can be written as u a = 0 , v a = −U a0 , and the unit is m · s −1 . Wind stress is calculated by the absolute wind stress (Eq. 3) and the relative wind stress (Eq. 2).
. Red arrows and blue arrows indicate eddy velocity using the equations of absolute wind stress and relative wind stress, respectively (Unit: m · s −1 ).

Response of an eddy to wind stress
The instantaneous work by the absolute (relative) wind stress on the above idealized eddy is plotted in Fig. 2a (2e), and it is suggested that the relative wind stress inhibits the total eddy kinetic energy growth, while the absolute wind stress does not change the energy. Here, although the eddy flow field and shape may change, these are not considered under the long-term act by wind stress. These changes may, however, also affect the magnitude and direction of wind stress, as well as the total power received from wind stress on the eddy. In this section, a steady, uniform wind is applied to the idealized eddy over a long period of time, and changes in the flow field are studied.
Considering that the water body is only affected by wind stress, it is assumed that the friction stress at the bottom of Ekman layer can be ignored, and the thickness of the Ekman layer (the deepest depth the wind stress can affect) is D e . The vertical average velocity of water body is To simplify the problem of obtaining the mathematical relationship between the surface velocity and the vertical-averaged velocity in the upper layer, the vertical profile of the velocity is assumed to be linear other than exponential form as Ekman drift in the present study. By only considering the energy input by the wind stress (assuming that the Coriolis effect and pressure gradient balance do not enter the dynamic balance at the first adjustment stage), the velocity change due to the wind stress can be expressed as Integrating over time, we have The spatial distribution of eddy velocity at different time steps can be obtained according to Eq. (5), as shown in Fig. 1. As can be seen from this figure, the speed on the left of the eddy increases and the speed on the right decreases. The eddy core center also gradually deviates to the right forced by both the relative and absolute wind stresses, and the velocity magnitudes forced by the absolute wind stress are slightly larger than those by the relative wind stress.
Under this condition U a0 = 10m · s −1 , U e0 = 1m · s −1 , using Eqs. (1), (2), and (5), the powers of the wind stress work at different times can be obtained. As shown in Fig. 2a, forced by the absolute wind stress, the area of negative work is equal to that of positive work at the beginning, while the area of negative work is larger than that of positive work by the relative wind stress. This result is consistent with the reports in the literature (e.g., Anderson et al., 2011;Renault et al., 2018;Flexas et al., 2019;Rai et al., 2021), in which changes in the eddy flow field and shape are not considered under the long-term act by wind stress. However, these changes may also affect the magnitude and direction of relative wind stress, as well as the total power received from the wind stress on the eddy. In this section, a steady, uniform wind is applied to the idealized eddy over a certain time, and changes in the flow field are studied. After 4 h, the area of positive power gradually increases forced by both the relative and absolute wind stresses, so that the wind stress on the whole eddy area does positive work on the eddy. Figures 2b-d gives results for the absolute wind stress, and Fig. 2f-h gives results for the relative wind stress. Shown in Fig. 2, due to considering the eddy velocity in the case with the relative wind stress, the zero line of the power is a curved line for the case of the relative wind stress, while it is a straight line for the case of the absolute wind stress. The temporal changes of the powers integrated in the initial eddy area forced by the relative and absolute wind stresses are plotted in Fig. 3a. Both the powers increase with time but the power forced by the absolute wind stress is zero or positive, while the power is negative before the time 1.58 h and then becomes positive for the case of the relative wind stress. The main reason for the increase of the positive wind work area is that the eddy velocity has changed in response to the wind stress and gradually follows the wind stress direction. The eddy velocity field adapts to the overlying wind field. Then, and only then, will wind stress do the positive work on the eddy.
As resulted from the power change with the time, Fig. 3b shows the changes of the total kinetic energy over the eddy area with the time. The KE increases with time forced by the absolute wind stress due to the power is always positive, while the KE decreases with time within the first 1.58 h and then increase with time forced by the relative wind stress corresponding to the power change with time in Fig. 3a. Figure 3a and b also demonstrate that the wind energy input by the relative wind stress is significantly smaller than that of the absolute wind stress.

Discussion
To examine the sensitivity of the inflection point to the wind intensity from the negative wind work to the positive wind work, the eddy KE under different wind speeds are calculated. Under different wind speed conditions (8, 10, and 12 m · s −1 ), the change of the total kinetic energy of the eddy with time within the initial eddy range is calculated, as shown in Fig. 4. Within 12 h, the total kinetic energy of the eddy decreases at first and then increases. The inflection points of the eddy kinetic energy is at 3.07, 1.58, and 0.92 h with wind speeds of 8, 10, 12 m · s −1 , respectively. It can be concluded that under the above assumption, the work of relative wind stress on the eddy is negative work first and then positive work.
From the above discussion, it can be speculated that when an eddy is generated by the wind forcing (i.e., wind curl), the eddy current flow will be aligned with the wind direction, and then the wind stress will do the positive work. The hypothesis is supported by Xu et al. (2016), where the authors quantifies the work of global relative wind stress on mesoscale eddies using satellite observations. They identified that work is negative in most areas but is surprisingly positive in the northeastern tropical Pacific, see Fig. 1 in Xu et al. (2016). The area is constantly forced by the wind jets over the Gulfs of Tehuantepec and Papagayo and the low-frequency wind force is suggested as the primary force for the generation of the eddies in the northeastern tropical Pacific (Liang et al. 2012). According to the theoretical analysis outlined in this paper, under the action of a long-term wind field, eddy structures have adjusted and adapted to the wind field, so that wind energy can input kinetic energy into the eddy, i.e., the wind stress does the positive work on the eddies.

Conclusions
Previous studies so far have suggested that the relative wind stress acts as a "killer" to oceanic eddies by suppressing the development of oceanic eddies. The  In this study, to explore the mechanism of the variation in wind stress work on an oceanic eddy, under an assumption that the wind stress plays a major role and other factors are ignored, we design a simplified ideal model by only considering the momentum transfer effect of atmospheric wind field on an ideal oceanic eddy. The effects of ocean mesoscale eddies on wind stress energy input are studied using two different equations of absolute and relative wind stress. It is found that at the beginning, the total power of wind stress integrated over the eddy area is 0 when the eddy velocity is not included and the total power of wind stress is negative when the eddy velocity is included. However, considering the response and adjustment of the eddy under the action of wind stress, the present research shows that the work of absolute wind stress is always positive and the work of relative wind stress is negative in the first few hours, after the inflection point, the work is positive. It implicates that when the eddy velocity field is adaptive to the wind field, the wind stress does the positive work on the eddy otherwise it will do negative work. In other words, if the eddy is generated by wind stress, such as the eddies in the southwest of Mexico in the Pacific Ocean, where the wind stress does positive work on eddies. Therefore, it can be concluded that the wind stress does not always do the negative work on an eddy. When an eddy structure is independent on the wind field, the wind stress does negative work integrated over the eddy area but when the eddy velocity field adjusts to the wind field, the wind stress does the positive work. This finding implicates that the energy input from the wind is sensitive to the mesoscale eddy velocity structure.
It should be noted that the KE of eddy cannot grow infinitely. In reality, the eddy development involves many physical processes: wind energy input, eddy and mean flow interaction, turbulent dissipation, et al. After the early stage of the adjustment, all the processes mentioned above will take effect and the eddy will finally reach an equilibrium state. The dynamic balance determines the time scale to the equilibrium. The investigation on the eddy equilibrium is beyond the scope of the present manuscript and will be left for the future study.
The present study, together with the previous relevent research, suggest that the signs (positive or negative) of relative wind stress work on oceanic eddies depend on the duration and variation frequency of the sea surface wind. The further exploration of the topic requires more extensive investigations including numerical modeling and in situ observations, which will be left for future studies.