Unexpected findings on the taxonomic status of Est Mediterranean Crocidura russula auct. (Mammalaia, Insectivora)

With an approach, based on chromosomal, biochemical and morphological analysis, we show that the animals wich are classically attributed to C. r. gueldenstaedti and/or to C. r. monacha belong in fact to the taxon C. suaveolens. The oriental forms, formely seen as C. russula, are of the same karyotype as C. suaveolens and are biochemically close to C. suaveolens in Central Europe. The origin of the taxonomical confusion might stem from the large morphological variation in and between populations of C. suaveolens in the Near East. Based on our results C. russula monacha Thomas, 1906 belongs to the polymorphic species C. suaveolens Pallas, 1811 and we propose to consider C. (russula) gueldenstaedti Pallas, 1881 as "incertae sedis".

Genetic homogeneity of dolphinfish (Coryphaena hippurus) in the western Mediterranean and the eastern Atlantic

Dolphinfish (Coryphaena hippurus) is an epipelagic, highly migratory species distributed worldwide in tropical and temperate waters including the Mediterranean Sea. Protein electrophoresis analyses can provide knowledge of the genetic population structure of the species and therefore be used as a tool for fishery management. Areas sampled include the islands of Majorca and Sicily in the western Mediterranean and the Canary Islands in the eastern Atlantic. The results of the protein electrophoresis reveal a level of genetic variability similar to other highly migratory species. No differences were found among locations, and it was not possible to reject the null hypothesis of one panmictic population in the area studied

Lipid reserves of red mullet (Mullus barbatus) during pre-spawning in the northwestern Mediterranean

Lipid reserves are a particularly important attribute of fishes because they have a large influence on growth, reproduction and survival. This study analyses the lipid content of red mullet (Mullus barbatus) pre-spawners in three different areas of the northwestern Mediterranean in relation to trawling activities and river runoff. The muscle lipid was considered as an indicator of the somatic condition of individuals whilst the gonad lipid was used as a proxy of the energy invested in reproduction. The results show that fish with the highest muscle lipid levels inhabited the area where fishing impact was lowest. Since the abundance and biomass of polychaetes, which represent the main food source for red mullet, were found to be lower in trawled zones than in unfished ones, we suggest that differences in the muscle lipid levels between areas might be attributed to variation in prey abundance in relation to fishing impact. However, no impact of river runoff on lipid reserves of red mullet was observed. The results also show that muscle and gonad lipid reserves are not related to each other during pre-spawning

Lipid (energy) reserves of european hake (Merluccius Merluccius) in the North-Western Mediterranean

This study analyses for the first time the lipid (energy) reserves of European hake (Merluccius merluccius) in the north-western Mediterranean from an ecophysiological perspective. Results show that there is a progressive accumulation of lipids in the liver of maturing hake -where the bulk of the fat is stored- as individuals grow. Results also indicate that female pre-spawners expend much energy on reproductive activities since they present lower liver lipid reserves than juveniles and maturing individuals. Furthermore, results show that female pre-spawners with higher lipid reserves in their livers had a higher amount of lipids in their ovaries, suggesting that maternal condition (spawner quality) may affect the reproductive potential of hake. Overall, the results of this study suggest that the analysis of liver lipid reserves during pre-spawning, along with the evaluation of the gonadosomatic index and the consideration of the reproductive stage, can contribute to improve the estimation of the reproductive potential of gadoid species such as hake

Changing landscapes of Switzerland: assessing effects of historic and future land uses and land covers on the occurrence of selected butterfly species

Summary Landscapes are continuously changing. Natural forces of change such as heavy rainfall and fires can exert lasting influences on their physical form. However, changes related to human activities have often shaped landscapes more distinctly. In Western Europe, especially modern agricultural practices and the expanse of overbuilt land have left their marks in the landscapes since the middle of the 20th century. In the recent years men realised that mare and more changes that were formerly attributed to natural forces might indirectly be the result of their own action. Perhaps the most striking landscape change indirectly driven by human activity we can witness in these days is the large withdrawal of Alpine glaciers. Together with the landscapes also habitats of animal and plant species have undergone vast and sometimes rapid changes that have been hold responsible for the ongoing loss of biodiversity. Thereby, still little knowledge is available about probable effects of the rate of landscape change on species persistence and disappearance. Therefore, the development and speed of land use/land cover in the Swiss communes between the 1950s and 1990s were reconstructed using 10 parameters from agriculture and housing censuses, and were further correlated with changes in butterfly species occurrences. Cluster analyses were used to detect spatial patterns of change on broad spatial scales. Thereby, clusters of communes showing similar changes or transformation rates were identified for single decades and put into a temporally dynamic sequence. The obtained picture on the changes showed a prevalent replacement of non-intensive agriculture by intensive practices, a strong spreading of urban communes around city centres, and transitions towards larger farm sizes in the mountainous areas. Increasing transformation rates toward more intensive agricultural managements were especially found until the 1970s, whereas afterwards the trends were commonly negative. However, transformation rates representing the development of residential buildings showed positive courses at any time. The analyses concerning the butterfly species showed that grassland species reacted sensitively to the density of livestock in the communes. This might indicate the augmented use of dry grasslands as cattle pastures that show altered plant species compositions. Furthermore, these species also decreased in communes where farms with an agricultural area >5ha have disappeared. The species of the wetland habitats were favoured in communes with smaller fractions of agricultural areas and lower densities of large farms (>10ha) but did not show any correlation to transformation rates. It was concluded from these analyses that transformation rates might influence species disappearance to a certain extent but that states of the environmental predictors might generally outweigh the importance of the corresponding rates. Information on the current distribution of species is evident for nature conservation. Planning authorities that define priority areas for species protection or examine and authorise construction projects need to know about the spatial distribution of species. Hence, models that simulate the potential spatial distribution of species have become important decision tools. The underlying statistical analyses such as the widely used generalised linear models (GLM) often rely on binary species presence-absence data. However, often only species presence data have been colleted, especially for vagrant, rare or cryptic species such as butterflies or reptiles. Modellers have thus introduced randomly selected absence data to design distribution models. Yet, selecting false absence data might bias the model results. Therefore, we investigated several strategies to select more reliable absence data to model the distribution of butterfly species based on historical distribution data. The results showed that better models were obtained when historical data from longer time periods were considered. Furthermore, model performance was additionally increased when long-term data of species that show similar habitat requirements as the modelled species were used. This successful methodological approach was further applied to assess consequences of future landscape changes on the occurrence of butterfly species inhabiting dry grasslands or wetlands. These habitat types have been subjected to strong deterioration in the recent decades, what makes their protection a future mission. Four spatially explicit scenarios that described (i) ongoing land use changes as observed between 1985 and 1997, (ii) liberalised agricultural markets, and (iii) slightly and (iv) strongly lowered agricultural production provided probable directions of landscape change. Current species-environment relationships were derived from a statistical model and used to predict future occurrence probabilities in six major biogeographical regions in Switzerland, comprising the Jura Mountains, the Plateau, the Northern and Southern Alps, as well as the Western and Eastern Central Alps. The main results were that dry grasslands species profited from lowered agricultural production, whereas overgrowth of open areas in the liberalisation scenario might impair species occurrence. The wetland species mostly responded with decreases in their occurrence probabilities in the scenarios, due to a loss of their preferred habitat. Further analyses about factors currently influencing species occurrences confirmed anthropogenic causes such as urbanisation, abandonment of open land, and agricultural intensification. Hence, landscape planning should pay more attention to these forces in areas currently inhabited by these butterfly species to enable sustainable species persistence. In this thesis historical data were intensively used to reconstruct past developments and to make them useful for current investigations. Yet, the availability of historical data and the analyses on broader spatial scales has often limited the explanatory power of the conducted analyses. Meaningful descriptors of former habitat characteristics and abundant species distribution data are generally sparse, especially for fine scale analyses. However, this situation can be ameliorated by broadening the extent of the study site and the used grain size, as was done in this thesis by considering the whole of Switzerland with its communes. Nevertheless, current monitoring projects and data recording techniques are promising data sources that might allow more detailed analyses about effects of long-term species reactions on landscape changes in the near future. This work, however, also showed the value of historical species distribution data as for example their potential to locate still unknown species occurrences. The results might therefore contribute to further research activities that investigate current and future species distributions considering the immense richness of historical distribution data. Résumé Les paysages changent continuellement. Des farces naturelles comme des pluies violentes ou des feux peuvent avoir une influence durable sur la forme du paysage. Cependant, les changements attribués aux activités humaines ont souvent modelé les paysages plus profondément. Depuis les années 1950 surtout, les pratiques agricoles modernes ou l'expansion des surfaces d'habitat et d'infrastructure ont caractérisé le développement du paysage en Europe de l'Ouest. Ces dernières années, l'homme a commencé à réaliser que beaucoup de changements «naturels » pourraient indirectement résulter de ses propres activités. Le changement de paysage le plus apparent dont nous sommes témoins de nos jours est probablement l'immense retraite des glaciers alpins. Avec les paysages, les habitats des animaux et des plantes ont aussi été exposés à des changements vastes et quelquefois rapides, tenus pour coresponsable de la continuelle diminution de la biodiversité. Cependant, nous savons peu des effets probables de la rapidité des changements du paysage sur la persistance et la disparition des espèces. Le développement et la rapidité du changement de l'utilisation et de la couverture du sol dans les communes suisses entre les années 50 et 90 ont donc été reconstruits au moyen de 10 variables issues des recensements agricoles et résidentiels et ont été corrélés avec des changements de présence des papillons diurnes. Des analyses de groupes (Cluster analyses) ont été utilisées pour détecter des arrangements spatiaux de changements à l'échelle de la Suisse. Des communes avec des changements ou rapidités comparables ont été délimitées pour des décennies séparées et ont été placées en séquence temporelle, en rendrent une certaine dynamique du changement. Les résultats ont montré un remplacement répandu d'une agriculture extensive des pratiques intensives, une forte expansion des faubourgs urbains autour des grandes cités et des transitions vers de plus grandes surfaces d'exploitation dans les Alpes. Dans le cas des exploitations agricoles, des taux de changement croissants ont été observés jusqu'aux années 70, alors que la tendance a généralement été inversée dans les années suivantes. Par contre, la vitesse de construction des nouvelles maisons a montré des courbes positives pendant les 50 années. Les analyses sur la réaction des papillons diurnes ont montré que les espèces des prairies sèches supportaient une grande densité de bétail. Il est possible que dans ces communes beaucoup des prairies sèches aient été fertilisées et utilisées comme pâturages, qui ont une autre composition floristique. De plus, les espèces ont diminué dans les communes caractérisées par une rapide perte des fermes avec une surface cultivable supérieure à 5 ha. Les espèces des marais ont été favorisées dans des communes avec peu de surface cultivable et peu de grandes fermes, mais n'ont pas réagi aux taux de changement. Il en a donc été conclu que la rapidité des changements pourrait expliquer les disparitions d'espèces dans certains cas, mais que les variables prédictives qui expriment des états pourraient être des descripteurs plus importants. Des informations sur la distribution récente des espèces sont importantes par rapport aux mesures pour la conservation de la nature. Pour des autorités occupées à définir des zones de protection prioritaires ou à autoriser des projets de construction, ces informations sont indispensables. Les modèles de distribution spatiale d'espèces sont donc devenus des moyens de décision importants. Les méthodes statistiques courantes comme les modèles linéaires généralisés (GLM) demandent des données de présence et d'absence des espèces. Cependant, souvent seules les données de présence sont disponibles, surtout pour les animaux migrants, rares ou cryptiques comme des papillons ou des reptiles. C'est pourquoi certains modélisateurs ont choisi des absences au hasard, avec le risque d'influencer le résultat en choisissant des fausses absences. Nous avons établi plusieurs stratégies, basées sur des données de distribution historique des papillons diurnes, pour sélectionner des absences plus fiables. Les résultats ont démontré que de meilleurs modèles pouvaient être obtenus lorsque les données proviennent des périodes de temps plus longues. En plus, la performance des modèles a pu être augmentée en considérant des données de distribution à long terme d'espèces qui occupent des habitats similaires à ceux de l'espèce cible. Vu le succès de cette stratégie, elle a été utilisée pour évaluer les effets potentiels des changements de paysage futurs sur la distribution des papillons des prairies sèches et marais, deux habitats qui ont souffert de graves détériorations. Quatre scénarios spatialement explicites, décrivant (i) l'extrapolation des changements de l'utilisation de sol tels qu'observés entre 1985 et 1997, (ii) la libéralisation des marchés agricoles, et une production agricole (iii) légèrement amoindrie et (iv) fortement diminuée, ont été utilisés pour générer des directions de changement probables. Les relations actuelles entre la distribution des espèces et l'environnement ont été déterminées par le biais des modèles statistiques et ont été utilisées pour calculer des probabilités de présence selon les scénarios dans six régions biogéographiques majeures de la Suisse, comportant le Jura, le Plateau, les Alpes du Nord, du Sud, centrales orientales et centrales occidentales. Les résultats principaux ont montré que les espèces des prairies sèches pourraient profiter d'une diminution de la production agricole, mais qu'elles pourraient aussi disparaître à cause de l'embroussaillement des terres ouvertes dû à la libéralisation des marchés agricoles. La probabilité de présence des espèces de marais a décrû à cause d'une perte générale des habitats favorables. De plus, les analyses ont confirmé que des causes humaines comme l'urbanisation, l'abandon des terres ouvertes et l'intensification de l'agriculture affectent actuellement ces espèces. Ainsi ces forces devraient être mieux prises en compte lors de planifications paysagères, pour que ces papillons diurnes puissent survivre dans leurs habitats actuels. Dans ce travail de thèse, des données historiques ont été intensivement utilisées pour reconstruire des développements anciens et pour les rendre utiles à des recherches contemporaines. Cependant, la disponibilité des données historiques et les analyses à grande échelle ont souvent limité le pouvoir explicatif des analyses. Des descripteurs pertinents pour caractériser les habitats anciens et des données suffisantes sur la distribution des espèces sont généralement rares, spécialement pour des analyses à des échelles fores. Cette situation peut être améliorée en augmentant l'étendue du site d'étude et la résolution, comme il a été fait dans cette thèse en considérant toute la Suisse avec ses communes. Cependant, les récents projets de surveillance et les techniques de collecte de données sont des sources prometteuses, qui pourraient permettre des analyses plus détaillés sur les réactions à long terme des espèces aux changements de paysage dans le futur. Ce travail a aussi montré la valeur des anciennes données de distribution, par exemple leur potentiel pour aider à localiser des' présences d'espèces encore inconnues. Les résultats peuvent contribuer à des activités de recherche à venir, qui étudieraient les distributions récentes ou futures d'espèces en considérant l'immense richesse des données de distribution historiques.

Anomalous rainfall and associated atmospheric circulation in the northeast Spanish Mediterranean area and its relationship to sediment fluidization events in a lake

Changes in the dynamics of sediment transport in a Mediterranean lake (sediment fluidization events) are linked to atmospheric circulations patterns (trough monthly precipitation). In the basins of Lake Banyoles, located in the northeast of Spain, water enters mainly through subterranean springs, and associated fluctuations in the vertical migration of sediment distribution (fluidization events) present episodic behavior as a result of episodic rainfall in the area. The initiation of the fluidization events takes place when the monthly rainfall is ∼2.7 times greater than the mean monthly rainfall of the rainiest months in the area, especially in spring (April and May), October, and December. The duration of these events is found to be well correlated with the accumulated rainfall of the preceding 10 months before the process initiation. The rainfall, in turn, is mainly associated with six atmospheric circulation patterns among the 19 fundamental circulations that emerged in an earlier study focused on significant rainfall days in Mediterranean Spain. Among them, accentuated surface lows over the northeast of Spain, general northeasterly winds by low pressure centered to the east of Balearic Islands and short baroclinic waves over the Iberian Peninsula, with easterly flows over the northeastern coast of Spain, are found the most relevant atmospheric circulations that drive heavy rainfall events

Integrating detection probabilities in species distribution models of amphibians breeding in Mediterranean temporary ponds

Aim  The imperfect detection of species may lead to erroneous conclusions about species-environment relationships. Accuracy in species detection usually requires temporal replication at sampling sites, a time-consuming and costly monitoring scheme. Here, we applied a lower-cost alternative based on a double-sampling approach to incorporate the reliability of species detection into regression-based species distribution modelling.Location  Doñana National Park (south-western Spain).Methods  Using species-specific monthly detection probabilities, we estimated the detection reliability as the probability of having detected the species given the species-specific survey time. Such reliability estimates were used to account explicitly for data uncertainty by weighting each absence. We illustrated how this novel framework can be used to evaluate four competing hypotheses as to what constitutes primary environmental control of amphibian distribution: breeding habitat, aestivating habitat, spatial distribution of surrounding habitats and/or major ecosystems zonation. The study was conducted on six pond-breeding amphibian species during a 4-year period.Results  Non-detections should not be considered equivalent to real absences, as their reliability varied considerably. The occurrence of Hyla meridionalis and Triturus pygmaeus was related to a particular major ecosystem of the study area, where suitable habitat for these species seemed to be widely available. Characteristics of the breeding habitat (area and hydroperiod) were of high importance for the occurrence of Pelobates cultripes and Pleurodeles waltl. Terrestrial characteristics were the most important predictors of the occurrence of Discoglossus galganoi and Lissotriton boscai, along with spatial distribution of breeding habitats for the last species.Main conclusions  We did not find a single best supported hypothesis valid for all species, which stresses the importance of multiscale and multifactor approaches. More importantly, this study shows that estimating the reliability of non-detection records, an exercise that had been previously seen as a naïve goal in species distribution modelling, is feasible and could be promoted in future studies, at least in comparable systems.

Predicting the potential for lightning activity in mediterranean storms based on the WRF model dynamic and microphysical fields

A new parameter is introduced: the lightning potential index (LPI), which is a measure of the potential for charge generation and separation that leads to lightning flashes in convective thunderstorms. The LPI is calculated within the charge separation region of clouds between 0 C and 20 C, where the noninductive mechanism involving collisions of ice and graupel particles in the presence of supercooled water is most effective. As shown in several case studies using the Weather Research and Forecasting (WRF) model with explicit microphysics, the LPI is highly correlated with observed lightning. It is suggested that the LPI may be a useful parameter for predicting lightning as well as a tool for improving weather forecasting of convective storms and heavy rainfall.

A methodology for the classification of convective structures using meteorological radar: Application to heavy rainfall events on the Mediterranean coast of the Iberian Peninsula

During the period 1996-2000, forty-three heavy rainfall events have been detected in the Internal Basins of Catalonia (Northeastern of Spain). Most of these events caused floods and serious damage. This high number leads to the need for a methodology to classify them, on the basis of their surface rainfall distribution, their internal organization and their physical features. The aim of this paper is to show a methodology to analyze systematically the convective structures responsible of those heavy rainfall events on the basis of the information supplied by the meteorological radar. The proposed methodology is as follows. Firstly, the rainfall intensity and the surface rainfall pattern are analyzed on the basis of the raingauge data. Secondly, the convective structures at the lowest level are identified and characterized by using a 2-D algorithm, and the convective cells are identified by using a 3-D procedure that looks for the reflectivity cores in every radar volume. Thirdly, the convective cells (3-D) are associated with the 2-D structures (convective rainfall areas). This methodology has been applied to the 43 heavy rainfall events using the meteorological radar located near Barcelona and the SAIH automatic raingauge network.

The cases of June 2000, November 2002 and September 2002 as examples of Mediterranean floods

Four flood events that affected three different countries are here described in terms of meteorological genesis and in terms of consequences on the population and on the territory. Each event is a good representative of a class of phenomena that produce important effects on the urban and extra-urban tissue and that must be taken into account in an optic of civil protection and risk evaluation. This is the subject of the HYDROPTIMET project, part of the Interreg IIIB program, which is collocated in the framework of the prevention of natural hazards and, in particular, those related to severe meteo-hydrological events. This paper aims at being a general introduction of the four events which are the subject of more detailed studies, already published or under submission.

Towards a database on societal impact of Mediterranean floods in the framework of the HYMEX project

The NW Mediterranean region experiences every year heavy rainfall and flash floods that occasionally produce catastrophic damages. Less frequent are floods that affect large regions. Although a large number of databases devoted exclusively to floods or considering all kind of natural hazards do exist, usually they only record catastrophic flood events. This paper deals with the new flood database that is being developed within the framework of HYMEX project. Results are focused on four regions representative of the NW sector of Mediterranean Europe: Catalonia, Spain; the Balearic Islands, Spain; Calabria, Italy; and Languedoc-Roussillon, Midi-Pyrenées and PACA, France. The common available 30-yr period starts in 1981 and ends in 2010. The paper shows the database structure and criteria, the comparison with other flood databases, some statistics on spatial and temporal distribution, and an identification of the most important events. The paper also provides a table that includes the date and affected region of all the catastrophic events identified in the regions of study, in order to make this information available for all audiences.

The role of population origin and microenvironment in seedling emergence and early survival in Mediterranean maritime pine (Pinus pinaster Aiton).

Understanding tree recruitment is needed to forecast future forest distribution. Many studies have reported the relevant ecological factors that affect recruitment success in trees, but the potential for genetic-based differences in recruitment has often been neglected. In this study, we established a semi-natural reciprocal sowing experiment to test for local adaptation and microenvironment effects (evaluated here by canopy cover) in the emergence and early survival of maritime pine (Pinus pinaster Aiton), an emblematic Mediterranean forest tree. A novel application of molecular markers was also developed to test for family selection and, thus, for potential genetic change over generations. Overall, we did not find evidence to support local adaptation at the recruitment stage in our semi-natural experiment. Moreover, only weak family selection (if any) was found, suggesting that in stressful environments with low survival, stochastic processes and among-year climate variability may drive recruitment. Nevertheless, our study revealed that, at early stages of recruitment, microenvironments may favor the population with the best adapted life strategy, irrespectively of its (local or non-local) origin. We also found that emergence time is a key factor for seedling survival in stressful Mediterranean environments. Our study highlights the complexity of the factors influencing the early stages of establishment of maritime pine and provides insights into possible management actions aimed at environmental change impact mitigation. In particular, we found that the high stochasticity of the recruitment process in stressful environments and the differences in population-specific adaptive strategies may difficult assisted migration schemes.

A preliminary catalogue of natural substances of opistobranch molluscs from Western Mediterranean and near Atlantic

A list of opisthobranch molluscs species from the western Mediterranean and nearby Atlantic is presented. These species have natural products that are of interest because of their chemical structure, origin and/or function in benthic ecosystems. This review contains data on the origin and activity of these molecules, collection sites of the animals, and their bibliographic references. A discussion of these subjects is also included.

Kleptoparasitism, disturbance and predation of yellow-legged gulls on Audouin's gulls in three colonies of the western Mediterranean

The impact of yellow-legged gulls on Audouin´s gulls was studied by means of observations from a blind, at the Ebro Delta, the Chafarinas Islands and the Columbretes Islands colonies, during different stages of Audouin´s Gulls breeding cycle. The rates of predation (upon eggs and chicks) and kleptoparasitism (aerial and on courtship and chick-feeding regurgitates) were recorded to evaluate this impact. Kleptoparasitism and predation rates recorded at the three colonies may be considered low when compared with similar studies dealing with other gull species. It seems that interactions did not represent a threat to the population dynamics of Audouin´s Gulls in any of the three colonies. The highest rates of aerial kleptoparasitism, courtship feeding and chick mortality were recorded at the Columbretes Islands, and the lowest at the Ebro Delta. The ratio of the number of yellow-legged gulls to the number of Audouin´s gulls seems to influence the rate of disturbances. Food availability was also likely to play an important role in the number of interactions, which increased when food was in shorter supply. The study at the Ebro Delta and Chafarinas Islands was restricted to some subcolonies and hence results may not be necessarily representative of the entire colonies. The existence of some yellow-legged gull individuals specialised as predators was recorded. Finally, we discuss the suitability of several conservation measures commonly applied in colonies where both species breed syntopically.

Spatial and temporal variation of natural toxicity in cnidarians, bryozoans and tunicates in Mediterranean caves

The natural toxicity of cnidarians, bryozoans and tunicates in two caves was assessed using the Microtox® technique in spring and autumn. One cave was located in the Cabrera Archipelago (Balearic Islands) and the other in the Medes Islands (Catalan littoral). The organisms analysed were good representatives of the coverage of each Phylum in the communities; however, these Phyla are less abundant than sponges which are the dominant group in these caves. Seventy-one percent of the species of cnidarians and bryozoans analysed were toxic in one of the caves, communities or seasons, which indicates the relevance of bioactive species in these groups. The tunicate Lissoclinum perforatum was the most toxic species. Although all three Phyla had some highly toxic species, a common pattern that related the caves, communities and seasons was not found. Seasonal variation of toxicity in cnidarians and bryozoans was higher in the Cabrera than in the Medes cave. Moreover, variation in toxicity either between communities or between seasons was a common trait for most cnidarians and bryozoans, whereas tunicates remained toxic throughout communities and seasons.