# Consequence of a sudden wind event on the dynamics of a coastal phytoplankton community: an insight into specific population growth rates using a single cell high frequency approach

Mathilde Dugenne, Melilotus Thyssen, David Nerini, Claude Manté, Jean-Christophe Poggiale, Nicole Garcia, Fabrice Garcia, Gérald Grégori

► To cite this version:

Mathilde Dugenne, Melilotus Thyssen, David Nerini, Claude Manté, Jean-Christophe Poggiale, et al. Consequence of a sudden wind event on the dynamics of a coastal phytoplankton community: an insight into specific population growth rates using a single cell high frequency approach. Frontiers in Microbiology, 2014, 5, 10.3389/fmicb.2014.00485. [https://doi.org/10.3389/fmicb.2014.00485](https://doi.org/10.3389/fmicb.2014.00485). [hal-01146092](https:// HAL Id: )

Submitted on 12 Oct 2018

**HAL** is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire **HAL**, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Consequence of a sudden wind event on the dynamics of a coastal phytoplankton community: an insight into specific population growth rates using a single cell high frequency approach

Mathilde Dugenne*, Melilotus Thyssen, David Nerini, Claude Mante, Jean-Christophe Poggiale, Nicole Garcia, Fabrice Garcia and Gérald J. Grégori

Aix Marseille Université, Université de Toulon, CNRS/INSU, IRD, Mediterranean Institute of Oceanography, UM 110, Marseille, France

**Edited by:**

Susann Müller, Helmholtz Centre for
Environmental Research Leipzig
Germany, Germany

**Reviewed by:**

Peter Bergholz, North Dakota State
University, USA

Steven Singer, Lawrence Berkeley
National Laboratory, USA

Peter Allan Thompson,
Commonwealth Scientific Industrial
Research Organisation, Australia

**\*Correspondence:**

Mathilde Dugenne, Aix Marseille
Université, CNRS/INSU, IRD,
Mediterranean Institute of
Oceanography, UM 110,
13288 Marseille, France
e-mail: mathilde.dugenne@
univ-amu.fr

Phytoplankton is a key component in marine ecosystems. It is responsible for most of the marine primary production, particularly in eutrophic lagoons, where it frequently blooms. Because they are very sensitive to their environment, the dynamics of these microbial communities has to be observed over different time scales, however, assessment of short term variability is often out of reach of traditional monitoring methods. To overcome these limitations, we set up a Cytosense automated flow cytometer (Cytobuoy b.v.), designed for high frequency monitoring of phytoplankton composition, abundance, cell size, and pigment content, in one of the largest Mediterranean lagoons, the Berre lagoon (South-Eastern France). During October 2011, it recorded the cell optical properties of 12 groups of pico-, nano-, and microphytoplankton. Daily variations in the cluster optical properties were consistent with individual changes observed using microscopic imaging, during the cell cycle. We therefore used an adaptation of the size-structured matrix population model, developed by Sosik et al. (2003) to process the single cell analysis of the clusters and estimate the division rates of 2 dinoflagellate populations before, during, and after a strong wind event. The increase in the estimated *in situ* daily cluster growth rates suggest that physiological changes in the cells can prevail over the response of abundance.

**Keywords:** phytoplankton, flow cytometry, high frequency, *in situ* growth rates, Berre lagoon, wind

## INTRODUCTION

Phytoplankton is responsible for about 50% of the annual global net primary production, yet its biomass represents only about 2% of the entire chlorophylous biomass (Field et al., 1998). This compartment is made up of mostly drifting autotrophic microorganisms, whose dynamics are important for carbon and nutrient fluxes in the ecosystem. With a high turnover compared to terrestrial plants (Margalef, 1978), phytoplankton plays an essential role in driving the biogeochemical cycles and the Redfield ratio in the ocean (Goldman et al., 1979). However, obtaining an accurate description of the dynamic distribution of phytoplankton can be limited by the methods available for measuring, on a short temporal scale, the processes which may influence it. Various methods have now been developed to investigate phytoplankton, from the single cell level (by microscopy and flow cytometry after *in situ* sampling), the global scale using bulk methods (*in situ* fluorimetry and metagenomics) up to the ocean scale using remote sensing. The ocean color obtained by satellite observations indicates bulk chlorophyll concentrations, used as a proxy for phytoplanktonic biomass, can be highly variable in coastal waters (Antoine et al., 1996). This variability is due to forcings, that control net population growth, grazing, and physical transport (Cloern, 1996), acting over hourly or daily scales. Since

conventional sampling strategies hardly reach these scales, a legitimate concern has been raised over the underestimation of the natural short-term dynamics of phytoplankton during oceanographic campaigns and routine surveys (Rantajärvi et al., 1998). Or it is well known that the short generation time of phytoplankton can explain their ability to coexist on the same limited resources (Hutchinson, 1961).

The most efficient method for quantifying phytoplankton productivity is the estimation of the dividing rate of a population (in relation to the cell cycle) indicative of NPP. This combined with the variation in abundance is the only method for quantifying loss rates which is a good indicator of the sensitivity of the phytoplankton to trophic interactions or abiotic forcings. Estimations of dividing rates, expressed as population growth rate, are generally based on DNA measurements with lags from an hour to a day (Binder et al., 1996; Vaulot and Marie, 1999). Unfortunately, this requires important manual handling (i.e., sampling, fixing, storing, staining, and analysis) of the samples. To address these limitations, priority has been given to autonomous flow cytometry (Peeters et al., 1989; Olson et al., 2003), like the Cytosense flow cytometer (Cytobuoy b.v., Netherlands) in order to observe fine-scale temporal dynamics of phytoplankton (Dubelaar and Gerritzen, 2000; Dubelaar and Jonker, 2000; Dubelaar et al.,
2004; Thyssen et al., 2008, 2014). Single cell analysis using several optical properties gives a good alternative for estimating growth rates, as they are linked to cell size. As previously reported in Durand (1995) and Binder et al. (1996), periodic increase and decrease in light scatter intensity can be interpreted as a response of cellular growth and division occurring during the cell cycle. Diel variations in the phytoplankton cell dimensions based on flow cytometry proxy therefore represent a necessary compromise for automated measurement requirements in order to calculate growth rates of natural phytoplanktonic groups during high frequency surveys (Sosik et al., 2003). Moreover, image processing allow to further improve the estimations of cell dimensions using mounted imaging systems implemented in the last generation of instruments (Sieracki et al., 1998).

In this study, we report on a high frequency *in situ* monitoring of phytoplankton in the Berre lagoon during October 2011. The Berre lagoon, located in South Eastern France is one of the largest brackish lagoons off the Mediterranean shore. Since 1995, the site has benefited from a special attention relying on a local ecological survey conducted by the GIPREB Joint Union (Groupement d'Intêret Public pour la Réhabilitation de l'Etang de Berre) to observe and report on the eutrophication of the lagoon. However, in this ecosystem, the influence of important seasonal events like natural floods or strong wind events, which may have an influence on light and nutrient availability, are poorly documented and may limit the full understanding of the mechanisms driving the phytoplanktonic variability (Gouze et al., 2008). If previous studies have reported a high eutrophic state characterized by the dominance of only a few species, more recent data highlights a high biodiversity, superior to that found in the nearby Bay of Marseille (Thyssen et al., 2011; Malkassian, 2012). Using a Cytosense automated flow cytometer, equipped with the “image-in-flow” device, abundances and optical properties were measured every hour for species ranging from pico- to small microphytoplankton. Pictures of microphytoplanktonic cells were used to estimate cell size from the sideward scatter signal and apply the full size distribution model developed by Sosik et al. (2003). *In situ* population growth rates were eventually determined and investigated with regards to the environmental variables monitored at the same frequency.

# MATERIALS AND METHODS

## STUDY AREA

The Berre lagoon is located on the SE French coast, close to the city of Marseille (800,000 inhabitants) and the Frioul archipelago where meteorological data (wind speed and direction) is permanently recorded by the meteorological station of the Mediterranean Institute of Oceanography (MIO). It is one of the largest brackish water pools off the Mediterranean shore (with a surface of 155 km²), receiving both seawater via the Caronte channel and freshwater inputs from natural and drifted tributaries (Arc, Touloubre and Durance rivers). It is divided into two shallow sub-basins: the “Grand Etang” (max. depth 9 m) and the “Etang de Vaïne” (max. depth 4 m) where the experiment took place. Every month, the GIPREB assesses the vertical structure of the water column at 10 hydrological stations, with one located in the “Etang de Vaïne.” Profiles of temperature, salinity, dissolved oxygen, pH, chlorophyll, turbidity and nutrients

($NO_3^-$, $NO_2^-$, $NH_4^+$, $PO_4^{3-}$) are measured at several depths to provide information on the stratification of the lagoon.

## SAMPLING SYSTEM

A multi-instrumental platform composed of a Cytosense automated flow cytometer and hydrological sensors was set up in the GIPREB laboratory located on the “Berre l’Etang” harbor (Figure 1). Water was pumped directly from the Vaïne lagoon (43°47 N, 5°17 E, 2.5 m depth) to the laboratory using a 250 m pipe (inner diameter 50 mm). Samples for flow cytometric analyses were taken every hour from a 1-dm³ sampling reservoir and hydrological sensors were placed in an 80-dm³ tank supplied after the subvolume. Every hour, water was pumped (JABASCO pump) at a flow rate of 30 dm³·min⁻¹ for 17 min so that the pipe, the sampling reservoir and the tank were flushed and renewed several times prior to measurements. Samples were analyzed using the Cytosense 2 min after the pump stopped, when all air bubbles had disappeared. An evacuation pipe led the water back to the lagoon (by gravity). The JABASCO pump uses flexible impeller technology so that it does not squeeze the water passing through, therefore avoiding the damage to the cells (Thyssen et al., 2009).

Between October 1st and 14th 2011, nitrate ($NO_3^-$) concentration was measured using a Atlantic's ISUS V3 nitrate sensor set up in the 80-dm³ tank with a CTD (Hydrolab) which recorded temperature, chlorophyll concentration, turbidity, and salinity (conductivity) of the water. Throughout the month, the phytoplanktonic assemblage was characterized using the Cytosense of the PRECYM flow cytometry platform of the MIO (http://precym.mio.univ-amu.fr). Additional discrete samples were manually collected to determine nutrient ($NO_3^-$, $NO_2^-$, $PO_4^{3-}$) concentration and calibrate the ISUS sensor using a Technicon Autoanalyser (detection limits: 50 nM for nitrate and 20 nM for nitrite and phosphate). Ammonium ions ($NH_4^+$) were measured by spectrophotometry. An *in situ* HOBO sensor, fixed near to the pipe inlet, measured incident light intensity and temperature.

## FLOW CYTOMETRY

The Cytosense is a flow cytometer optimized to analyse phytoplankton. It has been specifically designed to analyse particles (cells, chains, colonies) of 1–800 µm diameter and up to several mm in length (chain forming cells). The sample flows in a sheath fluid to an exciting light source (a 488 nm laser beam) at a rate of 9 mm³·s⁻¹ in order to optimize the processing of cells, separation, and alignment. The inherent optical properties of the cells were recorded in full pulse shape as they crossed the laser beam: sideward light scatter (SWS) and red (FLR, 668–734 nm), orange (FLO, 601–668 nm) and yellow (FLY, 536–601 nm) fluorescences were collected on photomultiplier tubes. Forward light scatter (FWS) was collected by a PIN photodiode. Data was acquired in Log scale. The cytoUSB software (Cytobuoy b.v.) was used to control the flow cytometer and acquire the data stored on a computer. Data acquisition was triggered on the chlorophyll-induced red fluorescence (FLR) of the phytoplankton cell. The threshold of the FLR signal was set at 9 mV in order to compute optical pulse shapes of the photosynthetic cells and get rid of the background noise (heterotrophs and detritic particles). Samples were
<figure><img src="image_3.png" /><figcaption>FIGURE 1 | Location of the sampling point (black circle) and the 250 m pipe (black line) to pump the water to the laboratory, settled on the Berre l'étang harbor.</figcaption></figure>

hydrodynamically focused in the flow cytometer by the sheath fluid: just prior to the experiment, the cytometer was filled with 1 µm filtered seawater preserved with a 1% formaldehyde solution. After each analysis, the sheath fluid and the samples were mixed together and the sheath fluid was recycled by filtering through 0.45 and 0.1 µm integrated Nucleopore filters. A mix of fluorescent beads (Polysciences microspheres of 1.00 ± 0.03 µm, 1.6 ± 0.1 µm, 2.9 ± 0.1 µm, and 10.3 ± 0.4 µm) was regularly analyzed to ensure quality control, reliable analyses and calibrate cell sizes. Cells sharing similar optical fingerprints were manually grouped in clusters on a set of 2D projections (cytograms) in the Cytoclus software (Cytobuoy b.v.).

## ESTIMATION OF CELL VOLUME

The Cytosense used in this study is equipped with an “image-in-flow” device controlled by the CytoUSB software (Cytobuoy b.v.). During this October 2011 survey, cells displaying the highest red fluorescence and scatter intensities (a.u.) were manually gated within a targeted window where a maximum of 60 pictures per analysis (current limitation of the software) were taken (Figure 2). These pictures were used to identify and estimate the dimensions of the cells in the different clusters. Area, length, and width were automatically calculated after image processing of given files containing the original pictures (Figure 3A) using the Fiji software, an open source version of the ImageJ software

(Schindelin et al., 2012). The first step of the image processing was to automatically segment the picture into foreground (cells) with pixels values of 255 (black) and background with pixels of 0 (white) (Figure 3B). Noise was filtered (Figure 3C) and edges between pixels' value were selected as illustrated by Figure 3D. Finally we used the external edge as the contour of the cell (Figure 3E) to compute the measurements. Biovolume of all the pictured cells was estimated considering either spherical or ellipsoidal model according to cell shapes. An exponential function explained 84% of the variance between cell volume and SWS (Figure 4).

## ESTIMATION OF IN SITU GROWTH RATES

To estimate phytoplankton *in situ* growth rates, we used the size-structured population model described in Sosik et al. (2003) by applying the regression between cells' volume and sideward scatter intensities to all the cells within a specific cluster. In this model, cells are classified into several size classes according to their dimensions at time *t*. The number of classes, *m*, was chosen in order to cover the entire observed biovolume spectrum from $v_{min}$ to $v_{max}$. Classes were logarithmically spaced as follows:

$$ \text{For } i = 1, 2, \dots, m \quad v_i = v_{min} \cdot 2^{(i-1)\Delta v} $$

with $\Delta v$ constant.
At any time t, the number of cells in a size ranging from $v_i$ to $v_{i+1}$, noted $N_i(t)$, was projected to $t + dt$ via matrix multiplication so that:

$$N(t + dt) = A(t) \cdot N(t)$$

The elements $a_{ij}(t)$ of the matrix A correspond to the fraction of cells of class j at time t that becomes cells of the class i at time $t + dt$. Changes in cell size are linked to two phases of the cell cycle: the size growth during the interphase and the size decrease after the mitosis. In $dt$, one cell can either go to the direct next size class, that is to say from $v_i$ to $v_{i+1}$ through cellular growth, or can become a cell of half its original size after division. If in $dt$ one cell neither grows nor divides, it remains in the same size class (i.e., stasis).

Adjustments have been made in order to apply the model described in Sosik et al. (2003) for growth rate estimations from the dynamics of the clusters in this study. Basic hypothesis were adapted to find the best fit with natural data. Elements of the -1 shifted main diagonal, corresponding to the fraction of cells that grew between $t$ and $t + dt$ and reached the next size class were directly calculated from:

<figure><img src="image_4.png" /><figcaption>FIGURE 2 | Cytogram total red fluorescence vs. total sideward scatter (a.u.) with 95% confidence ellipses for the various clusters optically resolved during this study. The "image-in-flow" window used to trigger the picture acquisition is represented by the black square.</figcaption></figure>

$$a_{i+1,i}(t) = \gamma(t)$$

However, elements of the $\Delta v$ shifted main diagonal, corresponding to the fraction of cells that divided between $t$ and $t + dt$ and reached the size class of half their size, were counted from the fraction of cells that were not currently growing:

$$a_{i,i+1/\Delta v}(t) = 2 \cdot \delta_{i+1/\Delta v} \cdot (1 - \gamma(t))$$

Elements of the main diagonal remained those fractions of cells that has neither grown nor divided between $t$ and $t + dt$:

$$a_{i,i}(t) = (1 - \delta_i) \cdot (1 - \gamma(t))$$

The probability of cells growing to the next size class was given by an asymptotic function, assuming that it only depends on the light intensity necessary to photosynthesis:

$$\gamma(t) = \gamma_{max} \cdot (1 - e^{-E(t)/E^*})$$

with $E$, the light intensity, $\gamma_{max}$ and $E^*$ diel constants.

<figure><img src="image_5.png" /><figcaption>FIGURE 4 | Exponential relationship between total sideward scatter and biovolume of the cells pictured by the "image-in-flow" camera (n = 693) mounted on the Cytosense.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>FIGURE 3 | Image processing of a Gymnodinium picture performed by the Fiji software in order to extract the cell dimensions. (A) Original picture (B) Binary transformed picture (C) Filtered binary picture (D) Edges detection between pixels (E) Contour selection of the cell.</figcaption></figure>
The probability of cells dividing was given by a log-normal density function, assuming that it depends on the size spectrum of the vegetative cells only:

$$ \delta_i = a \cdot e^{-(\log v_i - \log \bar{v})^2 / 2 \cdot \log \sigma^2} $$

with $a, \bar{v}$ and $\sigma$ diel constants.

This model was chosen to ensure that only the size distribution of vegetative cells will account for the estimation of the population growth rates. Even if dinoflagellates species are haplontics, the coexistence of asexual and sexual forms is well documented and expected after sediment resuspension induced by turbulent mixing (Kremp and Heiskanen, 1999). Since growth rates are manifestations of the asexual life cycle, the incorrect use of sexual stages could inaccurately estimate division rates. However, dinoflagellate stages are generally heteromorphic (Pfiester and Anderson, 1987). Vegetative and planozygote/cyst cells can display distinct size (Von Stosch, 1973; Walker, 1982; Anderson et al., 1983; Blackburn et al., 1989; Figueroa et al., 2008) so the Gaussian function should enable either to restrain the size spectrum if sexual or alternative (G0) forms are present or to take into account the entire spectrum if not.

In fine, growth rates can be calculated according to the following formula:

$$ \mu = \ln \left( \frac{\sum_{i=1}^{m} \hat{N}_i (24)}{\sum_{i=1}^{m} N_i (0)} \right) $$

with $N(0)$ the initial observed size distribution, at $t = 00:00$ and $\hat{N}(24)$ the fitted size distribution at $t = 24:00$.

The parameters of the projection matrix A were those that minimize the sum of squared differences between observed and projected proportion distributions, w, ever 24 h. Either the entire distribution or only the median of the distribution was used to optimize the function:

$$ \sum_{t=1}^{24} \sum_{i=1}^{m} (\mathbf{w}(t) - \hat{\mathbf{w}}(t))^2 $$

## STATISTICAL ANALYSIS

Statistical analysis and modeling codes were performed using the R freeware (http://www.r-project.org). The time series of mean wind speed have been fitted using a linear regression (1 m function, CRAN) in order to select the window within which the temporal trend explained most of the variation of wind speed. Significance of trends of hydrological variables in this window have been estimated by comparing linear mixed-effects regression models (lmer function, lme4 library, CRAN) adapted to pseudoreplication of observations in time series (Pinheiro and Bates, 2000). To account for the strong serial correlation between measurements, random effects have been modeled as deviations of intercepts between factors grouping measurements equivalent to one period in each time series. Periodicity has been studied using discrete Fourier transforms (fft function, CRAN) on detrended time series. The maximum amplitude associated with the harmonic frequency equivalent to one period determined by Fourier analysis have been applied to differences between observed variables and non-linear trends fitted by local polynomial regressions

(loess function, CRAN). We compared, with a likelihood ratio test (anova function, CRAN), models estimated by maximum likelihood with fixed effects formulated as seasonal cycles with and without a linear time trend:

$$ \begin{aligned} Y(t) &= a_0 + a_1 \text{index}(t) + a_2 \sin(2t\pi) + a_3 \cos(2t\pi) + \varepsilon(t) \text{ and} \\ Y(t) &= a_4 + a_5 \sin(2t\pi) + a_6 \cos(2t\pi) + \varepsilon(t) \end{aligned} $$

with (t) scaled to one period, and $\varepsilon(t)$ residuals iid~N(0,σ²). For these models, residuals did not reveal any deviations from homoscedasticity, normality, and temporal dependency as tested by the autocorrelation function (acf function, CRAN). Correlation coefficients between environmental variables and cluster abundances were represented with a PCA (FactoMineR library, CRAN). Cross-correlation between wind speed time series and all the clusters abundance time series were performed to estimate the lag time between respective maxima. Negative lag time correspond to cluster abundance maximum prior to the wind event and inversely for positive lag time.

## RESULTS

### OPTICAL RESOLUTION OF THE FLOW CYTOMETRY CLUSTERS

Up to 12 flow cytometry clusters of phytoplanktonic cells (arbitrary labeled C1–C12) have been resolved over the sampling period. Each cell contributed to a particular cluster on the basis of the distinct optical properties recorded. Clusters showed a proportional relationship between SWS (correlated to cell biovolume, 0.82, $p < 0.001$) and red fluorescence emitted by chlorophyll pigments. Almost all the clusters were well discriminated on red fluorescence vs. sideward scatter cytogram (Figure 2). With the same level of red fluorescence, clusters C5 and C8 emitted more orange fluorescence than clusters C6 and C9, respectively. The size spectrum of the cells ranged within the pico-, nano-, and microphytoplankton size classes. Smallest cells were observed in the cluster C12 (mean $0.9 \pm 0.1$ µm) and the biggest cells belonged to cluster C1 (mean $56.4 \pm 12.2$ µm). Clusters C1–C4 were composed of different dinoflagellates species identified as Akashiwo sanguinea (Hirasaka), Prorocentrum micans (Erhenberg), Scrippsiella sp. (Balech) and Gymnodinium sp. (Stein), respectively over the all sampling period. Clusters C5–C10 were assigned to nanophytoplankton with no visual identification possible by the “image-in-flow” pictures due to the lack of resolution of the pictures. Picophytoplankton was represented by clusters C11 and C12 with mean length estimated by the calibrated beads at $1.1 \pm 0.1$ µm and $0.9 \pm 0.1$ µm.

### PHYTOPLANKTONIC COMMUNITY DYNAMICS

#### Environmental variables

During October 2011, the wind was the main forcing on the hydrological state of the Berre lagoon and more particularly the intense Mistral blowing from the North (330–360°), with mean speeds exceeding 20 m/s (Figure 5). Turbulent mixing of the water column increased as gusts reached $18.7 \pm 3.2$ m/s between October 7th and 9th (Figure 6A). During these 2 days, all the hydrological variables showed significant trends with the exception of in situ light intensity [$\chi^2_{(1)}, n = 79, p < 0.05$]. Mean water temperature measured at the sampling point cooled by 5.2°C
<figure><img src="image_7.png" /><figcaption>FIGURE 5 | Distribution of wind direction during October 2011 with gusts exceeding 20 m/s.</figcaption></figure>

between October 05th and 11th (Figure 6B). Salinity distribu-
tion showed a wider range of variation during the Mistral event,
suggesting a gradient between surface and bottom layers prior
to the wind: difference between 0.05 and 0.95 percentiles was
0.19 on October 6th and 0.59 on October 8th. On September
16th and October 18th, the nutrient vertical profiles measured by
the GIPREB at surface and bottom depths displayed higher con-
centrations of recycled nutrients (NH<sub>4</sub><sup>+</sup>, PO<sub>4</sub><sup>3-</sup>) at the bottom
of the lagoon, suggesting the stratification of the water column
prior to and after the first Mistral event (Figure 7). Variation
in nitrate concentrations was similar to temperature (Figure 6C)
with a negative trend of -0.04 µM per hour [χ<sup>2</sup><sub>(1)</sub> = 17.2, n = 79,
p < 0.001] which led to a minimum of 0.29 µM on October 9th at
19:00. Chlorophyll concentration reached 20.96 µg/l on October
6th 13:00 (Figure 6D). After the maximal gust of 25 m/s, the mean
value dropped to 3.97 µg/l. Despite the decrease in chlorophyll
concentration the *in situ* light intensity was not significantly dif-
ferent before or during the Mistral event [χ<sup>2</sup><sub>(1)</sub> = 0.07, n = 79,
p = 0.8]. The photoperiod was between 08:00 and 19:00 and the
light intensity was always maximal at 12:00 (Figure 6E).

Community response to the wind event

A Principal Component Analysis was performed to analyze the
multivariate time series of hydrological variables and cluster
abundances. The main relationships between phytoplankton con-
centrations and the environment were well represented on the
two first components, which explained 59.1% of the total vari-
ance. High chlorophyll concentration was observed simultane-
ously with the occurrence of higher abundances of clusters C1–C9
(Figure 8). The increase of chlorophyll concentration and turbid-
ity in the Berre lagoon was phased with the sudden proliferation

<figure><img src="image_8.png" /><figcaption>FIGURE 6 | Dynamics of the hydrological variables before, during and after the strong wind (Mistral) event. (A) Wind speed, (B) Temperature, (C) NO₃⁻ concentration, (D) Chlorophyll concentration, and (E) Light intensity.</figcaption></figure>

of clusters C1–C3, C5, and C8 (correlation coefficient 0.36, 0.46,
0.31, 0.50, and 0.51, respectively, *p* < 0.001), as well as the pres-
ence of C9 cells. C9 was the most abundant cluster with mean
concentration 17,024 ± 10,880 cells/ml between October 05th
and 06th and a maximum of 51,142 cells/ml on October 05th
20:00. All of these clusters were more abundant before the intensi-
fication of the Mistral wind speed (Table 1). The lag time between
maximum wind speed and abundances mainly ranged between
24 and 48 h. Turbidity was correlated to high temperature (0.63,
*p* < 0.001) and high nitrate concentration (0.26, *p* < 0.001).
Conversely, picophytoplankton cluster (C11 and C12) concentra-
tions increased less than 1 h after the wind speed rose (correlation
coefficient 0.72 and 0.60, *p* < 0.001). Coefficients of cross corre-
lation were thus maximal for the 0 h lag time between wind speed
and abundance time series. Mean abundance was 1482 ± 1009
cells/ml for C11 and 923 ± 317 cells/ml for C12 between October
05th and 06th and reached 3787 ± 1557 cells/ml and 1690 ± 411
cells/ml during October 07th and 08th. The increase in these clus-
ter concentrations was reproduced during a second strong Mistral
event (Figure 9), with an incremented amplitude. The abundance
of the *Gymnodinium* cluster did not significantly change between
October 05th and 11th. The 3rd and 4th principal component
that explained 8.7% and 6.7% of the total variance respectively
(data not shown) showed the negative correlation between light
<figure><img src="image_9.png" /><figcaption>FIGURE 7 | Vertical distribution of (A) ammonium and (B) phosphates concentrations measured by the GIPREB at surface and bottom layers at the hydrological station in the Vaïne lagoon in 2011.</figcaption></figure>

<table><caption>Table 1 | Maximum cross correlation coefficient and associated lag time of wind speed and cluster abundance time series between October 5th 00:00 and October 11th 23:00.</caption><thead><tr><th>Clusters</th><th>C1</th><th>C2</th><th>C3</th><th>C4</th><th>C5</th><th>C6</th><th>C7</th><th>C8</th><th>C9</th><th>C10</th><th>C11</th><th>C12</th></tr></thead><tbody><tr><td>Coefficient</td><td>0.57</td><td>0.62</td><td>0.62</td><td>0.48</td><td>0.62</td><td>0.60</td><td>0.63</td><td>0.73</td><td>0.57</td><td>0.46</td><td>0.96</td><td>0.92</td></tr><tr><td>Lag (h)</td><td>-32</td><td>-30</td><td>-28</td><td>132</td><td>-30</td><td>-39</td><td>-14</td><td>-26</td><td>-48</td><td>49</td><td>0</td><td>0</td></tr></tbody></table>

intensity and turbidity (–0.14, *p* < 0.1) and the positive influence of light on chlorophyll concentration mostly before the Mistral intensification (0.36, *p* < 0.1).

## DYNAMICS OF THE PHYTOPLANKTONIC POPULATIONS’

### Gymnodinium dynamics

From the pictures taken by the “image-in-flow” device, cluster C4 has been identified as *Gymnodinium* sp. Its dynamics was quasi-periodic during the Mistral event (Figure 10A). The concentration was marked by several peaks of abundance with a mean of 49.9 ± 28.8 cells/ml. The highest concentration reached 142.4 cells/ml on October 6th at 19:00 and the abundance regularly dropped down below 50 cells/ml. Between October 05th and 09th, a high abundance was generally observed at the end of the day, but after the lull in the wind two distinct peaks were observed each day. On October 10th, maximal concentrations (of 122.8 and 79.2 cells/ml) were measured at 05:00 and 13:00, respectively. On the 11th, they were measured at 03:00 and 13:00, with 82.8 and 71.6 cells/ml, respectively. The same day, the maximum mean biovolume of the cells was 3 times greater than the minimum, passing from 4852–16,294 µm³ (Figure 10B). During this period, cells were highly phased within cluster C4, making the estimation of growth rates based on the biovolume distribution

more reliable (Figure 11). To apply the model and compare the observed (Figure 11A) and theoretical (Figure 11B) distribution of the biovolume, $v_{min}$ and $v_{max}$ were set to 3785 µm³ and 59,874 µm³, dt to 10 min, $\Delta v$ to 1/4 and $m$ to 17. The mean cell biovolume increased 1 h after dawn until approximately 20:00 (Figure 11D) before it progressively decreased due to mitotic divisions (Figure 11C). The amplitude of these two processes varied before, during and after the wind event so that growth rates ranged between 0.22 d⁻¹ (<1 division per day) and 0.85 d⁻¹ (> 1 division per day) (Figure 10C).

### Dynamics of cluster C5

Cells forming the cluster C5 could not be identified on the pictures. It was one of the clusters that showed a sudden peak of abundance on October 6th, reaching up to 838 cells/ml at 15:00 (Figure 12A). Two hours later, the concentration was 101 cells/ml and remained at a mean value of 85.9 ± 57.7 cells/ml. The cells were well phased with the light/dark period between October 5th and 11th as indicated by the mean cell biovolume that exhibited a clear diel variation in relation to the cell cycle (Figure 12B). Cell size began to increase 1 h after dawn (09:00), until a few hours after dusk (20:00–23:00). The ratio of max/min mean biovolume ranged from 1.26 on October 7th to 1.55 on October 11th.
<figure><img src="image_10.png" /></figure>

**FIGURE 8** | Principal component analysis performed on the dataset. Projection of the observations (gray points, i.e., dates from 10/05/11 00:00 to 10/11/11 23:00) and the variables (arrows, i.e., hydrological variables and clusters abundance) on the two first principal components of the PCA (59.1 % of the total variance).

Between these dates, growth rates progressively increased with regard to the variation of mean biovolume amplitude and reached 0.46 $d^{-1}$ on October 11th (Figure 12C). To apply the model to cluster C5 and compare the observed and theoretical distribution of the biovolume, $v_{min}$ and $v_{max}$ were set to 3707.6 $\mu m^3$ and 16,017.6 $\mu m^3$, $dt$ to 10 min, $\Delta v$ to 1/9 and $m$ to 20.

# DISCUSSION

## PHYTOPLANKTON IN THE BERRE LAGOON

Because of high nutrients inputs and regeneration (Gouze et al., 2008), the primary production in the Berre lagoon has always been high compared to others Mediterranean lagoons (Minas, 1976; Kim, 1983). High production rates leading to blooms combined with great diversity have made this site a unique area for phytoplankton study.

In October, these blooms have been observed due to persistent summer-like temperatures, causing a thermal stratification of the water column, and the resurgence of nutrients supply by the power plant activities. Before the freshwater discharges were regulated by the power plant, algal proliferation was almost monospecifically dominated by red tide dinoflagellates such as *Prorocentrum minimum*. However, the wind regimes occurring in this region play a major role in mixing and oxygenating the water column. This limits the negative impact of hypoxia that could be induced by the heterotrophic activity. The Mistral is a strong Northerly wind that generally drives the superficial layers of the Berre lagoon toward the Caronte channel to join the Mediterranean Sea (Ulises et al., 2005). It causes the homogenization of the water mass, mixing warm surface layer and colder

<figure><img src="image_11.png" /></figure>

**FIGURE 9** | Time series of the (A) C11 cluster and (B) C12 cluster abundances (black line) during the sampling period, characterized by two strong Mistral event (gray dotted line).

bottom layer (Nerini et al., 2001). All these forcings drove the phytoplanktonic community to evolve in order to develop within the range of temperature and salinity gradients induced by wind, rivers and freshwater discharges. However, since the regulation of discharge imposed on the hydroelectric plant, species richness has increased (Raimbault et al., 2013). Diatoms (*Chaetoceros* sp., *Pseudo-Nitzschia* sp., and *Thalassionema* sp.), dinoflagellates have been sampled at 10 hydrological stations and identified by microscopy and flow cytometry (data not shown), technique recently implemented in the monitoring strategy of the GIPREB. Flow cytometric analyses by an autonomous *in situ* flow cytometer have revealed more phytoplankton clusters in the lagoon than in the Bay of Marseille, in the NW Mediterranean Sea (Thyssen et al., 2011). In 2005–2006, 7 clusters were resolved in the Bay of Marseille assigned to pico (2 clusters) and nanophytoplankton (5 clusters) compared to 12 clusters identified in 2011 in the Berre lagoon at the sampling point. Two clusters were assigned to picophytoplankton, 6 clusters of nanophytoplankton and 4 of microphytoplankton. In the microphytoplanktonic clusters, cells were identified as red tide dinoflagellates species such as *Akashiwo sanguinea*, *Prorocentrum micans*, *Gymnodinium* sp., and *Scrippsiella* sp. The first two are eurythermale and euryhaline species, able to grow within temperature (10–30°C) and salinity (10–40) gradients (Matsubara et al., 2007; Dhib et al., 2013)

## COMMUNITY AND PRODUCTION CHANGES TO A STRONG MISTRAL EVENT IN THE BERRE LAGOON

Offsets of phytoplanktonic production are calculated from the loss and growth rates of natural populations. Loss in net productivity at a fixed point can result from trophic interactions
<figure><img src="image_12.png" /><figcaption>**FIGURE 10** | Time series of the *Gymnodinium* cluster during the strong Mistral event. (A) Population dynamics. (B) Variation of the mean cells biovolume. (C) Estimations of the *in situ* growth rate provided by the ratio</figcaption></figure>

between maximal and minimal mean diel biovolume (gray points) and the model (black points). Error intervals correspond to growth rate estimates after applying the biovolume regression model with standard errors.

such as grazing and viral lysis as well as physical transport (Riley et al., 1949). For our experiment, the sampling strategy for high frequency characterization of the phytoplankton community has been limited to one depth at a fixed station, and thus represents only a partial view of the response of phytoplankton clusters in the lagoon. However, several 3D hydrodynamical models simulate the circulation in the station location as confined within the Vaïne lagoon in case of a strong Mistral event (Leredde et al., 2002; Alekseenko et al., 2013). Moderate turbulent mixing is usually beneficial for large non-motile cells of phytoplankton because it counterbalances their sedimentation and prevents their exportation toward the aphotic layer (Margalef, 1978). The classical succession scheme induced by temporal or spatial ergoclines (Legendre and Demers, 1985) predicts the outbreak of diatoms followed by dinoflagellates when nutrients are depleted in the surface layers. Because of both inputs of freshwater and seawater, maintaining the haline gradient between surface and bottom layers, the deepest part of the Berre lagoon is almost permanently stratified (Nerini et al., 2001). Only strong gusts of wind (>10 m/s) allow the homonogeneization of the water column. In October 2011, automated flow cytometry analyses

prior the strong Mistral event revealed the sudden proliferation of red tide dinoflagellates as well as nanophytoplankton clusters. Dinoflagellate species regroup complex microorganisms. Their life cycle implying sexual, asexual, and resting stages as well as their trophic strategy can both justify the potential of dinoflagellates to outcompete other phytoplanktonic species (Schnepf and Elbrachter, 1992; Kremp, 2013). When inorganic compounds are not sufficient in their environment, dinoflagellates can swim in order to have access to the deeper nutrient pool or ingest small prey in food vacuoles if they are mixotrophs (Stoecker, 1999). The proliferation of dinoflagellates during our experiment thus supports the idea that they are K strategists relying on the microbial loop for the uptake of regenerated nutrients (Smayda, 1997; Pitcher et al., 1998; Kudela et al., 2008). On the contrary, high frequency sampling revealed that picophytoplankton clusters might show a reactive capacity superior to other groups such as diatoms, and so be defined as pioneer species. Autonomous sampling has already indicated that picophytoplankton could be the first to develop following a strong wind event (Thyssen et al., 2008). Higher assimilation rates give these small cells an advantage compared to larger species accounting for new production as well
<figure><img src="image_13.png" /><figcaption>**FIGURE 11** | Output of the size structured population model applied to cluster C4 on October 11th 2011: Diel variations of the (A) observed and (B) theoretical distributions of cell proportions in each size class. The</figcaption></figure>

theoretical distribution was estimated only on the basis of two cell cycle processes, division and size growth with probability depending on (C) the biovolume of the cells and on (D) the incident light intensity respectively.

(Eppley et al., 1969; Aknes and Egge, 1991). In theory, the concentration of such r strategists should increase within a short lag time after the perturbation. In this study, clusters C11 and C12 started to multiply less than 1 h after the wind speed increase. This behavior was reproduced during a second strong Mistral event, which also favored their proliferation soon after the weather lull, but likely prevented the reappearance of the following successors, accounting for regenerated production fuelled by the organic matter recycling (Legendre and LeFèvre, 1989).

## GROWTH RATE ESTIMATIONS

As sideward scatter is strongly correlated to cell size (Ackleson and Spinrad, 1988; Simon et al., 1994), determination of cell optical properties can provide information about the life cycle of micro-organisms (Kiefer et al., 1979). The model developed by Sosik et al. is based on the full size distribution that varies during the sequential phases of the cell cycle. Therefore, cells will grow with the support of photosynthesis and divide according to their size. In an attempt to accommodate these simple assumptions to clusters observed in this study, we tested an alternative model for division, relaxing the hypothesis of the strict increase of division rates with biovolume. However, caution should be obviously taken with *in situ* growth rates estimations of dinoflagellate populations, particularly with regard to life stages differentiations (Cetta and Anderson, 1990; Gisselson et al., 1999). No inhibitory effect of high light intensity or UV exposure was expressed in the growth model since the daily red fluorescence variation did not exhibit any depression during the day (data not shown), as can be observed when quenching occurs (Jacquet et al., 1998). This model fitted C4 and C5 clusters size distribution patterns that were, as for dinoflagellates (Sweeney, 1959), phased to the photocycle. Cellular growth began ~1 h after dawn and continued

briefly after dawn, whereas division mainly occurred during the night, in agreement with the literature (Weiler and Chisholm, 1976). As a greater proportion of cells was in the G1 phase, we changed the dependence relation between division and growth to reflect that division mainly occurred when cellular growth was no more dominant. After the Mistral event, cells reached higher dimensions during the photoperiod, suggesting a possible light limitation. Since duration of the interphase and mitotic divisions was unchanged, growth rates naturally increased. Cells in the *Gymnodinium* cluster underwent from less than one division to more than one division per day after October 9th. The maximal growth rate was at 0.85 d⁻¹ on October 10th. Growth rates of the C5 cluster were always inferior to 0.69 d⁻¹ but they progressively rose to 0.57 d⁻¹.

Clusters C1–C3 were not abundant enough relatively to the volume analyzed to ensure a reliable application of the model. The size variation did not show any clear pattern phased to the photoperiod. The shape of size histograms was constant during the day, giving evidence of asynchronous populations (Campbell and Yentsch, 1989). Asynchrony of dinoflagellates appears when spatial migration is triggered by physiological needs, local current, cyst bed or wind rather than by the light/dark cycle (Ralston et al., 2007). Nutritional migration strategy over photo/geotaxis in natural populations leads to deep chlorophyll maximums which coincide with the primary production, near the nutricline (Cullen, 1982). However, synergy between nutrient availability and irradiance needs depends not only on hydrodynamics but also on internal cells’ nutritional status, reinforcing the variability within populations (Ji and Franks, 2007). In some particular ecosystems, the vertical distribution of dinoflagellates is typically bimodal, with a second abundance peak at the surface (Townsend et al., 2001). Even if diel size variation can be persistent along
<figure><img src="image_14.png" /><figcaption>**FIGURE 12** | Time series of the cluster C5 during the strong Mistral event. (A) Population dynamics. (B) Variation of the mean cells biovolume. (C) Estimations of the *in situ* growth rate provided by the ratio between</figcaption></figure>

maximal and minimal mean diel biovolume (gray points) and the model (black points). Error intervals correspond to growth rate estimates after applying the biovolume regression model with standard errors.

the water column, growth rates of population necessarily vary within layers. For populations synchronized to the photoperiod, difference in incident irradiance with depth can induce the retardation of the cell cycle whereas nutrient stratification and pigment quenching due to high UV can also completely invert its timing between deep and surface layers (Vaulot and Marie, 1999). Water mixing could thus influence the growth rate of a patchy population, as it should homogenize the nutrient pool and light availability. The vertical gradient of cell size can lead to growth rate differences between model and mean biovolume estimates. Presence of non-vegetative larger cells (Blackburn et al., 1989), displaying smaller division probability, could explain the biovolume ratio being superior to the model ratio when mean size is used. Stratification of the water column, as demonstrated by the time series of hydrological variables recorded, should increase the differences between these two estimations. Although there is strong evidence of the inhibitory effects of turbulence on dinoflagellate cell division and growth (Berdalet, 1992), the rise of C4–C5 cluster growth rates at the sampling point with the wind was likely due to higher growth rates at the bottom, suggesting

the presence of a nutricline before the mixing, and a decrease in the light limitation due to turbidity dispersal. The influence of drop of temperature on growth rates is improper with no or positive correlation in the range 15–20°C for phytoplankton (Eppley, 1972) and particularly for species of the genus *Gymnodinium* (Thomas, 1975; Nielsen, 1996; Yamamoto et al., 2002).

## CONCLUSION

In this study abundances and optical properties of phytoplankton community were monitored in autumn in one of the biggest European brackish lagoons using an autonomous flow cytometer especially designed for photosynthetic microorganisms. This *in situ* monitoring carried out in the Berre lagoon and remotely operated (thanks to internet connection) enabled to follow at the single cell level the dynamics of a large diversity of phytoplanktonic functional groups, as defined by flow cytometry, in response to brief and sudden environmental forcings. In October, the concentrations of 2 groups of picophytoplankton, 6 groups of nanophytoplankton, and 4 groups of microphytoplankton were measured at high frequency before, during, and after a
strong wind event, providing in near real time the *in situ* struc-
ture of the phytoplankton community. Data analyses revealed
a high variability of the autotrophic community in response to
this natural forcing, that has been undetected by the routine
survey conducted monthly by a conventional manual sampling.
Phytoplankton responses to wind have also been reproduced dur-
ing a second event and have been previously described in the
literature, leading to the definition of the concept of response
functional groups (Thyssen et al., 2014).

Besides, a potential application to harmful algal bloom moni-
toring is conceivable with this technology combining flow cytom-
etry and image analysis (Campbell et al., 2010). In this study,
we managed to identify red tide species and their dynamics at
population and intrinsic levels. A stable hydrological environ-
ment triggered their proliferation before the wind event, result-
ing in a chlorophyll biomass peak, with the exception of the
*Gymnodinium* cluster. The turbulent mixing generally affects the
spatial distribution of phytoplankton by physical transport and
may hide influence on the net abundance of populations by phys-
iological changes that need to be measured at the individual
level. To address this issue, estimates of *in situ* population growth
rates have been calculated using the light scatter signal intensity
recorded during the passage of each cells through the 488 nm laser
beam. This signal was strongly correlated to the biovolume of the
cells pictured by the “image-in-flow” device. Growth rates were
calculated from the diel changes of size distribution monitored
by the instrument thanks to the high frequency measurements
that are difficult to obtain using more conventional methods.
The increase of asexual reproductive rates was predicted by the
model of Sosik et al. (2003) for the two phased clusters, proba-
bly as a result of the attenuation of light/nutrient limitation at the
sampling point in response to the homogenization of the water
column.

Although chlorophyll biomass dropped during the wind-
induced mixing, the evolution of population growth rates suggest
that in some cases, the calculation of NPP based on ocean color
might not match the estimation of photosynthetic carbon fixa-
tion. Observing marine microorganisms at the single cell level
with an automated flow cytometer operating at high frequency
(up to several sampling and analyses per hour) takes into account
the short-term variability of phytoplankton. Coupled to other
hydrological sensors also operating at high frequency it should
bring a new insight into the phytoplankton structure and func-
tioning. It should also fill the gap in the primary production
budget estimations based mostly on low frequency sampling (typ-
ically every month or twice a month) that cannot take into
consideration pulsed events. Not taking into account the fast
biomass pulses of phytoplankton could significantly impact the
estimations of biogeochemical fluxes and budgets on an annual
scale (Lomas et al., 2009).

ACKNOWLEDGMENTS

The instatement of the automated platform near the Berre lagoon
was supported by the project EC2CO-DRIL (CNRS) of “In situ
Monitoring of phytoplanktonic assemblages dynamics” and the
PRECYM flow cytometry platform of the MIO for provid-
ing the Cytosense flow cytometer. We thank the GIPREB for

logistic support and Beatriz Beker for the identification of phy-
toplankton species. We are grateful to Mathieu Fallet for his
contribution to the images treatment. M. Dugenne benefits from
a fellowship of the French Minister of High Education and
Research for a phD thesis supervised by Gérald Grégori and Denys
Pommeret.

REFERENCES

Ackleson, S. G., and Spinrad, R. W. (1988). Size and refractive index of individual marine particulates: a flow cytometric approach. *Appl. Opt.* 27, 1270–1277. doi: 10.1364/AO.27.001270

Aknes, D. L., and Egge, J. K. (1991). A theoretical model for nutrient uptake in phytoplankton. *Mar. Ecol. Prog. Ser.* 70, 65–72. doi: 10.3354/meps070065

Alekseenko, E., Roux, B., Sukhinov, A., Kotarba, R., and Fougere, D. (2013). Coastal hydrodynamics in a windy lagoon. *Comput. Fluids* 77, 24–35. doi: 10.1016/j.compfluid.2013.02.003

Anderson, D. M., Chisholm, S. W., and Watras, C. J. (1983). Importance of life cycle events in the population dynamics of *Gonyaulax tamarensis*. *Mar. Biol.* 76, 179–189. doi: 10.1007/BF00392734

Antoine, D., Andre, J. M., and Morel, L. (1996). Oceanic primary production. Estimation at global scale from satellite (coastal zone color scanner) chlorophyll. *Global Biogeochem. Cycles* 10, 57–69. doi: 10.1029/95GB02832

Berdalet, E. (1992). Effects of turbulence on the marine dinoflagellates *Gymnodinium nelsonii*. *J. Phycol.* 28, 267–272. doi: 10.1111/j.0022-3646.1992.00267.x

Binder, B. J., Chisholm, S. W., Olson, R. J., Frankel, S. L., and Worden, A. Z. (1996). Dynamics of picophytoplankton, ultraphytoplankton and bacteria in the central equatorial Pacific. *Deep-sea Res. II, Top. Stud. Oceanogr.* 43, 907–931. doi: 10.1016/0967-0645(96)00023-9

Blackburn, S. I., Hallegraeff, G. M., and Bolch, C. J. (1989). Vegetative repro-
duction and sexual life cycle of the toxic dinoflagellate *Gymnodinium cate-
natum* from Tasmania, Australia. *J. Phycol.* 25, 577–590. doi: 10.1111/j.1529-
8817.1989.tb00264.x

Campbell, J. W., and Yentsch, C. M. (1989). Variance within homogeneous phytoplankton populations, II: analysis of clonal cultures. *Cytometry* 10, 596–604. doi: 10.1002/cyto.990100515

Campbell, L., Olson, R. J., Sosik, H. M., Araham, A., Henrichs, D. W., Hyatt, C. J., et al. (2010). First harmful Dinophysis (Dinophyceae, Dinophysiales) bloom in the U.S is revealed by automated imaging flow cytometry. *J. Phycol.* 46, 66–75. doi: 10.1111/j.1529-8817.2009.00791.x

Cetta, C. M., and Anderson, D. M. (1990). Cell cycle studies of the dinoflagellates *Gonyaulax polydra* Stein and *Gyrodinium catenum* Hulburt during sexual and asexual reproduction. *J. Exp. Mar. Biol. Ecol.* 135, 69–84. doi: 10.1016/0022-0981(90)90199-M

Cloern, J. E. (1996). Phytoplankton bloom dynamics in coastal ecosystems: a review with some general lessons from sustained investigation of San Francisco Bay, California. *Rev. Geophys.* 34, 127–168. doi: 10.1029/96RG00986

Cullen, J. J. (1982). The deep Chlorophyll maximum: comparing vertical profiles of Chlorophyll a. *Can. J. Fish. Aquat. Sci.* 39, 791–803. doi: 10.1139/f82-108

Dhib, A., Frossard, V., Turki, S., and Aleya, L. (2013). Dynamics of harmful dinoflagellates driven by temperature and salinity in a northeastern Mediterranean lagoon. *Environ. Monit. Assess.* 185, 3369–3382. doi: 10.1007/s10661-012-2797-4

Dubelaar, G. B. J., Geerders, P. J. F., and Jonker, R. R. (2004). High frequency monitoring reveals phytoplankton dynamics. *J. Environ. Monit.* 6, 946–952. doi: 10.1039/b409350j

Dubelaar, G. B. J., and Gerritzen, P. L. (2000). CytoBuoy: a step forward towards using flow cytometry in operational oceanography. *Sci. Mar.* 64, 255–265.

Dubelaar, G. B. J., and Jonker, R. R. (2000). Flow cytometry as a tool for the study of phytoplankton. *Sci. Mar.* 64, 135–156. doi: 10.3989/scimar.2000.64n2135

Durand, M. D. (1995). *Phytoplankton Growth and Diel Variations in Beam Attenuation Through Individual Cell Analysis.* Thèse de doctorat, MIT/WHOI.

Eppley, R. W. (1972). Temperature and phytoplankton growth in the sea. *Fish. Bull.* 70, 1063–1085.

Eppley, R. W., Rogers, J. N., and McCarthy, J. J. (1969). Half-saturation constants for uptake of nitrate and ammonium by marine phytoplankton. *Limnol. Oceanogr.* 14, 912–920 doi: 10.4319/lo.1969.14.6.0912
Field, C. B., Behrenfield, M. J., Randerson, J. T., and Falkowski, P. G. (1998). Primary production of the biosphere: integrating terrestrial and oceanic components. *Science* 281, 237–240. doi: 10.1126/science.281.5374.237

Figueroa, R. I., Bravo, I., Ramilo, I., Pazos, Y., and Morono, A. (2008). New life-cycle stages of *Gymnodinium catenatum* (Dinophyceae): laboratory and field observations. *Aquat. Microb. Ecol.* 52, 13–23. doi: 10.3354/ame01206

Gisselson, L. A., Graneli, E., and Carlsson, P. (1999). Using cell cycle analysis to estimate *in situ* growth rates of the dinoflagellates *Dinophysis acuminate*: drawbacks of the DNA quantification method. *Mar. Ecol. Prog. Ser.* 184, 55–62. doi: 10.3354/meps184055

Goldman, J. C., McCarthy, J. J., and Peavey, D. G. (1979). Growth rate influence on the chemical composition of phytoplankton in oceanic waters. *Nature* 279, 210–215. doi: 10.1038/279210a0

Gouze, E., Raimbault, P., Garcia, N., Bernard, G., and Picon, P. (2008). Nutrient and suspended matter discharge by tributaries into the Berre Lagoon (France): the contribution of flood events to the matter budget. *C. R. Geosci.* 340, 233–244. doi: 10.1016/j.crte.2007.12.007

Hutchinson, G. P. (1961). The paradox of the plankton. *Am. Nat.* 95, 137–145. doi: 10.1086/282171

Jacquet, S., Lennon, J. F., Marie, D., and Vaulot, D. (1998). Picoplankton population dynamics in coastal waters of the Northwestern Mediterranean Sea. *Limnol. Oceanogr.* 43, 1916–1931.

Ji, R., and Franks, P. J. S. (2007). Vertical migration of dinoflagellates: model analysis of strategies, growth and vertical distribution patterns. *Mar. Ecol. Prog. Ser.* 344, 49–61. doi: 10.3354/meps06952

Kiefer, D. A., Olson, R. J., and Wilson, W. H. (1979). Reflectance spectroscopy of marine phytoplankton. Optical properties as related to age and growth rates. *Limnol. Oceanogr.* 24, 664–672. doi: 10.4319/lo.1979.24.4.0664

Kim, K. T. (1983). Production primaire pelagique de l'Etang de Berre en 1977 et 1978. Comparaison avec le milieu marin (Méditerranée nord-occidentale). *Mar. Biol.* 73, 325–341. doi: 10.1007/BF00392258

Kremp, A. (2013). "Diversity of dinoflagellates life cycles: facets and implications of complex strategies," in *Biological and Geological Perspectives of Dinoflagellates*, eds J. Lewis, F. Marret, and L. Bradley (London: The Micropalaeontological Society; Special Publications; Geological Society), 197–205.

Kremp, A., and Heiskanen, A. S. (1999). Sexuality and cyst formation of the spring-bloom dinoflagellate *Scrippsiella hangoei* in the coastal northern Baltic Sea. *Mar. Biol.* 134, 771–777. doi: 10.1007/s002270050594

Kudela, R. M., Lane, J. Q., and Cochlan, W. P. (2008). The potential role of anthropogenically derived nitrogen in the growth of harmful algae in California, USA. *Harmful Algae* 8, 103–110. doi: 10.1016/j.hal.2008.08.019

Legendre, L., and Demers, S. (1985). Auxiliary energy, ergoclines and aquatic biological production. *Nat. Can.* 112, 5–14.

Legendre, L., and LeFèvre, L. (1989). "Hydrodynamical singularities as controls of recycled versus export production in oceans," in *Productivity of the Ocean: Present and Past*, ed W. H. Berger (Chichester: John Wiley & Sons), 44–63.

Leredde, Y., Dekeyser, Y., and Devenon, J. L. (2002). T-S data assimilation to optimize turbulent viscosity: an application to the Berre lagoon hydrodynamics. *J. Coast. Res.* 18, 555–567.

Lomas, M. W., Roberts, N., Lipschultz, F., Krause, J. W., Nelson, D. M., and Bates, N. R. (2009). Biogeochemical responses to late-winter storms in the Sargasso Sea IV. Rapid succession of major phytoplankton groups. *Deep Sea Res. Part I Oceanogr. Res. Pap.* 56, 892–909. doi: 10.1016/j.dsr.2009.03.004

Malkassian, A. (2012). *Méthodes d'analyse Fonctionnelle et Multivariée Appliquées à l'étude du Fonctionnement Écologique des Assemblages Phytoplanctoniques de l'étang de Berre*. Thèse de Doctorat, Université d'Aix-Marseille, France.

Margalef, R. (1978). Life-forms of phytoplankton as survival alternatives in an unstable environment. *Oceanol. Acta* 1, 493–509.

Matsubara, T., Nagasoe, S., Yamasaki, Y., Shikata, T., Shimasaki, Y., Oshima, Y., et al. (2007). Effects of temperature, salinity, and irradiance on the growth of the dinoflagellate *Akashiwo sanguinea*. *J. Exp. Mar. Biol. Ecol.* 342, 226–230. doi: 10.1016/j.jembe.2006.09.013

Minas, M. (1976). Production organique primaire dans un milieu saumâtre eutrope (Etang de Berre). Effets d'une forte dilution (dérivation des eaux de la Durance). *Mar. Biol.* 35, 13–29. doi: 10.1007/BF00386671

Nerini, D., Mante, C., Durbec, J.-P., and Garcia, F. (2001). Une méthode statistique de détermination de séquences caractéristiques dans une série temporelle de plusieurs variables. Application à la physico-chimie des eaux de l'étang de Berre. *C. R. Acad. Sci. II.* 332, 457–464. doi: 10.1016/S1251-8050(01)01561-0

Nielsen, M. V. (1996). Growth and chemical composition of the toxic dinoflagellate *Gymnodinium galateanum* in relation to irradiance, temperature and salinity. *Mar. Ecol. Prog. Ser.* 136, 205–211. doi: 10.3354/meps136205

Olson, R. J., Shalapyonok, A., and Sosik, H. M. (2003). An automated submersible flow cytometer for analyzing pico- and nanophytoplankton: FlowCytobot. *Deep Sea Res. Part I Oceanogr. Res. Pap.* 50, 301–315. doi: 10.1016/S0967-0637(03)00003-7

Peeters, J. C. H., Dubelaar, G. B. J., Ringelberg, J., and Visser, J. W. M. (1989). Optical plankton analyser: a flow cytometer for plankton analysis, 1: design considerations. *Cytometry* 10, 522–528. doi: 10.1002/cyto.990100507

Pfiester, L. A., and Anderson, D. M. (1987). "Dinoflagellates reproduction," in *Botanical Monographs*, Chapter 14, Vol. 21, ed F. J. R. Taylor (Oxford: Blackwell Scientific Publications), 611–648.

Pinheiro, J. C., and Bates, D. M. (2000). *Mixed-Effects Models in S and S-Plus*. New York, NY: Springer-Verlag.

Pitcher, G. C., Boyd, A. J., Horstman, D. A., and Mitchell-Innes, B. A. (1998). Subsurface dinoflagellate populations, frontal blooms and the formation of red tide in the southern Benguela upwelling system. *Mar. Ecol. Prog. Ser.* 172, 253–264. doi: 10.3354/meps172253

Raimbault, P., Beker, B., Garcia, N., Fornier, M., and Gouze, E. (2013). "L'étang de Berre, cycle de la matière, production primaire et communautés phytoplanctoniques," in *LAGUN'R*, ed GIPREB (Éang de Berre).

Ralston, D. K., McGillicuddy, D. J., and Townsend, D. W. (2007). Asynchronous vertical migration and bimodal distribution of motile phytoplankton. *J. Plankton Res.* 29, 803–821. doi: 10.1093/plankt/fbm061

Rantajärvi, E., Olson, R., Hälfors, S., Leppänen, J. M., and Raateoja, M. (1998). Effect of sampling frequency on detection of natural variability of phytoplankton: unattended high-frequency measurements on board ferries in the Baltic Sea. *J. Mar. Sci.* 55, 697–704.

Riley, G. A., Stommel, H., and Bumpus, D. F. (1949). Quantitative ecology of the plankton of the Western-North Atlantic. *Bull. Bingham Oceanogr. Collect.* 12, 1–169.

Schindelin, J., Arganda-Carreras, I., Frise, E., Kaynig, V., Longair, M., Pietzsch, T., et al. (2012). Fiji: an open-source platform for biological-image analysis. *Nat. Methods* 9, 676–682. doi: 10.1038/nmeth.2019

Schnepf, E., and Elbrachter, M. (1992). Nutritional strategies in dinoflagellates. *Eur. J. Protist.* 28, 3–24.

Sieracki, K. K., Sieracki, M. E., and Yentsch, C. S. (1998). An imaginig-in-flow system for automated analysis of marine microplankton. *Mar. Ecol. Prog. Ser.* 168, 285–296. doi: 10.3354/meps168285

Simon, N., Barlow, R. G., Marie, D., Partensky, F., and Vaulot, D. (1994). Flow cytometry analysis of oceanic photosynthetic picoeukaryotes. *J. Phycol.* 30, 922–935. doi: 10.1111/j.0022-3646.1994.00922.x

Smyda, T. J. (1997). Harmful algal blooms: their ecophysiology and general relevance to phytoplankton blooms in the sea. *Limnol. Oceanogr.* 42, 1137–1153. doi: 10.4319/lo.1997.42.5_part_2.1137

Sosik, H. M., Olson, R. J., Neubert, M. G., and Shalapyonok, A. (2003). Growth rates of coastal phytoplankton from time-series measurements with a submersible flow cytometer. *Limnol. Oceanogr.* 48, 1756–1765. doi: 10.4319/lo.2003.48.5.1756

Stoecker, D. K. (1999). Mixotrophy among dinoflagellates. *J. Eukaryot. Microbiol.* 46, 397–401. doi: 10.1111/j.1550-7408.1999.tb04619.x

Sweeney, B. M. (1959). "Endogenous diurnal rhythms in marine dinoflagellates," in *International Oceanography Congress Preprints*, ed M. Sears (Washington, DC: America Association for the Advancement of science), 204–207.

Thomas, W. H. (1975). Effects of temperature and illuminance on cell division rates of three species of tropical oceanic phytoplankton. *J. Phycol.* 11, 17–22. doi: 10.1111/j.1529-8817.1975.tb04536.x

Thyssen, M., Beker, B., Ediger, D., Yilmaz, D., Garcia, N., and Denis, M. (2011). Phytoplankton distribution during two contrasted summers in a Mediterranean harbour: combining automated submersible flow cytometry with conventional techniques. *Environ. Monit. Assess.* 173, 1–16. doi: 10.1007/s10661-010-1365-z

Thyssen, M., Garcia, N., and Denis, M. (2009). Sub meso scale phytoplankton distribution in the North East Atlantic surface waters determined with an automated flow cytometer. *Biogeosciences* 6, 569–583. doi: 10.5194/bg-6-569-2009

Thyssen, M., Grégori, G. J., Grisoni, J., Pedrotti, M., Mousseau, L., Artigas, L. F., et al. (2014). Onset of the spring bloom in the northwestern Mediterranean Sea: influence of environmental pulse events on the *in situ* hourly-scale
dynamics of the phytoplankton community structure. *Front. Microbiol.* 5:387.
doi: 10.3389/fmicb.2014.00387

Thyssen, M., Mathieu, D., Garcia, N., and Denis, M. (2008). Short-term variation of phytoplankton assemblages in Mediterranean coastal waters recorded with an automated submerged flow cytometer. *J. Plankton Res.* 30, 1027–1040. doi: 10.1093/plankt/fbn054

Townsend, D. W., Pettygrew, N. R., and Thomas, A. C. (2001). Offshore of the red tide dinoflagellates Alexandrium sp. in the gulf of maine. *Cont. Shelf Res.* 21, 347–369. doi: 10.1016/S0278-4343(00)00093-5

Ulses, C., Grenz, C., Marsaleix, P., Schaaff, E., Estournel, C., Meulé, S., et al. (2005). Circulation in a semi-enclosed bay under influence of strong freshwater input. *J. Mar. Syst.* 56, 113–132. doi: 10.1016/j.jmarsys.2005.02.001

Vaulot, D., and Marie, D. (1999). Diel variability of photosynthetic picoplank-
ton in the equatorial Pacific. *J. Geophys. Res.* 104, 3297–3310. doi: 10.1029/
98JC01333

Von Stosch, H. A. (1973). Observations on vegetative reproduction and sex-
ual life cycles of two freshwater dinoflagellates, *Gymnodinium pseudopalus-
tre Schiller Woloszynskia* apiculata sp. Nov. *Br. Phycol. J.* 8, 105–134 doi:
10.1080/00071617300650141

Walker, L. M. (1982). Evidence for a sexual cycle in the Florida red tide dinoflagellate Ptychodiscus brevis (=Gymnodinium breve). *Bioscience* 32, 809–810. doi: 10.2307/1308977

Weiler, C. S., and Chisholm, S. W. (1976). Phased cell division in natural popula-
tions of marine dinoflagellates from shipboard cultures. *J. Exp. Mar. Biol. Ecol.*
25, 239–247. doi: 10.1016/0022-0981(76)90126-X

Yamamoto, T., OH, S. J., and Kataoka, Y. (2002). Effects of temperature, salinity and irradiance on the growth of the toxic dinoflagellates *Gymnodinium catenatum* (Dynophyceae) isolated from Hiroshima Bay, Japan. *Fish. Sci.* 68, 356–363. doi: 10.1046/j.1444-2906.2002.00433.x

**Conflict of Interest Statement:** The authors declare that the research was con-
ducted in the absence of any commercial or financial relationships that could be
construed as a potential conflict of interest.

Received: 18 April 2014; accepted: 28 August 2014; published online: 15 September 2014.

Citation: Dugenne M, Thyssen M, Nerini D, Mante C, Poggiale J-C, Garcia N, Garcia F and Grégori GJ (2014) Consequence of a sudden wind event on the dynamics of a coastal phytoplankton community: an insight into specific population growth rates using a single cell high frequency approach. *Front. Microbiol.* 5:485. doi: 10.3389/fmicb.2014.00485

This article was submitted to Systems Microbiology, a section of the journal Frontiers in Microbiology.

Copyright © 2014 Dugenne, Thyssen, Nerini, Mante, Poggiale, Garcia, Garcia and Grégori. 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.