## 1. Introduction

The ocean circulation and its water mass distribution and variability present different length and time scales. We can distinguish length scales from the order of the basin size to mesoscale (some hundred kilometers), and from a continuum of downscaling length scales down to microturbulence. Our knowledge about the spatial distribution of the different properties of the ocean comes mainly from in situ measurements and satellite data, and more recently from autonomous profilers. One of the main techniques used in the analysis of the data from oceanographic surveys is their interpolation into a regular grid for contouring and for further calculations. This provides information about the water mass distributions as well as estimates of water, heat, salt, etc., fluxes.

Some of the most frequently used interpolation methods are the distance-weighting and optimal interpolation (OI) schemes (Thiébaux and Pedder 1987 for a review of both methods). In both methods the length scales that can be solved are limited by the separation between observations. Obviously we cannot solve for wavelengths shorter than twice the separation between stations. Weighting schemes allow the filtering out of those wavelengths in which we are not interested, or simple separation thereof (Maddox 1980; Gomis and Alonso 1990; Tintoré et al. 1991). Optimal interpolation also allows some sort of filtering through the choice of the noise-to-signal ratio and the decaying distance of the used autocovariance function (Bretherton et al. 1976; Roemmich 1983; Pedder 1989; Fukimori and Wunsch 1991; Gomis et al. 2001; Lavín 1999). The length scales that are retained in the analysis depend on which is the phenomenon of study. Mesoscale, and even submesoscale, eddies will have great importance if the interest is in the physical–biological coupling in the upper hundred meters (Rodríguez et al. 2001; Barton et al. 1998; McGillicuddy and Robinson 1997; Fernández and Pingree 1996). On the other hand, if the interest is in some property average, the eddy-induced variability will simply be considered as a noise.

A good example of the eddy-induced variability around a mean value is found in the historical study of oceanographic sections. One of the best sampled of these sections is the one along the 24.5°N parallel in the Atlantic, from the African coast to the Bahamas, which has been sampled in 1957, 1981, 1992, and 1998 (Fuglister 1960; Roemmich and Wunsch 1984, 1985; Parrilla et al. 1994; Bryden et al. 1996; Lavín et al. 2003; Baringer and Mollinari 1999). In these works it has been shown that there is a warming in the intermediate layer (800–3000 m), which is represented by a trend in the potential temperature mean value (zonally averaged). The temperature change is obscured by the shorter time-scale variability that is induced by the eddy field; therefore, it was filtered using a Gaussian filter (Roemmich and Wunsch 1984; Parrilla et al. 1994). This eddy-induced noise also imposes a limitation to the significance of the temperature trends that are detected from 1957 to 1992.

It is usually assumed that the dataset obtained in an oceanographic survey is synoptic, independently of which the interpolation method is used, and the length scale of the phenomenon is of interest. Distributions of temperature, salinity, etc., usually exhibit strong eddy activity (Fig. 2 in Bryden et al. 1996 or Fig. 5 in Lavín et al. 2003), with the separation between stations or, as already mentioned, a particular interest being the only limitation for their resolution. Nevertheless, these structures can have time scales that are shorter than the time needed to accomplish the complete sampling. In this case the mesoscale field left after interpolation, and contouring would not give any real picture of the spatial domain covered by the survey. However, this eddy field can still be representative of certain statistical properties of the real one, such as the main length scales or the variance of the field.

The problem of determining whether or not a certain dataset can be considered synoptic is difficult, because oceanographic surveys are not frequently repeated. One exception can be seen in Gomis et al. (2001), where an oceanographic survey was repeated 4 times in a 14- day period and a synopticity test for each survey could be addressed. The area that was covered in that case was quite reduced (100 km × 80 km). A similar test, using traditional sampling techniques, is very difficult for large areas of the ocean.

In March 2002, the Instituto Español de Oceanografía (IEO) deployed (within the European project Gyroscope) 19 autonomous profilers in the Canary Basin. During the deployment, 38 CTD profiles were collected. One of the cruise legs was along 24.5°N, covering the eastern half of the section that was already sampled in 1957, 1981, 1992, and 1998. The temperature and salinity fields that are obtained from the CTD casts, and the first data transmitted by the profiling floats with a 5-day delay, provided an exceptional opportunity to check the time variability of a hydrologic transatlantic section. The main objective of this paper is to show that mesoscale structures can change in a 5-day period. Taking into account that this time is shorter than the time that is needed to accomplish a transatlantic section, it is implied that the data cannot be considered as synoptic. The direct comparison of temperature and salinity profiles obtained with a CTD probe and profiling floats at the same location shows that the differences observed are a result of the high temporal variability of the mesoscale field. This result also implies that it is not possible to compare different instruments, such as CTD probes and profiling floats, at isobaric levels if there is a time delay between both datasets. As a secondary objective it is shown that the CTD–profiler comparison should be made by analyzing the *θ*–*S* relationship that is obtained for the water masses sampled with the CTD probe and the profilers.

The paper is organized as follows. Section 2 presents the dataset. Section 3 contains a comparison of the two repetitions of the 24.5°N section and the direct comparison of the potential temperature and salinity profiles obtained at the same locations. Section 4 deals with the estimation of the mesoscale time variability using the time series provided by the profiling floats. In section 5 the performance of the profiling floats and the CTD probe are compared, considering the *θ*–*S* diagrams that are obtained from both datasets. A summary and conclusions are given in section 6.

## 2. Data

The cruise Gyroscope E 03–2002 (Grupo Gyroscope-España 2002), carried out by IEO in March 2002 within the homonymous European project, was planned to deploy 20 autonomous profilers, equipped with pressure, temperature, and conductivity sensors, along a box bounded by the parallels between 24.5° and 30°N and the meridians between 18° and 42°W. The southern leg of the box ran from the African coast to the Mid-Atlantic Ridge. It coincides with the eastern half of the sections that were previously sampled in 1957, 1981, 1992, and 1998.

The profiler set consisted of 15 PROVOR and 5 autonomous profiler explorer (APEX) types. The main difference is that the latter ones are equipped with a salinity module, providing salinity measurements directly, instead of conductivity. Both types of profilers are programed to drift at a 1500-m depth and sink to a 2000-m depth and emerge to the surface recording a CTD profile every 10 days. One of the PROVOR profilers failed because of a programming error, so the final set was of 19 profilers (Fig. 1).

Additionally, it was also planned to make CTD stations along the box. Unfortunately, the loss of the probe limited the planned number of casts to 38 (Fig. 1). The CTD probe was equipped with redundant conductivity and temperature sensors. Water samples were also obtained for salinity calibration and oxygen and nutrient determination. The redundant sensors and water samples were used to calibrate the CTD to World Ocean Circulation Experiment (WOCE) standards (see Grupo Gyroscope-España 2002 for technical details).

The only CTD section completed was that along 24.5°N, the one on which we will focus in section 3. Nevertheless we will also use CTD and profiler data along the western boundary of the box and the northern section (uncompleted) for comparison of the CTD and profiler data (section 5).

The CTD was lowered to the bottom depth, which varies between some 3300 m near the African slope and over the Mid-Atlantic Ridge, and more than 5600 m in the abyssal plane. CTD profiles along the 24.5°N section were taken from 4 to15 March 2002. The profiles from the first section were completed by the autonomous profilers, which correspond to 8–20 March 2002. The central time point in the CTD section is 1200 UTC 9 March 2002, and that in the profiler section is 14 March 2002. There is a delay close to 5 days between both of them. If there were differences in the length scales that the station separation would be able to solve, it would indicate that mesoscale structures can evolve in a 5-day period, which is shorter than the time needed to complete the whole section—12 days for the CTD and 13 for the profiler section. Note that the three previous cruises completed the section from the African coast to the Bahamas in 23, 24, and 26 days, respectively. Because the complete section has approximately double the length of that of the section of the Gyroscope cruise, the time needed to sample the sections is similar. If a 12-day period is too long for ignoring the time evolution of mesoscale structures, the situation is even worse for the 24-day-long sections.

## 3. CTD and profiler section comparison

*D*= 250 km and

_{x}*D*= 33 m;

_{y}**u,**are any couple of position vectors within the vertical section.

*υ*Figures 2a and 2c show the potential temperature and salinity fields plotted from the CTD data (solid line) and from the first transmission of the profiling floats (dashed line). Figures 2b and 2d show the differences between both datasets. These differences range between −0.7° and 0.7°C for potential temperature and from −0.14 to 0.15 for salinity. The background fields (not shown), which capture those structures with length scales of the order of the analyzed section, are almost identical for both datasets. Differences seem to be associated with the different shape and position of the mesoscale structures. To enhance these differences we present a zoom of the upper 200 m (Figs. 2e, f).

*θ*is the interpolated increment field

*j, i*denotes the horizontal and vertical positions within the grid,

*N*and

_{x}*N*are the number of grid points on the horizontal and vertical directions, and

_{y}*L*and

_{x}*L*the horizontal and vertical lengths of the grid. Each plane wave has a wavenumber

_{y}**k**= 2

*π*[(

*k*/

*L*

_{x}), (

*l*/

*L*)] and wavelength

_{y}*λ*

_{k}= 2

*π*/|

**k**|. The averaged power of the signal corresponding to each wave is |

*C*

_{k,l}|

^{2}. Figure 3b shows the power associated with each plane wave as a function of

*1*/

*λ*for the interpolated temperature field. Figure 3c is the same but selects only the horizontal plane waves, that is, those with

*l*= 0 [see expression (2)], because it is on the horizontal direction where both datasets have a different resolution. We see that in both cases the power of the signal is limited to values of

*1*/

*λ*lower than 0.0018 km

^{−1}, that is, there is no energy associated with the wavelengths shorter than 550 km. This first result indicates that the differences outlined in Fig. 2 are not a result of the different spatial resolution of both datasets.

Observed differences could also be a result of the interpolation error. The OI allows the estimation of the error variance, defined as 〈(*θ*_{g} − *θ̂*_{g})^{2}〉, where *θ*_{g} is the true value at the grid point **g**, *θ̂*_{g} is the estimated value, and the brackets denote expectation. The error variance is estimated at each grid point. Taking the square root of the error variance and averaging over the whole section we get an average standard deviation for the error of 0.14°C for temperature and 0.03 for salinity. Values as large as 0.7°C or 0.14 are not likely to be caused by the interpolation error, so we discard this cause.

The contouring of interpolated fields evidences variations in the mesoscale structures. For instance, Figs. 2e and 2f show that the position where the 22°C isotherm and the 37.2 isohaline outcrop the surface has shifted more than 100 km from one survey to the other. Nevertheless, another possibility is that all of the differences outlined along this section are an artifact of the analysis and interpolation process. The interpolated sections are simply smoothed versions of the coincident CTD and profiler measurements, and a direct comparison of these profiles should show these differences too. To perform this analysis we only use the PROVOR floats for which there are 10 coincident with CTD casts (only three APEX profilers are coincident). The differences between temperature and salinity profiles at coincident locations in the upper 1000-m range between −1° and 1°C and −0.5 and 0.2, respectively. As in the interpolation process the existence of errors in the observations is assumed, the interpolation process smoothes the field as is already evidenced by the spectral response of the interpolated fields (Figs. 3b,c). This explains that the differences are slightly larger in the direct comparison.

The mesoscale field, estimated from the interpolated sections, has a standard deviation that changes along the water column, with the maximum in the upper 200 m where the eddy activity is likely to be more intense. Averaging over the water column, we get the values of 0.18°C and 0.03 for the temperature and the salinity, respectively. If we estimate the standard deviation for the mesoscale field using the PROVOR profiles (once a second-order polynomial is subtracted to account for the background or large-scale field), we get the values of 0.21°C and 0.04 for the temperature and the salinity, respectively. As already stated, the differences between the CTD and profiler data range between −1° and 1°C. The standard deviations of these differences are 0.3°C and 0.05 for the temperature and salinity. Note that the standard deviation of the difference between the CTD and profilers is *σ*^{2}_{CTD} + *σ*^{2}_{P} − *γ*(CTD, *P*)*σ*^{2}_{CTD} is the variance of the CTD field and *σ*^{2}_{P} is the variance of the profiler data, which we assume is the same. Here, *γ* is the covariance between the CTD and profiler data. Both datasets must be correlated because they are only 5 days apart (we will see the time scale of the mesoscale field in the following section). Note that for the profiler data the standard deviation is 0.21, which, multiplied by

## 4. Autocorrelation and integral time scale

In the previous section it has been shown that the mesoscale field can change in a time, which is comparable to that needed to complete half a transatlantic section (in our case, around 12 days). Nevertheless some of the statistics that are used have been derived from a single CTD survey, which is usually the procedure in oceanography. To estimate the covariance between two stations separated a certain distance, let us say *d*, we search all of the pairs of measurements that are separated by *d*± *b*/2, where *b* is the size of the bin. This implies that the statistics of the signal are homogeneous, and we substitute averages over time by averages over a certain spatial domain. Most of the OI methods were initially developed by meteorologists. In meteorology, long time series at different fixed stations are available; the covariance between them can be obtained and we can use it in the OI schemes. Contrary to this situation, in oceanography the background field and the statistics of the signal are frequently estimated from one single oceanographic survey. The use of autonomous profilers can produce more accurate estimations of the signal statistics. Nevertheless, this method still has certain limitations when compared with meteorological stations, because profilers drift and it is more difficult to obtain the correlation between two different locations or the autocorrelation for a single site.

For profiler 4 we fit (3) to 25 time series corresponding to 25 depths from the surface to 2000 m. The fit is represented in Fig. 4a by the dashed lines. Because profiler 4 was deployed close to the eastern boundary of the section, a similar analysis was applied to profiler 10, in the western edge of it, and no differences were found in the following results.

Once the low-frequency variability was subtracted at the 25 depths selected, the autocovariance function for the residuals that represent the mesoscale field was estimated. Figure 4b shows the autocorrelation function at 5-, 50-, 100-, 250-, 500-, 1000-, 1500-, and 2000-m depths. We only represent these eight depths for clarity; nevertheless, the thick line represents the average of the autocorrelation function for the 25 depths that are analyzed. The lag is expressed in days with an interval of 10 days, which is the time between two consecutive transmissions. A striking result is that the correlation for a 10-day delay is quite low. To obtain another estimation of the decorrelation between hydrological profiles that are delayed in time, we estimated the correlation between the PROVOR temperature profiles over the 24.5°N section and the coincident CTD profiles (after detrending). The correlation obtained was 0.8. PROVOR and CTD profiles were delayed 3 days (APEX profiles were delayed 9 days, which were responsible for the average 5-day delay). Using the profiler time series, we obtained a correlation of 0.1 for a 10-day lag. Interpolating linearly from 1 (at 0 lag) to 0.1 (at 10-day lag), we estimate a value of 0.73 for a 3-day lag, which seems to be coherent with the 0.8 value obtained from the direct comparison of profiler and CTD data.

*ρs*) is the correlation coefficient for the lag

*s*. In practice, the estimated correlation function decays to zero and then remains close to zero, but with positive and negative oscillations that make its estimation uncertain. Freeland (1990) solves this problem by fitting an exponential function to the autocorrelation function that is estimated for the first lags. Lenschow and Stankov (1986) calculated

*τ*(

*r*) = ∫

^{r}

_{0}

*ρ*(

*s*)

*ds*, and then took the first maximum of the function

*τ*(

*r*). Figure 4d shows the integral time scale that is calculated following Lenschow and Stankov (1986) for each of the 25 depths studied. We also computed the integral time scale considering only the autocorrelation function until its first zero and no main changes were observed. The average value is 11 days.

To check the coherence of this result we consider that the length scale is related to the time scale by the mean speed of the flow, which can range between 10 and 30 cm s^{−1} for the eastern boundary and between 5 and 10 cm s^{−1} for the rest of the section (Lavín et al. 2003). For these values, the mesoscale length scale would range between 100–300 and 50–100 km respectively, coinciding in a rough way with an spectral analysis that is applied to the nonsmoothed CTD data (not shown), as well as with previous works (Lavín 1999).

According to these results, two consecutive transmissions of a profiler (10 days apart) are very weakly correlated (concerning the mesoscale field). But what is more interesting follows: if we need 12 days to complete the section, the mesoscale field that is obtained at the beginning and the end of the section is very poorly correlated and synopticity is not a good hypothesis. The situation would be even worse when comparing the eastern and western North Atlantic basins in those surveys that completed the section from Africa to the Bahamas.

## 5. Influences on the intercomparison of CTD and profiler performances

One of the objectives of the Gyroscope E 03–2002 cruise was to compare the performance of the profilers. Because there is a small delay (5 days) between both sections it is implied that the eddy variability is the main source of error as it happened in other cases, for example, the warming detection from 1957 to 1992 (Parrilla et al. 1994). The temperature sensors used in the profilers seem to perform properly, but the same could not be said of the conductivity sensors (Wong et al. 2003). It is clear that the eddy activity can produce large vertical excursions of isopycnals. Roemmich and Wunsch (1985) stated that a single eddy could change the temperature at a 1000-m depth by 2°C over a distance of 200 km. It is clear that it is not sensible to compare temperature or salinity at isobaric levels. On the other hand, if the effect of eddy activity is just to rise (cyclonic eddies) or sink (anticyclonic) the isopycnals, the *θ*–*S* relationships would remain constant. Figure 5a shows the *θ*–*S* diagrams for both the CTD casts and the first transmission of the profilers. Because there are two different types of profilers and each one has a different type of salinity sensor, we used only the PROVOR profilers (14), because a statistical analysis would be more difficult with APEX profilers (only 5). For comparison we used only those CTD casts at the same locations as the PROVOR profilers, that is, profilers 1–13 and the CTD casts coinciding with them.

As noted in previous sections, it is clear that the direct comparison of two coincident profiles makes no sense because of eddy variability. Nevertheless, because the total length of the spatial domain is much larger than the eddy length scale, fluctuations should be canceled out and the average *θ*–*S* relationship should be the same, obviously within the uncertainty imposed by the noise dispersion, which is a result of eddy variability.

Apparently, from Fig. 5a, the salinity values corresponding to the PROVOR profilers are slightly displaced toward lower values than those of the CTD. The mean *θ*–*S* curve has to be modeled by fitting a smooth function to the cloud of *θ*–*S* points. Another possibility is to average the salinity for each temperature value.

Figures 5b and 5c show the result of fitting a third- and fourth-order polynomial to the main thermocline (9°–22°C) and the intermediate (4°–8°C) layers. In the fitted polynomial the salinity is considered as a function of potential temperature, which is assumed to have no error. Though the fitted line is displaced to lower salinities for the profiler data, the differences in the parameters of the fit were not significantly different. This does not allow for the rejection of the null hypothesis that both the profiler and CTD give the same *θ*–*S* relation. Another factor to bear in mind is that the length of the section is 2200 km in the horizontal, which is around 7 times the considered mesoscale length scale of 300 km (Lavín 1999). This makes it easy to find differences in the mean values or in the least squares fit. Another factor to be considered is that confidence intervals were calculated using the usual expressions (see, i.e., Jenkins and Watts 1968, their appendix A4.1), where the autocorrelation of the residuals is not taken into account. The autocorrelation of residuals reduces the degrees of freedom and enlarges the confidence intervals, making it more difficult to reject the null hypothesis. This problem was faced by Roemmich and Wunsch (1984), Parrilla et al. (1994), and Lavín et al. (2003), considering that the number of degrees of freedom is the number of length scales (for the mesoscale) contained in the complete length of the section. This procedure was applied to the isobaric levels where the correlation between consecutive casts was considered. In the least squares fit the situation is not so simple. Two consecutive values of temperature in the cloud of points from Fig. 5a can correspond to nonconsecutive stations and different isobaric levels. For this reason we ignore the autocorrelation of residuals and simply call the attention of the reader to this point. The two factors discussed above support the fact that the differences found in both datasets are compatible with the eddy variability. Contrary to this, we notice that the displacement toward low salinity is observed in the whole water column, in particular in the main thermocline and the intermediate layer.

## 6. Summary and conclusions

Transatlantic sections typically need around 20 days to be completed. The contouring of different properties over these sections show an intense mesoscale activity with length scales of the order of 100 km. The study of the time variability of these structures is not usually addressed, because the repetition in time of these sections is difficult. For this reason it is not usually checked whether or not the datasets obtained in these sections are synoptic. In other words, the mesoscale structures that are sampled at the beginning of the section could have evolved when similar structures are sampled at the end of it. If this were the case, the contoured mesoscale field would not represent a real picture of the section at any time. The Gyroscope 03–2002 cruise was a unique opportunity to obtain two repetitions of such a section within around a 5-day delay. The analysis of two independent datasets that are obtained from the CTD survey and the first transmission of the profilers shows that the mesoscale structures can change in time during the time that is needed to complete the section. The analysis of the interpolated fields and the direct comparison of temperature and salinity profiles obtained at the same locations with some days of delay show differences that cannot be attributed to interpolation errors or the different resolution of both datasets. The range of these differences is coherent with the range of variability of the eddy field. According to our results, the eddy field has an integral time scale of 11 days. Consecutive transmissions of the profilers (10 days apart) are very weakly correlated, and the mesoscale field of the initial and final parts of a transatlantic section cannot be considered as synoptic. The mesoscale field that is typically contoured in this kind of section could represent the statistical properties of the true field (length scales, variance, etc.), but it would not represent any real situation. This would not be very important if we were interested in some average property, because in that case the mesoscale field should be removed, but we think that other studies, for instance, those related to the physical–biological coupling, could be seriously affected by this lack of synopticity.

Our results show that it is not possible to compare temperature and salinity fields that are obtained with different instruments if there is a certain time delay. For instance, for a 3-day delay, the temperature field is decorrelated from 1 to 0.7–0.8, and differences over the isobaric levels range between −1° and 1°C, which is several orders of magnitude higher than the accuracy of the instruments (0.002°C). Instead of this, we propose that the *θ*–*S* properties of the water masses that are sampled should be compared. The average *θ*–*S* obtained with different instruments can be compared, and the time variability imposed by the eddy field could be removed.

## Acknowledgments

We thank Norbert Cortes and David Marcote who deployed the profilers and the crew of *Vizconde de Eza* for their invaluable help during the whole cruise. We also thank Nicolás González, José Escánez, Rosario Carballo, Gerardo Casas, Eugenio Fraile, Venicio Pita, Elena Tel, and Daura Vega who contributed greatly to the acquisition and control of the data. This research has been funded by the European Commission project Gyroscope (EVK2-CT-2000-00087).

## REFERENCES

Baringer, M O., , and Mollinari R. , 1999: Atlantic Ocean baroclinic heat flux at 24 to 26°N.

,*Geophys. Res. Lett.***26****,**353–356.Barton, E D., and Coauthors, 1998: The transition zone of the Canary Current upwelling region.

*Progress in Oceanography*, Vol. 41, Pergamon, 455–504.Bretherton, F P., , Davis R E. , , and Fandry C B. , 1976: A technique for objective analysis and design of oceanographic experiments applied to MODE-73.

,*Deep-Sea Res.***23****,**559–582.Bryden, H., , Griffiths M J. , , Lavín A. , , Millard R C. , , Parrilla G. , , and Smethie W M. , 1996: Decadal changes in water mass characteristics at 24°N in the subtropical North Atlantic Ocean.

,*J. Climate***9****,**3162–3186.Fernández, E., , and Pingree R D. , 1996: Coupling between physical and biological fields in the North Atlantic subtropical front southeast of the Azores.

,*Deep-Sea Res.***43A****,**1369–1393.Freeland, H J., 1990: Sea surface temperatures along the coast of British Columbia: Regional evidence for a warming trend.

,*Can. J. Fish. Aquat. Sci.***47****,**346–350.Fuglister, F C., 1960:

*Atlantic Ocean Atlas of Temperature and Salinity Profiles and Data from the International Geophysical Year of 1957–1958.*Atlas Series, Vol. 1, Woods Hole Oceanographic Institution, 209 pp.Fukimori, I., , and Wunsch C. , 1991: Efficient representation of the North Atlantic hydrographic and chemical distributions.

*Progress in Oceanography*, Vol. 27, Pergamon, 111–195.Gomis, D., , and Alonso S. , 1990: Diagnosis of a cyclogenetic event in the western Mediterranean using an objective technique for scale separation.

,*Mon. Wea. Rev.***118****,**723–736.Gomis, D., , Ruíz S. , , and Pedder M A. , 2001: Diagnostic of the 3D ageostrophic circulation from a multivariate spatial interpolation of CTD and ADCP data.

,*Deep-Sea Res.***48A****,**269–295.Grupo Gyroscope-España, 2002: Informe de la campaña Gyroscope E 03-2002. Datos y resúmenes del Instituto Español de Oceanografía, 25 pp.

Jenkins, G M., , and Watts D G. , 1968:

*Spectral Analysis and Its Applications*. Holden Day, 522 pp.Lavín, A., 1999: Fluxes, trends and decadal changes in the subtropical North Atlantic. Ph.D. dissertation, Universidad de Cantabria, 222 pp.

Lavín, A., , Bryden H. , , and Parrilla G. , 2003: Mechanisms of heat, freshwater, oxygen and nutrient transports and budgets at 24.5°N in the subtropical North Atlantic.

,*Deep-Sea Res.***50A****,**1099–1128.Lenschow, D H., , and Stankov B B. , 1986: Length scales in the convective boundary layer.

,*J. Atmos. Sci.***43****,**1198–1209.Maddox, R A., 1980: An objective technique for separating macroscale and mesoscale features in meteorological data.

,*Mon. Wea. Rev.***108****,**1108–1121.McGillicuddy D J. Jr., , , and Robinson A R. , 1997: Eddy-induced nutrient supply and new production in the Sargasso Sea.

,*Deep-Sea Res.***44A****,**1427–1450.Parrilla, G., , Lavín A. , , Bryden H. , , García M. , , and Millard R. , 1994: Rising temperatures in the subtropical North Atlantic Ocean over the past 35 years.

,*Nature***369****,**48–51.Pedder, M A., 1989: Limited area kinematic analysis by a multivariate statistical interpolation method.

,*Mon. Wea. Rev.***117****,**1695–1708.Rodríguez, J., and Coauthors, 2001: Mesoscale vertical motion and the size structure of phytoplankton in the ocean.

,*Nature***410****,**360–363.Roemmich, D., 1983: Optimal estimation of hydrographic station data and derived fields.

,*J. Phys. Oceanogr.***13****,**1544–1549.Roemmich, D., , and Wunsch C. , 1984: Apparent changes in the climatic state of the deep North Atlantic Ocean.

,*Nature***307****,**447–450.Roemmich, D., , and Wunsch C. , 1985: Two transatlantic sections: Meridional circulation and heat flux in the subtropical North Atlantic Ocean.

,*Deep-Sea Res.***32****,**619–664.Thiébaux, H J., , and Pedder M A. , 1987:

*Spatial Objective Analysis: With Applications in Atmospheric Science*. Academic Press, 300 pp.Tintoré, J., , Gomis D. , , Alonso S. , , and Parrilla G. , 1991: Mesoscale dynamics and vertical motion in the Alborán Sea.

,*J. Phys. Oceanogr.***21****,**811–823.Wong, A. P. S., , Johnson G C. , , and Owens W B. , 2003: Delayed-mode calibration of autonomous CTD profiling float salinity data by

*θ*–*S*climatology.,*J. Atmos. Oceanic Technol.***20****,**308–318.