Temporal variations of food web in a marine bay ecosystem based on LIM-MCMC model

Climate change has led to significant fluctuations in marine ecosystems,including alterations in the structure and function of food webs and ecosystem status.Coastal ecosystems are critical to the functioning of the earth’s lifesupporting systems.However,temporal variations in most of these ecosystems have remained unclear so far.In this study,we employed a linear inverse model with Markov Chain Monte Carlo(LIM-MCMC)combined with ecological network analysis to reveal the temporal variations of the food web in Haizhou Bay of China.Food webs were constructed based on diet composition data in this ecosystem during the year of 2011 and 2018.Results indicated that there were obvious temporal variations in the composition of food webs in autumn of 2011 and 2018.The number of prey and predators for most species in food web decreased in 2018 compared with 2011,especially for Trichiurus lepturus,zooplankton,Amblychaeturichthys hexanema,and Loligo sp.Ecological network analysis showed that the complexity of food web structure could be reflected by comprehensive analysis of compartmentalized indicators.Haizhou Bay ecosystem was more mature and stable in 2011,while the ecosystem’s self-sustainability and recovery from disturbances were accelerated from 2011 to 2018.These findings contribute to our understanding of the dynamics of marine ecosystems and highlight the importance of comprehensive analysis of marine food webs.This work provides a framework for assessing and comparing temporal variations in marine ecosystems,which provides essential information and scientific guidance for the Ecosystem-based Fisheries Management.

Optimization of environmental variables in habitat suitability modeling for mantis shrimp Oratosquilla oratoria in the Haizhou Bay and adjacent waters

Habitat suitability index(HSI)models have been widely used to analyze the relationship between species abundance and environmental factors,and ultimately inform management of marine species.The response of species abundance to each environmental variable is different and habitat requirements may change over life history stages and seasons.Therefore,it is necessary to determine the optimal combination of environmental variables in HSI modelling.In this study,generalized additive models(GAMs)were used to determine which environmental variables to be included in the HSI models.Significant variables were retained and weighted in the HSI model according to their relative contribution(%)to the total deviation explained by the boosted regression tree(BRT).The HSI models were applied to evaluate the habitat suitability of mantis shrimp Oratosquilla oratoria in the Haizhou Bay and adjacent areas in 2011 and 2013–2017.Ontogenetic and seasonal variations in HSI models of mantis shrimp were also examined.Among the four models(non-optimized model,BRT informed HSI model,GAM informed HSI model,and both BRT and GAM informed HSI model),both BRT and GAM informed HSI model showed the best performance.Four environmental variables(bottom temperature,depth,distance offshore and sediment type)were selected in the HSI models for four groups(spring-juvenile,spring-adult,falljuvenile and fall-adult)of mantis shrimp.The distribution of habitat suitability showed similar patterns between juveniles and adults,but obvious seasonal variations were observed.This study suggests that the process of optimizing environmental variables in HSI models improves the performance of HSI models,and this optimization strategy could be extended to other marine organisms to enhance the understanding of the habitat suitability of target species.

Fishery biology of whitespotted conger Conger myriaster(Brevoort, 1856) in the Yellow Sea and East China Sea

Whitespotted conger Conger myriaster is a commercially important species in the seas around China, Korea and Japan. The coastal waters of China serve as an important feeding ground for congers, but the spatio-temporal variations in the fishery and biological characteristics of the population have been rarely evaluated and less well understood in this area. We studied the growth, spawning and feeding characteristics of C. myriaster on the basis of samples collected from October 2016 to April 2017 in the Yellow Sea and East China Sea. A total of 529 specimens were collected, with ages ranging from 1 to 6 years and total length ranging from 132 mm to 834 mm. The parameters of von Bertalanffy growth equation L∞ and k were 1 026 mm and 0.226 a^–1, respectively;the sex ratio was 88:0 (female: male) in the East China Sea and 2.67:1 in the South Yellow Sea;the development stage of ovary ranged from peri-nucleolus stage to secondary yolk globule stage, and the testis of two males was at midmeiotic stage;Crustacean was the major prey for conger of small length, and food source shift to fish with somatic growth. The results showed substantial differences from previous studies in Japan and Korean waters, as well as from China seas in the 1980s, suggesting potential spatiotemporal changes in the biological characteristics of C. myriaster. This study may improve current understanding of the fishery biology of C. myriaster in the Yellow Sea and East China Sea.

Impacts of sample size for stomach content analysis on the estimation of ecosystem indices

This study used Ecopath model of the Jiaozhou Bay as an example to evaluate the effect of stomach sample size of three fish species on the projection of this model. The derived ecosystem indices were classified into three categories:(1) direct indices, like the trophic level of species, influenced by stomach sample size directly;(2)indirect indices, like ecology efficiency(EE) of invertebrates, influenced by the multiple prey-predator relationships;and(3) systemic indices, like total system throughout(TST), describing the status of the whole ecosystem. The influences of different stomach sample sizes on these indices were evaluated. The results suggest that systemic indices of the ecosystem model were robust to stomach sample sizes, whereas specific indices related to species were indicated to be with low accuracy and precision when stomach samples were insufficient.The indices became more uncertain when the stomach sample sizes varied for more species. This study enhances the understanding of how the quality of diet composition data influences ecosystem modeling outputs. The results can also guide the design of stomach content analysis for developing ecosystem models.

Evaluating the eff ect of input variables on quantifying the spatial distribution of croaker Johnius belangerii in Haizhou Bay,China

A habitat model has been widely used to manage marine species and analyze relationship between species distribution and environmental factors.The predictive skill in habitat model depends on whether the models include appropriate explanatory variables.Due to limited habitat range,low density,and low detection rate,the number of zero catches could be very large even in favorable habitats.Excessive zeroes will increase the bias and uncertainty in estimation of habitat.Therefore,appropriate explanatory variables need to be chosen first to prevent underestimate or overestimate species abundance in habitat models.In addition,biotic variables such as prey data and spatial autocovariate(SAC)of target species are often ignored in species distribution models.Therefore,we evaluated the eff ects of input variables on the performance of generalized additive models(GAMs)under excessive zero catch(>70%).Five types of input variables were selected,i.e.,(1)abiotic variables,(2)abiotic and biotic variables,(3)abiotic variables and SAC,(4)abiotic,biotic variables and SAC,and(5)principal component analysis(PCA)based abiotic and biotic variables and SAC.Belanger’s croaker Johnius belangerii is one of the dominant demersal fish in Haizhou Bay,with a large number of zero catches,thus was used for the case study.Results show that the PCA-based GAM incorporated with abiotic and biotic variables and SAC was the most appropriate model to quantify the spatial distribution of the croaker.Biotic variables and SAC were important and should be incorporated as one of the drivers to predict species distribution.Our study suggests that the process of input variables is critical to habitat modelling,which could improve the performance of habitat models and enhance our understanding of the habitat suitability of target species.

Evaluating impacts of pulse fishing on the effectiveness of seasonal closure

Seasonal fishing closures are often used in fisheries management to conserve overfished stocks.As one of the unintended consequences,fishermen often contend for maximizing catches immediately after reopening fisheries.The resultant large catch landings in a short time period(i.e.,pulse fishing)may undermine the benefit of closure.We implemented an end-to-end model OSMOSE-JZB(Object-oriented Simulator of Marine ec OSystem Exploitation OSMOSE)modelling ecosystem in the Jiaozhou Bay located in China to evaluate the impact of pulse fishing on the effectiveness of seasonal closure at levels of fish community,population,and individual.Our study demonstrated that the three-month closure was successful in conserving fish stocks.There were small variations on ecological indicators(i.e.,total biomass of the community,mean trophic level of the community,mean trophic level of the catch,and Shannon-Wiener biodiversity index)when pulse fishing occurred.Pulse fishing seemed not to result in a great shift in community structure.Compared to other species,the biomass of two large predatory fishes were more susceptible to pulse fishing.Pulse fishing could change the pressure of predators to fish stocks via food webs,especially for young individuals.Our simulations indicate that we can improve the effectiveness of seasonal closure by managing pulse fishing.Although the results derived in this study may be specific to the target ecosystem,the general approach is applicable to other ecosystems when evaluating fishing impacts.

中国海洋生物地理学研究进展和热点:物种分布模型及其应用

海洋生物地理学是研究海洋生物地理分布时空格局、形成过程及其影响因素的学科,对于生物多样性保护、生态系统结构稳定以及海洋资源可持续利用具有重要支撑作用。物种分布模型是研究物种地理分布变化的重要手段,其发展对于海洋生物地理学的发展具有重要的推动作用。本文综述了中国海洋生物地理学的发展历程和现状;提供了物种分布模型的定义、类型、发展历程及相关重要数据库;整理分析了中国海洋生物地理学研究中的部分热点问题,主要涉及全球变化下海洋生物分布区变化、海洋生物谱系地理结构、海洋生物入侵、种群连通性、海洋保护规划、海洋生态修复、海洋生物对极端环境适应、海洋渔业资源管理和海水养殖规划等。本文同时对海洋生物地理学研究的未来发展进行了展望,强调了改进物种分布模型的重要性,提出要强化环境和生物数据库建设,加强海洋生物地理学与其他学科的交叉融合等。本综述对于今后我国海洋生物地理学的研究,尤其是物种分布模型在相关领域中的应用具有参考价值。

Feeding ecology of Japanese Spanish mackerel(Scomberomorus niphonius)along the eastern coastal waters of China

Feeding activities provide necessary nutrition and energy to support the reproduction and development of fish populations.The feeding ecology and dietary plasticity of fish are important factors determining their recruitment and population dynamics.As a top predator,Japanese Spanish mackerel(Scomberomorus niphonius)supports one of the most valuable fisheries in China.In this study,the feeding ecology and diet composition of Japanese Spanish mackerel spawning groups were analysed based on samples collected from six spawning grounds along the eastern coastal waters of China during spring(March to May)in 2016 and 2017.Both stomach contents and stable isotope analysis were conducted.Stomach content analysis showed that spawning groups of Japanese Spanish mackerel mainly fed on fish,consuming more than 40 different prey species.Diets were significantly different among sampling locations.The most important prey species were Stolephorus in Fuzhou,Japanese anchovy Engraulis japonicus in Xiangshan,Euphausia pacifica in Lüsi,sand lance Ammodytes personatus in Qingdao and Weihai,and Leptochela gracilis in Laizhou Bay.Stable isotope analysis showed that the trophic level of Japanese Spanish mackerel was relatively high and generally increased with latitude from south to north.In the 1980 s,the diet of Japanese Spanish mackerel was dominated solely by Japanese anchovies in the eastern coastal waters of China.The results in the present study showed that the importance of Japanese anchovies declined considerably,and this fish was not the most dominant diet in most of the investigated waters.Both the spatial variations in diet composition and changes in the dominant diet over the long term indicated the high adaptability of Japanese Spanish mackerel to the environment.Combining the results of stomach analysis and stable isotope analysis from different tissues provided more comprehensive and accurate dietary information on Japanese Spanish mackerel.The study provides essential information about the feeding ecology of Japanese Spanish mackerel and will benefit the management of its populations in the future.

Performance evaluation of fixed-station sampling design for a fishery-independent survey with multiple objectives

Fixed-station sampling design was widely used in fishery-independent surveys because of its characteristics of convenient sampling station setting,but the non-probabilistic(fixed)nature made it more uncertainty of drawing inferences on population.The performance of fixed-station sampling design for multispecies survey has not been evaluated,and we are uncertain if the design could detect the temporal trends of different populations in multispecies fishery-independent survey.In this study,spatial distribution of abundance indices for three species with different spatial distribution patterns including small yellow croaker(Larimichthys polyactis),whitespotted conger(Conger myriaster)and Fang’s blenny(Enedrias fangi)were simulated using ordinary kriging interpolation as the“true”population distribution.The performance of fixed-station sampling design was compared with simple random sampling design by resampling the simulated“true”populations in this simulation study.The results showed that the fixed-station sampling design had the power to detect the seasonal trends of species abundance.The effectiveness of fixed-station sampling design were different in different species distribution patterns.When the species had even distribution,fixed-station sampling design could get high quality abundance data;when the distribution was uneven with heterogeneity or patchiness,fixed-station sampling design tended to underestimate or overestimate the abundance.Evidently,the estimates of abundance index based on the fixedstation sampling design must be calibrated cautiously while applying them for fisheries stock assessment and management.This study suggested that fixed-station sampling design could catch the temporal dynamics of population abundance,but the abundance estimates from the fixed-station sampling design could not be treated as the absolute estimates of populations.

Evaluating multispecies survey designs using a joint species distribution model

Survey designs should be efficient as marine survey programs are usually expensive and time-consuming;however,surveys have rarely been evaluated for multiple species.In the present study,we evaluated multispecies fisheries surveys with respect to three influential factors,i.e.,sampling methods,estimation methods and sample size.A joint species distribution model(JSDM)developed in north Yellow Sea,China was used as the operating model to simulate the spatial distribution of multiple species simultaneously.We examined the precision of multispecies abundance estimation using diverse sampling methods[random sampling(RDS),systematic sampling(SYS),stratified random sampling(SRS),generalized random-tessellation stratified sampling(GRT)and spatial coverage sampling(SPC)],estimation methods[arithmetic mean(Arm),universal kriging(Ukr),multivariate distribution model(Mvd),and boral model(Brm)],and a range of sample sizes(from 30 to 300).The results showed significant differences in estimation among sampling methods,where GRT and SYS yielded less relative absolute bias(RAB)over all and RDS showed the least precision.Regarding estimation methods,Mvd and Arm showed the best performances and Brm yielded the least precision.Significant interactions existed between sampling and estimation methods.Arm worked best with GRT,likewise Mvd with SYS and Ukr with SPC.SPC and Mvd showed the best performances for a small sample size(N=30),and all sampling and estimation methods provided similar results for a large sample size(N=300).Generally,doubling sample size resulted in a decrease of RAB by 0.097 on average,a rate depending on species,sampling and estimation methods.This study contributed to an integrative framework for evaluating designs of multispecies fisheries surveys.

Comparison of sampling effort allocation strategies in a stratified random survey with multiple objectives

Stratified random survey is commonly used to estimate abundance indices of fish populations in multispecies survey,providing reliable data for stock assessment and fisheries management.In some cases,however,the sample size is relatively small because of the limitation of survey cost or other factors.The allocation methods of sampling efforts among strata in stratified random surveys with small sample size may need adjustment compared with traditional approaches.In this study,two sampling stations were allocated to each stratum first and then the remaining sampling units were allocated among strata using five traditional allocation methods.In order to distinguish them from traditional methods,we called them adjusted methods in this study.A simulation study was conducted to compare the performances of different allocation strategies of sampling efforts in a stratified random survey for estimating abundance indices of multiple target species.Relative estimation error(REE)and relative bias(RB)were used to measure the precision and accuracy of estimates of abundance indices under different allocation schemes of sampling efforts in the multispecies survey.The performances of different allocation schemes in estimating abundance indices varied greatly for different species over different seasons.The adjusted Neyman allocation scheme could significantly reduce the REE and RB of estimates of abundance index for single species survey.For multiple species surveys,the adjusted average-Neyman allocation method,the adjusted Yate allocation method,the adjusted proportional allocation method and current allocation method had relatively high accuracy and precision of estimates of abundance indices for four species in terms of the total_(REE) and total_(RB).Though the adjusted average-Neyman allocation scheme did not always have the best performance,it was the optimal one considering the accuracy and precision of estimates of abundance indices for all species simultaneously.The allocation of sampling efforts among strata in stratified random surveys targeting for estimating abundance indices of multiple species should comprehensively consider the variance of abundance of different species in stratum and the seasonal changes.

Special topic-Fisheries monitoring:Eyes in the ocean:Innovative and multidimensional solutions to monitor fisheries

In the era of climate change,the visibility of environmental changes dictates public attention.Pictures of untamable bushfires,intense hurricanes,collapsing ice sheets are all gripping images that alarm us and urge us to take action.The ocean,however,gives us less such visuals;the changes that are taking place there are often abstract and hidden.However,the environmental challenges in the ocean are less visible but no less grave;they come from multiple sources:pollution,plastic waste,ocean surface warming,ocean acidification(IPCC,2019).But the most direct and age-old impact we exert on the ocean and its ecosystems is overfishing(FAO,2018).