Bias in freshwater biodiversity sampling: the case of Iberian water beetles

Extensive distributional data bases are key tools in ecological research, and good-quality data are required to provide reliable conservation strategies and an understanding of biodiversity patterns and processes. Although the evaluation of data bases requires the incorporation of estimates of sampling effort and bias, no studies have focused on these aspects for freshwater biodiversity data. We used here a comprehensive data base of water beetles from the Iberian Peninsula and the Balearic Islands, and examine whether these data provide an unbiased, reliable picture of their diversity and distribution in the study area. Based on theoretical estimates using the Clench function on the accumulated number of records as a surrogate of sampling effort, about a quarter of the Iberian and Balearic 50 × 50 km Universal Transverse Mercator grid cells can be considered well prospected, with more than 70% of the theoretical species richness actually recorded. These well-surveyed cells are not evenly distributed across biogeographical and physicoclimatic subregions, reflecting some geographical bias in the distribution of sampling effort. Our results suggest that recording was skewed by relatively simple variables affecting collector activity, such as the perceived 'attractiveness' of mountainous landscapes and protected areas with recently described species, and accessibility of sampling sites (distance from main research centres). We emphasize the importance of these evaluation exercises, which are useful to locate areas needed of further sampling as well as to identify potential biases in the distribution of current biodiversity patterns.     

Butterfly species richness in mainland Portugal: predictive models of geographic distribution patterns

A three-step protocol described elsewhere is used to obtain a map of butterfly species density in Portugal on a 50×50 km grid. First, all available faunistic information was compiled and analysed to explore the historic patterns of butterfly sampling in Portugal, and to determine which grid cells are sufficiently prospected to produce reliable estimates of species richness. Then, we relate the estimated species richness scores from these areas to a set of environmental and spatial variables by means of General Linear Models, obtaining a function to extrapolate of species density scores to the rest of Portugal. Finally, the model is validated, results explored and outliers identified and deleted. Any spatial autocorrelation remaining in the residuals is examined. Lastly, model parameters are recalculated in absence of deleted outliers, and the resulting function is used to predict species richness scores throughout mainland Portugal. A highly-predictive function based on some variables previously related to butterfly composition at macro-scale, such as number of sunny days per year, temperature or environmental heterogeneity, was obtained. However, in Portugal those variables are highly spatially structured along a steep latitudinal gradient, leading to difficulty in ascertaining if the latitudinal gradient detected by our analysis is due to macroecological or historic effects. Information on European and Iberian butterfly assemblages and causal processes are discussed in the light of the patterns observed. Then, previous information obtained on Portuguese scarabs is added to identify conservation areas, biogeographically important for both insect groups. Finally, the main drawbacks and advantages of this approach to mapping biodiversity for conservation are discussed briefly.     

The ‘rotiferologist’ effect and other global correlates of species richness in monogonont rotifers

Global biodiversity patterns are often driven by different environmental variables at different scales. However, it is still controversial whether there are general trends, whether similar processes are responsible for similar patterns, and/or whether confounding effects such as sampling bias can produce misleading results. Our aim is twofold: 1) assessing the global correlates of diversity in a group of microscopic animals little analysed so far, and 2) inferring the influence of sampling intensity on biodiversity analyses. As a case study, we choose rotifers, because of their high potential for dispersal across the globe. We assembled and analysed a new worldwide dataset of records of monogonont rotifers, a group of microscopic aquatic animals, from 1960 to 1992. Using spatially explicit models, we assessed whether the diversity patterns conformed to those commonly obtained for larger organisms, and whether they still held true after controlling for sampling intensity, variations in area, and spatial structure in the data. Our results are in part analogous to those commonly obtained for macroorganisms (habitat heterogeneity and precipitation emerge as the main global correlates), but show some divergence (potential absence of a latitudinal gradient and of a large‐scale correlation with human population). Moreover, the effect of sampling effort is remarkable, accounting for >50% of the variability; this strong effect may mask other patterns such as latitudinal gradients. Our study points out that sampling bias should be carefully considered when drawing conclusions from large‐scale analyses, and calls for further faunistic work on microorganisms in all regions of the world to better understand the generality of the processes driving global patterns in biodiversity.     

Iberian–Balearic fern regions and their explanatory variables

This article delineates the compositional regions present in the Iberian–Balearic fern flora and compares these regions to previously proposed biogeographic units. It also assesses the extent to which environmental variables could explain the regions and the fern species richness gradients found within them. A combination of 40 previously published and new maps were used to compile the distribution of 123 pteridophytes on a 50 × 50 km UTM grid. Cluster analysis of the resulting 257 squares was used to classify 10 regions based on fern species assemblages. Discriminant function analysis identified the environmental variables that best explained these fern composition regions. Using generalized linear models; the number of species in each square was regressed against topography, climate, geology, environmental diversity, land use and spatial variables within each region. Two main latitudinal pteridophyte zones can be recognized in the Iberian Peninsula. These two zones are longitudinally subdivided into two sub zones. The 10 regions established significantly differ both in species richness and influential environmental variables. Climatic variables discriminate the most among regions, followed by topography, heterogeneity and geology. Pteridophyte richness varies, with richer areas being located along the coast and the main mountain ranges and the poorest areas being in the central plateaus and some north eastern and south western river basins. Species richness variation in Iberia is positively correlated with altitude range, precipitation, maximum altitude and area with siliceous soils. It is negatively correlated with the total annual days of sun, however. The fact that species richness is explained by different variables within each of the 10 regions indicates that the specific factors determining the spatial distribution of species richness vary from region to region. Some coastal regions are poorly explained by the model, and display a negative correlation with the selected causal factors. This finding suggests that persistent historic effects might play a local role in determining species assemblages in these regions. An erratum to this article can be found at     

Biogeographical analysis of species of the tribe Cytiseae (Fabaceae) in the Iberian Peninsula and Balearic Islands

Aim To analyse quantitatively the biogeographical distribution pattern of species of the Cytiseae Bercht. & J. Presl (= Genisteae Benth.) tribe in the Iberian Peninsula and the Balearic Islands, and to identify environmental variables related to the distributional patterns. Location Iberian Peninsula and Balearic Islands, using the 61 administrative provinces as operational geographical units. Methods In order to identify chorotypes (groups of species with similar geographical distribution), we performed a upgma classification based on the similarity index of Baroni‐Urbani & Buser. The method of McCoy et al. [Ecology 67 (1986), 749] enabled us to detect the significant groups and to differentiate them from those groupings that could be generated at random. Logistic regression analyses and environmental gradient analyses (DCA and CCA) were performed in order to find the relationships between the environmental variables and the observed distributional patterns. Results Sixteen chorotypes were obtained between the Cytiseae species of the Iberian Peninsula and the Balearic Islands. The thermal variables showed the greatest influence in species distribution. Specifically, temperatures (maximum, minimum and mean) of the coldest months were associated with the majority of the chorotypes. Main conclusions The species of the Cytiseae tribe were not randomly distributed in space, and can be classified in groups of species with common distributional patterns. The importance of cold tolerance in the distribution of these species, as well as their preference for acidic soils, was demonstrated. Certain general tendencies seem to exist with respect to the distribution of the biota in the Iberian Peninsula, and these seem to be independent of the taxonomic groups involved and to be determined by convergent macroclimatic factors.     

Spatial and environmental determinants of vascular plant species richness distribution in the Iberian Peninsula and Balearic Islands

Using an exhaustive data compilation, Iberian vascular plant species richness in 50 times 50 UTM grid cells was regressed against 24 explanatory variables (spatial, geographical, topographical, geological, climatic, land use and environmental diversity variables) using Generalized Linear Models and partial regression analysis in order to ascertain the relative contribution of primary, heterogeneous and spatially structured variables. The species richness variation accounted for by these variables is reasonably high (65% of total deviance). Little less than half of this variation is accounted for spatially structured variables. A purely spatial component of variation is hardly significant. The most significant variables are those related to altitude, and particularly maximum altitude, whose cubic response reflects the occurrence of the maximum number of species at the highest altitudes. This result highlighted the importance of Iberian mountains as hotspots of diversity and the relevance of large and small scale historical factors in contemporary plant distribution patterns. Climatic or energy-related variables contributed little, whereas geological (calcareous and acid rocks) and, to a lesser extent, environmental heterogeneity variables (land use diversity and altitude range) seem to be more important.     

Spatial distribution and habitats of useful plants: an initial assessment for conservation on an indigenous territory, Panama

In order to examine the spatial distribution of forest resources on local territories and to understand the factors controlling such distributions, we studied the spatial patterns of a group of 23 useful plant species on the territory of a Kuna community in the province of Darien, Panama. A stratified random sampling scheme was used to survey the distribution and abundance of the species across a 3500 ha area around the village. Data on the physical environment as well as the geographic coordinates of the sample plots were also obtained. A series of canonical analyses was conducted to evaluate the species–environment relationships and to identify spatial structures in the species distributions left unexplained by the environmental variables. Four distinct distribution patterns were identified among the species; these were most strongly explained by land-use, the degree of canopy closure and topography. Significant spatial structures, independent of the measured environmental variables, were related to anthropogenic pressure and an edaphic gradient, and the habitat associations of the individual species were described. The results obtained from this case study suggest that land-use dynamics may play a predominant role in structuring inhabited landscapes, and that diversity in distribution patterns and habitat associations will require a combination of spatially explicit management strategies to ensure the local resource base.     

Linkages among land-use, water quality, physical habitat conditions and lotic diatom assemblages: A multi-spatial scale assessment

We assessed the importance of spatial scales (catchment, stream network, and sample reach) on the effects of agricultural land-use on lotic diatom assemblages along a land-use gradient in the agricultural Willamette Valley Ecoregion of Oregon. Periphyton, water chemistry, and physical habitat conditions were characterized for 25 wadeable streams during a dry season (July to September, 1997). Additional water chemistry samples were collected in the following wet season (February 1998) to assess seasonal effects of land-use on stream water chemistry. Percent agricultural land-use in the study catchments ranged from 10% to 89% with an average of 52%. Partial canonical correspondence analysis (CCA) with the first axis constrained by % agricultural land-use showed that % agricultural land-use at 3 spatial scales explained between 3.7%–6.3% of variability in the diatom species dataset. Monte Carlo Permutation tests indicated that the variance explained by % agricultural land-use was only significant at the spatial scale of the stream network with 10- and 30-m band width (p<0.05, 999 permutations). In addition to the effects of % agricultural land-use, partial CCAs with a forward selection option showed that water chemistry (e.g., SiO₂), reach-scale stream channel dimensions (e.g., width, depth, and slope), reach-scale in-stream habitats (substrates and filamentous algal cover in stream beds), and riparian vegetative buffer were all important with relation to diatom species assemblages. Percent of obligately nitrogen-heterotrophic taxa was the only diatom autecological metric that showed a significant but weak correlation with % agricultural land-use along the stream network (r=0.50), but not at catchment or sample reach scale. Correlation between % agricultural land-use and water chemistry variables varied among the spatial scales and between seasons. Physical habitat variables (log₁₀ erodible substrate diameters and stream reach slope) were significantly correlated with % agricultural land-use along the stream network but not at catchment or sample reach scale. Our data suggest that spatial scales are important in assessing effects of land-use on stream conditions but the spatial scale effects may vary between seasons. Direct linkages between agricultural land-use and lotic diatom assemblages were weak during summer base-flow time regardless of the spatial scales. Summer sampling may underestimate the effects of catchment land-use on stream conditions in areas where seasonal patterns are so distinctive as in the Willamette Valley.     

Environmental and geographical constraints on common swift and barn swallow spring arrival patterns throughout the Iberian Peninsula

Aim Still poorly understood, the main migratory pathways for most trans‐Saharan species pass through the Iberian Peninsula, which acts as a gateway to the European–African migratory system. Arrival patterns in this region for the common swift (Apus apus) and barn swallow (Hirundo rustica), of similar morphology and flight capabilities, were described, and the environmental and geographical factors best explaining them were examined, in a search for common ecological constraints on these two migratory species. Location Latitude ranged from 36.02 to 43.68°N, longitude from 9.05°W to 3.17°E, and altitude from 0 to 1595 m a.s.l. for 482 common swift and 812 barn swallow Spanish localities spread widely over the Iberian breeding grounds of the two species. Methods Our data set, covering the years 1960–1990, consisted of 3206 first‐arrival dates for common swifts and 6036 for barn swallows. Forty topographical, climatic, river basin, geographical and spatial variables were used as explanatory variables in general regression models (GRMs). GRMs included polynomial terms up to cubic functions in all variables when they were significant. A backward stepwise selection procedure was applied in all models until only significant terms remained. GRMs were applied in two steps. First, we searched for the best model in each one of the five types of variables (topographical, climatic, river basin, geographical and spatial). To cope with the unavoidable correlation between explanatory variables, the relative importance of each type of variable was assessed by hierarchical variance partitioning. Secondly, we searched for that model able to explain the maximum amount of the observed variability in arrival date. To obtain this model all significant explanatory variables were subjected jointly to a GRM. Spatial variables were then added to this model to take any remaining spatial structure in the data into account. Moran's I autocorrelation coefficient was used to check for spatial autocorrelation. Results Both species arrived earlier in the south‐western Iberian Peninsula, where summers are warmer and drier. From there, both species followed the main southern Iberian river basins towards the north‐east; however, several mountainous regions impede the colonization of eastern Iberia. The best models for each type of variable explained 19–47% of the variability in common swift arrival dates and 14–44% in barn swallow arrival dates. Variance partitioning indicated that climatic and geographical variables best explained variability. The best predictive models built with all variables accounted for 52% of the variability in common swift arrival dates and 50% for the barn swallow. Residuals from both models were not spatially autocorrelated, an indication that all major spatially structured variation had been accounted for. Main conclusions Spring arrival patterns are highly dependent on the geographical configuration of the Iberian Peninsula. This spatial constraint forces both species to converge very closely in their spring migration, because common swifts and barn swallows are subject to a trade‐off between optimum migratory pathways and territories ecologically suitable for breeding.     

Choice of biota in stream assessment and monitoring programs in tropical streams: A comparison of diatoms,macroinvertebrates and fish

The objective of this study was to compare the responses of diatoms, macroinvertebrates and fish to agriculture, urbanization and mining in the Manyame River Basin. Water quality sampling and benthic diatom, macroinvertebrate and fish community data were collected in April (end of the rain season) and September (dry season) 2013 at 44 sampling stations spread out across three land-use categories: commercial agricultural, communal agricultural and urban-mining areas. Commercial agricultural areas were relatively pristine as they were characterized by mature deciduous riparian forest strips which acted as riparian buffers thus protecting water resources from nonpoint source pollution. In communal agricultural areas a combination of poor agricultural practices (stream bank cultivation, overgrazing, soil erosions) and high human population densities had negative effects on water quality of streams draining these areas. Streams in urban-mining areas were highly stressed, being impacted primarily by physical habitat degradation and both point and nonpoint sources of pollution. A suite of environmental variables that varied with land-use pattern was assessed to find the combination of variables that best explained patterns of biota community composition. Community metrics i.e. the Trophic Diatom Index (TDI) based on diatoms, the South African Scoring system version 5 (SASS 5) based on macroinvertebrates and the Fish Assemblage Integrity Index (FAII) were used to determine the ecological status of study streams in relation to human-induced stressors. Data were also subjected to multivariate statistical techniques; canonical correspondence analysis (CCA), mantel test and cluster analysis to determine environmental gradients along which the diatom, macroinvertebrate and fish assemblages were distributed as well as to elucidate hypothesized differences in response to stressors among communities per land-use type. Using CCA, we assessed the individual importance of a set of environmental variables on each biotic community structure. ANOVA, showed a significant difference (p < 0.05) in physical and chemical variables among commercial agricultural, communal agricultural and urban-mining sampling stations with no significant differences (p > 0.05) between the 2 sampling periods. Based on CCAs carried out using individual variables, the strengths of relationships between diatoms and macroinvertebrates was generally high for nutrient levels, organic and metal pollution and other variables. However, fish assemblages showed a relatively low association with all water quality variables in the study; this might be explained by the high abundance of omnivores and air breathers which are able to tolerate a variety of environmental conditions. These patterns were also confirmed by the mantel test as well as the other CCAs carried out to investigate the simultaneous effects of environmental variables. These findings indicate that diatoms are more powerful indicators in accessing ecological stream/river quality and have potential for application in routine monitoring programs in tropical streams.     

Masting in Buxus balearica: assessing fruiting patterns and processes at a large spatial scale

Masting consists of the synchronous highly variable seed production among years by a plant population. We studied spatiotemporal variation in fruit production in ten populations of Buxus balearica (six in the Balearic Islands and four in the Iberian Peninsula) from 2001 to 2004 in the light of masting. In some of them we assessed, by means of both observational and experimental data, the relationship between fruit production and some abiotic variables, the role of previous reproduction, the "pollination efficiency" and the "predator satiation" hypotheses, as well as the consequences for seedling density and survival. Fruit production in B. balearica showed substantial between-year variation, especially in island compared to mainland populations. Correlative evidence indicated that this variation and its geographic pattern were related to differences between regions in rainfall variability, cost of reproduction and the degree of ambophily. We found no indication of predator satiation. However, experimental tests failed to support many of our results, namely a negative effect of previous reproduction on future flowering in island populations and lower pollen limitation with increasing flower production. We therefore warn against exclusive reliance on correlations when testing hypotheses related to masting. In addition, seedling recruitment increased after some episodes of high fruit production but probably additional factors had a role in recruitment, suggesting that mast events not always translate into increased reproductive success. Although a limited time series only allows considering B. balearica as showing 'putative' masting, weather and pollination-related processes are good candidates for further exploration of fruiting patterns and processes at a large spatial scale.     

Revealing biogeographical patterns by nonlinear ordinations and derived anisotropic spatial filters

Aim Spatial floristic and faunistic data bases promote the investigation of biogeographical gradients in relation to environmental determinants on regional to continental scales. Our aim was to extract major gradients in the distribution of vascular plant species from a grid‐based inventory (the German FLORKART data base) and relate them to long‐term precipitation and temperature records as well as soil conditions. We present an ordination technique capable of coping with this complex data array. The goal was also to sort out the influence of spatial autocorrelation, assuming floristic autocorrelation is anisotropic. Location Germany, at a spatial resolution of 6′ × 10′. Methods Isometric feature mapping (Isomap) was applied as a nonlinear ordination method. Isomap was coupled to ‘eigenvector‐based filters’ for generating spatial reference models representing spatial autocorrelation. What is novel here is that the derived filters are not based on the assumption of equidirectional autocorrelation. Instead, the so‐called ‘principal coordinates of anisotropic neighbour matrices’ build filters to test the influence of geographical vicinity in directions of high similarity among observations. Results The Isomap ordination of floristic data explained more than 95% of the data variance in six dimensions. The leading two dimensions (representing about 80% of the FLORKART data variance) revealed clear spatial gradients that could be related to independent effects of temperature, precipitation and soil observations. By contrast, the third and higher FLORKART dimensions were dominated by an antagonism of anisotropic spatial autocorrelation and soil conditions. A subsequent cluster analysis of the floristic Isomap coordinates educed the spatial organization of the floristic survey, indicating a considerable sampling bias. Conclusions We showed that Isomap provides a consistent methodical framework for both ordination and derived spatial filters. The technique is useful for tracing the often nonlinear features of species occurrence data to environmental drivers, taking into account anisotropic spatial autocorrelation. We also showed that sampling biases are a conspicuous source of variance in a frequently used floristic data base.     

Balearic insular isolation and large continental spread framed the phylogeography of the western Mediterranean Cheirolophus intybaceus s.l. (Asteraceae)

Recent Quaternary geological and climate events have shaped the evolutionary histories of plant species in the Mediterranean basin, one of the most important hotspots of biodiversity. Genetic analyses of the western Mediterranean Cheirolophus intybaceus s.l. (Asteraceae) based on AFLP were conducted to establish the relationships between its close species and populations, to reconstruct the phylogeography of the group and to analyse potential unidirectional versus bidirectional dispersals between the Ibero‐Provençal belt and the Balearic Islands. AFLP data revealed two main genetic groups, one constituted by the Balearic populations and Garraf (NE Iberia) and the other formed by the remaining mainland populations that were further sub‐structured into two geographically separated subgroups (SE + E Iberia and NE Iberia + SW France). Genetic diversity and spatial structure analyses suggested a mid‐Pleistocene scenario for the origin of C. intybaceus in southern Iberia, followed by dispersal to the north and a single colonisation event of the Balearic archipelago from the near Dianic NE Iberian area. This hypothesis was supported by paleogeographic data, which showed the existence of terrestrial connections between the continent and the islands during the Middle–Late Pleistocene marine regressions, whereas the more recent single back‐colonisation of the mainland from Mallorca might be explained by several hypotheses, such as long‐distance dispersal mediated by migratory marine birds or sea currents.     

Regional Insect Inventories Require Long Time,Extensive Spatial Sampling and Good Will

Understanding how faunistic knowledge develops is of paramount importance to correctly evaluate completeness of insect inventories and to plan future research at regional scale, yet this is an unexplored issue. Aim of this paper was to investigate the processes that lead to a complete species inventory at a regional level for a beetle family. The tenebionid beetles of Latium region (Italy) were analysed as a case study representative of general situations. A comprehensive faunistic database including 3,561 records spanning from 1871 to 2010 was realized examining 25,349 museum specimens and published data. Accumulation curves and non-parametric estimators of species richness were applied to model increase in faunistic knowledge over time, through space and by collectors’ number. Long time, large spatial extent and contribution of many collectors were needed to obtain a reliable species inventory. Massive sampling was not effective in recovering more species. Amateur naturalists (here called parafaunists) were more efficient collectors than professional entomologists. Museum materials collected by parafaunists over long periods and large spatial extent resulted to be a cost effective source of faunistic information with small number of collected individuals. It is therefore important to valuate and facilitate the work of parafaunists as already suggested for parataxonomists. By contrast, massive collections by standardized techniques for ecological research seem to be of scarce utility in improving faunistic knowledge, but their value for faunistic studies may be enhanced if they are conducted in poorly surveyed areas.     

The Balearic Islands: a reservoir of cpDNA genetic variation for evergreen oaks

Aim To analyse the role of the Balearic Islands as a refuge area for evergreen Quercus (cork oak: Quercus suber L., holm oak: Q. ilex L., kermes oak: Q. coccifera L.), by using molecular, historical and palaeobotanical data. Location The Western Mediterranean Basin (Balearic Islands, eastern Iberia, Provence, Sardinia, Corsica, Sicily, Malta, Italy, Northern Africa). Methods We sampled 108 populations and used the PCR‐RFLP technique with five universal cpDNA primers to define haplotypes in the sampled populations. Diversity, differentiation parameters and spatial analysis of the populations, using a spatial version of amova , were linked to the geological history of the Western Mediterranean Basin in order to explain the present spatial pattern of the evergreen Quercus populations in the Balearics. Results Evergreen Quercus cpDNA shows a complex structure, with remnants of ancient diversity in the Balearics. Balearic populations of holm oak are related to Iberian populations, while for cork and kermes oaks, we found both Tyrrhenian and Iberian haplotypes. Main conclusions The complex spatial patterns of cpDNA in Balearic evergreen Quercus appears explicable in terms of a combination of physical (vicariance and long distance dispersal) and biological (introgressive hybridization) factors. The Balearics constitute a glacial refuge area and a reservoir of genetic variation with traces of ancient diversity from Messinian–Pliocene stages.     

Patterns of nuclear and mitochondrial DNA variation in Iberian populations of <Emphasis Type="Italic">Emys orbicularis</Emphasis> (Emydidae): conservation implications

The European pond turtle (Emys orbicularis) is threatened and in decline in several regions of its natural range, due to habitat loss combined with population fragmentation. In this work, we have focused our efforts on studying the genetic diversity and structure of Iberian populations with a fine-scale sampling (254 turtles in 10 populations) and a representation from North Africa and Balearic island populations. Using both nuclear and mitochondrial markers (seven microsatellites, ∼1048 bp nDNA and ∼1500 bp mtDNA) we have carried out phylogenetic and demographic analyses. Our results show low values of genetic diversity at the mitochondrial level although our microsatellite dataset revealed relatively high levels of genetic variability with a latitudinal genetic trend decreasing from southern to northern populations. A moderate degree of genetic differentiation was estimated for Iberian populations (genetic distances, F _ST values and clusters in the Bayesian analysis). The results in this study combining mtDNA and nDNA, provide the most comprehensive population genetic data for E. orbicularis in the Iberian Peninsula. Our results suggest that Iberian populations within the Iberian–Moroccan lineage should be considered as a single subspecies with five management units, and emphasize the importance of habitat management rather than population reinforcement (i.e. captive breeding and reintroduction) in this long-lived species.     

Protecting persistent dynamic oceanographic features: transboundary conservation efforts are needed for the critically endangered Balearic shearwater

The protection of key areas for biodiversity at sea is not as widespread as on land and research investment is necessary to identify biodiversity hotspots in the open ocean. Spatially explicit conservation measures such as the creation of representative networks of marine protected areas (MPAs) is a critical step towards the conservation and management of marine ecosystems, as well as to improve public awareness. Conservation efforts in ecologically rich and threatened ecosystems are specially needed. This is particularly urgent for the Mediterranean marine biodiversity, which includes highly mobile marine vertebrates. Here, we studied the at sea distribution of one of the most endangered Mediterranean seabird, the critically endangered Balearic shearwater Puffinus mauretanicus. Present knowledge, from vessel-based surveys, suggests that this species has a coastal distribution over the productive Iberian shelf in relation to the distribution of their main prey, small pelagic fish. We used miniaturised satellite transmitters to determine the key marine areas of the southern population of Balearic shearwaters breeding on Eivissa and spot the spatial connections between breeding and key marine areas. Our tracking study indicates that Balearic shearwaters do not only forage along the Iberian continental shelf but also in more distant marine areas along the North African coast, in particular W of Algeria, but also NE coast of Morocco. Birds recurrently visit these shelf areas at the end of the breeding season. Species distribution modelling identified chlorophyll a as the most important environmental variable in defining those oceanographic features characterizing their key habitats in the western Mediterranean. We identified persistent oceanographic features across time series available in the study area and discuss our results within the current conservation scenario in relation to the ecology of the species.     

Spatial predictability of juvenile fish species richness and abundance in a coral reef environment

Juvenile reef fish communities represent an essential component of coral reef ecosystems in the current focus of fish population dynamics and coral reef resilience. Juvenile fish survival depends on habitat characteristics and is, following settlement, the first determinant of the number of individuals within adult populations. The goal of this study was to provide methods for mapping juvenile fish species richness and abundance into spatial domains suitable for micro and meso-scale analysis and management decisions. Generalized Linear Models predicting juvenile fish species richness and abundance were developed according to spatial and temporal environmental variables measured from 10 m up to 10 km in the southwest lagoon of New Caledonia. The statistical model was further spatially generalized using a 1.5-m resolution, independently created, remotely sensed, habitat map. This procedure revealed that : (1) spatial factors at 10 to 100-m scale explained up to 71% of variability in juvenile species richness, (2) a small improvement (75%) was gained when a combination of environmental variables at different spatial and temporal scales was used and (3) the coupling of remotely sensed data, geographical information system tools and point-based ecological data showed that the highest species richness and abundance were predicted along a narrow margin overlapping the coral reef flat and adjacent seagrass beds. Spatially explicit models of species distribution may be relevant for the management of reef communities when strong relationships exist between faunistic and environmental variables and when models are built at appropriate scales.     

Low effectiveness of the Natura 2000 network in preventing land-use change in bat hotspots

One of the most important issues in biodiversity conservation is an exploration of the relationships among protected areas, land-use changes and biodiversity, so we aimed to assess the performance of the Natura 2000 network (N2000) in representing the bat conservation hotspots in peninsular Spain and the Balearic Islands and to compare the rates of land-use changes within these hotspots with those observed throughout the rest of the study area. First, we applied a Combined Index that integrates various biodiversity metrics (species richness, rarity and vulnerability) to identify hotspots, and once they were identified, we used null models to assess the performance of N2000 in representing them. Finally, also using null models, we tested whether the changes in land use (“anthropization” or “naturalization”) within the hotspots occurred at a significantly higher or lower frequency than in the rest of the study area; for this, we considered two temporal windows (1980–2006 and 2006–2012) corresponding with periods before and after the official designation of the N2000 sites. Our results show that bat hotspots are effectively represented in the Iberian N2000, but although land-use changes were generally higher in Spain before 2006, hotspots have not experienced lower rates of change compared to the remainder of the territory (regardless of the period under consideration). This suggests low effectiveness of the Iberian N2000 in preventing land-use change in bat hotspots, so to preserve the Iberian bat fauna, we encourage the urgent implementation of management plans to avoid intensive changes in land use both inside and around bat hotspots.