Ecological-niche factor analysis: how to compute habitat-suitability maps without absence data?

We propose a multivariate approach to the study of geographic species distribution which does not require absence data. Building on Hutchinson's concept of the ecological niche, this factor analysis compares, in the multidimensional space of ecological variables, the distribution of the localities where the focal species was observed to a reference set describing the whole study area. The first factor extracted maximizes the marginality of the focal species, defined as the ecological distance between the species optimum and the mean habitat within the reference area. The other factors maximize the specialization of this focal species, defined as the ratio of the ecological variance in mean habitat to that observed for the focal species. Eigenvectors and eigenvalues are readily interpreted and can be used to build habitat-suitability maps. This approach is recommended in Situations where absence data are not available (many data banks), unreliable (most cryptic or rare species), or meaningless (invaders). We provide an illustration and validation of the method for the alpine ibex, a species reintroduced in Switzerland which presumably has not yet recolonized its entire range.

A multi-temporal approach to model endangered species distribution in Europe: the case of the Eurasian otter in Italy

We analysed the relationship between changes in land cover patterns and the Eurasian otter occurrence over the course of about 20 years (1985-2006) using multi-temporal Species Distribution Models (SDMs). The study area includes five river catchments covering most of the otter's Italian range. Land cover and topographic data were used as proxies of the ecological requirements of the otter within a 300-m buffer around river courses. We used species presence, pseudo-absence data, and environmental predictors to build past (1985) and current (2006) SDMs by applying an ensemble procedure through the BIOMOD modelling package. The performance of each model was evaluated by measuring the area under the curve (AUC) of the receiver-operating characteristic (ROC). Multi-temporal analyses of species distribution and land cover maps were performed by comparing the maps produced for 1985 and 2006. The ensemble procedure provided a good overall modelling accuracy, revealing that elevation and slope affected the otter's distribution in the past; in contrast, land cover predictors, such as cultivations and forests, were more important in the present period. During the transition period, 20.5% of the area became suitable, with 76% of the new otter presence data being located in these newly available areas. The multi-temporal analysis suggested that the quality of otter habitat improved in the last 20 years owing to the expansion of forests and to the reduction of cultivated fields in riparian belts. The evidence presented here stresses the great potential of riverine habitat restoration and environmental management for the future expansion of the otter in Italy

Effect of spatial processes on the evolution of ecological niches and dispersal in metacommunity systems

RésuméLa coexistence de nombreuses espèces différentes a de tout temps intrigué les biologistes. La diversité et la composition des communautés sont influencées par les perturbations et l'hétérogénéité des conditions environnementales. Bien que dans la nature la distribution spatiale des conditions environnementales soit généralement autocorrélée, cet aspect est rarement pris en compte dans les modèles étudiant la coexistence des espèces. Dans ce travail, nous avons donc abordé, à l'aide de simulations numériques, la coexistence des espèces ainsi que leurs caractéristiques au sein d'un environnement autocorrélé.Afin de prendre en compte cet élément spatial, nous avons développé un modèle de métacommunauté (un ensemble de communautés reliées par la dispersion des espèces) spatialement explicite. Dans ce modèle, les espèces sont en compétition les unes avec les autres pour s'établir dans un nombre de places limité, dans un environnement hétérogène. Les espèces sont caractérisées par six traits: optimum de niche, largeur de niche, capacité de dispersion, compétitivité, investissement dans la reproduction et taux de survie. Nous nous sommes particulièrement intéressés à l'influence de l'autocorrélation spatiale et des perturbations sur la diversité des espèces et sur les traits favorisés dans la métacommunauté. Nous avons montré que l'autocorrélation spatiale peut avoir des effets antagonistes sur la diversité, en fonction du taux de perturbations considéré. L'influence de l'autocorrélation spatiale sur la capacité de dispersion moyenne dans la métacommunauté dépend également des taux de perturbations et survie. Nos résultats ont aussi révélé que de nombreuses espèces avec différents degrés de spécialisation (i.e. différentes largeurs de niche) peuvent coexister. Toutefois, les espèces spécialistes sont favorisées en absence de perturbations et quand la dispersion est illimitée. A l'opposé, un taux élevé de perturbations sélectionne des espèces plus généralistes, associées avec une faible compétitivité.L'autocorrélation spatiale de l'environnement, en interaction avec l'intensité des perturbations, influence donc de manière considérable la coexistence ainsi que les caractéristiques des espèces. Ces caractéristiques sont à leur tour souvent impliquées dans d'importants processus, comme le fonctionnement des écosystèmes, la capacité des espèces à réagir aux invasions, à la fragmentation de l'habitat ou aux changements climatiques. Ce travail a permis une meilleure compréhension des mécanismes responsables de la coexistence et des caractéristiques des espèces, ce qui est crucial afin de prédire le devenir des communautés naturelles dans un environnement changeant.AbstractUnderstanding how so many different species can coexist in nature is a fundamental and long-standing question in ecology. Community diversity and composition are known to be influenced by heterogeneity in environmental conditions and disturbance. Though in nature the spatial distribution of environmental conditions is frequently autocorrelated, this aspect is seldom considered in models investigating species coexistence. In this work, we thus addressed several questions pertaining to species coexistence and composition in spatially autocorrelated environments, with a numerical simulations approach.To take into account this spatial aspect, we developed a spatially explicit model of metacommunity (a set of communities linked by dispersal of species). In this model, species are trophically equivalent, and compete for space in a heterogeneous environment. Species are characterized by six life-history traits: niche optimum, niche breadth, dispersal, competitiveness, reproductive investment and survival rate. We were particularly interested in the influence of environmental spatial autocorrelation and disturbance on species diversity and on the traits of the species favoured in the metacommunity. We showed that spatial autocorrelation can have antagonistic effects on diversity depending on disturbance rate. Similarly, spatial autocorrelation interacted with disturbance rate and survival rate to shape the mean dispersal ability observed in the metacommunity. Our results also revealed that many species with various degrees of specialization (i.e. different niche breadths) can coexist together. However specialist species were favoured in the absence of disturbance, and when dispersal was unlimited. In contrast, high disturbance rate selected for more generalist species, associated with low competitive ability.The spatial structure of the environment, together with disturbance and species traits, thus strongly impacts species diversity and, more importantly, species composition. Species composition is known to affect several important metacommunity properties such as ecosystem functioning, resistance and reaction to invasion, to habitat fragmentation and to climate changes. This work allowed a better understanding of the mechanisms responsible for species composition, which is of crucial importance to predict the fate of natural metacommunities in changing environments

Fire management, climate change and their interacting effects on birds in complex Mediterranean landscapes: dynamic distribution modelling of an early-successional species - the near-threatened Dartford Warbler (Sylvia undata)

The current challenge in a context of major environmental changes is to anticipate the responses of species to future landscape and climate scenarios. In the Mediterranean basin, climate change is one the most powerful driving forces of fire dynamics, with fire frequency and impact having markedly increased in recent years. Species distribution modelling plays a fundamental role in this challenge, but better integration of available ecological knowledge is needed to adequately guide conservation efforts. Here, we quantified changes in habitat suitability of an early-succession bird in Catalonia, the Dartford Warbler (Sylvia undata) ― globally evaluated as Near Threatened in the IUCN Red List. We assessed potential changes in species distributions between 2000 and 2050 under different fire management and climate change scenarios and described landscape dynamics using a spatially-explicit fire-succession model that simulates fire impacts in the landscape and post-fire regeneration (MEDFIRE model). Dartford Warbler occurrence data were acquired at two different spatial scales from: 1) the Atlas of European Breeding Birds (EBCC) and 2) Catalan Breeding Bird Atlas (CBBA). Habitat suitability was modelled using five widely-used modelling techniques in an ensemble forecasting framework. Our results indicated considerable habitat suitability losses (ranging between 47% and 57% in baseline scenarios), which were modulated to a large extent by fire regime changes derived from fire management policies and climate changes. Such result highlighted the need for taking the spatial interaction between climate changes, fire-mediated landscape dynamics and fire management policies into account for coherently anticipating habitat suitability changes of early succession bird species. We conclude that fire management programs need to be integrated into conservation plans to effectively preserve sparsely forested and early succession habitats and their associated species in the face of global environmental change.

Using species richness and functional traits predictions to constrain assemblage predictions from stacked species distribution models

Aim: Modelling species at the assemblage level is required to make effective forecast of global change impacts on diversity and ecosystem functioning. Community predictions may be achieved using macroecological properties of communities (MEM), or by stacking of individual species distribution models (S-SDMs). To obtain more realistic predictions of species assemblages, the SESAM framework suggests applying successive filters to the initial species source pool, by combining different modelling approaches and rules. Here we provide a first test of this framework in mountain grassland communities. Location: The western Swiss Alps. Methods: Two implementations of the SESAM framework were tested: a "Probability ranking" rule based on species richness predictions and rough probabilities from SDMs, and a "Trait range" rule that uses the predicted upper and lower bound of community-level distribution of three different functional traits (vegetative height, specific leaf area and seed mass) to constraint a pool of environmentally filtered species from binary SDMs predictions. Results: We showed that all independent constraints expectedly contributed to reduce species richness overprediction. Only the "Probability ranking" rule allowed slightly but significantly improving predictions of community composition. Main conclusion: We tested various ways to implement the SESAM framework by integrating macroecological constraints into S-SDM predictions, and report one that is able to improve compositional predictions. We discuss possible improvements, such as further improving the causality and precision of environmental predictors, using other assembly rules and testing other types of ecological or functional constraints.

Impact of model complexity on cross-temporal transferability in Maxent species distribution models: An assessment using paleobotanical data

Maximum entropy modeling (Maxent) is a widely used algorithm for predicting species distributions across space and time. Properly assessing the uncertainty in such predictions is non-trivial and requires validation with independent datasets. Notably, model complexity (number of model parameters) remains a major concern in relation to overfitting and, hence, transferability of Maxent models. An emerging approach is to validate the cross-temporal transferability of model predictions using paleoecological data. In this study, we assess the effect of model complexity on the performance of Maxent projections across time using two European plant species (Alnus giutinosa (L.) Gaertn. and Corylus avellana L) with an extensive late Quaternary fossil record in Spain as a study case. We fit 110 models with different levels of complexity under present time and tested model performance using AUC (area under the receiver operating characteristic curve) and AlCc (corrected Akaike Information Criterion) through the standard procedure of randomly partitioning current occurrence data. We then compared these results to an independent validation by projecting the models to mid-Holocene (6000 years before present) climatic conditions in Spain to assess their ability to predict fossil pollen presence-absence and abundance. We find that calibrating Maxent models with default settings result in the generation of overly complex models. While model performance increased with model complexity when predicting current distributions, it was higher with intermediate complexity when predicting mid-Holocene distributions. Hence, models of intermediate complexity resulted in the best trade-off to predict species distributions across time. Reliable temporal model transferability is especially relevant for forecasting species distributions under future climate change. Consequently, species-specific model tuning should be used to find the best modeling settings to control for complexity, notably with paleoecological data to independently validate model projections. For cross-temporal projections of species distributions for which paleoecological data is not available, models of intermediate complexity should be selected.

Climate niche conservatism does not explain restricted distribution patterns in Tynanthus (Bignonieae, Bignoniaceae)

Studies on niche evolution allow us to establish how species niches have changed over time as well as to identify how long-term evolutionary processes have led to present-day species distributions. Here, we investigate the patterns of climatic niche evolution in Tynanthus (Bignonieae, Bignoniaceae), a genus comprising narrowly distributed species. We test the hypothesis that niche conservatism has played an important role in the diversification history of this group of Neotropical lianas. For that, we perform univariate and multivariate comparisons between species' climatic niches and associated environmental data with information on species' phylogenetic relationships. We encountered considerable divergence in niches among species, indicating that niche conservatism in climatic variables has does not seem to havenot played a key role in the diversification of the genus. Our results are used as a basis to discuss patterns of ecological niche evolution in the group and to suggest novel approaches for future analyses.

Social differences in avoidable mortality between small areas of 15 European cities: an ecological study

Health and inequalities in health among inhabitants of European cities are of major importance for European public health and there is great interest in how different health care systems in Europe perform in the reduction of health inequalities. However, evidence on the spatial distribution of cause-specific mortality across neighbourhoods of European cities is scarce. This study presents maps of avoidable mortality in European cities and analyses differences in avoidable mortality between neighbourhoods with different levels of deprivation. Methods: We determined the level of mortality from 14 avoidable causes of death for each neighbourhood of 15 large cities in different European regions. To address the problems associated with Standardised Mortality Ratios for small areas we smooth them using the Bayesian model proposed by Besag, York and Mollié. Ecological regression analysis was used to assess the association between social deprivation and mortality. Results: Mortality from avoidable causes of death is higher in deprived neighbourhoods and mortality rate ratios between areas with different levels of deprivation differ between gender and cities. In most cases rate ratios are lower among women. While Eastern and Southern European cities show higher levels of avoidable mortality, the association of mortality with social deprivation tends to be higher in Northern and lower in Southern Europe. Conclusions: There are marked differences in the level of avoidable mortality between neighbourhoods of European cities and the level of avoidable mortality is associated with social deprivation. There is no systematic difference in the magnitude of this association between European cities or regions. Spatial patterns of avoidable mortality across small city areas can point to possible local problems and specific strategies to reduce health inequality which is important for the development of urban areas and the well-being of their inhabitants

Diameter structure of stratum and ecological groups of a Semideciduous Forest area in Dionísio-MG

The objective of this study was to obtain homogeneous groups of species and information on their density, dominance and volume, in terms of ecological group and diameter structure of an area of Submontane Semideciduous forest (Mata do Mumbaça) in Dionísio, MG. This work was conducted with data of the diameter distribution per species from floristic and phytosociological (Mata do Mumbaça) survey of 120 plots with 10 x 10 m each one. The 120 plots were contiguous and corresponding to a total sample area of 12,000 m² distributed over the topographic units (Low Ramp, Lower Slope, Upper Slope and Hill Top). The topographic units Low Ramp, Lower Slope and Upper Slope were in the middle stage of succession as they presented incipient stratification into two strata (canopy and understory) i.e. canopy ranging from 5 to 12 m high. However, the stratum Hill Top was classified as intermediate/advanced succession because it had a total height equal to or greater than 12 m. The distribution of individual trees of the four strata on diameter classes showed a typical J-inverted pattern that is, high concentration of individuals in smaller diameter classes and a sharp reduction towards the larger classes. In relation to absolute dominance and total volume of species, the ecological group that stood out in the four strata (Low Ramp, Lower Slope, Upper Slope and Hill Top) was the initial secondary, which were in the intermediate stage of secondary, rapidly developing into the mature phase.

Influence of seed size and ecological factors on the germination and emergence of field bindweed (Convolvulus arvensis)

An understanding of seed germination ecology of weeds can assist in predicting their potential distribution and developing effective management strategies. Influence of environmental factors and seed size on germination and seedling emergence of Convolvulus arvensis (field bindweed) was studied in laboratory and greenhouse conditions. Germination occurred over a wide range of constant temperatures, between 15 and 40 ºC, with optimum germination between 20 and 25 ºC. Time to start germination, time to 50% germination and mean germination time increased while germination percentage and germination index decreased with an increase in temperature from 20 ºC, salinity and osmotic stress. However, germination was tolerant to low salt (25 mM) or osmotic stress (0.2 MPa), but as salinity and osmotic stress increased, germination percentage and germination index decreased. Seeds of C. arvensis placed at soil surface showed maximum emergence and decreased as seeding depth increased. Seeds of C. arvensis germinated over a wide range of pH (4 to 9) but optimum germination occurred at pH 6 to 8. Under highly alkaline and acidic pH, time to start germination, time to 50% germination and mean germination time increased while germination percentage and germination index decreased. Increase in field capacity caused decreased time to start germination, time to 50% germination and mean germination time but increased germination percentage and germination index. Bigger seeds had low time to start germination, time to 50% germination and mean germination time but high germination percentage and germination index. Smaller seeds were more sensitive to environmental factors as compared to larger or medium seeds. It can be concluded that except for pH, all environmental factors and seed sizes adversely affect C. arvensis as regards seed germination or emergence and germination or emergence traits, and larger seeds result in improved stand establishment and faster germination than small seeds, regardless of moisture stress or deeper seeding depth.

STUDY OF AUTUMN-WINTER WEEDS IN PORTUGAL AND ITS ECOLOGICAL PREFERENCES

ABSTRACTThe conclusion of the dam project located in Alqueva, in Southern Portugal, has resulted in a significant increase of new irrigated areas, since 2006. This has meant that, in recent years, there have been progressive flora changes in farming systems traditionally implemented in the Alentejo region. The present work has analyzed the weed flora in an early stage of these changes, and the impact of environmental factors on the distribution of natural vegetation under Mediterranean climate conditions in the influence area of Alqueva. In 2007, 105 floristic surveys were carried out in autumn-winter crop plots or other soil use, and 264 species were identified. Families with higher expression were: Asteraceae, Poaceae, and Fabaceae. Only three species have been identified in more than half of farms, Avena sterilis, Phalaris minor and Lolium rigidum and they were part of the 15 species that revealed high and very high infestation degrees. Soil texture and extractable phosphorus have been determined as active ecological factors, according to the method of ecological profiles and Mutual Information. Therefore, these factores were those with the greatest influence on the species distribution. L.rigidum distribution showed to be associated with medium soil texture and A.sterilis distribution also showed to be associated with medium and fine soil texture soils, without showing ecological preference by extractable phosphorus. The distribution of P.minorwas not related to the soil texture but showed preference for soils with medium phosphorus content.

Taxonomy and distribution of the green algal genus Halimeda (Bryopsidales, Chlorophyta) in Brazil

Halimeda is a genus of calcified coenocytic green algae with a well known ecological importance in some tropical areas. Bleached calcified segments of Halimeda may accumulate in large deposits of economic potential as is the case in the northeastern coast of Brazil. In a survey of the genus in Brazil based on recent collections and examination of abundant material deposited on Brazilian herbaria we identified seven species: Halimeda cuneata Hering, H. discoidea Decaisne, H. gracilis Harvey ex J. Agardh, H. incrassata (Ellis) Lamouroux, H. opuntia (Linnaeus) Lamouroux, H. simulans Howe and H. tuna (Ellis & Solander) Lamouroux. These species are described in detail, with emphasis on diagnostic characters. Our study has shown that the shape and size of the utricula in surface view, under scanning electron microscopy, can be utilized to discriminate some species. Fertile specimens of Halimeda cuneata and H. discoidea are reported for the first time in the region. Data on vertical and geographical distribution are presented for each species and the southern limit of the genus in the western Atlantic was extended.

Patterns of spatial distribution in macroalgal communities from tropical lotic ecosystems

Three sampling sites were analysed in each of the following tropical regions: 1) northwestern São Paulo State, representing a disturbed region; 2) Bonito, Mato Grosso do Sul State, representing a hard water region; and 3) Ubatuba, northern costal region of São Paulo State, a well preserved tropical rainforest region. The hard water region had the highest mean values for macroalgal species richness (6.3) and diversity index (H' = 0.62). Northwest and rainforest regions had the highest percent cover values (22.5% and 17.0%, respectively). All sites in the northwest region had one or two dominant species (percent cover significantly higher than the remaining species), characterizing the niche pre-emption distribution pattern. The same pattern was found in two sites of the Atlantic rainforest. The hard water region had dominance of one species in two out of the three sites, but differently from the northwest region, niche overlap values were lower, evidencing a patch distribution. Competition for space was one of the main factors to explain spatial distribution. Overall, sites characterized by niche pre-emption had lower species richness, higher values for niche width and overlap, dominance index and percent cover of dominant species. In contrast, sites characterized by patch distribution had higher species richness and lower values for niche overlap and width, dominance index and percent cover.

Some ecological effects of beaver upon the watersheds in the Porcupine Hills, Alberta /

This ecological investigation of the beaver (Castor canadensis canadensis Kuhl) was part of the Eastern Slopes (Alberta) Watershed Research program conducted during the summers of 1968 . 1969 . 1970. and 1971 . A soil particle size analysis indicated that soil heavy with clay was used by the beaver for construction purposes in the Porcupine Hills . Examples were given of the beaver controlling erosion with the construction of dams . and also causing erosion to occur . However . in general . the beaver slow down soil erosion and decrease the loss of soil from this region. The beaver utilized measureable amounts of herbaceous vegetation . However, the utilization of herbs by the beaver requires further investigation. A system of ex~sures and enclosures of herbaceous vegetation plots was used to study the utilization of sedges , grasses , and forbs . The beaver indicated stronger species preference for willow as compared to aspen . The size preference for willow was in the 1.0 to 10.0 centimeters basal diameter classes , while the beaver utilization of aspen indicated a preference for the largest trees with basal diameter greater than 20.0 centimeters. Willow was the most important plant in the low lying areas with regard to distribution, abundance and to produce sustained yield. The beaver used this plant for food and construction throughout the study area. The distribution of aspen was limited. and this species did not appear to produce more than one crop in the lifetime of a beaver colony. Nine out of 15 woody plant types were sampled by the beaver in this region. A plot-intercept transect technique was used for systematic vegetation sampling of the woody vegetation in six intensively studied watersheds. The beaver population of the Porcupine Hills region of SW-Alberta is believed to depend upon the chinooks. During the 1971 ground census, a total of 60 active beaver colonies were tallied on the 930 square kilometers large study area. The beaver of the region were not found to store large food caches during fall and winter, however, they are believed to collect feed periodically throughout the winter months. It was observed that the severe winter in 1968-69, reduced by 27 per cent the number of active beaver colonies within the study area. The Porcupine Hills region had 0.07 beaver colony per square kilometer in 1971, a low density of beaver colonies due to the rough topography of the area. However, the importance of the beaver ponds was somewhat clarified as they provide increased moisture, which lessens the fire damage, and store water for wildlife and cattle in the area. Meteorological data was collected by the author in collaboration with the Department of Transport.

Ecological interactions between zooplankton and photosynthetic bacteria in Crawford Lake, Ontario

The present study was carried out to test the hypothesis that photosynthetic bacteria contribute a large portion of the food of filter feeding zooplankton populations in Crawford Lake, Ontario. The temporal and spatial variations of both groups of organisms are strongly dependent on one another. 14 By using C-Iabelled photosynthetic bacteria. the ingestion and clearance rates of Daphnia pulex, ~. rosea, and Keratella spp were estimated during summer and fall of 1982. These quantitative estimations of zooplankton ingestion and clearence rates on photosynthetic bacteria comprised an original addition to the literature. Photosynthetic bacteria comprised a substantial portion of the diet of all four dominant zooplankton species. The evidence for this is based on the ingestion and clearance rates of the dominant zooplankton species. Ingestion rates of D. pulex and D. rosea ranged 5 5 -1 -1 - -- 5 - -- 5 from 8.3X10 -1 to 14.6XlO -1 cells.ind. hr and 8.1X10 to 13.9X10 cells.ind. hr • Their clearance rates ranged from 0.400 to 1.000 -1 -1 -1 -1 ml.ind. hr. and 0.380 to 0.930 ml.ind. hr • The ingestion and clearance -1 -1 -1 -1 rates of Keratella spp were 600 cell.ind. hr and 0.40 ul.ind. hr respectively. Clearance rates were inversely proportional to the concentration of food cells and directly proportional to the body size of the animals. It is believed that despite the very short reg~neration times of photosynthetic bacteria (3-8 hours) their population densities were controlled in part by the feeding rates of the dominant zooplankton in Crawford Lake. By considering the regeneration times of photosynthetic bacteria and the population clearance rates of zooplankton, it was estimated that between 16 to 52% and 11 to 35% of the PHotosynthetic bacteria were' consumed· by Daphnia· pulex. and Q.. rosea per day. The temporal and spatial distribution of Daphnia pulex, !.. rosea, Keratella quadrata, K. coChlearis and photosynthetic bacteria in Crawford Lake were also investigated during the period of October, 1981 to December, 1982. The photosynthetic bacteria in the lake, constituted a major food source for only those zooplankton Which tolerate anaerobic conditions. Changes in temperature and food appeared to correlate with the seasonal changes in zooplankton density. All four dominant species of zooplankton were abundant at the lake's surface (O-4m) during winter and spring and moved downwards with the thermocline as summer stratification proceeded. Photosynthetic bacteria formed a 2 m thick layer at the chemocline. The position of this photosynthetic bacterial J-ayer changed seasonally. In the summer, the bacterial plate moved upwards and following fall mixing it moved downwards. A vertical shift of O.8m (14.5 to 15.3m) was recorded during the period of June to December. The upper limit of the photosynthetic bacteria in the water column was controlled by dissolved oxygen, and sulfide concentrations While their lower limit was controlled by light intensity. A maximum bacterio- 1 chlorophyll concentration of 81 mg Bchl.l was recorded on August 9, 1981. The seasonal distribution of photosynthetic bacteria was controlledinpart' by ·theg.-"z1ai'_.Q;~.zoopl. ank:tCm;-.Qther -ciactors associated with zooplankton grazing were oxygen and sulfide concentrations.