Establishing species' environmental requirements to understand how the southernmost species of South American pitvipers (Bothrops,Viperidae) are distributed: A niche‐based modelling approach

Three Bothrops species are known to be present along an extensive and generally xeric band of the south‐eastern portion of South America. Yet, the environmental factors responsible for the maintenance of the structure of this community have remained undetermined. To have a better understanding of snakes' geographic range limits in Argentina, we modelled the ecological niche of Bothrops alternatus, B. ammodytoides and B. diporus to identify areas of their occurrence and sympatry. We used snakes presence records and pseudo‐absence data with topographic (altitude and slope), vegetation and 19 bioclimatic variables as predictors. These variables were evaluated for each species, using Multiple Logistic Regression Analysis. The distribution of B. alternatus was correlated with temperature, precipitation and vegetation; that of B. ammodytoides with altitude, slope, precipitation, temperature and vegetation; and that of B. diporus related only to precipitation and temperature variables. Five variables appeared in all three models (i.e. annual precipitation, precipitation seasonality, precipitation of the driest month, mean diurnal temperature range and temperature range), although with different effects in different species. Probable sympatry areas were identified, for B. alternatus – B. ammodytoides in a narrow area, for B. diporus – B. alternatus in a wide area and B. ammodytoides – B. diporus in a small area. While some factors constrain the distribution of the southernmost pitvipers community in South America to spatially non‐overlapping ranges, others facilitate sympatry by allowing spatial coexistence. We found high sympatry between phylogenetically distantly related pitvipers. Niche‐based Geographic Information Systems modelling allowed us to identify the environmental factors correlated with each pitviper range in Argentina and in its contact zones, and to understand species distributional limits.