ZHAO, ZHU, CHEN, ZHOU, WANG, and LI: Numerical Simulation of the Vertical Migration of Microcystis (Cyanobacteria) Colonies Based on Turbulence Drag

Numerical Simulation of the Vertical Migration of Microcystis (Cyanobacteria) Colonies Based on Turbulence Drag

Abstract

The vertical migration and accumulation of Microcystis is an important process in water blooms, and colony migration is influenced by colony size and wind-wave disturbance. The vertical migration of Microcystis colonies in turbulence can be simulated in a numerical model. In this study, we model such migration by coupling the colony size and hydrodynamics, including the gravity, colony buoyancy, and the viscous drag force of turbulence. The turbulence intensity was represented by the turbulent kinetic energy (KZ); the larger the KZ, the stronger the wind-wave disturbance. The simulated vertical distribution of Microcystis well agreed with the measured values in a laboratory experiment indicating that our model can simulate the vertical distribution of Microcystis under different hydrodynamic conditions. We also found a size-dependent critical turbulent kinetic energy (TKZ), such that if the turbulent kinetic energy of water exceeds the critical value (i.e., KZ > TKZ), the colonies sink under the drag forces of turbulence; conversely, if KZ < TKZ, the colonies can overcome the turbulent mixing and float. The TKZ of each colony was linearly related to colony diameter. The model is crucial for prediction and prevention of water blooms. The simulated threshold turbulent kinetic energy, at which water blooms disappear in Lake Taihu (a large freshwater lake in the Yangtze Delta, Jiangsu Province, China), was 55.5 cm2s–2.




INTRODUCTION

Cyanobacterial blooms are a major environmental problem in eutrophic lakes and Microcystis is the most common genus of bloom forming cyanobacterium. Water blooms can be effectively predicted in numerical models that simulate the occurrence and disappearance of water blooms. The cyclic appearance and disappearance of water blooms is mainly driven by the spatial aggregation and dispersion of Microcystis by wind waves (Wu and Kong, 2009; Zhu et al., 2014). Consequently, this process is assumed in numerical simulations of blooms.

Microcystis colonies achieved vertical migration with the help of buoyancy regulating system of gas vesicles and ballast (e.g., carbohydrate) (Oliver and Walsby, 1984; Kromkamp et al., 1984). These mechanisms were well described by Walsby (1994) that Microcystis colonies always rise during night due to synthesis of gas vesicles and carbohydrate utilization. In the contrary, they sink in the day because of the accumulation of carbohydrate and collapse of gas vesicles. In a laboratory study, Nakamura et al. (1993) found that the measured floating velocity of Microcystis is proportional to colony size. Wu and Kong (2009) investigated Microcystis blooms in Lake Taihu, a large freshwater lake in the Yangtze Delta of Jiangsu Province, China. They found that large Microcystis colonies (>120 µm) can overcome wind-wave disturbance and thus accumulate at the water surface, whereas small colonies (<36 µm) tend to mix with the turbulence. Chien et al. (2013) simulated the vertical migration of differently sized colonies in a trajectory model, and reported that large colonies (≥300 µm) can resist turbulent mixing and undertake diurnal vertical migration, whereas smaller colonies are more homogenously scattered in the vertical direction. According to Stokes’ law (V=2gr2(ρ–)A/(9Φn), where V and r are the velocity and radius of the colony, respectively, ρ and ρw are the densities of the colony and water, respectively, A is the proportion of cell volume relative to colony volume, Φ the form resistance and n the viscosity of the water), without any other interference, the vertical migration velocity of Microcystis colonies mainly depends on the cell density and the colony size (Reynolds et al., 1987), and the floating rate is proportional to the square of the colony diameter. All of these studies suggest that colony size directly affects the vertical migration of Microcystis.

Besides the internal factors, the vertical distribution of Microcystis is mainly driven by wind-wave disturbance. George and Edwards (1976) investigated the relationship between water blooms and wind speed in Lake Nynyod (South Wales), and found that algae gather into blooms at the water surface at wind speeds below 3.0 ms–1, but disperse at wind speeds of 3.0 ms–1. Cao et al. (2006) conducted a field investigation of Lake Taihu, and observed that at a wind speed of 2.0 m s–1 and wave height of 4.4 cm, approximately 37% of the whole Microcystis biomass gathers on the water surface, forming water blooms. Conversely, at a wind speed of 3.1 ms–1 and wave height of 6.2 cm, the Microcystis tend to distribute homogenously through the water column. Wu et al. (2013) found that disturbances mix the Microcystis in deeper layers. In a field study, Ding et al. (2012) investigated the effect of the tropical storm Morakot on Lake Taihu’s water blooms, and observed that chlorophyll a concentration homogenously distributed during the typhoon crossing, and water blooms reappeared on the surface after the typhoon. The above studies show that turbulent mixing determines the vertical distribution of Microcystis. However, although these studies relate wind speed to water blooms, they do not clarify the motion mechanism of Microcystis under different hydrodynamic conditions, nor the effect of turbulence on the migration of colonies.

Many researchers have established mathematical models to study the vertical migration and distribution of Microcystis colonies under different hydrodynamic conditions (Visser et al., 1997; Wallace et al., 2000; Huisman and Sommeijer, 2002; Rabouille et al., 2005; Serizawa et al., 2008; Chen et al., 2009). Visser et al. (1997) combined a model of cell density change with Stokes’ law, and simulated the vertical migration of Microcystis colonies in static water. Wallace et al. (2000) combined a buoyancy regulation model with a hydrodynamic model into a vertical migration model of Microcystis, and computed the vertical position of Microcystis in a turbulent shallow lake (depth 1.4 m). A similar vertical migration model of Microcystis colonies was developed by Chen et al. (2009). The fluid mechanics described particle migration in the flow in details and indicated that the viscous drag force was the driven factor making the particles to migrate in the flow (White et al., 1977; Van Rijn, 1984). This principle could also be applied to Microcystis colonies. The viscous drag was ignored in the models of the vertical migration of Microcystis colony. To more accurately simulate the vertical migration of Microcystis colonies, we require numerical models based on the most influential physical processes.

In this study, we investigate the effect of Microcystis colony size on resistance to disturbance. To this end, we develop a mathematical model that allows floating and sinking of Microcystis colonies under viscous drag as well as gravity and buoyancy. The mathematical model is verified in a laboratory experiment.

METHODS

Forces acting on Microcystis colonies

In the present model, the forces acting on Microcystis colonies are restricted to the vertical direction. The three main forces acting on Microcystis in water (where the z axis is positive downward) are listed below.

1) Gravity

(eq. 1)
limno-2017-1-1501-e001.jpg

2) Buoyancy

(eq. 2)
limno-2017-1-1501-e002.jpg

3) Viscous drag force of turbulence

(eq. 3)
limno-2017-1-1501-e003.jpg

These forces are summed to give the resultant force

(eq. 4)
limno-2017-1-1501-e004.jpg
ρw and ρs as the densities of water (998 kg m–3) and Microcystis, respectively. Reynolds et al. (1981) indicated that the lowest density of Microcystis colonies in Bleham Tarn, Cumbria was 985 kg m–3. In addition, we recently measured the density of Microcystis colonies taken from Lake Taihu and found that the average value was 985 kg m–3. Therefore, we set ρs as 985 kg m3 in the current study. We also require the gravitational acceleration g (9.8 m s–2), the volumes V of a single colony and Vf of the colony in water, the colony diameter D, and the drag coefficient CD. As the colony is assumed to be spherical, its volume is given by V = πD3/6; further, we assume Vf = V, and set CD to 0.45 in turbulence (Dong et al., 2007). We also define uw and vs as the vertical velocities of flow and vertical migration of Microcystis colonies respectively, and a as the acceleration of the colonies. The mass m of a single colony is calculated as m=ρsV.

Vertical position of Microcystis colonies

As turbulence varies in a transient manner, we divided the motion time t into n equal segments with a time step t of 10–4 s. After discretizing eq. (4) as eq. (5), we solved the vertical migration velocity of each Microcystis colony by a finite difference scheme implemented in MATLAB R2009a software (The MathWorks).

(eq. 5)
limno-2017-1-1501-e005.jpg

The initial vertical migration velocity of the colonies (regardless of their diameter) was set to 0 m s–1 (vs0 = 0). vsi was determined from the known vertical velocities of flow (uwi–1). Within the very short time step, the motion of the colonies is assumed uniform and their instantaneous displacement is given by Si=vsi × t. The total displacement (S) and position (z) of a Microcystis colony in water are described by the following flow equations:

(eq. 6)
limno-2017-1-1501-e006.jpg
(eq. 7)
limno-2017-1-1501-e007.jpg

Where z0 is the initial position of the Microcystis colony.

Turbulent kinetic energy

In this study, the disturbance intensity was represented by the turbulent kinetic energy (KZ), which is calculated by eq. (8):

(eq. 8)
limno-2017-1-1501-e008.jpg

Here, ui are the instantaneous vertical flow velocities and ū is their mean value, n is the number of ui. The larger the KZ, the stronger the disturbance.

Laboratory experiment (vertical wave-making simulator) and simulation

The effects of disturbance on the vertical migration and distribution of Microcystis colonies were investigated in a vertical wave-making simulator. Survival kinds of amplitude A (4, 6, 8 cm) with variable frequency (1/T) were set to obtain different hydrodynamic conditions. Microcystis samples taken from Lake Taihu, which were mainly composed of M. aeruginosa, were used in this experiment (Xiao et al., 2013).

Measurement of vertical flow velocity

The vertical flow velocity (u) in the simulator was measured by Particle Image Velocimetry (PIV), and the KZ was calculated by Eq. (8). Fig. 1a shows the working principle of PIV. Before measurement, the tracing particles are placed in the simulator with high flow, and one side of the simulator is illuminated by a pair of lasers. The trajectories of the particles are then tracked perpendicularly to the laser beam by a high-speed CCD camera. The particle trajectories are photographed over 100 μs (the time step of the lasers), and the two-dimensional velocity of the water particle is calculated. Because the laser beam must pass through regions containing liquid, the lowest water level must lie outside the light irradiation area. To ensure this condition, we placed the laser at 50 cm and 80 cm from the bottom of the simulator (Fig. 1b). The vertical flow velocities were calculated from 80 photographs under different hydrodynamic conditions.

Simulation of vertical flow velocity

Considering the wave patterns in the simulator and random turbulence, the vertical flow velocity in the simulator was estimated by Eq. (9), where A, T are amplitude and period set in experiment, respectively, w=2π/T. A random number R was employed here to indicate random fluctuation component of turbulence which was always described as u = ū +u’ (where u is the instantaneous velocity, ū is the average velocity, and u’ is the random fluctuation component).

(eq. 9)
limno-2017-1-1501-e009.jpg

The vertical flow of the wave-making simulator under different conditions was simulated by eq. (9). For each hydrodynamic condition, the flow velocities were calculated in one period, and the turbulent kinetic energy (KZ) was estimated by eq. (8).

Simulation of vertical distribution of Microcystis

Four groups of Microcystis colonies (un-sieved original mixed-size Microcystis samples, colony sizes of 100-300, 300-500 and >500 µm) were operated in the experiment of vertical distribution of Microcystis under different hydrodynamic conditions (Xiao et al., 2013).

Simulation was conducted in following steps. All Microcystis colonies were initially set at the surface (zo = 0). And colonies’ displacements at different time (1, 5, 10 min) were calculated by applying eq. (1) - eq. (7). Then, colonies at a particular depth (z) were counted (Nz) according to their displacements. Finally, the Chlorophyll a concentration at different depths was determined through multiplying ratio of Nz and N (the total number of the colonies in the system) by the whole Chlorophyll a concentration C0 (acquired from laboratory experiment). Thus, the Microcystis concentration at depth (z) is given by eq. (10):

(eq. 10)
limno-2017-1-1501-e010.jpg

The measured results obtained by Xiao et al. (2013) was used to assess the consistency of the predicted results and the measured results.

Critical hydrodynamic conditions

The condition under which the colonies happen to sink is the critical hydrodynamic disturbance, below this hydrodynamic condition, the Microcystis colonies can resist the water motions and float up. In this study, the critical turbulent kinetic energy (TKZ) was used to represent the disturbance intensity. If the turbulent kinetic energy of the water body exceeds the critical kinetic energy of the colonies (i.e., KZ> TKZ), the colonies will sink under the drag forces of turbulence. Conversely, if the turbulent kinetic energy of water is below the critical kinetic energy of the colonies (KZ < TKZ), the colonies will float.

Individual colony was placed on the water surface separately. The critical turbulent kinetic energy (TKZ) under varying hydrodynamic conditions was obtained when the colony happened to leave the water surface. The experiments were performed in three times.

In the numerical simulation, the colony displacements were investigated under increasing hydrodynamic disturbances until the displacement is just greater than zero. The critical hydrodynamic conditions of colonies with different diameters (100, 200, 300, 400, 500 µm) were obtained, hence their critical turbulent kinetic energies (TKZ). Linear regression was used to estimate the dependence of critical kinetic energy on colony size.

RESULTS

Turbulent kinetic energy under different hydrodynamic conditions

The simulated turbulent kinetic energy (KZ) under different hydrodynamic conditions is showed in Tab. 1 and Fig. 2. The turbulent kinetic energy (KZ) is a quadratic function of the hydrodynamic conditions (A/T) as equation 11. The measured and simulated values favorably agree, confirming that our model can simulate the flows in a vertical wave-making simulator.

(eq. 11)
limno-2017-1-1501-e011.jpg

Vertical distribution of Microcystis

After estimating the displacements of all Microcystis colonies under various experimental conditions, we calculated the chlorophyll a concentrations at different depths by eq. (10), and compared them with the experimental results. Fig. 3 showed the variation in vertical distribution of Microcystis colonies (300-500 µm). At an amplitude of 8 cm and period of 2 s, most of the Microcystis gathered at the surface after 1 min. The surface aggregation decreased after 5 min disturbance, and the Microcystis were homogeneously distributed after 10 min. The simulated values well agreed with the measured values, and exhibited the same behavior. Under less vigorous hydrodynamic conditions (A = 6 cm, T = 2 s; A = 6cm, T = 3 s), the simulated chlorophyll a concentration at the surface was higher than the measured one, and the simulated levels at lower depths agreed with the measured levels.

Critical turbulent kinetic energy

The critical turbulent kinetic energy (TKZ) was proportional to colony diameter (D) as equation 12. The strong agreement between the measured and simulated values is clarified in Fig. 4.

(eq. 12)
limno-2017-1-1501-e012.jpg

DISCUSSION

The 1-D wave simulator designed by Xiao et al. (2013) was used to simulate the water movement vertically only. It simplified the water movement of real lakes which varies temporally and spatially, and the simulation was reasonable. The vertical turbulence induced by the motor was strong at the deeper water column, which was consistent with the investigations in Lake Taihu that, deeper layers had more intensive turbulence. Therefore, the effect of vertical water movement on the vertical distribution of Microcystis sp. could be well studied with the device and the experimental set up designed by Xiao et al. (2013). In the present study, model results have been compared with detailed measurements in three aspects (vertical flow velocity, the vertical distribution of Microcystis, and the critical turbulent kinetic energy of different-sized colonies). Simulated values agreed with measured values.

Water motion affects the vertical migration and distribution of Microcystis in lakes, rivers and reservoirs (Ha et al., 2000; Wallace et al., 2000; Medrano et al., 2013; Cui et al., 2016). Microcystis will migrate to the surface forming water blooms under weaker disturbance, and sink to deeper depths as water mixing become stronger. These observations consistent with our experimental data and model results.

On the other hand, hydrological conditions affect water blooms directly or indirectly by means of affecting the sizes of Microcystis colonies, which can alter the floating capacity of colonies. Large colonies may be destroyed to smaller colonies under oscillation mixing (Robarts and Zohary, 1984; O’Brien et al., 2004), and these smaller colonies can’t overcome the drag force of turbulence and sink. Resulting in decrease of Microcystis biomass at the water surface. These may be a reasonable explanation for the discrepancy between measured and simulated Microcystis concentration at the surface under moderate mixing (A = 6 cm, T = 2 s; A = 6 cm, T = 3 s). Our simulation results showed that larger colonies (>400 µm) can resist the disturbance and remain at the water surface under these two conditions. Some of these large colonies may break up into smaller colonies under sustained oscillation mixing in the experiments. While these smaller colonies would sink for that they can’t resist turbulent mixing. Resulting in the decrease of Microcystis biomass at the water surface. But colonies break-up was ignored in our simulation, leading to the differences between simulated and measured Microcystis concentrations. The relationship between the measured and simulated Microcystis concentrations was showed in Fig. 5. The figure indicated that our model can accurately simulate the vertical distribution of Microcystis under different hydrodynamic conditions.

However, this kind of dispersion is small even for the highest turbulence intensities experienced by Microcystis in the field (O’Brien et al., 2004). Some studies discussed the influence of turbulence on the ecology of Microcystis (Reynolds, 1994; Wilkinson et al., 2016). Their results suggested that M. aeruginosa are powerful enough to tolerate intensive mixing for a few days (>96 h) (Regel et al., 2004). In summary, it is reasonable to ignore the effect of turbulence on colony viability and size in simulation process, and simulated results are suitable for field.

The phenomenon has been found that water blooms occurred under poor weather or nutrient conditions (Ozawa et al., 2005; Ma et al., 2015). The reason may be explained as the effect of water motion. It may not be suitable to predict the outbreaks of water bloom by simply relying on the nutrient, light or temperature. Establishing the relationship between water and Microcystis motion is necessary. While existing numerical models and studies reported consistent results. Namely, that Microcystis large colonies can resist turbulence and float to the surface whereas small colonies are easily mixed in the water column (Wallace et al., 2000, Medrano et al., 2013). But failed to quantify the relationship between disturbance resistance and colony diameter. In this study, our model computed the critical turbulent kinetic energies as eq. 12, above which the colonies will sink.

According to the long-term monitoring study of Meiliang and Gonghu Bays in Lake Taihu, a large shallow lake (Zhu et al., 2015), reported that during the water blooms season (October to December), Microcystis colonies are 300-400 µm wide in Lake Taihu (Zhu et al., 2016). The critical turbulent kinetic energy of 400 µm colony is 55.5 cm2 s–2, it means that under this condition, most of the colonies will sink and the water blooms will disappear. In addition, the turbulent kinetic energy (KZ) is the following quadratic function of the wind speed x in Lake Taihu:

(eq. 13)
limno-2017-1-1501-e013.jpg

Inserting the critical turbulent kinetic energy into this expression, the corresponding wind speed is determined as 3.3 ms–1, similar to the threshold wind speed (3.0 ms–1) at which the water blooms disappears in Lake Taihu (Zhu and Cai, 1997; Cao et al., 2006; Zhang et al., 2008). At a wind speed of 1 m s–1, the turbulent kinetic energy of water is 25.2 cm2 s–2, indicating that Microcystis colonies will sink if their diameters are below 200 µm. When the wind speed increases to 4.02 m s–1, the turbulent kinetic energy of water reaches 71.0 cm2 s–2, and colonies with diameters smaller than 500 µm will leave the water surface. At a wind speed of 5.61 ms–1, the turbulent kinetic energy is 117.7 cm2 s–2, ensuring that even large colonies (800 µm) are mixed throughout the water column. The critical turbulent kinetic energy is a useful parameter for predicting the occurrence and disappearance of water blooms and the depth at which Microcystis potentially occurs in different type of lakes.

Wind-waves frequently disturb phytoplankton and sediment in shallow lakes. Light, temperature and nutrients are homogeneous in the whole water column mainly affected by water disturbance (Nixdorf and Deneke, 1997). And colonies will be evenly mixed in the whole water column under intensive mixing conditions, assemble at the near surface in calm periods. However, in deep lakes, disturbance degree varies with depth under disturbance, thus different-sized colonies will gather at different depths (at which the turbulent kinetic energy of water no more than their critical energy). Microcystis will occur at this lower intensity mixing and deeper locations. Light, temperature, may not be favorable for Microcystis growth at deep layers, which will affect algal physiology of Microcystis. Thus affecting morphological and size of Microcystis colony. Shallow lakes may hold the potential to bear favorable conditions for the presence of Microcystis blooms, because of the homogeneous light, temperature, and turbulence mixing, as well as the frequently disturbance which can cause the release of nutrients from the sediment (Fan et al., 2004).

CONCLUSIONS

This paper introduced a vertical migration model of Microcystis colonies based on turbulence drag, revealing the mechanism by which turbulence and colony size affects the vertical migration of Microcystis. The critical turbulent kinetic energy (TKZ), above which the colony cannot resist the turbulence conditions and sinks, depends on the colony diameter. The critical turbulent kinetic energy is a useful parameter for predicting the occurrence and disappearance of water blooms in different type of lakes. The study is crucial for prediction and prevention of water blooms.

ACKNOWLEDGMENTS

The project was funded by the Program on Furtherance of Scientific Research of Japan, Fundament C (15K00630), the Lake Taihu water pollution control special funds (Tenth Phase) Scientific research topic (JSZC-G2016-198), the Fundamental Research Funds for the Central Universities (2015B36214).

REFERENCES

1 

HS Cao, FX Kong, LC Luo, XL Shi, Z Yang, XF Zhang, Y Tao, 2006. Effects of wind and wind-induced waves on vertical phytoplankton distribution and surface blooms of Microcystis aeruginosa in Lake Taihu. J. Freshwater Ecol. 21:231-238.

2 

Y Chen, X Qian, Y Zhang, 2009. Modelling turbulent dispersion of buoyancy regulating cyanobacteria in wind-driven currents, p. 1-4. In: Proc. 3rd Int. Conf. on IEEE Bioinformatics and Biomedical Engineering, 2009. ICBBE.

3 

YC Chien, SC Wu, WC Chen, CC Chou, 2013. Model simulation of diurnal vertical migration patterns of different-sized colonies of Microcystis employing a particle trajectory approach. Environ. Eng. Sci. 30:179-186.

4 

YJ Cui, DF Liu, J Zhang, ZJ Yang, ST Khu, LH Long, 2016. Diel migration of Microcystis during an algal bloom event in the Three Gorges Reservoir, China. Environ. Earth Sci. 75:1-9.

5 

Y Ding, B Qin, G Zhu, T Wu, Y Wang, L Luo, 2012. Effects of typhoon Morakot on a large shallow lake ecosystem, Lake Taihu, China. Ecohydrology 5:798-807.

6 

CY Dong, WL Luan, ST Zhou, Q Zhang, 2007. [Analysis and application of model for solid particle movement in Newton fluid].[Article in Chinese with English abstract]. J. China Univ. Petrol. (Edition of Natural Science) 31:56-63.

7 

CX Fan, L Zhang, BQ Qin, S Wang, W Hu, C Zhang, 2004. Estimation on dynamic release of phosphorus from wind-induced suspended particulate matter in Lake Taihu. Sci. Chin. Ser. D. 47:710.719.

8 

D George, R Edwards, 1976. The effect of wind on the distribution of chlorophyll a and crustacean plankton in a shallow eutrophic reservoir. J. Appl. Ecol. 13:667-690.

9 

K Ha, HW Kim, KS Jeong, GJ Joo, 2000. Vertical distribution of Microcystis population in the regulated Nakdong River, Korea. Limnology 1:225-230.

10 

J Huisman, B Sommeijer, 2002. Population dynamics of sinking phytoplankton in light-limited environments: simulation techniques and critical parameters. J. Sea. Res. 48:83-96.

11 

JC Kromkamp, LR Mur, 1984. Buoyant density changes in the cyanobacterium Microcystis aeruginosa due to changes in the cellular carbohydrate content. FEMS Microbiol. Lett. 25:105-109.

12 

J Ma, B Qin, HW Paerl, JD Brookes, NS Hall, K Shi, 2015. The persistence of cyanobacterial (Microcystis spp.) blooms throughout winter in Lake Taihu, China. Limnol. Oceanogr. 61:711-722.

13 

EA Medrano, RE Uittenbogaard, LMD Pires, BJH Van De Wiel, HJH Clercx, 2013. Coupling hydrodynamics and buoyancy regulation in Microcystis aeruginosa for its vertical distribution in lakes. Ecol. Model. 248:41-56.

14 

T Nakamura, Y Adachi, M Suzuki, 1993. Flotation and sedimentation of a single Microcystis floc collected from surface bloom. Water Res. 27:979-983.

15 

B Nixdorf, R Deneke, 1997. Why ‘very shallow’ lakes are more successful opposing reduced nutrient loads. Hydrobiologia 342/343:269-284.

16 

KR O’Brien, DL Meyer, AM Waite, GN Ivey, DP Hamilton, 2004. Disaggregation of Microcystis aeruginosa colonies under turbulent mixing: laboratory experiments in a gridstirred tank. Hydrobiologia 519:143-152.

17 

RL Oliver, AE Walsby, 1984. Direct evidence for the role of light-mediated gas vesicle collapse in the buoyancy regulation of Anabaena flos-aquae (cyanobacteria). Limnol. Oceanogr. 29:879-886.

18 

K Ozawa, H Fujioka, M Muranaka, 2005. Spatial distribution and temporal variation of Microcystis species composition and microcystin concentration in Lake Biwa. Environ. Toxicol. 20:270-276.

19 

S Rabouille, M-J Salençon, J-M Thébault, 2005. Functional analysis of Microcystis vertical migration: a dynamic model as a prospecting tool: I—Processes analysis. Ecol. Model. 188:386-403.

20 

RH Regel, JD Brookes, GG Ganf, RW Griffiths, 2004. The influence of experimentally generated turbulence on the Mash01 unicellular Microcystis aeruginosa strain. Hydrobiologia 517:107-120.

21 

CS Reynolds, 1994. The role of fluid motion in the dynamics of phytoplankton in lakes and rivers. p. 141-187. In: P.S. Giller, A.G. Hildrew, D.G. Raffaelli (eds.), Aquatic ecology: scale, pattern and process. Blackwell Scientific Publications.

22 

CS Reynolds, RL Oliver, AE Walsby, 1987. Cyanobacterial dominance: the role of buoyancy regulation in dynamic lake environments. New Zeal. J. Mar. Fresh. 21:379-390.

23 

CS Reynolds, G Jaworski, H Cmiech, G Leedale, 1981. On the annual cycle of the blue-green alga Microcystis aeruginosa Kutz. emend. Elenkin. Philos. T. R. Soc. B. 293:419-477.

24 

RD Robarts, T Zohary, 1984. Microcystis aeruginosa and underwater light attenuation in a hypertrophic lake (Hartbeespoort Dam, South Africa). J. Ecol. 72:1001-1017.

25 

H Serizawa, T Amemiya, AG Rossberg, K Itoh, 2008. Computer simulations of seasonal outbreak and diurnal vertical migration of cyanobacteria. Limnology 9:185-194.

26 

LC Van Rijn, 1984. Sediment transport, part I: bed load transport. J. Hydraul. Eng. 110:1431-1456.

27 

PM Visser, J Passarge, LR Mur, 1997. Modelling vertical migration of the cyanobacterium Microcystis. Hydrobiologia 349:99-109.

28 

BB Wallace, MC Bailey, DP Hamilton, 2000. Simulation of vertical position of buoyancy regulating Microcystis aeruginosa in a shallow eutrophic lake. Aquat. Sci. 62:320-333.

29 

AE Walsby, 1994. Gas vesicles. Microbiol. Rev. 58:94-144.

30 

BR White, JC Schulz, 1977. Magnus effect in saltation. J. Fluid Mech. 81:497-512.

31 

A Wilkinson, M Hondzo, M Guala, 2016. Effect of small-scale turbulence on the growth and metabolism of Microcystis aeruginosa. Adv. Microbiol. 6:351.

32 

T Wu, B Qin, G Zhu, L Luo, Y Ding, G Bian, 2013. Dynamics of cyanobacterial bloom formation during short-term hydrodynamic fluctuation in a large shallow, eutrophic, and windexposed Lake Taihu, China. Environ. Sci. Pollut. Res. 20:8546-8556.

33 

X Wu, F Kong, 2009. Effects of light and wind speed on the vertical distribution of Microcystis aeruginosa colonies of different sizes during a summer bloom. Int. Rev. Hydrobiol. 94:258-266.

34 

M Xiao, 2014. [The influence of vertical turbulence on the vertical dispersion and accumulation of Microcystis colonies].[Article in Chinese with English Abstract]. Nanjing: Hohai University

35 

M Xiao, W Zhu, M Li, QQ Sun, PN Nkrumah, X Tan, 2013. The influence of water oscillation on the vertical distribution of Microcystis colonies of different sizes. Fresen. Environ. Bull. 22:3511-3518.

36 

M Zhang, F Kong, X Wu, P Xing, 2008. Different photochemical responses of phytoplankters from the large shallow Taihu Lake of subtropical China in relation to light and mixing. Hydrobiologia 603:267-278.

37 

W Zhu, M Li, X Dai, M Xiao, 2015. Differences in vertical distribution of Microcystis morphospecies composition in a shallow hypertrophic lake (Lake Taihu, China). Environ. Earth Sci. 73(9): 5721-5730.

38 

W Zhu, M Li, Y Luo, X Dai, L Guo, M Xiao, J Huang, X Tan, 2014. Vertical distribution of Microcystis colony size in Lake Taihu: Its role in algal blooms. J. Great Lakes Res. 40:949-955.

39 

W Zhu, X Zhou, H Chen, L Gao, M Xiao, M Li, 2016. High nutrient concentration and temperature alleviated formation of large colonies of Microcystis: evidence from field investigations and laboratory experiments. Water Res. 101:167-175.

40 

Y Zhu, Q Cai, 1997. [The dynamic research of the influence on wind field on the migration of algae in Taihu Lake].[Article in Chinese with English Abstract]. J. Lake Sci. 9:152-158.

Fig. 1.

a) The working principle of PIV. b) Area of flow measurement.

limno-2017-1-1501-g001.jpg
Fig. 2.

Turbulent kinetic energy under different hydrodynamic conditions.

limno-2017-1-1501-g002.jpg
Fig. 3.

The vertical distribution of Microcystis under different hydrodynamic conditions. The continuous lines are simulated results; the dots are measured results from Xiao et al. (2013).

limno-2017-1-1501-g003.jpg
Fig. 4.

The critical turbulent kinetic energy of different-sized Microcystis colonies. Error bars represent the standard deviation of measured values.

limno-2017-1-1501-g004.jpg
Fig. 5.

The relationship between the measured and simulated Microcystis concentrations.

limno-2017-1-1501-g005.jpg
Tab. 1.

The simulated turbulent kinetic energy in different hydrodynamic conditions.

Amplitude A/period T (cm s–1) 8/2 6/2 8/3 6/3 6/4 5/4 4/4
Hydrodynamic condition A/T (cm s–1) 4 3 2.67 2 1.5 1.25 1
Turbulent kinetic energy KZ (cm2 s–2) 160 90 71 40 23 16 10
Abstract views:
427

Views:
PDF
142
HTML
210

Article Metrics

Metrics Loading ...

Metrics powered by PLOS ALM


Copyright (c) 2016 Hongru Zhao, Wei Zhu, Huaimin Chen, Xiaohua Zhou, Ruochen Wang, Ming Li

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.
 
© PAGEPress 2008-2017     -     PAGEPress is a registered trademark property of PAGEPress srl, Italy.     -     VAT: IT02125780185