SciELO - Scientific Electronic Library Online

 
vol.87 número3Diversidad genética en grupos de monos aulladores de manto (Alouatta palliata mexicana) en la Reserva de la Biosfera Los Tuxtlas (Veracruz, México)¿Dónde cultivar el árbol milagro, Moringa oleifera, en México? Un análisis de su distribución potencial índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista mexicana de biodiversidad

versión On-line ISSN 2007-8706versión impresa ISSN 1870-3453

Rev. Mex. Biodiv. vol.87 no.3 México sep. 2016

https://doi.org/10.1016/j.rmb.2016.07.001 

Ecología

Sobre la relación entre idoneidad del hábitat y la abundancia poblacional bajo diferentes escenarios de dispersión

On the relationship between habitat suitability and population abundance under different dispersal scenarios

Luis Alfredo Osorio-Olveraa  * 

Manuel Falconib 

Jorge Soberónc 

a Facultad de Ciencias, Universidad Nacional Autónoma de México, Circuito exterior s/n, Cd. Universitaria, 04510 , Ciudad de México, México

b Departamento de Matemáticas, Facultad de Ciencias, Universidad Nacional Autónoma de México, Circuito exterior s/n, Cd. Universitaria, 04510 , Ciudad de México, México

c Biodiversity Institute, University of Kansas, Dyche Hall, 1345 Jayhawk, Boulevard, Lawrence, KS, 66045, Estados Unidos


Resumen:

Trabajos recientes basados en modelación del nicho para estimar la distribución de las especies han aportado evidencias de que la abundancia de una especie decrece conforme aumenta la distancia al centroide del nicho climático de la especie. Esta relación empírica se ha determinado utilizando modelos correlativos que relacionan datos de presencia/ausencia de una especie con factores abióticos, bajo el supuesto de que las especies se encuentran en equilibrio con su ambiente, pero sin tomar en cuenta explícitamente procesos poblacionales. La falta de un mecanismo que relacione explícitamente procesos poblacionales con el nicho de la especie dificulta la interpretación de los resultados y la obtención de conclusiones generales. En este trabajo se explora la relación entre la idoneidad del hábitat y la abundancia bajo diferentes escenarios de dispersión, utilizando una especie virtual. Se encontró que, en un escenario sin dispersión, hay una clara correlación entre la abundancia y la idoneidad del hábitat, pero en aquellos en los que se incluye el proceso de dispersión, a medida que se incrementa la proporción de individuos que se dispersan, la correlación entre abundancia e idoneidad del hábitat disminuye. Estos resultados apoyan la idea de que existe una relación entre idoneidad del hábitat y las abundancias, sin embargo, también muestran que esta se encuentra modulada por la dispersión.

Palabras clave: Modelos de distribución de especies; Modelos híbridos; Núcleo de dispersión; Especie virtual

Abstract:

Recent studies based on niche modeling to estimate the distribution of species have provided evidence that the abundance of a species decreases as the distance to the centroid of the climatic niche of the species increases. This empirical relationship is determined using correlative models that relate the presence/absence of a species with abiotic factors under the assumption that species are in equilibrium with its environment, but without taking into account population processes explicitly. The lack of a mechanism that explicitly relates population processes to the niche hinders to interpret results and obtain general conclusions. In this paper the relationship between habitat suitability and abundance under different dispersal scenarios is explored using a virtual species. A clear correlation between abundance and habitat suitability in scenarios without dispersal was found, on the other hand, this correlation decreased as the proportion of individuals who are dispersed increases. These results support the idea that there is a relationship between habitat suitability and abundances, however, they also show that this relationship is modulated by dispersal.

Keywords: Species distribution models; Hybrid models; Dispersal Kernel; Virtual species

Introducción

Entender los patrones de distribución de las especies es uno de los tópicos más importantes de la ecología (Sagarin, Gaines y Gaylord, 2006; VanDerWal, Shoo, Johnson y Williams, 2009). Los patrones de distribución están determinados por factores ecológicos, evolutivos y geográficos, que hacen que el estudio de la distribución de las especies sea un problema complejo (Cain, 1994; Gaston, 2003; Udvary, 1969). Desde el trabajo clásico de Grinnell (1917), se ha considerado que las preferencias ambientales de una especie, denominadas «nicho ecológico», determinan, a través de los parámetros demográficos de cada población, las posibilidades de que una especie colonice exitosamente una localidad. En conjunto con el nicho, parámetros demográficos como la reproducción, la mortalidad y la dispersión, determinan la dinámica metapoblacional y por ende la abundancia y la distribución a diferentes escalas espaciales (Guisan y Thuiller, 2005; Lira-Noriega, Soberón y Milller, 2013; Marion et al., 2012; Schurr et al., 2007; Thuiller et al., 2008).

La forma en la que se relaciona el nicho ecológico con la abundancia es un tema que ha llamado la atención (Siqueira, Bini, Cianciaruso, Oliveira y Trivinho-Strixino, 2009; VanDerWal et al., 2009). En términos generales, se espera que la reproducción y la sobrevivencia sean altas en localidades ubicadas al interior del nicho de una especie y por ende, se alcancen abundancias poblacionales relativamente altas, y lo contrario en localidades periféricas (de Moraes y Viveiros, 2012; Maguire, 1973; Martínez-Meyer, Díaz-Porras, Townsend y Yáñez-Arenas, 2013; Siqueira et al., 2009; VanDerWal et al., 2009). Esta idea fue sugerida por Hutchinson (1957), quien expresa la posibilidad de que exista una región «óptima» dentro del nicho, es decir, con condiciones óptimas para la sobrevivencia de la especie, y regiones con condiciones menos favorables en la periferia. Maguire (1973) es el primero en proponer explícitamente la idea de que diferentes regiones del espacio de nicho corresponden a diferentes valores de la tasa intrínseca de crecimiento y de la abundancia poblacional. Una idea derivada de lo anterior es la hipótesis del «centroide del nicho», la cual plantea que las mayores abundancias se alcanzan en el centro del nicho y conforme una población esté más alejada de este, tenderá a presentar menores abundancias (Martínez-Meyer et al., 2013).

Se ha empezado a acumular evidencia empírica de que la posición de una población en el espacio de nicho determina su abundancia (de Moraes y Viveiros, 2012; Escalante y Martínez-Meyer, 2013; Martínez-Meyer et al., 2013; Siqueira et al., 2009; Torres et al., 2012; VanDerWal, 2009; Yáñez-Arenas, Martínez-Meyer, Mandujano y Rojas-Soto, 2012; Yáñez-Arenas, Guevara, Martínez-Meyer, Mandujano y Lobo, 2014). En estos trabajos se utilizan algoritmos de estimación de nichos para obtener una aproximación al nicho fundamental (medido como «idoneidad de hábitat»), y luego se establece una correlación entre la idoneidad de hábitat con datos de abundancia tomados en campo o de la literatura. Otra forma en la que se ha analizado la relación entre el nicho y la abundancia es identificando el centroide del nicho, calculado a partir de los valores promedio de distintas variables climáticas, y estableciendo una correlación entre la idoneidad de hábitat y la distancia a ese centroide. En ninguno de estos dos casos se incluyen explícitamente procesos poblacionales como natalidad, mortalidad o migración, etc. Por tal motivo, si bien estos procedimientos pueden revelar la existencia de una relación entre abundancia y posición en el espacio de nicho, los resultados son correlaciones que no permiten profundizar en el entendimiento de los mecanismos que determinan esta relación (Dormann, 2009; Pagel y Schurr, 2012).

Se han propuesto diversos métodos para integrar aspectos demográficos en el análisis de los patrones de distribución y abundancia. Los modelos espacio-temporales, por ejemplo, combinan medidas del nicho con parámetros demográficos como la dispersión y el crecimiento poblacional. Se les denomina modelos acoplados o híbridos y pueden dar estimaciones de abundancia en un espacio geográficamente explícito (Catterall, Cook, Marion, Butler y Hulme, 2012; Nenzén, Swab, Keith y Araújo, 2012; Pagel y Schurr, 2012). Estos modelos se han utilizado frecuentemente para estimar la distribución de las especies en escenarios de cambio climático (i.e., Franklin, 2010), y para estudiar la relación entre demografía y la dinámica del área de distribución de algunas especies invasoras (Franklin, 2010; Gallien, Münkemüller, Albert, Boulangeat y Thuiller, 2010; Jongejans et al., 2008; Marion et al., 2012; Pagel y Schurr, 2012; Thuiller et al., 2008).

Trabajos como el de Jongejans et al. (2008) utilizan un modelo híbrido que incluye mecanismos poblacionales y dispersión para estimar la dispersión de una planta, pero sin referencia alguna al concepto de nicho. A su vez, Kearny y Porter (2009) explican cómo los modelos mecanísticos basados en procesos biofísicos ecofisiológicos, se pueden utilizar para estimar la distribución de una especie, pero no muestran cómo relacionarlos con los procesos demográficos y de dispersión. Buckley et al. (2010) proponen un modelo híbrido en el cual la tasa de crecimiento es función del ambiente, sin incluir al proceso de dispersión. La relación entre abundancia y nicho se menciona en un par de líneas, pero no se analiza, y finalmente, Pagel y Schurr (2012) utilizan un modelo híbrido sin estructura de edades para analizar cómo la tasa de crecimiento, la capacidad de carga y la densidad poblacional se relacionan con la probabilidad de presencia de una especie.

El objetivo de este trabajo es analizar la relación entre la abundancia y la posición en el espacio de nicho utilizando un modelo híbrido, que incorpora la estructura de edades de la población y que integra algunos parámetros poblacionales básicos, así como el proceso de dispersión. Con esta aproximación novedosa se pretende contribuir a la teoría que estudia la relación entre abundancia poblacional y la posición en el espacio del nicho, de tal forma que permita profundizar en el entendimiento de los mecanismos que determinan esta relación.

Materiales y métodos

Para estudiar la relación entre la idoneidad del hábitat y las abundancias poblacionales bajo diferentes escenarios de dispersión se utilizó una especie virtual. Se asumió que el nicho fundamental de la especie se mide en términos de la idoneidad del hábitat, la cual está definida por una función f(x, y) donde x y y son las variables ambientales que definen las condiciones donde la especie puede sobrevivir. Se considera que la tasa óptima de reproducción de la especie se obtiene en un cierto punto x-,y- del nicho fundamental de modo tal que fx,y es menor que fx-,y- para todo x,yx-,y- .

En este trabajo se utilizó como función fx,y a la función de densidad de probabilidad normal bivariada (ecuación 1) (Broennimann et al., 2012; Saupe et al., 2012), en donde la tasa óptima de reproducción fx-,y- se alcanza en las medias de las variables independientes x y y. La función de distribución multinormal bivariada es la siguiente:

fx,y=12πσxσy1-ρ2exp-121-ρ2x-μx2σx2+y-μy2σy2-2ρx-μxy-μyσxσy, ()1

donde ρ es la correlación entre la variable x y la variable y; μx y μy son las medias de cada una de las variables independientes; σx y σy son las varianzas de ambas. Este método es equivalente a utilizar un método envolvente tipo distancia de Mahalanobis (Farber y Kadmon, 2003).

La función de la ecuación 1 es una función de densidad de probabilidad y por lo tanto esta se reescaló para que sus valores se encuentren en el intervalo cerrado [0,1]. La función de reescalamiento está dada por la expresión:

Sx,y=fx,ymaxfx,yǀx,yR. ()2

Dicha función de reescalamiento (ecuación 2) representa la denominada idoneidad del hábitat (Nenzén et al., 2012; Helaouet, Beaugrand y Edwards, 2013), en la que el valor de S(x, y) = 1 ocurre en los niveles de las variables ambientales donde se alcanza la tasa máxima de crecimiento poblacional.

Modelado de la dinámica espacial de crecimiento de la especie

Para analizar la relación entre abundancia poblacional y posición en el espacio de nicho bajo diferentes escenarios de dispersión, se utilizó un modelo híbrido de dinámica discreta de la forma dada por:

Nt+1lon,lat=Sxlon.lat, ylon,latMNtlon,lat+Δlon,lat, ()3

donde Nt +1 es el vector que contiene las abundancias de la población en cada uno de los estadios en las coordenadas (lon,lat) al tiempo t + 1. La matriz M = (ai,j ) es la matriz de proyección poblacional de tipo Lefkovich (Caswell, 2000) y contiene: 1), las proporciones que pasan de un estadio a otro (elementos ai , j con ij); 2), la proporción de individuos que se mantiene en el mismo estadio (elementos de la diagonal principal) y 3), las fertilidades de los estadios (primer renglón de la matriz). Cada elemento ai , j es una variable aleatoria con una distribución normal N(μi,j). El vector Δ representa el flujo neto de inmigrantes de cada estadio al parche con coordenadas (lon,lat). Con Δk,rNlon,lati se denota a los inmigrantes del estadio i al parche (lon,lat) provenientes del parche (k,r); es decir, Δk,rNlon,lati=Δk,rlon,latNk,ri-Δlon,latk,rNlon,lati , donde ΔγνNγi es el número de individuos de la clase i que migran del parche γ al parche ν; por tanto la inmigración neta Δ está dada por:

Δlon,lat=k,rΔk,rNlon,lat1,k,rΔk,rNlon,lat2,,k,rΔk,rNlon.latn. ()4

S(x, y) es una matriz diagonal de n × n dimensiones con las idoneidades de cada estadio (ecuación 5). Las variables x, y dependen de la posición (lon, lat), por lo tanto,

Sxlon,lat,ylon,lat=Slon,lat1.00Slon,lat2000Slon,latn. ()5

De la estructura del modelo (ecuación 3) se sigue que hay un valor umbral mínimo de idoneidad S 0 por debajo del cual la supervivencia de la especie es nula; es decir, en ausencia de dispersión, si Si (x, y) ≤ S 0 para todo i = 1, ..., n entonces la población en (x, y) tiende a extinguirse debido a que su tasa intrínseca de crecimiento λ será menor a 1. El valor S 0 se calcula como el inverso de la tasa máxima intrínseca de crecimiento

S01λmax, ()6

Nótese que el valor máximo de la tasa intrínseca de crecimiento λmax se alcanza en el valor de idoneidad Sx-,y-=1 .

Parametrización del modelo híbrido

La matriz M que se utilizará en las simulaciones es la correspondiente a la especie Astragalus scaphoides (tabla 1) tomada de la base de datos demográficos COMPADRE (COMPADRE Plant Matrix Database, 2015). Los estadios son latente, pequeña, grande y reproductiva. Para la dispersión se considerarán dos kernels de dispersión que indican el flujo neto de inmigrantes (esto corresponde al término Δ de las ecuaciones 3 y 4). El primer kernel se refiere a una dispersión a celdas cercanas (kernel SDD) y el segundo a la dispersión a celdas lejanas (kernel LDD). El porcentaje de la población de un parche que se dispersa a los parches contiguos SDD, lo denotamos con α y la que se dispersa a los más lejanos LDD, lo denotamos con β. Se considera que la dispersión ocurre solo en el estadio de semilla. El kernel SDD se da entre las 8 celdas contiguas a cada parche, con una proporción de 0.2α colonizadores hacia cada una de las 4 celdas laterales y una proporción de 0.05α hacia cada celda diagonal (fig. 1a). La dispersión a celdas no contiguas está determinada por un kernel (LDD) de dispersión exponencial del tipo siguiente:

Pd=aexp-dcbcon a,b,c ∈0,1,xDmax, ()7

donde Dmax es la distancia máxima a la que puede viajar un individuo y P(d) es la fracción de migrantes que viajan a una distancia d. La proporción máxima de la población que se desplaza está determinada por a. La razón entre c y b define que tan grande es la proporción de individuos que se desplazan a una distancia d. Los valores de los parámetros elegidos fueron a = 0.8, b = 0.8, c = 0.1 y D = 2. La distribución del kernel dados estos parámetros se muestra en la figura 1.

Tabla 1 Matriz de proyección poblacional de Astragalus scaphoides (NR=no reproductiva). 

  Latente Pequeña (NR) Grande (NR) Reproductiva
Latente 0.50 0.03 0.04 5.00
Pequeña (NR) 0.40 0.45 0.04 0.00
Grande (NR) 0.00 0.24 0.37 0.16
Reproductiva 0.10 0.06 0.52 0.80

Figura 1 Kerneles de dispersión. Imagen izquierda, dispersión a las celdas contiguas (SDD). Imagen derecha, dispersión a celdas distantes (LDD) como función de la distancia d al parche inicial S. 

Para definir el nicho fundamental de la especie virtual, se eligieron como variables ambientales la temperatura media anual y la precipitación en el mes más húmedo, por ser las más importantes para el desarrollo y supervivencia de varias especies (Bartlein, Prentice y Webb III, 1986; Woodward, 1987). Para analizar el proceso de dispersión se eligió como «espacio geográfico» una zona al noroeste de México. En este espacio la especie virtual inicia el proceso de dispersión, y se definió tomando como base el nicho de la especie, es decir los sitios del espacio geográfico donde la especie puede subsistir. Los valores de temperatura media anual (bio1) y de precipitación en el mes más húmedo (bio13) fueron obtenidos de la página de Worldclim (Hijmans, Cameron, Parra, Jones y Jarvis, 2005) a una resolución de 30 segundos de arco. Aquí suponemos que la tasa óptima de reproducción de la especie virtual se alcanza en las medias (μx , μy ) de bio1 y bio13 en el espacio geográfico, los cuales fueron 24.2 °C y 260 mm (este es el centroide del nicho). La correlación de las variables nichos fue de ρ = −0.77 con una desviación de σt = 3.89 para la bio1 y σP = 9.85 para bio13. El valor de la tasa máxima intríseca de crecimiento (λmax ) de la matriz M en condiciones de idoneidad máxima S(μxy )=1 es de λmax = 1.6

Utilizando la función rasterToPoints del paquete «raster» de R (Hijmans, 2015), se extrajeron los valores de bio1 y bio13 para el espacio geográfico en cuestión, y se aplicaron las ecuaciones (1) y (2) para obtener la idoneidad S(x, y) en cada una de las celdas (fig. 2). Se eligieron aleatoriamente 10 celdas cuya idoneidad del hábitat fuera mayor al valor de S 0 (en este caso S 0 = 0.625). En estas celdas se colocaron las poblaciones iniciales de la especie virtual y constituyeron el área inicial de distribución (puntos negros en el espacio geográfico, fig. 3).

Figura 2 Nicho fundamental de la especie virtual. La función de distribución f(x, y) que se utilizó para definir al nicho fundamental de la especie virtual, fue la función multinormal bivariada; se eligió como variable x a la temperatura media anual (bio1) y como variable y a la precipitación en el mes más húmedo (bio13). 

Figura 3 Mapa de idoneidad del hábitat. Los puntos negros muestran los sitios donde se colocaron a las 10 poblaciones iniciales de la especie virtual. 

Los escenarios

Se puso a prueba la hipótesis del centroide del nicho (la mayor abundancia se encuentran más cercanos al centroide del nicho) bajo distintos escenarios de dispersión: a) sin dispersión de la especie, i.e. α = 0 = β; b) con dispersión; se simularon un total de 36 escenarios de dispersión, los cuales son resultado de combinar diferentes proporciones de α y β. Los valores posibles de α y β son [0.01, 0.04, 0.06, 0.8, 0.11, 0.14]. Para realizar las simulaciones se utilizó el paquete demoniche de R (Nenzén et al., 2012).

Se realizaron regresiones lineales entre el índice de idoneidad y la abundancia poblacional, usando los coeficientes de determinación. La comparación de las r 2 permitió comparar los resultados entre los diferentes escenarios de dispersión.

Resultados

Modelo sin dispersión

Las mayores abundancias se encontraron en los valores de idoneidad de hábitat S cercanos al óptimo de la especie (i.e. centroide del nicho). La correlación r fue positiva y la r 2 fue de 0.95 (fig. 4). La variación no explicada se debe a la estocasticidad introducida en los parámetros demográficos del modelo.

Figura 4 Relación entre idoneidad del hábitat y la abundancia para el escenario donde no hay dispersión. Imagen izquierda, la correlación entre estas variables fue de r2 =0.95, la banda alrededor de la línea de tendencia representa el intervalo de confianza de la estimación a un nivel de confianza del 95%. Imagen derecha, proyección de las abundancias en el espacio geográfico. 

Modelo con dispersión

En general, los valores que tomó la abundancia fueron más altos conforme la proporción de individuos con dispersión a celdas cercanas SDD fue mayor. Por ejemplo, el número máximo de individuos que alcanzó un parche en el escenario LDD = 0.01 SDD = 0.01 fue de 3,500 individuos, mientras que en el escenario de mayor dispersión (LDD = 0.01 SDD = 0.08 fue de 120,000 individuos (fig. 5).

Figura 5 Proyección al espacio geográfico de las abundancias poblacionales para diferentes escenarios de dispersión a cortas distancias SDD (0.01≤α≤0.08) y valores de dispersión a largas distancias LDD fijos (β=0.01). 

La relación idoneidad del hábitat S y abundancia fue positiva para todos los escenarios. Sin embargo, se observó que conforme incrementó la proporción de α individuos que viajaban a distancias cortas (SDD), los coeficientes de determinación r 2 resultaron considerablemente más bajos (fig. 6). En estos escenarios, conforme incrementó la dispersión, la relación entre abundancia e idoneidad se hace más débil, por lo tanto fue posible encontrar abundancias altas en sitios con idoneidad baja (fig. 6).

Figura 6 Correlaciones entre idoneidad del hábitat y las abundancias poblacionales para diferentes escenarios de dispersión a cortas distancias SDD y a largas distancias LDD. La banda alrededor de las líneas de tendencia representa el intervalo de confianza de la estimación a un nivel de confianza del 95%. 

Las relaciones que se encontraron entre la idoneidad del hábitat y la abundancia poblacional para el escenario donde la proporción de individuos susceptibles a dispersarse a las largas distancias es de 0.01 ≤ β = 0.06, con proporciones de 0.01 ≤ α ≤ 0.14 individuos susceptibles dispersarse a cortas distancias, muestran una correlación r positiva en todos los casos, sin embargo también muestran que tanto r como r 2 disminuyen al incrementar las proporciones α y β (fig. 6).

En términos generales las r 2 obtenidas en las simulaciones para todos los escenarios de dispersión disminuyeron conforme se incrementaron las proporciones de individuos con dispersión a celdas cercanas SDD (α) y de individuos dispersión a celdas no contiguas LDD (β) (fig. 7).

Figura 7 Valores de r2 para todos los escenarios de dispersión a grandes y a cortas distancias. 

Discusión

En este trabajo analizamos la relación entre la abundancia de una especie virtual y su posición en el espacio de nicho, utilizando un modelo «híbrido» con estructura de edades que integra procesos demográficos y de dispersión. A pesar de que en la literatura hay trabajos en los que se utilizan modelos híbridos y mecanicistas para estudiar la distribución y abundancia de las especies, la diferencia principal radica en que el nuestro integra diferentes escenarios de dispersión.

El principal resultado que encontramos es que existe una correlación positiva entre abundancia poblacional y la idoneidad de hábitat, es decir, se confirma la hipótesis del centroide del nicho. Otro resultado fundamental es que esta correlación está regulada por la relación que guardan los parámetros demográficos respecto a los de la dispersión.

Los resultados obtenidos muestran que la intensidad de la correlación positiva entre abundancia e idoneidad del hábitat depende de la capacidad de dispersión de la especie en cuestión: la correlación entre abundancia e idoneidad del hábitat es casi perfecta en un escenario de no dispersión; sin embargo, cuando se incluye el factor dispersión, la correlación declina abruptamente pudiéndose encontrar poblaciones con abundancias mayores en sitios donde la calidad del hábitat es baja. Este resultado concuerda con lo señalado por Pulliam (2000), quien afirmó que las especies frecuentemente se encuentran en sitios donde su tasa de crecimiento es negativa (Pulliam, 2000; Schurr et al., 2012). A este fenómeno se le ha denominado dinámica fuente-sumidero (Pagel y Schurr, 2012; Pulliam, 1988, 2000; Schurr et al., 2012) y depende en gran medida de si las respuestas densodependientes afectan la dinámica de dispersión de la especie.

Nuestros resultados refuerzan las predicciones teóricas informales derivadas de la teoría del nicho de Hutchinson (Maguire, 1973; Schurr et al., 2012) y da elementos teóricos a los resultados empíricos como los de VanDerWal et al. (2009), Escalante y Martínez-Meyer (2013) y Martínez-Meyer et al. (2013), ya que dan una explicación plausible en términos poblacionales a las correlaciones que se obtienen entre abundancia y distancia en el espacio de nicho. Además, sugieren algunas hipótesis acerca de las condiciones en las que los procesos demográficos y de dispersión pueden hacer que la correlación entre abundancia y distancia al centroide del nicho sea baja, por ejemplo, se espera que las especies cuyas tasas de migración son altas, muestren una correlación baja entre su abundancia y su nicho.

Asimismo, nuestro trabajo complementa estudios como el de Holt, Lawton, Gaston y Blanchette (1997) en el que se analiza la relación entre la abundancia y el rango de distribución de la especie. A diferencia del enfoque de Holt et al., que se basa en la idea de que el proceso de dispersión de la especie ha concluido y no consideran el efecto de la migración-inmigración sobre las abundancias locales, nosotros mostramos cómo el proceso de dispersión y el nicho pueden determinar la abundancia de las poblaciones en el rango geográfico.

Finalmente, si bien el enfoque de los modelos híbridos es prometedor, existe una limitante práctica para incluir el término de dispersión de estos modelos, ya que la estimación de los kernels de dispersión es difícil, los datos son escasos y kerneles con colas largas son matemáticamente difíciles de manejar (Barve et al., 2011). A pesar de que en este trabajo se enfatizó el papel de la dispersión en la relación entre abundancia y posición en el espacio de nicho, otra limitante importante para la comprensión de esta relación son las interacciones, de manera que la ecuación de abundancia de una especie en un espacio geográficamente explícito debe ser una función de su nicho fundamental, su capacidad de dispersión y las interacciones que tenga con las demás especies (Soberón, 2010). El trabajo que presentamos aquí pretende ser una base teórica sobre la cual avanzar incluyendo este otro tipo de factores, lo cual permitirá entender mejor la compleja red de factores que determinan la geometría y la estructura interna del área de distribución de una especie.

Agradecimientos

Se agradece al Posgrado en Ciencias Biológicas por la formación académica de Luis Osorio. Asimismo se agradece al PAPIIT-UNAM IN112715 (2015) por el apoyo parcial brindado para la realización de esta investigación. Jorge Soberón agradece al proyecto NSF 1208472. Agradecemos a Andrés Lira-Noriega y a los revisores anónimos cuyos comentarios nos ayudaron a mejorar la presentación del trabajo.

Referencias

Bartlein et al., 1986 Bartlein P.J, Prentice I.C, Webb III T. Climatic response surfaces from pollen data for some eastern North American Taxa. Journal of Biogeography. 1986; 13:35-57 [ Links ]

Barve et al., 2011 Barve N, Barve V, Jiménez-Valverde A, Lira-Noriega A, Maher S.P, Peterson A.T, et al. The crucial role of the accessible area in ecological niche modeling and species’ distribution modeling. Ecological Modelling. 2011; 222:1810-9 [ Links ]

Broennimann et al., 2012 Broennimann O, Fitzpatrick M.C, Pearman P.B, Petitpierre B, Pellissier L, Yoccoz N.G, et al. Measuring ecological niche overlap from occurrence and spatial environmental data. Global Ecology and Biogeography. 2012; 21:481-97 [ Links ]

Buckley et al., 2010 Buckley L, Urban M, Angilletta M, Crozier L, Rissler L, Sears M. Can mechanism inform species’ distribution models?. Ecology Letters. 2010; 13:1041-54 [ Links ]

Cain, 1944 Cain S.A. Foundations of plant geography. New York: Harper & Brothers; 1944. [ Links ]

Caswell, 2000 Caswell H. Matrix population models: construction, analysis and interpretation. Sunderland, Massachusetts: Sinauer Associates; 2000. [ Links ]

Catterall et al., 2012 Catterall S, Cook A.R, Marion G, Butler A, Hulme P.E. Accounting for uncertainty in colonization times: a novel approach to modelling the spatio-temporal dynamics of alien invasions using distribution data. Ecography. 2012; 35:901-11 [ Links ]

COMPADRE, 2015 COMPADRE Plant Matrix Database (2015). Max Planck Institute for demographic research (Germany). (consultado 19 Sep 2015) Disponible en: www.compadre-db.org. [ Links ]

de Moraes y Viveiros, 2012 De Moraes M, Viveiros C. Does environmental suitability explain the relative abundance of the tailed tailless bat, Anoura caudifer?. Natureza a Conservacao. 2012; 10:221-7 [ Links ]

Dormann, 2009 Dormann C.F. Response to comment on methods to account for spatial autocorrelation in the analysis of species’ distributional data: a review. Ecography. 2009; 32:379-81 [ Links ]

Escalante y Martínez-Meyer, 2013 Escalante T, Martínez-Meyer E. Ecological niche modeling and wildlife management units (UMAS): an application to deer in campeche, México. Tropical and Subtropical Agroecosystems. 2013; 16:183-91 [ Links ]

Farber y Kadmon, 2003 Farber O, Kadmon R. Assessment of alternative approaches for bioclimatic modeling with special emphasis on the Mahalanobis distance. Ecological Modeling. 2003; 160:115-30 [ Links ]

Franklin, 2010 Franklin J. Moving beyond static species’ distribution models in support of conservation biogeography. Diversity and Distributions. 2010; 16:321-30 [ Links ]

Gallien et al., 2010 Gallien L, Münkemüller T, Albert C.H, Boulangeat I, Thuiller W. Predicting potential distributions of invasive species: where to go from here?. Diversity and Distributions. 2010; 16:331-42 [ Links ]

Gaston, 2003 Gaston K.J. The structure and dynamics of geographic ranges. Oxford series in ecology and evolution. New York: Oxford University Press; 2003. [ Links ]

Grinnell, 1917 Grinnell J. The niche-relationships of the California thrasher. The Auk. 1917; 34:427-33 [ Links ]

Guisan y Thuiller, 2005 Guisan A, Thuiller W. Predicting species’ distribution: offering more than simple habitat models. Ecology Letters. 2005; 8:993-1009 [ Links ]

Helaouet et al., 2013 Helaouet P, Beaugrand G, Edwards M. Understanding long-term changes in species’ abundance using a niche-based approach. PloS one. 2013; 8:e79186 [ Links ]

Hijmans, 2015 Hijmans, R. J. (2015). Raster: geographic data analysis and modeling. R package version 2.4-18. (consultado 4 Nov 2015). Disponible en: http://CRAN.R-project.org/package=raster. [ Links ]

Hijmans et al., 2005 Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. G. y Jarvis, A. (2005). The WorldClim interpolated global terrestrial climate surfaces. (consultado el 01 Oct 2014). Disponible en: http://www.worldclim.org/. [ Links ]

Holt et al., 1997 Holt R.D, Lawton J.H, Gaston K.J, Blanchette C.A. On the relationship between range size and local abundance: back to basics. Oikos. 1997; 78:183-90 [ Links ]

Hutchinson, 1957 Hutchinson G.E. Concluding remarks. Cold Spring Harbor Symposia on Quantitative Biology. 1957; 22:415-27 [ Links ]

Jongejans et al., 2008 Jongejans E, Shea K, Skarpaas O, Kelly D, Sheppard A.W, Woodburn T.L. Dispersal and demography contributions to population spread of Carduus nutans in its native and invaded ranges. Journal of Ecology. 2008; 96:687-97 [ Links ]

Kearny y Porter, 2009 Kearny M, Porter W. Mechanistic niche modelling: combining physiological and spatial data to predict species’ ranges. Ecology Letters. 2009; 12:334-50 [ Links ]

Lira-Noriega et al., 2013 Lira-Noriega A, Soberón J, Milller C. Process-based and correlative modeling of desert mistletoe distribution: a multiscalar approach. Ecosphere. 2013; 4:1-23 [ Links ]

Maguire, 1973 Maguire, B. Jr. Niche response structure and the analytical potentials of its relationship to the habitat. The American Naturalist. 1973; 107:213-46 [ Links ]

Marion et al., 2012 Marion G, McInerny G.J, Pagel J, Catterall S, Cook A.R, Hartig F, et al. Parameter and uncertainty estimation for process-oriented population and distribution models: data, statistics and the niche. Journal of Biogeography. 2012; 39:2225-39 [ Links ]

Martínez-Meyer et al., 2013 Martínez-Meyer E, Díaz-Porras D, Townsend P, Yáñez-Arenas C. Ecological niche structure and range wide abundance patterns of species. Biology Letters. 2013; 9:1-5 [ Links ]

Nenzén et al., 2012 Nenzén H.K, Swab R.M, Keith D.A, Araújo M.B. Demoniche - an R-package for simulating spatially-explicit population dynamics. Ecography. 2012; 35:577-80 [ Links ]

Pagel y Schurr, 2012 Pagel J, Schurr F.M. Forecasting species’ ranges by statistical estimation of ecological niches and spatial population dynamics. Global Ecology and Biogeography. 2012; 21:293-304 [ Links ]

Pulliam, 1988 Pulliam H.R. Sources, sinks, and population regulation. The American Naturalist. 1988; 132:652-61 [ Links ]

Pulliam, 2000 Pulliam H.R. On the relationship between niche and distribution. Ecology Letters. 2000; 3:349-61 [ Links ]

Sagarin et al., 2006 Sagarin D.R, Gaines S, Gaylord B. Moving beyond assumptions to understand abundance distributions across the ranges of species. Trends in Ecology and Evolution. 2006; 21:524-30 [ Links ]

Saupe et al., 2012 Saupe E.E, Barve V, Myers C.E, Soberón J, Barve N, Hensz C.M, et al. Variation in niche and distribution model performance: the need for a priori assessment of key causal factors. Ecological Modelling. 2012; 237-238:11-22 [ Links ]

Schurr et al., 2007 Schurr F.M, Midgley G.F, Rebelo A.G, Reeves G, Poschlod P, Higgins S.I. Colonization and persistence ability explain the extent to which plant species fill their potential range. Global Ecology and Biogeography. 2007; 16:449-59 [ Links ]

Schurr et al., 2012 Schurr F.M, Pagel J, Cabral J.S, Groeneveld J, Bykova O, O’Hara R.B, et al. How to understand species’ niches and range dynamics: a demographic research agenda for biogeography. Journal of Biogeography. 2012; 39:2146-62 [ Links ]

Siqueira et al., 2009 Siqueira T, Bini L, Cianciaruso M, Oliveira F, Trivinho-Strixino S. The role of niche measures in explaining the abundance-distribution relationship in tropical lotic chironomids. Hydrobiologia. 2009; 1-10 [ Links ]

Soberoń, 2010 Soberoń J. Niche and area of distribution modeling: a population ecology perspective. Ecography. 2010; 33:159-67 [ Links ]

Thuiller et al., 2008 Thuiller W, Albert C, Araujo M.B, Berry P.M, Cabeza M, Guisan A, et al. Predicting global change impacts on plant species’ distributions: future challenges. Perspectives in Plant Ecology Evolution and Systematics. 2008; 9:137-52 [ Links ]

Torres et al., 2012 Torres M.N, Júnior P.D.M, Santos T, Silveira L, Jácomo-de Almeida A.T, Diniz-Filho A.F. Can species’ distribution modelling provide estimates of population densities? A case study with jaguars in the Neotropics. Diversity and Distributions. 2012; 18:615-27 [ Links ]

Udvary, 1969 Udvary M.D.F. Dynamic zoogeography: with special reference to land animals. New York: Van Nostrand Reinhold Comp; 1969. [ Links ]

VanDerWal et al., 2009 VanDerWal J, Shoo L, Johnson C, Williams S. Abundance and the environmental niche: environmental suitability estimated from niche models predicts the upper limit of local abundance. The American Naturalist. 2009; 174:282-91 [ Links ]

Woodward, 1987 Woodward F.I. Climate and plant distribution. Cambridge: Cambridge University Press; 1987. [ Links ]

Yáñez-Arenas et al., 2014 Yáñez-Arenas C, Guevara R, Martínez-Meyer E, Mandujano S, Lobo J.M. Predicting species’ abundances from occurrence data: effects of sample size and bias. Ecological Modelling. 2014; 294:36-41 [ Links ]

Yáñez-Arenas et al., 2012 Yáñez-Arenas C, Martínez-Meyer E, Mandujano S, Rojas-Soto O. Modelling geographic patterns of population density of the white-tailed deer in central Mexico by implementing ecological niche theory. Oikos. 2012; 121:2081-9 [ Links ]

** La revisión por pares es responsabilidad de la Universidad Nacional Autónoma de México.

Recibido: 26 de Febrero de 2015; Aprobado: 07 de Abril de 2016

* Autor para correspondencia. Correo electrónico: luismurao@comunidad.unam.mx (L.A. Osorio-Olvera).

Creative Commons License Este es un artículo publicado en acceso abierto bajo una licencia Creative Commons