News Column

Surface Roughness Imaging of Currents Shows Divergence and Strain in the Wind Direction

August 1, 2014

Klein, Patrice


Images of sea surface roughness-for example, obtained by synthetic aperture radars (SAR) or by radiometers viewing areas in and around the sun glitter-at times provide clear observations of meso- and submesoscale oceanic features. Interacting with the surface wind waves, particular deformation properties of surface currents are responsible for those manifestations. Ignoring other sources of surface roughness variations, the authors limit their discussion to the mean square slope (mss) variability. This study confirms that vortical currents and currents with shear in the wind direction shall not be expressed in surface roughness images. Only divergent currents or currents with no divergence but strained in the wind direction can exhibit surface roughness signatures. More specifically, nondivergent currents might be traced with a 458 sensitivity to the wind direction. A simple method is proposed in order to interpret high-resolution roughness images, where roughness variations are proportional to ?u/?x + a?v/?y, a linear combination of the along-wind and crosswind current gradients. The polarization parameter a depends upon the sensor look direction and the directional properties of the surface waves selected by the sensor. The use of multiple look directions or possible acquisitions with different wind directions shall thus help to retrieve surface currents from surface roughness observations.

(ProQuest: ... denotes formulae omitted.)

1. Introduction

Images of sea surface roughness are routinely obtained using active or passive sensors-such as microwave radars including synthetic aperture radar (SAR) (e.g., Beal et al. 1981) or optical radiometers viewing areas in and around the sun glitter (e.g., Apel et al. 1975)-that operate at different wavelengths, polarizations, and viewing angles. Under low to moderate wind conditions, those images provide clear observations of a wide range of oceanic phenomena at scales of about 1 to 30 km, including internal waves, mesoscale, and submesoscale features such as fronts, filaments, and spiraling eddies (e.g., Fu and Holt 1983; Alpers 1985). Figure 1 shows an example of such an image. Because surface currents can strongly affect small-scale (capillaries and short gravity) waves under moderate wind conditions, it is possible to track and quantify current gradients (Kudryavtsev et al. 2005; Johannessen et al. 2005).

Spurred partly by improvements in computing power, numerical models have recently confirmed the impor- tance of ocean features at the submesoscale for ex- changes within the ocean and between the ocean and the atmosphere (e.g., Fox-Kemper et al. 2011; LÉvy et al. 2012; Ferrari 2011). There is however a lack of system- atic observations at these scales. Efforts to remedy this situation are under way, as illustrated by the upcoming Surface Water and Ocean Topography (SWOT) satel- lite mission (Fu et al. 2012), but at the moment obser- vations at these scales come mainly from images of infrared-based sea surface temperature and ocean color reflectance. The use of surface roughness images has been largely overlooked in that context, whereas with large footprints and high-resolution capabilities, surface roughness images can already be valuable to further improve qualitative and quantitative understandings of upper-ocean submesoscale processes.

Surface roughness variations in sun glitter images are controlled by specularly reflected sunlight (Cox and Munk 1954) and related to changes of wave slope dis- tribution. The latter mostly correspond to variations of mean square slope (mss) (Kudryavtsev et al. 2012b). For radar images, the interpretation is more complex as ra- dar backscatter intensity is affected by the mss but also by scattering effects and geometrical properties of the surface (Kudryavtsev et al. 2005). Details on how to estimate mss from glitter brightness or from radar backscatter are beyond the scope of the present paper, and the reader is referred to Kudryavtsev et al. (2012b). For simplicity in this paper, we will assume an un- equivocal relationship between surface roughness and mss, which is an adequate assumption at least for optical observations.

Local variations of surface roughness have been at- tributed to three main processes: direct wave-current interactions (e.g., Johannessen et al. 1996), attenuation by surfactants accumulated in current convergence zones (e.g., Espedal et al. 1998), and varying wind field resulting from transformations of the atmospheric boundary layer across sea surface temperature fronts (e.g., Beal et al. 1997). Surface slicks have specific roughness signatures and mainly occur in coastal waters under light wind conditions (e.g., McWilliams et al. 2009). In this paper, we do not consider the effects of surfactants. Atmospheric boundary layer adjustments depend on the atmosphere stability induced by changes in the sea surface temperature.1 Variations of atmo- spheric stability seemingly occur at scales larger than those shown in Fig. 1 (see, e.g., Kudryavtsev et al. 2005, their section 4), therefore we do not consider these ef- fects either.

Thus, not considering the processes associated with surfactants and atmospheric stability, the primary pro- cess that helps imaging the underlying flow dynamics is governed by interactions between surface wind waves and surface currents. Accordingly, questions arise. What information on ocean turbulence at meso- and sub- mesoscale is contained in surface roughness variations? Surface current gradients can be described with four types of deformations: divergence, vorticity, strain, and shear2 (Fig. 2, sketches 1 to 4). Which type of current deformation is mostly visible in roughness images? Can we develop a method that inverts surface roughness into current information?

The surface roughness images that are most easily interpreted are those depicting currents induced by bottom topography or by internal waves (e.g., Alpers 1985), where clear patterns of convergence/divergence of surface currents are detected as bright/dark contrasts. More generally, an individual surface wave traveling over a varying current might be deflected by the gradient of the cross-propagation current (sometimes inadequately called current shear, Fig. 2, sketch 5, if we suppose that the wave travels in the x direction), but the gradient of the along-propagation current (sometimes inade- quately called divergence, Fig. 2, sketch 6) can more strongly alter the steepness properties (e.g., Phillips 1984). As such, remote sensing experts often refer the current deformation as current divergence, implicitly assuming that it refers to the gradient of the along- propagation current (Phillips 1984; Alpers and Hennings 1984).

Based on the conservation of action (e.g., Phillips 1984), modulation not only of an individual wave but of a full spectrum of wind waves can be determined. Such spectrum with waves propagating in different directions is modulated by different current deformation types, including current divergence but also current shears. Dulov and Kudryavtsev (1990) demonstrated the latter effect by considering fluctuations of wave breaking in- tensity. Yet, subsequent works reported the effect of current shear to be negligible compared to that of cur- rent divergence (e.g., Jansen et al. 1998). As such, it has become customary to assume that divergence was the dominant deformation type (e.g., Johannessen et al. 2005; Kudryavtsev et al. 2005, 2012a). This is indeed a good approximation for currents with significant di- vergence, as those induced by bottom topography (Alpers and Hennings 1984), by internal waves, or along sharp density fronts (Jansen et al. 1998). On the other hand, for meso- and submesoscale turbulence, currents are constrained by Earth rotation and consequently exhibit weak divergence and larger shears. Thus, it might be somehow misleading to immediately assume that current divergence is into play whenever effects of current deformation are explored.

In this paper, we wish to clarify the interpretation of surface roughness variations in terms of current defor- mation. In section 2, we recall the modulation of waves by currents, extend it from a single wave to a spectrum symmetrical about the wind direction, and calculate resulting mss variations. Out of the four types of sur- face current deformations, two can exhibit a roughness signature, the divergence and the strain in the wind direction, whereas the vorticity and the shear in the wind direction have no signature. The wind direction and the sensor look direction are irrelevant for di- vergence detection. This is not the case for strain de- tection. In section 3, we describe the roughness responses to divergence and strain for different spec- tral angular spreads of wind waves. This provides a simple method to interpret roughness images, de- pending on the wave directional spread and the sensor look direction. In section 4, results are illustrated using a forward numerical simulation including mss varia- tions. Conclusions and perspectives are finally drawn in section 5.

2. Surface roughness and current deformation types

a. Surface waves in currents

In the presence of a current, the conservation of wave action N(x, k, t) reads as (e.g., Phillips 1984)

... (1)

where cg is the group velocity, u is the current, x is the horizontal position, k is the wavenumber, v is the in- trinsic frequency, and S is the energy source that con- tains input from the wind, dissipation, and wave-wave nonlinear interactions. Repeated indices i, j 5 1 to 2 in- dicate summation over horizontal components.

Following a relaxation approach (e.g., Keller and Wright 1975; Hughes 1978; Alpers and Hennings 1984), the action is written N (x, k, t) 5 N0(k) 1 N~ (x, k, t), where it is supposed to experience a small disturbance N~ with respect to a background value N0 corresponding to the state undisturbed by currents. Sources rapidly re- store equilibrium and are represented by a linear re- laxation S/v 52N~ /tc ,wherethetimescaletc(k)depends on the wavenumber. Assuming that the typical length scale L of the current is larger than the relaxation length scale lr 5 tcjcgj, the advection term on the left-hand side of (1) can be ignored. In a steady state, the action anomaly N~ due to current variations then is

... (2)

And in polar coordinates

... (3)

where k and f are the wavenumber magnitude and di- rection, where mk 5 ? lnN0/? lnk and mf 5 ? lnN0/?f, and where notations x, y and u, y are used in place of x1, x2 and u1, u2.

b. Surface roughness changes

As stated in the introduction, we suppose that the detected variations of surface roughness directly relate to local mss variations. From the wave action anomaly N~ , the mss anomalies are

... (4)

where mssx and mssy are the x and y projection of the total mss anomaly msst .

The wavenumber interval where the integration has to be performed is sensor dependent. For instance, at high zenith incidence angles, SARs in the C band are rather sensitive to scattering by waves of the order of 0.1-m wavelength. On the contrary, sun glitters in- volve specular reflections by waves from the whole spectrum.

Additionally, a sensor might not be sensitive to the total slope distribution msst but rather to the distribution of slopes in the azimuthal look direction of the sensor (e.g., mssx if the sensor looks in the x direction). That sensitivity is not well established yet fsee Kudryavtsev et al. [2005, their (10)] for the SAR and see Kudryavtsev et al. [2012b, their (6)] for the sun glitterg. Awaiting further improvements, one might consider that a sensor detects a weighted average of the total mss and the mss projection in the look direction.

c. Current tensor decomposition

To ease the analysis, the deformation tensor ?uj/?xi is separated in divergence D, vorticity V, strain St, and shear Sh:

... (5)

... (6)

That decomposition is sketched in Fig. 2.

The action anomaly in (3) is then rewritten

... (7)

d. Wind-wave spectral symmetry

The proposed decomposition is valid in any direction. Yet, wind-wave spectral energy is concentrated in the wind direction. Without loss of generality, we will as- sume that the wind direction is aligned with the x axis.

We will also assume that the angular spread of the wave spectrum is symmetrical about the wind direction. The wave spectrum is then an even function of f, and the same property is transferred to the relaxation time tc(k, f). Note that mkN0 is even whereas mfN0 is odd. Any integral over the wave directions f of an odd function will cancel out in the mss response (Kudryavtsev et al. 2005). Consequently, this results in the following:

* The mss (any mss, i.e., mss x , mss y , or mss t ) response to vortical currents (V ) is zero, since mf is an odd function of f. This physically means that mss anomaly due to the modulation of waves traveling in the direction f is exactly compensated by that of waves traveling in the direction 2 f.

* The mss response to sheared currents (Sh) is also zero, since sin(2f)mk 1 cos(2f)mf is also odd.

* The mss response to divergent currents (D) is non- zero, since mk is even.

* The mss response to strained currents ( St ) is nonzero, since cos(2f)mk 2 sin(2f)mf is also even.

Vortical and sheared currents will thus not appear as surface roughness variations. Only currents with di- vergence or with strain in the wind direction will have a signature on the mss.

These results are translated onto the canonical base of the tensor [(3)] as follows: ?u/?y and ?y /?x will have no signature in terms of mss, while ?u/?x and ?y/?y will have a nonzero signature. A polarization in the wind direction (i.e., the nonzero response to strain in that direction) appears because the responses to ?u/?x and to ?y /?y are generally not equal. This will be quantified in section 3a by defining a polarization index.

e. Wind direction dependency

Divergence D and vorticity V of the currents are in- dependent of the axis orientation. This is not the case for the strain St and shear Sh. For a change of coordinates, if u is the rotation to the new coordinates,3 the strain and shear are not conserved and transform as follows:

... (8)

This explicitly shows that the strain in one direction is related to shear in a direction rotated by 458. The two orthogonal directions for which S0h 5 0 define the major and minor axes of the deformation (e.g., Kundu 1990, p. 60).

To illustrate our purpose, consider an infinite front with alongfront current [u 5 0; y 5 y(x); see sketch 5 in Fig. 2]. This current has no divergence and no strain in the x direction, D 5 St 5 0, but has nonzero vorticity and nonzero shear in the x direction, V 5 Sh 5 ?y/?x.If the wind blows at an angle u, the strain in the wind direction reads St0 52sin(2u)Sh according to (8). Such a front will therefore not appear in roughness images when the wind blows in the north-south or east-west directions. But, if the wind blows obliquely (the maxi- mum being at 458 angle), the current will have a non- zero strain in the wind direction and will thus appear in roughness images.

In summary, divergent currents with no strain and no shear will appear in surface roughness images in- dependently of the wind direction, with a roughness anomaly independent of the wind direction (see Fig. 2). For nondivergent currents, the relative direction of the wind is crucial. Only currents strained in the wind di- rection will exhibit an mss signature, while currents strained 458 to the wind will not appear. If the wind turns by 908, the signature of strained currents will further change sign.

3. Toward a practical interpretation of surface roughness variations

The analysis presented so far demonstrates that some current deformations appear in surface roughness im- ages. Analysis of such images in the sun glitter spot (Cox and Munk 1954), or use of ancillary sources (e.g., models or other satellite data), can provide additional in- formation about the wind speed and direction. It is then of interest what information we can extract about the current deformation field when we combine all this information.

a. Polarization index

In this section the wind is set again in the x direction. Because the mss responses to ?u/?x and to ?y /?y are generally not equal, we introduce a new parameter, the polarization index, to quantify the relative contribution of each current gradient. The total mss response to current can be rewritten from (3) and (4) as

... (9)

where we define the polarization index as

... (10)

We can define similarly the polarization indexes ax and ay for the x and y components of the mss.

Note that a combination ?u/?x 1 a?y/?y corresponds to a divergence and strain combination in the form D 1 bSt, with b 5 (1 2 a)/(1 1 a):

(i) If a ' 1 (i.e., b ' 0), the roughness equally adjusts to ?u/?x and ?y /?y. The roughness patterns represent divergence, and this is independent of the wind/ current angle. Nondivergent currents do not appear even if they are strained in one direction. The imaging mechanism, that is, the mechanism to detect currents from roughness variations, is thus unpolarized.

(ii) If a ^ 1 (i.e., b ' 1), the roughness adjusts to ?u/?x and is insensitive to ?y /?y. The imaging mechanism is completely polarized in the wind direction, that is, roughness patterns solely represent current gradients in the wind direction ?u/?x.

(iii) If a ^ 1 (i.e., b ' 21), the imaging mechanism is completely polarized in the crosswind direction.

b. Polarization and wave angular spread

The polarization index mainly depends on the angular spread of the wave spectrum to be dependent on the wavenumber (e.g., Elfouhaily et al. 1997; Yurovskaya et al. 2013). The present subsection illustrates this by considering an idealized spectrum of background wind sea of the form

... (11)

where the wavenumber exponent mk is 24.5 for a saturation spectrum (Phillips 1958) and where the directionality cosine exponent n takes typical values from around 1 for a broad spectrum (i.e., with wide- spread directions of wave propagations) to 9 for a narrow spectrum (i.e., with waves propagating mainly along the wind direction).

The angular spread of the relaxation time t c has to be prescribed. To do so we argue that, related to the growth rate parameters, waves traveling in the wind direction are more closely related to the local wind and adjust quicker than waves near the crosswind direction. For simplicity, the relaxation time is assumed to be separa- ble, that is, tc(k, f) 5 t 1(k)t2(f). By analogy with the relaxation time of Kudryavtsev et al. [2005, their (40)] we express

... (12)

Physically, (12) means that waves have a relaxation time independent of their direction except for waves within 458 of the wind that have a shorter relaxation time.

Polarization parameters ax, ay, and at (for the upwind, crosswind, and total mss, respectively) are shown in Fig. 3a for our idealized wave spectrum [(11-12)], as a function of the directional spread. Prominent features are the following:

* Very directional wave spectra ( n ' 9) mainly adjust to ?u/?x, with any a (ax, ay,orat) of the order of 0.1-0.5. This is expected as there are very few waves traveling in the crosswind direction y in that case. On the contrary, less directional spectra (smaller n) are more sensitive to ?y /?y.

* The upwind mss, as expected, mainly adjusts to ?u/?x , with ax of the order of 0.1-0.6 (thin solid line).

* The crosswind mss principally responds to ? u / ? x for narrow spectra (n ' 9), whereas for broad spectra (n ' 1) it principally responds to ?y/?y.

c. Strain versus divergence signature

Written in the form D 1 bSt, with b 5 (1 2 a)/(1 1 a), the polarization index gives information on the relative mss response to strain and divergence of similar ampli- tude. That is, a current divergence of magnitude D (in s21) produces an mss response similar to that of a current strain of magnitude St 5 D / b .

The polarization index b is shown in Fig. 3b from which two important points are made:

* The magnitude of b is of the order of 0.5, which means that the response to a current strain is smaller but has a similar order of magnitude as that of divergence. Quasigeostrophic surface currents have a small di- vergence compared to strain and are therefore ex- pected to be observed due to their strain in the wind direction.

* Waves traveling in all directions become less steep when currents are divergent. On the contrary, when currents are strained in the wind direction, waves traveling in the wind direction become less steep, whereas waves traveling in the crosswind direction steepen. Under certain conditions (especially for crosswind-looking sensors), we thus expect b to be- come negative, that is, we expect that the mss response to strain takes an opposite sign to that of the response to divergence. As shown in Fig. 3, such conditions occur when the wave spectrum is broad (e.g., n , 4) for crosswind mss and when it is very broad (e.g., n , 1) for the total mss. For upwind mss, such conditions never occur and b is always positive.

d. Polarization for a complete spectrum

For more realistic spectra, the energy-containing waves are quite directional, while smaller waves tend to be more omnidirectional (Donelan et al. 1985). In turn, short gravity waves with parasitic capillaries are often reported to gain directionality (Elfouhaily et al. 1997; Kudryavtsev et al. 1999). For our purpose, we use the background spectral form proposed by Kudryavtsev et al. (2005) constrained by recent ob- servations (Yurovskaya et al. 2013). The polarization parameters a computed with this spectrum are shown in Fig. 4a. For most wavenumbers (k . 2m21), the polar- ization corresponds to a case of broad spectrum (n ' 1) in (11). The reason is that a cosh(u)22 (Donelan et al. 1985) is flatter than a cosn(u) angular distribution. The crosswind mss response is then always dominated by ?y /?y (Fig. 4a, dashed line) and thus the imaging mech- anism for a crosswind-looking sensor is strongly polarized in the crosswind direction (ay . 1). On the contrary, the imaging mechanism for an upwind-looking sensor is solely polarized in the upwind direction (ax , 1).

e. Generation of surface roughness by breaking waves

Kudryavtsev et al. (2005) introduced an additional source term to take into account short-scale generation by longer breaking waves. The mechanism was shown to help in reproducing observed large variations of mss. Waves of intermediate scales, which are indeed experi- encing the largest modulations in the presence of cur- rent, can further enhance the overall roughness through their breaking modulation.

Accounting for the effect of breaking waves on sur- face roughness leads to polarization parameters lower than 1 (a ' 0.44), for both upwind and crosswind mss (see Fig. 4b). Short waves created by larger, more di- rectional breaking waves become indeed more sensitive to ?u/?x. Note that polarization parameters for the up- wind, crosswind, and total mss are then very similar for all short waves with k . 5m21.

f. Practical summary

The typical mss response to different current de- formations is sketched on Fig. 2. To simplify, we as- sume that for the upwind (crosswind) mss, the imaging mechanism is polarized in the upwind (crosswind) direction.

Figure 2 is therefore not valid for the case corre- sponding to section 3e (for small waves, k . 5m21; see Fig. 4b), where for the crosswind mss the imaging mechanism is polarized in the upwind direction. In that case the crosswind mss would have the same response as that of the upwind mss.

In the special case of an unpolarized imaging mecha- nism, the responses to strain and shear would be zero.

4. Illustration

The polarization sensitivity is illustrated next with a numerical simulation of mesoscale and submesoscale oceanic turbulence in a beta-plane channel performed with a primitive equation model. The simulation is forced by an unstable westerly zonal flow. The turbu- lence is intensified at the surface where the sub- mesoscale dynamics are active and the spectrum of horizontal kinetic energy is shallow (slope in (2p/L)22). The averaged horizontal kinetic energy is about 10 cm s21 at the surface. The horizontal resolution of the model is 2 km, and the vertical resolution is about 3 m near the surface. More details can be found in Klein et al. (2008). Surface currents are shown in Fig. 5a.

The forward simulation of the mss variations uses the model of Kudryavtsev et al. (2005) and includes wave propagation [left-hand side of (1)]. The horizontal res- olution of the wave model is also 2 km. The wind is westerly 6 m s21, and background waves are supposed fully developed. We consider here the mss produced by waves in the wavenumber range 5 , k , 28 m21.

The resulting upwind mss is shown in Fig. 5b. The current convergence (52D) is shown in Fig. 5c. It is not correlated to the upwind mss (corr ' 20.08). The av- erage value of the polarization parameter ax in this range of short waves is 0.44 (Fig. 4b, thin solid line), which corresponds to bx 5 0.39. The field 2?u/?x 2 0.44?y/?y is shown in Fig. 5d. As can be visually checked, it is correlated with the upwind mss (corr ' 0.56). The latter correlation is not perfect, but improves when the long-wave propagation is omitted (corr . 0.9). This re- sult validates the decomposition of the right-hand side of (1) with our polarization parameter.

Yet, even without propagation effects, the correla- tion between upwind mss and convergence remains very low (corr ' 0.2). The upwind mss anomaly mssx can thus be considered as a good tracer of ?u/?x 1 0.44?y/?y. This field corresponds to D 1 0.39St.Inthis submesoscale turbulence simulation, strain is roughly 3 to 4 times larger than divergence. Strain and di- vergence thus make equivalent contributions to the upwind mss anomaly and cannot be distinguished from asingleimage.

One can exploit properties of oceanic turbulence in order to advance further. Namely, the magnitude of divergence compared to that of strain varies with the horizontal scale L.InFig. 6, we show the energy spectra of D and of 0.39St. Strain contribution is larger at scales L . 25 km, whereas divergence contribution is larger at scales L , 25 km. For this particular oceanic turbulence, the upwind mss variations shall pro- vide a reasonable source of information related to the divergence at high resolution (L , 25 km) and to thestraininthewinddirectionatlargerscales(L . 25 km).

5. Conclusions

Surface wind waves interact with ocean currents to reveal specific deformation properties of surface cur- rents. Considering surface roughness images, we limited our discussion to variations of wave mean square slope (mss). We ignored other mechanisms that potentially affect surface roughness, such as surfactants and atmo- spheric boundary layer modifications. Under a classical relaxation approach, the surface current gradient per- turbs the background equilibrium sea state. We focused on waves of short and intermediate scales, where the propagation can be neglected compared to the current length scale. We also supposed that the wave spectrum is symmetrical about the wind direction.

Accordingly, the mss is modified by the along-wind current gradient ?u/?x or by the crosswind current gra- dient ?y/?y. The other terms of the deformation field, ?u/?y and ?y /?x, have no impact on the mss. This shows that currents with pure vorticity or with pure shear in the wind direction will not appear in roughness images.

The mss responses to ?u/?x and to ?y /?y are not equal in general, and their ratio defined a polarization pa- rameter a.Ifa 6 1, the imaging mechanism is called ''polarized'' in the sense that current detection depends on the relative direction between the wind and the current feature.

A convenient way to look at the imaging mechanism when polarized is through divergence and strain of the current. Divergent currents appear in any wind di- rection. On the contrary, only currents strained in the wind direction exhibit an mss variation, while currents strained 458 to the wind do not appear. As such, the wind direction must be taken into account.

Additionally, the polarization parameter depends on whether the observation geometry favors upwind or crosswind components of the mss. Divergent (convergent) currents shall always appear as negative (positive) mss anomalies, while currents strained in the wind direction shall be sensor look (upwind or crosswind) dependent. This latter aspect further depends upon the directional spread of short wind waves and the interactions between short waves and longer breaking waves.

Thus, a very useful application of surface roughness images shall be to help diagnose upper-ocean vertical velocities due to submesoscale turbulence. As discussed, inversion of a single roughness image to obtain current divergence is certainly not straightforward. Especially when the imaging mechanism is polarized, the signature of the strain can be comparable in magnitude (around 30%) to that of divergence of similar amplitude. Yet, when using observations with very high resolution, the analysis can concentrate on horizontal scales for which divergence makes the dominant contribution. For larger scales, additional observations are needed to assess and possibly remove the strain signatures. The present analysis suggests the potential use of multiple observa- tions with different azimuthal look angles of the sensor (for conditions and sensors for which the crosswind mss imaging is polarized in the crosswind direction) or multiple delayed observations with 458 or 908 different wind directions. The use of multiple look directions or possible acquisition with different wind directions can then help to retrieve surface currents from surface roughness observations.

Nevertheless, if individual images of surface rough- ness are to be used to characterize turbulence regimes in the upper ocean, a suitable diagnostic variable related to the wind direction, that is, ?u/?x 1 a?y /?y including the polarization parameter a, shall be used rather than the isotropic current divergence ?u/?x 1 ?y/?y.

Finally, in the perspective of future wide swath radar altimetry (e.g., the SWOT mission) for which the sensor physics is principally dominated by the mss strength, unique opportunities shall be offered to combine high- resolution (1 km) detected surface roughness changes with sea surface elevation measurements at about 5-km resolution.

Acknowledgments. We acknowledge the financial support of the French Agence Nationale pour la Recherche (ANR) through the REDHOTS project, of the European Space Agency (ESA) through the STSE MESO3D project and through the GlobCurrent project, and of the French Laboratoire d'Excellence (LabexMer) Axe1. We also acknowledge V. Kudryavtsev, RSHU, SPb, for the fruitful discussions and help to implement the numerical forward simulation and F. Collard, Ocean- DataLab, to motivate the analysis and interpretation of numerous observations.

1 Under moderate wind conditions, there is also a local coupling between the atmosphere flow and surface roughness. This effect shall follow roughness variations and further enhance them, as the form drag scaled by the squared friction velocity is directly linked to the mss (Plant 1982).

2 In the literature, ''strain'' is also called ''stretch'' (e.g., Okubo 1970) or ''normal strain'' (e.g., Hua 1994) and ''shear'' is also called ''shear strain.''

3 For a change of coordinates, Q 5 cosu 2sinu is the u ro- sinu cosu tation tensor to the new coordinates. A deformation tensor T changes in^the new^ coordinates to T0 5 QtTQ.^ The dive^rgence tensor D 5 D0 and the vorticity tensor V 5 2V of the 0D 0 current are then conserved, that is, D0 5 D and V0 5 V, because t^hey comm^ute with Q. On the cont^rary, the^strain tensor St 5 St 0 0 Sh 0 St and the shear tensor Sh 5 Sh 0 are not conserved.


Alpers, W., 1985: Theory of radar imaging of internal waves. Na- ture, 314, 245-247, doi:10.1038/314245a0.

_____, and I. Hennings, 1984: A theory of the imaging mechanism of underwater bottom topography by real and synthetic aper- ture radar. J. Geophys. Res., 89, 10 529-10 546, doi:10.1029/ JC089iC06p10529.

Apel, J. R., H. M. Byrne, J. R. Proni, and R. L. Charnell, 1975: Observations of oceanic internal and surface waves from the earth resources technology satellite. J. Geophys. Res., 80, 865- 881, doi:10.1029/JC080i006p00865.

Beal, R. C., P. S. DeLeonibus, and I. Katz, Eds., 1981: Spaceborne Synthetic Aperture Radar for Oceanography. Vol. 7. Johns Hopkins University Press, 215 pp.

_____, V. N. Kudryavtsev, D. R. Thompson, S. A. Grodsky, D. G. Tilley, V. A. Dulov, and H. C. Graber, 1997: The influence of the marine atmospheric boundary layer on ERS 1 synthetic aperture radar imagery of the Gulf Stream. J. Geophys. Res., 102, 5799-5814, doi:10.1029/96JC03109.

Cox, C., and W. Munk, 1954: Measurement of the roughness of the sea surface from photographs of the sun's glitter. J. Opt. Soc. Amer., 44, 838-850, doi:10.1364/JOSA.44.000838.

Donelan, M. A., J. Hamilton, and W. H. Hui, 1985: Directional spectra of wind-generated waves. Philos. Trans. Roy. Soc. London, A315, 509-562, doi:10.1098/rsta.1985.0054.

Dulov, V., and V. Kudryavtsev, 1990: Imagery of the inhomo- geneities of currents on the ocean surface state. Sov. J. Phys. Oceanogr., 1, 325-336, doi:10.1007/BF02196830.

Elfouhaily, T., B. Chapron, K. Katsaros, and D. Vandemark, 1997: A unified directional spectrum for long and short wind-driven waves. J. Geophys. Res., 102, 15 781-15796, doi:10.1029/ 97JC00467.

Espedal, H. A., O. M. Johannessen, J. A. Johannessen, E. Dano, D. Lyzenga, and J. Knulst, 1998: COASTWATCH'95: ERS 1/2 SAR detection of natural film on the ocean surface. J. Geophys. Res., 103, 24969-24982, doi:10.1029/98JC01660.

Ferrari, R., 2011: A frontal challenge for climate models. Science, 332, 316-317, doi:10.1126/science.1203632.

Fox-Kemper, B., and Coauthors, 2011: Parameterization of mixed layer eddies. III: Implementation and impact in global ocean climate simulations. Ocean Modell., 39, 61-78, doi:10.1016/ j.ocemod.2010.09.002.

Fu, L.-L., and B. Holt, 1983: Some examples of detection of oceanic mesoscale eddies by the SEASAT synthetic-aperture radar. J. Geophys. Res., 88, 1844-1852, doi:10.1029/JC088iC03p01844.

_____, D. Alsdorf, R. Morrow, E. Rodriguez, and N. Mognard, 2012: SWOT: The Surface Water and Ocean Topography mission: Wide-swath altimetric measurement of water elevation on Earth. Jet Propulsion Laboratory JPL Publ. 12-5, 228 pp. [Available online at]

Hua, B. L., 1994: Skewness of the generalized centrifugal force divergence for a joint normal distribution of strain and vor- ticity components. Phys. Fluids, 6, 3200, doi:10.1063/1.868101.

Hughes, B., 1978: The effect of internal waves on surface wind waves 2. Theoretical analysis. J. Geophys. Res., 83, 455-465, doi:10.1029/JC083iC01p00455.

Huot, J., M. Rast, S. Delwart, J. Bezy, G. Levrini, and H. Tait, 2001: The optical imaging instruments and their applications: AATSRand MERIS.ESABull.,106,56-66.[Available online at Bulletin-106/.]

Jansen, R. W., C. Y. Shen, S. R. Chubb, A. L. Cooper, and T. Evans, 1998: Subsurface, surface, and radar modeling of a Gulf Stream current convergence. J. Geophys. Res., 103, 18723-18743, doi:10.1029/98JC01195.

Johannessen, J. A., R. A. Shuchman, G. Digranes, D. Lyzenga, C. Wackerman, O. M. Johannessen, and P. Vachon, 1996: Coastal ocean fronts and eddies imaged with ERS 1 synthetic aperture radar. J. Geophys. Res., 101, 6651-6667, doi:10.1029/ 95JC02962.

_____, V. Kudryavtsev, D. Akimov, T. Eldevik, N. Winther, and B. Chapron, 2005: On radar imaging of current features: 2. Mesoscale eddy and current front detection. J. Geophys. Res., 110, C07017, doi:10.1029/2004JC002802.

Keller, W., and J. Wright, 1975: Microwave scattering and the straining of wind-generated waves. Radio Sci., 10, 139-147, doi:10.1029/RS010i002p00139.

Klein, P., B. Hua, G. Lapeyre, X. Capet, S. Le Gentil, and H. Sasaki, 2008: Upper ocean turbulence from high-resolution 3D sim- ulations. J. Phys. Oceanogr., 38, 1748-1763, doi:10.1175/ 2007JPO3773.1.

Kudryavtsev, V. N., V. K. Makin, and B. Chapron, 1999: Coupled sea surface-atmosphere model: 2. Spectrum of short wind waves. J. Geophys. Res., 104, 7625-7639, doi:10.1029/1999JC900005.

Kudryavtsev, V., D. Akimov, J. Johannessen, and B. Chapron, 2005: On radar imaging of current features: 1. Model and comparison with observations. J. Geophys. Res., 110, C07016, doi:10.1029/2004JC002505.

_____, A. Myasoedov, B. Chapron, J. A. Johannessen, and F. Collard, 2012a: Joint sun-glitter and radar imagery of sur- face slicks. Remote Sens. Environ., 120, 123-132, doi:10.1016/ j.rse.2011.06.029.

_____, _____, _____, _____, and _____, 2012b: Imaging mesoscale upper ocean dynamics using synthetic aperture radar and optical data. J. Geophys. Res., 117, C04029, doi:10.1029/ 2011JC007492.

Kundu, P., 1990: Fluid Mechanics. Academic Press, 638 pp.

LÉvy, M., R. Ferrari, P. J. Franks, A. P. Martin, and P. RiviÈre, 2012: Bringing physics to life at the submesoscale. Geophys. Res. Lett., 39, L14602, doi:10.1029/2012GL052756.

McWilliams, J., F. Colas, and M. Molemaker, 2009: Cold filamen- tary intensification and oceanic surface convergence lines. Geophys. Res. Lett., 36, L18602, doi:10.1029/2009GL039402.

Okubo, A., 1970: Horizontal dispersion of floatable particles in the vicinity of velocity singularities such as convergences. Deep-Sea Res. Oceanogr. Abstr., 17, 445-454, doi:10.1016/ 0011-7471(70)90059-8.

Phillips, O. M., 1958: The equilibrium range in the spectrum of wind-generated waves. J. Fluid Mech., 4, 426-433, doi:10.1017/ S0022112058000550.

_____, 1984: On the response of short ocean wave components at a fixed wavenumber to ocean current variations. J. Phys. Oceanogr., 14, 1425-1433, doi:10.1175/1520-0485(1984)014,1425: OTROSO.2.0.CO;2.

Plant, W. J., 1982: A relationship between wind stress and wave slope. J. Geophys. Res., 87, 1961-1967, doi:10.1029/ JC087iC03p01961.

Yurovskaya, M., V. Dulov, B. Chapron, and V. Kudryavtsev, 2013: Directional short wind wave spectra derived from the sea surface photography. J. Geophys. Res., 118, 4380-4394, doi:10.1002/ jgrc.20296.


Laboratoire d'OcÉanographie Spatiale, Institut Francais de Recherche pour l'Exploitation de la Mer, PlouzanÉ, France,

and Nansen-Tutu Center, Department of Oceanography, University of Cape Town, Rondebosch, South Africa


Laboratoire d'OcÉanographie Spatiale, Institut Francais de Recherche pour l'Exploitation de la Mer, PlouzanÉ, France


Laboratoire de Physique des OcÉans, Institut Francais de Recherche pour l'Exploitation de la Mer, PlouzanÉ, France


Laboratoire d'OcÉanographie Spatiale, Institut Francais de Recherche pour l'Exploitation de la Mer, PlouzanÉ, France


Laboratoire de Physique des OcÉans, Institut Francais de Recherche pour l'Exploitation de la Mer, PlouzanÉ, France

(Manuscript received 16 December 2013, in final form 25 April 2014)

Corresponding author address: Nicolas Rascle, Laboratoire d'OcÉanographie Spatiale, IFREMER, ZI Pointe du Diable, 29280 PlouzanÉ, France.


For more stories covering the world of technology, please see HispanicBusiness' Tech Channel

Source: Journal of Physical Oceanography

Story Tools Facebook Linkedin Twitter RSS Feed Email Alerts & Newsletters