Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Freshwater exchanges and surface salinity in the Colombian basin, Caribbean Sea

  • Emilio Beier ,

    Contributed equally to this work with: Emilio Beier, Gladys Bernal, Mauricio Ruiz-Ochoa, Eric Desmond Barton

    Roles Conceptualization, Formal analysis, Funding acquisition, Methodology, Supervision, Writing – original draft, Writing – review & editing

    Affiliation Centro de Investigación Científica y de Educación Superior de Ensenada–Unidad La Paz, La Paz, Baja California Sur, México

  • Gladys Bernal ,

    Contributed equally to this work with: Emilio Beier, Gladys Bernal, Mauricio Ruiz-Ochoa, Eric Desmond Barton

    Roles Conceptualization, Formal analysis, Funding acquisition, Methodology, Supervision, Writing – original draft, Writing – review & editing

    gbernal@unal.edu.co

    Current address: Departamento de Geociencias y Medio Ambiente, Universidad Nacional de Colombia, Medellín, Colombia.

    Affiliation Departamento de Geociencias y Medio Ambiente, Universidad Nacional de Colombia, Medellín, Colombia

  • Mauricio Ruiz-Ochoa ,

    Contributed equally to this work with: Emilio Beier, Gladys Bernal, Mauricio Ruiz-Ochoa, Eric Desmond Barton

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

    Affiliations Departamento de Geociencias y Medio Ambiente, Universidad Nacional de Colombia, Medellín, Colombia, Programa de Ingeniería Ambiental, Facultad de Ingeniería, Universidad Manuela Beltrán, Bucaramanga, Colombia

  • Eric Desmond Barton

    Contributed equally to this work with: Emilio Beier, Gladys Bernal, Mauricio Ruiz-Ochoa, Eric Desmond Barton

    Roles Conceptualization, Investigation, Methodology, Supervision, Writing – original draft, Writing – review & editing

    Affiliation Departamento de Oceanografía Instituto Investigaciones Marinas (CSIC), Vigo, España

Abstract

Despite the heavy regional rainfall and considerable discharge of many rivers into the Colombian Basin, there have been few detailed studies about the dilution of Caribbean Surface Water and the variability of salinity in the southwestern Caribbean. An analysis of the precipitation, evaporation and runoff in relation to the climate variability demonstrates that although the salt balance in the Colombian Basin overall is in equilibrium, the area south of 12°N is an important dilution sub-basin. In the southwest of the basin, in the region of the Panama-Colombia Gyre, Caribbean Sea Water is diluted by precipitation and runoff year round, while in the northeast, off La Guajira, its salinity increases from December to May by upwelling. At the interannual scale, continental runoff is related to El Niño Southern Oscillation, and precipitation and evaporation south of 12°N are related to the Caribbean Low Level Jet. During El Niño years the maximum salinification occurs in the dry season (December-February) while in La Niña years the maximum dilution (or freshening), reaching La Guajira Coastal Zone, occurs in the wet season (September-November).

Introduction

The Caribbean Sea is considered a concentration basin (i.e. evaporation exceeds precipitation) even though the increase of salinity is not enough to cause deep vertical convection [13]. Schmidt et al. [4] concluded that Caribbean salinity is not significantly affected by freshwater runoff and primarily reflects the evaporation/precipitation ratio. On the other hand, Yoo and Carton [2] noted that the river discharge and the annual average rate of advection of freshwater into the Caribbean Sea from the Atlantic differ by just over 10%.

The transport of water through the succession of Grenada, Venezuelan, Colombian, Cayman and Yucatan basins that constitute the Caribbean forms part of the western boundary current system of the North Atlantic Ocean and therefore is associated with the global meridional overturning circulation. Surface currents in the upper 50-to-100 m layer transport Caribbean Surface Water (CSW), which is formed by dilution of Atlantic Water by freshwater discharge mainly from the Amazon and Orinoco Rivers as it is transported into the Caribbean Sea by the Guiana Current [57]. Recently, Grodsky et al. [8] investigated the propagation of interannual changes in the Amazon plume as salinity anomalies through the central Caribbean on into the Gulf Stream system. CSW is transported by the Caribbean Current, which accelerates toward the Yucatan Channel, and has an average speed greater than 25 cm s-1 [911]. The Caribbean Current is made up of two cores that flow westward through the Venezuelan Basin to join in the Colombian Basin. As in all Caribbean sub-basins, the flow is strongly influenced by coherent eddy motions with horizontal scales of 100 km [7, 10, 1214]. Salinity anomalies propagate westward with these eddies, at an average speed of 11 cm s-1 [8]. In the southwest Caribbean, roughly south of 12°N, a recirculation takes place in the cyclonic Panama-Colombia Gyre as described by Mooers and Maul [15] and Richardson [10], among others. The southern limb of the Panama-Colombia Gyre constitutes the nearshore Panama-Colombia Countercurrent, which near 75°W mainly turns back with the general westward flow, though a proportion continues eastward along the continental slope as an undercurrent [12, 16].

The water flux budgets strongly affect the sea surface salinity, which in combination with temperature defines the surface water masses as well as the density field, which in turn influences ocean and coastal dynamics. Salinity is also relevant to biological processes, e.g. the life cycles of a number of organisms from phytoplankton to larger predators are affected by haline concentration and its variability. Up to now, there has been no detailed investigation of the salinity regime of CSW in the Colombian Basin (Fig 1). The region south of 12°N receives the discharge of the Magdalena River [1718], freshwater runoff from Central American rivers, and considerable precipitation over the ocean near Central America. The combined contribution of runoff and precipitation to the Colombian Basin results in strong surface salinity gradients between coastal and offshore waters [17, 1921]. However, the variability of these salinity contrasts has not been studied in any detail.

thumbnail
Fig 1. Limits of Colombian Basin with the location of the river discharges (dots).

The black line at 12°N marks the boundary of the southwestern Caribbean containing the Panama-Colombia Gyre. River outflows into the zone are marked by colored dots grouped into four regions: Central American Rivers (blue), Darien Gulf Zone (green), Magdalena River Zone (red) and La Guajira Coastal Zone (yellow).

https://doi.org/10.1371/journal.pone.0182116.g001

The surface water masses entering the Caribbean Sea through the passages and channels of the Greater and Lesser Antilles are transported by the Caribbean Current but may be modified by local processes including the freshwater inflow and coastal upwelling. In this study, we use data of evaporation, precipitation, and river discharge into the Colombian Basin and water mass analysis to show that the freshwater contributions and introduction of higher salinity water to the surface layers by upwelling modify CSW along the Colombian Coast.

Materials and methods

Precipitation over the basin was obtained from the monthly long-term mean data of the Climate Prediction Center Merged Analysis of Precipitation (http://www.esrl.noaa.gov/psd/data/gridded/data.cmap.html) with a resolution of 2.5° x 2.5° [22]. The data of Skin Temperature from the NCEP/NCAR Reanalysis with a spatial resolution of 1.86° x 1.86° (http://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.derived.surfaceflux.html) [23] and Net Latent Heat Flux from North American Regional Reanalysis, 0.25° x 0.25° (ftp://ftp.cdc.noaa.gov/Datasets/NARR/Monthlies/monolevel/apcp.mon.mean.nc) [23].

Evaporation during the common period of the three series (January 1979 to December 2012) was computed using the method described by Gill [24]: (1) where E is the evaporation rate (m s-1), Qe is the latent heat flux (W m-2), ρ is the density of water (kg m-3), Lv = 2.5008x106 –(2.3x103T) is the latent heat of evaporation (J kg-1), and T is the surface temperature (°C).

NCEP/NCAR and NARR reanalyses are among the most dependable, have a reasonably high resolution, and are widely used by the geophysical community. As a check on data quality, our evaporation series estimated from NCEP/NCAR data at (10.625°N, 79.375°W) was compared to that in the ERA-Interim reanalysis (http://www.ecmwf.int/en/research/climate-reanalysis/era-interim) for a location in the Caribbean representative of the area south of 12° (10.75°N, 78.25°W). The two series had similar characteristics over the period of comparison 1979–2014, with high correlation r = 0.93 and p-value = 0.001, which suggests good agreement between the two series. The NCEP/NCAR estimates are used for consistency with the other analyses in this paper.

Data from 36 rivers discharging into the Colombian Basin (Table 1) were obtained from different institutions of the countries around the basin. Long-term monthly records of 17 rivers were obtained from the Colombian Hydrology, Meteorology and Environmental Institute (Instituto de Hidrología, Meteorología y Estudios Ambientales de Colombia, IDEAM, https://www.sivirtual.gov.co/memoficha-tramite/-/tramite/T209, http://institucional.ideam.gov.co/jsp/info/institucional/solicitud/solicitud.pdf); and runoff climatology of rivers from Central America were obtained from the Panama Electrical Transmission Enterprise (Empresa de Transmisión Eléctrica S.A. de Panamá, ETESA, http://www.hidromet.com.pa/hidro_historicos.php), the 21th Hydrological Bulletin of the Costa Rica Electricity Institute (Instituto Costarricense de Electricidad, http://www.cfia.cr/descargas/informe1.pdf) [25] and the River Discharge Database Version 2.0 for Nicaragua (http://nelson.wisc.edu/sage/data-and-models/riverdata/index.php). Data in Table 1 were compared and completed with data of river discharges to the Colombian Basin including data from Restrepo et al. [18], Restrepo and Kjerfve [20] and Milliman and Farnsworth [26]. There are no reports of the estimation error of discharges. The data are separated into four groups in the table and figures for convenience of later reference.

thumbnail
Table 1. Annual average of the Rivers draining to the Colombian Basin.

https://doi.org/10.1371/journal.pone.0182116.t001

High-resolution (0.25° x 0.25°) gridded ocean climatology was obtained from the National Oceanographic Data Center Web site (http://www.nodc.noaa.gov/OC5/woa13/woa13data.html), hereafter WOA13v2 data [2728]. Among the few world databases of observed hydrography, this has the best resolution and is the most widely used. WOA13v2 hydrographic data provide objectively analyzed climatological mean fields, with greater vertical and spatial resolution than previously available [2729].

Climatology of temperature was derived from 1 to 1 32 casts per quadrant (grid square) and salinity from 1 to 37 casts. Though salinity observations in WOA13v2 are fewer than those of temperature, it is the best available source of information on in situ salinity. Observations were available for between 1 and 11 months of the year in each quadrant.

Monthly mean Sea Surface Salinity (SSS) from Aquarius/SAC-D observatory was downloaded with spatial resolution of 1°x1° [3031] from the Physical Oceanography Distributed Active Archive Center (PODAAC, http://podaac.jpl.nasa.gov/aquarius). We used Combined Active-Passive (CAPv3) Algorithm for the sea surface roughness correction to enable the retrieval of SSS, from September 2011 to April 2015 [31]. To compare the salinity data from WOA13v2 and Aquarius/SAC-D we linearly interpolated the 1-degree Aquarius data onto a 0.25-degree grid.

The climate of Colombia has been described by Mesa et al. [32]. The major dry season from December to February (DJF) is followed by a minor wet season from March to May (MAM), minor dry season (“midsummer drought”) from June to August (JJA), and major wet season from September to November (SON). We used the same seasons in our analysis. Seasonal total discharge in the rivers was calculated by regions (Table 1): Central America Rivers, Darien Gulf Zone with Atrato and Sinu as the most significant, Magdalena Rivers Zone including Dique Channel and Magdalena flow, and La Guajira Coastal Zone. As can be seen in Fig 1 and Table 1, river discharge in the Colombian Basin occurs mainly south of 12°N and can be represented by two rivers (Atrato and Magdalena) which combined contribute over 60% of the total runoff. Both rivers have long records and were used for further analysis.

For each forcing agent of SSS (evaporation, precipitation, and runoff) we have performed a seasonal fitting using Least Squares Fitting and estimations of uncertainties [3335]: (2) (3) where A0 is the mean value, Aα and As are the annual and semiannual amplitudes for each time series; ω = 2π/365.25 is the annual radian frequency; φα and φs are the phases of annual and semiannual harmonics referred to the beginning of the year; and t is the time. The term Fnonseason represents the anomalies and contains the non-seasonal variability (interannual and mesoscale). It was used to understand the interannual freshwater input into the basin, through comparisons with climatic indexes. Of the different El Niño Southern Oscillation (ENSO) indexes, the Oceanic El Niño Index appeared the most representative of the salinity forcing agents and is the one used here (ONI, http://www.cpc.ncep.noaa.gov/products/analysis_monitoring/ensostuff/ensoyears.shtml). We used the Caribbean Low Level Jet Index (CLLJ Index, Wang, 2007), defined by taking the negative of the 925hPa zonal wind anomalies in the region of 12.5°–17.5°N, 80o–70oW. Since the lower tropospheric winds in the region are easterly, the definition indicates that when the CLLJ Index is positive, the CLLJ is anomalously strong.

All time series used in this work are long enough (>25 years) to represent properly the interannual variability. To account for the effects of autocorrelation within each time series, the Effective Degrees of Freedom (Nef) were calculated following Trenberth [36].

To further investigate the forcing of interannual freshwater inputs in the Colombian Basin on water masses, we used information by seasons (three-monthly periods) of modeled SSS and currents near the coast during typical El Niño and La Niña years. These results correspond to outputs for the Colombian Basin over the period from 1979 to 1998 of a global run of the Parallel Oceanic Circulation Model (POCM-4C, courtesy of Dr. Robin Tokmakian from the Naval Postgraduate School, Monterey, CA, USA). The spatial and temporal resolution of the model was 0.25°x0.25° and one month, respectively. The POCM model [3738] reproduces properly in qualitative terms the seasonal and interannual variability of large scale circulation [3940].

Results and interpretation

Fresh water inputs

The climatological seasonal maps of evaporation minus precipitation (E-P) (mm month-1) are shown in Fig 2. From the central Colombian Basin toward the Mosquito and Darien Gulfs, (Fig 1) there is a gradient of E-P during the whole year. The high E-P values in the eastern, central half of the basin coincide with the core of the CLLJ. During the main dry season (DJF) evaporation dominates all over the basin from a maximum 180 mm month-1 in the central basin to a minimum of 40 mm month-1 in the Mosquito and Darien Gulfs. In MAM precipitation exceeds evaporation in the Mosquito Gulf. E-P values continue decreasing during the rest of the year (JJA and SON), with dilution in the southwest and concentration in the northeast. During the “midsummer drought” (JJA) the lowest values of E-P occur close to Central America, and during the main wet season (SON) the area of dilution extends as far as 76°W. The spatial pattern of dilution and concentration basins during the second half of the year reflects the Walker-type convection cell and the divergence of winds over the Colombian Basin [4143].

thumbnail
Fig 2. Evaporation minus precipitation (E-P) over the Colombian Basin by seasons.

Heavy contours indicate positive values, where evaporation exceeds precipitation (concentration conditions) and light contours indicate negative values (dilution conditions).

https://doi.org/10.1371/journal.pone.0182116.g002

The total annual mean river discharge into the Colombian Basin is estimated at 15397 m3 s-1 (Table 1). This is equivalent to about ~50% of the annual mean discharge of the Orinoco River (31000 m3 s-1) [7] into the eastern Caribbean. Over 60% of total water discharge directly into the Colombian basin arises from only two rivers, the Magdalena (7100 m3 s-1) and the Atrato (2400 m3 s-1). Rivers north of 12°N make a smaller contribution to the overall budgets. The climatological monthly mean of rivers draining to the Colombian Basin has a strongly seasonal signal, as can be seen in Fig 3 where the climatological values (thick lines) are very close to the seasonal fit (thin lines). Seasonal fitting explains from 93% to 99% of the climatological variance. Water discharges have a bimodal seasonal cycle (similar annual and semiannual amplitudes) with maximums in October–November and May–June, and minimums in February–March and July–August, consistent with the wet and dry seasons. The sum of all river runoff in the Colombian Basin has the minimum (maximum) runoff values in March (November), e.g., Magdalena River contributes around 10000 m3 s-1 in November and diminishes to 4000 m3 s-1 in March.

thumbnail
Fig 3. Seasonal cycle of river discharge into the Colombian Basin.

The regions as shown in Table 1 and Fig 1. The thin lines represents the seasonal fitting, and the respective percentage of explained variance is included.

https://doi.org/10.1371/journal.pone.0182116.g003

Taking into account the E-P and runoff shown in Figs 2 and 3, most of the surface water dilution in the Colombian Basin occurs south of 12°N. Climatological monthly means of evaporation (E), precipitation (P) and runoff (R), as well as the resulting budget between them, E-(P+R), are depicted in Fig 4A and 4C for the entire Colombian Basin and for the area south of 12°N Fig 4B, 4A and 4D, respectively.

thumbnail
Fig 4. Monthly means of evaporation, precipitation and runoff.

(a, b) and the resulting budget (c, d) over the entire Colombian Basin (a, c) and the area south of 12°N (b, d).

https://doi.org/10.1371/journal.pone.0182116.g004

Overall, in the Colombian Basin, precipitation and evaporation are up to 5 times higher than runoff (Fig 4A) while south of 12°N all are more similar in magnitude (Fig 4B). From May to November, the Colombian Basin is a dilution basin while from December to April it is a concentration basin (Fig 4C). The region south of 12°N suffers net dilution from April to December and only weak concentration in the remaining months (Fig 4D). Total evaporation in the Colombian basin is 57965 m3 s-1 and total precipitation 49184 m3 s-1.

The net evaporation for the Colombian Basin E-(P+R) is -6616 m3 s-1, an order of magnitude smaller than E, P or R. But the freshwater budget for the area south of 12°N is determined by evaporation of 12959 m3 s-1, precipitation of 17105 m3 s-1 and runoff at 11845 m3 s-1, and so the net balance of E-(P+R) = -15991 m3 s-1. Given that, the sum and its constituent terms all have the same order of magnitude; this sub-region appears to be a zone of net dilution. However, relative errors in each of the terms may be large, and more observations are needed to confirm this result.

Surface water properties

Using the WOA13v2 data in the Colombian Basin, the main Caribbean water masses reported in previous works [5, 15, 4446] can be identified. The mean potential temperature and salinity profiles as a function of depth, averaged over the whole Colombian Basin with their respective minimums and maximums (Fig 5), show a wide variation of the salinity of the CSW. Centered at 150 m, the Subtropical Underwater (SUW) is seen as a typical salinity maximum of 36.5 above the deeper Tropical Atlantic Central Water (TACW) and Antarctic Intermediate Water (AAIW) layers.

thumbnail
Fig 5. Mean vertical profiles of temperature (thick black) and salinity (thin blue) based on data from WOA13v2.

The lighter curves indicate range of the data, while horizontal lines indicate water mass limits.

https://doi.org/10.1371/journal.pone.0182116.g005

The spatial and seasonal changes can be seen in detail in the θ-S diagrams of surface waters down to 100 m depth (Fig 6) averaged by depths and regions as follows: north of 12°N; south of 12°N and east of 74°W (La Guajira upwelling zone); and south of 12°N and west of 74°W (the intense precipitation area of the Panama-Colombia Gyre). Surface waters are within the intervals reported for CSW (temperature between 25 and 29°C; salinity higher than 35.1) but, in the area of the Panama-Colombia region south of 12°N, waters above 30 m are fresher and warmer all year than those north of 12°N. Meanwhile, in La Guajira region, surface waters are significantly colder than in the other regions from December through May because of strong upwelling [4748]. During these months, a thin (~10 m) low salinity surface layer is present, but waters between 10 and 100 m are significantly saltier than in the other regions (also because of the upwelling). In JJA, surface water in la Guajira is colder but with the same salinity as the Panama–Colombia region (lower than the Central Caribbean). During SON, when upwelling is weaker, water masses of the three regions become similar, except for the shallowest surface waters (above 5 m) of Panama-Colombia and La Guajira whose salinity remains low.

thumbnail
Fig 6. Seasonal θ-S diagrams in the Colombian Basin for the upper 100 m layer.

Data from WOA13 hydrographic climatologic data, averaged by depths in three areas. (a) North of 12°N, (b) La Guajira (South of 12°N and East 74°W) and (c) Panamá-Colombia (South of 12°N and West 74°W).

https://doi.org/10.1371/journal.pone.0182116.g006

The seasonal distributions of WOA13v2 SST (Fig 7) show the strong surface cooling produced by the upwelling in La Guajira Coastal Zone particularly during DJF and MAM. During JJA and SON, although La Guajira region remains locally cooler, surface heating is evident throughout the basin, including the regions of anomalously high temperatures in the Darien and Mosquitos Gulfs (see Fig 1). These anomalies were reported and correlated with surface forcing by Ruiz-Ochoa et al. [41] based on Pathfinder satellite SST. The strong similarity of the temporal and spatial patterns from WOA13v2 and Pathfinder SST confirms the capability of the former to characterize the climatology of hydrography in the basin, although there are fewer salinity than temperature observations in the database.

thumbnail
Fig 7. Seasonal variability of sea surface temperature [°C] from WOA13 climatological data.

(a) DJF, (b) MAM, (c) JJA, and (d) SON.

https://doi.org/10.1371/journal.pone.0182116.g007

On the basis of WOA13v2 salinity surface quarterlies (Fig 8), the central basin salinities are generally highest in MAM and JJA, and lowest in SON, the period of most rain and advection of waters diluted by the Orinoco from the Atlantic [5, 11]. The cooler waters of La Guajira region are also generally saltier than the rest of the basin. In contrast with the saltier waters off La Guajira, fresher waters of salinities as low as 34.2 occur in the Darien and Mosquito Gulfs, coincident with higher temperatures. These low salinity waters extends eastwards along the Colombian coast in JJA reaching their maximum extension in SON, when precipitation and runoff are maximum. The eastward extension of low salinity waters is minimal in MAM.

thumbnail
Fig 8. Seasonal variability of sea surface salinity from WOA13 climatological data.

(a) DJF, (b) MAM, (c) JJA, and (d) SON.

https://doi.org/10.1371/journal.pone.0182116.g008

The spatial distributions of sea surface salinity and temperature suggest that in the Colombian Basin the CSW is diluted south of 12°N, where it acquires lower salinities and higher temperatures. Waters in the basin are more homogenous in MAM, after the main dry season, when they are saltier, and in SON, the main rainy season, when they are fresher. Off La Guajira, the colder and saltier water associated with the upwelling of SUW in DJF and MAM becomes fresher and warmer in SON because of the advection of waters from the west. From DJF through MAM, the intensity of the CLLJ and wind stress curl along the La Guajira Coast is most favorable to upwelling [41].

The recent availability of remotely sensed sea surface salinity from the Aquarius mission provides an independent check on the salinity field for the basin. Although Aquarius data are currently available only from September 2011 to April 2015 and are known to be unreliable close to shore, they still provide useful information on the region. We calculated seasonal means of these data using only pixels flagged as valid and free from land contamination. The seasonal composites over this period (2011–2015) showed higher salinity values (between 0.25 and 0.5) but a similar seasonal pattern to the WOA13v2 data (r = 0.72, p-value = 0.001) despite the different record lengths and the low spatial resolution of Aquarius.

Interannual variability of freshwater inputs and water masses

The seasonal cycles explain 61% of the total variance of E, 80% of P, and 53% and 30% of river outflow (R) from Magdalena and Atrato rivers, respectively. Fitting errors of amplitudes and phases are low. Monthly means of Magdalena runoff, Mnonseason, from 1950 to 2004, 648 observations, correlate with the ONI Index with an r = -0.52, p-value<10−10, Nef = 168 and a lag of one month. This suggests that during El Niño (La Niña) events, the Colombian Basin receives less (more) runoff. Furthermore, the Magdalena anomalies of discharge (Mnonseason) are correlated with the Niño1+2 Index, which represents the SST in the extreme Eastern Tropical Pacific, with r = -0.64, p-value<10−11, Nef = 172 and a lag of one month. Similar values of correlation coefficient and statistical significance were found between the Atrato river outflow and ONI index.

Monthly means of Pnonseason from 1980 to 2004, 288 observations, correlate well with the CLLJ, r = 0.55, p-value< 10−10, Nef = 160 with zero lag. Finally, Enonseason from 1979 to 2008, 348 observations, correlates with the CLLJ with r = -0.31, p-value< 10−7, Nef = 271 with no lag. Time series of anomalies of surface forcing of freshwater balance south of 12°N during the common period between 1980 and 2004 are shown in Fig 9, with El Niño and La Niña periods (notice that r, p-value, and Nef changes with time series length). For each time series Fig 9 also shows the relation with the interannual variability through the CLLJ Index (Fig 9A and 9B), and ONI (Fig 9C). The sum of P+R in Fig 9B and 9C is one order of magnitude higher than E. The figure shows that E and R decrease during El Niño, and increase during La Niña periods. The pattern for precipitation in relation with El Niño or La Niña is not clear, whenever it is not well correlated with ENSO.

thumbnail
Fig 9. Anomalies of evaporation and precipitation south of 12°N and water discharge of Magdalena River.

In dotted lines their best-correlated index (CLLJ Index and ONI, respectively). El Niño (dark gray) and La Niña (light gray) periods according to ONI are shaded. Red lines show 5-month moving averages of E and P.

https://doi.org/10.1371/journal.pone.0182116.g009

Output from the POCM-4C was used to test qualitatively if ENSO phases affect the river discharges sufficiently to modify the coastal surface water masses. Here Caribbean SSS maps were extracted from the global run over the period 1980 to 1998 as explained in the Methods section. Results for contrasting years show that during El Niño of 1983 (Fig 10) the low salinity coastal water in the Panama-Colombia Gyre was less extended and the high salinity off La Guajira was more extended than in the climatology of POCM, for both dry and wet seasons (Fig 11). Conversely, during La Niña of 1988 (Fig 10), the low salinity coastal water was more extended than in the climatology of POCM (Fig 11) and extended as far as La Guajira. Although the model captures the seasonality only qualitatively, it represents the general oceanography of the Colombian Basin and the spatial distribution of the water masses. Therefore, it serves as a first approximation to their interannual behaviour and represents the ENSO variability reasonably well.

thumbnail
Fig 10. Parallel Oceanic Circulation Model (POCM-4C) outputs of sea surface salinity (12.5 m).

(a) Wet season of El Niño 1982–83, (b) Dry season of El Niño 1982–83, (c) Wet season of La Niña 1988–89, and (d) Dry season of La Niña 1988–89.

https://doi.org/10.1371/journal.pone.0182116.g010

thumbnail
Fig 11. Parallel Oceanic Circulation Model (POCM-4C) outputs of sea surface salinity (12.5 m).

(a) January climatology, (b) October climatology, and (c) December climatology.

https://doi.org/10.1371/journal.pone.0182116.g011

Discussion and conclusions

The Caribbean Sea is generally considered a concentration basin where evaporation exceeds precipitation [14]. Remotely sensed SSS and numerical modelling indicates that, on an interannual basis, the salinity of water entering the Caribbean through the Lesser Antilles can vary by 0.5 psu [8]). This results from variability of the Amazon plume caused by changes in discharge rates and in the ocean and atmospheric forcing that govern its extent. Thus interannual anomalies of salinity unrelated to local conditions propagate westward through the central Caribbean and on into the Gulf Stream system. On the other hand, our results show that variability in the southwestern Caribbean is related to the local conditions. Unlike the rest of the Caribbean, the Colombian Basin is not a concentration basin, and south of 12°N it is a strong dilution basin in which the contribution of runoff to the water balance is as important as evaporation and precipitation (Fig 5D).

The net dilution in the southwest and upwelling in the northeast of the southern Colombian Basin modifies CSW, which shows wide variation of its salinity, as noted by Andrade and Barton [49]. Furthermore, Corredor [47], reported alternation of surface waters off La Guajira Peninsula between properties of the Panama-Colombia Countercurrent and wind-induced coastal upwelling. In situ observations by Andrade and Barton [48] and simulations by Lonin et al. [50] found that warm and low salinity waters from the Panama-Colombia Gyre block the southwestward extension of La Guajira coastal upwelling that normally reaches 10.5° N. During upwelling, we found a thin surface layer (0–10 m) of low-salinity water with very salty water between 10 and 100 m. Presumably when upwelling is weaker off La Guajira, the western water mass penetrates further to the east.

Interannual variability of the evaporation and precipitation in the south of the Colombia Basin correlates better with the local forcing of the CLLJ Index than with ENSO (ONI). This permanent low-level wind jet over the basin [5155] enhances evaporation, transports humidity and forms strong convection over topography in Central America. Conversely, runoff is better correlated with ENSO (ONI) because it originates principally in continental Colombia where, during El Niño events, precipitation and runoff decrease, while La Niña events have the contrary effects [19, 32, 5660]. Even though the annual cycles of precipitation and runoff are related, the deseasonalized time series are poorly correlated (r = 0.19, p-value < 0.03, Neff = 129 over 288 observations), because of the isolation of the catchment area of the Colombian rivers in the northern Andes from the Caribbean. Although both spatial patterns of precipitation over the Colombian mainland and their temporal variability are complicated, it has been recognized that ENSO is the greatest single cause of interannual variability within the region (e.g. Poveda et al. [57]) as indicated here by the significant correlation between the ONI and Magdalena River outflow.

Amador [54] reported that during El Niño (warm) phases, stronger than normal CLLJ winds occur, while during La Niña (cool) phases the CLLJ is weaker, in both cases especially during summer. A stronger CLLJ was found to be associated with greater than normal precipitation in the western Caribbean and weaker precipitation in the central Caribbean, while a weaker CLLJ produced the opposite. Our results for that portion of the Colombian basin south of 12°N show that precipitation is indeed positively correlated with the CLLJ Index, as expected. However, the precipitation anomalies were not well correlated with ENSO. We found that during the El Niño dry season (boreal summer) lower runoff results in less dilution of CSW in the Panama-Colombia Gyre and greater presence of upwelling waters (enhanced also by stronger winds of the CLLJ). In contrast, during the wet season of La Niña, runoff extends dilution beyond its normal range to reach as far as La Guajira Coastal Zone.

The interplay between warm and low salinity waters from the Panama-Colombia Gyre and cold and high salinity waters from La Guajira is related to the zonal wind circulation that plays a role in the Caribbean–Pacific connection [4143]. The CLLJ is the surface component of the cell with upward convective processes at the Caribbean side of Central America, and low-level divergence and downward drying over La Guajira. According to Andrade and Barton [42] and Hidalgo et al. [43], the convergent zone is connected with another cell in the west, with the sinking of dry air over the Eastern Tropical Pacific, where the drying favors the southward displacement of the Inter Tropical Convergence Zone (ITCZ). This could be the link between ENSO and CLLJ that produces the seasonal and interannual patterns of water masses observed here.

La Guajira upwelling is part of the South Caribbean Upwelling System [61]. It is characterized by strong physical forcing [41, 48, 6263], strong mixing of subsurface water with surface water, but low chlorophyll concentration [6364]. Different explanations reported for the low productivity in La Guajira include oligotrophic source waters [16], the upwelling of Panama-Colombia undercurrent water [63], a greater depth of the SUW in front of La Guajira than off Venezuela, or reduced efficiency of wind energy transfer [64]. Our findings indicate that from June to November, surface waters off La Guajira (see Figs 6 and 8) are similar to the western diluted CSW; especially during La Niña years (see Fig 10). This low salinity water, originating in the southwest of the basin, favors a strong surface stratification that inhibits vertical mixing and enhances surface heating and thus, plays a role in the low productivity off La Guajira.

Based on observed surface temperature, salinity and runoff, a first evaluation of the seasonal surface water balance of the Colombian Basin has been made. Results demonstrate that south of 12°N in the Colombian Basin, the surface water balance indicates a net dilution, responsible for producing a strong dilution of CSW. This water reaches the La Guajira region from June to November, where it dilutes saltier water produced by local upwelling.

At interannual scale, precipitation and evaporation are related to variability of the CLLJ, while runoff is related to ENSO. In turn, the CLLJ can change in relation to ENSO. Since south of 12°N runoff has the same order of magnitude as precipitation, the result is an anomalous decrease in the area of dilution and persistence of saltier water off La Guajira during the dry season of El Niño. In contrast, during the wet season of La Niña, the area of dilution extends beyond its normal range to reach La Guajira Coastal Zone. The intensified seasonal invasion of La Guajira upwelling by warm and freshwater from the Panama-Colombia Gyre during La Niña events favors surface stratification and can contribute to the low productivity of this upwelling system.

Supporting information

Acknowledgments

Thanks are due to IDEAM, ETESA, the Costa Rican Electricity Institute and the Nelson Institute for Environmental Studies of University of Wisconsin—Madison, for runoff data. To Dr. Robin Tokmakian from the Naval Postgraduate School, Monterey, CA, USA, for the results of POCM-4C model. And to NOAA, JPL Physical Oceanography Distributed Active Archive Center (PODAAC) and others mentioned in the text for the free data provided.

References

  1. 1. Etter PC, Lamb PJ, Portis DH. Heat and freshwater budgets of the Caribbean Sea with revised estimates for the Central American Seas. J Phys Oceanogr. 1987;17: 1232–1248.
  2. 2. Yoo JM, Carton JA. Annual and interannual variation of the fresh-water budget in the Tropical Atlantic-Ocean and the Caribbean Sea. J Phys Oceanogr. 1990;20(6): 831–845.
  3. 3. Tomczak M, Godfrey J. Regional Oceanography: An introduction. 2nd ed. New Delhi: Daya Publishing House; 2003.
  4. 4. Schmidt MW, Spero HJ, Lea DW. Links between salinity variation in the Caribbean and North Atlantic thermohaline circulation. Nature. 2004;428: 160–163. pmid:15014495
  5. 5. Kelly PS, Lwiza KMM, Cowen RK, Goni GJ. Low salinity pools at Barbados, West Indies: The origin, frequency, and variability. J Geophys Res Oceans. 2000;105(C8): 19699–19708.
  6. 6. Morrison JM, Smith OP. Geostrophic transport variability along the Aves Ridge in the eastern Caribbean Sea during 1985–1986. J Geophys Res Oceans. 1990;95(C1): 699–710.
  7. 7. Hu C, Montgomery ET, Schmitt RW, Muller-Karger FE. 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 2 Top Stud Oceanogr. 2004;51: 1151–1171.
  8. 8. Grodsky SA, Johnson BK, Carton JA, Bryan FO. Interannual Caribbean salinity in satellite data and model simulations. J Geophys Res Oceans. 2015;120: 1375–1387.
  9. 9. Fratantoni DM. North Atlantic surface circulation during the 1990’s observed with satellite-tracked drifters. J Geophys Res Oceans. 2001;106(C10): 22,067–22,093.
  10. 10. Richardson PL. Caribbean Current and eddies as observed by surface drifters. Deep Sea Res Part 2 Top Stud Oceanogr. 2005;52: 429–463.
  11. 11. Chérubin LM, Richardson PL. Caribbean Current variability and the influence of the Amazon and Orinoco freshwater plumes. Deep Sea Res Part 1 Oceanogr Res Pap. 2007; 54: 1451–1473.
  12. 12. Jouanno J, Sheinbaum J, Barnier B, Molines J-M, Debreu L, Lemarié F. The mesoscale variability in the Caribbean Sea. Part I: Simulations and characteristics with an embedded model. Ocean Model. 2008;23: 82–101.
  13. 13. Centurioni LR, Niiler PP. On the surface currents of the Caribbean Sea. Geophys Res Lett. 2003;30(6): 1279.
  14. 14. Alvera-Azcárate A, Barth A, Weisberg RH. The surface circulation of the Caribbean Sea and the Gulf of Mexico as inferred from satellite altimetry. J Phys Oceanogr. 2009;39: 640–657.
  15. 15. Mooers CNK, Maul GA. Intra-Americas Sea circulation. In: Robinson AR, Brink KH, editors. The Sea, vol. 11: Regional Studies and Syntheses. New York: John Wiley; 1998. pp. 183–208.
  16. 16. Andrade CA, Barton ED, Mooers CNK. Evidence for an eastward flow along the Central and South American Caribbean Coast. J Geophys Res Oceans. 2003;108(C6): 3185–3202.
  17. 17. Restrepo JD, Zapata P, Díaz J, Garzón-Ferreira J, García C. Fluvial fluxes into the Caribbean Sea and their impact on coastal ecosystems: The Magdalena River, Colombia. Global Planet. Change. 2006;50(1–2): 33–49.
  18. 18. Restrepo JC, Ortiz JC, Pierini J, Schrottke K, Maza M, Otero L, et al. Freshwater discharge into the Caribbean Sea from the rivers of Northwestern South America (Colombia): Magnitude, variability and recent changes. J Hydrol. 2014;509: 266–281.
  19. 19. Restrepo JD, Kjerfve B. Magdalena River: Interannual variability (1975–1995) and revised water discharge and sediment load estimates. J Hydrol. 2000;235(1–2): 137–149.
  20. 20. Restrepo JD, Kjerfve B. The Pacific and Caribbean Rivers of Colombia: Water discharge, sediment transport and dissolved loads. In: Lacerda L, Santelli R, Duursma E, Abrao J, editors. Environmental geochemistry in tropical and subtropical environments. Berlin: Springer Verlag; 2004. pp. 169–187.
  21. 21. Cañón M, Santamaría-del-Ángel E. Influencia de la pluma del río Magdalena en el Caribe colombiano. Bol Cient CIOH. 2003;21: 66–84.
  22. 22. Xie P, Arkin PA. Global precipitation: A 17-year monthly analysis based on gauge observations, satellite estimates, and numerical model outputs. Bull Amer Meteor Soc. 1997;78(11): 2539–2558.
  23. 23. Kalnay E, Kanamitsu M, Kistler R, Collins W, Deaven D, Gandin L, et al. The NCEP/NCAR 40-year reanalysis project. Bull Amer Meteor Soc. 1996;77: 437–471.
  24. 24. Gill A. Atmosphere-ocean dynamics. Vol. 30. San Diego: Academic press. International Geophysics Series; 1982.
  25. 25. Grant A, Oreamuno R, Serrano A, Vargas O. Comisión sobre la problemática de inundaciones en la vertiente Atlántica. Costa Rica: Colegio Federado de Ingenieros y de Arquitectos de Costa Rica; 2004.
  26. 26. Milliman JD, Farnsworth KL. River Discharge to the Coastal Ocean: A Global Synthesis. Cambridge: Cambridge University Press; 2011.
  27. 27. Boyer T, Locarnini RA, Zweng MM, Mishonov AV, Reagan JR, Antonov JI, et al. Changes to calculations of the World Ocean Atlas 2013 for version 2. Silver Spring: National Oceanographic Data Center; 2015.
  28. 28. Locarnini RA, Mishonov AV, Antonov JI, Boyer TP, Garcia HE, Baranova OK, et al. World Ocean Atlas 2013, Volume 1: Temperature. In: Levitus S, Mishonov A, editors. NOAA Atlas NESDIS 73. Silver Spring: National Oceanographic Data Center; 2013. pp. 1–40.
  29. 29. Zweng MM, Reagan JR, Antonov JI, Locarnini RA, Mishonov AV, Boyer TP, et al. World Ocean Atlas 2013, Volume 2: Salinity. In: Levitus S, Mishonov A, editors. NOAA Atlas NESDIS 74. Silver Spring: National Oceanographic Data Center; 2013. pp. 1–39.
  30. 30. Lagerloef G, Colomb R, Le Vine D, Wentz F, Yueh S, Ruf C, et al. The Aquarius/SAC-D mission: Designed to meet the salinity remote-sensing challenge. Oceanogr. 2008;21(1): 68–81.
  31. 31. Yueh S, Tang W, Fore A, Hayashi A. Aquarius CAP Algorithm and Data User Guide. Version: 3.0. Jet Propulsion Laboratory. California: California Institute of Technology; 2014.
  32. 32. Mesa O, Poveda G, Carvajal L. Introducción al clima de Colombia. 1era ed. Medellín: Universidad Nacional de Colombia; 1997.
  33. 33. Draper N, Smith H. Applied regression analysis. 3rd Ed. New York: John Wiley; 1981.
  34. 34. Beron-Vera FJ, Ripa P. Three-dimensional aspects of the seasonal heat balance in the Gulf of California. J Geophys Res Oceans. 2000;105(C5): 11,441–11,447.
  35. 35. Ripa P. Ajuste de datos por cuadrados mínimos. Cienc. Mar. 2002;28(1): 79–105.
  36. 36. Trenberth K. Signal versus noise in the Southern Oscillation. Mon. Weather Rev. 1984;112(2): 326–332.
  37. 37. Semtner AJ, Chervin RM. Ocean general-circulation from a global eddy-resolving model. J Geophys Res Oceans. 1992;97(C4): 5493–5550.
  38. 38. Tokmakian R. A high resolution Ocean model with variable forcing of wind, heat, and freshwater: Initial evaluation. Int. WOCE News. 1998;32: 26–28.
  39. 39. Matano R, Beier E, Strub P, Tokmakian R. Large scale forcing of the Agulhas variability: The seasonal cycle. J Phys Oceanogr. 2008;32(4): 1228–1241.
  40. 40. Matano R, Beier E, Strub P. The seasonal variability of the circulation in the South Indian Ocean: Model and observations. J Marine Syst. 2008;74(1–2): 315–328.
  41. 41. Ruiz-Ochoa M, Beier E, Bernal G, Barton ED. Sea surface temperature variability in the Colombian Basin, Caribbean Sea. Deep Sea Res Part 1 Oceanogr Res Pap. 2012;64: 43–53.
  42. 42. Andrade CA, Barton ED. Sobre la existencia de una celda de circulación atmosférica sobre el Caribe y su efecto en las corrientes de Ekman del Caribe suroccidental. Bol Cient CIOH. 2013;31: 73–94.
  43. 43. Hidalgo HG, Durán-Quesada AM, Amador JA, Alfaro EJ. The Caribbean Low-Level Jet, the Inter-Tropical Convergence Zone and precipitation patterns in the Intra-Americas Sea: A proposed dynamical mechanism. Geogr Ann Ser A. 2015;97: 41–59.
  44. 44. Morrison JM, Nowlin WD Jr. General distribution of water masses within the eastern Caribbean Sea during the winter of 1972 and fall of 1973. J Geophys Res Oceans. 1982;87(C6): 4207–4229.
  45. 45. Rhein M, Kirchner K, Mertens C, Steinfeldt R, Walter M, Fleischmann-Wischnath U. Transport of South Atlantic water through the passages south of Guadeloupe and across 16° N, 2000–2004. Deep Sea Res Part 1 Oceanogr Res Pap. 2005;52(12): 2234–2249.
  46. 46. Smith RH. Atlantic-Caribbean exchange through Windward Passage. M.Sc. Thesis, University of Miami. 2010. Available from: http://scholarlyrepository.miami.edu/oa_theses/24/.
  47. 47. Corredor JE. Apuntes sobre la circulación costera en el Caribe noroccidental colombiano. Bol Cient CIOH. 1981;3: 3–8.
  48. 48. Andrade CA, Barton ED. The Guajira upwelling system. Cont. Shelf Res. 2005;25: 1003–1022.
  49. 49. Andrade CA, Barton ED. The atmospheric Low Level Jet and the surface mesoscale circulation of the Caribbean Sea. GEOS. 2009;29(1): 85.
  50. 50. Lonin SA, Hernández JL, Palacios DM. Atmospheric events disrupting coastal upwelling in the southwestern Caribbean. J Geophys Res Oceans. 2010;115(C06030).
  51. 51. Wang C. Variability of the Caribbean Low-Level Jet and its relations to climate. Clim Dyn. 2007;29(4): 411–422.
  52. 52. Muñoz E, Busalacchi A, Nigam S, Ruiz-Barradas A. Winter and summer structure of the Caribbean Low-Level Jet. J Clim. 2008;21(6): 1260–1276.
  53. 53. Whyte FS, Taylor MA, Stephenson TS, Campbell JD. 2008. Features of the Caribbean Low Level Jet. Int J Climatol. 2008; 28(1): 119–128.
  54. 54. Amador J. The Intra-Americas Sea Low-Level Jet: Overview and Future Research. Ann NY Acad Sci. 2008;1146(1): 153–188. pmid:19076415
  55. 55. Cook KH, Vizy EK. Hydrodynamics of the Caribbean Low-Level Jet and Its Relationship to Precipitation. J Clim. 2010;23: 1477–1494.
  56. 56. Poveda G, Jaramillo A, Gil MM, Quiceno N, Mantilla RI. Seasonality in ENSO-related precipitation, river discharges, soil moisture, and vegetation index in Colombia. Water Resour Res. 2001;37(8): 2169–2178.
  57. 57. Poveda G, Waylen P, Pulwarty R. Annual and inter-annual variability of the present climate in northern South America and southern Mesoamerica. Palaeogeogr Palaeoclimatol Palaeoecol. 2006;234(1): 3–27.
  58. 58. Poveda G, Álvarez DM, Rueda OA. Hydro-climatic variability over the Andes of Colombia associated with ENSO: A review of climatic processes and their impact on one of the Earth’s most important biodiversity hotspots. Clim Dyn. 2010;36(11–12), 2233–2249.
  59. 59. Poveda G. La hidroclimatología de Colombia: una síntesis desde la escala inter-decadal hasta la escala diurna. Rev. Acad. Col. Cien. Exact. Fís. Nat. 2004;28(107): 201–222.
  60. 60. Tootle GA, Piechota TC, Gutiérrez F. The relationships between Pacific and Atlantic Ocean sea surface temperatures and Colombian streamflow variability. J Hydrol. 2008;349(3–4): 268–276.
  61. 61. Gordon A. Circulation of the Caribbean Sea. J Geophys Res Oceans. 1967;72(24): 6207–6223.
  62. 62. Petus C, Garcia-Valencia C, Thomas Y, Cesaracio M. 2007. Étude de la variabilitésaisonnière et interanuelle de la résurgence de la Guajira (Colombie) par analyse de donnéessatellitaires AMI-Wind, SeaWinds et AVHRR. Rev Télédétection. 2007;7(1-2-3-4): 143–156.
  63. 63. Paramo J, Correa M, Nuñez S. Evidencias de desacople físico-biológico en el sistema de surgencia en La Guajira, Caribe colombiano. Rev. Biol. Mar. Oceanogr. 2011;46(3): 421–430.
  64. 64. Rueda-Roa DT, Muller-Karger FE. The southern Caribbean upwelling system: Sea surface temperature, wind forcing and chlorophyll concentration patterns. Deep Sea Res Part 1 Oceanogr Res Pap. 2013;78: 102–114.