The temporal and spatial changes of zooplankton and chlorophyll a concentration were studied during the warm stratification period (early June) at three stations whose traits corresponded to the coastal, frontal, and offshore-dome water conditions described for the Catalan Sea. We sampled the stations for 12 days at a frequency ranging from less than 10 to 102 h, with a spatial resolution ranging from 10 to 104 m. The objective was to determine the variability of mesozooplankton and phytoplankton (chlorophyll a) biomass, and average individual size (mass) across a coast-offshore transect in relation to the stratification conditions prevailing in the NW Mediterranean during summer. The vertical distribution of phytoplankton biomass displayed a clear deep maximum at 60 m depth except at the coastal station. This maximum exists during most of the year and is especially important during the density stratification period. It was accompanied during daylight hours by a coherent zooplankton maximum. At sunset mesozooplankton ascended and dispersed, with larger organisms from deeper layers joining the migrating community and increasing the average individual mass. The highest variability of mesozooplankton biomass, individual mass and chlorophyll a concentration occurred at the front station due to the coupling between the vertical migration of zooplankton and the particular characteristics of the front. According to the data shown, the highest variability was observed at the lowest scales.
RESUMEN
Durante el periodo de estratificación térmica (principios de junio) se estudiaron los cambios temporales y espaciales de concentración de zooplancton y clorofila a en tres estaciones cuyos rasgos corresponden a las condiciones de aguas costeras, frontal y domo central descritas para el Mar Catalán. Muestreamos las estaciones a lo largo de 12 días con una frecuencia que varió entre <10 h y 102 h, y con una resolución espacial entre 10 y 104 m. El objetivo fue determinar la variabilidad de la biomasa del fitoplancton (clorofila a) y del mesozooplancton, así como de su tamaño individual promedio a lo largo de un transecto costa-mar abierto en relación con las condiciones de estratificación térmica típicas del Mediterráneo del Noroccidental desde finales de primavera a finales de verano. La distribución vertical de la biomasa de fitoplancton mostró la existencia de un máximo profundo a 60 m de profundidad, excepto en la estación costera. Aunque cierta heterogeneidad vertical en la biomasa de fitoplancton suele observarse durante la mayor parte del año, es especialmente importante durante el período de estratificación de densidad. El máximo de fitoplancton estuvo acompañado durante las horas de luz por un máximo de zooplancton, que en la crepúsculo ascendía y se dispersaba, y al que se unían organismos de mayor tamaño procedentes de zonas más profundas, contribuyendo a aumentar la masa individual promedio. La mayor variabilidad de biomasa mesozooplanctónica, de la masa individual y de la concentración de clorofila a se observó en la estación frontal debido al acoplamiento entre la migración vertical del zooplancton y las características particulares del frente. Según los datos presentados la mayor variabilidad fue observada en las escalas espaciales y temporales más bajas.
KEYWORDSzooplankton biomassNW Mediterraneanspatial distributionshort-term variabilityPALABRAS CLAVEbiomasa de zooplanktonMediterráneo noroccidentaldistribución espacialvariabilidad de corto plazoINTRODUCTION
The study of the effects of physical singularities (i.e., mesoscale features and density stratification) on the characterisitics of plankton distribution, structure and function in the Catalan Sea started about four decades ago, in the framework of a series of research projects led by Dr Marta Estrada (see Marta Estrada bibliography in this volume). The projects merged different research groups of the former Instituto de Investigaciones Pesqueras (Institute of Fishery Research of Barcelona, now the Marine Sciences Institute) and the University of Barcelona, and received the invaluable advice and participation of Prof. Ramon Margalef. The area has thus been the object of successive studies almost unenterrumptedly, and this paper is a late consequence of one of them.
A better understanding of the role of zooplankton in food web dynamics and carbon balance in pelagic systems requires not only simultaneous estimates of the structural parameters and metabolic rates, but also to characterize their reliability in relation to the spatial and temporal scales considered. While rate processes seem to exhibit fast responses to physical or biological external forcing (Alcaraz 1988, Calbet et al. 1996), extensive properties such as biomass and taxonomic composition are considered to follow more gradual changes through cumulative processes (i.e. growth and reproduction, Alcaraz et al. 2007). However, their rate of change depends not only on demographic parameters. Behavioural characteristics of zooplankton and advection processes strongly determine their patchy nature and therefore the limits of what can be accepted as synoptic in terms of data representativeness and ecosystem dynamics (Steele and Henderson 1992).
In Stommel’s classic diagram about the time-space-intensity scales of variability in plankton, and of the factors contributing to them (Stommel 1963, Haury et al. 1978, Vance and Doel 2010), frequency and intensity of changes are inversely related. However, at the highest frequencies (small-scale variability) the magnitude of the changes is by no means negligible. Physical, chemical and social patterns (Stavn 1971), coupled with the effects of size and nature of patch formation and diffusion, among other parameters, impose severe limits on the acquisition of synoptic information of plankton structures (Longhurst et al. 1966, Haury 1976a,b). As a matter of fact, the problem is reduced to assessing how far in time and space a particular zooplankton sample can be considered representative (Winsor and Clarke 1940). In most cases, part of the variability derived from predictable sources such as tidal currents and light-induced activity rhythms can be avoided or at least significantly reduced by adopting suitable sampling strategies. However, as pointed out by Haury et al. (1978): “…if the trade-off is between economics, impatience and ecological convenience, usually impatience and economics win”.
In the NW Mediterranean, organismal behaviour and physical heterogeneity are important as variability sources, specifically for zooplankton (Pinot and Jansá 2001). A sub-surface maximum of phytoplankton productivity and biomass occurs most of the year (Estrada 1985) and is especially strong during the summer stratification period. This deep phytoplankton biomass and production layer is usually accompanied by a zooplankton biomass maximum (Alcaraz 1985, 1988, Alcaraz et al. 2010). Although the zooplankton biomass contributes a significant proportion of the short-term temporal and spatial variability due to the daily vertical migration of the organisms (Alcaraz 1985, Saiz and Alcaraz 1990, Alcaraz et al. 2010), highly dynamic hydrographic structures such as density fronts could be further sources of daily and seasonal zooplankton changes (Boucher et al. 1987, Fernández de Puelles et al. 2004, Mazzocchi et al. 2014).
The main objective of the Spanish research project VARIMED (AMB94-1019) was to improve our understanding of the matter and energy flux and carbon balance in the pelagic domain in the Catalan Sea (NW Mediterranean). To this end, we studied the short-term changes in the structural properties and rate processes for the different plankton components, from bacteria and phytoplankton (Acinas et al. 1997, Estrada 1996, Agawin and Agustí 1997) to zoooplankton and seabirds (Saiz et al. 1997, Abelló and Oro 1998). In this paper the variability in time and space scales ranging from hours to weeks and from metres to tens of kilometres of the structural properties of zooplankton is analysed in relation to phytoplankton (chlorophyll a, hereafter Chl a) and the physical characteristics of coastal-frontal-offshore conditions in a transect in the Catalan Sea. Accordingly, the two main objectives were: 1) to determine the variability pattern of biomass and size (individual mass) of zooplankton at the beginning of the summer stratification period, and 2) to estimate the possible interaction between the physical and biological (Chl a) and the space and time variability patterns observed.
METHODOLOGY
The study took place in the Catalan Sea (NW Mediterranean) during the oceanographic cruise VARIMED-95 (2-14 June 1995) on board the R/V Hespérides. The study was done in a 120-km transect approximately perpendicular to the coastline (Fig. 1), running from Barcelona to the middle of the Majorca–Menorca channel. This transect had been systematically sampled in previous oceanographic cruises (for a complete list, see Alcaraz et al. 2007).
Map of the study area with the position of the hydrographic stations in the 120 km transect between Barcelona and the Majorca-Menorca channel (dots). Large black dots: coastal (C), frontal (F) and offshore (O) stations. Shaded rectangle: area occupied by the grid of CTD casts made to characterize the area hydrographically.
Three structures of different characteristics define the hydrography and ecological peculiarities of the Catalan Sea (Font et al. 1988, Salat 1996, Alcaraz et al. 2007): the coastal zone, influenced by waters of continental origin, with an outer current running NE-SW parallel to the coast, the Liguro-Provençal-Catalan Current; the Catalan density front at the outer margin of the current, usually occluded during the summer stratification period but clearly marked below the thermocline; and a dome-like structure in the central part of the Catalan-Balearic Sea (Alcaraz 1988, Estrada and Margalef 1988). These structures are corresponding analogues of those described for the Nice-Calvi transect by Boucher (1984).
Sampling strategy
Previous to the biological study, we obtained quasi-synoptic information on the hydrography of the area (temperature, salinity and fluorescence profiles) with a Neil-Brown MARK-V CTD provided with a fluorometer, in a grid of stations (shaded rectangle in Fig. 1) including the study stations. According to the physical characteristics and Chl a distribution pattern, three biological stations were chosen as representative of the above-mentioned coastal, frontal and offshore conditions (hereafter C (41°20’N, 2°17’E), F (41°07’N, 2°27’E) and O (40°40’N, 2°53’E), located at 10, 35.9 and 89.4 km from the shore and above 80, 900 and 2000 m depth, respectively (big dots in Fig. 1).
The sampling schedule consisted of three visits (legs) to the three stations at regular intervals during a 12-day period (Table 1). Every 4 days, each station was sampled for 24 h at 4- to 8-h intervals, resulting from 3 to 5 zooplankton hauls per day, usually 1 to 2 during dusk or dawn hours, one at noon and 1to 2 during the night. For Station O we had only partial samples on the first day due to a failure of the net. In total, we made 29 hauls (10 at Station C, 11 at Station F, 8 at Station O), completing 198 depth-strata samples. Between successive legs and on the way back to the first sampled station (Station O), we spent 24 h doing fast CTD casts along the transect in order to detect changes in the physical structure that would make it necessary to re-consider the geographic position of the stations.
Sampling schedule during the VARIMED-95 cruise. The grey shadowed columns correspond to days dedicated to obtaining quick hydrographic surveys along the transect shown in Figure 1. The black zones correspond to biological sampling, generally repeated four times per day: 1) early morning, from 6 to 8.30 h; 2) midday, from 11.30 to 13.30 h; 3) evening, from 17 to 19 h; 4) night, around 22 h.
Julian day
153
154
155
156
157
158
159
160
161
162
163
164
Repet
1
2
3
4
1
2
3
4
1
2
3
4
1
2
3
4
1
2
3
4
1
2
3
4
1
2
3
4
1
2
3
4
1
2
3
4
C
F
O
Chlorophyll concentration
The vertical profiles of Chl a concentration were obtained by relating the fluorescence at selected depths obtained with the fluorometer installed in the CTD to the Chl a concentration analysed in water samples taken simultaneously at the depths of the corresponding fluorescence readings. Chl a was analysed on board fluorimetrically in acetone extracts according to Yentsch and Menzel (1963). The equation relating Chl a to fluorescence was µg Chl a L–1=0.132+14.26 Fl, r=0.88, n=77.
Mesozooplankton
After an initial CTD cast, according to the vertical structure of density and fluorescence, 6 to 8 depths intervals between 0 and either 60 m (St. C) or 100 m depth (Stations F and O) were chosen depending on the station depth. The mesozooplankton samples were obtained by descending oblique net tows made with a Longhurst-Hardy Plankton Recorder® net (hereafter LHPR) fitted with 200 µm-mesh netting and equipped with a flowmeter and a Seabird® CTD. These tows provided data on temperature, salinity and in vivo fluorescence, as well as average, minimum and maximum depth interval, and filtered water volume (m3) for each sampled depth interval. The water volume filtered for each depth interval ranged from 11.94 to 68.97 m3 (average=34.86 m3, standard deviation=15.67 m3). The zooplankton samples were preserved in formalin buffered with hexamethylene-tetramine to a final concentration of 4%.
Mesozooplankton biomass as organic carbon (hereafter CZOO) was estimated after biovolume determinations by applying the corresponding conversion factor between mesozooplankton biovolume (VZ) and organic C contents for mixed zooplankton communities of the NW Mediterranean (Alcaraz et al. 2003). Sample VZ was estimated by integrating the volumes of individual organisms corresponding to sample aliquots after excluding macroplankton (a few jellyfish, salps and siphonophores) by screening the sample through a 1-cm-mesh plastic grid. Images of the fixed samples were obtained with an EDI CCD video camera installed in a stereomicroscope connected to a computer. The images were analysed with a semi-automatic free image analysis software (NIH-Image®, National Institute of Health, Bethesda, Md., USA) that counted the individuals in the sample and measured the area contained by their silhouette, transforming it into an ellipse of the same area (Alcaraz et al. 2003). The volume of the organism was assimilated to that of the ellipsoid resulting from the rotation of the ellipse about its major axis. We used the ZV–CZOO transform factor given by Alcaraz et al. (2003) for formalin-preserved mixed zooplankton (gelatinous components excluded) of the western Mediterranean: 1 mm3 ZV=0.078 mg CZOO.
Aliquots of the samples corresponding to the different depth intervals of each haul were analysed, the biomass being calculated according to the fraction of the sample analysed and the water volume filtered by the net along the depth stratum:
CZOO-i = CZOO-a (F Vi)–1
(1)
where CZOO-i is the averge biomass (µgCZOO m–3) in the depth interval I, CZOO-a is the biomass of the aliquot analysed, F is the fraction of the total sample corresponding to the analysed aliquot, and Vi is the water volume filtered by the net along the depth interval. For graphic and statistical purposes, the sampling depth of each interval was considered to be the average interval depth. The abundance of mesozooplankton in individuals m–3 (IND) was calculated similarly, as
IND-i = IND-a (F Vi)–1
(2)
where IND-i is the number of individuals corresponding to the depth interval I, IND-a is the number of individuals in the aliquot analysed, and F and Vi are as in Equation 1.
The average zooplankton biomass for each haul was calculated as the integrated biomass for the whole haul divided by the haul depth,
(3)
where CZOO-h is the average haul mesozooplankton biomass, CZOO-i is the average biomass in the layer I, Zi+1 and Zi are the lower and upper depths of the considered water layers, Zt is the total depth of the haul and n is the total number of depth intervals sampled.
Average individual zooplankton biomass (IM) was calculated from the quotient between the average biomass and the average IND. CZOO, IM and IND are expressed in µgCZOO L–1, µgC ind–1 and individuals m–3, respectively. The variability was calculated as the coefficient of variation CV=100 (σ/μ). Statistical analyses were made using the JMP® software for Macintosh, the graphs about average depth distribution of mesozooplankton biomass (CZOO) using Systat® statistical and graphical software, and the mapping of the variables with the Surfer® software for PC.
RESULTSHydrography and chlorophyll a
The hydrographic structure along the transect corresponded to the typical thermal stratification situation of late spring-early summer in the NW Mediterranean (Estrada 1985, Estrada and Margalef 1988, Saiz and Alcaraz 1990), with the signature of the front occluded at surface by the mixing layer above the thermocline. However, the position of the front was identifiable at depth by the salinity and density gradients (Fig. 2). The deep chlorophyll layer (DCL) was well developed below the thermocline at around 60 m depth, except at the coast (station C), where in addition to the deep layer a surface maximum was visible (Fig. 2). The DCL was a persistent structure, the maximum concentration of Chl a increasing with the distance from the coast from 0.8 to 1.8 µg Chl a L–1 (Fig. 2). During the 12 study days, according to the three fast hydrographic surveys made along the transect (the initial one on Julian day [JD] 153 plus two between legs, JD 157 and 161, respectively), there were no significant changes in the physical structure of the study area (i.e. temperature or salinity profiles), in the position of the front, or in the depth of the DCL.
Distribution of temperature (A), salinity (B) and chlorophyll a (C) from 0 to 100 m depth along the transect between Barcelona and the Majorca-Menorca channel during the first study day (JD 153). In the horizontal axis, the distance to the coast refers to Barcelona.
CZOO, IM and IND
The vertical distribution of CZOO showed a maximum around 60 m depth (noon samples) in association with the DCL at the base of the thermocline. This pattern was clear at Stations F and O, whereas at Station C the zooplankton maximum occurred at 20 m depth, with a minimum at 40 m and a new secondary maximum at 60 m depth (Fig. 3). Overall, the depth-integrated average CZOO and IND for the three stations indicated a positive coast-offshore gradient, with the highest values offshore (Table 2). At time and space resolutions ranging from 4 to 6 h and from 10 to 20 m depth, respectively (Fig. 4), the relative maximum of CZOO was at around 60 m depth for noon samples, coinciding with the DCL (stations F and O), while at Station C it was less conspicuous. IM appeared to be inversely related to CZOO at Stations C and F. At Station O the loss of part of the data made it impossible to confirm this observation.
Box-plot of the average depth distribution of CZOO (µgC L–1). C, F and O as in Figure 1. At Station C the maximum depth stratum sampled was 60 m, as the station depth was only 80 m. *: Outlier.
Average values and variability (CV, in parentheses) of CZOO (µg C L–1), IM (µg C ind.–1) and IND (individuals m–3) at coastal (C), frontal (F) and offshore (O) stations. The time- and space scales are 102 h and 102 m (every 4 days, 100 m water column integrations). Station C, 10 hauls, 63 depth strata; Station F, 11 hauls, 75 depth strata; Station O, 8 hauls, 60 depth strata.
Variable
C
F
O
CZOO
1.34 (39)
1.81 (61)
2.19 (31)
IM
7.29 (43)
11.06 (62)
8.46 (32)
IND
370 (74)
533 (97)
696 (65)
Chl a and CZOO depth distribution at noon along a 24 h series for one of the legs (JD 162-164). The trends were similar for legs JD 154-156 and JD158-160. Horizontal axis, GMT+2; vertical axis, depth in m. C, F and O as in Figure 1.
In a space-time-CZOO diagram (time-space scales: 104 m and 102 h, respectively), the daily average CZOO followed the above-mentioned coast-offshore positive gradient, which was more visible in the final part of the cruise (Fig. 5). At the same time, for Stations C and F there was a decrease in CZOO over time, while at Station O the tendency seemed to be the opposite. No significant changes were observed between day and night samples. The IM (Fig. 6) showed a tendency to decrease during the study in day samples, while the opposite pattern was observed in night samples, especially at Station F.
Depth-integrated average values of mixed day and night samples of mesozooplankton biomass (CZOO, ) and its variability (CV) along the studied transect during the period analysed. Horizontal axis, sampling day (JD); vertical axis, distance from the coast (km). The integrated depth was 60 m for Station C and 100 m for Stations F and O.
Average values of biomass (CZOO) and individual size (IM) and their variability (CV) for day and night samples. Horizontal and vertical axes as in Figure 5.
Variability
At the time and space scales of 102 h and 102 m (successive legs for a given station and for the whole 0-100 m water column, respectively), the average values and variability of CZOO, IM and IND at Stations C, F and O are indicated in Table 2. The highest variability corresponded to Station F and the lowest to Station O).
Table 3 shows the variability of CZOO, IM and Chl a for all the time and space scales considered. The variability of CZOO was inversely related to the time and space scales, with CVs ranging from 59 (4-6 h, 10 m) to 18 (102 h, 104 m). IM showed a similar trend, but there was a significant increase in the CV to 81 at time scales coinciding with the day-night cycle (about 10 h, 102 m), and similar results were found for Chl a (CV = 47).
Variability (CV, %) of zooplankton biomass, CZOO (µg C L–1), individual mass, IM (µg C ind.–1), and Chl a, (µg chlorophyll a L–1) at the time- and space scales explored during the survey.
Scales
Parameters
Time
Space
CZOO
IM
Chl a
(h)
(m)
(µg C L–1)
(µg C ind–1)
(µg Chl a L–1)
4-6
101
59
70
36
4-6
102
41
67
30
101
102
38
81
47
102
102
24
58
12
102
104
18
35
16
In Figure 5, the time-space-variability diagram of zooplankton biomass indicates the higher variability associated with Station F (CV>60) and a tendency to increase towards the end of the cruise. The highest variability of IM was directly associated with the average individual size in daylight samples (Fig. 6), while during the night the tendency was the opposite.
The variability associated with the combined time and space scales included in the study for CZOO, IM and Chl a is represented in Figure 7. According to this diagram, the highest variability of CZOO was associated with the smallest time scales (<10 h, and ≈10 m; repeated samples and depth interval, CV≈60), decreasing monotonically according as the time and space scales increased. The maximum variability of IM took place at time scales of ≈10 h and 102 m, equivalent to day-night changes along the whole water column (vertical migration). For Chl a the highest variability (CV≈50) was associated with similar time and space scales to those of IM. The lowest variability for mesozooplankton and chlorophyll occurred at scales of 102 h (about 4 days, between successive legs) and 104 m, a scale determined by the distance between stations. The variability for IND has not been represented, as it was very similar to that of CZOO.
Contours of the variability (CV) of CZ, IM, and Chl a at the spatial and temporal scales considered in the study, from <10 (depths strata) to 104 m (between stations), and from <10 h (between successive daily hauls) to 102 h (between successive legs). Horizontal axis, log space scale (m); vertical axis, log time scale (h).
DISCUSSION
Most of the problems when studying variability in plankton systems arise from difficulties in obtaining synoptic or quasi-synoptic samples in physically unstable environments at the suitable spatial resolution. In the case of zooplankton, the problem is increased by the capacity of most of its components to conduct relatively important movements in the vertical plane which, combined with horizontal shear, can lead to major advection (Steele and Henderson 1992).
Although the resolution of most of the existing sampling techniques for zooplankton is not absolutely satisfactory, it is possible to obtain quantitative zooplankton samples separated into spatial (horizontal or vertical) strata fast enough to be compatible with the rate of change of the community (Donaghay et al. 1992). The displacement velocity of the gravity centre of zooplankton populations ranges from less than 1 m h–1 to 7 m h–1 (Krause and Radach 1989), although for large, fast swimming organisms it can be significantly higher (Enright 1977, Smith at al 1989).
Although Fasham et al (1974) and Fasham (1978) concluded that the spatial variability (patchiness) of plankton cannot be completely apprehended through LHPR hauls, according to Haury (1976a) the samples thus obtained should provide reasonable snapshots of the vertical structure of zooplankton biomass. In our case, the 0-100 m depth layer was sampled in less than 30 min, so the sampling of the different depth strata can be considered as quasi-synoptic.
Regarding the estimates of mesozooplankton biomass through the volume-carbon relationships and volume estimations by semi-automatic image analysis, the method has proven to be reliable. The overall variability (CV), including sub-sampling plus individual volume and C analytical errors, ranges from 9% to 20% (Alcaraz et al. 2003, Alcaraz et al. 2010). The obtained values of mesozooplankton carbon (CZOO) were significantly lower (from 50% to 75%) than previous values obtained by low-pressure reverse-flow water filtration of discrete water volumes in the same area and time of the year (Alcaraz 1980, Alcaraz et al. 1985, Calbet et al. 1996). This is due to the higher retention efficiency of zooplankton organisms by low-pressure filtration methods in comparison with hauling plankton nets with the same mesh-size (Alcaraz 1977). Calbet et al. (2001) obtained similar results when comparing samples taken with a 200-µm mesh Juday-Bogorov net with simultaneously obtained microplankton samples that were low-pressure filtered through a 200-µm gauze. Further discussions on the problem can be found in Turner (2004), and comments on the sampling bias of the LHPR net in Haury (1973). However, the quantitative underestimation of biomass by oblique plankton hauls should not affect the relative spatial pattern and temporal variability of the extensive properties of zooplankton.
According to the physical structure observed, the three stations corresponded to the typical coastal, frontal and offshore conditions described in Font et al. (1988), Saiz et al. (1992) and Calbet et al. (1996, 2002) for the same area and time of the year. While the coast-offshore gradient in zooplankton biomass was opposed to previous observations (Alcaraz et al. 1994, Calbet et al. 2002), the pattern of vertical distribution at noon repeated the typical deep zooplankton biomass maximum (DZM) associated with the phyto- and chlorophyll deep maximum (Alcaraz 1985, Saiz and Alcaraz 1990). No previous data on organism size allow us to generalize the average larger IM of zooplankton organisms associated with the density front, although it agrees with its nature as a nursery area for copepods (Boucher et al. 1987), where euphausiids accumulate and fish larvae are retained (Sabatés and Olivar 1996), thus enlarging the size-range of zooplankton organisms gathered in the area.
The higher variability of zooplankton biomass at small time- and space scales could probably result, as pointed out by Haury (1976a), from the combination of the scale size of zooplankton microstructure with the constant vertical water motion due to internal wave activity. For depth-integrated (102 m) and day-night (≈10 h) scales the variability is lower, or of the same order as the values found in single LHPR hauls by Haury (1976b) for organism numbers in the Pacific Central Gyre and the California current, and the average IM shows the highest variability at the same time and space scales. This finding could be explained by the high number of samplings coinciding with the time of maximum ascension speed at dusk of large, deep zooplankton components, and their descent at dawn. The higher backscatter strength recorded by moored ADCP in the same area and time of the year during the morning and evening migrations of pelagic organisms (Pinot and Jansá 2001), and off Southern California (Smith et al 1989), seem to confirm this hypothesis. In general, the higher variability of CZOO and IND as compared with Chl a is the consequence of the non-motile nature of most of the phytoplankton components, and of the variability imposed by vertical migration on zooplankton.
Although the spatial-temporal diagram of haul-integrated zooplankton biomass and its variability indicate a persistence of the coast-offshore positive gradient, it also indicates a considerable degree of variability, especially at the front. All the data identify it as a hotspot of zooplankton biomass of larger size, during both day and night. The front is thus an ecotone between hydrographically distinct systems: the Liguro-Provençal-Catalan Current on the coast side, the dome waters offshore and, occluding everything during the summer stratification period, an almost uniform mixing layer, each system having different zooplankton assemblages (Boucher 1984, Calbet et al. 1996).
Further studies in contrasting hydrographic situations (i.e. the winter mixing period or the spring outbreak and autumn collapse of the thermocline) would help to identify the relevant variability scales for extensive properties of the different plankton components and improve the parameterization in plankton production models.
ACKNOWLEDGEMENTS
We are deeply grateful to Marta Estrada for her friendship, stimulating inspiration and collaboration over almost five decades. This paper is intended as a small homage to her kind and inspired leadership.
The VARIMED-95 cruise took place in the framework of the Spanish research project AMB94-1019 from the Ministry of Science to M. Alcaraz. Financial support was also provided by the projects MAR98-0854 and CTM2006-12344-c02 to M.A., TOPCOP-CTM2011-23480 to E. Saiz and FERMI (CGL2014-59227-R) to AC and ES. The authors are indebted to all the participants in the VARIMED-95 cruise, and to the crew of the RV Hespérides. This is a contribution of the Marine Zooplankton Ecology Research Group of the Generalitat de Catalunya (2014 SGR 498).
REFERENCES Abelló P. Oro D. Offshore distribution and assemblages of breeding seabirds in the Catalano-Balearic Sea (northwestern Mediterranean). 199821422426 Acinas S.G. Rodríguez-Valera F. Pedrós-Alió C.Spatial and temporal variation in marine bacterioplankton diversity as shown by RFLP fingerprinting of PCR amplified 16S rDNA 1997242740http://dx.doi.org/10.1111/j.1574-6941.1997.tb00420.x Alcaraz M. Muestreo cuantitativo de zooplancton: análisis comparativo de la eficacia de mangas y botellas en un sistema estuárico 197741258294 Alcaraz M. Evolución y distribución vertical de la biomasa de zooplancton expresada en carbono y nitrógeno orgánicos, relación C/N y carbono detrítico, en una zona marina afectada por aguas residuales de la ciudad de Barcelona198044265274 Alcaraz M. Vertical distribution of zooplankton biomass during summer stratification in the Western Mediterranean Gibbs P.E.1985 Cambridge Cambridge Univ. Press 135143 Alcaraz M. Summer zooplankton metabolism and its relation to primary production in the Western Mediterranean Minas H.J. Nival P.1988 Oceanol. Acta Nº SP 9185191 Alcaraz M. Estrada M. Flos J. Particulate carbon and nitrogen and plankton biomass in oligotrophic and upwelling systems Bas C. Margalef R. Rubies P.1985BarcelonaCSIC435448 Alcaraz M. Saiz E. Estrada M. Excretion of ammonia by zooplankton and its potential contribution to nitrogen requirements for primary production in the Catalan Sea (NW Mediterranean). 19941196976http://dx.doi.org/10.1007/BF00350108 Alcaraz M. Saiz E. Calbet A. Estimating zooplankton biomass through image analysis.2003143307315http://dx.doi.org/10.1007/s00227-003-1094-8 Alcaraz M. Calbet A. Estrada M. Physical control of zooplankton communities in the Catalan Sea200774294312http://dx.doi.org/10.1016/j.pocean.2007.04.003 Alcaraz M. Almeda R. Calbet A. The role of arctic zooplankton in biogeochemical cycles: respiration and excretion of ammonia and phosphate during summer 20103317191731http://dx.doi.org/10.1007/s00300-010-0789-9 Agawin N.S.R. Agustí S. Abundance, frequency of dividing cells and growth rates of Synechococcus sp. (cyanobacteria) in the stratified Northwestern Mediterranean Sea 19971915991615http://dx.doi.org/10.1093/plankt/19.11.1599 Boucher J. Localization of zooplankton populations in the Ligurian marine front: role of ontogenic migration. 198431469484http://dx.doi.org/10.1016/0198-0149(84)90097-9 Boucher J. Ibanez F. Prieur L. Daily and seasonal variations in the spatial distribution of zooplankton populations in relation to the physical structure in the Ligurian sea front 198745133173http://dx.doi.org/10.1357/002224087788400891 Calbet A. Alcaraz M. Saiz E. Planktonic herbivorous food webs in the Catalan Sea (NW Mediterranean): temporal variability and comparison of indices of phyto-zooplankton coupling based in state variables and rate processes 19961823292347 Calbet A. Garrido S. Saiz E. Annual zooplankton succession in coastal NW Mediterranean waters: the importance of the smaller size fractions200123319331 Calbet A. Saiz E. Alcaraz M. Copepod eggproduction in the NW Mediterranean: effects of winter environmental conditions 2002273173184 Donaghay P.L. Rines H.M. Sieburt J.McN Simultaneous sampling of fine scale biological, chemical and physical structure in stratifed waters 19923697108 Enright J.T. Copepods in a hurry: sustained high-speed upward migration 197722118125 Estrada M. Deep phytoplankton and chlorophyll maxima in the Western Mediterranean Moraitou-Apostoulou M. Kiortsis V.1985 New York Plenum Press 247277 Estrada M. Primary production in the northwestern Mediterranean 1996605564 Estrada M. Margalef R. Océanographie pélagique méditerranéenne Minas H.J. Nival P. Supply of nutrients to the Mediterranean photic zone along a persistent front 19889133142 Fasham M.J.R. The statistical and mathematical analysis of plankton patchiness. 1978164379 Fasham M.J.R. Angel M.V. Roe S.J. An investigation of the spatial pattern of zooplankton using the Longhurst-Hardy plankton recorder 19741693112 Fernández de Puelles M.L. Valencia J. Jansá J. Hydrographical characteristics and zooplankton distribution in the Mallorca channel (Western Mediterranean): spring 2001200461654666 Font J. Salat J. Tintoré J. Permanent features of the circulation in the Catalan sea Minas H.J. Nival P. Pelagic Mediterranean oceanography 198895157 Haury L.R. Sampling biass of a Longhurst-Hardy Plankton Recorder 197318500506 Haury L.R. Small-scale pattern of a California current zooplankton assemblage 197637137155http://dx.doi.org/10.1007/BF00389124 Haury L.R. A comparison of zooplankton patterns in the California current and North Pacific Central Gyre 197637159167http://dx.doi.org/10.1007/BF00389125 Haury L.R. Mc Gowan J.A. Wiebe P.H. Patterns and processes in the time-space scales of plankton distributions Steele J.H.1978NY Plenum Press 277327http://dx.doi.org/10.1007/978-1-4899-2195-6_12 Krause M. Radach G. On the relations of vertical distribution, diurnal migration and nutritional state of herbivorous zooplankton in the northern North Sea during FLEX 1976.198974371417 Longhurst A.R. Reith A.D. Bower R.E.A new system for the collection of multiple serial plankton samples 196613213222http://dx.doi.org/10.1016/0011-7471(66)91101-6 Mazzocchi M.G. Siokou I. Tirelli V. Regional and seasonal characteristics of epipelagic mesozooplankton in the Mediterranean Sea based on an artificial neural network analysis 20141356480http://dx.doi.org/10.1016/j.jmarsys.2013.04.009 Pinot J.M. Jansá J. Time variability of acoustic backscatter from zooplankton in the Ibiza channel (western Mediterranean)20014816511670http://dx.doi.org/10.1016/S0967-0637(00)00095-9 Sabatés A. Olivar P. Variation of larval fish distributions associated with variability in the location of a shelf-slope front 19961351120http://dx.doi.org/10.3354/meps135011 Saiz E. Alcaraz M. Pigment gut contents of copepods and deep phytoplankton maximum in the Western Mediterranean 199012665672http://dx.doi.org/10.1093/plankt/12.3.665 Saiz E. Rodriguez V. Alcaraz M. Spatial distribution and feeding rates of Centropages typicus in relation to frontal structures in the Catalan Sea (Western Mediterranean). 19921124956http://dx.doi.org/10.1007/BF00349727 Saiz E. Calbet A. Trepat I. Food availability as a potential source of bias on the egg production method for copepods 199719114http://dx.doi.org/10.1093/plankt/19.1.1 Salat J. Review of hydrographic environmental factors that may influence anchovy habitats in northwestern Mediterraneaan 1996602132 Smith P.E. Ohman M.D. Eber L.E. Analysis of the patterns of distribution of zooplankton aggregations from an accoustic doppler current profiler 19893088103 Stavn R.H.
Stavn R.H. 1971. The horizontal-vertical distribution hypothesis: Langmuir circulations and Daphnia distributions 197116453446http://dx.doi.org/10.4319/lo.1971.16.2.0453 Steele J.H. Henderson E.W. A simple model for plankton patchiness 19921413971403http://dx.doi.org/10.1093/plankt/14.10.1397 Stommel H. Varieties of oceanographic experience 1963139572576http://dx.doi.org/10.1126/science.139.3555.572 Turner J.T. The importance of small planktonic copepods and their roles in pelagic marine food webs 200443255266 Vance T.C. Doel R.E. Graphical methods and cold war scientific practice: the Stommel diagram's intriguing journey from the physical to the biological sciences 201040147http://dx.doi.org/10.1525/hsns.2010.40.1.1 Winsor C.P. Clarke G.L. A statistical study of variation in the catch of plankton nets 19403134 Yentsch C.S. Menzel D.W. A method for the determination of phytoplankton chlorophyll and pheophytin by fluorescence 196310221231