山东半岛南部产卵场鳀鱼的产卵生态Ⅰ.鳀鱼鱼卵和仔稚幼鱼的数量与分布特征

2000-2004年5-7月,"北斗"号生物资源调查船对山东半岛南部产卵场进行了11个航次多学科的综合调查,鳀鱼鱼卵和仔稚幼鱼的数量与分布及其与环境因氐墓系为重点调查内容.调查结果表明:1)山东半岛南部产卵场5月中、下旬鳀鱼才进入产卵初期;6月上旬产卵个体逐渐增多,即将进入产卵盛期;6月中、下旬进入产卵盛期; 2)山东半岛南部产卵场中,35°N断面及其邻近海区的鳀鱼鱼卵密集分布区相当稳定,这一带海域基本上每年都有鱼卵的密集分布中心形成;北部海域的密集分布区的分布范围和分布位置变化比较大,但是,密集分布区一旦形成,鱼卵的平均密度相当可观;西南部沿岸一带海域的鱼卵数量很少.鳀鱼鱼卵的分布与海流、潮汐以及水温分布等物理海洋学要素的关系密切,南、北两个鱼卵密集分布区正是潮锋所在的位置,锋区的辐聚和卷夹作用,使鳀鱼鱼卵聚集成密集的斑块,从而形成了鱼卵的密集分布中心; 3)鳀鱼鱼卵在水体中的空间分布主要集中在上混合层,温跃层中的数量次之,下混合层的数量最少;统计学检验, 上混合层、温跃层和下混合层之间的鱼卵数量差异显著; 4)鳀鱼仔稚幼鱼的数量明显少于鱼卵,水平分布特征与鱼卵也有明显的差异.     

黄海中部近岸春夏季鱼卵、仔稚鱼群落结构特征

2013年4月下旬至7月上旬,采用浅水I型网垂直拖网的调查方式对黄海中部20m以浅海域(34°30'—35°30'N,119°20'—120°20'E)的鱼卵、仔稚鱼种类组成与数量分布进行了6个航次的综合调查。春夏季共采集鱼卵21种(未定种3种),共计2247粒,隶属于9科15属。仔稚鱼22种(未定种1种),共计286尾,隶属于13科17属。鱼卵优势种为日本鳀Engraulis japonicus、赤鼻棱鳀Thryssa kammalensis、江口小公鱼Stolephorus commersonnii、皮氏叫姑鱼Johnius belengerii和属Callionymus spp.。仔稚鱼优势种为大银鱼Protosalanx chinensis、白姑鱼Argyrosomus argentatus、日本鳀、赤鼻棱鳀和皮氏叫姑鱼。鱼卵丰度和分布与经度呈负相关、与SST呈正相关,但与其他环境因子如纬度、表层盐度(SSS)和水深等无显著相关性。仔稚鱼丰度和分布与上述环境因子均未表现出显著相关性。基于Bray-Curtis指数的聚类分析表明,鱼卵优势种包括集中分布于研究水域34°48'N以南海域(江口小公鱼、皮氏叫姑鱼和青鳞小沙丁Sardinella zunasi)、10m等深线以内(赤鼻棱鳀),以及北部水域(日本鳀和属)3种主要聚集类型;仔稚鱼优势种包括集中分布于34°48'N—35°12'N海域(赤鼻棱鳀和日本鳀)、连云港外海浅水区(白姑鱼、皮氏叫姑鱼和梭鱼Liza haematocheilus,以及南部水域(褐菖鲉Sebastiscus marmoratus和鲬Platycephalus indicus)3种主要聚集类型。研究结果揭示了黄海中部近岸海域鱼卵、仔稚鱼群落的时空分布格局动态,为科学认识环境变动下该海区渔业种群数量动态和资源变化提供了科学依据。     

长江口水域夏季鱼卵和仔稚鱼年间变化

基于2005年、2008年、2009年和2011年8月(夏季)在长江口水域(30°30'—31°45'N,121°15'—123°10'E)4个航次的浮游生物拖网资料,分析了长江口水域鱼卵和仔稚鱼的种类组成、数量分布特征及其年间变化。结果表明:4个航次采集的鱼卵和仔稚鱼鉴定到种的种类有17种,隶属于8目13科,以鲈形目种类最多,11种,其次是鲱形目,5种,其他各目种类均小于5种;种类数存在明显年间差异,2005年种类数最多(鱼卵3种,仔稚鱼8种),其次是2009年和2011年,2008年种类数最少(鱼卵1种,仔稚鱼5种)。优势种年间更替明显,长蛇鲻(Saurida elongata)、虾虎鱼(Gobiidae spp.)和中华小公鱼(Stolephorus chinensis)在2005年是优势种,2008年优势种为鳀鱼(Engraulis japonicus),2009年优势种为鳀鱼、寡鳞飘鱼(Pseudolaubuca engraulis)、虾虎鱼等,2011年虾虎鱼和小公鱼(包括小公鱼属未定种Stolephorus spp.和中华小公鱼Stolephorus chinensis)成为优势种。2005年鱼卵和仔稚鱼数量分布的密集区在嵊泗列岛附近水域,2008年鱼卵和仔稚鱼出现较少,未出现明显的数量密集区;2009年鱼卵数量较少,仔稚鱼数量较多,密集区主要分布在在长江口以外123°E附近水域;2011年鱼卵主要分布在在长江北支口门外附近水域,仔稚鱼在调查区内分布相对均匀。     

西太平洋热带水域长体翻车鲀鱼卵和仔稚鱼的空间分布及其形态特征

1978年10月利用向阳红5号调查船对西太平洋14°00′-17°00′S、171°00′-174°00′E热带水域进行鱼类浮游生物调查。调查期间共采集到长体翻车鱼卵27粒和仔稚鱼11尾。长体翻车鱼卵集中分布于调查海域东南部的15°28.4′-16°17.3′S、172°43.8′-173°29.0′E海区,而仔稚鱼主要分布于调查海域东南部的15°28.4′-16°17.3′S、172°43.8′-173°29.0′E和西北部的14°29.0′S、171°48.4′E海域。鱼卵和仔稚鱼的密度分别为1-14粒/100m3和1-3尾/100m3。出现鱼卵和仔稚鱼分布海域的水温和盐度分别为25.72-26.22℃、35.10-35.32和25.72-27.02℃、34.80-35.50。鱼卵和仔稚鱼均出现在海洋锋面处。长体翻车鱼卵为圆球型、浮性、卵膜光滑、透明,卵子直径1.40-1.60mm(1.50±0.054mm,n=27),多油球,20-30个,其直径为0.06-0.16mm(0.10±0.064mm,n=27),卵黄均匀无龟裂。当胚胎发育到尾芽期,可数肌节20对,胚体背部和油球表面出现黑色素细胞分布,皮质体壳已初步形成。孵化前,皮质体壳的表面布满了浓密、细微的小泡,唯尾部和胸鳍外,整个胚体被体壳包住。仔鱼的肛前长为1.20mm(全长1.80mm,体长1.68mm)时,体壳已退化脱尽,肛前长为3.20mm(体长4.00mm)时,鱼体上的部分棘开始退化、变小,尾鳍完全退化、消失,尾锥体萎缩成瘤状,腹部的棘与棘之间出现许多锯齿状小刺。卵子的形状和大小、卵膜的结构、油球的数量、黑色素的分布以及体棘的形状和数量等是区别于其它形目鱼类的卵子和仔稚鱼的依据。     

黄海中南部越冬鳀鱼空间分布及其与水温年际变化的关系

利用1986-2010年的底拖网调查数据和遥感海表温度,基于GIS技术,通过空间分析和相关性分析等方法,研究黄海中南部水域越冬鳀鱼空间分布的年际变化,并探讨其与水温变化的关系.结果表明:1986-2010年鳀鱼空间分布的年际变化明显,2004年其资源丰度下降至最低且向岸分布,2010年资源丰度有所回升并在东部水域集中分布,鳀鱼捕获站位、资源密度在经纬度方向上的分布集中区域年际变化明显;鳀鱼资源密度重心纬度与代表性等温线平均纬度的变化趋势相一致,表明水温变化影响越冬鳀鱼纬度方向上的分布,但其向岸分布趋势则由黄海暖流的年际变动决定;在影响鳀鱼分布的众多因子中,资源丰度变化主要受捕捞压力的影响,而空间上的变化主要由海水温度决定.     

黄海鳀鱼时空分布季节差异分析

利用1986-2010年的底拖网调查数据,基于GIS的空间分析和地统计分析技术,从不同角度阐明在捕捞压力和气侯变化的双重胁迫下,黄海鳀鱼资源密度时空分布及其变动的季节差异.结果表明:不同季节鳀鱼空间分布趋势变化明显,冬季,资源密度从调查区域东北部向岸边和南部水域递减,春季,与冬季相反,资源密度由岸边向东部水域递减,夏季,整体上资源密度呈北高南低的趋势,但高值点散乱分布,秋季,与冬季相似,资源密度由东部水域向岸边环形递减,但变化梯度较冬季缓和;鳀鱼资源密度重心季节变化明显,但不同年代变化规律并不相同;1986-2010年,春、夏和秋3个季节,鳀鱼资源密度重心均向北移动,冬季,资源密度重心年间变化明显,但无整体迁移规律;不同季节,鳀鱼资源密度高值聚集区空间分布位置不同.     

中街山列岛海域鱼卵、仔稚鱼群落结构特征及其与环境因子的关系

为了解中街山列岛海域鱼卵仔稚鱼空间分布、种类组成及与环境因子的关系,于2010年5月—2012年2月在中街山列岛海域利用水平和垂直拖网开展了8个航次的鱼卵、仔稚鱼调查。运用非度量多维标度排序、聚类分析和典范对应分析(CCA)等方法对该海域鱼卵、仔稚鱼的空间分布进行了综合分析。结果表明:两年水平拖网共采集到鱼卵1783粒,仔稚鱼67尾;垂直拖网共采集到鱼卵576粒,仔稚鱼59尾。2010年平拖优势种为短吻红舌鳎(Cynoglossus joyneri),垂拖重要种为鲉形目未定种(Scorpaeniformes sp.)和扁舵鰹(Auxis thazard)。2011年平拖重要种为短吻红舌鳎、鰏属未定种(Leiognathuss sp.)、龙头鱼(Harpodon nehereus)和鳀(Engraulis japonicus);垂拖重要种为短吻红舌鳎、龙头鱼和中华小公鱼(Stolephorus chinensis)。经单因子相似性分析(ANOSIM),鱼卵仔稚鱼的优势种和重要种群落物种组成存在显著差异(P0.01)。研究结果表明,温度和盐度对鱼卵仔稚鱼的影响较大,夏季是中街山列岛海域鱼类产卵的重要时期。该海域以恋礁的小型底层鱼类为主,是褐菖鲉(Sebastiscus marmoratus)、龙头鱼等小型底层经济鱼类的重要产卵场和育幼场。     

热带大西洋黄鳍金枪鱼垂直分布空间分析

为了解热带大西洋黄鳍金枪鱼(Thunnus albacares)延绳钓适宜渔获水温的等温线时空分布,分析黄鳍适宜的垂直和水平空间分布范围,采用Argo浮标剖面温度数据重构热带大西洋13℃和距海洋表层水温8℃(Δ8℃)的月平均等温线场,网格化计算了13℃和Δ8℃等温线深度值和温跃层下界深度差,并结合大西洋金枪鱼会委员(International Commission for the Conservation of Atlantic Tunas ICCAT)的黄鳍金枪鱼延绳钓渔业数据,绘制了13℃和Δ8℃等温线深度与月平均CPUE的空间叠加图,用于分析热带大西洋黄鳍金枪鱼中心渔场单位捕捞努力渔获量(Catch per unit effort CPUE)时空分布和次表层环境季节性变化关系。结果表明,13℃等温线,在高值CPUE出现的海域深度值大多小于250 m,主要在170—249 m,深度值超过250 m的海域CPUE普遍较小。5°S—9°N区域,Δ8℃等温线高值CPUE出现的海域深度值大多小于150 m,主要在50—139 m;7—10月份在南半球的非洲西海岸,在Δ8℃等温线深度值为150—350 m的海域也会出现中心渔场。全年在低纬度区域,高渔获率的垂直分布深度更加集中。13℃等温线影响热带大西洋黄鳍金枪鱼的空间分布,温跃层下界温度影响黄鳍金枪鱼的垂直分布。采用频次分析和经验累积分布函数计算其适宜次表层环境因子分布,13℃等温线180—240 m;Δ8℃等温线50—139 m;与下界深度差:13℃等温线-70—29 m;海表以下8℃等温线30—149 m。文章初步得出热带大西洋黄鳍金枪鱼适宜的水平、垂直深度分布区间。结果可以辅助渔情预报,为热带大西洋黄鳍金枪鱼实际生产作业和资源管理提供参考依据。     

黄海近岸海域短吻红舌鳎夏季产卵场的空间分布及其年际变化

根据2014年夏季、2016年夏季与2017年夏季3个航次的黄海近岸海域产卵场调查数据,采用服从Tweedie分布的广义加性模型(GAM),对黄海近岸海域短吻红舌鳎产卵场的空间分布及其年际变化进行了研究。结果表明: 短吻红舌鳎的产卵场主要分布在海州湾(34°00′—35°18′ N,119°30′—121°30′ E)、苏北浅滩(32°18′—34°00′ N,120°18′—122°00′ E)海域范围内,而山东半岛南部海域(34°42′—36°48′ N,119°30′—122°00′ E)鱼卵数量较少。短吻红舌鳎产卵场的空间分布与经纬度、水深及海水表层温度均显著相关,与海水表层盐度、表层叶绿素a相关性不显著;其最适水深、海水表层温度、海水表层盐度、表层叶绿素a浓度范围分别为15～26 m、29～32 ℃、22～25和0.10～3 mg·m^-3。短吻红舌鳎产卵场的空间分布及其影响因子存在一定的年际变化,总体上其产卵场分布较为稳定,但表层温度较高的年份可能存在北移的趋势。     

基于水声学方法的天目湖鱼类季节和昼夜空间分布研究

采用水声学方法系统地对天目湖春夏秋冬四季鱼类水平和垂直空间分布进行了调查研究,同时对昼夜探测的差异性进行了比较,并对不同季节间鱼类聚群形态进行了探讨。研究结果表明,天目湖鱼类在不同季节间存在规律性水平迁移,从全湖角度分析,天目湖鱼类资源昼夜间的水平空间分布无明显差异特征,但在春末和夏季的局部区域里,昼夜间鱼类存在一定近岸—远岸的水平迁移;不同季节和昼夜间,鱼类垂直分布差异明显,且在存在温跃层的夏季7月昼间,鱼类密度垂直分布与水温的垂直分布关系密切,温跃层以下的鱼类密度基本为0;天目湖鱼类在四个季节都属于成群分布类型,但季节间鱼类聚群形态不同,在冬季的1月呈现出典型的聚集,相对于昼间,春夏季鱼类在夜间分布更为均匀、分散,且水平探测表明在夜间水体表层存在大量鱼类分布。     

水体热分层对万峰湖水环境的影响

2009年9月(秋季)和2010年1月(冬季)对南盘江流域峡谷型水库--万峰湖的水温和水化学(DO、pH、总磷)进行监测,结果表明:万峰湖水体在秋、冬季节均存在明显分层,秋季.水体分为3层,0～10 m为混匀层,10～50 m为斜温层,50 m以下为滞水层,这种温度分层阻止了水的对流混合,引起显著的水化学分层,形成底部厌氧层.冬季,水温下降,在水下50 m左右分层,表层随深度增加水温下降,下层为滞水层,水温较为均匀.相关性分析表明,万峰湖在2009年9月(秋季)和2010年1月(冬季)水温和DO、pH、总磷之间均为极显著相关,水化学分层与温度分层同步.     

Vertical distribution of eggs and larvae of the Baltic Sprat <Emphasis Type="Italic">Sprattus sprattus balticus</Emphasis> (Clupeidae) in relation to seasonal and diurnal variation

Vertical distribution of eggs and larvae of the Baltic sprat Sprattus sprattus balticus was studied. Diurnal variation was revealed in the distribution of eggs of the first stage of development and larvae with a length of more than 5 mm, which is related to dynamics of spawning and vertical migrations of sprat larvae. Seasonal variation in the distribution of eggs manifested itself in: (1) widening of the range of occurrence from the zone of halocline at the beginning of spawning (winter) to the entire water column in spring-summer and (2) emergence in the surface layer in July of all stages of egg development unlike the presence from May to June mainly of the first stage of development. It is suggested that the ability of early ontogenetic stages of sprat to maintain in the surface layer depended on the extent of development of vertical gradients of density in the thermocline.     

Catastrophic transition in vertical distributions of phytoplankton: alternative equilibria in a water column

We showed a catastrophic transition between a surface maximum pattern and a subsurface maximum pattern of phytoplankton in a water column by a mathematical model considering the vertical distribution of phytoplankton and two resources, light and a limiting nutrient. In our model, we assumed that a water column consists of two layers: a complete mixing layer above a seasonal thermocline and an incomplete mixing layer below it. From numerical calculation of the model, we obtained that there are two stable vertical patterns of phytoplankton over a certain range of parameters of the model: a pattern having its maximum below the thermocline and another having its maximum above the thermocline. As other models having multiple stable equilibria, our model also exhibits a hysteresis effect and catastrophic transition when one of the parameters of the model changes continuously. These results indicate the possibility of the existence of alternative equilibria of vertical patterns of phytoplankton even if the trophic status and physical condition of the water column are similar. Moreover, the catastrophic transition between the steady states suggests one of the possible mechanisms of autumn algal blooms.     

The non-linear response of plankton to wind mixing events -- implications for survival of larval northern anchovy

A one-dimensionsal (z, t) plankton model incorporating mixedlayer dynamics suggests die following: wind events can be detrimentalto larval northern anchovy, Engraulis mordax, which need tofeed on concentrations of plankton localized within the watercolumn. Wind mixing dissipates vertical structure in prey distributions.Interacting biological and physical processes determine thetime interval before high concentrations of prey are re-established,i.e. the starvation period endured by the anchovy larvae. Phytoplanktongrowth, herbivore grazing and reproduction, and plankton verticalmigration govern the rate of re-establishment of vertical structurein plankton distributions, once wind conditions allow turbulencein the upper water column to dissipate. If there is a net nutrientflux into the euphotic zone, local concentrations of prey afterthe wind event may be higher than before the event. This mayenhance larval anchovy survival, providing the starvation periodis not too long.     

Biogeochemical response to physical forcing in the water column of a warm monomictic lake

Based on microstructure measurements of temperature and horizontalcurrent velocity the physical structure in the water column of Lake Kinneretwascharacterized as a five layer system consisting of a surface mixed layer, lowerepilimnion, metalimnion, upper hypolimnion and benthic boundary layer. Usingoxygen and hydrogen sulfide as natural chemical tracers, the time scale ofchemical change was identified in relation to advection, mixing and biologicalprocesses. Rapid changes due to advection that took place on an hourly timescale were removed by referring the data back to the temperature of the water.Biological activity dominated the hydrochemical changes observed in the meta-and upper hypolimnion. These were expressed by DO depletion rates of 2.0 and0.4g m^–2 d^–1, respectively.Verticaland horizontal mixing were shown to occur on a seasonal time scale. Once thechemical stratification process was completed the slow mixing through thebenthic boundary layer became the limiting factor for subsequent reactions inthe water column.     

Effects of thermal stratification and mixing on reservoir water quality

In this study, the effect of thermal stratification on water quality in a reservoir has been investigated by field observations and statistical analysis. During the summer period, when stratification is evident, field observations indicate that the observed dissolved oxygen concentrations drop well below the standard limit of 5 mg l⁻¹ at the thermocline, leading to the development of anoxia. The reasons for variations in the dissolved oxygen concentrations were investigated. Variations of air temperature and other meteorological factors and lateral flows from side arms of the lake were found to be responsible for the increase of dissolved oxygen concentrations. It was also observed that turbidity peaked mostly in the thermocline region, closely related to the location of the maximum density gradient and thus low turbulence stabilizing the sediments in the vertical water column. Relatively cold sediment-laden water flowing into the lake after rain events also resulted in increased turbidity at the bottom of the lake. Nondimensional analysis widely used in the literature was used to identify the strength of the stratification, but this analysis alone was found insufficient to describe the evolution of dissolved oxygen and turbidity in the water column. Thus correlation of these parameters was investigated by multivariate analysis. Fall (partial mixing), summer (no mixing), and winter (well mixed) models describe the correlation structures between the independent variables (meteorological parameters) and the dependent variables (water-quality parameters). Statistical analysis results indicate that air temperature, one day lagged wind speed, and low humidity affected variation of water-quality parameters.     

Basin-scale internal waves in the bottom boundary layer of ice-covered Lake Müggelsee, Germany

We performed high-resolution temperature measurements under ice cover in Lake Müggelsee, Germany, during the winter of 2005–2006. Intense seiche-like temperature oscillations developing after the ice-on have been encountered in a thin water layer above the sediments. The oscillations were initiated immediately after lake freezing by the release of the potential energy of the thermocline slope and existed for several weeks without appreciable external forcing. The oscillations were associated with a basin-scale internal waves existing in the lower stratified part of the water column. The weakness of the density stratification under ice ensured the long wave periods, exceeding the period of geostrophic inertial oscillations at the lake’s latitude. As a result, two frequency peaks were present in the oscillations corresponding to two rotational waves, one of Kelvin-wave type and another of Poincaré type wave. The rotational character provided long dissipation times of the waves and allowed the oscillations to persist in lake several weeks. Temperature measurements in the upper several centimeters of the sediment demonstrated that oscillations of the near-bottom temperature produced vertical density instability and pore-water convection in the upper sediments.     

Vertical distribution and abundance of mesoplanktonic medusae and siphonophores from the Weddell Sea,Antarctica

The composition, abundance and vertical distribution of mesoplanktonic cnidarians collected along a transect across the Weddell Sea have been analysed. The transect was characterized by a thermocline, approximately between 200 and 100 m, which deepened significantly towards the shelf edges. In total, 10 species of medusae and 18 species of siphonophores were identified. The most abundant medusae were Pantachogon scotti (up to 11,671 specimens/1,000 m³) and Arctapodema ampla (up to 960 specimens/1000 m³). The most abundant siphonophores were Muggiaea bargmannae (up to 1,172 nectophores/1,000 m³) and Dimophyes arctica (up to 230 nectophores/1,000 m³). Five assemblages of planktonic cnidarians were distinguished: (a) epipelagic species located in and above the thermocline; (b) epi- and upper mesopelagic species located in, above and just below the thermocline; (c) epi- and mesopelagic species located in and below the thermocline; (d) mesopelagic species; (e) lower mesopelagic species. Differences in the depth distribution of the various species gave rise to a clear partitioning of the mesoplanktonic cnidarian population throughout the water column. This vertical partitioning was related to the existence of a thermocline, the structure of the water column and the vertical distribution of prey.     

Daily patterns of mixing and nutrient concentrations during early autumn circulation in a small sheltered lake

1. Autumn circulation in lakes is currently conceived to occur very rapidly, being controlled mainly by wind‐power dynamics, decreasing irradiance and heat flux. In addition, autumn mixing is usually related to nutrient redistribution in the vertical column, resulting in its overall increase. To test these assumptions, mixing and nutrient dynamics in a Spanish small, wind‐sheltered, mesotrophic, seepage lake were studied daily during autumn circulation. 2. The seasonal erosion of the pycnocline in Las Madres Lake was the outcome of vertical and horizontal exchanges of heat and matter. The overall mixing of the water column lasted 3 months, which was an unexpected period for a rather shallow lake. Two periods of mixing could be envisaged until full circulation was attained. First, a slightly faster period of pycnocline deepening than that predicted by the heat flux and wind stress model of Fisher et al. (1979) occurred for 41 days, mixing most of the water column down to within two meters of the bottom. Then a much slower process took place promoting frequent instability of the bottom layer and resulting in entire mixing in a further 52 days. 3. Vertically, the whole mixing process was a response to weak surface cooling, resulting from the mild air temperatures of the semiarid climate of the area, and weak wind stress, because of low wind fetch and high shelter. Horizontally, a gravity current transporting cold, denser water from western shallower areas of the lake and materials produced by the decomposition of organic matter of littoral origin may produce a bottom layer of increased density, thus impinging on vertical stability. Seepage inputs of water of roughly constant temperature might also have increased bottom density. Bottom density enhancement resulted in a double diffusion process. 4. Only in‐lake nitrogen content increased until full circulation was attained, whereas carbon showed no trend and phosphorus declined. External processes, such as seepage exchange and atmospheric deposition, coupled to internal processes, such as nitrification, oxidised phosphorus precipitation and complexation with organic carbon, might have been responsible for the areal nutrient patterns observed. 5. Our study demonstrates that current models of water column mixing and nutrient redistribution in lakes during autumn circulation must be improved to encompass the effects of external inputs, including horizontal heat and matter exchange.     

Computer simulations of seasonal outbreak and diurnal vertical migration of cyanobacteria

Algal blooms caused by cyanobacteria are characterized by two features with different time scales: one is seasonal outbreak and collapse of a bloom and the other is diurnal vertical migration. Our two-component mathematical model can simulate both phenomena, in which the state variables are nutrients and cyanobacteria. The model is a set of one-dimensional reaction-advection-diffusion equations, and temporal changes of these two variables are regulated by the following five factors: (1) annual variation of light intensity, (2) diurnal variation of light intensity, (3) annual variation of water temperature, (4) thermal stratification within a water column and (5) the buoyancy regulation mechanism. The seasonal change of cyanobacteria biomass is mainly controlled by factors, (1), (3) and (4), among which annual variations of light intensity and water temperature directly affect the maximum growth rate of cyanobacteria. The latter also contributes to formation of the thermocline during the summer season. Thermal stratification causes a reduction in vertical diffusion and largely prevents mixing of both nutrients and cyanobacteria between the epilimnion and the hypolimnion. Meanwhile, the other two factors, (2) and (5), play a significant role in diurnal vertical migration of cyanobacteria. A key mechanism of vertical migration is buoyancy regulation due to gas-vesicle synthesis and ballast formation, by which a quick reversal between floating and sinking becomes possible within a water column. The mechanism of bloom formation controlled by these five factors is integrated into the one-dimensional model consisting of two reaction-advection-diffusion equations.