©IUCN/SCC Otter Specialist Group
A Review of Otter Distribution Modeling: Approach, Scale, and Metrics
Ana L. Quiñónez C., Todd K. Fuller*, and Timothy O. Randhir
Department of Environmental Conservation, University of Massachusetts,
Amherst, MA 01003, USA
* Corresponding Author e-mail: tkfuller@eco.umass.edu
INTRODUCTION
All of 13 otter species in the world are on the International Union for Conservation of Nature Red List (IUCN Otters Specialist Group, 2013). They often are used as flagship species for conservation and are considered indicators of healthy watershed habitats (Kruuk, 2011; Stevens et al., 2011). Other than the two marine species (Enhydra lutris and Lontra felina), otters live in a variety of freshwater habitats, though some species may also be found along marshes, rocky coasts, and mangroves (Kruuk, 2006). In general, aquatic habitat requirements for otters include rivers that contain deep pools that retain water during the winter and dry seasons, dense vegetation surrounding the river for protection, rivers with sandy banks, and a substantial prey density (Pardini, 1998; Ruiz-Olmo et al., 2001; Kruuk, 2006; Sánchez, 2007). As apex predators, otters have an important role within their local food chain, feeding on fish, crustaceans, amphibians, and even birds, reptiles, and mammals (Pardini, 1998; Kruuk, 2006). Unfortunately, as semi-aquatic top predators, they also are highly vulnerable to habitat degradation, as well as to direct removal for the protection of other species and human livelihoods (e.g., fisheries and domestic fowl) and because of the value of their furs (Kruuk 2006, 2011; Scorpio et al., 2016).
Studying otters can be challenging because they are scarce, elusive, can be nocturnal, sometimes live in difficult-to-access habitats, and have fairly large territories and home ranges (Kruuk, 2011). Because of these factors, data collection can be expensive and labor-intensive (Kruuk, 2011). Consequently, the presence of some otter species is largely unknown over large geographic areas, population declines are nearly impossible to detect (Kruuk, 2011), and many questions about their biology, population sizes, and distribution remain unanswered (Foster-Turley et al., 1990; Kruuk, 2006). Thus, it is very challenging to try to ameliorate threats that affect otter persistence, since factors such as population and range size can affect their vulnerability (Brodie et al., 2013).
In the face of a biodiversity crisis, efforts are increasing to improve our understanding of species declines and conservation efforts (Marcelli and Fusillo, 2009). Understanding the factors that drive the distribution of species is important for their conservation and determining their ecological requirements (Lopes Rheingantz et al., 2014). Scientific studies can answer several questions about potential species distributions and species conservation, but can be prohibitive for mammals with large ranges due to cost and amounts of effort required (Lopes Rheingantz et al., 2014). Species distribution modeling is becoming an essential tool for the management of ecosystems and species conservation, as it gives a geographical perspective that can be used as context for future studies (Barbosa et al., 2003). By studying over time, these models can allow tracking of occurrence patterns and changes in population in order to focus conservation efforts in areas that require it the most.
Many modeling approaches exist that accommodate a varied number of data types to estimate species distributions. Therefore, we classified the modeling approaches found, based on the type of data that was used for the model’s development and the statistical process used. When there is lack of information of otter presence or when studying over a large range, common environmental factors related to the species, in general, can be used. Deductive habitat suitability models do not require otter presence data as they use conceptual knowledge about the species-habitat relationships based on expert opinion, literature, and research (Ali et al., 2010). Presence-only models use incomplete information (presence-only data) to represent the ecological niche of a species from the analysis of several variables to define their distribution across an area (Santiago-Plata, 2013; Gomez et al., 2014; Bieber, 2016). In an occupancy model analysis, presence/absence data is used with different variable combinations to choose the model that best accounts for the probability of an individual occupying a site and being detected in a survey (Bennett, 2014). Finally, phenomenological models used presence/absence data to find relationships or correlation between species presence and a variety of factors, sometimes defining whether these relationships are negative or positive. Due to information gaps and the diverse habitat and resource use among the different otter species, much effort has been invested in modeling otter distribution and discovering more about the correlations between their presence and their surrounding habitats and conditions (Barbosa et al., 2001; Park et al., 2002; Nel and Somers, 2007; Sepúlveda et al., 2009; Ali et al., 2010; Jeffress et al., 2011; Lopes Rheingantz et al., 2014; de Oliveira et al., 2015).
In this paper, we assess models constructed for otter species in order to inform future efforts. We compile information on the different variables used to model the distribution of otters, what information is commonly available and what may be required, what different results can be obtained with different models, and model limitations. This should make future modeling of otter distribution easier and more efficient by providing insights as to what is necessary to make successful and useful models, thus improving our conservation efforts. This will also help disseminate knowledge to non-scientists about important factors that can affect the distribution of otter species.
METHODS
In order to find articles related to otter distribution and the different modeling methods, we searched Google Scholar and Web of Science with a series of combination of the following terms: “otter”, “species distribution models”, “semi-aquatic mammals”, “occupancy models”, “habitat suitability”, “distribution”, and “distribution models”. From the publications found, we made a list of the different methods of analysis and the many variables used in order to have a better understanding of what environmental factors were considered and why.
RESULTS
We found 29 publications containing different modeling methods for assessing otter distributions and correlating distribution with environmental factors. Most of the articles found were related to Lutra lutra, and were from studies across its range in the Tyne catchment (England), Abruzzo region (Italy), southern Italy, Molise region (Italy), Soraksan Nacional Park (Korea), Spanish provinces, the Iberian Peninsula (Spain and Portugal), Hungary, Switzerland, Italy and across Europe. For Lontra canadensis, we identified studies in Maine, the Midwest, New Jersey and Nebraska in the U.S. For Lontra longicaudis, we located studies in Ibera Lake (Argentina), central México, Rio San Juan (Costa Rica), Parana River Delta (Argentina), Pueblo Nuevo (México) and across its geographical range, and for Lontra provocax in Nahuel Huapi National Park (Argentina) and Chile. Studies of Aonyx capensis were from South Africa, Pteronura brasiliensis from the northern Brazilian Amazon, Lutrogale perspiccillatea from the Indus plains of Pakistan, and Enhydra lutris from Glacier Bay, Alaska (see Appendix A1 for Tables A1-A8 that summarize the variables used for each species). Species such as Lontra provocax, Aonyx capensis, Pteronura brasiliensis, and Lutroale perspicillata are relatively underrepresented. These results also highlight unrepresented species, such as Aonyx cinereus, Aonyx congicus, Lontra feline, Hydrictis maculicollis, and Lutra sumatrana.
Among these studies there were four main modeling approaches for defining otter distribution across a given area or the relationship between otter species and certain environmental factors (Table 1). These included deductive habitat suitability models, presence-only models, occupancy model analysis, and presence-absence phenomenological models/variable correlation. The use of these different methods could have been influenced by the availability of modeling techniques at the time, the availability of information, the objectives of the study or analysis, and the extent and/or type of area being used or described (see Appendix A2 for information on advantages and disadvantages for each modeling approach, according to articles reviewed).
Table 1: Model types used to assess otter distribution | ||
Model type | No of References | References |
Presence-only models | 8 | Sepúlveda et al., 2009 Cianfrani et al., 2010; 2011 Santiago-Plata, 2013 Cirelli and Sánchez-Cordero, 2009 Gomez et al., 2014 Lopes Rheingantz et al., 2014 Bieber, 2016 |
Occupancy model analysis | 4 | Santiago-Plata, 2013 Jeffress et al., 2011 Bennett, 2014 Bieber, 2016 |
Deductive Habitat Suitability | 5 | Ottino et al., 1995 Loy et al., 2009 Ottaviani et al., 2009 Ali et al., 2010 Gomez et al., 2014 |
Presence-Absence phenomenological models/variable correlation | 19 | Dubuc et al., 1990 Kemenes and Demeter, 1995 Thom et al., 1998 Barbosa et al. 2001; 2003 Park et al., 2002 Aued et al., 2003 Gori et al., 2003 Nel and Somers, 2007 Sepúlveda et al., 2009 Marcelli and Fusillo, 2009 Cianfrani et al., 2010; 2011; 2013 Gomez et al., 2014 Carone et al., 2014 de Oliveira et al., 2015 Cruz et al., 2017 Williams et al., 2017 |
There were also differences in the scales of the modeling efforts (Table 2), including local, regional, country, continental, and at the geographic extent of the species. The use of different scales was likely due to the information need/gaps for the different species and the objectives of the study or analysis, such as species status in an area or reintroduction efforts (see Appendix A3 for Tables A9-A13 that summarize the metrics used at each scale).
Table 2: Scales used to model otter distributions | ||
Scale | No of References | References |
Local | 9 | Dubuc et al., 1990 Thom et al., 1998 Park et al., 2002 Aued et al., 2003 Gori et al., 2003 Santiago-Plata, 2013 Gomez et al., 2014 Cruz et al., 2017 Williams et al., 2017 |
Regional | 10 | Ottino et al., 1995 Cirelli and Sánchez-Cordero, 2009 Loy et al., 2009 Cianfrani et al., 2010 Ali et al., 2010 de Oliveira et al., 2015 Jeffress et al., 2011 Bennett, 2014 Carone et al., 2014 Bieber, 2016 |
Country | 8 | Kemenes and Demeter, 1995 Barbosa et al. 2001; 2003 Nel and Somers, 2007 Sepúlveda et al., 2009 Marcelli and Fusillo, 2009 Ottaviani et al., 2009 Cianfrani et al., 2013 |
Continent | 1 | Cianfrani et al., 2011 |
Geographic range | 1 | Lopes Rheingantz et al., 2014 |
Anthropogenic Disturbance variables
Anthropogenic disturbance variables are significant contributing factors to negative effects on otter presence and habitat quality, and are the result of human population growth and human behaviors (de Oliveira et al., 2015). Freshwater habitats are highly impacted by anthropogenic activity, such as pollution and water diversion and use, which affect water quality and quantity, and riparian vegetation (Kemenes and Demeter, 1995; Barbosa et al., 2003; Nel and Somers, 2007; Sepúlveda et al., 2009; Cianfrani et al., 2010). Roads have an effect on otters via habitat fragmentation, high sedimentation of watercourses, and increased human disturbance due to greater access to otter habitat (Barbosa et al., 2003).
The variables most commonly used were roads, population density/distribution, and land use. Metrics significantly affecting otter presence were distance to roads (Park et al., 2002), number of visitors to a park (Park et al., 2002), human settlements (Aued et al., 2003), water use and pollution (Nel and Somers, 2007), roads (Sepúlveda et al., 2009), agriculture/livestock adjacent areas (Marcelli and Fusillo, 2009), proportion of urban areas (Marcelli and Fusillo, 2009), distance from industrial areas (Marcelli and Fusillo, 2009), human population density (Lopes Rheingantz et al., 2014; Marcelli and Fusillo, 2009), distance from surface excavations (Cianfrani et al., 2010), distance from productive areas (Cianfrani et al., 2010), proportion of area comprised of cropland (Jeffress et al., 2011), location of fishing nets (de Oliveira et al., 2015), location of homes (de Oliveira et al., 2015), and distance to the nearest otter release site (km)( Bieber, 2016).
Sometimes, using anthropogenic disturbance variables can be complex since their effect on otter presence is variable. Otters have been found to have high resistance to disturbance factors and can be found in areas that we normally consider too disturbed to be ideal for their use (e.g. Kemenes and Demeter, 1995; Bennett, 2014). Human density or development may have a negative effect on otter distribution, but this may vary at a regional scale or with habitat quality (Bennett, 2014). The effect of disturbance factors may be direct, causing otters to avoid certain areas, or indirect through changes inhabitant conditions (Bennett, 2014).
Terrestrial variables
For otters, vegetation could be important as a source of refuge (Gori et al., 2003), as resting and breeding sites, for providing water quality, and for increasing fish productivity (Cianfriani et al., 2010, 2011, 2013; Carone et al., 2014). Riparian vegetation is commonly correlated with high water quality, high primary productivity, high fish biomass, and high availability of alternative prey species (Ottaviani et al., 2009). Altitude/elevation is limiting when considering that there is more food availability at lower and medium river sections than in headwaters (Barbosa et al., 2003). Acclivity (upward slope) may be considered important because very steep river banks have been considered good indicators of areas inaccessible to humans, and as optimal sites for otter holts/dens and couches/resting sites (Ottaviani et al., 2009). Mean altitude has been used as a surrogate of otter habitat quality and its variation (Marcelli and Fusillo, 2009). Slope and topographic convexity are variables that may influence the hunting opportunities (Kruuk, 2006; Cianfrani et al., 2013). Soil permeability is a factor that may affect otter presence in a negative way due to its effect on superficial freshwater availability (Barbosa et al., 2003).
Vegetation and elevation were the variables most commonly used. Metrics found significant for otter presence were percent of forested land composed of birch-aspen (Dubuc et al., 1990), percent of forested land composed of mixed hardwood-softwood (Dubuc et al., 1990), sum of the areas of all water bodies characterized by emergent herbaceous vegetation (Dubuc et al., 1990), density of bank vegetation (Kemenes and Demeter, 1995), soil permeability (Barbosa et al., 2001), mean longitude (Barbosa et al., 2001), coarse scale extra-riparian CORINE (Coordination of Information on the Environment, Land Cover database developed by project of Commission European of the European Union) land cover (Park et al., 2002), vegetation type of stream bank zone (Gori et al., 2003), vegetation complexity (Aued et al., 2003), elevation (Aued et al., 2003), semi-dense riparian vegetation (Sepúlveda et al., 2009), proportion of survey area buffer comprised of woodland (Jeffress et al., 2011), and proportion of survey area comprised of grassland (Jeffress et al., 2011)
Sometimes using terrestrial variables can be complex, since their effect on otter presence is variable. For example, Lopes Rheingantz et al. (2014) describe how elevation was not found to influence their model as it did in other studies, and that vegetation cover had little influence on the model. Variables such as elevation, slope, and density of bank were identified as significant in some studies, but not in others.
Aquatic variables
Water availability is crucial for otters (Cianfrani et al., 2011), a semi-aquatic species that spends a large part of its time in aquatic environments. Water bodies are also a source of fish, which is the most common otter food. Water availability and water quality should have an influence on otter distribution and presence (Kemenes and Demeter, 1995; Cianfrani et al., 2013). Otters appear sensitive to reduction of water depth (Kemenes and Demeter, 1995), as well as stream order and its variation (Marcelli and Fusillo, 2009). Hierarchy of tributaries is used as a proxy of water flow (Ottaviani et al., 2009). If otters are forced to find food sources out of the water or too close to the shore, they may become more vulnerable to terrestrial predators, in turn affecting their survival rates (Ruiz-Olmo and Jimenez, 2009).
In modeling efforts, water body characteristics and river hierarchy were the variables most commonly used. Metrics found significant for otter presence included mean shoreline diversity index (shape; Dubuc et al., 1990), total stream length, over all stream orders (Dubuc et al., 1990), water depth (Kemenes and Demeter, 1995; Nel and Somers, 2007), river/stream width (Park et al., 2002), bottom structure of stream (Park et al., 2002), bank (shore) type (Gori et al., 2003), current type (Nel and Somers, 2007), anastomosed (two or more interconnected channels that enclose flood basins) river length (Sepúlveda et al., 2009), sum of the waterbody perimeters/sum of waterbody areas for entire watershed (Jeffress et al. 2011), sum of stream (3rd order) km within the watershed/watershed area (Jeffress et al., 2011), number of waterbodies within the watershed/watershed area (Jeffress et al., 2011), river water level (de Oliveira et al., 2015), long-term median flow rate of the river (ft3/s) (Bieber, 2016), flow zone (Cruz et al., 2017), total dissolved solids (Cruz et al., 2017) and pH (Cruz et al., 2017).
Sometimes using aquatic variables can be a challenge, since the characteristics that need to be measured for their effect on otter presence are not easily defined. Water fluctuations can have a negative effect on fish abundance and size; therefore, floods and droughts can cause otters to abandon areas (Ruiz-Olmo et al., 2001), though this reaction seems to vary from species to species, and from area to area. In some areas, species such as Lutra lutra are able to live in dry rivers during the summer, as long as there are pools that provide enough fish to eat throughout the season (Ruiz-Olmo et al., 2001; Prenda et al., 2001). For other species, like Aonyx capensis in South Africa, freshwater availability is more important than prey availability (Van Niekerk et al., 1998). Also, variability in the metrics presented shows how difficult it is to define which characteristics of a water body can affect otters. Sometimes the scale of the study could be what affects the effect of the metrics; river/stream width was considered a significant variable by Park et al. (2002; local scale), but in Nel and Somers (2007) it was not (country scale) (Appendix A3, Tables A9-A13).
Climate variables
Climate mostly influences distributions of species at macroscales (Cianfrani et al., 2011). Climate factors at large scales have high potential as surrogates for local freshwater availability, and water warming could affect fish species diversity and abundance (Cianfrani et al., 2011). Floods can increase the deposit of suspended solids, which tend to bury potential denning areas, as well as decreasing food availability for fish and otters (Ruiz-Olmo et al., 2001). Droughts may also increase mortality, because with a drought comes diminishing food availability which may trigger an increase in territoriality among individuals (Prenda et al., 2001).
In the models we looked at, temperature and precipitation where the climate variables most commonly used. Significant metrics affecting otter presence included relative humidity in January (Barbosa et al., 2001) and annual temperature (Lopes Rheingantz et al., 2014). Lopes Rheingantz et al. (2014) also found that annual precipitation was the most relevant climatic metric for neotropical otter (Lontra longicaudis) distribution within its geographic range.
Using climatic variables can have its difficulties, such as finding the information at an appropriate scale for the study. The common use of macroscale global climatic data, even in local studies, is a clear example of this. The effect of climatic variables on otters has not been directly measured, so many studies use variables that are assumed to be most significant to the species and/or were used in previous mammal studies (e.g., Lopes Rheingantz et al., 2014).
Biological interaction variables
Food availability has been found to be the factor of most importance for otter presence (Kruuk, 2006; Nel and Somers, 2007; Cianfrani et al., 2013). Mink (Mustela vison) are considered a potential competitor for resources (Aued et al., 2003), and beaver (Castor canadensis) presence has been found to be a predictor of otter presence (Dubuc et al., 1990; Bennett, 2014). The metrics found significant for otter presence were: percent of all wetlands with active or inactive beaver sign (Dubuc et al., 1990), food availability (Nel and Somers, 2007), and freshwater crab and crayfish abundance (Sepúlveda et al., 2009).
Disadvantages of using biological interaction variables include how costly and time consuming it is to define food availability, which is considered as the most important factor affecting otter distribution. In many cases, it is easier to use surrogate variables for food availability rather than measure actual food presence and abundance, but this could cause overprediction, may not be as accurate, and makes it riskier to interpret the data (Sepulveda et al., 2009).
DISCUSSION
The reliability of modeling efforts may be influenced by the modeling scale (i.e., ecological scale or the extent of the landscape under consideration) used, measuring issues, and several other shortcomings that we are forced to face as we work with such complex ecosystems (cf. Elith and Graham 2009). Some of the complications found or mentioned in the reviewed articles are as follows:
Influence of modeling scale
The importance of scale is very often underestimated or not accounted for in many ecological studies (Thom et al., 1998). Using course resolution in an analysis can make it complicated to assess land use and connectivity during the analysis (Cianfrani et al., 2011). In addition, scale mismatch between data and ecological process is a major problem in ecological modeling. Fine-scale data is important for modeling some species, including otters, because characteristics such as riparian vegetation cover may not be well represented in coarser data layer such as land cover (Loy et al., 2009). Habitat variables that might be effective to predict species response at one scale might not be as effective other scales (Ali et al., 2010); unfortunately, the resolution of variables usually depends on data available. HIS models are used under the assumption that habitat-wildlife relationships are consistent at all scales (Ali et al., 2010). Environmental data for freshwater bodies (water temp, depth, water velocity, etc.) is usually not spatially accurate to be used in models (Cianfrani et al., 2013).
Otter habitat is complex, consisting of a narrow strip of an aquatic and riparian ecosystem, and though individuals may move several hundred meters from this area, their activity mostly occurs close to this strip (Ruiz-Olmo et al., 1998; Ottaviani et al., 2009). Therefore, fine-scale modeling is appropriate to measure decreasing habitat suitability as one considers habitats away from a waterway, as well as the effects of land use as it moves towards riparian habitat (Ottaviani et al., 2009). It is not easy to obtain useful large-scale information of habitat suitability based on fine scales habitat linearity (Ottaviani et al., 2009). Relevant information such as fish abundances, water flow, hunting pressure, and water pollution are rarely available or reliable at large spatial scales, particularly since they fluctuate a lot within time and space. Also, survey techniques are difficult to standardize at large scales for some species in some systems; therefore, proxies are commonly used (Ottaviani et al., 2009; Cianfrani et al., 2013). Large-scale efforts need to be refined with local data such as pollution, food availability, and human disturbance when you want to apply them in local conservation plans, or else only be used for large scale conservation strategies (Ottaviani et al., 2009).
Fine-scale models are limited in their application and evaluation of potential habitat at a larger scale (Park et al., 2002). Sometimes one has to use a larger, less accurate scale to identify the fine scale of microhabitat (Lopes Rheingantz et al., 2014); this is because the fine-scale data are usually not available over large areas.
Measurement issues
Sparse otter occurrence may lead to overestimated range sizes when including areas of sporadic occurrences (Marcelli and Fusillo, 2009). Some areas where otters are observed are used only to move from one site to a better one, and they do not necessarily represent an area that the otter regularly inhabits. Using the characteristics of these rarely used areas as basis for defining habitat suitability can cause for overestimating the areas that are used by the species. Also, insufficient numbers of data points may not provide enough information to use independently for validation and calibration, therefore performance could be overestimated (Sepúlveda et al., 2009). In general, overprediction and underprediction can affect results for distribution ranges and for conservation efforts; therefore, making different types of models and then overlaying them may be a good way of reducing this effect, but there are definitely pitfalls to avoid (Cade, 2015). Optimization methods in modeling fitting and testing for such deviations in predictive ability should be useful way forward in model selection (e.g., Elith and Leathwick, 2009; Merow et al., 2014)
Sometimes it is difficult to measure the time in which a track was made or when spraints were deposited, usually due to their location; therefore, there may be discrepancies with the actual environmental conditions in which they were made (Kermenes and Demeter, 1995). Weather (snow and rain) and water level variations may affect the ability to detect indirect signs of otter presence, and therefore affect the results of our models, so monitoring should be done consistently during drier seasons (de Oliveira et al., 2014; Bieber, 2016).
The use of more detailed variables will allow a better understanding of actual relationships between otter and their surrounding habitat (species or types of vegetation for example). Fish population estimates at each site might not be representative of the population of the full stretch of habitat/river (Thom et al., 1998), though obtaining a more representative estimate might be impractical due to financial and time constraints.
In most of the work related to absence/presence data there is a possibility of being biased due to “false absences”. It is difficult to know a priori which absences are reliable and which ones are not since species distribution is usually a snapshot in time of a system that is dynamic (Cianfrani et al., 2010). Sometimes the species may be considered as absent in an area, but in reality, it was just not detected (Jeffress et al., 2011; Ruiz-Olmo et al., 2001). There is also always a possibility of errors when using multiple/inexperienced observers.
Sometimes when using data collected over time, methods might not be comparable between the surveys (e.g. difference in grid system or lack there off, surveying one or both bank sides; Marcelli and Fusillo, 2009). Information (scat census or density and vegetation cover) used may not always be from the same time periods, and may have changed since the time the information was produced (outdated), affecting the relationships found within the model and therefore the accuracy the results. Usually when using data from different time frames and projects, the geographic extent of the surveys is not the same (Marcelli and Fusillo, 2009).
The information available for modeling efforts can come from different sources such as scientific collections, museums, herbariums and online databases (Table 4). This information might have several issues: there may only be presence data, the species might not have been classified, the data might not be correctly georeferenced, and data are usually collected for different reasons, without a standardized methodology, therefore representing a biased distribution of the species (Santiago-Plata, 2014).
Shortcomings
It is common to lack enough data points for the range being modeled, or for there to be areas that are unrepresented by available data. When working with data from different time frames, usually one of the observation times has a limited dataset (Carone et al., 2014). Fish density (when available) is usually not measured evenly within the otter’s range.
The level of productivity of rivers varies, being usually low in the headwater, increasing in the middle reaches and peaking in the lower reaches (Nel and Somers, 2007). Therefore, considering all of a river as suitable or with the same suitability is not adequate or realistic. Temporal, spatial, or quantitative variation in negative or positive effects of factors on the species may be case specific (Marcelli and Fusillo, 2009). Lack of observed influence of some variables could be due to their low variability through study area (Jeffress et al., 2011, Cruz et al., 2017). Using different climate scenario for future predictions can cause discrepancies in the models (Cianfrani et al., 2011).
Spraint numbers are not sensitive to changes in otter distribution in relation to changes in prey distributions (Thom et al., 1998). Location of spraint and spraint sites might sometimes interfere with their relationship to the periods of feeding activity, and in some sites, they might last longer (Thom et al., 1998).
Sometimes the lack of information on local water fluctuations will force one to downgrade the suitability of certain areas such as smaller streams and areas located at certain altitudes (Ottaviani et al., 2009). It is common to ignore water regimes of different watercourses in the modeling process, even though it affects their carrying capacity at several levels (Ottaviani et al., 2009).
It is not a simple task to determine if you have been able to choose all relevant factors that affect otter presence and distribution (Barbosa et al., 2003), particularly because of surrogate variables and the extent to which they may correlate to other variables used. Using regional bioclimatic variables (representing annual tendencies, stationarity, and extreme factors) without taking into account that some of them might not have a relationship with the species being modeled may generate instability in the models generated (Santiago-Plata, 2014). Selection of variables may lead to randomness in the predictions (Cianfrani et al., 2010). Reliability of models also depends on the species ability to adaptation and the environment’s temporal and spatial variability (Cianfrani et al., 2010). Pseudo-absence data will affect distribution modeling efforts, no matter how minimized their effect is (Carone et al., 2014).
Additional considerations
Sometimes, as in the case of Aued et al. (2003), the lack of otter presence in an area within the study cannot be defined quantitatively, therefore you can infer what factors can qualitatively be causing this absence. These inferences cannot be statistically demonstrated, but perhaps further research will provide the information needed.
There are also cases in which there are too many variables to consider or factors that cannot be accounted for directly, and in order to make the analysis less complex or to include other possibilities, a single variable of “other” is used. This is a complicated decision to make, since the effect of these compound variables is not being clearly defined, and since any of the many options included could be responsible for it. Also, the extent to which each variable is responsible for affecting the otter species is difficult, if not impossible, to measure. This highlights the need for more careful attention to model choice and model fitting, for which there is a vast literature (e.g., Burnham and Anderson, 2003; Johnson and Omland, 2004; Guisan et al., 2017).
Many of the variables that are important in determining otter abundance are affected by climate change. Temperature and precipitation ranges and distribution are the main factors affected directly by climate change. Though their direct effect on species is usually unknown, their effect through freezing, drought and flooding could have negative consequences on otter populations. These two variables have a huge impact on other variables such as vegetation, food availability, and waterbody-related variables. Vegetation community assemblages are expected to vary due to climate change (Brodie et al., 2013), and vegetation is considered an important variable because it is commonly used as a proxy for refuge and food resources.
Though climate change studies are more common for terrestrial species, we found one paper related to climate change and otters (Cianfrani et al., 2011). The main focus of this research was identifying the effects of temperature and precipitation on otter distribution. Precipitation is expected to have an important role in water availability and distribution. Temperature on the other hand, is expected to affect fish assemblages as water warms up. Results indicated that climate change may cause a profound reshuffling in the potential otter distribution across Europe, though there was some variation in outcomes across the range. Even when vulnerability to climate change and conservation status seem to be correlated, their relationship is not perfect, as other factors may affect their degree of correlation (Brodie et al., 2013).
CONCLUSIONS AND RECOMENDATIONS
As apex predators, otters have an important role in ecosystems, but their dependence on aquatic habitats makes them vulnerable. This dependency on water sources and the food and shelter they provide makes freshwater species more vulnerable, with higher extinction rates, than terrestrial species (Scorpio et al., 2016). Otters are considered among the most threatened mammals in the world (Kruuk, 2006; Scorpio et al., 2016) and prey and water availability seem to be the two most important factors that limit otters (Prenda and Granado-Lorencio, 1996; Prenda et al., 2001; Ruiz-Olmo et al., 2001).
Otters are difficult to study, and thus large information gaps exist for many of the species (see Kruuk, 2011 for species-specific research recommendations based on information gaps). Defining otter distribution, even with the advances in software and technology today, is a complicated process given the gaps in important information. There are several factors that have been deemed important for otter distribution, based on the studies reviewed, such as anthropogenic disturbance, climatic, terrestrial, aquatic and biological interaction variables. Unfortunately, we are still uncertain on how many of the factors are directly or indirectly affecting otters or whether they have the adaptability to deal with changes within them. The relationship between otters and some of these factors are simple to interpret, but other relationships are still unclear. Sometimes a variable deemed as important in one study is considered unimportant in another. Even when a variable is known to be important, it may be hard to measure. Most research had focused on Lutra lutra, though based on the differences between their life histories and geographic ranges with other otters, the results can often not be extrapolated to predict impacts on other species. Even if the results of species-specific studies cannot often be extrapolated, the processes and data types used, the data gathering methods, and the variables/metrics considered can still provide guidance for future research. Researchers should still be careful in considering what really applies to their study areas and species of interest.
We still need to have a better understanding of the relationships among many factors and otter distribution, as well as their varying effects on different otter species, and articles cited in this review provide many useful suggestions for improved modeling. Modeling should be considered a dynamic process in order to progressively improve the quality of the predictions, and adequate evaluation indexes should be used when evaluating model quality (Cianfrani et al., 2010). Robust spatially explicit models for identifying and hierarchically assessing areas for otter conservation and restoration can be achieved with sequential implementation of methods combining species modeling and place prioritization (Cirelli and Sánchez-Cordero, 2009). The use of multiple survey methods, data sets and analysis methods to allow a better representation of the areas of interest and the direct comparison between the methods being used (Bieber 2016).
With respect to the type of model being used, accuracy of presence data is important for calibrating habitat suitability models (Cianfrani et al., 2013), as is definition of the cut-off point above which the presence of a species is more likely than expected at random, since it can be used to correct the established thresholds that are used to separate unsuitable from suitable areas (Ali et al., 2010). When the species-environment equilibrium assumption (this assumption presumes that a species occupies all suitable habitat that is available) is not met (e.g. recolonization and expansion), habitat suitability models’ predictions should be assessed carefully (Cianfrani et al., 2010). Prior to employing environmental niche models, an important step is to test for environmental similarity (Cianfrani et al., 2013). When dealing with species with unstable spatial equilibrium, presence-only models may be a better option than presence–absence methods for making reliable predictions of suitable areas for expansion (Cianfrani et al., 2010).
Using the same variables at different scales may have different effects on populations, and therefore should be analyzed the most appropriate scale (Thom et al., 1998). Unfortunately, there are many variables whose information is not usually available at different scales. Studies of regional-scales processes are an important complement for local-scale studies, providing a broader geographic perspective that can be seen as context in local studies, and allowing us to take into account factors that have an effect on a larger scale (Barbosa et al. 2003).
Better research on direct factors affecting otter distribution is needed, since proxy variables, and even seemingly direct habitat characteristics, sometimes indicate the opposite of expectations (Kemenes and Demeter, 1995). In considering the shifting dynamics over time and how factors effect on species distributions may vary (Sépulveda et al., 2009), radio-telemetry may provide the best data for the analysis of habitat choice and use by otters (Nel and Somers, 2007). Despite being difficult to quantify, additional emphasis should be placed on water quality and prey availability, given their importance. Finding better ways of integrating these factors into analyses will allow results to be more reliable. Distribution studies could be directed towards areas where there is previous information available regarding other important factors, such as water quality assessment (Bennett, 2014) and food availability, in order to help include this critical but hard to obtain information within the studies. PCB contamination and fish density are necessary in a spatially explicit way, and therefore should be a priority for future research (Cianfrani et al., 2013). When trying to build models on the effect of environmental variables, the results from field surveys should always be compared with water quality data (Kemenes and Demeter, 1995).
Also, it is important to collect long-term site occupancy data and to use modeling procedures that account for imperfection in detectability (Marcelli and Fusillo, 2009). Avoiding the concentration of data points in a particular area over another (i.e., having an equal distribution of data along the area of study) is needed to prevent bias in the probability of distribution of the models (Santiago-Plata, 2014). It also is important to obtain data on trends over time from periodic survey of factors such as water quality, land use, anthropogenic disturbance changes, and vegetation.
We need to have a better understanding of the relationship between otters and the variables that are commonly used to describe their habitats, especially water quality and prey availability, as these seem to be the most important but the hardest to quantify. Ultimately, to gain a more accurate and meaningful understanding of otter survival, we need to focus on finding better ways of integrating measurable variables into our analyses.
After considering the modeling approaches used in the papers reviewed, a few suggestions come to mind. Williams et al. (2017) applied a Bayesian approach for the first time in these modeling efforts while calculating sea otter occupancy and abundance; it would be interesting to try these methods on other otter species, as well as other methods such as simulations and mission learning tools. Advances in technology and science open the door for the use of new tools for our conservation efforts; thus, we need to keep an eye on these advances and an open mind to distribution studies made for other species.
Appendix A1 – Variables used for different otter species
Appendix A2 - Advantages and disadvantages of modeling approaches
Appendix A3 - Appendix A3 – Metrics used at each scale
Acknowledgements: The authors would like to thank Brian Gerber for very helpful review comments and insights.
REFERENCES
Ali, H., Saleem, R., Qamer, F. M., Khan, W. A., Abbas, S., Gunasekara, K., Hazarika, M., Ahmed, M. S., Akhtar, M. (2010). Habitat evaluation of smooth-coated otter (Lutrogale perspicillata) in Indus plains of Pakistan using remote sensing and GIS. Remote Sensing and Spatial Information Science 38: 127-132.
Aued, M., Chehebar, C., Porro, G., Macdonald, D., Cassini M. (2003). Environmental correlates of the distribution of southern river otters Lontra provocax at different ecological scales. Oryx 37: 413.
Barbosa, A. M., Real, R., Marquez, A. L., Rendon, M. A. (2001). Spatial, environmental and human influences on the distribution of otter (Lutra lutra) in the Spanish provinces. Divers. Distrib. 7: 137-144.
Barbosa, A. M., Real, R., Olivero, J., Vargas, J. M. (2003). Otter (Lutra lutra) distribution modeling at two resolution scales suited to conservation planning in the Iberian Peninsula. Biol. Conserv. 114: 377-387.
Bennett, C. H., III (2014). A predictive habitat occupancy model for North American river otter (Lontra canadensis) along low order streams in inland New Jersey. M.S. Thesis, University of Delaware, Delaware, USA.
Bieber, N. R. (2016). River Otter (Lontra canadensis) Distribution and Habitat Suitability in Nebraska.University of Nebraska - Lincoln, Nebraska, USA.
Brodie, J. F., Post, E. S., Doak, D. F. (2013). Wildlife conservation in a changing climate. The University of Chicago Press, Chicago.
Burnham, K.P., Anderson, D.R. (2003). Model selection and multimodel inference: a practical information-theoretic approach (2nd ed.). Springer-Verlag, New York, U.S.A.
Cade, B.S. (2015). Model averaging and muddled multimodel inferences. Ecology 96: 2370-2382.
Carone, M., Guisan, A., Cianfrani, C., Simoniello, T., Loy, A., Carranza, M. (2014). A multi-temporal approach to model endangered species distribution in Europe. The case of the Eurasian otter in Italy. Ecol. Mode.l 274: 21–28.
Cianfrani, C., Le Lay, G., Hirzel, A. H., Loy, A. (2010). Do habitat suitability models reliably predict the recovery areas of threatened species? J. Appl. Ecol. 47: 421-430.
Cianfrani, C., Le Lay, G., Maiorano, L., Satizábal, H. F., Loy, A., Guisan, A. (2011). Adapting global conservation strategies to climate change at the European scale: the otter as a flagship species. Biol. Conserv .144: 2068-2080.
Cianfrani, C., Maiorano, L., Loy, A., Kranz, A., Lehmann, A., Maggini, R., Guisan, A. (2013). There and back again? Combining habitat suitability modelling and connectivity analyses to assess a potential return of the otter to Switzerland. Anim. Conserv. 16: 584-594.
Cirelli, V., Sánchez-Cordero V. (2009). Selection of restoration and conservation areas using species ecological niche modeling: A case study of the neotropical river otter Lontra longicaudis annectens in central México. Universidad Nacional Autónoma de México, México D.F., México.
Cruz García, F., Contreras Balderas, A. J., Nava Castillo, R., Gallo Reynoso, J. P. (2017). Habitat and abundance of the Neotropical otter (Lontra longicaudis annectens) in Pueblo Nuevo, Durango, Mexico. Therya 8: 123-130.
de Oliveira, I. A. P., Norris, D., Michalski, F. (2015). Anthropogenic and seasonal determinants of giant otter sightings along waterways in the northern Brazilian Amazon. Mamm. Biol. – Z. Saugetierkd 80: 39-46.
Dubuc, L. J., Krohn ,W. B., Owen Jr, R. B. (1990). Predicting occurrence of river otters by habitat on Mount Desert Island, Maine. J. Wildlife. Manage. 54: 594-599.
Elith, J., Graham, C.H. (2009). Do they? How do they? WHY do they differ? On finding reasons for differing performances of species distribution models. Ecography 32: 66–77
Elith, J., Leathwick, J.R. (2009). Species Distribution Models: Ecological explanation and prediction across space and time. Annu. Rev. Ecol. Evol. Syst. 40: 677–97.
Foster-Turley, P., Macdonald, S. M., and Mason, C. F. (1990). Otters: an action plan for their conservation. IUCN Gland, Switzerland.
Gomez, J. J., Túnez, J. I., Fracassi, N., Cassini, M. H. (2014). Habitat suitability and anthropogenic correlates of Neotropical river otter (Lontra longicaudis) distribution. J. Mammal. 95: 824-833.
Gori, M., Carpaneto, G. M., Ottino, P. (2003). Spatial distribution and diet of the Neotropical otter Lontra longicaudis in the Ibera Lake (northern Argentina). Acta Theriol. 48: 495-504.
Guisan, A., Thuiller, W., Zimmermann, N.E. (2017). Habitat suitability and distribution models: with applications in R. Cambridge University Press, Cambridge, U.K.
IUCN Otter Specialist Group (2015). The 13 species of otters. <http://www.otterspecialistgroup.org/Species.html>. Accessed 11 January 2018.
Jeffress, M. R., Paukert, C. P., Whittier, J. B., Sandercock, B. K., Gipson, P. S. (2011). Scale-dependent factors affecting North American river otter distribution in the Midwest. The American Midland Naturalist 166: 177-193.
Johnson, J.B., Omland, K.S. (2004). Model selection in ecology and evolution. Trends in Ecol Evol 19: 101-108.
Kemenes, I., Demeter, A. (1995). A predictive model of the effect of environmental factors on the occurrence of otters (Lutra lutra L.) in Hungary. Hystrix 7: 209-218.
Kruuk, H. (2006). Otters: ecology, behaviour, and conservation. Oxford University Press, New York.
Kruuk, H. (2011). Ecological research and conservation management of otters
<www.otterspecialistgroup.org/Library/OSG_Research_Guidelines.pdf> Accessed 11 January 2018.
Lopes Rheingantz, M., Saraiva de Menezes, J. F., de Thoisy, B. (2014). Defining Neotropical otter Lontra longicaudis distribution, conservation priorities and ecological frontiers. Trop. Conserv. Sci. 7: 214-229.
Loy, A., Carranza, M. L., Cianfrani, C., D'Alessandro, E., Bonesi, L., Di Marzio, P., Minotti, M., Regiani, G. (2009). Otter Lutra lutra population expansion: assessing habitat suitability and connectivity in southern Italy. Folia Zool. 58: 309-326
Marcelli, M., Fusillo, R. (2009). Assessing range re-expansion and recolonization of human-impacted landscapes by threatened species: a case study of the otter (Lutra lutra) in Italy. Biodivers. Conserv. 18: 2941-2959.
Merow, C., Smith, M.J., Edwards Jr, T.C., Guisan, A., McMahon, S.M., Normand, S., Thuiller, W., Wüest, R.O., Zimmermann, N.E., Elith, J. (2014). What do we gain from simplicity versus complexity in species distribution models? Ecography 37: 1267–1281.
Nel, J. A., Somers, M. J. (2007). Distribution and habitat choice of Cape clawless otters, in South Africa. S. Afr. J. Wildl. Res. 37: 61-70.
Ottaviani, D., Panzacchi, M., Lasinio, G. J., Genovesi, P., Boitani, L. (2009). Modelling semi-aquatic vertebrates’ distribution at the drainage basin scale: The case of the otter Lutra lutra in Italy. Ecol. Model. 220: 111-121.
Ottino, P., Prigioni, C., Taglianti, A. V. (1995). Habitat suitability for the otter (Lutra lutra) of some rivers of Abruzzo Region (Central Italy). Hystrix 7: 265-268.
Pardini, R. (1998). Feeding ecology of the neotropical river otter Lontra longicaudis in an Atlantic Forest stream, south-eastern Brazil. J. Zool .245: 385-391.
Park, C., Joo, W., Seo C. (2002). Eurasian otter (Lutra lutra) habitat suitability modeling using GIS; A case study on Soraksan National Park. Journal of Korea Spatial Information Society 10: 501-513.
Prenda, J., Granado-Lorencio, C. (1996). The relative influence of riparian habitat structure and fish availability on otter Lutra lutra L. Sprainting activity in a small Mediterranean catchment. Biol. Conserv. 76: 9–15.
Prenda, J., López Nieves, P., Bravo, R. (2001). Conservation of otter (Lutra lutra) in a Mediterranean area: The importance of habitat quality and temporal variation in water Aquat. Conserv. 11: 343-355.
Ruiz-Olmo, J. (1998). Influence of altitude on the distribution, abundance and ecology of the otter (Lutra lutra). In behaviour and ecology of riparian mammals. Dunstone, N., Gorman, M.L. (Eds). Sym. Zool. S. 71: 159–176.
Ruiz-Olmo, J., López-Mart, J. M., Palazón, S. (2001). The influence of fish abundance on the otter (Lutra lutra) populations in Iberian Mediterranean habitats. J. Zoo.l, 254: 325–336.
Ruiz-Olmo, J., Jiménez J. (2009). Diet diversity and breeding of top predators are determined by habitat stability and structure: A case study with the Eurasian otter (Lutra lutra L.). Eur. J. Wildlife. Res. 55: 133–144.
Sánchez, O. (2007). Método de evaluación del riesgo de extinción de las especies silvestres en México (MER): Secretaria de Medio Ambiente y Recursos Naturales: Instituto Nacional de Ecología: Instituto de Ecología de la Universidad Nacional Autónoma de México: Comisión Nacional para el Conocimiento y Uso de la Biodiversidad, México.
Santiago-Plata, V.M. (2013). Ocupación y distribución potencial de la nutria neotropical (Lontra longicaudis) asociada a variables ambientales en la cuenca del río San Juan, Costa Rica. Centro Agronómico Tropical de Investigación y Enseñanza, Turrialba, Costa Rica. [In Spanish]
Scorpio, V., Loy, A., Di Febbraro, M., Rizzo, A., Aucelli P. (2016). Hydromorphology Meets Mammal Ecology: River Morphological Quality, Recent Channel Adjustments and Otter Resilience. River Res Appl 32: 267-279.
Sepúlveda, M., Bartheld, J., Meynard, C., Benavides, M., Astorga, C., Parra, D., and Medina‐Vogel, G. (2009). Landscape features and crustacean prey as predictors of the Southern river otter distribution in Chile. Anim. Conserv. 12: 522-530.
Stevens, S., Organ J. F., Serfass T. L. (2011). Otters as flagships: social and cultural considerations. Proceedings of Xth International Otter Colloquium. IUCN Otter Spec. Group Bull. 28: 150–161.
Thom, T., Thomas, C., Dunstone, N., Evans, P. (1998). The relationship between riverbank habitat and prey availability and the distribution of otter (Lutra lutra) signs: An analysis using a geographical information system. Sym Zool S. 71: 135-157.
Van Niekerk, C. H., Somers, M. J., Nel, J. A. (1998). Freshwater availability and distribution of Cape clawless otter spraints and resting places along the south-west coast of South Africa. S. Afr. J. Wildl. Res., 28: 68–72
Williams, P. J., Hooten, M. B., Womble, J. N., Esslinger, G. G., Bower, M. R, Hefley, T. J. (2017). An integrated data model to estimate spatiotemporal occupancy, abundance, and colonization dynamics. Ecology 98: 328-336
Synthèse Bibliographique des Modèles de la Distribution de la Loutre: Démarche, Échelle et Métriques
Toutes les espèces de loutre sont concernées par la conservation et sont utilisées comme espèce parapluie pour la conservation et comme indicateur de la santé des bassins versants; En conséquence, l’identification et la compréhension de leur distribution est une nécessité élémentaire. Nous avons passé en revue la littérature publiée afin d’identifier un modèle de distribution de la loutre dans le monde et la compilation de données sur différentes métriques/variables utilisées: quelle information est communément disponible et peut être requise? Quels résultats distincts peuvent être obtenus avec différents modèles et limitations de modèle? Nous avons identifié 29 études sur 8 espèces qui utilisaient principalement 4 méthodes de modélisation de distribution de la loutre à travers une région déterminée ou la relation entre les espèces de loutre et certains facteurs environnementaux. Les études ont été modélisées sur base d’une distribution à différentes échelles, incluant l’extension géographique locale, régionale, nationale et continentale. Nous avons catalogué 301 métriques environnementales différentes utilisées dans les modélisations que nous avons ensuite classées en 6 catégories principales : perturbation anthropogénique, climat, interactions terrestre, aquatique et biologique. L’alimentation, la disponibilité en eau et sa qualité, ainsi que les influences anthropogéniques sont régulièrement identifiées comme variables importantes corrélées avec la distribution de la loutre, mais sont souvent mesurées de différentes façons, ou identifiées dans des modèles par des variables indirectes ou de substitution parce que la disponibilité des données pertinentes est insuffisante ou absente. L’échelle, la démarche, et la sélection des métriques, tout cela demande à être envisagé avec précaution pour chaque étude. Cependant, la compréhension des problèmes de dimensionnement et des lacunes de modélisation identifiées par d’autres devraient permettre d’améliorer ce type de modélisation dans le futur. L’examen du contenu de cet article de synthèse bibliographique peut fournir des indications sur les efforts ultérieurs dans les processus de modélisation, les types de données utilisées, les méthodes de collecte des informations, et des variables/métriques à inclure. Cette information doit encore être évaluée avec précaution pour une utilisation sur des zones d’études spécifiques, des espèces dignes d’intérêt, et comme base de développement de méthodes innovantes et plus efficaces.
Revenez au dessus
Resumen: Review del Modelado de Distribución en Nutrias: Enfoque, Escala y Métrica
Todas las especies de nutrias son de preocupación de conservación, y son utilizadas tanto como especies-bandera para la conservación, así como indicadores de salud de cuencas; consecuentemente, identificar y entender su distribución es una necesidad básica. Revisamos la literatura publicada para identificar esfuerzos de modelado de distribución, en todo el mundo, y luego compilamos información sobre distintas métricas/variables usadas, qué información está comúnmente disponible y qué se puede requerir, qué distintos resultados se pueden obtener con diferentes modelos, y limitaciones de los modelos. Identificamos 29 estudios sobre 8 especies que utilizaron 4 métodos principales para modelar distribución de nutrias en una determinada área, o la relación entre las especies de nutria y ciertos factores ambientales. Los estudios modelaron distribución a través de una variedad de escalas, incluyendo la local, regional, nacional, continental, y de toda la distribución de la especie. Catalogamos 301 diferentes métricas ambientales usadas en los modelos, que luego clasificamos en seis categorías principales: disturbios antropogénicos, clima, interacciones terrestres, acuáticas y biológicas. El alimento, disponibilidad y calidad del agua, y las influencias antropogénicas todas fueron regularmente identificadas como variables importantes que se correlacionan con la distribución de nutrias, pero son a menudo medidas en una variedad de maneras, o identificadas en los modelos mediante variables proxy ó sucedáneas, porque la disponibilidad de los datos relevantes fue escasa o ausente. La escala, el enfoque, y la selección de métricas, todas necesitan ser cuidadosamente consideradas para cada estudio, pero entender los temas de medición y las limitaciones identificadas en los modelos, deberían ayudar a mejorar el modelado de nutrias en el futuro. La revisión de información de este paper de review puede ayudar a futuros esfuerzos respecto de procesos de modelado, tipos de datos usados, métodos de recolección de datos, y variables/métricas a incluir. Esta información debería de todos modos ser cuidadosamente evaluada antes de usarse en áreas de estudio o especies específicas, y como una base para desarrollar métodos innovadores y más efectivos.
Vuelva a la tapa