**APPENDIX I. HEAT FLUXES ESTIMATION**

The surface heat flux computation
was made based on bulk formulae. The net heat flux *Q* through the sea
surface is given by

where *Q _{s}, Q_{b},
Q_{e},* and

*1. Incoming short wave
radiation*

The short wave radiation was estimated following Reed (1977)

where *C* is the fraction
of sky covered by clouds in tenths,
is the altitude of the sun from the horizontal at noon (in degrees), *A =*
0.06 is the albedo, is computed with
the relation sin = sin *l*
sin[23.87 sin(2*π*(*t*–82)/365)] + cos *l*
cos[23.87 sin(2*π*(*t*-82)/365)], where *l*
is the latitude (Reed, 1977). The radiation under clear sky (*Q _{c}*)

from Seckel and Beaudry (1973),
cited in Reed (1977). In equation (21), *ɸ* = (2*π*/365)(*t*–21)
is a function related to the day of the year in which *t* is the Julian
day. The coefficients in equation (24) were computed following Reed (1983).

*2. Long wave radiation*

The long wave radiation was computed as Reed (1983):

In this equation *σ*
= 5.67 x 10^{-8} *Wm*^{–2}*K*^{–4}
is the Stefan Boltzman constant, is
the emissivity of the sea surface, *T* is the temperature in degrees
Kelvin, *e _{a}* is the water vapor pressure (in millibars) computed
following Gill (1982) and is given by

with *H* the percentage
of relative humidity, and *e _{w}* is the saturation vapor pressure
computed as

where *P* is the atmospheric
pressure at sea level in millibars and γ = (0.7859 + 0.03477*t*_{s})/(1
+ 0.00412*t*_{s}), in which *t*_{s} is the temperature
in Celsius.

*3. Latent heat flux*

The latent heat flux was computed with the relation

In this relation *ρ _{a}*
is the air density,

In this work, a constant value
of *C _{E}*

*4. Sensible heat flux*

The sensible heat flux was calculated with the formulae

with *C _{H}* the
sensible heat turbulent coefficient,