The Sample Size Necessary to Assess Changes in Fish Biomass – A Reply

Abstract. In this article, we compare biomass results for three fish species (Symphodus ocellatus, Serranus scriba and Diplodus annularis) between May 1993 and May 2000 in the Bonifacio Strait Marine Reserve. The data were obtained from seagrass meadows using the visual census method. A randomised procedure was used to test the null hypothesis of biomass equality. To this end, decreasing sample sizes were used for the May 2000 data to determine the sample size necessary to reject the null hypothesis (5 % level). For Symphodus ocellatus, 10 samples (fixed points) were necessary to reject the null hypothesis, whereas 15 samples were required for Diplodus annularis. The null hypothesis was never rejected for Serranus scriba regardless of the sample size considered.
These results are discussed, as is the power statistical analysis proposed by Ortiz (2002; this issue, pp. 1 – 9) , to determine the optimal sample size necessary to detect perturbation effects in the marine environment. Although the power analysis is certainly a useful method to determine the sampling design in assessment the effect of marine reserves on fish biomass, the magnitude of the marine reserve effect remains problematic in ecology.     

Precision in systematic trawl surveys as assessed from replicate sampling by parallel trawling off Namibia

Following independence of Namibia in 1990, abundance of that country's hake stocks was monitored by trawl surveys conducted by the Norwegian F.R.V. Dr Fridtjof Nansen using a systematic survey design. Precision from such designs is considered to be better than with a random design, but it cannot be quantified through standard analysis of variance. In 1998 and 1999, the trawl surveys were duplicated in full using a commercial trawler that operated in parallel with the Dr Fridtjof Nansen. Both vessels had the same fishing gear and rigging. In both years, biomass and distribution patterns estimated by the vessels were similar. The paired datasets collected from the surveys were used to analyse variability in the point samples. Absolute differences between pairs of catches were, in general, proportional to the catch, but they varied randomly following a normal distribution around the catch level. Local variability was analysed as random noise, modelled and later reapplied on single vessel data series to evaluate the effect of sample size and replicates on the survey mean. In a survey with 200 stations, local sample noise accounted for about 4% of the variability in the survey mean. Alternatively, running series of simulated surveys by bootstrapping on pairs of catch data gave similar results, when a small systematic vessel effect was adjusted for. The main statistical techniques applied were less susceptible to outlier catches than straightforward correlations or regressions and could therefore, perhaps with some advantage, also be used to estimate the vessel factor when intercalibrating trawl survey vessels.     

Landscape ecology of interactions between seagrass and mobile epifauna: The matrix matters

There is increasing interest among ecologists about how the type of matrix surrounding a habitat patch influences the organisms living in that patch. This question is virtually unstudied in marine systems. In this paper I show that the mobile faunal assemblage in seagrass patches does depend on the surrounding matrix. Faunal assemblages in patches of Posidonia surrounded by sand are different than in those surrounded by Heterozostera, another seagrass, having more than double the abundance of both amphipods and polychaetes. However, the differences are not simply due to spillover from the matrix habitat, but rather are an emergent property of the patch context that cannot be predicted. Posidonia surrounded by sand actually has an assemblage that is intermediate between Heterozostera and Posidonia surrounded by Heterozostera. Differences in habitat structure do not account for this pattern, as seagrass biomass did not vary, and the same result was found in artificial seagrass. The faunal assemblage did not vary depending on the location within the patch (edge or centre) for Heterozostera, Posidonia or artificial seagrass. Patch size, however, did have an effect for Heterozostera, with smaller patches having 2–3 times as many isopods per sample as large, but less than half the number of some amphipod families. These results suggest that the landscape context is as important in marine systems as it is known to be in terrestrial systems.     

Estimating periphyton standing crop in streams: a comparison of chlorophyll a sampling and visual assessments

Periphyton standing crop is often measured as chlorophyll a but there is increasing interest in using visual assessments of periphyton cover. Visual methods are rapid and require no sample analysis, but can be regarded as subjective with high inter-operator variability. We investigated variability in periphyton standing crop across operators, rivers and time, as determined visually and from chlorophyll a measurements. We found that visual assessments (comprising percentage cover estimates of up to eight periphyton categories) distinguished sites and occasions as effectively as chlorophyll a. Furthermore, an estimate of chlorophyll a could be derived from the visual assessments. Because our surveys were conducted in only three rivers, general applicability of the derivation of estimated chlorophyll a from the visual assessments warrants further investigation. The current recommendation of 20 views was sufficient for realistic visual assessments of average cover. Overall, our results indicated that inter-operator variability in visual assessments need not be a major concern, given adequate training.     

The latitudinal and bathymetric ranges of marine fishes: a global analysis to test the application of Rapoport′s Rule

We evaluated the applicability of Rapoport’s Rule (RR) to the marine bony fishes of the world. The biogeographical pattern predicted by RR has been the subject of a large number of studies, some supporting it and some not. In this exercise, we attempted to generate results free of biases from taxonomy or geographic scale. The study area encompassed all world oceans. Our analysis was based on secondary data. We tested the relationship of the geographic range and the bathymetric range to the latitude gradient, using the method of Stevens (American Naturalist 1989, 133 , 240–256) . We compared all known species of marine bony fishes together, and performed a second analysis limited to the non‐pelagic species. Our results were generated from a databank including 13,957 species. The results indicated that RR (latitudinal range) is valid for 11 regions. However, there were exceptions for the northern part of the Eastern Atlantic and Pacific, the northern part of the Western Indian, and the Arctic Oceans. The analyses of RR (bathymetric range) were supported for 13 regions. We found only one exception, for the northern part of the Western Indian Ocean. The applicability of RR and its relationship to energy supply, represented here by the temperature variability of the marine water masses, was clear. The exceptions appeared when some feature associated with the distribution of water masses was superposed on the latitudinal gradient, or when some life trait characteristic of a taxon, for instance body size, was confounded with the effect of climate variability.     

Sediment–water interactions of natural oestrogens under estuarine conditions

The natural human female hormones oestrone and 17β-oestradiol have been implicated in the disruption of endocrine systems in some wildlife adjacent to sewage effluents. The sorption behaviour of these two compounds under estuarine conditions was studied by spiking either 2.55 μg of oestrone or 2.65 μg of 17β-oestradiol in kinetic experiments. In equilibrium experiments, 3 ng of oestrone or 3.2 ng of 17β-oestradiol was added in each of the centrifuge tubes. Sorption onto sediment particles was relatively slow, with sorption equilibrium being reached in about 70 and 170 h for oestrone and 17β-oestradiol, respectively. The effects of a variety of environmental parameters on sorption were studied including salinity, sediment concentration (SC), the presence of a third phase, particle size and, also, surfactant concentrations. Results show that although salinity did not induce any statistically significant effect on the sorption of 17β-oestradiol, it did statistically enhance the sorption of oestrone, and a salting constant of 0.3 l mol⁻¹ was derived. The partition coefficient for both compounds decreased with increasing sediment concentration, a phenomenon that has been widely reported and attributed to the presence of colloids (which could enhance dissolved concentrations). In this paper, the true partition coefficients for sediment particles (K_p^true) and colloidal particles (K_c^true) have been calculated, and a K_p^true value of 141 and 102 ml g⁻¹ was obtained for oestrone and 17β-oestradiol, respectively. In addition, K_c^true values for oestrone (222×10² ml g⁻¹) and 17β-oestradiol (135×10² ml g⁻¹) were two orders of magnitude higher than their respective K_p^true values, suggesting that the colloidal particles are significantly stronger sorbents for natural oestrogens than sediment particles. Particles of different sizes were found to have different partition coefficients due to the strong relationships between partition coefficients for the two compounds and particulate organic carbon (POC) contents and specific surface areas (SSAs). The presence of a surfactant was shown to reduce the partition coefficients for the two compounds, although its concentrations being used were higher than those normally found in the natural environment.     

High concentrations of marine snow and diatom algal mats in the North Pacific Subtropical Gyre: Implications for carbon and nitrogen cycles in the oligotrophic ocean

A video plankton recorder (VPR) and a remotely operated vehicle (ROV) were utilized on three cruises in the oligotrophic North Pacific Subtropical Gyre (NPSG) between 1995 and 2002 to quantify the size and abundance of marine snow and Rhizosolenia diatom mats within the upper 305 m of the water column. Quantitative image analysis of video collected by the VPR and an ROV-mounted particle imaging system provides the first transect of marine snow size and abundance across the central North Pacific Gyre extending from 920 km NW of Oahu to 555 km off Southern California. Snow abundance in the upper 55 m was surprisingly high for this oligotrophic region, with peak values of 6.0–13.0×10³ aggregates m⁻³ at the western- and easternmost stations. At stations located in the middle of the transect (farthest from HI and CA), upper water column snow abundance displayed values of 0.5–1.0×10³ aggregates m⁻³. VPR and ROV imagery also provided in situ documentation of the presence of nitrogen-transporting, vertically migrating Rhizosolenia mats from the surface to >300 m with mat abundances ranging from 0 to 10 mats m⁻³. There was clear evidence that Rhizosolenia mats commonly reach sub-nutricline depths. The mats were noted to be a common feature in the North Pacific Gyre, with the lower salinity edge of the California Current appearing to be the easternmost extent of their oceanic distribution. Based on ROV observations at depth, flux by large (1.5 cm) mats is revised upward 4.5-fold, yielding an average value of 40 μmol N m⁻² d⁻¹, a value equaling previous estimates that included much smaller mats visible only to towed optical systems. Our results suggest that the occurrence across a broad region of the NPSG of particulate organic matter production events represented by high concentrations of Rhizosolenia mats, associated mesozooplankton, and abundant detrital marine aggregates may represent significant stochastic components in the overall carbon, nitrogen, and silica budgets of the oligotrophic subtropical gyre. Likewise, their presence has important implications for the proposed climate-driven, ecosystem reorganization or domain shift occurring in the NPSG.     

Trade-off between safety and feeding in the sea anemone Anthopleura aureoradiata

After a physical disturbance, such as a predatory stimulus, several marine organisms withdraw into a refuge to ensure survival. Consequently, they cannot engage in other essential activities. Re-emergence reflects a trade-off between avoidance of disturbance and lost feeding opportunities. The New Zealand sea anemone (Anthopleura aureoradiata) avoids predators by withdrawing its tentacles into its upper body. In a field experiment, responses of 480 A. aureoradiata to a physical disturbance were studied at three locations to assess the influence of the following factors on hiding time: physical contact with a neighbour; density of nearby conspecifics; total density per pool; and body size. Withdrawal times were significantly shorter for anemones in contact with a conspecific, occurring in high-density patches or of larger size. These patterns are consistent with rapid re-emergence being stimulated by competition for food from immediate neighbours and/or a lower risk of predation when in groups (i.e. a dilution effect), with larger size also potentially reducing vulnerability to predation.     

Do fecundity‐time models really predict extreme optimal egg sizes?

In 1973, Vance published a mathematical model that explains differences in reproductive strategies such as planktotrophy and lecithotrophy found in marine invertebrates. The original model Vance (The American Naturalist, 107 , 339–352) traded fecundity against development time and resulted in a U‐shaped fitness function suggesting that only extreme levels of egg provisioning (referred to as egg size) are evolutionarily stable. Several later models independently predicted evolutionary stability of extreme egg sizes, and this has been commonly interpreted as a proof that the Vance model, albeit oversimplified, reflects the fundamental evolutionary principle of disruptive egg size evolution in marine invertebrates. However, empirical studies of the last several decades have illustrated that intermediate levels of egg provisioning are commonly observed in nature. As a result, the recent modifications of the original Vance model have attempted to explain evolution of such intermediate strategies. Here we re‐examine the earlier fecundity‐time models of the ‘Vance family’ and show that only Vance and its minor modification by Grant (The American Naturalist, 122 , 549–555) invariably predict evolutionary stability of extreme egg sizes. Most modifications of the original Vance model in fact can predict at least some intermediate optimal egg sizes and reproductive strategies. There is no reason to assume that selection towards extreme egg sizes is the driving force in the evolution of marine invertebrate life histories.     

Stable Nitrogen Isotopes in Fish: Literature Synthesis on the Influence of Ecotonal Coupling

Stable nitrogen isotopes have customarily been used to delineate trophic position with only infrequent regard to source variability. A compilation of literature data supports a previous tentative hypothesis that marine fish are enriched in¹⁵N relative to those inhabiting freshwaters. Estuarine and anadromous fish were also confirmed to have intermediate δ¹⁵N values depending on their respective time spent feeding in either fresh- or salt water. The use of fish δ¹⁵N as a measure of continental-marine coupling in complex coastal environments will therefore provide additional support for analyses based on the more traditionally utilized isotopes δ³⁴S or δ¹³C.     

Trophic ecology of Pomatoschistus microps within an intertidal bay (Roscoff,France), investigated through gut content and stable isotope analyses

The diet of Pomatoschistus microps has been studied using both gut content and stable isotope analyses. In the Roscoff Aber Bay (Brittany, France), this fish is commonly found on sandy muddy intertidal flats. Gut content analyses were also interpreted using trophic indices. Owing to the large diversity of prey consumed, these indices emphasised the opportunistic feeding behaviour of P. microps. Here, this species fed mainly on endofauna with meiofauna being of high relative importance. The main biotic components of its trophic habitat, characterized by δ¹³C and δ¹⁵N, provided evidence of a major trophic pathway based on drift Enteromorpha sp. Trophic positions estimated by both diet analyses and isotopic analyses led to similar results. In this bay, P. microps is a first‐order predator with a low degree of omnivory. Despite a preferential consumption of the amphipod Corophium arenarium, we assumed that this goby behaves as a generalist feeding on a uniform variety of endofauna taxa.     

What do egg size distributions in marine invertebrates tell us about validity of fecundity‐time models?

An early mathematical model of egg size evolution in marine invertebrates (American Naturalist, 1973a, 107 , 353) traded fecundity against development time and resulted in a U‐shaped optimality curve suggesting that only extremes of egg sizes should be evolutionary stable. It was later suggested that the model predicts a bimodal distribution of egg sizes and therefore can be tested by examining histograms of egg sizes. More recent fecundity‐time models (American Naturalist, 1997, 150 , 48; American Naturalist, 2000, 156 , 175) came up with predictions of intermediate optimal egg sizes; however, comparisons of optimality curves with histograms of egg sizes remained central to model testing. Here the question whether the bimodality assumption has any basis has been addressed theoretically. Simulation modeling and analysis show that egg size distributions generally do not match the optimality curves produced by fecundity‐time models and that the shape of predicted frequency distributions depends on both the optimality curve and other model parameters. Therefore, egg size frequency distributions alone cannot be used for model testing.     

Analyses of stomach contents and stable isotopes reveal food sources of estuarine detritivorous fish in tropical/subtropical Taiwan

Detritivorous fish generally refers to fish that primarily ingest unidentified organic detritus. We analyzed stomach contents in combination with stable isotopes to trace and compare the food sources of the large-scale mullet Liza macrolepis and other detritivorous fish species in subtropical mangrove creeks and a tropical lagoon in Taiwan. The volume of organic detritus always contributed >50% of the stomach content of L. macrolepis in the two habitats. However, consumed items were distinct between the two habitats and corresponded to the types in which they reside. The consumed items in the lagoon were more diverse than those observed in the mangroves. In the mangroves, the diet composition of L. macrolepis was primarily determined by season, not by body size. In the lagoon, there were no clear seasonal or size-dependent grouping patterns for the diet composition. There were significant seasonal and spatial variations in δ¹³C and δ¹⁵N values of potential food sources and L. macrolepis. However, neither δ¹³C nor δ¹⁵N values of L. macrolepis were correlated with fish body size. Joint analyses of stomach contents and stable isotopes indicated that benthic microalgae on sediments were the most important assimilated food in both seasons for the dominant detritivorous fish in the mangroves, whereas a greater reliance on microalgal and macroalgal periphyton on oyster-culture pens was observed in the lagoon. Mangrove and marsh plants and phytoplankton, which are mostly locally produced within each habitat, were of minor importance in the assimilated food.     

Spatial and temporal variations in sediment grain size in tidal wetlands, Yangtze Delta: On the role of physical and biotic controls

To examine the spatial and temporal variability of sediment grain size in exposed tidal wetlands with ample sediment supply, we sampled sediments and measured hydrodynamics, accretion/erosion rates, and vegetation characteristics in the Yangtze Delta. Sediment grain size exhibited a landward/upward decreasing trend. This trend is mainly attributed to attenuation of hydrodynamics. A 630-day series of daily surface sediment sampling at a fixed site on an unvegetated intertidal flat revealed significant seasonal and storm-cyclic changes in grain size. This temporal variability was related to alternating accretion/erosion events, with erosion associated with coarser grain size. Such temporal dynamics were not present in vegetation, where sediment remained fine grained throughout the year. In the marsh, vegetation cover enables the trapping of fine-grained sediments in the following ways: (a) adherence of suspended sediments onto plants; (b) deposition of suspended sediments stimulated by attenuation of hydrodynamics through plant obstruction; and (c) prevention of resuspension of fine-grained deposits due to the protection of the plant canopy. The influence of vegetation on sediment grain size was clearly seen when comparing sediment trapped by different vegetation types and seasonal patterns of trapped sediment on different vegetation canopy densities. The relatively high plant biomass of the recently introduced Spartina alterniflora enhanced the trapping effect, whereas plant degradation due to buffalo grazing reduced the trapping effect. We conclude that for exposed tidal wetlands with ample sediment supply such as the Yangtze Delta, the spatial and temporal variability of sediment grain size is governed predominantly by physical controls on the unvegetated flat and predominantly by biophysical interaction of hydrodynamics and vegetation in the salt marsh, rather than by sediment supply.     

Effects of pulsed nutrient inputs on phytoplankton assemblage structure and blooms in an enclosed coastal area

The response of phytoplankton assemblage structure to terrestrial nutrient inputs was examined for the Gulf of Kalloni in the Northern Aegean Sea, a productive semi-enclosed coastal marine ecosystem. The study was focused on a typical annual cycle, and emphasis was placed on the comparative analysis between blooms developing after significant nutrient inputs from the watershed, and naturally occurring blooms. Baseline information was collected on a monthly basis from a network of stations located in the oligotrophic open sea and the interior and more productive part of the embayment. Intensive sampling was also carried out along a gradient in the vicinity of a river which was the most important source of freshwater and nutrient input for the Gulf. Phytoplankton assemblage structure was analyzed from 188 samples using diversity indices (Shannon and Average Taxonomic Distinctness), multivariate plotting methods (NMDS), multivariate statistics (PERMANOVA), and canonical correspondence analysis (CCA). Three characteristic assemblages were recognized: (1) an autumn assemblage developed under nutrient depleted conditions, having low diversity due to the dominance of two small diatoms, (2) a winter bloom of the potentially toxic species Pseudo-nitzschia calliantha occurring immediately after a nutrient peak and characterized by very low diversity, and (3) a naturally occurring early summer bloom of centric diatoms with relatively high diversity. The results of the study support the view that moderate nutrient inputs may have a beneficial effect on the functioning of coastal ecosystems, stimulating the taxonomic diversity through the growth of different taxonomic groups and taxa. On the other hand, a sudden pulse of high nutrient concentrations may greatly affect the natural succession of organisms, have a negative effect on the diversity through the dominance of a single species, and can increase the possibility of a harmful algal bloom development.     

A note on the relationships between organic matter and some geotechnical properties of a marine sediment

Abstract

An analysis of the relationship between organic matter and liquid and plastic limits, and grain‐specific gravity of a marine sediment was accomplished by making a stepwise adjustment in the organic content of that sediment. The sample used was from Santa Barbara Basin (off southern California) and is typical of fine‐grained marine sediments: it is a clayey silt with a common suite of minerals and other constituents. During the experiment, texture and composition (except organic content) were constant; only the quantity of natural, indigenous organic matter was changed.

A strong linear relationship exists between the independent variable, the amount of organic carbon present in the sample, and the dependent variables. Liquid limit, plastic limit, and plasticity index all increased with increasing organic content over the range studied (0.57–3.20% organic carbon). Grain‐specific gravity decreased. All had linear correlation coefficients (r) greater than |0.90| and r ² values greater than 90%, except the plasticity index (83%).

Comparing the results of regression analyses from this and several similar studies shows that although there is good qualitative agreement, there are quantitative inconsistencies. In particular there is considerable overall variability in the regression coefficients. Among studies on marine sediments the inconsistencies are less pronounced, yet still evident. The increase in liquid limit as organic carbon increased by 1 % sediment dry weight ranged from 9 to 28% water content; in the plastic limit the range was from 4 to 18%. However, in these marine studies regression coefficients are relatively close in value in some cases, levels of significance of the regressions are high in most cases, and in all cases the relationships appear to be linear over the range of organic carbon percentage studied. Finally, we believe that a relatively clear relationship between plasticity and organic carbon begins to emerge when the latter exceeds a value of 2%.     

Influence of sample size on ecological status assessment using marine benthic invertebrate‐based indices

Routine monitoring of the quality status of water bodies demands the best cost‐benefit relation and sample‐size reduction is therefore welcomed. However, great caution is needed because such reduction affects the accuracy and variation of the results. In the present study we tested the influence of sample size (number of replicate samples) on reference condition values and within‐sample ecological quality ratio (EQR) variability of six commonly used ecological indices (taxa richness, Shannon–Wiener diversity, AMBI, Medocc index, Bentix and M‐AMBI). Analysis of soft‐bottom benthic invertebrate data from Slovenian coastal waters showed that sample size influenced the reference condition values of richness/diversity indices (taxa richness and Shannon–Wiener diversity) but not of the sensitivity/tolerance indices (AMBI, Medocc index, Bentix). Increased sample size decreased the within‐sample EQR variability and, concomitantly, increased the accuracy of site ecological status classification for all indices. The size of EQR variability differed depending on the index used. EQR variability of M‐AMBI, an index composed of metrics with different within‐sample EQR variability, was statistically the same as that of the metric with the lowest within‐sample EQR variability. Whether this is a common principle for multimetric indices remains to be confirmed. Based on these results, the use of at least three replicates is suggested to obtain reliable measures of reference condition and EQR for the assessment of ecological status. This level of replication is particularly necessary in areas with high diversity and environmental patchiness, and when richness/diversity measures and indices that include these measures are used.     

Variability of the apparent quantum efficiency of CO photoproduction in the Gulf of Maine and Northwest Atlantic

The photochemical oxidation of colored, dissolved organic matter (CDOM) is important for carbon cycling in the ocean. This oxidation process produces a number of products, including carbon monoxide (CO). While the photochemical production efficiency of CO (apparent quantum yield, AQY, defined in terms of CDOM absorbance) has been reported to be similar for many water types, a full evaluation of the observed natural variability in CO AQY requires additional study. Here we use a polychromatic irradiation system to determine twenty AQY spectra at sea on fresh samples ranging from the near coastal waters of the Gulf of Maine to the offshore waters of the Northwest Atlantic. Despite the geographic variability of these marine samples the AQY of CO production in the Gulf of Maine and Northwest Atlantic exhibited only a small degree of variability, none of which was not correlated with measured environmental parameters. Consequently, a single aggregate AQY spectrum _λ = e^{(−(9.134+0.0425(λ−290)))+e(−(11.316+0.0142(λ−290)))} was found to adequately represent the entire data set. Significantly, the accuracy of an AQY spectrum determined using this multispectral/statistical technique was confirmed with data obtained from a monochromatic irradiation technique on a single open ocean sample. Taken together, the AQY spectra determined in this study were similar in magnitude and shape to those previously published for marine samples and, overall, were somewhat lower than those previously reported for freshwater studies.     

Food web structure of two Mediterranean lagoons under varying degree of eutrophication

The food web structure and functioning of two north-western Mediterranean lagoons exhibiting contrasting degrees of eutrophication and marine influences were compared through δ¹³C and δ¹⁵N analysis of major potential food sources and consumers. The Lapalme Lagoon is well preserved and has kept a natural and temporary connection with the open sea. Conversely, the Canet Lagoon is heavily eutrophicated and its water exchange with the open sea has been artificially reduced. In Lapalme, all potential food sources and consumers exhibited δ¹⁵N values indicative of pristine coastal areas. Suspended particulate organic matter (POM) and sediment organic matter (SOM) pools seemed to constitute the main food sources of most primary consumers. Both primary producers and all consumers were much more ¹⁵N-enriched (by  10‰) and more ¹³C-depleted in Canet than in Lapalme. This reflected: (1) the assimilation of important amounts of anthropogenic nitrogen in the food web, and (2) a marked and uniform influence of ¹³C-depleted allochtonous sources of carbon. Based on the mean δ¹⁵N of primary consumers, we found rather similar food web lengths in both lagoons with top consumers at trophic levels 3.6 and 4.0 in Canet and Lapalme, respectively. However, the eutrophication of the Canet Lagoon resulted in a simplification of the food web structure (i.e., a single trophic pathway from a ¹⁵N-enriched fraction of the SOM pool to top predators) compared to what was observed in Lapalme Lagoon where additional ¹³C-enriched food sources played a significant trophic role. Moreover, some consumers of Canet tended to exploit primary producers to a larger extent (and thus to exhibit lower trophic levels) than in Lapalme.     

Implementing geostatistical analysis to study spatio-temporal distribution patterns of swimming crabs(Portunus trituberculatus)

This study attempted to compare the performance of local polynomial interpolation,inverse distance weighted interpolation,and ordinary kriging in studying distribution patterns of swimming crabs.Cross-validation was used to select the optimum method to get distribution results,and kriging was used for making spatial variability analysis.Data were collected from 87 sampling stations in November of 2015(autumn) and February(winter),May(spring) and August(summer) of 2016.Results indicate that swimming crabs widely distributed in autumn and summer:in the summer,they were more spatially independent,and resources in each sampling station varied a lot;in the winter and spring,the abundance of crabs was much lower,but the individual crab size was bigger,and they showed the patchy and more concentrative distribution pattern,which means they were more spatially dependent.Distribution patterns were in accordance with ecological migration features of swimming crabs,which were affected by the changing marine environment.This study could infer that it is applicable to study crab fishery or even other crustacean species using geostatistical analysis.It not only helps practitioners have a better understanding of how swimming crabs migrate from season to season,but also assists researchers in carrying out a more comprehensive assessment of the fishery.Therefore,it may facilitate advancing the implementation in the pilot quota management program of swimming crabs in northern Zhejiang fishing grounds.