1. Introduction
The northern Colombian Caribbean has one of the most productive ocean systems, where oceanographic and meteorological processes are key to the understanding of the functioning of marine communities [1,22,23].
This region is influenced by the upwelling event, in which masses of cold deep waters enriched with nutrients used by phytoplankton emerge, increasing the sea’s productivity [2] and sustaining fisheries [4,39]. This event is present and frequent off the La Guajira Peninsula and near the Tayrona National Park [10,16,20,29].
In socioeconomic terms, activities such as ports, industry (e.g., mining, oil, coal), tourism and fishing converge in this area. Actions associated with these operations such as construction, ballast water, coal particles and wastewater, among others, may have an impact on marine environmental conditions and the water quality of the region [7,11,24,30]. Therefore, the evaluation of the dynamics of physicochemical variables facilitates the control and monitoring of water quality, guarantying the continuity of natural resources and of all trophic dynamics [6,42]. Additionally, it provides important information for environmental emergency responses. The characteristics of oceanographic conditions in the northern Colombian Caribbean have been widely described [2,4,8,20,28,29,32,37,38] but few investigations have been carried out on a wide spatial scale. Also, most studies have focused solely on the department of Magdalena [5,7,21,22,23,25,30] and the information has not recently been updated. Information on the dynamics of the most important physicochemical variables within the northern Colombian Caribbean is presented to a maximum depth of 150 m.
2. Materials and methods
2.1. Study area
The northern Colombian Caribbean coast covers the departments of Magdalena and La Guajira, between the mouth of the Magdalena River (11° 06’ 44.0’’ N; 74° 50’ 27.5’’ W) and Punta Gallinas (12° 33’ 00.7’’ N; 72° 13’ 20.1’’ W). The department of Magdalena has a narrow continental shelf, with rocky coasts, due to the mountainous relief of the Sierra Nevada de Santa Marta. On the other hand, La Guajira has a wide continental shelf of between 8 and 17 km, where a sandy bottom predominates [2,9,17,18].
The climate of the region is seasonal, and determined by the north-south movement of the Intertropical Convergence Zone (ITCZ) and the presence of trade winds from the NE [10,21]. Two climatic periods have been defined: (i) the dry season (December-April), in which the trade winds from the NE have the greatest influence, there is a decrease in water temperature (20-25 °C) and an increase in salinity (36) and wave action; and (ii) the rainy season (May-November) when trade winds weaken, Panama’s counter current dominates, water is warmer (29 °C), rains are frequent and continental discharges increase, mainly from the Magdalena River and from the Ciénaga Grande de Santa Marta (CGSM) [3,9,19,21,34].
The upwelling event varies throughout the year along with the wind velocity; both increase during the dry season [3,34].
2.2. Sample collection
In 2018, four scientific voyages were undertaken (May, August, November and December) and samples were taken from the 17 stations located in the water column between the mouth of the Magdalena River and Manaure (La Guajira) (Fig. 1). Places near the coast (4 km) at a maximum depth of 50 m and those further from the coast (8-10 km) at a maximum depth of 150 m were evaluated. Due to the fact that a sailing permit was denied, in August only five stations (E2A-E6A) located in the department of Magdalena were evaluated. Considering the oceanographic characteristics of the Caribbean Sea described by Chollett et al. [12], the study area includes two regions: 1) The first (Region 1) has heterogeneous, murky waters with low salinity and direct influence of rivers and the CGSM, corresponding to the stations located from the mouth of the Magdalena River to Santa Marta Bay (E1-E6A), and 2) the second (Region 2) is influenced by upwelling events, so the average sea temperature is colder, corresponding to the stations located from Tayrona National Park to Manaure (E7-E21).
At each station, salinity, dissolved oxygen, dissolved oxygen saturation and sea temperature were measured at every meter of depth, using a CTDO EXO1 (YSI) probe from the surface down to 150 m or 1 m from the bottom when water was shallower. Water transparency was determined using a Secchi disk. The pH was measured with a WTW 3210 probe, but due to a device malfunction, the data for this variable was not considered.
2.3. Analysis of the information
Hydrographic profiles of the variables measured were constructed for each voyage, showing the behavior of the water column. For water transparency, the light extinction coefficient k was calculated with equation (1):
where d is the depth of the Secchi disk [40]. The lower the value, the greater the light penetration. The characterization of the water mass was carried out by calculating measures of central tendency (arithmetic mean) and measures of variability (standard error, maximum, minimum, coefficient of variation). The principal component analysis (PCA) was performed with the software PAST v. 3.25, prior to standardization of the variable data [13]. Due to the different depths of the stations, this analysis used the average of the data down to the depth recorded at all of the stations each voyage: May down to 21 m, August down to 55 m and November and December down to 7 m. For each voyage, TS curves were made with the information from the stations located at a depth of more than 80 m. To determine the climate in which each sample was taken, the study used rainfall information and the ambient temperature measured by the 2018 IDEAM report from its coastal stations located in the Magdalena and La Guajira departments.
3. Results
The hydrographic profiles of the variables measured with the CTDO probe (Fig. 2-5) were taken into account in order to establish, with greater certainty, the climate in which each sample was carried out. In this way, samples from May and December were taken in the dry season. Samples from November and August were taken in the rainy season.
For the hydrographic profiles it is necessary to take into account that stations E15 and E18A, located in Region 2 in La Guajira, are the shallowest, with a depth of 7 and 10 m respectively. Table 1 summarizes the descriptive statistics of the variables evaluated by region during the study period. The temperature, dissolved oxygen and oxygen saturation values were higher on the surface and decreased slightly with depth. The salinity showed an opposite pattern, these values increased with depth until they reached a maximum between 50 and 100 m, and from there the gradient of salinity was reversed, decreasing slightly.
The temperature ranged between 18.9 and 28.5 °C, being colder during the dry season, especially in May. In general, in the water column the warmest layers were found in the first 40 m of depth and from 100 m down, colder waters below 22 °C (Fig. 2) were found. Spatially, Region 2 recorded the coldest surface waters, particularly between stations E16 and E21, whose values were lower than 26.1 °C. In Region 1 there was greater variability, with warmer surface waters in November.
Salinity varied between 34.8 and 38.2. In the first 10 m of depth the lowest values were recorded, particularly in Region 1, in August and November, with values between 34 and 36. Spatially, the values of surface salinity were lower in Region 1. In Region 2, the values were more stable (CV <0.5%) and were always higher than 36 (Fig. 3, Table 1).
The dissolved oxygen was between 4.4 and 6.6 mg / L. Below 80 m deep, the lowest values of oxygen concentration were found. Spatially, Region 2 registered the greatest dynamics of this variable over time (Table 1). In all voyages there were high saturation levels on the surface, with values ranging between 86.2 and 102.0%. The dynamics of this variable are similar to those observed for dissolved oxygen, in which subsatured waters predominated in both vertical and spatial distribution, except in May at stations E11 and E12 (Region 2) where supersaturated waters were detected (Figs. 4 and 5, Table 1).
The values of the light extinction coefficient k ranged from 0.04 m-1 in May (E2) to 0.97 m-1 in December (E11). The lowest values were recorded in the dry season and the highest in the rainy season (Fig. 6). In general, the stations closest to the river mouth such as E2A, E5, E10 and E11 registered higher values.
All of the TS curves presented a high salinity of 37.5, located between 80 and 90 m deep, between the surface and intermediate water layers. Likewise, the curves showed the presence of two bodies of water (Fig. 7). The first was found to a depth of close to 60 m, with temperatures between 27 and 29 °C and salinity between 34.8 and 37.0. The second body of water had a temperature of between 21 and 24 °C and maximum salinity values of close to 37.6. Temporary changes were observed in the first body of water; in the rainy season lower salinity was measured in a layer approximately 10 m deep.
According to the PCA, the first component (PC1) explained more than 48.4% of the variance, with temperature being the most characteristic variable in the first 20 m of depth, followed by salinity. In Fig. 8, station segregation of the first factor is seen, due to the variations in both variables. The stations were distributed into three broad groups, although this may change somewhat throughout the year according to the influence of upwelling and continental discharges: the first group was characterized by warmer waters, with temperatures higher than 26.9 °C, and were located closer to freshwater runoff; this applies to most of the stations in Region 1. The second group showed intermediate temperature and salinity values; stations from both regions were categorized in this group. The third group showed colder waters, with temperatures of less than 25.9 °C and with higher salinity (37.4). This group brought together stations located further east (Region 2) and further away from freshwater sources.
4. Discussion
Oceanographic conditions in the northern Colombian Caribbean are associated with the wind regime and the movement of the ITCZ [4, 10, 21]. In general, it was evidenced that the dynamics and hydrographic structure of the region were linked to upwelling events in the dry season and to the greater fluvial runoff due to heavy rainfall in the rainy season.
May and December showed features of the dry season, in response to less intense upwelling reflected in lower surface temperatures, higher salinities, clearer waters and low rainfall, especially in Region 1. Also, August and November presented features of the rainy season, with warmer surface waters and more frequent and intense rainfall that lead to a greater supply of fresh water (Fig. 2-6). All recorded values were similar to those found in other studies for the dry and rainy season [4,7,8,21,25,30,34,38].
The dynamics of the variables in the hydrographic profiles showed the typical structure of the tropical oceans, measuring larger values on the surface and decreasing slightly with depth, except for salinity, which showed an opposite pattern [27,39]. This behavior demonstrates the effect of depth on physicochemical variables, especially on temperature and dissolved oxygen.
The upwelling in the dry season usually implies a lower concentration of dissolved oxygen [9,21,23]. Contrary to this, in the dry season, high oxygen saturation values were recorded on the surface (88.9-102.0%), coinciding with values of higher transparency (Figs. 4, 6). According to Mancera et al. [30], this could be associated with an incursion of oceanic waters that are rich in oxygen. Similarly, the action of the winds could have a positive effect on the dissolved oxygen saturation. However, it must be considered that the concentration of this gas depends on several factors such as temperature, salinity, oxidation activities, photosynthesis, respiration, decomposition of organic matter and winds [21,33,35].
The results of the PCA (Fig. 8), showed a general spatial sectorization, grouping the western, central and eastern stations of the northern Colombian Caribbean together, as all presented dynamics that were regulated by variations in temperature and salinity, with temperature being most relevant. Similar results were obtained by Mancera-Pineda et al. [30] in Santa Marta bay, where both variables explained more than 50% of the data behavior. According to Knaus [27], these two variables can vary locally depending on continental runoff, evaporation and precipitation. The western stations (Region 1) recorded the strongest surface variations, particularly salinity (34-36) and higher temperatures of above 27 °C. The sea off the coast of the department of Magdalena is influenced by the discharges of the Magdalena River and the CGSM [29]. Therefore, its fluctuation is due to the proximity of these rivers’ contributions, which are greater during the rainy season [3,9,21]. The stations located in the east which are further away from the river runoff, especially between E16 and E21 (Region 2), recorded the coldest surface waters, which may indicate that a permanent upwelling condition occurs in this area. This condition coincides with other studies [4,38] which mention that colder surface waters in La Guajira can persist throughout the year.
The central sector, which brings together stations from both regions, recorded intermediate temperature and salinity values, since they receive waters with features recorded at the eastern and western stations. They are also influenced by river discharges from some of smaller rivers such as Manzanares, Mendihuaca Don Diego and Buritaca in Magdalena, and the Jerez, Tapias, Palomino and Ranchería rivers in La Guajira, whose flow rates are less than 0.06 L / s [36].
According to the TS curves, stations with depths greater than 80 m are constituted of two water masses (Fig. 7). The characteristics of the first water mass, only for the August and November assessments, coincide with the Caribbean Surface Waters (CSW) that have temperatures between 27 and 30 °C and salinities between 34.5 and 36.0 [4,26,29,31,32,41]. These waters with high temperatures and low salinities are the result of the seasonal displacement of the ITCZ, ocean-atmosphere events and the contribution of river runoff in the rainy season [8,29,31].
However, in the dry season, a surface mass which is slightly colder (26-27 °C) and considerably more saline (36.4-37.0) than the CSW was observed, which may be associated with a major upwelling of subsurface water mixed with surface waters. This noticeable variation in salinity was also found by Beier et al. [8] and Correa-Ramírez et al. [14] in this area. It should also be considered that in the dry season, the northern Colombian Caribbean presents colder and more saline waters in relation to the rest of the Colombian basin [8]. These results demonstrate the high variability of surface water bodies.
The second water mass has low temperatures (21-24 °C) and the maximum subsurface salinity value (37.6), which corresponds to subtropical subsurface waters (SSW). These waters have a core of maximum salinity (> 37) and derive from the waters of the Central Tropical Atlantic, where evaporation exceeds precipitation [26,31]. This body of water emerges in the upwelling events in the northern Colombian Caribbean [15], at which time the exchange of nutrients and even organisms [29] takes place.
5. Conclusions
In the northern Colombian Caribbean, the water features of the May and December 2018 voyages coincided with the dry season and those of August and November 2018 with the rainy season. Spatially, a general sectorization of the waters regulated by variations in temperature and salinity was evident. Therefore, the dynamics of the stations located further west (Region 1) were linked to their proximity to river discharges, while the variations of the stations farther east of Region 2 were linked to upwelling, which is persistent throughout the year.
In the stations with depths greater than 80 m, the presence of two types of water mass was identified: the first corresponds to the Caribbean Surface Waters, but with the influence of the upwelling in the dry season; and the second to subtropical subsurface waters. The behavior of physicochemical variables shows the great influence that the upwelling and continental discharges have on the hydrographic structure of the region.