Zooplankton and Karenia brevis in the Gulf of Mexico

Blooms of the toxic dinoflagellate Karenia brevis are common in the Gulf of Mexico, yet no in situ studies of zooplankton and K. brevis have been conducted there. Zooplankton abundance and taxonomic composition at non-bloom and K. brevis bloom stations within the Ecology of Harmful Algal Blooms (ECOHAB) study area were compared. At non-bloom stations, the most abundant species of zooplankton were Parvocalanus crassirostris, Oithona colcarva, and Paracalanus quasimodo at the 5-m isobath and P. quasimodo, O. colcarva, and Oikopleura dioica at the 25-m isobath. There was considerable overlap in dominance of zooplankton species between the 5 and 25-m isobaths, with nine species contributing to 90% of abundance at both isobaths. At stations within K. brevis blooms however, Acartia tonsa, Centropages velificatus, Temora turbinata, Evadne tergestina, O. colcarva, O. dioica, and P. crassirostris were dominant. Variations in abundance between non-bloom and bloom assemblages were evident, including the reduction in abundance of three key species within K. brevis blooms.     

Direct effects of Daphnia-grazing, not infochemicals, mediate a shift towards large inedible colonies of the gelatinous green alga Sphaerocystis schroeteri

The influence of Daphnia galeata×hyalina grazing and of infochemicals released by the daphnids on the colony size and growth rate of the colonial gelatinous green alga Sphaerocystis schroeteri (Chlorococcales) was investigated in laboratory batch experiments run for 96 h. High zooplankton grazing pressure was exerted by a final concentration of 100 daphnids L⁻¹ in the Daphnia treatments. Infochemicals were obtained by filtration (0.2 μm) of water from D. galeata×hyalina cultures (200 ind. L⁻¹ exposed for 24 h). This filtrate was added to the S. schroeteri cultures in two concentrations corresponding to 7 and 50 daphnids L⁻¹, respectively. The growth rate of S. schroeteri was neither affected significantly by direct Daphnia grazing nor by the presence of Daphnia infochemicals, in comparison to the control. However, the portion of inedible S. schroeteri colonies (diameter>50 μm) increased under direct grazing pressure, whereas the Daphnia infochemicals did not influence the colony size significantly. We conclude that the shift in colony size by direct zooplankton grazing denotes an effective defence mechanism against size selective feeding for colonial gelatinous green algae. This effective defence in combination with unchanged growth rates of the larger colonies (under non-limiting nutrient and light conditions) falsifies the assumption of a trade-off between minimising grazing losses and maximising growth by optimising the colony size.     

新疆布伦托海轮虫群落结构及其与环境因子的关系

本文于2006年秋至2008年夏,对新疆布伦托海轮虫的群落组成、时空分布及其与环境因子的相关性进行了系统的调查.结果表明,本次调查采集到布伦托海轮虫30种,主要优势种有针簇多枝轮虫(Polyarthra trigla)、长三肢轮虫(Filina longiseta)、奇异六腕轮虫(Hexarthra mira)、方形臂...     

Blooms of the toxic dinoflagellate Alexandrium fundyense in the western Gulf of Maine in 1993 and 1994: A comparative modeling study

Blooms of the toxic dinoflagellate Alexandrium fundyense commonly occur in the western Gulf of Maine but the amount of toxin observed in coastal shellfish is highly variable. In this study, a coupled physical–biological model is used to investigate the dynamics underlying the observed A. fundyense abundance and shellfish toxicity in 1993 (a high toxicity year) and 1994 (low toxicity year). The physical model simulates the spring circulation, while the biological model estimates the germination and population dynamics of A. fundyense based on laboratory and field data. The model captures the large-scale aspects of the initiation and development of A. fundyense blooms during both years, but small-scale patchiness and the dynamics of bloom termination remain problematic. In both cases, the germination of resting cysts accounts for the magnitude of A. fundyense populations early in the spring. Simulations with low net A. fundyense growth rates capture the mean observed concentration during the bloom peak, which is of similar magnitude during both years. There is little evidence that large-scale changes in biological dynamics between 1993 and 1994 were a primary driver of the differences in shellfish toxicity. Results instead suggest that the persistent southwesterly flow of the western Maine Coastal Current led to A. fundyense populations of similar alongshore extent by late May of both years. This period coincides with peak cell abundance in the region. Variations in wind forcing (downwelling favorable in 1993, upwelling favorable in 1994) and subsequent cell transport (inshore in 1993, offshore in 1994) in early June then provides a plausible explanation for the dramatic mid-June differences in shellfish toxicity throughout the western Gulf of Maine.     

Comparison of life-history parameters and effectiveness in competition with Moina macrocopa between two Keratella tropica morphs

Competition and predation are the two major biotic interactions in zooplankton community. Various rotifer species exhibit polymorphism in lorica and spine morphology in the presence of predators and competitors, which effectively increases their defence competition capacity. This study compared life-history parameters between long-spined morph (LSM) and short-spined morph (SSM) of Keratella tropica at three temperatures. The results showed that both morphs showed significantly higher intrinsic rate of population increase at 25 °C than those at 20 °C and 30 °C. The intrinsic rate of population increase was significantly lower in LSM than SSM at 25 °C and 30 °C, and the net reproductive rate was significantly lower in LSM than SSM at 30 °C. These results suggested that K. tropica had cost on development of long spines, which, however, was dependent on temperature. Moreover, the competition outcome of LSM and SSM with Moina macrocopa were investigated at three densities of Scenedesmus obliquus (0.5 × 10⁶, 1.0 × 10⁶ and 2.0 × 10⁶ cells/mL) and three initial densities of M. macrocopa (0.1, 0.3 and 0.9 ind./mL), revealing that both LSM and SSM populations decreased significantly at all tested algal densities and Moina levels. LSM showed higher population growth rate than SSM at lower algal density (0.5 × 10⁶ cells/mL) and higher Moina densities (0.3 and 0.9 ind./mL), which suggested that LSM could be more effective to defend M. macrocopa. Three-way analysis of variance showed no significant effects of algal density on population growth rate of K. tropica, but interaction between algal density and M. macrocopa density significantly affected its population growth rate.     

Relationship between differential retention of Escherichia coli and Enterococcus faecalis and variations in enzyme activity in the scallop Patinopecten yessoensis

Uptake of Escherichia coli and Enterococcus faecalis and variations of trypsin amylase activity acid phosphatase and superoxide dismutase in tissue of the scallop Patinopecten yessoensis were detected. The results showed that P. yessoensis accumulated E. faecalis in larger numbers and more rapidly than E. coli, both with the highest concentration in the digestive tract and lowest in hemolymph. Compared to E. coli, all scallops exposed to E. faecalis showed significantly higher trypsin and AMS activity. SOD activity in hemocytes and ACP activity in hemolymph was significantly higher in the treatments with 5 log₁₀CFU/ml E. coli than with E. faecalis. But no significant differences in ACP activity of P. yessoensis exposed to a 3 log₁₀CFU/ml inoculum of both bacteria were recorded. In conclusion, the mass retention of gut microflora in P. yessoensis is positively correlated with digestive enzymes activity and negatively correlated with ACP activity in the hemocyte.     

PCB exposure in sea otters and harlequin ducks in relation to history of contamination by the Exxon Valdez oil spill

Exposure to contaminants other than petroleum hydrocarbons could confound interpretation of Exxon Valdez oil spill effects on biota at Prince William Sound, Alaska. Hence, we investigated polychlorinated biphenyls (PCBs) in blood of sea otters and harlequin ducks sampled during 1998. PCB concentrations characterized by lower chlorinated congeners were highest in sea otters from the unoiled area, whereas concentrations were similar among harlequin ducks from the oiled and unoiled area. Blood enzymes often elevated by xenobiotics were not related to PCB concentrations in sea otters. Only sea otters from the unoiled area had estimated risk from PCBs, and PCB composition or concentrations did not correspond to reported lower measures of population performance in sea otters or harlequin ducks from the oiled area. PCBs probably did not influence limited sea otter or harlequin duck recovery in the oiled area a decade after the spill.     

The Determination of Relationship between Zooplankton and Abiotic Factors Using Canonical Correspondence Analysis (CCA) in the Ova Stream (Ankara/Turkey)

In this study totally 25 rotifer species, 3 cladocerus species, and 2 copepoda species were identified in Ova Stream. Among zooplankton rotifer species were dominant. Population densities of planktonic organisms were calculated as individual per cubic meter and the relationship of the planktonic organisms with physicochemical parameters was determined applying canonical correspondence analysis (CCA). Also diversity index has been calculated according to the sampling stations. According to the results, the diversity index changes between stations and sampling time. The CCA results show that the rotifer species Keratella, Notholca showed negative correlation with the increasing chemical parameters and temperature but Brachionus, Mytilina, Colurella, and Testudinella have positive correlation with increasing temperature.     

转轮虫(Rotaria rotatoria)生活史策略对温度和食物密度变化的响应

温度和食物密度一直是浮游动物生长和繁殖过程中的重要影响因子,而有关两者对蛭态类轮虫生活史策略的影响研究鲜有报道.本研究基于生命表,探究了温度（24、28和32℃）和食物密度（0.5×10⁶、1.0×10⁶、2.0×10⁶和4.0×10⁶ cells/ml）对蛭态类转轮虫（Rotaria rotatoria）生长、存活和繁殖的影响.结果表明,温度对转轮虫的所有生活史参数均具有极显著影响,食物密度以及温度和食物密度的交互作用显著地影响轮虫的主要发育阶段历时和生命期望.无论食物密度如何,随着温度的升高,轮虫的发育速度加快,繁殖期缩短,寿命减小,特定年龄存活率开始下降的时间逐渐前移,下降的速度也逐渐加快,但内禀增长率增加,这是转轮虫在不利条件下的生活史策略——能量分配权衡.在较低温度下（24和28℃）,转轮虫的各主要发育阶段历时随着食物密度的增加而减小,但在高温下（32℃）高食物密度却对轮虫有利,轮虫的繁殖期和寿命随食物密度的增加而增加;可见,在适宜培养温度下,较低的藻浓度有利于轮虫的生长和繁殖,而高温下则需要高食物密度以保证轮虫对物质和能量的需求.在较低食物密度下,温度对转轮虫生活史参数的影响更加显著,这表明食物密度是重要的限制因子,其在量上的限制或不足会增加轮虫对基础代谢生理调节的代价,从而对其他因子的耐受能力可能随之缩小.     

Modelling the effects of temperature and food on individual growth and reproduction of Daphnia and their consequences on the population level

Individual based simulations of population dynamics require the availability of growth models with adequate complexity. For this purpose a simple-to-use model (non-linear multiple regression approach) is presented describing somatic growth and reproduction of Daphnia as a function of time, temperature and food quantity. The model showed a good agreement with published observations of somatic growth (r² = 0.954, n = 88) and egg production (r² = 0.898, n = 35). Temperature is the main determinant of initial somatic growth and food concentration is the main determinant of maximal body length and clutch size. An individual based simulation was used to demonstrate the simultaneous effects of food and temperature on the population level. Evidently, both temperature and food supply affected the population growth rate but at food concentrations above approximately 0.4 mg Cl⁻¹ Scenedesmus acutus temperature appeared as the main determinant of population growth.

Four simulation examples are given to show the wide applicability of the model: (1) analysis of the correlation between population birth rate and somatic growth rate, (2) contribution of egg development time and delayed somatic growth to temperature-effects on population growth, (3) comparison of population birth rate in simulations with constant vs. decreasing size at maturity with declining food concentrations and (4) costs of diel vertical migration. Due to its plausible behaviour over a broad range of temperature (2–20 °C) and food conditions (0.1–4 mg Cl⁻¹) the model can be used as a module for more detailed simulations of Daphnia population dynamics under realistic environmental conditions.     

Concentrations of PAHs (polycyclic aromatic hydrocarbons) in various marine organisms in relation to those in sediments and to trophic level

Sediments and marine organisms such as mussels (Mytilus galloprovincialis), fish (Mullus barbatus, Serranus scriba), crabs (Polybius henslowi) and shrimp-like crustaceans (mysids, euphausiids) were sampled on European coasts, France, Spain) and analysed for their PAH content by gas chromatography/mass spectrometry (GC-MS). The PAH concentrations determined in the organisms were significantly negatively correlated with logK_ow values, but to different extents according to the species. This was explained by the greater exposure of the organisms to the lower molecular weight compounds with differences resulting from different feeding habits, habitat, and biotransformation capacities of the organisms in relation to trophic levels. The influence of these observations on the use of different sentinel species in biomonitoring programmes is discussed.     

Relationships between surface-active organic substances, chlorophyll a and nutrients in the northern Adriatic Sea

Surface-active organic substances (SAS), nutrients, chlorophyll a (Chl a), dissolved oxygen (O₂), salinity (S) and temperature (T) were measured approximately monthly in the northern Adriatic Sea (NA) during two years. Results were elaborated for two stations of different trophic status. Exhaustive statistical examinations of measured variables were performed to contribute to better understanding of the processes and interdependence of the measured parameters. The results of those analyses allowed the region to be described based on several groups of data relevant for particular processes. Chl a appear to be independent parameter; orthosilicate, nitrite and ammonium represent the parameter group most relevant for regeneration processes; the second group includes nutrients dominant in freshwater inputs (orthophosphate and nitrate), while in the third group are parameters related to primary production processes (T, O₂, organic phosphorus and SAS). Probably due to time lags between nutrient uptake and phytoplankton growth as well as between cells growth/division and SAS release, both, the correlation between nutrients and Chl a and the correlation between Chl a and SAS were not statistically significant. Although significantly higher nutrient and Chl a concentrations were found in the mesotrophic western part of the region in comparison to the oligotrophic eastern part, the SAS concentration differences were less marked. It is assumed that more OM is produced, but also at the same time remineralized, at the mesotrophic western part, leading to higher concentrations of regenerated nutrients there. Presence of different phytoplankton taxa and/or fractions (micro and nano) at the two stations may have also contributed to minimize the difference in SAS contents. Higher SAS acidities were noticed in the mesotrophic than in the oligotrophic part.     

A model that helps explain Sr-isotope disequilibrium between feldspar phenocrysts and melt in large-volume silicic magma systems

Feldspar phenocrysts of silicic volcanic rocks are commonly in Sr-isotopic disequilibrium with groundmass. In some cases the feldspar is more radiogenic, and in others it is less radiogenic. Several explanations have been published previously, but none of these is able to accommodate both senses of disequilibrium. We present a model by which either more- or less-radiogenic feldspar (or even both within a single eruptive unit) can originate. The model requires a magma body open to interaction with biotite- and feldspar-bearing wall rock. Magma is incrementally contaminated as wall rock melts incongruently. Biotite preferentially melts first, followed by feldspar. Such melting behavior, which is supported by both field and experimental studies, first contaminates magma with a relatively radiogenic addition, followed by a less-radiogenic addition. Feldspar phenocrysts lag behind melt (groundmass of volcanic rock) in incorporating the influx of contaminant, thus resulting in Sr-isotopic disequilibrium between the crystals and melt. The sense of disequilibrium recorded in a volcanic rock depends on when eruption quenches the contamination process. This model is testable by isotopic fingerprinting of individual feldspar crystals. For a given set of geologic boundary conditions, specific core-to-rim Sr-isotopic profiles are expectable. Moreover, phenocrysts that nucleate at different times during the contamination process should record different and predictable parts of the history. Initial results of Sr-isotopic fingerprinting of sanidine phenocrysts from the Taylor Creek Rhyolite are consistent with the model. More tests of the model are desirable.     

Potential link between exposure to fungicides chlorothalonil and mancozeb and haemic neoplasia development in the soft-shell clam Mya arenaria: A laboratory experiment

The aetiology of haemic neoplasia (HN) is unknown, so far but many causative factors are suggested such as viral, pollution and genetics. The aim of this study was to determine if, under chronic exposure, two major pesticides (chlorothalonil and mancozeb) which are used in potato production could induce HN in soft-shell clams (Mya arenaria). Short-term experiments with acute exposure were also performed. Clams were collected from an epizootic site (North River, PEI) and from a site free of the disease (Magdalen Islands, Quebec). The tetraploid level of haemocytes was assessed by flow cytometry for each clam to determine the HN status. The bioaccumulation of pesticides in tissues was quantified by gas chromatography/mass spectrometry (GC/MS) for chlorothalonil while mancozeb and manganese were quantified by inductively coupled plasma-mass spectrometer (ICP/MS). Long term exposure to fungicide Bravo 500^® did not induce high tetraploid levels on negative calm from North River and the analysis of the digestive gland and the mantle did not reveal any detectable level of chlorothalonil. In the Manzate 200 DF^®, some clams revealed high level of tetraploid cells but no difference were observed between the treatments and the control. The analysis of the digestive gland and the mantle for manganese did not highlight any significant difference in tissue concentration (p = 0.05). For the acute exposure, chlorothalonil analysis showed that the active ingredient is distributed between four chlorinated compounds: 99.5% for chlorothalonil isomers, 0.4% for pentachlorothalonil and 0.1% for trichlorothalonil isomers. For a 72 h experiment, the accumulation was within 4 h; the higher tissue concentration of chlorothalonil was 59.2 μg g⁻¹ in the mantle after 48 h, following by a decrease to an undetectable level at the end. For the manganese, the accumulation was detected after 4 h; the higher tissue concentration was 48.8 μg g⁻¹ in the mantle after 24 h and, over the following 48 h, the accumulation decreased until the end of the trial. Based on the data, the accumulation of these fungicides seems to be transitory. Chlorothalonil and mancozeb are both oxidative-stress promoters and could have induced cell dysfunction while in the tissue. Study on the effect of these fungicides on the p53 protein system is an example of strategy that would provide information on cellular events promoting neoplasia.     

The relationship between the growth process of the ferromanganese crusts in the pacific seamount and Cenozoic ocean evolvement

Base on the Os isotope stratigraphy together with the empirical growth rate models using Co concentrations, the growth ages of the ferromanganese crusts MHD79 and MP3D10 distributed in the seamount of Pacific are confirmed. Through the contrast and research on the previous achievements including ODP Leg 144 and the crusts CD29-2, N5E-06 and N1–15 of the seamount of the Central Pacific, the uniform five growth and growth hiatus periods of them are found, and closely related to the Cenozoic ocean evolvement process. In the Paleocene Carbon Isotope Maximum (PCIM), the rise of the global ocean productivity promoted the growth of the seamount crust; the first growth hiatus (I) of the ferromanganese crust finished. In the Paleocene-Eocene Thermal Maximum (PETM), though the vertical exchange of seawater was weakened, the strong terrestrial chemical weathering led to the input of a great amount of the terrigenous nutrients, which made the bioproductivity rise, so there were no crust hiatuses. During 52–50 Ma, the Early Eocene Optimum Climate (EECO), the two poles were warm, the latitudinal temperature gradient was small, the wind-driven sea circulation and upwelling activity were weak, the terrestrial weathering was also weakened, the open ocean bioproductivity decreased, and the ferromanganese crust had growth hiatus again (II). From early Middle Eocene-Late Eocene, Oligocene, it was a long-term gradually cooling process, the strengthening of the sea circulation and upwelling led to a rise of bioproductivity, and increase of the content of the hydrogenous element Fe, Mn and Co and the biogenous element Cu, Zn, so that was the most favorable stage for the growth of ferromanganese crust (growth periods III and IV) in the studied area. The hiatus III corresponded with the Eocene- Oligocene boundary, is inferred to relate with the global climate transformation, celestial body impact event in the Eocene-Oligocene transition. From the early to the middle Miocene, a large-scale growth hiatus (hiatus period IV) of the ferromanganese crust in the studied area is inferred to relate with temporary warm up climate and ephemeral withdrawal of Antarctic bottom water in the early Miocene. After that, the Antarctic ice sheets extended, the bottom water circumfluence strengthened, the ocean fertility increased, and the once interrupted crust continued to grow in the late Miocene (growth period V). Supported by China Ocean Mineral Resources Research and Development Association “10th Five Year” Topic (Grant No. DY105-01-04-14)     

Climatic features of atmospheric heat source/sink over the Qinghai-Xizang Plateau in 35 years and its relation to rainfall in China

Using the 1961–1995 monthly averaged meteorological data from 148 surface stations in the Qinghai-Xizang Plateau (QXP) and its surrounding areas, calculation of the 35-year atmospheric heat source/sink (<Qi>) and an analysis on its climatic features and relation to rainfall in China have been made. It is found that on the average, the atmospheric heat source over the QXP is the strongest in June (78 W / m²) and cold source is the strongest in December (−72 W/m²). The sensible heat of the surface increases remarkably over the southwest of the QXP, causing the obvious increase of <Qi> there in February and March, which makes a center of the atmospheric heat source appear over the north slope of the Himalayas. Afterwards, this center continues to intensify and experiences noticeable migration westwards twice, separately occurring in April and June. The time when the atmosphere over the east of the QXP becomes heat source and reaches strongest is one month later than that over the southwest of the QXP. In summer, the latent heat of condensation becomes a heating factor as important as the sensible heat and is also a main factor that makes the atmospheric heat source over the east of the QXP continue growing. On the interdecadal time scale, (Q₁) of the QXP shows an abrupt change in 1977 and a remarkable increase after 1977. The atmospheric heat source of the spring over the QXP is a good indicator for the subsequent summer rainfall over the valleys of the Changjiang and Huaihe rivers and South China and North China. There is remarkable positive correlation between the QXP heat source of summer and the summer rainfall in the valleys of the Changjiang River.     

胜利油田油水井动态与华北地震关系的初步分析

近些年来,在收集总结大量复杂多样的地震前兆资料中,油水井动态异常以其变化量显著而引人注目。本文以近几年华北地区所发生的地震为实例,旨在分析研究胜利油田油水井的产油量、产液量、油水井压力(地层压力)在震前震后的变化规律。     

Sensitivity of soil parameters in unsaturated zone modelling and the relation between effective,laboratory and in situ estimates

Simulation of soil moisture content requires effective soil hydraulic parameters that are valid at the modelling scale. This study investigates how these parameters can be estimated by inverse modelling using soil moisture measurements at 25 locations at three different depths (at the surface, at 30 and 60 cm depth) on an 80 by 20 m hillslope. The study presents two global sensitivity analyses to investigate the sensitivity in simulated soil moisture content of the different hydraulic parameters used in a one‐dimensional unsaturated zone model based on Richards' equation. For estimation of the effective parameters the shuffled complex evolution algorithm is applied. These estimated parameters are compared to their measured laboratory and in situ equivalents. Soil hydraulic functions were estimated in the laboratory on 100 cm³ undisturbed soil cores collected at 115 locations situated in two horizons in three profile pits along the hillslope. Furthermore, in situ field saturated hydraulic conductivity was estimated at 120 locations using single‐ring pressure infiltrometer measurements. The sensitivity analysis of 13 soil physical parameters (saturated hydraulic conductivity (K_s), saturated moisture content (θ_s), residual moisture content (θ_r), inverse of the air‐entry value (α), van Genuchten shape parameter (n), Averjanov shape parameter (N) for both horizons, and depth (d) from surface to B horizon) in a two‐layer single column model showed that the parameter N is the least sensitive parameter. K_s of both horizons, θ_s of the A horizon and d were found to be the most sensitive parameters. Distributions over all locations of the effective parameters and the distributions of the estimated soil physical parameters from the undisturbed soil samples and the single‐ring pressure infiltrometer estimates were found significantly different at a 5% level for all parameters except for α of the A horizon and K_s and θ_s of the B horizon. Different reasons are discussed to explain these large differences. Copyright © 2004 John Wiley & Sons, Ltd.     

Comparison of hydraulic geometry between sand‐ and gravel‐bed rivers in relation to channel pattern discrimination

A comparison has been made between the hydraulic geometry of sand‐ and gravel‐bed rivers, based on data from alluvial rivers around the world. The results indicate a signi?cant difference in hydraulic geometry among sand‐ and gravel‐bed rivers with different channel patterns. On this basis, some diagrams for discrimination of meandering and braided channel patterns have been established. The relationships between channel width and water discharge, between channel depth and water discharge, between width–depth ratio and water discharge and between channel slope and water discharge can all be used for channel pattern discrimination. The relationship between channel width and channel depth can also be used for channel pattern discrimination. However, the accuracy of these relationships for channel pattern discrimination varies, and the depth–discharge relationship is a better discriminator of pattern type than the classic slope–discharge function. The cause for this difference has been explained qualitatively. To predict the development of channel patterns under different natural conditions, the pattern discriminator should be searched on the basis of independent or at least semi‐independent variables. The relationship between stream power and bed material grain size can be used to discriminate channel patterns, which shows a better result than the discriminator using the slope–discharge relationship. Copyright © 2004 John Wiley & Sons, Ltd.     

Factors contributing to heavy metal accumulation in sediments and in the intertidal mussel Perna perna in the Gulf of Annaba (Algeria)

This paper presents the results of a seasonal survey of heavy metals accumulated in sediments and in the soft parts of the body of the mussel Perna perna at four stations in the Gulf of Annaba (Algeria). Pooled soft tissues from 10 mussels representing the entire range of sizes were digested in nitric acid. Statistical analysis reveals a significant seasonal effect on all the measured metals, the highest values being recorded in winter. With the exception of Cr, the levels for all metals were significantly higher in the east, at the outlet of the Seybouse River, than at all other monitoring stations. The study also shows that north-western waters are subject to a significantly lower degree of heavy metal pollution than elsewhere in the gulf. Levels were nevertheless within the limits of public health standards.