Introducción

El conocimiento del uso y de la selección del hábitat posee importantes implicaciones para explicar la distribución de los ungulados silvestres y nos permite diferenciar entre la disponibilidad y la calidad de diferentes hábitat (^{Manly et al. 2002}). Dicha Información resulta esencial para establecer acciones de manejo, como reintroducciones, translocaciones, conservación, aprovechamiento y establecimiento de áreas protegidas (^{Araujo y Williams 2000}).

El borrego cimarrón es un buen ejemplo para investigar los mecanismos de selección del hábitat, debido a que se trata de una especie oportunista en la búsqueda de su forraje y requiere una diversidad topográfica y paisajes escarpados para detectar y evadir a sus depredadores (^{Valdez y Krausman 1999}). La relación de estos componentes del hábitat con el borrego cimarrón también se ha utilizado para explicar la segregación sexual. Algunos estudios encuentran una correlación directa entre la disponibilidad del terreno de escape y el tamaño poblacional de esta especie (^{Holl 1982}; ^{Bleich et al. 1997}; ^{Divine et al. 2000}; ^{Mckinney et al. 2003}). De tal forma, la evaluación de la selección que hace el borrego cimarrón del hábitat es fundamental para determinar la distribución y el número de organismos que puede soportar un área (^{Brewer et al. 2014}).

A pesar de la importancia de evaluar las preferencias de hábitat del borrego cimarrón, son pocos los estudios realizados en México, y la mayoría se ha enfocado en monitoreos poblacionales (^{Sandoval et al. 2014}). La carencia de estudios sobre el hábitat del borrego cimarrón en Baja California no permite pronosticar los posibles efectos del cambio climático (e. g. la reducción en la disponibilidad de alimento como consecuencia de periodos largos sin lluvias, reducción que ha sido una de las causas de extinción local y extirpación de poblaciones de borrego cimarrón en zonas áridas; ^{Epps et al. 2004}). Los objetivos de esta investigación fueron: 1) Determinar la disponibilidad del terreno de escape, para lo que se utilizó un modelo de rugosidad del terreno en cuatro categorías (plano, sierras bajas, cañones y sierras altas). 2) Evaluar la disponibilidad de la cobertura vegetal mediante un Índice de Vegetación de Diferencia Normalizada (NDVI), que en estudios recientes ha sido útil para explicar la distribución de ungulados silvestres en zonas áridas, e. g. el venado bura (^{Marshal et al. 2006}), el impala (^{VanBommel et al. 2006}) y el borrego cimarrón (^{Sesnie et al. 2012}; ^{Hoglander et al. 2015}), y 3) Analizar las preferencias de uso y la selección que hace el borrego cimarrón de las variables del hábitat. Esta es la primera evaluación de preferencia de hábitat del borrego cimarrón en Baja California. Se espera que los resultados de este estudio sirvan para la identificación de sitios a fin de manejar y conservar el Ovis canadensis y su hábitat en Sierra Santa Isabel, Baja California.

Material y métodos

Área de estudio. La investigación se realizó en la Sierra Santa Isabel, localizada en la región central de Baja California (30.460° N, -115.191° W y 29.950° N, -114.757° W; Figura 1). La Sierra se caracteriza por contar con el mayor hábitat continuo para el borrego cimarrón en el estado (2,072.45 km²). Asimismo, es la que presenta el mejor grado de conservación. Se ha estimado un tamaño poblacional entre 300 - 400 ejemplares (^{DeForge et al. 1993}; ^{Martínez-Gallardo y Eaton-González 2008}; ^{Martínez 2011}; ^{Lee et al. 2012}).

Figura 1 Localidades de estudio en Sierra Santa Isabel, Baja California. Los cuadros en negro representan las áreas de monitoreo del borrego cimarrón. Los Hemes (1) y El Zamora (2). 

En la Sierra, se localizan dos ecorregiones: por un lado, la parte baja del Desierto de Colorado en los valles y lomeríos cercanos al Golfo de California, donde las principales especies vegetales son Larrea tridentata, Fouquieria splendens, Agave deserti, Psorathamnus spinosus, Olneya tesota y Prosopis sp.; por otro lado, el Desierto Central, que corresponde a la mayor parte de la Sierra; allí, las principales especies son Ambrosia dumosa, Prosopis microphylla, Bursera microphylla, Fouquieria columnaris, Cercidium microphyllum y Pachycereus discolor (^{González-Abraham et al. 2010}). El clima es semicálido, con 28 °C de temperatura media anual. La precipitación media anual en las serranías cercanas al Golfo de California es menor a 50 mm; en cambio, en la parte central y norte de la sierra, la precipitación oscila entre los 100 y 150 mm (^{SMN 2013}). Los principales depredadores del borrego cimarrón en la Sierra Santa Isabel son el puma (Puma concolor), el gato montés (Lynx rufus), el coyote (Canis latrans) y el águila real (Aquila chysaetos).

Localización de borregos cimarrones. Se seleccionaron dos localidades para observar a los borregos cimarrones, con un área de 10 x 10 km cada una. La primera, Los Hemes, se localiza en la región central de Sierra Santa Isabel. Dicha localidad se caracteriza por ser accesible y cuenta con una tinaja principal con agua la mayor parte del año. La altitud oscila entre 600 y 900 metros, y los rasgos topográficos son cañones y planicies en las partes altas. La segunda, El Zamora, se ubica en el sur de la Sierra. Es la localidad con mayor disponibilidad de agua para la fauna silvestre. Se caracteriza por tener cañones y cumbres a las que solo se puede acceder a pie (Figura 1). Estas localidades fueron seleccionadas debido a que son frecuentadas por los borregos cimarrones, dada su posibilidad de acceder a ellas en diferentes épocas del año y porque los dueños de los terrenos lo autorizan.

Mediante recorridos terrestres de 10 km de longitud efectuados en diferentes estaciones del año (julio del 2011 a febrero del 2012). Se realizaron observaciones directas de borregos cimarrones individuales o en grupos, en tres periodos: sequía (julio a septiembre), lluvias (octubre a diciembre) y poslluvias (diciembre a febrero). Los tres periodos corresponden aproximadamente a las temporadas de segregación sexual, apareamiento y crianza (^{Rubin et al. 2000}). Durante ocho días, en cada periodo y recorriendo los mismos transectos (cuatro en total), dos equipos de tres personas realizaron por separado búsquedas y observaciones directas de borregos cimarrones. Utilizaron binoculares 10 x 40. La colecta de datos se realizó a lo largo de cinco horas en la mañana (7:00 a 12:00) y cuatro en la tarde (14:00 a 18:00).

En cada observación individual o de grupos de borregos cimarrones, los observadores acudieron al sitio para ubicarlo espacialmente con un GPS. La edad y género de los borregos se estableció con base en la descripción propuesta por ^{Monson y Summer (1980}), quienes definen cuatro clases de edad para los machos, considerando el tamaño del cuerpo y el tipo de cuernos: Clase I (cuernos con un cuarto de curvatura) con edades entre 2 y 3 años. Clase II (cuernos con un medio curvatura), con edades entre tres y seis años. Clase III (el cuerno con tres cuartos de curvatura), y Clase IV (el cuerno con curvatura completa y la punta generalmente desgastada; puede superar la órbita ocular) para borregos cimarrones con edades superiores a los seis años. Las hembras se clasificaron en jóvenes (cuernos con inicio de curvatura), con edades entre uno y dos años, y adultas (cuernos largos con curvatura en forma de L), con edades mayores a los dos años. Por último, se nombró añero de ambos sexos a borregos cimarrones con edades menores a un año. En ellos era fue visible el inicio del crecimiento del cuerno, pero no resultó claro definir su sexo.

Terreno de escape. Por medio del portal del Servicio Geológico de Estado Unidos (USGS), se obtuvo un modelo digital de elevación (MDE). El modelo fue adquirido durante la misión SRTM (Shuttle Radar Topography Mission) y se ubica en el path 38 y row 39, con una resolución espacial de 30 metros para cada píxel. El MDE se procesó en el programa ArcGis 10 (^{ESRI 2010}) y se implementó un algoritmo para calcular el índice de rugosidad del terreno (IRT) en lenguaje Phyton, siguiendo la metodología propuesta por ^{Riley et al. (l999)}.

Los valores de IRT son equivalentes a la cantidad de terreno de escape disponible en el área de distribución del borrego cimarrón. Se obtuvieron a partir de un análisis de vecindario en una matriz de renglones y columnas. Cada celda (píxel) dispone de un valor de altitud procedente del MDE. A partir de la celda central, se calcula la dispersión de elevaciones con sus celdas vecinas. Las unidades de medida del IRT se expresan en metros de diferencia de cota entre puntos distantes entre sí (^{Riley et al. 1999}).

Con base en los valores obtenidos del IRT y en clasificaciones realizadas previamente por ^{Álvarez-Cárdenas et al. (2009}) y ^{Gilad et al. (2013}), se definieron cuatro categorías de rugosidad del terreno con los siguientes intervalos: plano con cerros bajos 0 - 15; poco rugoso con sierras bajas 16 - 30; sierras medianas con cañones 31 - 45, y sierras altas con cañones para valores mayores a 45. Los números corresponden a valores del IRT.

Cobertura vegetal. Mediante el portal USGS, se obtuvo tres imágenes de satélite Landsat-TM (path 38 y row 39) con una resolución espacial de 30 metros. Las fechas de adquisición fueron las siguientes: 7 de septiembre de 2011 (periodo de sequía), 16 de noviembre de 2011 (periodo de lluvias) y 12 de febrero de 2012 (periodo de poslluvias). Las imágenes fueron procesadas y co-registradas, utilizando el programa ENVI 4.7 (^{EXELIS 2012}).

Para cada localidad se calculó un Índice de Vegetación de Diferencia Normalizada (NDVI), considerado como uno de los mejores estimadores de cobertura vegetal mediante sensores remotos (^{Purevdorj et al.1998}; Pettoreli 2013). A los valores de NDVI se les aplicó la corrección radiométrica propuesta por ^{Guyot y Gu (1994}), recomendable para imágenes del sensor TM en ecosistemas áridos. Esta corrección se describe mediante la siguiente ecuación: NDVI = IRC - 0.801 × R ÷ IRC + 0.801 × R.

IRC y R indican la reflectividad del píxel en la banda del infrarojo cercano y del rojo para el sensor TM. Los valores positivos cercanos a 1 indican una mayor cobertura vegetal, mientras que los valores de 0 y negativos se relacionan con una baja cobertura y con suelo desnudo.

Se establecieron tres categorías de cobertura vegetal para cada imagen, a partir del valor de NDVI de cada píxel calculado con el módulo statistics en ENVI 4.7 (^{EXELIS 2012}). El tercio inferior corresponde a la categoría baja; el tercio intermedio, a la categoría media, y el tercio superior, a la categoría alta (^{VanBommel et al. 2006}).

Las categorías de cobertura vegetal para cada imagen se importaron al programa ArcGis 10 (^{ESRI 2010}). Con la herramienta de estadística focal (Focal Statistics), se obtuvo la categoría de cobertura vegetal en aproximadamente una hectárea, con una matriz de 90 m x 90 m. Esta dimensión espacial representa el área mínima que utilizan los borregos cimarrones para alimentarse (^{Valdez y Krausman 1999}; ^{Sandoval et al. 2014}).

Análisis estadísticos. Se realizó un análisis de uso-disponibilidad del hábitat, cuyo objetivo fue relacionar los registros de uso de las variables del hábitat por el borrego cimarrón con su disponibilidad en un área específica. Con las herramientas de análisis espacial (Spatial Analysis) en ArcGis 10, se determinó la proporción disponible de las categorías rugosidad y cobertura vegetal. Posteriormente, la información de registros de observación de borregos cimarrones de cada periodo se exportó como un vector de puntos al programa ArcGis 10. Con los módulos de tabulación de área y un análisis aleatorio en el mismo programa, se obtuvo la frecuencia de registros observados y esperados de acuerdo con la disponibilidad de cada categoría del hábitat.

Con la información de uso (frecuencia de observación) y disponibilidad (categorías de rugosidad y NDVI), se exploró la hipótesis de selección de hábitat, la cual establece que cuando el recurso tiene una frecuencia de utilización proporcional a su disponibilidad, se utiliza al azar y se define como uso. En cambio, la selección se refiere a la preferencia de un recurso; es decir, su frecuencia de utilización está por encima de su disponibilidad.

Para probar la hipótesis, se realizó una prueba Chi Cuadrada entre los valores de las frecuencias de uso observadas y esperadas, de acuerdo con la proporción disponible de las categorías de hábitat. Si no se encuentran diferencias significativas entre estos dos valores, se acepta la hipótesis estadística de utilización al azar. En contraste, si se obtienen diferencias significativas, se rechaza la hipótesis de utilización al azar y se continúa con la construcción de intervalos de confianza de Bonferroni al 95 % para determinar el patrón de utilización del recurso en términos de uso, rechazo o selección (^{Neu et al. 1974}; ^{Marcum y Loftsgaarden 1980}).

Resultados

Se obtuvo un total de 271 observaciones de borregos cimarrones. En el periodo de sequía, se registró 40.2 %; en el periodo de lluvias, 36.5 %, y en poslluvias, el menor registro: 23.2 % (Tabla 1). En la localidad El Zamora, la observación de hembras y añeros representó 79.6 % del total de 162 observaciones. La observación de machos fue de 15.4 %, y de los corderos, de 5.0 %. La mayor proporción añero/hembra fue observada durante poslluvias y fue de 70.5 %, lo que indicó que por cada 10 hembras, siete tenían añero. La mayor proporción macho/hembra se registró durante el periodo de lluvias, con 35.2 %, es decir, por cada macho se observaron tres hembras. La observación de corderos fue baja, con ocho registros durante el periodo de poslluvias (Tabla 1).

Tabla 1 Resultados de los conteos de borrego cimarrón en Sierra Santa Isabel, Baja California. 

En la localidad de Los Hemes, se observaron menos borregos cimarrones que en El Zamora, con un total de 109 registros. En el periodo de sequía, se obtuvo 43.1 % (n = 47) de los registros, seguido de las lluvias, con 38.5 % (n = 42), y el periodo de poslluvias, con 18.3 % (n = 20). Se registraron más hembras (49.5 %), seguidas de machos (31.1 %) y añeros (17.4 %). Solo se observaron dos corderos. Esta localidad se caracterizó por un mayor registro de machos clase III y IV, en comparación con la localidad El Zamora (Tabla 1).

Selección del terreno de escape. La selección que hicieron los borregos cimarrones de las categorías de rugosidad resultó significativamente diferente (χ² = 9.83, P < 0.05). Los terrenos planos fueron los más disponibles para las dos localidades: en Los Hemes, con 45.5 % (4,752 ha), y en El Zamora, con 30.4 % (3,045 ha). Sin embargo, esta categoría fue evitada por los borregos cimarrones en las dos localidades de estudio. En El Zamora, los rebaños de hembras (n = 86) y añeros (n = 43) seleccionaron las categorías de media y alta rugosidad, que tuvieron disponibles en 40.8 % (4,088 ha; Tabla 2). En Los Hemes, los borregos cimarrones machos jóvenes y adultos (n = 59) seleccionaron las categorías poco y medianamente rugosas, que tuvieron disponibles en 43.4 % (4,348 ha, Tabla 2). La categoría de alta rugosidad fue la de menor disponibilidad, con 9.0 % (900 ha). No se observó a ningún borrego cimarrón en esta categoría (Figura 2).

Tabla 2 Selección espacial que realizó el borrego cimarrón de las categorías de rugosidad en Sierra Santa Isabel, Baja California. * Valores cuyas diferencias son significativas (P < 0.05). Interpretación: "+" los borregos cimarrones prefieren el hábitat; "-" evitan el hábitat; "=" uso de acuerdo con su disponibilidad. Uso observado (UO). Proporción de uso esperado (E). Proporción de uso observada (Pi). Intervalo de Bonferroni (IB). Interpretación estadística (IE). 

Figure 2 Índice de rugosidad del terreno (IRT) en Sierra Santa Isabel, Baja California. Los círculos en color rojo son las observaciones de los borregos cimarrones en las localidades Los Hemes (1) y El Zamora (2). 

Selección de la cobertura vegetal. Los valores de NDVI obtenidos en la localidad de El Zamora fueron más altos durante los periodos de lluvias (entre 0.20 y 0.40) y poslluvias (entre 0.16 y 0.35). Los valores más bajos se encontraron en el periodo de sequía y fueron menores a 0.10. En Los Hemes, los valores de NDVI más altos se presentaron en poslluvias (0.12 - 0.24) y lluvias (0.11 - 0.24). Los valores más bajos se encontraron en el periodo de sequía (0.01 - 0.12).

En la localidad El Zamora, durante el periodo de sequía, la selección de las categorías de NDVI por los borregos cimarrones fue significativamente diferente (χ² = 11.58, P < 0.05) (Tabla 3). La categoría de NDVI bajo estuvo disponible en 12.7 % (1,277 ha) y fue seleccionada por los borregos cimarrones (Tabla 3). La categoría de NDVI medio fue seleccionada en proporción a su disponibilidad, y la categoría NDVI alto fue evitada. Estas categorías, en conjunto, representaron 87.2 % (8,723 ha; Figura 3).

Tabla 3 Uso espacial que hizo el borrego cimarrón de las categorías de NDVI en la localidad El Zamora, Baja California. * Valores cuyas diferencias son significativas (P < 0.05). Interpretación: "+" los borregos cimarrones prefieren el hábitat; "-" evitan el hábitat; "=" uso de acuerdo con su disponibilidad. Uso observado (UO). Proporción de uso esperado (E). Proporción de uso observada (Pi). Intervalo de Bonferroni (IB). Interpretación estadística (IE). 

Figura 3 Imagen Landsat TM, bandas 4, 3, 2 en monocromático. Fecha: 7 de septiembre de 2011, Sierra Santa Isabel, Baja California. Los puntos amarillos son los registros de borregos cimarrones en el periodo de sequía en las localidades 1) Los Hemes y 2) El Zamora. En color rojo, se representan las categorías de NDVI bajo que seleccionaron los borregos cimarrones. 

En el periodo de lluvias, también se encontraron diferencias significativas en la selección que realizaron los borregos cimarrones de las categorías de cobertura vegetal (χ² = 8.72, P < 0.05). La categoría NDVI alto fue seleccionada por los borregos y tuvo una disponibilidad de 40.8 % (4,083.66 ha). La categoría NDVI medio tuvo la mayor disponibilidad, con 47.1 % (4,717 ha) pero fue evitada. La categoría NDVI bajo representó 11.8 % (1,189 ha) y fue seleccionada en proporción a su disponibilidad (Tabla 3). En el periodo de poslluvias no se encontraron diferencias estadísticas significativas (χ² = 2.67, P > 0.05). Las categorías NDVI bajo, con 25.3 % (2,538 ha), y NDVI alto, con 28.3 % (2,836 ha), fueron similares en su disponibilidad (Figura 4).

Figura 4 Imagen Landsat TM, bandas 4, 3, 2 monocromático. Fecha: 16 de noviembre de 2011, Sierra Santa Isabel, Baja California. Los puntos amarillos son los registros de borregos cimarrones en el periodo de lluvias en las localidades 1) Los Hemes y 2) El Zamora. En color verde se representa las categorías de NDVI alto que seleccionaron los borregos cimarrones en la localidad El Zamora. 

En Los Hemes, durante el periodo de sequía, también se encontraron diferencias significativas en la selección que hizo el borrego cimarrón de las categorías de cobertura vegetal (χ² = 13.88, P < 0.05). La categoría de NDVI bajo estuvo disponible en 15.4 % (1,540 ha) y fue seleccionada por los borregos cimarrones. Las categorías evitadas por estos borregos fueron: NDVI medio, con 45.6 % (4,560 ha), y NDVI alto, con 39.0 % (3,900 ha; Tabla 4).

Tabla 4 Selección espacial que realizó el borrego cimarrón de las categorías de NDVI en la localidad Los Hemes, Baja California. * Valores cuyas diferencias son significativas (P < 0.05). Interpretación: "+" los borregos cimarrones prefieren el hábitat; "-" evitan el hábitat; "=" uso de acuerdo con su disponibilidad. Uso observado (UO). Proporción de uso esperado (E). Proporción de uso observada (Pi). Intervalo de Bonferroni (IB). Interpretación estadística (IE). 

En el periodo de lluvias también se encontraron diferencias significativas en la selección que hicieron los borregos de las categorías de cobertura vegetal (χ² = 17.10, P < 0.05). La categoría de NDVI bajo fue seleccionada por los borregos cimarrones y tuvo una disponibilidad de 18.0 % (1,800 ha). Por el contrario, la categoría NDVI medio se usó en proporción a su disponibilidad, y estuvo disponible en 37.0 % (3,700 ha). La categoría de NDVI alto fue evitada y estuvo disponible en 45.0 % (4,500 ha). En el periodo de poslluvias no se encontraron diferencias estadísticamente significativas en la selección que hicieron los borregos cimarrones de las categorías de cobertura vegetal (χ² = 1.96, P > 0.05; Tabla 4).

Discusión

En Sierra Santa Isabel, los borregos cimarrones seleccionaron las categorías del índice de rugosidad del terreno (IRT), que correspondieron a sierras medianas y altas con cañones. En cambio, las categorías de rugosidad plana fueron evitadas en las dos localidades de estudio. Estas diferencias en la selección de la rugosidad del terreno sugieren que la distribución de los borregos cimarrones se relaciona con la disponibilidad del terreno de escape. Los sitios con mayor rugosidad aparentemente son más seguros para las hembras y añeros, mientras que los sitios planos, con menor rugosidad, son evitados.

Estos resultados son consistentes con los registrados en la Sierra del Viejo, en Sonora (^{Tarango 2000}), y en la Sierra del Mechudo, en Baja California Sur (^{Álvarez-Cárdenas et al. 2001}). Se observó un mayor número de hembras y añeros en los sitios con mayor rugosidad, en comparación con sitios planos o con baja rugosidad. Dichas sierras son similares en topografía y elevación a la Sierra Santa Isabel. Previamente se registró que la disponibilidad de terreno de escape es el principal factor que determina la distribución del borrego cimarrón (^{Bleich et al. 1997}; ^{Bangs et al. 2005}; ^{Sappington et al. 2007}; ^{Schroeder et al. 2010}; ^{Sandoval et al. 2014}; ^{Hoglander et al. 2015}).

Aunque en estudio no se realizaron estimaciones poblacionales, la mayoría de las observaciones de borregos cimarrones se presentó en las categorías de mayor rugosidad. ^{Wilson (1968}) menciona que a medida en que se incrementa la cobertura rocosa, resulta mayor la probabilidad de avistamiento de borregos. ^{Divine et al. (2000}) y McKineey et al. (2003) concluyen que existe una correlación directa entre el tamaño de la población y el terreno de escape disponible. Sin embargo, en la localidad de Los Hemes, los borregos cimarrones evitaron los terrenos con mayor rugosidad. En consecuencia, no puede generalizarse que en todos los sitios se encontrará una correlación entre terreno de escape y el tamaño de la población.

La selección de los sitios usados por el borrego cimarrón varió estacionalmente en relación con la cobertura vegetal. En el periodo de sequía, estos animales seleccionaron los sitios que aparentemente poseen menor disponibilidad de alimento (menor cobertura vegetal de acuerdo con los valores del NDVI), pero con alta rugosidad en su terreno. Por lo tanto, son más seguros para los borregos cimarrones (^{Hansen 1980}; ^{Krausman et al. 1999}). Lo anterior sugiere que, para el borrego cimarrón, en el periodo de sequía es más importante seleccionar sitios que ofrezcan mayor posibilidad para escapar y con buena visibilidad para detectar a los depredadores en la distancia, aunque se tenga menor disponibilidad de alimento. Estos resultados son similares a los encontrados en otros estudios de selección del hábitat, donde se reporta que durante el periodo de sequía en los terrenos rugosos y escarpados hay una reducción del forraje disponible como consecuencia de altas temperaturas y ausencia de precipitación (^{Festa-Bianchet 1988}; ^{Andrew et al. 1999}; ^{Hoglander et al. 2015}).

Los resultados obtenidos coinciden con lo reportado para otras especies de ungulados silvestres, como el Odocoileus hemionus (^{Marshal et al. 2006}), Damaliscus lunatus (^{Bro- Jørgensen et al. 2008}, Procapra gutturosa (^{Mueller et al. 2008}) y Aepyceros melampus (^{VanBommel et al. 2006}). Estas especies también dependen de la visibilidad para la detección de sus depredadores y seleccionan hábitats con valores similares de NDVI a los obtenidos en este estudio. No obstante, algunos estudios sobre hábitos alimentarios y preferencias del borrego cimarrón por asociaciones vegetales en zonas áridas, reportan que durante el periodo de sequía se presentan lluvias ocasionales que promueven el crecimiento de especies herbáceas, como Lupinus excubitus, Lotus scoparius, Cyperus squarrosus, Dalea bicolor y Drymaria holosteoides, preferidas por el borrego cimarrón (^{Fish y Wildlife Service 2000}; ^{Martínez y Galindo 2001}). Sin embargo, se ha registrado que durante el periodo de sequía los borregos prefieren el consumo de suculentas (^{Warrick y Kraumsan 1989}; ^{Tarango et al. 2002}).

A causa de su tamaño y distribución dispersa, las especies vegetales de herbáceas y cactáceas no son detectadas por el sensor TM del satélite Landsat. En consecuencia, se trata de una limitante del método para investigar aspectos de ecología del forrajeo, e. g. la identificación de especies palatables para el borrego cimarrón.

En el periodo de lluvias, en la localidad de El Zamora, los grupos de añeros y hembras seleccionaron la categoría de NDVI alto, que indica mayor disponibilidad de forraje. Se ha considerado como posible que la disponibilidad del forraje (estimada por NDVI) sea el factor que más influye en la distribución espacial del borrego cimarrón durante el periodo de lluvias (^{Hoglander et al. 2015}).

En Los Hemes sucedió lo contrario: los borregos cimarrones volvieron a seleccionar las categorías de baja cobertura vegetal y evitaron los sitios con alta cobertura vegetal que estaban localizados en sitios planos. Lo anterior contrasta con lo reportado en las poblaciones de borrego cimarrón en Sonora, donde los machos, durante el periodo de lluvias, seleccionan los sitios con mayor cobertura vegetal en lomeríos y zonas poco escarpadas. Lo anterior aparentemente los expone a sus depredadores (^{Tarango 2000}).

En el periodo de poslluvias no encontramos una relación entre la distribución de los borregos cimarrones y la cobertura vegetal estimada por NDVI. Esto puede deberse a que el alimento se halla disponible de forma homogénea y, por consiguiente, no hay una selección por alguna de las categorías de NDVI; e. g. la distribución del venado bura en los desiertos de California tras el periodo de lluvias no se correlaciona con una mayor disponibilidad de forraje estimada por el NDVI (^{Marshal et al. 2006}).

Al parecer, la selección que realiza el borrego cimarrón de la cobertura vegetal en Sierra Santa Isabel no se relaciona con los sitios con alta cobertura vegetal, que indican una mayor disponibilidad de alimento. Lo anterior tiene sentido, ya que se ha encontrado que los borregos tienden a evitar sitios con una alta cobertura, pues obstruye su visibilidad y no permite detectar a sus depredadores en la distancia (^{Valdez y Krausman 1999}; ^{Tarango et al. 2002}).

El NDVI, a diferencia de otros índices, tiene la ventaja de ser compatible con diferentes sensores, lo que facilita la investigación en diferentes escalas (^{Pettorelli 2013}). Por ejemplo, el sensor AVHRR (Advanced Very High Resolution Radiometer), con píxeles de resolución espacial de 1 km², permite explorar la selección que hace un animal del hábitat en una escala geográfica. Con otros sensores, como MODIS (píxeles de 250 m) y Landsat (píxeles de 30 m), es posible analizar la selección del hábitat de una especie en un ámbito local (e. g. la sierra de este estudio).

Esta metodología de evaluación de la preferencia del borrego cimarrón por un hábitat puede replicarse en las diferentes áreas de distribución natural de la especie; principalmente, por la disponibilidad de modelos digitales de elevación provenientes de la misión SRTM para Norteamérica (^{USGS 2010}). La resolución espacial de estos modelos (píxeles de 20 y 30 m) ha sido útil para estudiar el uso-selección que realiza el borrego cimarrón del hábitat en diferentes investigaciones (^{Bangs et al. 2005}; ^{Sappington et al. 2007}; ^{Gilad et al. 2013}; ^{Villepique et al. 2015}). Asimismo, estos modelos y las imágenes de satélite de diferentes sensores se hallan disponibles de forma gratuita (www.glovis.usgs.gov) y es posible obtener imágenes históricas de los últimos 25 años, o bien actuales, como las provenientes del sensor Landsat 8 OLI (Operational Land Image).

Agradecimientos

Agradecemos a las siguientes personas e instituciones: a F. Gomis, E. Espinoza, R. Duarte, J. Sánchez, A. Guevara, I. Tafoya y A. Robles, por su ayuda en campo; al laboratorio de Sistemas de Información Geográfica del CIBNOR, en particular a J. Rivera Rosas; a CONACYT, por la beca (205965) otorgada a J. G. Escobar Flores; al Desert Bighorn Council, por el apoyo económico mediante la beca Hansen-Wells; a los propietarios de las UMAS El Volcán, Santa María y Rancho Tafoya, por permitirnos el acceso a sus predios; a M. Delgado y P. Peña, por sus sugerencias al manuscrito; al editor del resumen en inglés, I. L. Fogel, y a los tres revisores anónimos por sus valiosos comentarios, que ayudaron a mejorar la calidad de este trabajo.