Spanish (pdf)
Article in xml format
Article references
Automatic translation
Send this article by e-mail
Cited by SciELO 
Cited by Google 
Similars in
SciELO 
Similars in Google 
Permalink
Introducción
La simulación de hábitat físico es una herramienta numérica que permite cuantificar dicho hábitat en términos de la profundidad del agua (m), la velocidad (m s-1) y el tipo de sustrato (m), en una descarga específica para un caudal dado [1]. Además, predice el impacto de un cambio en el caudal en la disponibilidad de hábitat para las especies seleccionadas como objetivo [2]-[5]. Estas simulaciones de hábitat físico son un requisito esencial para evaluar la salud [6]-[8], el manejo, el inventario y la restauración de los ríos [8]-[11].
Los modelos de simulación de hábitat físico se desarrollan en la interfaz de dos disciplinas históricamente diferentes: la ecología y la hidráulica, combinando los procesos físicos y su influencia en la ecología acuática y la biología [12]. La necesidad de resolver problemas en el manejo del agua ha impulsado su desarrollo como un campo distinto de investigación aplicada con un interés creciente en todo el mundo [12]-[15]. La simulación de hábitat físico tiene como objetivo lograr una evaluación exhaustiva; los científicos han desarrollado centenares de métodos y técnicas de evaluación del caudal del rio, para diversas configuraciones y con diferentes bases teóricas [3], [16].
Entre los modelos más usados está el unidimensional de simulación hidráulica y de hábitat acuático, Physical Habitat Simulation System (phabsim) [17]. Este ha sido el más usado en estudios de simulación de hábitat físico en el mundo, en especial en Norte América [18]-[20]. Junto a este, existe el modelo en dos dimensiones, River2D [21], que se usa para los estudios de simulación de hábitat acuático [22]-[25].
PHABSIM y River2D son empleados en la Instream Flow Incremental Methodology (IFIM) [26], y quizás son los procedimientos más conocidos y documentados para la planificación de los recursos hídricos [3]. Estos modelos simulan el hábitat idóneo en las corrientes de agua, y los resultados se muestran en valores de hábitat potencial útil (HPU), que representa el área utilizable de hábitat idóneo para las especies acuáticas en los tramos de interés [27].
En México, en 2012, se publicó la Norma Mexicana NMX-AA-159-SCFI-2012, que establece el procedimiento para la determinación del caudal ecológico en cuencas hidrológicas [28]. La norma considera la simulación de hábitat fisico para definir el caudal mínimo ecológico, a partir de un estudio de todos los factores y condicionantes del hábitat para una o varias especies objetivo del ecosistema fluvial [29]. Sin embargo, la relación entre los conocimientos de hábitat físico (profundidad, anchura del cauce, pendiente de las orillas, velocidad del agua, tipo de sustrato, vegetación acuática y sombreado), los biotopos (el cauce, la ribera y la llanura de inundación), los parámetros físicos-químicos (temperatura, pH, conductividad, sólidos disueltos y oxígeno disuelto) y los ecológicos (alimentación, reproducción, competencia, crecimiento, entre otros), o la caracterización de los hábitats funcionales en México y Latinoamérica, presentan un desarrollo incipiente.
El río Ayuquila-Armería es uno de los ríos más importantes en el occidente mexicano; el agua que fluye en él es almacenada en cuatro grandes embalses: Tacotán, General Ramón Corona Madrigal (Trigomil), Basilio Badillo (Las Piedras) y El Nogal. Esta retención sirve para el aprovechamiento de los distritos de riego 094 en el Estado de Jalisco y 053 en el Estado de Colima, donde, durante la temporada de lluvias, se minimiza el flujo natural y se almacena para el riego del siguiente año [30]. El desvío y la retención han afectado la biota del río [29], [31] y [32].
Considerando las obligaciones normativas mencionadas y la gran importancia del río Ayuquila-Armería, el objetivo de este estudio es generar las simulaciones de hábitat físico, utilizando las curvas de idoneidad de hábitat de las especies: Nandopsis istlanum (Jordan y Schneider, 1899), Astyanax cf. aeneus (Günther, 1860), Poecilia butleri (Jordan, 1889) y Ilyodon furcidens (Jordan y Gilbert, 1882), que se encuentran dentro del río Ayuquila-Armería. La combinación de la hidráulica del río y las curvas de idoneidad determinarán el HPU para cada caudal simulado, además de relacionar las simulaciones con las series de tiempo del régimen hidrológico actual, utilizando el software International River Interface Cooperative (iRIC 1.5) en su plataforma FaSTMECH, y siguiendo la metodología IFIM. Los resultados conseguidos se utilizan en la estimación del caudal mínimo ecológico dentro del río para las especies objetivo seleccionadas. La información lograda proporcionará datos confiables a los gestores en la toma de decisiones, para un mejor uso y manejo de los recursos hídricos del río.
Materiales y métodos
Tramo de estudio
El tramo para la simulación de hábitat físico se ubicó a 30 km aguas abajo de la confluencia de los ríos Ayuquila y Tuxcacuesco: en las coordenadas 19° 28' 56,68'' (N) y 103° 50' 08,64'' (O) (Figura 1), y a una distancia aguas abajo de las presas: Trigomil y Las Piedras con 110 y 85 km, respectivamente. La distancia de las presas y la dinámica hidrológica que presenta en el tramo seleccionado se asocian a un régimen natural [26]. El tramo seleccionado limita entre los estados de Jalisco y Colima en el occidente de México, la cuenca vertiente drena una superficie aproximada de 9,864 km2 con un trayecto de 321 km desde su nacimiento hasta su desembocadura [29]. El caudal es extremadamente variable, con 10 m3 s-1 en la época seca del año y valores por arriba de los 160 m3 s-1 en la época de lluvias, con un flujo permanente durante todo el ano. El clima es semiárido, con una precipitación promedio anual de 836 mL y una temperatura promedio de 21 °C [33].
Selección del tramo de estudio
El tramo seleccionado tiene una longitud de 706 metros, incluyendo las diferentes unidades morfodinámicas (pozas, rápidos, rápidos someros, remansos, tablas e isletas) características del río. Con el manejo del equipo topográfico (estación total Sokkia 630 R y GPS Astech Promark 3), se levantaron 62 secciones de río, con una distancia aproximada entre secciones de 10 metros (Figura 1), con fecha de referencia del 12 al 14 de junio 2013. Cada sección consistió en un transecto perpendicular a la dirección de flujo, partiendo del punto de mayor inundación de un lado del río hasta el otro [29]. En cada una de las secciones del río se instauraron puntos de medición de los parámetros: profundidad de la columna de agua (m), velocidad (m s-1) y granulometría del sustrato (mm); este último utiliza la técnica Pebble Count en transectos [34]. Estas mediciones se establecieron considerando los cambios trasversales de la pendiente, la profundidad, la velocidad y el tipo de sustratos, cada punto quedó en coordenadas utm WGS84 13 Norte.
Para estimar el caudal circulante dentro del sitio de estudio, en cada punto donde se estimó la altura de la lámina de agua, a diferentes distancias en la sección transversal, se definieron unas celdas donde se midió la velocidad del agua, a una profundidad de la superficie de 0,6 veces del calado [29]. Para el correcto ajuste del modelo hidráulico, se realizaron dos mediciones del caudal circulante: una sección aguas arriba del sitio de trabajo y otra en la salida la sección aguas abajo del sitio del tramo en estudio) [29]. Se usó el método de velocidad-área para estimar el caudal circulante [35].
Simulación hidráulica
El software iRIC 1.5 implementa la dinámica general de fluidos, y fue desarrollado originalmente por la National Aeronautics and Space Administration (nasa), para el modelado estable y rápido en 2D, más la dinámica de fluidos complejos [23] y [36]. Este software hidrodinámico, en su plataforma conocida como FaSTMECH (Flow and Sediment Transport with Morphological Evolution of Channels), utiliza el modelado del caudal en dos dimensiones (2D), sus tiempos de ejecución son rápidos y su uso es fácil; las ecuaciones que utiliza se emiten en un sistema de coordenadas curvilíneas ortogonales, que se ajustan al canal del río y se describen solo por un trayecto de línea central, ancho y con resolución espacial; la ruta se elige en una línea recta y la cuadrícula o malla es cartesiana [37] y [38].
Parámetros de entrada y condiciones límite
Como dato de partida, iRIC 1.5 usa la representación del Modelo Digital del Terreno (mdt), que se obtiene de la topografía detallada del sitio de estudio. El mdt es almacenado en un formato de red irregular triangulada (tin); a este se incorporaron los datos topográficos y batimétricos en el mdt, para complementar los datos originales y verificar la “superficie” resultante (topográfica y batimétrica) [38] y [39].
La plataforma FaSTMECH dentro de iRIC 1.5 requiere definir tres condiciones para obtener una solución válida y estable. Estos valores dependen de las características hidráulicas del cauce, que son representadas por el dominio del modelo. Los valores óptimos se resuelven por lo general a través de un proceso de prueba y error [40].
La primera condición fue la estimación del caudal circulante en m3 s-1 aguas arriba, donde inicia el tramo de estudio, y se obtuvo de los datos observados en las campañas de campo. La segunda condición fue la elevación de la superficie de la lámina de agua a la salida del tramo de estudio (aguas abajo), que se consiguió del levantamiento topográfico del río a una elevación en metros sobre el nivel del mar.
La tercera condición de entrada corresponde al valor n Manning, que generalmente varía entre 0,030 a 0,038 [41]. En este estudio la rugosidad del canal fue representada por coeficientes de arrastre adimensionales y se calculó con base en el material de lecho, como el tamaño medio del sustrato encontrado en el tramo de estudio; por tanto, no se ajustó durante la calibración [38] y [39].
Calibración y validación
Con la información del levantamiento topográfico, fueron definidos 62 puntos de elevación de la superficie de la lámina de agua, identificados como el Thalweg a lo largo del tramo de estudio (Figura 1). Con la información del caudal circulante observado y asociado a la elevación de la superficie de la lámina de agua del Thalweg es posible comparar el modelo predicho [42] y los parámetros del modelo: caudal circulante observado (m3 s-1), elevación de la superficie de la lámina de agua (m) encontrada a la salida del tramo de estudio y la rugosidad constante (m). Dentro de FaSTMECH en el módulo Relaxation Parameters, para calibrar el modelo se ajustó la respuesta de la elevación de la superficie de la lámina de agua observada vs. la calculada por el programa, y que sean coincidentes para los 62 puntos del Thalweg, hasta obtener una solución satisfactoria.
Para validar el modelo de simulación, se usó el módulo Verification dentro de la plataforma FaSTMECH, y se obtuvo el error cuadrático medio (ECM) de la elevación de la superficie de la lámina de agua (m) calculada por el programa vs. observada en las campañas de campo. Al mismo tiempo, se empleó el coeficiente de determinación (R2) para medir la asociación del modelo calculado con el observado. Con los valores resultantes del ECM y R2, de acuerdo con la calidad de asociación, el modelo será validado. Después de aprobarse el modelo hidráulico, se generan las diferentes simulaciones a diferente caudal a lo largo del tramo de estudio, para las especies seleccionadas y para un horizonte de tiempo de 31 años (1971-2001) del régimen hidráulico actual.
Selección de las curvas de idoneidad
Una vez medidas las características hidráulicas del cauce y calibrado el modelo de simulación hidráulico, fue necesario seleccionar las especies indicadoras del tramo de estudio. Para esto se llevó a cabo un análisis de las bases de datos del programa de monitoreo bimensual de largo plazo (1998-2011), de calidad de agua y peces dentro del río Ayuquila-Armería [33], del proyecto regional Manejo y conservación de los recursos naturales del río Ayuquila-Armería, Departamento de Ecología y Recursos Naturales de la Universidad de Guadalajara (comunicación personal, Luis Manuel Martínez). Las especies seleccionadas como indicadoras del tramo de estudio fueron: N. istlanum, A. cf. aeneus, P. butleri y I. furcidens.
Las curvas de idoneidad se consiguieron a través de estadísticas con la técnica de histograma de frecuencia para cada variable de interés y para cada especie; por tanto, los datos de cada variable se clasificaron en grupos de intervalos siguiendo los pasos establecidos por [26] y [43] (Figura 2).
Fig. 2. Curvas de idoneidad de hábitat acuático, para las variables de profundidad, velocidad y tipo de sustratos
Los tres factores escogidos, velocidad, profundidad y tipo de sustratos, fueron integrados para relacionar las características físicas del tramo fluvial seleccionado, con las necesidades de hábitat de las especies que se eligieron.
Simulación del hábitat
Para simular el hábitat físico se empleó la aplicación Habitat Builder, que se encuentra dentro de la plataforma FaSTMECH, y se usó la información de las curvas de idoneidad mencionadas anteriormente.
El módulo Habitat Builder permite el cálculo de las estimaciones de área utilizables para las variables independientes: profundidad, velocidad y tipo de sustratos. Para cada variable se usó el método de interpolación tin (red de triángulos irregulares). La idoneidad final de esta celda se obtiene promediando (p.e. media geométrica) la idoneidad interpolada de las curvas de idoneidad [38]. Habitat Builder define una función de hábitat compuesta por uno o más predictores (curvas de idoneidad) para mapear el valor de la función de hábitat a la celda del modelo, y calcula un HPU que consiste en la suma para cada caudal evaluado del área o pixel de simulación hidráulica ponderado por la idoneidad conjunta [38].
Criterios de idoneidad del hábitat
Para producir una medida representativa del hábitat físico disponible en función de la descarga (caudal), los parámetros (profundidad, velocidad y tipo de sustratos) fueron combinados en la plataforma FaSTMECH, con la preferencia de las especies. Esta combinación genera el Índice de Aptitud Combinada (iac) y deduce que los parámetros hidráulicos, mencionados anteriormente, son variables independientes que caracterizan el hábitat, los cuales se usan comúnmente para describir la calidad del hábitat a escala del tramo de río seleccionado, y tienen la expresión de la ecuación (1) [40], [44] y [45]:
Donde n es el operador del producto, los subíndices i y j indican el i-th celda y el j-th el parámetro físico. El valor de iac varía espacialmente con la descarga.
El Hábitat Potencial Útil (HPU) y el Hábitat Potencial Hidráulico (hph), se deriva del iac y comúnmente utilizado para describir la cantidad del hábitat [46-47]. Tanto HPU como hph son funciones de descarga y utilizan las siguientes expresiones (2), (3) respectivamente:
Donde p es el número de celdas dentro del área mojada Aw del tramo del río y Ai es el área de la celda i-th. El HPU se expresa en unidades de área de microhábitat por distancia unitaria a lo largo del río (p.e., m2 por 1000 m) [40], [48].
Caudal simulado/hábitat potencial útil
Después de calibrarse el modelo hidráulico, se llevaron a cabo las simulaciones de caudal a: 0,2. 0,3. 0,4. 0,5. 0,8 y 1, 2, 3, 4 hasta alcanzar los 20 m3 s-1. Para las simulaciones de caudal se obtuvo el HPU para cada una de las especies seleccionadas.
Simulaciones del régimen hidráulico actual
Para generar el modelamiento hidráulico de caudales se introdujo una serie de caudales observados y se modeló de forma dinámica un horizonte de tiempo para 31 años (1971-2001); esto se hizo a través de la información de la estación hidrométrica 16032 Coliman, localizada en las coordenadas geográficas 18° 56' 25,26'' N y 103° 56' 39,72” O. Los valores promedio a nivel mensuales de la serie de caudales estudiados fueron simulados para obtener el HPU para cada una de las especies seleccionadas. Los datos fueron extraídos del Banco Nacional de Aguas Superficiales (Bandas) de la Comisión Nacional de Agua (Conagua) [49].
Resultados y discusión
El resultado del levantamiento topográfico en el río fue de 1,055 puntos distribuido en los 706 metros del tramo de estudio (Figura 1), mientras que el del caudal circulante observado para el tramo de estudio fue de 11,50 (m3 s-1), con una elevación de la superficie de la lámina de agua (ESLA) de 618,65 m s. n. m. a la salida del tramo de estudio (aguas abajo). Las características encontradas del lecho del río fueron las siguientes: madera en descomposición 1 %, limos y arcillas 1 %, arenas 16 %, gravas 12 %, cantos rodados 11 % y rocas 59 %.
La clasificación del tamaño del sustrato resultante fue D16 0,98 mm, D35 120 mm, D50 680 mm, D65 1200 mm, D84 1600 mm, D95 1900 mm, con un promedio de 39,6 mm. La rugosidad del canal, representada por el coeficiente de arrastre adimensionales, fue estimada en 0,0396 m, con un sesgo (skewness) de -0,79. Este valor negativo indica que los valores estarán más reunidos en niveles superiores a la media aritmética para el tramo de estudio. El tamaño del sustrato fue un aspecto importante para la simulación de hábitat, ya que definió el hábitat físico para los peces. En el tramo las partículas >1000 mm se presentaron en los transectos T1 a T20 con pendientes promedio de 1,17 %, en los transectos T21 a T40 con pendientes de 0,05 % dominados por limos, arcillas y arenas, y los transectos T41 a T60 con pendientes 0,01 % con presencia de arenas, gravas y gravillas (Figura 1).
Calibración
El resultado de la relación entre ESLA (puntos del Thalweg) y el caudal circulante presentó un buen ajuste entre los valores calculados por el modelo y los observados en las campañas de campo (Figura 3a), con un coeficiente de determinación de 0,9992 (Tabla 2). En el análisis de asociación el ECM entre ESLA calculada y observada fue de 0,040 metros (Tabla 2), lo que muestra un error de ESLA calculada cercano a cero metros (Figura 3b).
Los resultados de evaluación de las diferentes simulaciones mostraron que al calibrar el modelo a 11,50 (m3 s-1), las simulaciones generadas a diferente caudal (altos y bajos) expusieron un buen ajuste de acuerdo con R2, con valores por arriba de 0,90 en la mayor parte de las simulaciones (Tabla 1). Las menores asociaciones encontradas fueron para los caudales inferiores a 1 (m3 s-1) y los superiores a 200 (m3 s-1).
Tabla 1. Resultados de las simulaciones de caudal, el ECM de ESLA y el resultado de correlación entre ESLA calculada y observada
|
Caudal simulado (m3 s-1) |
ECM de ESLA (m) |
R2 |
Caudal simulado (m3 s-1) |
ECM de ESLA (m) |
R2 |
|---|---|---|---|---|---|
|
0,1 |
0,789 |
0,6201 |
11,49 caudal observado aguas arriba |
0,044 |
0,9990 |
|
0,2 |
0,768 |
0,6242 |
11,50 caudal observado aguas abajo |
0,040 |
0,9992 |
|
0,3 |
0,732 |
0,6306 |
12 |
0,088 |
0,9981 |
|
0,4 |
0,721 |
0,6348 |
14 |
0,109 |
0,9989 |
|
0,5 |
0,716 |
0,6511 |
16 |
0,135 |
0,9989 |
|
0,8 |
0,712 |
0,6525 |
18 |
0,119 |
0,9922 |
|
1 |
0,697 |
0,6540 |
20 |
0,141 |
0,9930 |
|
3 |
0,557 |
0,8209 |
40 |
0,532 |
0,9611 |
|
4 |
0,556 |
0,8489 |
60 |
0,731 |
0,9315 |
|
6 |
0,138 |
0,9836 |
80 |
0,937 |
0,9073 |
|
8 |
0,081 |
0,9936 |
100 |
1,210 |
0,9310 |
|
10 |
0,283 |
0,9928 |
200 |
1,937 |
0,7648 |
Fig. 3. a. asociación entre esia calculada por el software /s. observada en las campanas de campo, y b. ECM de ESIA entre el modelo calculado por el software y el observado en las campanas de campo
Es importante mencionar que el buen ajuste encontrado entre los caudales calculados con los observados también exhibió un aumento del ECM de ESLA, conforme aumentan o disminuyen las simulaciones con respecto al caudal calibrado. Las simulaciones inferiores a 4 (m3 s-1) mostraron un ECM de ESLA por arriba de los 0,5 m, mientras que las simulaciones superiores a 100 y 200 (m3 s-1) presentaron un ECM de ESLA por arriba de 1,2 m (Tabla 1). Es necesario para futuras investigaciones disenar campanas de campo para validar los caudales muy bajos y altos, dentro del modelado hidráulico.
Simulación de hábitat físico
La relación de caudal y hábitat presenta diferentes formas y valores absolutos de HPU entre las diferentes especies (N. istlanum; A. cf. aeneus; P. butleri y I. furcidens, Figura 5). La cantidad de aptitud del hábitat disminuye con el aumento de los caudales. La distribución horizontal calculada mostró las mayores velocidades y profundidades en el centro del canal y tiende a disminuir hacia las orillas del río. Además, se encontró que en los caudales simulados el aumento HPU se produce cerca del banco derecho e izquierdo del río, y solo para la especie I. furcidens exhibió un aumento de HPU que se extiende en el tramo seleccionado, a caudales > 4 (m3 s-1) (Figura 4). La superficie de HPU para las especies de peces modelados fue sobre-predichos a flujos más bajos de 4 (m3 s-1).
La distribución del HPU obtenida de las simulaciones a lo largo del tramo de estudio mostró que en aguas arriba donde inicia el tramo de estudio con pendientes de 1,60% se tiene la condición de hábitat más pobre y el aumento de HPU para las especies se mueve hacia aguas abajo del tramo de estudio, donde se mejorará la condición de aptitud del hábitat con pendientes de 0,04 y 0,01 % para las especies seleccionadas (Figura 4).
Un hallazgo importante en la simulación de hábitat físico para las cuatro especies fue la evidencia de distintas selecciones de hábitat en términos de profundidad, velocidad y tipo de sustratos.
Hábitat Potencial Útil vs. caudal simulado
Las curvas (caudal/HPU) mostraron un gradiente más elevado de HPU para P. butleri y I. furcidens, pero tal diferencia fue menos notable para N. istlanum y A. cf. aeneus. Las curvas producidas presentaron, en general, una misma tendencia (Figura 5). En los resultados se observó que los caudales simulados <1 (m3 s-1) no producen suficiente representación de HPU; por su parte, las simulaciones entre 1 y 6 (m3 s-1) evidenciaron un crecimiento rápido de HPU para las cuatro especies, y las simulaciones de 6 a 16 (m3 s-1), las especies N. istlanum; A. cf. aeneus y I. furcidens, muestran un crecimiento lineal con relación a los caudales.
La especie P. butleri alcanzó su máximo HPU en los caudales simulados a 6 y 7 (m3 s-1) y I. furcidens a 6 (m3 s-1). Estas dos especies de la familia Poeciliidae presentaron similitudes de caudal a la especie G. atripinnis, gracias a un estudio realizado en la cuenca baja del río Duero en México por [50]. Las especies N. istlanum alcanzaron su máximo HPU a 8 (m3 s-1) y A. cf. aeneus a 15 (m3 s-1) (Figura 5). Las variables hidráulicas modeladas (profundidad, velocidad y tipo de sustrato) cuando poseen un alto HPU para una determinada especie se debe al tamano de la celda o pixel, resultado de las secciones transversales del levantamiento topográfico. Las cuatro especies estudiadas mostraron diferentes variaciones en HPU que dependen del caudal.
Simulaciones del régimen hidráulico actual
La serie de tiempo (1971-2001) entre el caudal mes/HPU para las cuatro especies seleccionadas muestra que entre enero y mayo es el periodo de caudales decrecientes dentro del río, marcado por la estación seca. Las condiciones de aptitud de hábitat para las cuatro especies se analizan desde una evaluación general, considerando que, ordinariamente, durante el periodo seco en la primavera, el HPU es óptimo porque representa la condición de hábitat adecuada para la fecundación, el desarrollo embrionario, el nacimiento y el crecimiento de las especies [51] y [52].
Para las cuatro especies, estas fases las llevan a cabo entre marzo y mayo; por lo tanto, la primavera será un periodo de hábitat sensible para la etapa de la vida de estas especies [53]-[55] (Figura 6). Solo la especie N. istlanum presentó un HPU muy pobre para esta época, ocasionado por la retención de agua por las presas, lo que acerca los caudales del régimen hidrológico actual a un valor cero de caudal. Lo contrario ocurre en la época seca, en la época estival (lluvias), temporada cuando se presentan los mayores caudales de julio a octubre; en estos meses se muestra una disminución del HPU para las cuatro especies seleccionadas (Figura 6), hecho que confirma que un mayor caudal no representa un aumento de HPU. Estudios anteriores del río describen este comportamiento hidrológico de acuerdo con [26]. Es importante destacar que la toma de data para la construcción de las curvas de idoneidad de las especies seleccionadas fue conseguida, en gran parte del tiempo, durante la temporada de estiaje, cuando el río presenta menor caudal; esto ocasiona una incertidumbre al simular el HPU de los caudales altos.
Fig. 6. Distribución promedio mensual del régimen hidrológico actual (1971-2001) y su relación con las simulaciones de caudal y su Hábitat Potencial Útil para las especies
Esta serie temporal de julio a octubre muestra el decaimiento de HPU e ilustra que habrá un estado de hábitat muy pobre para las cuatro especies [29]. Las relaciones entre las series temporales y el HPU son una parte fundamental para los estudios hidrobiológicos, y es necesario establecer las series temporales de HPU en condición de referencia de caudal natural en comparación con diferentes escenarios, como la regulación del caudal ocasionada por la retención de las presas, la extracción de agua y el cambio climático [26] y [56]. Los caudales simulados de las series temporales presentados en este estudio servirán para apoyar a los gestores en la toma de decisiones para administrar de mejor manera los recursos hídricos del río Ayuquila-Armería.
El sistema analítico iRIC 1.5 constituye una herramienta provechosa para el ordenamiento y manejo de los recursos hídricos. El modelado bi-dimensional es potencialmente capaz de simular patrones de flujos complejos de acuerdo con [18], [57] y [58] y proporciona información detallada con respecto al patrón de flujo, la velocidad y la tensión de arrastre del lecho del río. Pero también presentó una serie de sesgos en la simulación de hábitat, por ejemplo, una vez calibrado el software a un caudal circulante (obtenido de la campaña de campo), se intuye (como verdad) que con la misma calibración de este caudal circulante se pueden estimar los caudales circulantes (inferiores 1 m3 s-1) y los extremos (mayor 100 m3 s-1); esto ocasiona una serie de incertidumbre en los resultados de la simulación de hábitat. La toma de datos extremos está relacionada con aspectos de seguridad personal (mayor profundidad y velocidad en el río, que lo hace inaccesible), lo cual imposibilita la toma de datos.
Las simulaciones de hábitat físico en este estudio, de acuerdo con la IFIM, mostró una eficaz relación entre los parámetros físicos del hábitat en diferentes escalas y simulaciones de caudal. Esta relación funcional entre el caudal y el HPU permitió la integración del hábitat acuático, así como la evaluación del efecto potencial del régimen de caudal, utilizando series temporales [59]. Estos resultados dentro de la IFIM no son extrapolables, ya que se han utilizado como indicador del tramo seleccionado, de manera que no es válido para otros tramos o estudios de simulación de hábitat físico dentro del río. Es necesario llevar a cabo estudios de transferibilidad de los criterios de idoneidad, cuantificando el hábitat hidráulico y utilizando el método de [60], el cual puede emplearse para determinar la calidad y la cantidad de hábitat en otra corriente o tramos seleccionados.
En la actualidad, en México y en Latinoamérica la aplicación de este tipo de métodos no es común y no existen simulaciones de hábitat físico para las especies seleccionadas, que podrían ser utilizadas para establecer comparaciones.
Conclusión
Este es uno de los primeros estudios que identifica especies de peces neotropicales para su uso en estudios relacionados con la simulación de hábitat físico. Esto amplía el uso de modelos de predicción y la capacidad de predecir cuantitativamente los impactos biológicos de los cambios en las variables hidráulicas en ambientes acuáticos.
La simulación de hábitat físico aparece en México y en Latinoamérica como un campo nuevo de investigación que debe fomentarse, ya que proporciona soluciones razonadas para salvaguardar la destacada biodiversidad de sus ecosistemas acuáticos, que se encuentran en proceso acelerado de degradación. Esta propuesta de simulación de hábitat físico permitirá a los gestores de decisiones en el río Ayuquila-Armería —como lo indica la Norma Mexicana NMX-AA-159-SCFI— proteger y conservar los ecosistemas acuáticos desde un punto de vista más ecológico.
El software iRIC 1.5 en su plataforma FaSTMECH fue adecuado en este estudio para simular las características del hábitat físico en dos dimensiones, lo que permitió cambiar el tamaño de la malla y variar la parametrización del modelo, y así mejorar los resultados de la simulación dentro del río Ayuquila-Armería y presentar un ajuste a las variables hidráulicas en el proceso de calibración. Las elevaciones de la lámina de agua observada de las 62 secciones transversales, que se obtuvieron del levantamiento topográfico y la combinación de las curvas de idoneidad propuestas, permitieron con mayor detalle describir las simulaciones del hábitat físico dentro del tramo seleccionado.
Un aporte importante en este estudio fueron las curvas de idoneidad de las especies N. istlanum; A. cf. aeneus; P. butleri y I. furcidens; estas no limitarán el uso de modelos de predicción y la capacidad de predecir de manera cuantitativa los impactos biológicos por la carencia de información, y sí servirán para futuras simulaciones de hábitat físico a lo largo de las vertientes (Atlántico y Pacifico). N. istlanum se localiza en el occidente mexicano, A. cf. aeneus utiliza ambas vertientes de las Américas y las Indias Occidentales; desde Veracruz en México hasta Costa Rica y Colombia. El P. butleri utiliza la vertiente del Pacífico en México y Mesoamérica, desde el río Fuerte, en Sonora, hasta la boca del río Comasagua, en El Salvador. I. furcidens se localiza en la vertiente del Pacífico, en los ríos Ameca, Armería, Marabasco y Coahuayana.
Para finalizar, estudios como este exponen las ventajas de seleccionar la técnica de simulación de hábitat físico, que se establece en la metodología de la Norma Mexicana NMX-AA-159-SCFI-2012 en su apéndice normativo E, para salvaguardar los sistemas de ríos en México.
Agradecimentos
Un agradecimiento al proyecto Manejo y conservación de los recursos naturales del río Ayuquila-Armería del Departamento de Ecología y Recursos Naturales de la Universidad de Guadalajara, por su acceso a la información. Se agradece a las aportaciones críticas del referee encargado de la evaluación de este artículo.
