Journal cover Journal topic
Ocean Science An interactive open-access journal of the European Geosciences Union
Journal topic
Ocean Sci., 15, 1419–1437, 2019
https://doi.org/10.5194/os-15-1419-2019
Ocean Sci., 15, 1419–1437, 2019
https://doi.org/10.5194/os-15-1419-2019

Research article 07 Nov 2019

Research article | 07 Nov 2019

# The impact of upwelling on the intensification of anticyclonic ocean eddies in the Caribbean Sea

The impact of upwelling on the intensification of anticyclonic ocean eddies in the Caribbean Sea
Carine G. van der Boog1,*, Julie D. Pietrzak1, Henk A. Dijkstra2, Nils Brüggemann5, René M. van Westen2, Rebecca K. James3, Tjeerd J. Bouma3, Riccardo E. M. Riva4, D. Cornelis Slobbe4, Roland Klees4, Marcel Zijlema1, and Caroline A. Katsman1 Carine G. van der Boog et al.
• 1Environmental Fluid Mechanics, Civil Engineering and Geosciences, Delft University of Technology, Stevinweg 1, 2628 CN Delft, the Netherlands
• 2Institute for Marine and Atmospheric Research, Utrecht University, Princetonplein 5, 3584 CC Utrecht, the Netherlands
• 3Estuarine and Delta Systems, Royal Netherlands Institute for Sea Research and Utrecht University, Korringaweg 7, 4401 NT Yerseke, the Netherlands
• 4Geoscience and Remote Sensing, Delft University of Technology, Stevinweg 1, 2628 CN Delft, the Netherlands
• 5University of Hamburg, Bundesstrasse 53, 20146 Hamburg, Germany
• *Invited contribution by Carine G. van der Boog, recipient of the EGU Ocean Sciences Outstanding Student Poster and PICO Award 2018.

Correspondence: Carine van der Boog (c.g.vanderboog@tudelft.nl)

Abstract

The mesoscale variability in the Caribbean Sea is dominated by anticyclonic eddies that are formed in the eastern part of the basin. These anticyclones intensify on their path westward while they pass the coastal upwelling region along the Venezuelan and Colombian coast. In this study, we used a regional model to show that this westward intensification of Caribbean anticyclones is steered by the advection of cold upwelling filaments. Following the thermal wind balance, the increased horizontal density gradients result in an increase in the vertical shear of the anticyclones and in their westward intensification. To assess the impact of variations in upwelling on the anticyclones, several simulations were performed in which the northward Ekman transport (and thus the upwelling strength) is altered. As expected, stronger (weaker) upwelling is associated with stronger (weaker) offshore cooling and a stronger (weaker) westward intensification of the anticyclones. Moreover, the simulations with weaker upwelling show farther advection of the Amazon and Orinoco River plumes into the basin. As a result, in these simulations the horizontal density gradients were predominantly set by horizontal salinity gradients. The importance of the horizontal density gradients driven by temperature, which are associated with the upwelling, increased with increasing upwelling strength. The results of this study highlight that both upwelling and the advection of the river plumes affect the life cycle of mesoscale eddies in the Caribbean Sea.

1 Introduction

The circulation in the Caribbean Sea is characterized by the throughflow of the wind-driven subtropical North Atlantic gyre (Gordon1967), known as the Caribbean Current. This flow is part of the upper branch of the Meridional Overturning Circulation and is highly variable, which manifests itself in meanders and mesoscale eddies. The temporal and spatial variability of the Caribbean Current is influenced by regional differences in wind , freshwater inflow from the Amazon and Orinoco River plumes , resonating Rossby waves , and, at interannual timescales, by the El Niño–Southern Oscillation .

Surface drifter data show that the mesoscale variability in this region is dominated by anticyclonic eddies . These anticyclones are formed in the eastern part of the basin and transport low-salinity anomalies originating from the Amazon and Orinoco River plumes westward . They can be generated from the interaction of the flow with the topography , from the meandering current , from instabilities due to the presence of the river plume , and from perturbations caused by the interaction of North Brazil Current Rings (NBC Rings) with topography . used a model to clarify the dominant generation mechanism and concluded that the mean flow in the Caribbean Sea is intrinsically unstable, which means that any perturbation could trigger the formation of Caribbean anticyclones.

After formation, the Caribbean anticyclones are advected westward with the mean flow and propagate along the wind-driven upwelling regions at the South American coast. Based on a hydrographic time series of the upwelling in the Cariaco Basin, found that the interannual variability of temperatures in the upwelling region is affected by the anticyclones that advect cold filaments of upwelled waters offshore. This advection affects the ecosystem as it transports larvae and nutrients offshore . The advection of these cold filaments also leads to cooling of the interior of the Caribbean Sea . The anticyclones leave the Caribbean Sea through the Yucatán Channel. Model studies examining the Yucatán Channel region have shown that Caribbean anticyclones could influence eddy-shedding events of the Loop Current .

During their propagation, the anticyclones become more energetic . Although this intensification is clearly present in observations , only a few studies elaborate (briefly) upon the dynamics of this intensification. Based on surface drifter data, suggested that the anticyclonic shear of the Caribbean Current could amplify the anticyclones. In contrast, found, based on satellite altimetry data, a direct relationship between the maximum curl of the wind stress and the westward intensification of anticyclones. This relationship highlights that wind stress alters the life cycle of the anticyclones. used a regional model to study the life cycle of Caribbean anticyclones and computed the mechanical energy balance of the flow in this region. Although this balance shows that baroclinic instabilities provide the energy necessary for the westward intensification of the anticyclones, it does not explain what drives the westward intensification of the anticyclones.

In this study, we hypothesize that this intensification is steered by the offshore advection of cold upwelling filaments that cool the interior of the basin. The advection of the filaments results in denser surface waters in the western part of the basin compared to the eastern part of the basin, suggesting that the density difference between the relatively light anticyclones and the surrounding waters will increase. Following the thermal wind balance, it can be expected that the vertical shear of the anticyclone, and consequently its strength, will increase.

To test this hypothesis, we use a regional model in which we only vary the upwelling strength. To isolate the impact of upwelling on the life cycle of the anticyclones, we kept all other forcing parameters constant. The coastal upwelling is adjusted by altering the magnitude of the zonal wind stress with a constant. The curl of the wind stress is kept constant by applying the adjustment over the full model domain. A description of the model is provided in Sect. 2 and is followed by a comparison between the modeled flow and observations (Sect. 3). Sections 4 and 5 contain the analysis of the westward intensification of the anticyclone and how this is related to the advection of upwelling filaments. The sensitivity of the mean flow and eddy variability to changes in upwelling strength is discussed in Sect. 6.

2 Model configuration and methods

## 2.1 Model configuration

The numerical simulations were performed with the hydrostatic configuration of the Massachusetts Institute of Technology (MIT) primitive equation model . The computational domain extended from 99 to 55 W and from 6 to 33 N (Fig. 1a); it was set up with a horizontal resolution of 1/12, which is well below the internal Rossby radius of deformation in this region (60–80 km, Chelton et al.1998). The same topography as the Operational Mercator global ocean analysis (Mercator) of the EU Copernicus Marine Information Service was used wherein the topography is based on ETOPO1 for the deep ocean and near the coast and on GEBCO8 for the coast and slopes. In this topographic setup, the majority of the islands of the Lesser Antilles are represented (Fig. 1). The islands that are too small to be captured are modeled as shallow ridges. In the vertical, the model contained 50 levels in z coordinates, increasing in depth from 1 m at the surface towards 459 m at the lower levels.

The time step was 240 s. All simulations had a total duration of 25 years, the first 5 years of which were considered spin-up and were excluded from the analysis. The model output was saved as 5 d averaged fields. The simulations were initialized with time-averaged fields for sea-surface height, temperature, salinity, and velocity. These fields were obtained from the years 2007–2017 of Mercator. The vertical diffusion of tracers was parameterized with the GGL90 mixed layer parameterization . The horizontal diffusion was parameterized with the Redi scheme (diffusion coefficient of 125 m2 s−1Redi1982). The sub-grid-scale mixing was parameterized with Smagorinsky viscosities .

Figure 1Setup of the regional model in MITgcm. (a) Bathymetry in meters. The sponge layer is indicated with the dashed line. (b) Magnitude of the wind stress. The direction in indicated with the vectors. Not all velocity vectors are shown for clarity purposes. (c) Surface temperature field (C) used for restoration. (d) Surface freshwater flux (m s−1). Positive values correspond to a net surface buoyancy loss.

All simulations were forced with time-averaged surface and lateral boundary conditions. A sponge layer with a meridional width of 1.25 (15 grid cells) was applied at the lateral boundaries. There, the velocity fields were relaxed towards the time-averaged fields of Mercator (years 2007–2017) with a relaxation time that varied linearly from a month at the offshore side of the domain towards a day at the domain boundary. At these boundaries, temperature, salinity, and velocities in the zonal and meridional directions were prescribed. At the surface, a freshwater flux, temperature restoration, and a wind stress were applied (Fig. 1b–d). The temperature was restored towards sea-surface temperature (SST) averaged over the years 2007–2017 of Mercator with a relaxation timescale of 1 month. The surface freshwater flux was obtained from the diagnosed freshwater flux of a 250-year simulation described in . The Orinoco River, Magdalena River, and Mississippi River are prescribed as stationary freshwater fluxes at the open boundaries with discharges based on .

The reference simulation (Ekman100) is forced by a time-averaged zonal wind forcing computed from the wind fields of years 2007–2017 of ERA-Interim . Prescribing stationary values implied that the NBC Rings, which can trigger the formation of the anticyclones , were not represented at the boundaries. However, and showed that a realistic eddy field in the Caribbean Sea can be obtained without the presence of NBC Rings. Moreover, in Sect. 3 we will show that a realistic eddy field is obtained with these boundary conditions. Therefore, we considered these boundary conditions sufficient for the purpose of this study.

To investigate the effect of wind-driven upwelling on the westward intensification of anticyclones, only this zonal wind forcing (τx) was altered between simulations. The magnitude of the zonal wind stress was reduced or increased in each simulation by the same constant over the entire domain (Fig. 2). With this approach, we ensured that we only change the upwelling strength and not the curl of the wind stress. The magnitude of the reduction or increase was determined based on the zonal wind stress in the upwelling region in Ekman100. In Ekman150, the zonal wind stress in the upwelling region was 50 % stronger than the wind stress in Ekman100, resulting in a theoretical increase in the northward Ekman transport of 50 %. The wind stress along the coast was weaker in Ekman50, leading to a theoretically weaker upwelling (50 %) in this simulation compared to Ekman100. The same principle was applied in Ekman75 and Ekman125.

Figure 2Average zonal wind stress (N m−2) in the southern Caribbean Sea (10.5–12.5 N, 65–70 W) in the simulations. The seasonal variation of the zonal wind stress at this location, obtained from ERA-Interim , is shown for reference. Negative values correspond to easterly winds.

The adjustment of the zonal wind forcing resulted in changes in mesoscale variability over the full domain. Because the applied constant was optimized to alter the northward Ekman transport in the Cariaco and Guajira upwelling regions, it resulted in unrealistic magnitudes of the zonal wind stress (and thus of the northward Ekman transport) in some other regions. Therefore, we will only analyze the impact of the changes in wind forcing on the mesoscale variability close to the Cariaco and Guajira upwelling regions (Fig. 3c) and disregard other coastal upwelling regions.

The range of τx applied in the simulations covers the seasonal variations (Fig. 2). The magnitude of the stationary zonal wind stress in Ekman50 was similar to the observed zonal wind stress in fall, while the stronger stationary zonal wind stress in Ekman125 was comparable to the wind stress during winter and summer months. The zonal wind stress in Ekman150 was stronger than the seasonal cycle but was similar to zonal wind stresses that are observed during years with anomalously high wind velocities .

## 2.2 Methods

To analyze the behavior of the mesoscale eddies, we analyzed both the eddy kinetic energy (EKE) and the individual eddies. Following , we calculated the EKE from 5 d averaged velocity fields that were high-pass filtered with a 125 d running mean. This 125 d period was long enough to capture the variability of the eddies, which have a characteristic period of 50–100 d , but was shorter than the interannual variability.

To gain insight into the westward intensification of Caribbean anticyclones, we used the py-eddy tracker to follow them . The py-eddy tracker used 5 d averaged sea-level anomaly fields to identify near-circular features. Anomalies were computed with respect to 20-year averaged fields. Negative anomalies were identified as cyclones and positive anomalies as anticyclones. More detailed information about the numerics of the py-eddy tracker can be found in . We set a minimum lifetime of 75 d. Taking into account the average westward propagation velocity of 0.13 m s−1 , the anticyclones propagate approximately 840 km during this minimum lifetime.

At the center of the eddies provided by the eddy tracker, we extracted the amplitude (Aeddy), swirl velocity (uswirl), radius (Reddy), and properties from the model output to assess their characteristics. The amplitude was defined as the sea-surface height difference between the eddy and the 20-year average sea-surface height at the center of the eddy. Since the py-eddy tracker identifies circular anomalies as eddies, we could define the swirl velocity as the average of the maximum northward and maximum southward velocity of the eddy. The location of both velocities was used to obtain a robust estimate of the eddy radius (Reddy), which we defined as half the distance between these locations. To characterize the local properties (T, S, σ) of the background and the eddies, these variables were averaged over the upper 50 m. This depth corresponds to the average depth of the pycnocline in the Caribbean Sea. The latter ensures that these properties not only reflect variations in surface forcing. The background properties were defined as the 125 d averaged values of the upper 50 m at the location of the eddy obtained from the py-eddy tracker. The differences between the properties of the eddies and the background are computed as $\mathrm{\Delta }=\mathrm{eddy}-\mathrm{background}$.

To gain insight into the geostrophic part of the westward intensification of the anticyclones in each simulation, we computed the strength of the horizontal density gradients ($|\mathrm{\nabla }\mathit{\sigma }|$). These gradients were calculated over the upper 50 m of the water column for 5 d averaged density fields as follows:

$\begin{array}{}\text{(1)}& |\mathrm{\nabla }\mathit{\sigma }|=\sqrt{{\left(\frac{\partial \mathit{\sigma }}{\partial x}\right)}^{\mathrm{2}}+{\left(\frac{\partial \mathit{\sigma }}{\partial y}\right)}^{\mathrm{2}}},\end{array}$

where σ is the density. The contribution of temperature ($|\mathrm{\nabla }{\mathit{\sigma }}_{T}|$) and salinity ($|\mathrm{\nabla }{\mathit{\sigma }}_{S}|$) to the horizontal density gradients was computed in a similar manner, whereby the density differences were calculated assuming a linear equation of state as $\mathrm{\Delta }{\mathit{\sigma }}_{T}={\mathit{\rho }}_{\mathrm{0}}\mathit{\alpha }\left({T}_{\mathrm{1}}-{T}_{\mathrm{0}}\right)$ and $\mathrm{\Delta }{\mathit{\sigma }}_{S}={\mathit{\rho }}_{\mathrm{0}}\mathit{\beta }\left({S}_{\mathrm{1}}-{S}_{\mathrm{0}}\right)$, respectively. Here, α, β, and ρ0 are constants: $\mathit{\alpha }=-\mathrm{3.1}×{\mathrm{10}}^{-\mathrm{4}}$C−1, $\mathit{\beta }=\mathrm{7.2}×{\mathrm{10}}^{-\mathrm{4}}$ psu−1, and ρ0=999.8 kg m−3.

Since the focus of this study is to analyze the westward intensification of the anticyclones, we not only calculate the EKE and strength of the horizontal density gradients over the full domain, but we also calculate their contribution associated with the anticyclones only. We estimate the latter by considering the EKE and density gradients around the core of an eddy at each 5 d averaged field over a spatial extent of 1.5×Reddy. This procedure allows us to study only the westward intensification of the anticyclones.

3 Performance of the regional model

## 3.1 Mean flow

The modeled surface current in Ekman100 enters the Caribbean Sea through the southern passages of the Lesser Antilles (Fig. 3a), which is similar to what is seen in observations . The flow through the Lesser Antilles is highly variable and depends on variations of the flow upstream and the wind forcing . Based on scarce data, estimated the transport at 66 W as 18.4±4.7 Sv. In our simulation, the transport is more stable due to the stationary forcing and is on average 13.6 Sv at 66 W, which is close to the lower range of the estimate of . Further westward at 80 W, the modeled flow accelerates over shallow topography at 17 N, where it continues northwestward towards the Yucatán Channel into the Gulf of Mexico.

Figure 3The 20-year near-surface average properties of the Caribbean Sea in Ekman100. (a) Mean sea-level anomaly (SLA) in meters with velocity vectors. Surface current vectors are only shown every eighth grid cell for clarity purposes. (b) Density ($\mathit{\sigma }=\mathit{\rho }-\mathrm{1000}$; kg m−3), (c) temperature (C), and (d) salinity (psu). Properties in panels (b, c, d) are averages over the upper 50 m.

In the southwest of the Caribbean Sea, the modeled surface currents display a cyclonic recirculation (Fig. 3a). This recirculation is known as the Panama–Colombia Gyre (PCG) and is wind-driven . Part of the PCG continues as an eastward subsurface countercurrent, which results from the Sverdrup circulation in the North Atlantic tropical cell . The countercurrent flows below the upwelling region, and its depth is related to the wind strength . The model is able to reproduce the subsurface countercurrent at a depth of approximately 100 m, which is slightly higher in the water column than at the 200 m observed by and the 150 m observed by ; the modeled strength (0.11 m s−1) is comparable to their observations.

Ekman100 displays a strong meridional density gradient in the mixed layer that varies between σ=25.1 kg m−3 in the south (11 N) and σ=22.7 kg m−3 in the north (18 N, Fig. 3b). The strongest meridional gradients are close to the surface and colocated with the Caribbean Current. The location of the dense waters coincide with the two major upwelling regions: Cariaco and Guajira. The upwelled waters are colder and more saline than the surrounding surface waters (Fig. 3c, d).

The Cariaco upwelling region in the southeast of the Caribbean Sea is located at 63–65 W, 10–12.5 N . The modeled minimum (24.9 C) and maximum temperatures (27.4 C) in the Cariaco Basin are less extreme than the observed temperatures (minimum: 20.3 C, maximum: 30.1 C;  Rueda-Roa and Muller-Karger2013). Although the model is not able to capture this seasonal variability, the modeled average SST (25.4 C) is similar to observations (25.2 C; Rueda-Roa and Muller-Karger2013). Note that these modeled temperatures are lower than the restoring temperature (Fig. 1c), which indicates that the model is able to reproduce coastal upwelling. The Guajira upwelling region is located west of the Cariaco upwelling region between 69 and 74 W . Here, the observed average SST is slightly higher (25.5 C) than in the Cariaco upwelling region . Similar to these observations, the model displays a temperature difference between the two upwelling regions (26.1 C in Guajira; Fig. 3c). Furthermore, the modeled temperatures are in line with , who found surface temperatures varying between 25.6 and 28 C in the Guajira upwelling region.

## 3.2 Eddy kinetic energy

In line with observations , we find that the flow in the Caribbean Sea is highly variable (Fig. 4). In the eastern part of the basin, the modeled surface EKE is relatively low (100–300 cm2 s−2; Fig. 4a) and similar to observations of . The EKE increases westward towards a maximum >900 cm2 s−2 at 78 W. The modeled magnitude of EKE is higher than found in satellite altimetry (>600 cm2 s−2Andrade and Barton2000; Jouanno et al.2012), but it is more similar to estimates obtained from surface drifters (>900 cm2 s−2Richardson2005). This is in line with other modeling studies , and this discrepancy is mainly attributed to the coarse resolution (0.25) of the gridded altimetry data products . The modeled spatial variability of EKE matches analyses of satellite altimetry well . Corresponding to observations , we find that the eddy kinetic energy is surface intensified (Fig. 4b). In line with the modeling results of , the magnitude of EKE at depth also increases towards the west (Fig. 4b).

Figure 4Eddy kinetic energy (EKE; cm2 s−2) of Ekman100 obtained from 5 d averaged velocity fields that are high-passed filtered with a 125 d mean. (a) The 20-year average of the EKE in the upper 50 m of the water column and (b) a cross section of EKE at 15 N over the upper 400 m.

The variability in the Caribbean Sea manifests itself in the presence of mesoscale eddies (Fig. 5). In line with observations , the mesoscale eddies are predominantly anticyclonic (Fig. 5a). As expected, the surface densities of these anticyclones are lighter than those of the surrounding surface waters (Fig. 5b). The density differences are due to both temperature and salinity (Fig. 5c, d). Figure 5c also displays the northward advection of a cold filament. Similar filaments have been observed, and it is known that these can be advected several hundred kilometers from the upwelling region . This advection results in the offshore cooling of surface waters .

Figure 5Near-surface properties of the Caribbean Sea averaged over a 5 d period during the final year of the simulation Ekman100. (a) Mean sea-level anomaly (SLA, in meters) with velocity vectors. Surface current vectors are only shown every eighth grid cell for clarity purposes. (b) Density ($\mathit{\sigma }=\mathit{\rho }-\mathrm{1000}$; kg m−3), (c) temperature (C), and (d) salinity (psu). Properties averages over the upper 50 m.

## 3.3 Eddy characteristics

From the py-eddy tracker , we found that on average 8.55 anticyclones and 5.70 cyclones are formed in the reference simulation Ekman100 east of 75 W between 12.5 and 17.5 N every year. These numbers are in line with the 8–12 anticyclones per year estimated from surface drifter data . To our knowledge, there are no previous estimates for the formation rate of cyclones in the Caribbean Sea.

The anticyclones have an average amplitude of Aeddy=0.17 m and swirl velocities of uswirl=0.60 m s−1 between 65–75 W and 12.5–17.5 N (Table 1). These amplitude and swirl velocity values are similar to those found in hydrographic surveys and similar to velocities obtained from surface drifters (uswirl=0.5 m s−1Richardson2005). In general, the core of the anticyclones is warmer ($\mathrm{\Delta }T=+\mathrm{0.15}$C) and fresher ($\mathrm{\Delta }S=-\mathrm{0.18}$ psu) than the surrounding surface waters (Fig. 5c–d).

Table 1Number of eddies per year and average Rossby number with standard deviations in each simulation between 12.5–17.5 N and 65–75 W. Amounts are obtained with the py-eddy tracker. The formation rate is computed as the number of eddies that are formed each year east of 75 W. The long tracks represents the percentage of the eddies that were formed east of 65 W and could be tracked west of 75 W. The Rossby number is calculated as $Ro={u}_{\mathrm{swirl}}/\left({R}_{\mathrm{eddy}}*f\right)$, where uswirl is the swirl velocity of the eddy, Reddy is the radius, and f the Coriolis parameter. Density, temperature, and salinity are taken at the core of the anticyclones; 1 standard deviation is added to indicate the interannual variability.

In the central Caribbean Sea (65–75 W and 12.5—17.5 N), the cyclones are less energetic than the anticyclones and have an average amplitude of ${A}_{\mathrm{eddy}}=-\mathrm{0.16}$ m and swirl velocity of uswirl=0.50 m s−1. Observations indicate that cyclones are generated near topographic features in the Caribbean Sea and can contain upwelling waters . The modeled properties of the cyclones are in agreement with these observations. It appears that the mesoscale eddies are close to geostrophic balance, with an average Rossby number of 0.15±0.01 (anticyclones) and 0.14±0.03 (cyclones), where the Rossby number is calculated as $Ro={u}_{\mathrm{swirl}}/\left({R}_{\mathrm{eddy}}*f\right)$, in which uswirl is the swirl velocity of the eddy, Reddy is the radius, and f the Coriolis parameter.

There is a strong and significant correlation between the amplitude of the tracked mesoscale eddies and the surface EKE. This suggests that the westward increase in EKE is related to the strength of the eddies. We find that 57 of the anticyclones that are formed in the Caribbean Sea during the 20 years of the simulation propagate from a region with weak EKE (<65 W) towards a region with high EKE (>75 W). The paths of these anticyclones are indicated with the black lines in Fig. 6a. The other 96 anticyclones that are formed in the Caribbean Sea during the 20 years of the simulation are either generated west of 65 W or do not pass 75 W (grey lines in Fig. 6a). In contrast to the anticyclones, cyclones display relatively short tracks: only one cyclone passes both 65 and 75 W (black line in Fig. 6b). This is similar to the observations of , who showed that the cyclones follow a different path than anticyclones.

Figure 6Paths of (a) anticyclones and (b) cyclones that are identified by the py-eddy tracker in Ekman100 during 20 years of model simulation. The black lines indicate tracks that pass both cross sections at 65 and 75 W (green lines). All other tracks are indicated in grey. (c) Estimate of EKE (cm2 s−2) of the long tracks (black lines in panels a and b). (d) Estimate of EKE of all other tracks (grey lines in panels a and b).

To assess the contribution of the anticyclones with the long tracks to the total EKE variability, we calculated their EKE from the zonal and meridional velocity fields by taking into account the EKE within a region of 1.5×Reddy around each eddy as described in Sect. 2. The spatial distribution of EKE due to these long-lived eddies (Fig. 6c) is similar to the spatial distribution of the total EKE (Fig. 4a). Notably, the magnitude of the EKE of these selected eddies is approximately 55 % of the total EKE (Fig. 4), while it is computed from only 37 % of the eddies (57 out of 153 anticyclones in 20 years). This shows that the modeled westward increase in EKE is dominated by the westward increase in EKE of a small number of long-lived eddies, which is also confirmed by the weaker eddy kinetic energy of the eddies with shorter tracks (Fig. 6d). Because our results show that the eddies with long tracks dominate the mesoscale variability in the Caribbean Sea, we will focus on these 57 energetic anticyclones in the remainder of this study.

4 Westward intensification

Figure 6 shows that the EKE of the anticyclones with long tracks increases towards the west. We hypothesized that the intensification of the anticyclones was governed by a westward strengthening of the horizontal density gradients. In this section, we evaluate the evolution of the horizontal density gradients of the anticyclones and the validity of the thermal wind balance.

Figure 7 shows the different components of the thermal wind balance ($\frac{\partial V}{\partial z}$,$\frac{g}{{\mathit{\rho }}_{\mathrm{0}}f}\mathrm{\nabla }\mathit{\sigma }$) computed from high-pass-filtered density and velocity fields. As in Fig. 6c, only the anticyclones with long tracks are taken into account. The vertical shear of the anticyclones (Fig. 7a) increases from the eastern part of the basin towards the west. The horizontal density gradients of the anticyclones (Fig. 7b) have a similar magnitude and spatial distribution as the vertical shear (Fig. 7a). This indicates that, on average, the anticyclones were close to geostrophy. For this case, the horizontal density gradients due to salinity (Fig. 7d) are slightly stronger than the horizontal density gradients due to temperature (Fig. 7c).

Figure 7(a) Vertical shear of the anticyclones averaged over the upper 50 m. (b) Horizontal density gradients of the anticyclones, scaled with g∕(ρ0f) according to the thermal wind equation. Horizontal density gradients due to (c) temperature gradients and (d) salinity gradients of the anticyclones.

To gain insight into the westward intensification of the anticyclones, we computed the zonal variation of the EKE, vertical shear, density gradients, and properties of the anticyclones (Fig. 8). Overall, the meridional maximum EKE that is contained in the anticyclones increases from approximately 200 cm2 s−2 at 65 W towards 530 cm2 s−2 at 75 W (Fig. 8a). These values are similar to those observed in .

Similar to the EKE, the horizontal density gradients and vertical shear of the anticyclones increase towards the west (Fig. 8b). In line with the low Rossby number of the anticyclones, there is a small difference in magnitude between the vertical shear and density gradients, indicating the presence of small ageostrophic velocities. As in Fig. 7c and d, the horizontal density gradients due to temperature and salinity gradients both increase towards the west. This increase in both temperature and salinity can be explained by the properties of the upwelled waters, which are both colder and more saline than the background environment. Note that the density gradients induced by salinity are stronger than those induced by temperature differences (dashed and dotted lines in Fig. 7c–d, respectively). Although the anticyclones become slightly denser on their path westward through mixing with the surrounding waters (dots in Fig. 8c), this density increase is much weaker than the westward density increase of the surrounding waters (solid line in Fig. 8c), suggesting that changes in the properties of the background environment are indeed the dominant factor resulting in the westward increase in the horizontal density gradients.

Because the westward increase in the horizontal density gradients and vertical shear of the anticyclones (Fig. 8b) is not constant, we computed the variations with longitude in the westward direction of these quantities (Fig. 8d). From this, three regions can be identified as locations of more rapid intensification (64.6, 66.7, 72 W). Up to 64.6 W, the westward increase in horizontal density gradients is not fully balanced by the vertical shear, indicating that the westward intensification of the anticyclones is not in geostrophic balance at this stage. The westward intensification becomes more geostrophic towards the second peak of rapid intensification at 66.7 W (Fig. 8d). This peak is located near the Cariaco upwelling region. As the anticyclones move closer to the Guajira upwelling region, they intensify more rapidly again (at 72 W in Fig. 8d). Although the third rapid increase is located eastward of the Beata Ridge at 73 W, it is possible that this topographic feature has some impact on the westward intensification, as was previously proposed by . However, our model results suggest that all three regions of more rapid increase are located close to preferred locations of the shedding of upwelling filaments (Fig. 5c).

Figure 8(a) Zonal variation of the (a) meridional maximum eddy kinetic energy, as well as (b) vertical velocity shear and shear computed from horizontal density gradients between 12.5 and 17.5 N of the anticyclones with long tracks. The shaded regions indicate the 25th and 75th percentile. All properties are averaged over the upper 50 m of the water column. (c) Density anomaly of the background (σbg) and anticyclones (σAC) averaged over the upper 50 m. The density anomaly of the anticyclones is averaged per degree longitude. (d) The variations of vertical shear and horizontal density gradients with longitude in the westward direction of the anticyclones with long tracks.

5 Offshore advection of dense upwelling filaments

Figure 8 shows that the westward increase in EKE coincides with a westward increase in the vertical shear and density gradients of the anticyclones. This westward density increase of the surrounding surface waters is possibly governed by the advection of cold upwelling filaments offshore . According to , these filaments, characterized by strong temperature gradients, are advected by the anticyclones. Since we found regions with a more rapid increase in horizontal density gradients of the anticyclones (Fig. 8d), we evaluate in this section whether that increase is related to the advection of cold upwelling filaments.

We defined a measure for the strength of the filaments to study if there is a relation between the filaments and the eddy kinetic energy of the anticyclones. Because of their sense of rotation, the filaments are expected on the western side of the anticyclones, and hence we define this measure as a combination of the asymmetry of the density gradients around the anticyclone and the density gradient (induced by temperature) on the western side of the anticyclone. The asymmetry is defined as the ratio between the strength of the horizontal density gradient on the western and eastern side of the anticyclone ($\frac{|\mathrm{\nabla }{\mathit{\sigma }}_{\mathrm{west}}|}{|\mathrm{\nabla }{\mathit{\sigma }}_{\mathrm{east}}|}$). An asymmetry larger than 1 corresponds to a situation in which the density gradients on the western side of the anticyclone are stronger than on the eastern side, which suggests that an anticyclone advects a filament. Furthermore, we expect the anticyclones that advect filaments to have strong horizontal temperature gradients on the western side. These two properties are shown in Fig. 9a. It displays a positive correlation between the asymmetry and the strength of the temperature gradients so that the most asymmetric anticyclones have the strongest temperature gradients on their western side, corresponding to our hypothesis that the anticyclones with strong temperature gradients advect cold upwelling filaments.

Figure 9(a) Asymmetry of the anticyclones $\left(\frac{|\mathrm{\nabla }{\mathit{\sigma }}_{\mathrm{west}}|}{|\mathrm{\nabla }{\mathit{\sigma }}_{\mathrm{east}}|}\right)$ versus the strength of the horizontal temperature gradients on the western side of the anticyclones ($|\mathrm{\nabla }{\mathit{\sigma }}_{T,\mathrm{west}}|$) in Ekman100. The asymmetry is defined as the ratio between the strength of the horizontal density gradient on the western and eastern side of the anticyclone between 65 and 75 W. Colors indicate the EKE of the eastern side of each eddy. Each dot represents one anticyclone, and the inlay contains a schematic of the locations where the properties are calculated. (b) The total lifetime versus latitude of each anticyclone. Anticyclones with long tracks (black lines in Fig. 6a) are indicated with dots, and anticyclones with short tracks (grey lines in Fig. 6a) are indicated with triangles. Colors denote the asymmetry of the anticyclones. All properties are calculated as average values between 65 and 75 W from individual anticyclones in all simulations.

The relation between the eddy kinetic energy and the advection of filaments is also shown in Fig. 9a. The average eddy kinetic energy is computed on the eastern side of the anticyclone to exclude the kinetic energy associated with the filament. The EKE of the anticyclone is positively correlated with the strength of the western temperature gradient of the anticyclone, which indicates that the anticyclones with strong filaments are more energetic than the anticyclones with weaker filaments. This means that the westward intensification of individual anticyclones can be affected by the advection of cold filaments.

To understand this positive correlation between the advection of filaments and the strength of the anticyclones, we compared the anticyclones with long tracks to the other anticyclones that are present in the Caribbean Sea between 65 and 75 W (black and grey lines in Fig. 6, respectively). This comparison between the asymmetry of the anticyclones, their lifetime, and their meridional location reveals two interesting aspects (Fig. 9b). First, the anticyclones with long tracks (dots in Fig. 9b) are located in the southern part of the basin (<16 N). Second, the anticyclones with long tracks are in general more asymmetric than the other anticyclones (triangles in Fig. 9b).

The combination of these two aspects suggests a relation between the advection of filaments and the lifetime of Caribbean anticyclones: anticyclones that propagate close to the upwelling region are more likely to advect a cold filament that increases their western horizontal temperature gradients. Such offshore advection of dense filaments increases the horizontal density gradients, leading to the westward intensification of the anticyclones. Consequently, these anticyclones become more energetic and their lifetime increases. In contrast to the anticyclones with long tracks, the anticyclones with short tracks (triangles in Fig. 9b) are in general more symmetric and located farther north, and thus they apparently do not advect cold filaments. This suggests that anticyclones that do not advect cold filaments dissipate, deform, or merge rather than intensify. This is in line with the earlier result that the anticyclones with long tracks are more energetic than the anticyclones with short tracks (Fig. 6c, d).

6 Impacts of varying upwelling

In the previous section, we found that the westward intensification of Caribbean anticyclones is characterized by an increase in the horizontal density gradients of the anticyclones due to the advection of cold upwelling filaments by the anticyclones themselves. We expect changes in upwelling strength to result in changes in the properties of the filaments, which would affect the westward intensification of these anticyclones. To study the impacts of upwelling on this westward intensification of the anticyclones, we performed sensitivity simulations in which we altered the upwelling strength (Sect. 2). This resulted in differences in both the mean flow and mesoscale variability.

## 6.1 Changes in upwelling

The zonal wind stress in Ekman50 and Ekman75 was reduced compared to Ekman100 such that the wind stress in the upwelling region was 50 % and 25 % less than Ekman100, respectively. The decrease in zonal wind stress in Ekman50 and Ekman75 results in weaker upwelling and consequently warmer SST in the Cariaco Basin (Fig. 10b, c). The highest coastal temperatures are present in Ekman50. Sea-surface salinity decreased in both Ekman75 and Ekman50 compared to Ekman100 (Fig. 10b, c). This freshening is related to the presence of a subsurface salinity maximum in the Caribbean Sea due to the presence of a water mass referred to as Subtropical Underwater. This water mass is located at approximately 100 m of depth and leads to more saline upwelled waters compared to the fresher surface waters . Weaker upwelling thus results in warmer and fresher surface waters (Fig. 10b, c), and stronger upwelling results in colder and more saline surface waters (Fig. 10d, e).

Figure 10Cross section at 66 W of the 20-year averaged salinity in (a) Ekman100, (b) Ekman50, (c) Ekman75, (d) Ekman125, and (e) Ekman150 (shading). The solid red line indicates the position of the subsurface countercurrent (contour line at an eastward velocity of 0.07 m s−1). The Caribbean Current is indicated with the dashed red line (contour line at a westward velocity of 0.2 m s−1). The black contour shows the 20-year average temperature in each simulation; the white contour indicates the 20-year average density in each simulation.

Besides differences in the mean properties (σ, T, S), the width of the upwelling region also changes in the simulations: in Ekman125 and Ekman150, the upwelling is confined to a narrower region than in Ekman75 and Ekman50. Moreover, we find that the core of the subsurface countercurrent is displaced upwards in Ekman50 and Ekman75 compared to Ekman100 (red contour in Fig. 10). This upward change in line with observations in which the core was also displaced upwards in weaker wind conditions . The increased northward Ekman transport in Ekman125 and Ekman150 leads to colder and more saline surface waters. In these simulations, the core of the subsurface countercurrent is positioned lower in the water column (red contour in Fig. 10). The changing position of the subsurface countercurrent and the changes in sea-surface properties in the upwelling regions indicate that the model is able to reproduce realistic upwelling properties, as also shown in Sect. 3.

## 6.2 Changes in EKE

The upwelling strength impacts the offshore EKE (Fig. 11). In the interior of the Caribbean Sea, the surface EKE and subsurface EKE decrease (increase) for weaker (stronger) zonal wind strength. For example, the magnitude of the mean EKE in the interior of the Caribbean Sea decreases by 27 % in Ekman50 compared to Ekman100 (Fig. 11a). In Ekman75, a slightly weaker decrease of 15 % is present (Fig. 11b). In Ekman125 and Ekman150, the surface EKE increases by 13 % and 28 %, respectively (Fig. 11c, d). Overall, the EKE variations indicate a positive correlation between the coastal wind stress and the strength of the mesoscale eddies. A similar correlation is found in the subsurface EKE. This is in line with the finding that the upwelling was shallower in the simulations with weaker upwelling and deeper in simulations with stronger upwelling.

Figure 11EKE anomaly compared to Ekman100 (Fig. 4) at the near surface and at a cross section at 15 N in (a) Ekman50, (b) Ekman75, (c) Ekman125, and (d) Ekman150. Positive values correspond to enhanced EKE compared to Ekman100. The black box indicates the interior of the Caribbean Sea (12.5–17.5 N, 65–75 W). The dashed line shows the location of the cross sections in the lower panels.

To gain insight into the sensitivity of the westward increase in EKE to variations in upwelling, the maximum of the total EKE between 12.5 and 17.5 N was computed as a function of longitude (Fig. 12a). The largest changes in EKE are seen in the western part of the basin. In Ekman100, the EKE increased by 123 % from 65 to 75 W, while the weaker upwelling in Ekman50 and Ekman75 resulted in a smaller westward increase of the EKE by 62 % and 97 %, respectively. Ekman125 and Ekman150 have a stronger westward increase in EKE compared to Ekman100. Overall, there is a clear positive correlation between the upwelling strength and the westward increase in EKE.

Figure 12(a) Zonal variation of the maximum EKE between 12.5 and 17.5 N from 65 to 75 W for each simulation. (b) As (a), but for the anticyclones with long tracks only. The EKE was averaged over the upper 50 m. Note the difference in scale between the two panels.

## 6.3 Changes in eddy characteristics

Next, we assessed whether the changes in the westward increase in EKE can be attributed to changes in the eddy field and to changes in the westward intensification of anticyclones. To this end, we studied the behavior of individual eddies in each simulation with the py-eddy tracker. The number of anticyclones that are formed each year in the interior of the Caribbean Sea is highly variable and varies between 8.55 (Ekman100) and 10.85 per year (Ekman50, Table 1). These numbers are in line with a formation rate of 8–12 anticyclones per year as observed by . In all simulations, between 34 % and 43 % of the anticyclones that are formed in the Caribbean Sea have long tracks (Table 1). The paths of these anticyclones with long tracks are, similar to Ekman100 (Fig. 6a), located in the southern part of the basin (not shown).

Figure 13Zonal variation of the meridional average (12.5–17.5 N) of the (a) average vertical shear of the anticyclones, (b) horizontal density gradients, (c) horizontal density gradients induced by temperature gradients, and (d) horizontal density gradients induced by salinity gradients. All properties are computed over the upper 50 m of the water column.

Similar to Ekman100, we computed the EKE associated only with the long-lived anticyclones that were formed in the eastern part of the basin (<65 W) and propagate beyond 75 W (Fig. 12b). Similar behavior is visible as for the total EKE in Fig. 12a: the EKE increases towards the west in all simulations, and stronger (weaker) upwelling corresponds to a larger (smaller) westward increase. In the simulations with strong upwelling (Ekman100, Ekman125, and Ekman150), the anticyclones with long tracks (Fig. 12b) are responsible for more than half of the total EKE (Fig. 12a), even though they constitute only 34 %–43 % of the total number of anticyclones in this region (Table 1). In Ekman50 and Ekman75, the anticyclones with long tracks become less dominant. In these simulations, the westward increase in EKE is also less pronounced than in the stronger wind conditions. Overall, this shows that in all simulations, a substantial part of the EKE and the longitudinal variations therein are governed by the evolution of a small number of anticyclones and that this effect is stronger in simulations with stronger upwelling.

In general, the anticyclones are fresher and warmer than the surrounding waters in all simulations. The freshest and warmest anticyclones are present in Ekman50, while the most saline and coldest anticyclones are found in Ekman150 (Table 1). It is interesting to note that the variation of the properties of the anticyclones is very small in each simulation (see the standard deviation of the anticyclone properties in Table 1). This suggests that, like shown for Ekman100 in Fig. 8c, for the simulations in which the upwelling was varied the properties of the anticyclones also remain approximately constant during their propagation westward.

In line with the observations of , we found fewer cyclones than anticyclones in each simulation (Table 1). The lowest average formation rate of cyclones is found in Ekman150 with 4.15 cyclones per year, and the highest average formation rate is present in Ekman50 with 6.15 cyclones per year. These formation rates were highly variable, and differences in formation rates between simulations were not significant. In none of the simulations was the py-eddy tracker able to track multiple cyclones from east to west (65–75 W). This implies that the cyclones are either deformed or dissipated too much such that the py-eddy tracker could not track their sea-level anomaly. Overall, the behavior of the mesoscale eddies is similar in all simulations, and the spatial pattern and magnitude of the surface EKE is governed by the anticyclones with long tracks. The westward intensification of these anticyclones is discussed in the next section.

## 6.4 Westward intensification

We showed in Sect. 4 that the westward increase in EKE in Ekman100 coincided with a westward strengthening of the vertical shear of the long-lived anticyclones (Figs. 6c and 7a). This increase in vertical shear was related via the thermal wind relation to the strengthening of the horizontal density gradients at the western edge of the anticyclone caused by the advection of upwelling filaments (Fig. 9a). Here, we analyze the thermal wind relation in the sensitivity simulations and its relation to the westward intensification of the anticyclones (Fig. 13).

In all simulations, the vertical shear of the anticyclones increases from east to west (Fig. 13a). This westward increase in vertical shear is similar to the westward strengthening of the horizontal density gradients (Fig. 13b). Similar to Ekman100, the magnitude of the horizontal density gradients is slightly smaller than the vertical shear in all simulations, but their variation is similar. In combination with the low Rossby number of the anticyclones (Table 1), this implies that the anticyclones are in near-geostrophic balance during their evolution.

The anticyclones display the strongest westward intensification in the simulations with stronger upwelling (orange and red lines in Fig. 13a, b). These simulations have relatively low vertical shear at 65 W. In the simulations with weaker upwelling, Ekman50 and Ekman75, the anticyclones intensify less than in Ekman100. Moreover, in Ekman50, the vertical shear even weakens between 70 and 73 W. At the same location, the EKE also does not increase in Ekman50 (Fig. 12b).

To understand why the salinity gradients in the eastern part of the basin are stronger for weaker upwelling, we analyzed the average spatial variation of the total density gradients. Figure 14 shows the magnitude of the time-averaged horizontal density gradients $|\mathrm{\nabla }\mathit{\sigma }|$ of Ekman100 (middle column) and the strengthening or weakening of these gradients in Ekman50 (right column) and Ekman150 (left column) compared to Ekman100. The strongest horizontal density gradients are located at the Guajira and Cariaco upwelling regions (Fig. 14b). These density gradients are weaker in Ekman50 than in Ekman100 (Fig. 14a) and stronger in Ekman150 (Fig. 14c).

Figure 14(a–c) Average strength of the horizontal density gradients ($|\mathrm{\nabla }\mathit{\sigma }|$) in Ekman50, Ekman100, and Ekman150 (see Sect. 2). The middle column shows the magnitude of the density gradients in Ekman100. The left and right columns show the difference between Ekman50 and Ekman100 and between Ekman150 and Ekman100, respectively. Red indicates an increase in the horizontal density gradient. (d–f) The average strength of the horizontal density gradients due to temperature variations ($|\mathrm{\nabla }{\mathit{\sigma }}_{T}|$). (g–i) The average strength of horizontal density gradients due to salinity variations ($|\mathrm{\nabla }{\mathit{\sigma }}_{S}|$). The solid black line indicates where S=35.5 psu. (j–l) The relative magnitude of density gradients due to salinity compared to density gradients due to temperature ($\frac{|\mathrm{\nabla }{\mathit{\sigma }}_{S}|}{|\mathrm{\nabla }{\mathit{\sigma }}_{T}|+|\mathrm{\nabla }{\mathit{\sigma }}_{S}|}$). All properties are averaged over 20 years and over the upper 50 m of the water column.

The horizontal density gradients at the upwelling region are due to both horizontal temperature and salinity gradients (Fig. 14e–h). Because upwelling brings relatively cold and saline waters towards the surface (Fig. 10), the strength of these density gradients increases for stronger upwelling (Fig. 14c) and decreases for weaker upwelling (Fig. 14a). In Ekman150, the strengthening of the horizontal density gradients is mainly due to temperature-induced density gradients in the upwelling region (Fig. 14f, i). In contrast, in Ekman50, the weakening of the temperature gradients (Fig. 14d) and salinity gradients (Fig. 14g) contributes roughly equally to the weakening of the total horizontal density gradients (Fig. 14a).

In the interior of the basin, the horizontal density gradients strengthen in Ekman50 compared to Ekman100 (Fig. 14a), while they are weaker in Ekman150 compared to Ekman100 (Fig. 14c). It appears that this response results from the increasing importance of the Amazon and Orinoco River plumes for weaker wind conditions. The wind affects both the propagation and spreading of the river plume. First of all, the northward Ekman transport advects the river plumes towards the north as previously observed by . This northward Ekman transport is weaker in Ekman50, and consequently the river plumes follow the mean flow towards the west (Fig. 14g). The stronger Ekman transport in Ekman150 results in a more northward path of the river plumes (Fig. 14i). The second consequence of changing wind conditions is the impact of the wind-driven vertical mixing. The strong zonal winds in Ekman150 induce more wind-driven mixing of the surface waters than the winds in Ekman50. As a result, salinity anomalies propagate less far westward in Ekman150 than in Ekman50. The farther advection of the river plume in Ekman50 explains why the anticyclones have stronger salinity gradients in Ekman50 than in the other simulations (Fig. 13d).

From this analysis we conclude that a weakening of the zonal wind in Ekman50 leads to an increased influence of the Amazon and Orinoco River plume compared to the other simulations. More specifically, in Ekman50, the upper ocean density gradients in the interior of the Caribbean Sea are dominated by salinity gradients (Fig. 14j), while the density gradients are predominantly temperature-driven in Ekman150 (Fig. 14l). So, depending on the zonal wind strength, the total horizontal density gradients in the Caribbean Sea are either driven by temperature or salinity gradients, which suggests that both the river plumes and the upwelling affect the variability in the basin.

7 Summary and discussion

In this study, we used a regional model of the Caribbean Sea to investigate the interaction between mesoscale anticyclones and the wind-driven coastal upwelling along the South American coast. We showed that the westward intensification of Caribbean anticyclones is driven by an increase in the horizontal density gradient between their cores and their surroundings (Fig. 7). Notably, the increase is governed by density changes in the surroundings and not by density changes within the anticyclone (Table 1). More specifically, the density of the surroundings increases through the offshore advection of cold filaments of upwelled water by the anticyclones, which in turn is governed by the passage of the anticyclones themselves. This increases the horizontal density gradient between the anticyclone and the surroundings on the western side of the anticyclone (Fig. 9a). As a result, its vertical shear increases, and the anticyclone becomes more energetic. Approximately two to four anticyclones per year showed this behavior. The anticyclones that are not associated with the advection of cold filaments strengthen this view that the westward intensification of Caribbean anticyclones is facilitated by the advection of filaments of upwelled water. These anticyclones were less energetic and had shorter life spans than the anticyclones that intensified towards the west.

Further support for this view is obtained from the series of simulations in which the strength of the upwelling was altered by adjusting the zonal wind stress (and thus the northward Ekman transport). Stronger upwelling leads to an increase in offshore cooling and thus to a stronger increase in the horizontal density gradient between the anticyclones and their surroundings. As a result, the anticyclones intensify more on their path westward compared to the case of weaker upwelling. We also found that a decrease in the northward Ekman transport in the basin allows the river plumes to advect further into the basin. This influences the early development of the anticyclones. This nonlinear response of the evolution of anticyclones to changes in upwelling highlights the fact that the horizontal density gradients in the Caribbean Sea are predominantly set by salinity gradients in the case of weak wind forcing and by temperature gradients in the case of strong wind forcing.

Previous studies in the Caribbean Sea mainly focused on either the influence of the river plumes (e.g., Hu et al.2004; Chérubin and Richardson2007) or on the influence of the wind (e.g., Oey et al.2003; Astor et al.2003; Jouanno et al.2009). However, based on our results, we argue that it is important to take the interaction between the pathway of the river plume and the upwelling into account when studying mesoscale variability in the Caribbean Sea.

Based on the different wind forcing in the simulations, we can speculate on the seasonal eddy variability in the Caribbean Sea: the strong zonal wind stress in winter is similar to Ekman125, and the weaker winds in fall are represented by Ekman50 and Ekman75. In Ekman150, the anticyclones intensified more than in Ekman50. This implies that during the strong wind forcing in winter the anticyclones intensify more. Furthermore, during weak wind forcing, we expect the salinity gradients of the river plume to become dominant. However, the dispersal of the Amazon and Orinoco River plumes also has a distinct seasonal cycle , which is absent in our model and would be an interesting topic for further study.

The results of this study also highlight some aspects of the interannual variability of the eddy field in the Caribbean Sea. In our simulations, stronger wind forcing resulted in a higher EKE in the center of the Caribbean Sea. Assuming that this relation holds on interannual timescales as well, these results suggest a positive correlation between the wind forcing and eddy variability in the interior of the basin. used a model with seasonally varying boundary conditions and identified a similar relationship. This is also found in observations (Fig. 15), which show that sea-surface variance is higher in years with stronger zonal winds. Figure 15 also suggests that the response of the sea-surface variance to the wind stress is nonlinear: although 2010 and 2011 were years with weak zonal winds, the sea-surface variance was relatively high. It is interesting to note that the salinity in the Cariaco Basin was anomalously low in these years . Taking into account the shallow depth of this salinity anomaly, it is plausible that it is related to the farther westward propagation of the river plumes as seen in Ekman50. This supports our view that both upwelling and the dispersal of the river plumes affect the mesoscale eddy field in the Caribbean Sea.

Figure 15Year-averaged sea-level variance and coastal zonal wind velocities. Sea-surface spatial variance is obtained from daily sea-level anomaly fields from the years 1993–2017 from satellite altimetry (CMEMS) and averaged over the interior of the Caribbean Sea (12–17.5 N, 65–75 W). Coastal zonal winds in the upwelling region (10–12.5 N, 65–75 W) are obtained from the years 1993–2017 from the ECMWF (ERA-Interim; Dee et al.2011). Each dot represents one year.

Over the past 2 decades, a weak decreasing trend has been observed in the coastal wind stress and upwelling strength . Not only does this have severe ecological impacts , but it also impacts mesoscale variability. Our results suggest that weaker upwelling will lead to less westward intensification of Caribbean anticyclones and less mesoscale variability in the interior of the basin. Furthermore, in weaker wind conditions the offshore advection of (nutrient-rich) cold filaments is weaker, which results in less offshore advection of nutrients and biota.

Overall, in this study we showed how the strength of the zonal wind stress in the Caribbean Sea impacts the eddy variability in this basin. We showed how two processes induced by the northward Ekman transport impact the development of the most energetic anticyclones. First, the river plume sets the properties of the anticyclones during their early development in the eastern Caribbean Sea. Second, the anticyclones intensify themselves by the advection of cold and saline upwelling filaments. Together these two processes explain some key aspects of the mesoscale variability in the Caribbean Sea.

Code and data availability
Code and data availability.

For information about the regional configuration of MITgcm and about the processing code of all simulations, contact Carine van der Boog. The boundary conditions for the regional setting were obtained from the Operational Mercator global ocean analysis product provided by the Copernicus Marine Environment Monitoring System.

Author contributions
Author contributions.

CvdB designed the model and computational framework and analyzed the data. CK, JP, and HD supervised the work and helped shape the research, analysis, and paper. NB provided valuable support on the configuration of the model. All other authors discussed the results and contributed to the final paper.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

We would like to thank Sabine Rijnsburger for her suggestions to analyze the effect of the river plume. We would also like to thank Sotiria Georgiou, Adam Candy, Steffie Ypma, and Juan Manuel Sayol for valuable discussions on this work.

Financial support
Financial support.

This research has been supported by the Delft University of Technology (Delft Technology Fellowship grant) and the Netherlands Organisation for Scientific Research (NWO) (ALW-Caribbean grant with project 858.14.061).

Review statement
Review statement.

This paper was edited by A. J. George Nurser and reviewed by two anonymous referees.

References

Alvera-Azcárate, A., Barth, A., and Weisberg, R. H.: The surface circulation of the Caribbean Sea and the Gulf of Mexico as inferred from satellite altimetry, J. Phys. Oceanogr., 39, 640–657, https://doi.org/10.1175/2008JPO3765.1, 2009. a

Andrade, C. A. and Barton, E. D.: Eddy development and motion in the Caribbean Sea, J. Geophys. Res.-Oceans, 105, 26191–26201, https://doi.org/10.1029/2000JC000300, 2000. a, b, c, d, e, f, g, h, i, j, k, l

Andrade, C. A. and Barton, E. D.: The Guajira upwelling system, Cont. Shelf Res., 25, 1003–1022, https://doi.org/10.1016/j.csr.2004.12.012, 2005. a, b, c, d, e, f, g

Andrade, C. A., Barton E. D., Mooers, C. N. K.: Evidence for an eastward flow along the Central and South American Caribbean Coast, J. Geophys. Res., 108, 3185, https://doi.org/10.1029/2002JC001549, 2003. a, b, c, d, e

Astor, Y., Müller-Karger, F., and Scranton, M. I.: Seasonal and interannual variation in the hydrography of the Cariaco Basin:implications for basin ventilation, Cont. Shelf Res., 23, 125–144, https://doi.org/10.1016/S0278-4343(02)00130-9, 2003. a, b

Baums, I. B., Miller, M. W., and Hellberg, M. E.: Geographic variation in clonal structure in a reef-building Caribbean coral, Acropora palmata, Ecol. Monogr., 76, 503–519, https://doi.org/10.1890/0012-9615(2006)076[0503:GVICSI]2.0.CO;2, 2006. a

Beier, E., Bernal, G., Ruiz-Ochoa, M., and Barton, E. D.: Freshwater exchanges and surface salinity in the Colombian basin, Caribbean Sea, PLoS ONE, 12, 1–19, https://doi.org/10.1371/journal.pone.0182116, 2017. a, b

Campbell, J. D., Taylor, M. A., Stephenson, T. S., Watson, R. A., and Whyte, F. S.: Future climate of the Caribbean from a regional climate model, Int. J. Climatol., 31, 1866–1878, https://doi.org/10.1002/joc.2200, 2011. a

Candela, J.: Yucatan Channel flow: Observations versus CLIPPER ATL6 and MERCATOR PAM models, J. Geophys. Res., 108, 3385, https://doi.org/10.1029/2003JC001961, 2003. a

Cariaco project: Variation in the salinity between 1995 and 2015 at the CARIACO Time-series station, available at: http://imars.marine.usf.edu/cariaco/temperature-salinity-variation, last access: 9 May 2019. a

Carton, J. A. and Chao, Y.: Caribbean Sea eddies inferred from TOPEX/POSEIDON altimetry and a 1/6 Atlantic Ocean model simulation, J. Geophys. Res., 104, 7743, https://doi.org/10.1029/1998JC900081, 1999. a, b, c, d

Centurioni, L. R. and Niiler, P. P.: On the surface currents of the Caribbean Sea, Geophys. Res. Lett., 30, 10–13, https://doi.org/10.1029/2002GL016231, 2003. a, b, c, d

Chang, Y.-L. and Oey, L.-Y.: Coupled response of the trade Wind, SST gradient, and SST in the Caribbean Sea, and the potential impact on Loop Current's interannual variability*, J. Phys. Oceanogr., 43, 1325–1344, https://doi.org/10.1175/JPO-D-12-0183.1, 2013. a

Chelton, D. B., DeSzoeke, R. A., Schlax, M. G., El Naggar, K., and Siwertz, N.: Geographical variability of the first baroclinic Rossby radius of deformation, J. Phys. Oceanogr., 28, 433–460, https://doi.org/10.1175/1520-0485(1998)028<0433:GVOTFB>2.0.CO;2, 1998. a

Chérubin, L. M. and Richardson, P. L.: Caribbean current variability and the influence of the Amazon and Orinoco freshwater plumes, Deep-Sea Res. Part I, 54, 1451–1473, https://doi.org/10.1016/j.dsr.2007.04.021, 2007. a, b, c, d, e

Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., Berg, L. V. D., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., and Dee, D. P.: The ERA-Interim reanalysis: configuration and performance of the data assimilation system, Q. J. Roy. Meteorol. Soc., 137, 553–597, https://doi.org/10.1002/qj.828, 2011. a, b, c

Ducet, N., Le Traon, P. Y., and Reverdin, G.: Global high-resolution mapping of ocean circulation from TOPEX/Poseidon and ERS-1 and -2, J. Geophys. Res.-Oceans, 105, 19477–19498, https://doi.org/10.1029/2000JC900063, 2000. a

Fekete, B., Vorosmarty, C., and Grabs, W.: Global composite runoff fields based on observed discharge and simulated water balance, Rep. 22, Tech. rep., Global Runoff Data Cent, Koblenz, Germany, https://doi.org/10.1029/1999GB001254, 2000. a

Gaspar, P., Goris, Y. G. R. I., and Lefevre, J.-M.: A simple eddy kinetic energy model for simulations of the oceanic vertical mixing: Tests at station Papa and long-term upper ocean study site, J. Geophys. Res.-Oceans, 95, 179–193, https://doi.org/10.1029/JC095iC09p16179, 1990. a

Goni, G. J. and Johns, W. E.: Synoptic study of warm rings in the North Brazil Current retroflection region using satellite altimetry, Elsevier Oceanogr. Ser., 68, 335–356, https://doi.org/10.1016/S0422-9894(03)80153-8, 2003. a, b

Gordon, A. L.: Circulation of the Caribbean Sea, J. Geophys. Res., 72, 6207, https://doi.org/10.1029/JZ072i024p06207, 1967. a, b

Hellweger, F. L. and Gordon, A. L.: Tracing Amazon River water into the Caribbean Sea, J. Mar. Res., 60, 537–549, https://doi.org/10.1357/002224002762324202, 2002. a

Hernandez-Guerra, A. and Joyce, T. M.: Water masses and circulation in the surface layers of the Caribbean at 66 W, Geophys. Res. Lett., 27, 3497–3500, https://doi.org/10.1029/1999GL011230, 2000. a

Hu, C., Montgomery, E. T., Schmitt, R. W., and Muller-Karger, F. E.: The dispersal of the Amazon and Orinoco River water in the tropical Atlantic and Caribbean Sea: Observation from space and S-PALACE floats, Deep-Sea Res. Part II, 51, 1151–1171, https://doi.org/10.1016/j.dsr2.2004.04.001, 2004. a, b

Hughes, C. W., Williams, J., Hibbert, A., Boening, C., and Oram, J.: A Rossby whistle: A resonant basin mode observed in the Caribbean Sea, Geophys. Res. Lett., 43, 7036–7043, https://doi.org/10.1002/2016GL069573, 2016. a

Jochumsen, K., Rhein, M., Hüttl-Kabus, S., and Böning, C. W.: On the propagation and decay of North Brazil Current rings, J. Geophys. Res.-Oceans, 115, 1–19, https://doi.org/10.1029/2009JC006042, 2010. a

Johns, W. E., Townsend, T. L., Fratantoni, D. M., and Wilson, W. D.: On the Atlantic inflow to the Caribbean Sea, Deep-Sea Res. Part I, 49, 211–243, https://doi.org/10.1016/S0967-0637(01)00041-3, 2002. a, b, c, d, e

Jouanno, J. and Sheinbaum, J.: Heat balance and eddies in the Caribbean upwelling system, J. Phys. Oceanogr., 43, 1004–1014, https://doi.org/10.1175/JPO-D-12-0140.1, 2013. a, b, c, d, e, f

Jouanno, J., Sheinbaum, J., Barnier, B., Molines, J. M., Debreu, L., and Lemarié, F.: The mesoscale variability in the Caribbean Sea. Part I: Simulations and characteristics with an embedded model, Ocean Modell., 23, 82–101, https://doi.org/10.1016/j.ocemod.2008.04.002, 2008. a, b, c, d

Jouanno, J., Sheinbaum, J., Barnier, B., and Molines, J. M.: The mesoscale variability in the Caribbean Sea, Part II: Energy sources, Ocean Modell., 26, 226–239, https://doi.org/10.1016/j.ocemod.2008.10.006, 2009. a, b, c, d

Jouanno, J., Sheinbaum, J., Barnier, B., Marc Molines, J., and Candela, J.: Seasonal and interannual modulation of the eddy kinetic energy in the Caribbean Sea, J. Phys. Oceanogr., 42, 2015–2055, https://doi.org/10.1175/JPO-D-12-048.1, 2012. a, b, c, d

Le Bars, D., Viebahn, J. P., and Dijkstra, H. A.: A Southern Ocean mode of multidecadal variability, Geophys. Res. Lett., 43, 2102–2110, https://doi.org/10.1002/2016GL068177, 2016. a

Lems-de Jong, H.: Ocean current patterns and variability around Curacao for Ocean Thermal Energy Conversion, Master thesis, Delft University of Technology, 2017. a

Lima, F. P. and Wethey, D. S.: Three decades of high-resolution coastal sea surface temperatures reveal more than warming, Nat. Commun., 3, 1–13, https://doi.org/10.1038/ncomms1713, 2012. a

Lin, Y., Sheng, J., and Greatbatch, R. J.: A numerical study of the circulation and monthly-to-seasonal variability in the Caribbean Sea: The role of Caribbean eddies, Ocean Dynam., 62, 193–211, https://doi.org/10.1007/s10236-011-0498-0, 2012. a

Marshall, J., Adcroft, A., Hill, C., Perelman, L., and Heisey, C.: A finite-volume, incompressible Navier Stokes model for studies of the ocean on parallel computers, J. Geophys. Res.-Oceans, 102, 5753–5766, https://doi.org/10.1029/96JC02775, 1997. a

Mason, E., Pascual, A., and McWilliams, J. C.: A new sea surface height-based code for oceanic mesoscale eddy tracking, J. Atmos. Ocean. Technol., 31, 1181–1188, https://doi.org/10.1175/JTECH-D-14-00019.1, 2014. a, b, c, d

Molinari, R. L., Spillane, M., Brooks, I., Atwood, D., and Duckett, C.: Surface currents in the Caribbean Sea as deduced from Lagrangian observations, J. Geophys. Res., 86, 6537–6542, https://doi.org/10.1029/JC086iC07p06537, 1981. a, b

Molleri, G. S., Novo, E. M. M., and Kampel, M.: Space-time variability of the Amazon River plume based on satellite ocean color, Conti. Shelf Res., 30, 342–352, https://doi.org/10.1016/j.csr.2009.11.015, 2010. a

Murphy, S. J., Hurlburt, H. E., and O'Brien, J. J.: The connectivity of eddy variability in the Caribbean Sea, the Gulf of Mexico, and the Atlantic Ocean, J. Geophys. Res., 104, 1431, https://doi.org/10.1029/1998JC900010, 1999. a

Nystuen, J. A. and Andrade, C. A.: Tracking mesoscale ocean features in the Caribbean Sea using Geosat altimetry, J. Geophys. Res., 98, 8389–8394, https://doi.org/10.1029/93JC00125, 1993. a

Oey, L., Lee, H., and Schmitz Jr., W.: Effects of winds and Caribbean eddies on the frequency of Loop Current eddy shedding: A numerical model study, J. Geophys. Res., 108, 1–25, https://doi.org/10.1029/2002JC001698, 2003. a, b

Pauluhn, A. and Chao, Y.: Tracking eddies in the subtropical north-western atlantic ocean, Phys. Chem. Earth, Part A, 24, 415–421, https://doi.org/10.1016/S1464-1895(99)00052-6, 1999. a, b

Redi, M.: Oceanic isopycnal mixing by coordinate rotation, J. Phys. Oceanogr., 12, 1154–1158, 1982. a

Richardson, P. L.: Caribbean Current and eddies as observed by surface drifters, Deep-Sea Res. Part II, 52, 429–463, https://doi.org/10.1016/j.dsr2.2004.11.001, 2005. a, b, c, d, e, f, g, h, i, j, k, l, m

Rudzin, J., Shay, L., Jaimes, B., and Brewster, J. K.: Upper ocean observations in eastern Caribbean Sea reveal barrier layer within a warm core eddy, J. Geophys. Res.-Oceans, 119, 7123–7138, https://doi.org/10.1002/2016JC012339, 2017. a, b

Rueda-Roa, D. T. and Muller-Karger, F. E.: The southern Caribbean upwelling system: Sea surface temperature, wind forcing and chlorophyll concentration patterns, Deep Sea Res. Part I, 78, 102–114, https://doi.org/10.1016/j.dsr.2013.04.008, 2013. a, b, c, d, e

Silander, M. F.: On the three-dimensional structure of Caribbean mesoscale eddies, Master thesis, University of Puerto Rico, 2005. a, b, c

Simmons, H. L. and Nof, D.: The squeezing of eddies through gaps, J. Phys. Oceanogr., 32, 314–335, https://doi.org/10.1175/1520-0485(2002)032<0314:TSOETG>2.0.CO;2, 2002. a

Smagorinsky, J., Galperin, B., and Orszag, S.: Large eddy simulation of complex engineering and geophysical flows, Evolut. Phys. Oceanogr., 3–36, https://doi.org/10.1002/qj.49712052017,1993. a

Torres, R. R. and Tsimplis, M. N.: Sea-level trends and interannual variability in the Caribbean Sea, J. Geophys. Res.-Oceans, 118, 2934–2947, https://doi.org/10.1002/jgrc.20229, 2013. a

Van der Boog, C. G., de Jong, M. F., Scheidat, M., Leopold, M. F., Geelhoed, S. C. V., Schultz, K., Pietrzak, J. D., Dijkstra, H. A., and Katsman, C. A.: Hydrographic and biological survey of a surface-intensified anticyclonic eddy in the Caribbean Sea, J. Geophys. Res.-Oceans, 124, 6235–6251, https://doi.org/10.1029/2018JC014877, 2019.  a, b, c, d

Van Westen, R., Dijkstra, H., Klees, R., Riva, R., Slobbe, D., van der Boog, C., Katsman, C. A., Candy, A., Pietrzak, J., Zijlema, M., James, R., and Bouma, T.: Mechanisms of the 40–70 day variability in the Yucatan Channel volume transport, J. Geophys. Res.-Oceans, 123, 1286–1300, https://doi.org/10.1002/2017JC013580, 2018. a

Villamizar, E. Y. and Cervigón, F.: Variability and sustainability of the Southern Subarea of the Caribbean Sea large marine ecosystem, 22, 30–41, https://doi.org/10.1016/j.envdev.2017.02.005, 2017. a

Whyte, F. S., Taylor, M. A., Stephenson, T. S., and Campbell, J. D.: Features of the Caribbean low level jet, Int. J. Climatol., 28, 119–128, https://doi.org/10.1002/joc.1510, 2008. a