Empirical modelling of lake water-temperature relationships: a comparison of approaches

1.  As a result of the role that temperature plays in many aquatic processes, good predictive models of annual maximum near-surface lake water temperature across large spatial scales are needed, particularly given concerns regarding climate change. Comparisons of suitable modelling approaches are required to determine their relative merit and suitability for providing good predictions of current conditions. We developed models predicting annual maximum near-surface lake water temperatures for lakes across Canada using four statistical approaches: multiple regression, regression tree, artificial neural networks and Bayesian multiple regression.
2.  Annual maximum near-surface (from 0 to 2 m) lake water-temperature data were obtained for more than 13 000 lakes and were matched to geographic, climatic, lake morphology, physical habitat and water chemistry data. We modelled 2348 lakes and three subsets thereof encompassing different spatial scales and predictor variables to identify the relative importance of these variables at predicting lake temperature.
3.  Although artificial neural networks were marginally better for three of the four data sets, multiple regression was considered to provide the best solution based on the combination of model performance and computational complexity. Climatic variables and date of sampling were the most important variables for predicting water temperature in our models.
4.  Lake morphology did not play a substantial role in predicting lake temperature across any of the spatial scales. Maximum near-surface temperatures for Canadian lakes appeared to be dominated by large-scale climatic and geographic patterns, rather than lake-specific variables, such as lake morphology and water chemistry.     

Evaluating model transferability for a threatened species to adjacent areas: Implications for rock‐wallaby conservation

When modelling the distribution of a species, it is often not possible to comprehensively sample the whole distribution of the species and managers may have habitat models based on data from one area that they want to apply in other areas. Hence, an important question is: how accurate are models of the distributions of species when applied beyond the areas where they were developed? A first step in measuring model transferability could be testing models in adjacent areas. We predicted the habitat associations of the brush‐tailed rock‐wallaby (Petrogale penicillata) across two spatial scales in two neighbouring study areas in eastern Australia, south‐east Queensland and north‐east New South Wales. We used classification trees for exploratory data analysis of habitat relationships and then applied logistic regression models to predict species occurrence. We assessed the within‐area discriminative ability of the habitat models using cross‐validation and threshold plots, and tested the predictive ability of the models for adjacent areas using the receiver operating characteristic statistic to determine the area under the curve. We found that models performed well within an area and extrapolating them to adjacent areas resulted in good predictive performance at the site scale but substantially poorer predictive performance at the landscape scale. We conclude that distribution models for wildlife species should only be extrapolated to neighbouring areas with caution when using landscape‐scale environmental variables. Alternatively, only key habitat associations predicted by the models at this scale should be transferred across adjacent areas once verified against local knowledge of the ecology of the study species.     

Validation of species–climate impact models under climate change

Increasing concern over the implications of climate change for biodiversity has led to the use of species–climate envelope models to project species extinction risk under climate‐change scenarios. However, recent studies have demonstrated significant variability in model predictions and there remains a pressing need to validate models and to reduce uncertainties. Model validation is problematic as predictions are made for events that have not yet occurred. Resubstituition and data partitioning of present‐day data sets are, therefore, commonly used to test the predictive performance of models. However, these approaches suffer from the problems of spatial and temporal autocorrelation in the calibration and validation sets. Using observed distribution shifts among 116 British breeding‐bird species over the past ～20 years, we are able to provide a first independent validation of four envelope modelling techniques under climate change. Results showed good to fair predictive performance on independent validation, although rules used to assess model performance are difficult to interpret in a decision‐planning context. We also showed that measures of performance on nonindependent data provided optimistic estimates of models' predictive ability on independent data. Artificial neural networks and generalized additive models provided generally more accurate predictions of species range shifts than generalized linear models or classification tree analysis. Data for independent model validation and replication of this study are rare and we argue that perfect validation may not in fact be conceptually possible. We also note that usefulness of models is contingent on both the questions being asked and the techniques used. Implementations of species–climate envelope models for testing hypotheses and predicting future events may prove wrong, while being potentially useful if put into appropriate context.     

An evaluation of methods for modelling species distributions

Aim Various statistical techniques have been used to model species probabilities of occurrence in response to environmental conditions. This paper provides a comprehensive assessment of methods and investigates whether errors in model predictions are associated to specific kinds of geographical and environmental distributions of species. Location Portugal, Western Europe. Methods Probabilities of occurrence for 44 species of amphibians and reptiles in Portugal were modelled using seven modelling techniques: Gower metric, Ecological Niche Factor Analysis, classification trees, neural networks, generalized linear models, generalized additive models and spatial interpolators. Generalized linear and additive models were constructed with and without a term accounting for spatial autocorrelation. Model performance was measured using two methods: sensitivity and Kappa index. Species were grouped according to their spatial (area of occupancy and extent of occurrence) and environmental (marginality and tolerance) distributions. Two‐way comparison tests were performed to detect significant interactions between models and species groups. Results Interaction between model and species groups was significant for both sensitivity and Kappa index. This indicates that model performance varied for species with different geographical and environmental distributions. Artificial neural networks performed generally better, immediately followed by generalized additive models including a covariate term for spatial autocorrelation. Non‐parametric methods were preferred to parametric approaches, especially when modelling distributions of species with a greater area of occupancy, a larger extent of occurrence, lower marginality and higher tolerance. Main conclusions This is a first attempt to relate performance of modelling techniques with species spatial and environmental distributions. Results indicate a strong relationship between model performance and the kinds of species distributions being modelled. Some methods performed generally better, but no method was superior in all circumstances. A suggestion is made that choice of the appropriate method should be contingent on the goals and kinds of distributions being modelled.     

Community versus single‐species distribution models for British plants

Aim Species distribution models are increasingly used to predict the impacts of global change on whole ecological communities by modelling the individualistic niche responses of large numbers of species. However, it is not clear whether this single‐species ensemble approach is preferable to community‐wide strategies that represent interspecific associations or shared responses to environmental gradients. Here, we test the performance of two multi‐species modelling approaches against equivalent single‐species models. Location Great Britain. Methods Single‐ and multi‐species distribution models were fitted for 701 native British plant species at a 10‐km grid scale. Two machine learning methods were used – classification and regression trees (CARTs) and artificial neural networks (ANNs). The single‐species versions are widely used in ecology but their multivariate extensions are less well known and have not previously been evaluated against one another. We compared their abilities to predict species distributions, community compositions and species richness in an independent geographical region reserved from model‐fitting. Results The single‐ and multi‐species models performed similarly, although the community models gave slightly poorer predictive accuracy by all measures. However, from the point of view of the whole community they were much simpler than the array of single‐species models, involving orders of magnitude fewer parameters. Multi‐species approaches also left greater residual spatial autocorrelation than the individualistic models and, contrary to expectation, were relatively less accurate for rarer species. However, the fitted multi‐species response curves had lower tendency for pronounced discontinuities that are unlikely to be a feature of realized niche responses. Main conclusions Although community distribution models were slightly less accurate than single‐species models, they offered a highly simplified way of modelling spatial patterns in British plant diversity. Moreover, an advantage of the multi‐species approach was that the modelling of shared environmental responses resolved more realistic response curves. However, there was a slight tendency for community models to predict rare species less accurately, which is potentially disadvantageous for conservation applications. We conclude that multi‐species distribution models may have potential for understanding and predicting the structure of ecological communities, but were slightly inferior to single‐species ensembles for our data.     

Fish distribution predictions from different points of view: comparing associative neural networks,geostatistics and regression models

Accurate prediction of species distributions based on sampling and environmental data is essential for further scientific analysis, such as stock assessment, detection of abundance fluctuation due to climate change or overexploitation, and to underpin management and legislation processes. The evolution of computer science and statistics has allowed the development of sophisticated and well-established modelling techniques as well as a variety of promising innovative approaches for modelling species distribution. The appropriate selection of modelling approach is crucial to the quality of predictions about species distribution. In this study, modelling techniques based on different approaches are compared and evaluated in relation to their predictive performance, utilizing fish density acoustic data. Generalized additive models and mixed models amongst the regression models, associative neural networks (ANNs) and artificial neural networks ensemble amongst the artificial neural networks and ordinary kriging amongst the geostatistical techniques are applied and evaluated. A verification dataset is used for estimating the predictive performance of these models. A combination of outputs from the different models is applied for prediction optimization to exploit the ability of each model to explain certain aspects of variation in species acoustic density. Neural networks and especially ANNs appear to provide more accurate results in fitting the training dataset while generalized additive models appear more flexible in predicting the verification dataset. The efficiency of each technique in relation to certain sampling and output strategies is also discussed.     

Logistic regression models for predicting occurrence of terrestrial molluscs in southern Sweden – importance of environmental data quality and model complexity

We compared the performance of four logistic regression models of different complexity with different environmental data quality, in predicting the occurrence of 49 terrestrial mollusc species in southern Sweden. Performance of models derived from an explanatory data set was evaluated on a confirmatory data set. The overall predictive success of our models (>80% for the three best model approaches), is as good as in other studies, despite the fact that we had to transform a text database into quantitative habitat variables. Simple models (no variable interactions), with forward selection, and detailed habitat data (from field visits) showed the best overall predictive success (mean=84.8%). From comparisons of model approaches, we conclude that data quality (map‐derived data vs habitat mapping) had a stronger impact than model complexity on model performance. However, most of these models showed relatively low values (mean=0.29) for Kappa (statistic for model evaluation), suggesting that the models need to be improved before they would be applied. Predictive success was strongly associated with species incidence but also Kappa was positively correlated with species incidence in univariate tests. Predictive success for true absences was negatively correlated with predictive success for true presences (R²=0.69) and most models failed to give a good prediction of both categories. Models for species with a high incidence in “Open dry sites” or “Mesic interior forests” had a better performance than expected, suggesting that occurrences of species with preference for “narrow” habitats are most easy to predict. Tree layer variables (openness and species abundance) were included in 48 of the 49 final predictive models, suggesting that these variables were good “indicators” of habitat conditions for ground‐living molluscs. Twenty‐four species models included distance to coast and altitude, and we interpret these associations as partly being related to differences in climate. In the final models, true presences (36.9% correctly classified) were much more difficult to predict than true absences (89.7% correct). Possible explanations might be that important habitat variables (e.g. chemical variables and site history) were not included. On the other hand, all suitable sites would not be expected to be occupied due to dynamics in local extinctions (meta‐population theory).     

Genetic inference as a method for modelling occurrence: A viable alternative to visual surveys

Management and conservation require a comprehensive understanding of species distributions and habitat requirements. Reliable species occurrence data are critical in the face of climate change and other anthropogenic activity, but are often difficult to obtain, particularly for wide ranging species. This directly affects ecological models of occurrence and habitat suitability and, in turn, conservation and management decisions. We used generalized linear mixed‐effects models to identify ecological determinants of occurrence for four macropod species (across a region of tropical northern Australia) using a non‐invasive genetic scat approach with and without additional observation records from visual surveys. We show that genetically derived occurrence data, alone, can be used to develop informative ecological models that describe the inter‐specific habitat requirements of macropods. Furthermore, we show that genetic scat surveys of macropods are cheaper and less time consuming to conduct, and tend to provide more occurrence records (and less false absences) than visual surveys. We conclude that indirect surveys using molecular approaches have an important role to play in modelling species' occurrence, and developing future management practices and guidelines to aid species conservation.     

Pale Rock Sparrow Carpospiza brachydactyla in the Mount Lebanon range: modelling breeding habitat

The Pale Rock Sparrow Carpospiza brachydactyla is a poorly known species with very little documentation of its breeding habitat preferences. Using detailed surveys and habitat modelling for a population in southern Lebanon's Bekaa valley, we have identified aspects of the habitat associated both temporally and spatially with breeding. Static spatial habitat modelling using three fundamentally different statistical techniques (multiple linear regression, regression trees and artificial neural networks) agreed on positive associations of breeding density with 'pebbliness' of ground cover and the quantity of available habitat, and a negative association with trees. Temporal associations were found between breeding and a rise in temperature and peaks in grasshopper and beetle abundance, the two main prey that we observed the birds taking. These findings are discussed in the context of Pale Rock Sparrow conservation and implications for species-directed habitat assessment more generally.     

Method selection for species distribution modelling: are temporally or spatially independent evaluations necessary?

To assess the realism of habitat projections in the context of climate change, we conduct independent evaluations of twelve species distribution models, including three novel ecosystem‐based modelling techniques. Habitat hindcasts for 24 western North American tree species were validated against 931 palaeoecological records from 6000, 11000, 14000, 16000 and 21000 yr before present. In addition, we evaluate regional extrapolations based on geographic splits of >55000 sample plots. Receiver operating characteristic analyses indicated excellent predictive accuracy for cross‐validations (median AUC of 0.90) and fair accuracy for independent regional and palaeoecological validations (0.78 and 0.75). Surprisingly, we found little evidence for over‐parameterisation in any method. Also, given high correlations found between model accuracies in non‐independent and independent evaluations, we conclude that non‐independent evaluations are effective model selection tools. Ecosystem‐based modelling approaches performed below average with respect to model sensitivity but excelled in specificity statistics and robustness against extrapolations far beyond training data, suggesting that they are well suited to reconstruct historical biogeographies and glacial refugia.     

Modeling the establishment of invasive species: habitat and biotic interactions influencing the establishment of <Emphasis Type="Italic">Bythotrephes longimanus</Emphasis>

Bythotrephes longimanus is an invasive pelagic crustacean, which first arrived in North America from Europe in early 1980s and can now be found throughout the Great Lakes and in many inland lakes and waterways. Determining the suitability of lakes to Bythotrephes establishment is an important step in quantifying its potential habitat range and environmental risk. Lake environmental conditions, planktivorous fishes, sport fishes and Bythotrephes occurrence data from 179 south-central Ontario lakes were used in this study to model lake characteristics suitable for its establishment. The performance of principal component analysis and different predictive models was used to determine the habitats that are suitable for the survival of Bythotrephes and the factors that may regulate its spread. Four modeling approaches were employed: linear discriminant analysis; multiple logistic regression; random forests; and, artificial neural networks. Ensemble prediction based on the four modeling approaches was also used as an indicator for predicting Bythotrephes occurrence. Bythotrephes appears to establish more readily in larger, deeper lakes with lower elevation, that have more sport fishes. Bythotrephes occurrence can be best predicted by artificial neural networks when including the measures of fish data, in addition to lake environmental data. Lake elevation, surface area and sport fish occurrence were ranked as the most important predictors of Bythotrephes invasion. The inclusion of biotic variables (occurrence or diversity of sport or planktivorous fishes) enhanced cross-validated models relative to analyses based on environmental data alone.     

Can niche‐based distribution models outperform spatial interpolation?

Aim Distribution modelling relates sparse data on species occurrence or abundance to environmental information to predict the population of a species at any point in space. Recently, the importance of spatial autocorrelation in distributions has been recognized. Spatial autocorrelation can be categorized as exogenous (stemming from autocorrelation in the underlying variables) or endogenous (stemming from activities of the organism itself, such as dispersal). Typically, one asks whether spatial models explain additional variability (endogenous) in comparison to a fully specified habitat model. We turned this question around and asked: can habitat models explain additional variation when spatial structure is accounted for in a fully specified spatially explicit model? The aim was to find out to what degree habitat models may be inadvertently capturing spatial structure rather than true explanatory mechanisms. Location We used data from 190 species of the North American Breeding Bird Survey covering the conterminous United States and southern Canada. Methods We built 13 different models on 190 bird species using regression trees. Our habitat‐based models used climate and landcover variables as independent variables. We also used random variables and simulated ranges to validate our results. The two spatially explicit models included only geographical coordinates or a contagion term as independent variables. As another angle on the question of mechanism vs. spatial structure we pitted a model using related bird species as predictors against a model using randomly selected bird species. Results The spatially explicit models outperformed the traditional habitat models and the random predictor species outperformed the related predictor species. In addition, environmental variables produced a substantial R² in predicting artificial ranges. Main conclusions We conclude that many explanatory variables with suitable spatial structure can work well in species distribution models. The predictive power of environmental variables is not necessarily mechanistic, and spatial interpolation can outperform environmental explanatory variables.     

Predicting species distributions: a critical comparison of the most common statistical models using artificial species

Aim To test statistical models used to predict species distributions under different shapes of occurrence–environment relationship. We addressed three questions: (1) Is there a statistical technique that has a consistently higher predictive ability than others for all kinds of relationships? (2) How does species prevalence influence the relative performance of models? (3) When an automated stepwise selection procedure is used, does it improve predictive modelling, and are the relevant variables being selected? Location We used environmental data from a real landscape, the state of California, and simulated species distributions within this landscape. Methods Eighteen artificial species were generated, which varied in their occurrence response to the environmental gradients considered (random, linear, Gaussian, threshold or mixed), in the interaction of those factors (no interaction vs. multiplicative), and on their prevalence (50% vs. 5%). The landscape was then randomly sampled with a large (n = 2000) or small (n = 150) sample size, and the predictive ability of each statistical approach was assessed by comparing the true and predicted distributions using five different indexes of performance (area under the receiver‐operator characteristic curve, Kappa, correlation between true and predictive probability of occurrence, sensitivity and specificity). We compared generalized additive models (GAM) with and without flexible degrees of freedom, logistic regressions (general linear models, GLM) with and without variable selection, classification trees, and the genetic algorithm for rule‐set production (GARP). Results Species with threshold and mixed responses, additive environmental effects, and high prevalence generated better predictions than did other species for all statistical models. In general, GAM outperforms all other strategies, although differences with GLM are usually not significant. The two variable‐selection strategies presented here did not discriminate successfully between truly causal factors and correlated environmental variables. Main conclusions Based on our analyses, we recommend the use of GAM or GLM over classification trees or GARP, and the specification of any suspected interaction terms between predictors. An expert‐based variable selection procedure was preferable to the automated procedures used here. Finally, for low‐prevalence species, variability in model performance is both very high and sample‐dependent. This suggests that distribution models for species with low prevalence can be improved through targeted sampling.     

A new method for dealing with residual spatial autocorrelation in species distribution models

Species distribution modelling (SDM) is a widely used tool and has many applications in ecology and conservation biology. Spatial autocorrelation (SAC), a pattern in which observations are related to one another by their geographic distance, is common in georeferenced ecological data. SAC in the residuals of SDMs violates the ‘independent errors’ assumption required to justify the use of statistical models in modelling species’ distributions. The autologistic modelling approach accounts for SAC by including an additional term (the autocovariate) representing the similarity between the value of the response variable at a location and neighbouring locations. However, autologistic models have been found to introduce bias in the estimation of parameters describing the influence of explanatory variables on habitat occupancy. To address this problem we developed an extension to the autologistic approach by calculating the autocovariate on SAC in residuals (the RAC approach). Performance of the new approach was tested on simulated data with a known spatial structure and on strongly autocorrelated mangrove species’ distribution data collected in northern Australia. The RAC approach was implemented as generalized linear models (GLMs) and boosted regression tree (BRT) models. We found that the BRT models with only environmental explanatory variables can account for some SAC, but applying the standard autologistic or RAC approaches further reduced SAC in model residuals and substantially improved model predictive performance. The RAC approach showed stronger inferential performance than the standard autologistic approach, as parameter estimates were more accurate and statistically significant variables were accurately identified. The new RAC approach presented here has the potential to account for spatial autocorrelation while maintaining strong predictive and inferential performance, and can be implemented across a range of modelling approaches.     

Neural hierarchical models of ecological populations

Neural networks are increasingly being used in science to infer hidden dynamics of natural systems from noisy observations, a task typically handled by hierarchical models in ecology. This article describes a class of hierarchical models parameterised by neural networks – neural hierarchical models. The derivation of such models analogises the relationship between regression and neural networks. A case study is developed for a neural dynamic occupancy model of North American bird populations, trained on millions of detection/non‐detection time series for hundreds of species, providing insights into colonisation and extinction at a continental scale. Flexible models are increasingly needed that scale to large data and represent ecological processes. Neural hierarchical models satisfy this need, providing a bridge between deep learning and ecological modelling that combines the function representation power of neural networks with the inferential capacity of hierarchical models.     

Updating bird species distribution at large spatial scales: applications of habitat modelling to data from long-term monitoring programs

Mapping of species distributions at large spatial scales has been often based on the representation of gathered observations in a general grid atlas framework. More recently, subsampling and subsequent interpolation or habitat spatial modelling techniques have been incorporated in these projects to allow more detailed species mapping. Here, we explore the usefulness of data from long-term monitoring (LTM) projects, primarily aimed at estimating trends in species abundance and collected at shorter time intervals (usually yearly) than atlas data, to develop predictive habitat models. We modelled habitat occupancy for 99 species using a bird LTM program and evaluated the predictive accuracy of these models using independent data from a contemporary and comprehensive breeding bird atlas project from the same region. Habitat models from LTM data using generalized linear modelling were significant for all the species and generally showed a high predictive power, albeit lower than that from atlas models. Sample size and species range size and niche breadth were the most important factors behind variability in model predictive accuracy, whereas the spatial distribution of sampling units at a given sample size had minor effects. Although predictive accuracy of habitat modelling was strongly species dependent, increases in sample size and, secondarily, a better spatial distribution of sampling units should lead to more powerful predictive distribution models. We suggest that data from LTM programs, now established in a large number of countries, has the potential for being a major source of good quality data suitable for the estimation and regularly update of distributions at large spatial scales for a number of species.     

Cost effective prediction of the eutrophication status of lakes and reservoirs

1. Eutrophication is a serious threat in many parts of the world, and identifying the environmental factors that determine the spatial distribution of eutrophicated waterbodies as well as the development of management tools is a challenge. 2. In this study, data from the Ile‐de‐France region were analysed to determine if catchment scale environmental variables could predict concentrations of chlorophyll a (used as a proxy for eutrophication status) of artificial lakes and reservoirs. 3. General additive models (GAM) and random forest models (RF) displayed greater predictive power than generalised linear models, indicating the importance of non‐monotonic relationships. Using RF modelling, very high predictive accuracy was achieved for both continuous and binomial (eutrophic or not) response variables (continuous: R² = 0.715; binomial: kappa = 0.764, 89% of waterbodies were accurately predicted). The better predictive power and robustness of RF versus GAM was attributed to the formers ability to better handle complex interactions between predictors and to account for threshold effects. 4. Our results confirmed the close link between the water quality of lakes and reservoirs and the characteristics of their catchments. Moreover, we also showed that (i) simple (e.g. linear and/or monotonic) relationships between catchment land use and water quality were only found for sub‐regional datasets, and (ii) land use needs to be considered in association with complementary environmental variables (hydromorphological variables) to best assess its impact on water quality.     

Predicting introduction, establishment and potential impacts of smallmouth bass

Aim The introduction of non‐indigenous species has resulted in wide‐ranging ecological and economic impacts. Predictive modelling of the introduction and establishment of non‐indigenous species is imperative to identify areas at high risk of invasion to effectively manage non‐indigenous species and conserve native populations. Smallmouth bass (Micropterus dolomieu), a warm water fish species native to central North America has negatively impacted native fish communities, including cyprinids and salmonid populations, as a result of intentional introductions. We predicted the introduction risk; species establishment based on habitat suitability; identified lakes at high risk of invasion; and finally assessed the consequential impacts on native salmon, trout and cyprinid populations. Location Ontario and British Columbia, Canada. Methods Classification tree and logistic regression models were developed and validated to predict the introduction and establishment of smallmouth bass for thousands of lakes. Results Densely human populated areas and larger lake surface areas successfully identify lakes associated with the introduction of smallmouth bass (introduction model) in British Columbia. Climate, lake morphology and water chemistry variables were the driving environmental parameters to define suitable smallmouth bass habitat (establishment model). A combination of the introduction and establishment model identified 138 lakes that are currently at risk in British Columbia to the introduction and establishment of smallmouth bass. Of these 138 high‐risk lakes, 95% of them contain at least one species of salmon, trout or cyprinid, thereby increasing the potential impact of an invasion by smallmouth bass. Main conclusions Our framework can be applied to other terrestrial and aquatic species to obtain a better understanding of the potential risk posed by a non‐indigenous species to an ecosystem. Furthermore, our methodology can be used to focus management efforts on areas at higher risk (e.g. number of potential releases, more favourable habitats) to control future introductions of non‐indigenous species, thereby conserving native populations.     

Community‐level vs species‐specific approaches to model selection

A topic of particular current interest is community‐level approaches to species distribution modelling (SDM), i.e. approaches that simultaneously analyse distributional data for multiple species. Previous studies have looked at the advantages of community‐level approaches for parameter estimation, but not for model selection – the process of choosing which model (and in particular, which subset of environmental variables) to fit to data. We compared the predictive performance of models using the same modelling method (generalised linear models) but choosing the subset of variables to include in the model either simultaneously across all species (community‐level model selection) or separately for each species (species‐specific model selection). Our results across two large presence/absence tree community datasets were inconclusive as to whether there was an overall difference in predictive performance between models fitted via species‐specific vs community‐level model selection. However, we found some evidence that a community approach was best suited to modelling rare species, and its performance decayed with increasing prevalence. That is, when data were sparse there was more opportunity for gains from “borrowing strength” across species via a community‐level approach. Interestingly, we also found that the community‐level approach tended to work better when the model selection problem was more difficult, and more reliably detected “noise” variables that should be excluded from the model.     

Geographic range determinants of two commercially important marine molluscs

Aim We modelled the spatial abundance patterns of two abalone species (Haliotis rubra Donovan 1808 and H. laevigata Leach 1814) inhabiting inshore rocky reefs to better understand the importance of current sea surface temperature (SST) (among other predictors) and, ultimately, the effect of future climate change, on marine molluscs. Location Southern Australia. Methods We used an ensemble species distribution modelling approach that combined likelihood‐based generalized linear models and boosted regression trees. For each modelling technique, a two‐step procedure was used to predict: (1) the current probability of presence, followed by (2) current abundance conditional on presence. The resulting models were validated using an independent, spatially explicit dataset of abalone abundance patterns in Victoria. Results For both species, the presence of reef was the main driver of abalone occurrence, while SST was the main driver of spatial abundance patterns. Predictive maps at c. 1‐km resolution showed maximal abundance on shallow coastal reefs characterized by mild winter SSTs for both species. Main conclusions Sea surface temperature was a major driver of abundance patterns for both abalone species, and the resulting ensemble models were used to build fine‐resolution predictive range maps (c. 1 km) that incorporate measures of habitat suitability and quality in support of resource management. By integrating this output with structured spatial population models, a more robust understanding of the potential impacts of threatening human processes such as climate change can be established.