ORIGINAL RESEARCH article
Horizontal and Vertical Distributions of Transparent Exopolymer Particles (TEP) in the NW Mediterranean Sea Are Linked to Chlorophyll a and O2 Variability
- Biologia Marina i Oceanografia, Consejo Superior de Investigaciones Científicas, Institut de Ciències del Mar, Barcelona, Spain
Transparent Exopolymer Particles (TEP) are relevant in particle and carbon fluxes in the ocean, and have economic impact in the desalination industry affecting reverse osmosis membrane fouling. However, general models of their occurrence and dynamics are not yet possible because of the poorly known co-variations with other physical and biological variables. Here, we describe TEP distributions in the NW Mediterranean Sea during late spring 2012, along perpendicular and parallel transects to the Catalan coast. The stations in the parallel transect were sampled at the surface, while the stations in the perpendicular transect were sampled from the surface to the bathypelagic, including the bottom nepheloid layers. We also followed the short-term TEP dynamics along a 2-day cycle in offshore waters. TEP concentrations in the area ranged from 4.9 to 122.8 and averaged 31.4 ± 12.0 μg XG eq L−1. The distribution of TEP measured in transects parallel to the Catalan Coast correlated those of chlorophyll a (Chla) in May but not in June, when higher TEP-values with respect to Chla were observed. TEP horizontal variability in epipelagic waters from the coast to the open sea also correlated to that of Chla, O2 (that we interpret as a proxy of primary production) and bacterial production (BP). In contrast, the TEP vertical distributions in epipelagic waters were uncoupled from those of Chla, as TEP maxima were located above the deep chlorophyll maxima. The vertical distribution of TEP in the epipelagic zone was correlated with O2 and BP, suggesting combined phytoplankton (through primary production) and bacterial (through carbon reprocessing) TEP sources. However, no clear temporal patterns arose during the 2-day cycle. In meso- and bathypelagic waters, where phytoplanktonic sources are minor, TEP concentrations (10.1 ± 4.3 μg XG eq l−1) were half those in the epipelagic, but we observed relative TEP increments coinciding with the presence of nepheloid layers. These TEP increases were not paralleled by increases in particulate organic carbon, indicating that TEP are likely to act as aggregating agents of the mostly inorganic particles present in these bottom nepheloid layers.
Transparent Exopolymer Particles (TEP) are defined as a subclass of gel-like organic particles, mainly composed by acidic polysaccharides, that are stainable with Alcian Blue (Alldredge et al., 1993). These particles are widespread in aquatic ecosystems, and their study in the ocean has biogeochemical and applied interests. Due to their high stickiness, TEP act as gluing agents for other particles to form larger aggregates susceptible to sink in the water column, hence stimulating the biological carbon pump (Passow et al., 2001; Burd and Jackson, 2009). However, TEP themselves have low density, and when unballasted, they can ascend through the water column (Azetsu-Scott and Passow, 2004) and accumulate in the sea surface microlayer (Wurl et al., 2011) where they can constitute a major source of primary aerosols (Orellana et al., 2011). The study of TEP has also gained interest in the water desalination industry since they are major agents of reverse osmosis membrane fouling (Berman, 2013). Given the ecological and economic relevance of TEP, there is a need to improve the knowledge about how these substances are distributed in the field and what factors affect their dynamics.
TEP were first observed, and have mostly been described, associated with phytoplankton blooms, in the field (Alldredge et al., 1993; Van Oostende et al., 2012) or in mesocosms and controlled chambers (Engel et al., 2015). From these studies we know that phytoplankton are a major source of TEP and TEP precursors in the sea. However, the relationship between phytoplankton and TEP varies depending on phytoplankton composition and physiology (Passow, 2002; Klein et al., 2011), and environmental variables such as nutrient availability (Mari et al., 2001), turbulence (Pedrotti et al., 2010), or UV irradiation (Ortega-Retuerta et al., 2009a). Therefore, even when phytoplankton are the likely main source for TEP, this is not necessarily translated into predictable relationships in the field between TEP and chlorophyll a (Chla), the most used proxy for phytoplankton biomass or production. Elucidating the sources of variability in the TEP-Chla relationships would help predicting the occurrence and dynamics of TEP in the ocean.
In addition to the role of phytoplankton, there are other sources of TEP in the sea, such as macroalgae (Thornton, 2004) or zooplankton (Prieto et al., 2001). Also bacteria are known to modify TEP distributions in the sea in various ways: They colonize and degrade TEP that are released by other organisms (Bar-Zeev et al., 2011; Taylor et al., 2014), thus acting as TEP sinks. Bacteria can also directly release TEP (Ortega-Retuerta et al., 2010) so they constitute TEP sources themselves. Finally, bacterial interactions with phytoplankton mediate TEP release (Van Oostende et al., 2013) and induce changes on their formation rates and properties such as their stickiness (Rochelle-Newall et al., 2010). The relative importance of these mechanisms governing the TEP dynamics in aquatic habitats, specifically in the Mediterranean Sea, remains unexplored.
The published information on TEP distributions in the Mediterranean Sea is particularly scarce (Prieto et al., 2006; Ortega-Retuerta et al., 2010; Bar-Zeev et al., 2011). The few published studies, however, concur in that TEP stocks are high when compared to other oceans. For instance, maximum TEP concentrations (up to 11,000 μg Xeq. L−1 in surface waters) were observed in Adriatic Sea samples (Passow, 2002). Exceptionally high in the Mediterranean Sea are the relative TEP concentrations with respect to Chla concentrations; higher TEP/Chla ratios than in other ocean basins have been taken to suggest that TEP are an important fraction of the particulate organic matter pool, and likely important drivers of carbon and particle fluxes in this oligotrophic sea.
Here, we report for the first time TEP distributions in the Catalan Sea (NW Mediterranean). Our specific goals were: (1) to determine the potential drivers of TEP from a wide range of physicochemical and biological variables and (2) to examine the variability in the TEP-Chla relationship across multiple spatial and temporal scales.
Materials and Methods
Study Site and Sampling
Samples were taken during the cruises NEMO1, NEMO2, and SUMMER2 in Mediterranean waters between the Catalan Coast and north of Majorca Island on board the Spanish RV “García del Cid” (Figure 1). Transects parallel to the Catalan coast (following the bathymetry line at 40 m bottom depth) between Barcelona and Blanes were conducted in May 10th (transect 1) and June 11th (transect 2, Figure 1). During these transects, surface (2 m) samples were taken every hour from the underway continuous flow with the ship moving at ~7 knots, so that each sample was taken approximately at every 12 km. A coast-to-offshore transect was performed during NEMO1, from May 11th to 20th, including stations located in the shelf (Stations 1 and 2), slope (stations 3, 4, and 5) and basin (stations 6, 8, and 9). Station 7 was sampled during NEMO2, 1 month later (June 12th). Water samples in these transects were collected using a rosette (12 Niskin bottles with external spring, 12 L each) coupled to a Sea-Bird Conductivity-Temperature-Depth profiler, a WET Labs C-Star transmissometer and a SeaPoint optical backscatter sensor. Up to six depths were sampled from each station, from surface to bottom (down to 2300 m) waters including the surface, the O2 maximum, the deep chlorophyll maximum (DCM) when present, mesopelagic waters, and bottom nepheloid layers.
Figure 1. Study area. Blue symbols, Coastal transects. symbols, Coast-to-open sea transect. Green symbol, 2-day cycle.
The 2-day lagrangian study (SUMMER2 cruise) was conducted aboard R/V “García del Cid” from 22nd to 24th May at ca. 45 nautical miles from the coast, within the core of a cyclonic eddy over a water-column depth of ca. 2000 m. A Lagrangian drifter was deployed to track the movement of the upper 15-m water layer. Each drifter consisted of a spherical floatable enclosure that contained a GPS and an emitter, from which 10 m cylindrical drogues hanged 5 m below the sphere. The drifters sent their position every 30 min, and all ship operations were conducted next to them. Samples were taken with the rosette every 4 h at six depths from surface to 200 m.
Chemical and Biological Analyses
TEP were analyzed following the colorimetric method proposed by Passow and Alldredge (1995). Samples (250–500 mL) were filtered through 25 mm diameter 0.4 μm pore size Polycarbonate filters (DHI) at low pressure (100 mm Hg). The filters were stained with 500 μL of Alcian Blue (0.02%, pH 2.5) for 5 s and rinsed with MilliQ water. The filters were soaked in 80% sulfuric acid for 3 h and the absorbance of the extract was determined at 787 nm in a Varian Cary spectrophotometer. Duplicates were taken for each sample. Previous analyses have shown a CV of 13% between TEP replicated measurements with this method (details not shown). We have calculated an average range of ±30.9% between duplicates in our dataset. Duplicate blanks (empty filters stained with alcian blue) were also taken at every station. The Alcian Blue dye solution was calibrated just before the cruise using a standard solution of xanthan gum processed with a tissue grinder and subsequently filtered through two sets of filters (five points in triplicate).
Chla concentration was determined by filtering 150 mL of seawater on GF/F filters (Whatman), extracting the pigment in acetone (90% v:v) in the dark at 4°C for 24 h, and measuring fluorescence with a Turner Designs fluorometer.
Analyses of dissolved inorganic nutrient concentrations [nitrate (NO3), nitrite (NO2), phosphate (PO4), and silicate (SiO2)], were done by standard segmented flow analyses with colorimetric detection (Hansen and Grasshoff, 1983) using an Seal Analytical AA3 High Resolution AutoAnalyzer.
Particulate organic carbon (POC) was measured by filtering 1000 mL of seawater on pre-combusted GF/F glass fiber filters (4 h, 450°C). The filters were frozen in liquid nitrogen and kept at −80°C until analysis. Prior to analysis, the filters were dried at 60°C for 24 h. Then the filters were dried again and analyzed with a C:H:N autoanalyser (Perkin-Elmer 240).
For bacterial abundance samples, 1.8 ml were preserved with 1% paraformaldehyde + 0.05% glutaraldehyde (final conc.) and frozen in liquid nitrogen until processed in the lab. Bacterial abundance (BA) was analyzed by flow cytometry (FACSCalibur cytometer, Becton and Dickinson) after staining with SYBRGreen I (Molecular probes). Bacteria were detected by their signature in a plot of side scatter vs. FL1 (green fluorescence) as explained in Gasol and del Giorgio (2000).
Bacterial Production (BP) was estimated using the 3H-leucine incorporation method described by Kirchman et al. (1985). Three 1.2-mL aliquots and two trichloroacetic acid (TCA)-killed controls (5% final concentration) of each sample were incubated with 40 nmol L−1, 3H-leucine (epipelagic samples) or 80 nmol L−1, 3H-leucine (meso- and bathy-pelagic samples). The incubations were carried out in a water bath at in situ temperature in the dark. The incorporation was stopped by adding cold TCA (5% final concentration) to the vials, and samples were kept at −20°C until processing as described by Smith and Azam (1992). Radioactivity was then counted on a Beckman scintillation counter. Leucine incorporation rates were converted into carbon production using the conversion factor of 1.55 kg C produced per mole of leucine incorporated and considering no isotope dilution (Simon and Azam, 1989).
Given that TEP are frequently enriched in fucose (Zhou et al., 1998), we determined fucosidase activity using a fluorogenic substrate, as in Sala et al. (2016). Each sample (350 μl) was pipetted in quadriplicates into 96 black well-plates, with 50 μl of the substrate 4-methylumbelliferyl β-D-fucoside (Sigma-Aldrich) at a final concentration of 125 μM. Fluorescence was measured immediately after addition of the substrate and after incubations, at in situ temperature and in the dark, for 15, 30 min, 1, 3, and 5 h. The measurements were done with a Modulus Microplate (DISMED, Turner BioSystems) at 450 nm excitation and 365 nm emission wavelengths. The increase of fluorescence units during the period of incubation was converted into enzymatic activity with a standard curve prepared with 4-methylumbelliferone (MUF, Sigma-Aldrich).
We used the Statistica 7.0 software package to test the potential drivers of TEP distributions across the different spatial and timescales. We performed pairwise Pearson correlations between TEP concentrations and the following physico-chemical and biological variables: Temperature, salinity, turbidity, O2, nutrients (NO3, PO4, SiO4), particulate organic carbon (POC) and nitrogen (PON), chlorophyll a (Chla), bacterial abundance (BA), bacterial production (BP), and extracellular fucosidase activity. Data were log10-transformed and Bonferroni-corrected when needed.
Horizontal TEP Distribution along the Catalan Coast
In the coastal transects, TEP concentrations ranged from 27.4 to 122.8 μg XG eq L−1, and were overall higher in June (average 83.7 ± 23.9 μg XG eq L−1) than in May (average 56.1 ± 25.8 μg XG eq L−1; Table 1). In contrast, Chla concentrations were overall higher in May (0.65 ± 0.55 μg L−1) than in June (0.24 ± 0.14 μg L−1) and BP rates were similar in the two transects (0.18 ± 0.16 μg C L−1 h−1, ranging from 0.91 to 10.97, in May, 0.14 ± 0.21 μg C L−1 h−1, ranging from 0.50 to 14.79 in June). Dissolved inorganic nitrogen (DIN, nitrate+nitrite+ammonia) averaged 0.65 μM in May, ranging from 0.19 to 1.28 μM, and averaged 0.41 μM in June, ranging from 0.21 to 1.02 μM. Dissolved phosphate concentrations averaged 0.07 μM in May, ranging from 0.05 to 0.10 μM, and averaged 0.06 μM in June, ranging from 0.05 to 0.10 μM (Supplementary Figure 2). In May, TEP showed maxima in waters near Barcelona and north of the outflow of the Tordera River (Figure 2A). In these locations DIN concentration was 1.2 μM and phosphate concentration was 0.093 μM, two to eight-fold higher than in the rest of the stations. TEP-values were significantly correlated to Chla concentration (r = 0.93, p = 0.0003, n = 7) and marginally correlated to BP (r = 0.72, p = 0.06, n = 7, Figure 3). In June, TEP distributions showed maxima south of the outflow of the Tordera River, and were uncorrelated to Chla (Figure 2B) nor to BP. The TEP/Chla ratios were markedly higher in June (434.0 ± 197.5) than in May (136.7 ± 91.0, Table 1).
Table 1. Ranges of TEP concentration and TEP/Chla ratios in the different transects, depth profiles and diel cycles presented here.
Figure 2. Variations of TEP (blue symbols) and chlorophyll a (green symbols) concentration in the coastal transects performed in May (A) and June (B).
Figure 3. Scatterplot between Chla and TEP concentrations in the coastal transects performed in May (filled symbols) and June (open symbols). TEP were related to Chla in May (r = 0.93, p < 0.001, n = 7).
Horizontal TEP Distribution from Coastal to Open Sea Waters
The concentration of TEP in the studied coast-to-open sea transect ranged from 4.9 to 54.2 μg XG eq L−1 with a mean concentration of 18.7 ± 11.4 μg XG eq L−1.
We calculated depth-averaged TEP concentrations in epipelagic waters (0–200 m) in an attempt to look at horizontal distribution patterns. Depth-averaged epipelagic TEP ranged from 9.9 to 24.9 μg XG eq L−1 (Figure 4A). TEP concentrations were highest near the coast and 60 km offshore, at the slope to basin transition (station 6, Figure 4A). The same horizontal patterns were observed for Chla, O2, BP (Figure 4B), rendering significant correlations with TEP (r = 0.7, p < 0.05, n = 9). TEP were likewise related to POC (r = 0.9, p < 0.01, n = 9) and to the ratio between BP and O2(r = 0.7, p < 0.05, n = 9), which can be considered a proxy of bacterial reprocessing of photosynthetically fixed carbon.
Figure 4. Depth-averaged TEP concentration (μg XG eq L−1) and Chla (μg L−1) (A) and BP (μg C L−1 h−1) and O2 (mg L−1) (B) in epipelagic waters of the coast to open sea transect.
TEP Vertical Distribution in Epipelagic Waters
TEP vertical distribution patterns were variable among stations. Shelf waters exhibited mixed temperature and salinity profiles (Figure 5) and other chemical and biological variables were also quite uniform in the vertical profile. In these stations, Chla ranged from 0.10 to 0.38 μg L−1, BP ranged from 0.012 to 0.078 μg C L−1 h−1, and POC ranged from 3.7 to 7.2 μM. TEP vertical distributions were also quite homogenous from surface to the bottom, averaging 24.9 ± 3.0 and 19.7 ± 2.5 μg XG eq L−1 in stations 1 and 2, respectively (Figure 5). Conversely, well-developed DCM were detected in slope and basin waters between 50 and 60 m, with Chla concentrations ranging from 0.41 (station 5, slope) to 1.73 (station 6, basin) μg L−1. The TEP/Chla ratios ranged from 16.8 (station 6, DCM) to 316.0 (station 3, surface). They were generally higher at the surface and lower at the DCM. Bacterial production ranged from 0.005 to 0.065 μg C L−1 h−1 and the vertical distribution varied between stations: at the slope stations, BP was highest at the surface and subsurface, while in basin stations BP showed a bimodal profile, with peaks at the surface and at the DCM. POC concentrations, that ranged from 2.9 to 11.5 μM, showed similar distributions than TEP in the slope stations (with surface or subsurface peaks) but covaried with Chla, with maxima at the DCM, in the basin stations. TEP also showed marked vertical changes in slope and basin waters (Figure 5): they generally peaked at the surface (slope stations 3 and 5) or at the subsurface (slope station 4 and basin stations). These subsurface maxima were located between 25 and 55 m and showed values from 33.0 to 54.2 μg XG eq L−1 (mean concentration 39.6 ± 10.3 μg XG eq L−1). TEP maxima were always located shallower than the DCM and coincided with O2 maxima, and nutrient minima (Figure 5). At the slope, TEP maxima were also coincident with BP and POC maxima, while the Chla maximum was always deeper. By contrast, in the basin, BP and POC coincided with Chla while TEP and O2 peaked at shallower depths.
Figure 5. Vertical profiles of Temperature (°C), salinity (practical salinity units), turbidity (FTU) and NO3 (μM) (left panels) and O2 (mg L−1), TEP concentration (μg XG eq L−1), Chla (μg L−1) and BP (μg C L−1 h−1) (right panels) in epipelagic waters of station 2 (shelf, upper panels), station 4 (slope, middle panels), and station 8 (basin, bottom panels).
Drivers of TEP Vertical Distribution
We assessed which were the environmental drivers of TEP vertical distributions in epipelagic waters performing Pearson correlation tests between TEP concentrations and a number of physical (temperature, salinity, turbulence), chemical (nutrients, O2), and biological (Chla, BP, bacterial abundance, bacterial fucosidase activity) variables (Table 2). Since vertical distribution patterns of TEP and biological variables such as Chla and BP differed among shelf, slope and basin stations, we separated shelf stations from the others for the analysis. TEP was never significantly correlated to Chla. In shelf waters, TEP was not related to any other physicochemical or biological parameter. In slope and basin waters, TEP was significantly and positively correlated to O2, BP and the ratio BP/O2, and significantly negatively correlated to nutrients and N/P ratios (Table 2).
Table 2. Results of Pearson correlations between TEP concentration and different physical, chemical, and biological variables measured in epipelagic waters during the NEMO cruises at the slope and basin stations.
TEP Vertical Distribution in Mesopelagic and Bathypelagic Waters
TEP concentrations were half lower in meso- and bathypelagic waters with respect to the epipelagic, showing minima between 4.9 and 11.2 μg XG eq L−1. By contrast, relative TEP increases were observed in waters near the bottom in all slope and basin stations, parallel to increases in turbidity (proxy of total particle concentration).
TEP in meso- and bathypelagic waters were only significantly related to turbidity (r = 0.54, p < 0.01 n = 23), indicative of the relevance of bottom nepheloid layers (BNL). Remarkably, TEP were uncorrelated with POC in these layers. No significant correlation was found between TEP and BP, but a correlation with bacterial fucosidase activity (r = 0.45, p < 0.05 n = 23) could be observed.
Diel TEP Variations
The same surface water mass was sampled over time during the 2-day lagrangian study in the SUMMER2 cruise, as confirmed by plotting depth-averaged temperature and salinity values of each CTD cast (Supplementary Figure 1), except for the last two casts were a warmer water mass was likely sampled. We did not detect a recurrent diel pattern of TEP or of any other biological variables such as Chla or BP (Figure 6), even though these variables varied highly during the cycle. Chla showed subsurface maxima of 0.75–1.33 μg L−1 between 40 and 54 m. TEP concentrations ranged ten-fold, from 5.7 to 55.9 μg XG eq L−1 (average 34.1 ± 5.7 μg XG eq L−1). TEP maxima were situated between 25 and 47 m with highest values during the night (10 p.m. and 2 a.m.). These maxima were always shallower than the DCM, and were again vertically coincident with O2 and BP maxima. TEP concentrations were not significantly correlated to Chla, but a significant correlation was observed with O2 (r = 0.70, p = 0.000, n = 78), BP (r = 0.48 p = 0.000, n = 78), and the ratio BP/O2 (r = 0.44, p = 0.0001, n = 78). We looked at temporal variations of TEP, Chla, O2, and BP using depth-averaged values for all the epipelagic depths. Diel dynamics of TEP, Chla, O2, and BP were not coupled, but we detected Chla and O2 increases that were followed by a TEP increase after 4–12 h (Figure 6). However, lagged correlations between these variables, using either maxima or integrated values, were not significant.
Figure 6. Diel variations of depth-averaged TEP (blue circles), Chlorophyll a (green triangles), O2 (purple triangles) concentration.
TEP dynamics can affect particle aggregation rates because of stickiness, and particle sinking rates due to their low density. Thus, determining and predicting TEP dynamics is crucial if we want to accurately estimate carbon and particle fluxes. The TEP-values observed in this study are within the range of data published for Mediterranean Sea waters (Prieto et al., 2006; Ortega-Retuerta et al., 2010; Bar-Zeev et al., 2011) and at the lower range of other coastal ocean areas (Klein et al., 2011; Van Oostende et al., 2012). By contrast, the TEP/Chla ratios were in the upper range of those previously published for other ocean basins (Prieto et al., 2006), except for those measured in the Mediterranean Sea (Ortega-Retuerta et al., 2010; Bar-Zeev et al., 2011).
TEP Variation in the Coastal Transects
TEP maxima were observed near the city of Barcelona and near the outflow of the Tordera River, areas with higher nutrient concentrations. Coincidentally, there are two desalination plants located at the mouth of the Llobregat (next to Barcelona) and Tordera Rivers, whose feedwater intake pipes are located between 800 and 2200 m from the coast and at 30 m depth, although they are not always operative. Given that TEP are directly linked to membrane fouling, it is important to be able to predict TEP occurrence in these locations. In May, TEP variations in these coastal waters could be explained by variations in Chla, suggesting a direct linkage between phytoplankton and TEP at this geographical scale. However, this relationship was absent in June (Figure 2). The higher TEP/Chla ratios in June were due to both higher TEP concentrations than in May and lower Chla concentrations than in May. Indeed, primary productivity (PP) is on average two-fold higher in May than in June in coastal NW Mediterranean Sea waters (Gasol et al., 2016). Bacteria did not seem to be a significant TEP source in June as no significant correlations between TEP and BP or TEP and the BP/Chla ratios were observed. Hence, TEP may accumulate in the sea surface at the beginning of the stratification period, similar to other phytoplankton-derived organic matter (Avril, 2002; Vila-Reixach et al., 2012; Romera-Castillo et al., 2013). These substances may be not taken up immediately by bacterioplankton due to nutrient defficiency, which is common in summer (Sala et al., 2002; Pinhassi et al., 2006). What we could resolve in this analysis was that the distribution of TEP in coastal waters could be predicted from Chla, a variable that is frequently monitored, but only at certain periods of the year.
TEP Horizontal Distribution from Coastal to Offshore Waters
Our dataset also allowed concluding that Chla is a good predictor of TEP concentrations at the horizontal scale from the coast to the open sea, since significant positive relationships were observed between these two variables. This information is important because it would allow estimating TEP concentrations, for instance, using remote sensing Chla values. It is worthy to mention that, if we restrict the analysis to surface values only instead of depth-averaged values, TEP and Chla were also related (r = 0.63 p = 0.06, n = 9), while not remarkable covariations were found between TEP and other biological variables. This reinforces the possibility of using remote sensing Chla data to estimate the geographical distributions of TEP in the area.
At the horizontal scale, TEP was also positively related to O2, POC, and BP. The covariance between all these variables at the horizontal scale suggests physical forces that drive them all. We observed maxima of TEP in epipelagic waters near the coast and at the interface between slope and basin waters. In this area, which is located next to the Catalan front, a salinity doming is frequently observed (Font et al., 1988; Estrada and Salat, 1989), and relatively high Chla and PP values are found due to increased deep ocean nutrient availability at shallower depths (Estrada, 1996; Pedrós-Alió et al., 1999). This is, however, not evident in our dataset since some basin stations (i.e., stations 8 and 9) were sampled during stormy conditions that rendered more mixed temperature and salinity profiles, while the rest were sampled in sunny days and well-stratified profiles.
TEP Vertical Distribution in Epipelagic Waters
In contrast to horizontal TEP distributions, vertical TEP distributions could not be predicted from Chla at the vertical scale. In this case, TEP maxima were always located shallower than the DCM. TEP maxima above the DCM and near the surface have been reported in previous studies in the Alboran Sea (García et al., 2002), the Eastern Mediterranean Sea (Bar-Zeev et al., 2011), and along a West to East transect in the Mediterranean Sea (Ortega-Retuerta et al., 2010). Furthermore, in a recent study (Kodama et al., 2014), TEP maxima were associated to layers of maximum O2 and nutrient minima, similar to what we observed in our study. Correlation analyses confirmed the Chla-TEP decoupling. From the whole set of variables, those that explained TEP distribution in both slope and basin waters were O2 (positively related to TEP) and nitrate (negatively related to TEP, Table 2). This suggests that TEP are a direct product of PP, which is frequently decoupled from Chla concentration and whose highest rates are located at depths shallower than the DCM (Estrada et al., 1993). This suggests that Chla is not the best proxy of phytoplankton biomass or PP at the vertical scale. Indeed, increases in Chla concentration at the DCM reflect photoacclimation to low light levels through increases in the Chla/carbon ratio, and do not necessarily match highest phytoplankton biomass (Delgado et al., 1992; Gernez et al., 2011).
Given that we lacked direct PP estimates, in this study we consider O2 as a proxy of PP in this area. We are aware, however, that O2 concentration in the ocean is a result of biological processes as well as physical processes, such as ventilation in the upper mixed layer. But a non-significant correlation between O2 and temperature indicates that O2 distribution in the area majorly reflects biological processes. Furthermore, in surface waters of our study, TEP and O2 were not significantly correlated, which supports our view that the O2-TEP correlation is through PP with O2 being a proxy of the latter. Similarly, the absence of negative correlations between O2 and bacterial abundance or production suggests that respiration by heterotrophs is not the main driver of O2 distributions.
The observed significant positive correlation with bacterial production would suggest one of the following mechanisms: (1) bacteria act as a source of TEP, (2) bacterial colonization and utilization of TEP, or (3) dependence of both TEP and bacteria on other factors, namely organic compounds released by phytoplankton during primary production. We performed multiple regression analyses with TEP as the dependent variable and O2 and BP as independent variables. Both variables significantly explained TEP variability (r2 = 0.71, p = 0.000, n = 25), with partial coefficients of 0.72 for O2 and 0.46 for BP, respectively. Additionally, we correlated the raw residuals of the O2-TEP regression against BP, resulting significant and positive (r = 0.38, p < 0.05, n = 25). Both analyses lead us to conclude that both PP (with O2 as surrogate) and bacterial activity (with BP as surrogate) has a significant influence on vertical TEP distributions. Finally, the correlation between TEP and the BP/O2 ratio was significant and positive. The positive correlation indicates that the higher the bacterial processing of organic matter originated in primary production, the higher TEP concentrations we observed, suggesting then a synergy between carbon fixation and bacterial reutilization of this fixed C to generate TEP.
Nutrients may also have an impact on TEP distributions as suggested by the negative correlations with NO3, PO4, and SiO4. TEP are enriched in carbon respect to N and P (Mari et al., 2001), and different experimental studies have demonstrated a higher TEP release rate when nutrients are limiting (Mari et al., 2005; Pedrotti et al., 2010). Also, the N/P ratios are likely important in determining TEP production and degradation, as suggested by a significant correlation between the N/P ratio and TEP concentration in epipelagic waters. Although how nutrient stoichiometry affects TEP dynamics is unclear and probably dependent on the composition of the initial microbial community (Gärdes et al., 2012); our negative correlation suggests that the lower the relative proportion of P, the higher the TEP; which contrasts with previous experimental results (Engel et al., 2015). However, the likely limitation by phosphorus in our system (e.g., Sala et al., 2002; Pinhassi et al., 2006) in contrast to the Atlantic area studied by Engel et al. (2015) could differently affect TEP production. Indeed, it has been proposed that extracellular organic carbon production is highest under P limitation (Mauriac et al., 2011). TEP properties also vary depending on whether they are released during active growth or during bloom senescence (Mari et al., 2001), with implications for the fate of these particles in seawater (degradation vs. export). In the NW Mediterranean, PP usually peaks at the end of winter and spring (Gasol et al., 2016), so we sampled during the transition of spring to summer, which coincides with the beginning of nutrient depletion and associated decreases in PP. Thus, we expected TEP to accumulate and be prone to further export.
TEP Vertical Distribution in Meso- and Bathy-Pelagic Waters
In meso- and bathypelagic waters TEP distributions were only explained by turbidity changes. Specifically, we could detect TEP increases associated with near-bottom particle-rich layers (BNL). Surprisingly, these TEP increases were not paralleled by increases in POC. This suggests that these BNL were composed mostly of mineral particles that could be coated or aggregated by TEP. A previous study in the NW Mediterranean (Puig et al., 2013) observed the presence of fine particles in the BNL. Via microscope visualizations, they showed that organic matter in the BNL was mainly composed by “amorphous aggregates,” and suggested that these aggregates had lower sinking rates than phytoplankton cells or other solid organic particles. This is in line with our findings, where TEP, which are low density particles, may have a longer residence time in the BNL than the rest of POC compounds. In bathypelagic waters, TEP cannot have a direct phytoplanktonic source as light is absent. However, a bacterial source was not evident either since no significant correlations were observed between TEP and BA or BP. Interestingly, TEP concentrations in that layer were positively correlated to bacterial fucosidase activity. Since TEP are enriched in fucose (Zhou et al., 1998), this may reflect bacterial degradation of TEP in deep waters. Therefore, bacteria would act as a sink instead of a source for TEP, and a probable non-local TEP source must exist, material either sunk from epipelagic waters, resuspended from the sediment, or advected off the shelf.
TEP Diel Variations
To our knowledge this is the first time that high frequency (every 4 h) and short-term (i.e., 2-day) TEP changes have been monitored in the field during a lagrangian study. However, we did not find a recurrent pattern of any of the variables measured. In our case, this lack of diel patterns of microbial biomass/activity likely explained the absence of recurrent TEP diurnal or nocturnal maxima. Additionally, although we confirmed in this study the vertical decoupling between Chla and TEP and the better coincidence with O2 concentrations, the results of short-term variations of these variables were less conclusive. We got some hints about short-term temporal decoupling, where TEP peaks followed Chla and O2 peaks, but further work, with longer diel sampling, is needed to explore this issue.
We showed that the TEP-Chla relationship in the ocean is variable and mainly depends on the time and spatial scale studied. TEP can be predicted from Chla distributions at the horizontal scale, which opens the possibility to estimate surface TEP using remote sensing Chla; but this relationship is not evident at the vertical scale, nor at a short timescale, and also likely varies seasonally. Since our dataset is limited, in an attempt to compare our results to other areas, we compiled information on the various TEP-Chla relationships published in the literature (Table 3). The previous results mainly concur with our observations: TEP patterns mimic Chla patterns horizontally, but they are vertically decoupled, and TEP concentration maxima are frequently found at depths shallower than the DCM. Hence, we propose O2 concentrations and bacterial production as predictive variables for vertical TEP distributions in the ocean. Further, we suggest looking at the spatial and temporal variations of TEP together with primary productivity measurements and orienting further work to elucidate what is the specific role of bacterioplankton at explaining geographical and vertical TEP distributions in the ocean.
EO, MS, FP, RS, CM, and JG designed the work. EO, CM, FA, CA, RG, EB, and MM sampled and performed laboratory analyses and processed the data. EO wrote the manuscript with the help and inputs of all co-authors.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was funded by projects funded by the Spanish Ministry of Science STORM (CTM2009-09352/MAR), SUMMER (CTM2008-03309/MAR), DOREMI (CTM2012-34294), REMEI (CTM2015-70340-R), ANIMA (CTM2015-65720-R), PEGASO (CTM2012-37615), and Grup consolidat de Recerca de la Generalitat de Catalunya (2014SGR/1179). We thank the captain and crew of the R/V “García del Cid” and the Marine Technology Unit for their assistance in the field. We also thank the reviewers for helping improve the manuscript. EO was supported by Beatriu de Pinós (Generalitat de Catalunya) and Juan de la Cierva (Spanish Ministry of Economy) postdoctoral fellowships.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2016.02159/full#supplementary-material
Supplementary Figure 1. Bivariate plot between depth-averaged (upper mixed layer) temperature (abscises) and salinity (ordinates) in the 13 sampling stations along the diel cycle.
Supplementary Figure 2. Variations of dissolved inorganic nitrogen (DIN, orange circles), dissolved inorganic phosphorus (DIP, green-yellow triangles), and estimated bacterial production with factor 1.55 kgC mol leucine (BP, purple triangles) in the coastal transects in May (A) and June (B).
Alldredge, A. L., Passow, U., and Logan, B. E. (1993). The abundance and significance of a class of large, transparent organic particles in the ocean. Deep Sea Res. I 40, 1131–1140. doi: 10.1016/0967-0637(93)90129-Q
Azetsu-Scott, K., and Passow, U. (2004). Ascending marine particles: significance of transparent exopolymer particles (TEP) in the upper ocean. Limnol. Oceanogr. 49, 741–748. doi: 10.4319/lo.2004.49.3.0741
Bar-Zeev, E., Berman, T., Rahav, E., Dishon, G., Herut, B., Kress, N., et al. (2011). Transparent exopolymer particle (TEP) dynamics in the eastern Mediterranean Sea. Mar. Ecol. Progr. Ser. 431, 107–118. doi: 10.3354/meps09110
Beauvais, S., Pedrotti, M. L., Villa, E., and Lemée, R. (2003). Transparent exopolymer particle (TEP) dynamics in relation to trophic and hydrological conditions in the NW Mediterranean Sea. Mar. Ecol. Progr. Ser. 262, 97–109. doi: 10.3354/meps262097
Berman, T. (2013). Transparent exopolymer particles as critical agents in aquatic biofilm formation: implications for desalination and water treatment. Desal. Water Treat. 51, 1014–1020. doi: 10.1080/19443994.2012.713585
Bhaskar, P. V., and Bhosle, N. B. (2006). Dynamics of transparent exopolymeric particles (TEP) and particle-associated carbohydrates in the Dona Paula bay, west coast of India. J. Earth Sys. Sci. 115, 403–413. doi: 10.1007/BF02702869
Delgado, M., Latasa, M., and Estrada, M. (1992). Variability in the size-fractionated distribution of the phytoplankton across the Catalan front of the north-west Mediterranean. J. Plankton Res. 14, 753–771. doi: 10.1093/plankt/14.5.753
Engel, A., Borchard, C., Loginova, A., Meyer, J., Hauss, H., and Kiko, R. (2015). Effects of varied nitrate and phosphate supply on polysaccharidic and proteinaceous gel particle production during tropical phytoplankton bloom experiments. Biogeosciences 12, 5647–5665. doi: 10.5194/bg-12-5647-2015
Estrada, M., Marrasé, C., Latasa, M., Berdalet, E., Delgado, M., and Riera, T. (1993). Variability of deep chlorophyll maximum characteristics in the Northwestern Mediterranean. Mar. Ecol. Prog Ser. 92, 289–300. doi: 10.3354/meps092289
García, C. M., Prieto, L., Vargas, M., Echevarría, F., García-Lafuente, J., Ruiz, J., et al. (2002). Hydrodynamics and the spatial distribution of plankton and TEP in the Gulf of Cádiz (SW Iberian Peninsula). J. Plankton Res. 24, 817–833. doi: 10.1093/plankt/24.8.817
Gärdes, A., Ramaye, Y., Grossart, H. P., Passow, U., and Ullrich, M. S. (2012). Effects of Marinobacter adhaerens HP15 on polymer exudation by Thalassiosira weissflogii at different N:P ratios. Mar. Ecol. Prog. Ser. 461, 1–14. doi: 10.3354/meps09894
Gasol, J., and del Giorgio, P. (2000). Using flow cytometry for counting natural planktonic bacteria and understanding the structure of planktonic bacterial communities. Sci. Mar. 64, 197–224. doi: 10.3989/scimar.2000.64n2197
Gasol, J. M., Cardelús, C., Morán, X. A. G., Balagué, V., Massana, R., Pedrós-Alió, C., et al. (2016). Seasonal patterns in phytoplankton photosynthetic parameters and primary production in a coastal NW Mediterranean site. Sci. Mar. 80, 63–77. doi: 10.3989/scimar.04480.06E
Gernez, P., Antoine, D., and Huot, Y. (2011). Diel cycles of the particulate beam attenuation coefficient under varying trophic conditions in the northwestern Mediterranean Sea: observations and modeling. Limnol. Oceanogr. 56, 17–36. doi: 10.4319/lo.2011.56.1.0017
Harlay, J., Borges, A. V., Van Der Zee, C., Delille, B., Godoi, R. H. M., Schiettecatte, L. S., et al. (2010). Biogeochemical study of a coccolithophore bloom in the northern Bay of Biscay (NE Atlantic Ocean) in June 2004. Progr. Oceanogr. 86, 317–336. doi: 10.1016/j.pocean.2010.04.029
Kirchman, D., Knees, E., and Hodson, R. (1985). Leucine incorporation and its potential as a measure of protein-synthesis by bacteria in natural aquatic systems. Appl. Environ. Microbiol. 49, 599–607.
Klein, C., Claquin, P., Pannard, A., Napoléon, C., Le Roy, B., and Véron, B. (2011). Dynamics of soluble extracellular polymeric substances and transparent exopolymer particle pools in coastal ecosystems. Mar. Ecol. Progr. Ser. 427, 13–27. doi: 10.3354/meps09049
Kodama, T., Kurogi, H., Okazaki, M., Jinbo, T., Chow, S., Tomoda, T., et al. (2014). Vertical distribution of transparent exopolymer particle (TEP) concentration in the oligotrophic western tropical North Pacific. Mar. Ecol. Progr. Ser. 513, 29–37. doi: 10.3354/meps10954
Mari, X., Beauvais, S., Lemée, R., and Pedrotti, M. L. (2001). Non-Redfield C: N ratio of transparent exopolymeric particles in the northwestern Mediterranean Sea. Limnol. Oceanogr. 46, 1831–1836. doi: 10.4319/lo.2001.46.7.1831
Mari, X., Rassoulzadegan, F., Brussaard, C. P. D., and Wassmann, P. (2005). Dynamics of transparent exopolymeric particles (TEP) production by Phaeocystis globosa under N- or P-limitation: a controlling factor of the retention/export balance. Harmful Algae 4, 895–914. doi: 10.1016/j.hal.2004.12.014
Mauriac, R., Moutin, T., and Baklouti, M. (2011). Accumulation of DOC in Low Phosphate Low Chlorophyll (LPLC) area: is it related to higher production under high N:P ratio? Biogeosciences 8, 933–950. doi: 10.5194/bg-8-933-2011
Orellana, M. V., Matrai, P. A., Leck, C., Rauschenberg, C. D., Lee, A. M., and Coz, E. (2011). Marine microgels as a source of cloud condensation nuclei in the high Arctic. Proc. Natl. Acad. Sci. 108, 13612–13617. doi: 10.1073/pnas.1102457108
Ortega-Retuerta, E., Duarte, C. M., and Reche, I. (2010). Significance of bacterial activity for the distribution and dynamics of transparent exopolymer particles in the Mediterranean Sea. Microb. Ecol. 59, 808–818. doi: 10.1007/s00248-010-9640-7
Ortega-Retuerta, E., Passow, U., Duarte, C. M., and Reche, I. (2009a). Effects of ultraviolet B radiation on (not so) transparent exopolymer particles. Biogeosciences 6, 3071–3080. doi: 10.5194/bg-6-3071-2009
Ortega-Retuerta, E., Reche, I., Pulido-Villena, E., Agustí, S., and Duarte, C. M. (2009b). Uncoupled distributions of transparent exopolymer particles (TEP) and dissolved carbohydrates in the Southern Ocean. Mar. Chem. 115, 59–65. doi: 10.1016/j.marchem.2009.06.004
Passow, U., and Alldredge, A. L. (1995). A dye-binding assay for the spectrophotometric measurement of transparent exopolymer particles (TEP). Limnol. Oceanogr. 40, 1326–1335. doi: 10.4319/lo.1922.214.171.1246
Passow, U., Shipe, R. F., Murray, A., Pak, D. K., Brzezinski, M. A., and Alldredge, A. L. (2001). The origin of transparent exopolymer particles (TEP) and their role in the sedimentation of particulate matter. Cont. Shelf Res. 21, 327–346. doi: 10.1016/S0278-4343(00)00101-1
Pedrós-Alió, C., Calderon-Paz, J. I., Guixa-Boixereu, N., Estrada, M., and Gasol, J. M. (1999). Bacterioplankton and phytoplankton biomass and production during summer stratification in the northwestern Mediterranean Sea. Deep Sea Res. I 46, 985–1019. doi: 10.1016/S0967-0637(98)00106-X
Pedrotti, M. L., Peters, F., Beauvais, S., Vidal, M., Egge, J., Jacobsen, A., et al. (2010). Effects of nutrients and turbulence on the production of transparent exopolymer particles: a mesocosm study. Mar. Ecol. Progr. Ser. 419, 57–69. doi: 10.3354/meps08840
Pinhassi, J., Gómez-Consarnau, L., Alonso-Sáez, L., Sala, M., Vidal, M., Pedrós-Alió, C., et al. (2006). Seasonal changes in bacterioplankton nutrient limitation and their effects on bacterial community composition in the NW Mediterranean Sea. Aquat. Microb. Ecol. 44, 241–252. doi: 10.3354/ame044241
Prieto, L., Navarro, G., Cózar, A., Echevarría, F., and García, C. M. (2006). Distribution of TEP in the euphotic and upper mesopelagic zones of the southern Iberian coasts. Deep Sea Res. II 53, 1314–1328. doi: 10.1016/j.dsr2.2006.03.009
Prieto, L., Sommer, F., Stibor, H. N., and Koeve, W. (2001). Effects of planktonic copepods on transparent exopolymeric particles (TEP) abundance and size spectra. J. Plankton Res. 23, 515–525. doi: 10.1093/plankt/23.5.515
Puig, P., Madron, X. D. D., Salat, J., Schroeder, K., Martín, J., Karageorgis, A. P., et al. (2013). Thick bottom nepheloid layers in the western Mediterranean generated by deep dense shelf water cascading. Progr. Oceanogr. 111, 1–23. doi: 10.1016/j.pocean.2012.10.003
Radic, T., Kraus, R., Fuks, D., Radic, J., and Pecar, O. (2005). Transparent exopolymeric particles' distribution in the northern Adriatic and their relation to microphytoplankton biomass and composition. Sci. Total Environ. 353, 151–161. doi: 10.1016/j.scitotenv.2005.09.013
Ramaiah, N., and Furuya, K. (2002). Seasonal variations in phytoplankton composition and transparent exopolymer particles in a eutrophicated coastal environment. Aquat. Microb. Ecol. 30, 69–82. doi: 10.3354/ame030069
Rochelle-Newall, E. J., Mari, X., and Pringault, O. (2010). Sticking properties of transparent exopolymeric particles (TEP) during aging and biodegradation. J. Plankton Res. 32, 1433–1442. doi: 10.1093/plankt/fbq060
Romera-Castillo, C., Álvarez-Salgado, X. A., Galí, M., Gasol, J. M., and Marrasé, C. (2013). Combined effect of light exposure and microbial activity on distinct dissolved organic matter pools. A seasonal field study in an oligotrophic coastal system (Blanes Bay, NW Mediterranean). Mar. Chem. 148, 44–51. doi: 10.1016/j.marchem.2012.10.004
Sala, M. M., Aparicio, F., Balagué, V., Boras, J. A., Borrull, E., Cardelús, C., et al. (2016). Contrasting effects of ocean acidification on the microbial food web under different trophic conditions. ICES J. Mar. Sci. 73, 670–679. doi: 10.1093/icesjms/fsv130
Sala, M. M., Peters, F., Gasol, J. M., Pedrós-Alió C., Marrasé C., and Vaqué D. (2002). Seasonal and spatial variations in the nutrient limitation of bacterioplankton growth in the Northwestern Mediterranean. Aquat. Microb. Ecol. 27, 47–56. doi: 10.3354/ame027047
Scoullos, M., Plavšić M., Karavoltsos, S., and Sakellari, A. (2006). Partitioning and distribution of dissolved copper, cadmium and organic matter in Mediterranean marine coastal areas: the case of a mucilage event. Estuarine Coast. Shelf Sci. 67, 484–490. doi: 10.1016/j.ecss.2005.12.007
Sun, C. C., Wang, Y. S., Li, Q. P., Yue, W. Z., Wang, Y. T., Sun, F. L., et al. (2012). Distribution characteristics of transparent exopolymer particles in the Pearl River estuary, China. J. Geophys. Res. Biogeosci. 117. doi: 10.1029/2012JG001951
Taylor, J. D., Cottingham, S. D., Billinge, J., and Cunliffe, M. (2014). Seasonal microbial community dynamics correlate with phytoplankton-derived polysaccharides in surface coastal waters. ISME J. 8, 245–248. doi: 10.1038/ismej.2013.178
Van Oostende, N., Harlay, J., Vanelslander, B., Chou, L., Vyverman, W., and Sabbe, K. (2012). Phytoplankton community dynamics during late spring coccolithophore blooms at the continental margin of the Celtic Sea (North East Atlantic, 2006–2008). Progr. Oceanogr. 104, 1–16. doi: 10.1016/j.pocean.2012.04.016
Van Oostende, N., Moerdijk-Poortvliet, T. C., Boschker, H. T., Vyverman, W., and Sabbe, K. (2013). Release of dissolved carbohydrates by Emiliania huxleyi and formation of transparent exopolymer particles depend on algal life cycle and bacterial activity. Environ. Microbiol. 15, 1514–1531. doi: 10.1111/j.1462-2920.2012.02873.x
Vila-Reixach, G., Gasol, J. M., Cardelús, C., and Vidal, M. (2012). Seasonal dynamics and net production of dissolved organic carbon in an oligotrophic coastal environment. Mar. Ecol. Progr. Ser. 456, 7–19. doi: 10.3354/meps09677
Keywords: transparent exopolymer particles, chlorophyll a, bacteria, carbon, Mediterranean Sea
Citation: Ortega-Retuerta E, Sala MM, Borrull E, Mestre M, Aparicio FL, Gallisai R, Antequera C, Marrasé C, Peters F, Simó R and Gasol JM (2017) Horizontal and Vertical Distributions of Transparent Exopolymer Particles (TEP) in the NW Mediterranean Sea Are Linked to Chlorophyll a and O2 Variability. Front. Microbiol. 7:2159. doi: 10.3389/fmicb.2016.02159
Received: 12 September 2016; Accepted: 22 December 2016;
Published: 31 January 2017.
Edited by:Kai Ziervogel, University of New Hampshire, USA
Reviewed by:Lars-Eric Heimbürger, Mediterranean Institute of Oceanography (CNRS), France
Daniel Conrad Ogilvie Thornton, Texas A&M University, USA
Astrid Anne-Marie Gärdes, Leibniz Center for Tropical Marine Ecology, Germany
Copyright © 2017 Ortega-Retuerta, Sala, Borrull, Mestre, Aparicio, Gallisai, Antequera, Marrasé, Peters, Simó and Gasol. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Eva Ortega-Retuerta, firstname.lastname@example.org