Las lagunas costeras del Estado de Veracruz en el Golfo de México se caracterizan por su alta biodiversidad y valor biológico, e incluyen una gran variedad de recursos pesqueros, ocupando el cuarto lugar en la producción pesquera nacional (Lara-Domínguez et al., 2011). En particular, el sistema lagunar La Mancha declarado sitio RAMSAR (Ficha Informativa de los Humedales de Ramsar, 2003) es pequeño y somero; presenta una conexión efímera con el mar la cual se cierra durante un periodo del año y se abre durante la época de lluvias, lo que permite el intercambio de organismos, nutrientes y agua. Este ecosistema lagunar, destaca por su alta producción primaria y riqueza de recursos biológicos, y es utilizado por una gran variedad de peces de manera temporal o permanente para llevar a cabo parte de sus ciclos de vida, como área de alimentación, crianza, desove y refugio, al mismo tiempo, representan especies para consumo local (Lara-Domínguez, Day, Yáñez-Arancibia, & Sainz-Hernández, 2006; Díaz-Ruiz, Aguirre-León, Juárez, Matsumoto, & Hernández, 2008). De esta manera, la utilización de los hábitats en estos sistemas por los peces, está relacionado con las variaciones de salinidad y temperatura, la vegetación circundante, la descarga de los ríos y el intercambio de agua con el mar, entre otras características, lo que influye en la estacionalidad de la diversidad, distribución, abundancia y frecuencia de las especies que habitan este sistema (Able, 2005; Barletta, Barletta-Bergan, Saint-Paul, & Hubold, 2005; Díaz-Ruiz, Aguirre-León, Calva-Benítez, & Barba-Macías, 2012; Zárate-Hernández, Castillo-Rivera, Sanvicente-Añorve, & Ortiz-Burgos, 2012; Aguirre-León, Pérez-Ponce, & Díaz-Ruiz, 2014).
Adicionalmente, los peces funcionan como reguladores energéticos en la comunidad, al transformar, intercambiar y almacenar la energía, a través de complejas interacciones biológicas (e.g.reproducción, desove, alimentación, crecimiento desde el mar hacia los hábitats lagunares), que están relacionadas con los patrones de migración y el uso diferencial de los hábitats, como respuesta a la heterogeneidad ambiental del sistema (Deegan et al., 1986; Castillo-Rivera, Zavala-Hurtado, & Zárate-Hernández, 2002; Arceo-Carranza, Vega-Cendejas, Montero, & Hernández de Santillana, 2010; Aguirre-León et al., 2014).
No obstante, a pesar de la importancia ecológica y económica de los recursos pesqueros de la laguna La Mancha, la información sobre la estructura ecológica de la comunidad de peces en este sitio es escasa. Por lo tanto, el conocimiento sobre los procesos ambientales y su influencia en la estructura ecológica de los peces de este sistema lagunar es relevante, toda vez que no existen antecedentes de estudios a largo plazo en esta región. Los estudios en este sistema sólo incluyen listas de especies, como el de Reséndez-Medina y Kobelkowsky (1991) y sobre evaluación ecológica de Díaz-Ruiz et al. (2008). Por lo tanto, el presente estudio representa la primera información sobre la ictiofauna del sistema lagunar con un enfoque de estructura ecológica. Para este fin, se plantearon los siguientes objetivos: 1) analizar el comportamiento fisicoquímico del sistema en escala espacial y temporal, 2) determinar la composición de la ictiofauna, 3) evaluar los principales patrones de diversidad, abundancia y conjuntos de peces, con relación a las principales variables fisicoquímicas, y 4) definir la dominancia, componentes comunitarios y niveles tróficos de las especies de peces.
Materiales y métodos
Área de estudio: La Laguna de La Mancha se localiza en la planicie costera central de Veracruz, en el Golfo de México (19°34’-19°36’ N & 96°22’-96°24’ W). Tiene una superficie de 1.5 km2y está integrada por dos cuerpos de agua unidos por un canal estrecho (Figura 1). Es una laguna somera con profundidad promedio de 0.60 m en el norte a 2.5 m en el sur. Tiene influencia marina a través de una boca de conexión efímera que se abre en época de lluvias, y se cierra durante la época seca. Por la boca cruza un gasoducto, lo que ha modificado el intercambio de agua con el mar (Lara-Domínguez et al., 2006). Recibe la descarga de los ríos Caño Grande y Sábalo en el extremo sur (Figura 1). El clima es cálido subhúmedo con lluvias en verano; la temperatura media anual fluctúa entre 22 y 26 °C y la precipitación entre 1 200 y 1 500 mm anuales (García, 2004). Se reconocen tres épocas climáticas: secas (marzo a junio), lluvias (julio a octubre) y nortes (noviembre a febrero) (Contreras-Espinosa, Castañeda-López, & Rivera-Guzmán, 2006). Existen diferentes ambientes acuáticos en el sistema con vegetación circundante de manglar (Rhizophora mangle), sumergida (Rupia maritima, Siringodium filiforme) y áreas de alta sedimentación (Lara-Domínguez et al., 2006).
Actividades de campo y laboratorio: La recolección de peces se realizó en siete estaciones de muestreo durante dieciséis meses, entre abril 2005 y enero 2008 (Figura 1). Los muestreos se realizaron cada tres meses en el 2005, y cada 2 meses en 2006 y 2007, incluyendo enero del 2008. Las capturas fueron diurnas (09:00-15:00 horas), y en cada estación se efectuó un arrastre con un chinchorro playero de 50 m de largo, (caída de 2 m, con luz de malla de 1.9 cm en el copo y 1.5 cm en las alas), y un área de 400 m2. Cada arrastre se llevó a cabo entre 20 a 30 minutos de maniobra, con un total de 95 capturas. Todas las muestras fueron fijadas en formaldehído al 4 % neutralizado con borato de sodio, y posteriormente preservadas en etanol al 70 % para su análisis en el laboratorio.
En cada muestreo, se registró la salinidad, temperatura y oxígeno disuelto, con un multianalizador YSI-3800; el pH con un potenciómetro Hanna-9812; la profundidad con una sonda graduada (cm) y la transparencia con un disco de Secchi (cm).
En el laboratorio se midió la longitud total (LT) de los peces con un ictiómetro graduado (cm) y se pesaron con una balanza digital Sartorius-Entris 2202-1S (precisión 0.01g). La identificación taxonómica de los peces se realizó con claves del Atlántico occidental y del Golfo de México (Fischer, 1995; Castro Aguirre, Espinoza-Pérez, & Schmitter-Soto, 1999). El estatus taxonómico se basó en el criterio de Nelson (2006) y para corroborar la autoridad y nomenclatura actual de las especies, se utilizó FishBase (Froese, & Pauly, 2017).
Se calcularon los parámetros ecológicos de la comunidad: diversidad H´ (Shannon, & Weaver, 1963), riqueza de especies D (Margalef, 1969) y equitatividad J´ (Pielou, 1966), así como la densidad (ind./m2), biomasa (g/m2) y peso promedio de los peces (g/ind.) (Aguirre-León et al., 2014). La dominancia de las especies se determinó mediante el Índice de Valor de Importancia (IVI) utilizando el programa ANACOM (De la Cruz-Agüero, 1994). Para determinar el tipo de componente comunitario al cual pertenece cada especie, se siguió el criterio de Yáñez-Arancibia (1986), quien las ubica según su permanencia temporal en ocasional, cíclica o permanente en los sistemas estuarinos: E= estuarino, ME= marino eurihalino, ES= marino estenohalino, D= dulceacuícola. Se definió la categoría trófica de cada especie, analizando la información publicada sobre la alimentación y de acuerdo con el criterio propuesto por Yáñez-Arancibia et al. (1985): a) Consumidores de primer orden, son herbívoros exclusivos, fitoplanctófagos y/o zooplanctófagos, peces detritívoros y omnívoros, b) Consumidores de segundo orden, peces predominantemente carnívoros, pueden incorporar en su dieta restos vegetales y detritus, c) Consumidores de tercer orden, peces exclusivamente carnívoros, donde los vegetales y el detritus es un alimento accidental.
Para estimar las diferencias estadísticas espaciales (entre estaciones) y temporales (entre meses y años) de los parámetros físicoquímicos y parámetros ecológicos de la comunidad, se evaluó previamente la homogeneidad de varianzas entre grupos, con la prueba de Levene (1960); en los casos necesarios, se hicieron transformaciones logarítmicas. Cuando se cumplió la homogeneidad de varianzas se utilizó una prueba ANOVA-F-, cuando no se cumplió, fue utilizada la prueba Kruskal-Wallis-H- (Zar, 1999). Para todas las pruebas anteriores se consideró un nivel de significancia p < 0.05 (Zar, 1999; Hair, Anderson, Tatham, & Black, 2001). Los análisis se realizaron con ayuda del paquete estadístico STATISTICA, versión 10. Se evaluó la relación entre las especies de peces y las variables ambientales del sistema, con un análisis de correspondencia canónica (ACC) (Ter Braak, & Verdonschot, 1995). Para este fin, se utilizó una matriz de abundancia de las especies y de datos físicoquímicos por estación y por mes. Las correlaciones inter-set determinaron las variables ambientales que fueron más importantes en el comportamiento de las abundancias de las especies (McGarigal, Cushman, & Stafford, 2000). El ACC fue generado usando el programa Multivariate Statistical Package (MSVP) versión 3.22.
Resultados
Ambiente-escala espacial: El análisis espacial de los físicoquímicos de 2005 a 2007-2008, mostró gradientes con valores bajos en la zona sur (estaciones 5 a 7), y altos en la zona norte del sistema (estaciones 1 a 4). El valor promedio más bajo de salinidad (11.8 ups ± 12.2) se registró en la estación 7 del 2005 y el más alto (31.3 ups ± 9.1) en la estación 1 del 2007-2008 (Figura 2A a C). Este parámetro sólo indicó diferencias significativas entre estaciones en 2007-2008 (ANOVA-F = 2.7, p < 0.05). La temperatura del agua presentó el valor promedio más bajo (26.0 °C ± 4.10) y más alto (31.6 °C ± 3.9) en la estación 6 durante 2005 y 2007-2008, respectivamente (Figura 2D a F). El ANOVA no señaló diferencias significativas entre estaciones (p > 0.05). Un comportamiento similar a la salinidad lo presentó el oxígeno disuelto, con el valor promedio más bajo (2.4 mg/L ± 1.21) en la estación 7 del 2005 y el más alto (5.7 mg/L ± 0.97) en la estación 1 del 2006 (Figura 2G a I). Este parámetro mostró diferencias significativas entre sitios en 2006 (ANOVA- F = 1.6, p < 0.05) y en 2007 (Kruskal Wallis-H = 7.8, p < 0.05). El valor promedio más bajo de pH (7.1 ± 0.7) se registró en la estación 6 del 2006, mientras que el más alto (7.7 ± 0.6) en la estación 1 del 2005 (Figura 2J a L), sin diferencias significativas entre sitios (p > 0.05). Un gradiente inverso a los parámetros anteriores se observó en la profundidad con el valor promedio más bajo (0.25 m ± 0.1) en la estación 1 y el promedio más alto (1.3 m ± 0.14) en la estación 7 del 2005 (Figura 2M a O). Así, la transparencia más baja (31 %, valor promedio 0.8 m ± 0.29) se registró en zonas profundas (estación 6) del 2005, mientras que la más alta (88 %, valor promedio 0.19 m± 0.05) en zonas someras (estación 1) del 2007-2008 (Figura 2M a O). La profundidad presentó diferencias significativas entre sitios en 2005 (ANOVA-F = 0.05, p < 0.05), 2006 (ANOVA-F = 0.08, p < 0.05) y 2007-2008 (ANOVA-F = 0.03, p < 0.05) y la transparencia en 2005 (Kruskal Wallis-H = 10.3, p < 0.05) y 2007-2008 (ANOVA-F=0.02, p < 0.05).
Ambiente-escala temporal: El valor promedio de salinidad más bajo (3.9 ups ± 2.4) se presentó en julio 2005, mientras que el más alto (34.3 ups ± 3.7) en septiembre del 2007 (Figura 3A a C). Este parámetro mostró diferencias significativas entre meses en 2005 (Kruskal Wallis-H =17.9, p < 0.0001), 2006 (Kruskal Wallis-H = 19.0, p < 0.001) y 2007-2008 (Kruskal Wallis-H = 20.6, p < 0.0001). La temperatura del agua registró el valor más bajo (20.9 °C ± 0.6) en diciembre 2005 y el más alto (33.0 °C ± 0.2) en mayo 2007 (Figura 3D a F). Este parámetro indicó diferencias significativas entre meses en 2005 (Kruskal Wallis-H = 21.9, p < 0.0001), 2006 (Kruskal Wallis-H = 35.5, p < 0.0001) y 2007-2008 (Kruskal Wallis-H = 25.9, p < 0.0001). El oxígeno disuelto mostró un comportamiento inverso a la temperatura con el valor promedio más bajo (3.02 mg/L ± 0.95) en julio 2005 y el más alto (6.6 mg/L ± 1.8) en enero 2008 (Figura 3G a I). Se presentaron diferencias significativas en 2005 (Kruskal Wallis-H = 9.8, p < 0.05) y 2007-2008 (Kruskal Wallis-H = 23.7, p < 0.0001). El valor promedio de pH más bajo se registró en junio 2006 (6.4 ± 0.3), mientras que valores más alcalinos en octubre 2006 (8.3 ± 0.2) (Figura 3J a L). Este parámetro mostró diferencias significativas en 2005 (Kruskal Wallis-H = 22.9, p < 0.0001), 2006 (ANOVA-F = 0.05, p < 0.0001) y 2007 (Kruskal Wallis-H = 38.1, p < 0.00001). Cambios en los niveles de agua se observaron en los meses de estudio con profundidades mínimas promedio (0.6 m ± 0.3) en abril 2006 y máximas promedio (1.3 m ± 0.4) en febrero 2006 (Figura 3M a O). La transparencia presentó el valor promedio más bajo (0.21 m ± 0.09) en julio 2005, mientras que el valor promedio más alto (0.56 m ± 0.26) en febrero 2006 (Figura 3M a O). Sólo se presentaron diferencias significativas entre meses del 2006 para valores promedio de profundidad (ANOVA-F = 0.07, p < 0.001) y transparencia (ANOVA-F = 0.02, p < 0.05). El Cuadro 1 muestra los valores promedios anuales de los parámetros físicoquímicos en el sistema.
Parámetro | 2005 | 2006 | 2007-2008 | Estimador | p |
---|---|---|---|---|---|
Salinidad (ups) | 17.18±9.9 | 22.58±10.3 | 26.86±7.90 | H=16.54 | 0.0003* |
Temperatura (°C) | 27.37±4.2 | 28.50±2.70 | 29.94±3.04 | H =7.08 | 0.030* |
Oxígeno disuelto (mg/L) | 3.85±1.57 | 4.30±1.67 | 5.42±1.32 | F=10.30 | <0.0001* |
pH | 7.45±0.62 | 7.26±0.70 | 7.47±0.48 | H =2.67 | 0.26 |
Profundidad (m) | 0.86±0.39 | 0.78±0.35 | 0.81±0.31 | F =0.43 | 0.65 |
Transparencia (m) | 0.35±0.25 | 0.35±0.20 | 0.42±0.19 | F =1.37 | 0.26 |
Composición ictiofaunística: En el sistema lagunar La Mancha, se recolectaron un total de 5 974 peces pertenecientes a 11 órdenes, 25 familias, 39 géneros y 50 especies (Cuadro 2). Las familias que presentaron mayor número de especies fueron Gerreidae con siete especies: Diapterus auratus, D. rhombeus, Eucinostomus argenteus, E. gula, E. melanopterus, Eugerres plumieri y Gerres cinereus; Engraulidae con cinco especies: Anchoa hepsetus, A. lamprotaenia, A. mitchilli, Anchovia clupeoides y Cetengraulis edentulus; Carangidae con cuatro especies: Hemicaranx amblyrhynchus, Oligoplites saurus, Selene vomer y Trachinotus falcatus; Centropomidae con tres especies: Centropomus parallelus, C. pectinatus y C. undecimalis; y Scianidae con tres especies: Bairdiella chrysoura, B. ronchus y Micropogonias undulatus. El número de especies varió de 37 (74 %) en 2005, 33 (66 %) en 2006 y 32 (64 %) en 2007-2008. Cinco especies resultaron registros nuevos para el sistema lagunar La Mancha: Brevoortia gunteri (2005), M. undulatus (2005), A. lamprotaenia (2006), Atherinomorus stipes (2007) y O. aenigmaticum (2007).
Familia | Especies | AB | Número (%) | Peso (%) | Frecuencia (%) | IVI (%) | CC | CT |
Dasyatidae | Hypanus sabinus | Hs | 0.05 | 5.93 | 0.48 | 6.46 | ME | 3° |
Elopidae | Elops saurus | Es | 0.07 | 0.50 | 0.64 | 1.20 | ME | 3° |
Megalopidae | Megalops atlanticus | Ma | 0.02 | 0.77 | 0.16 | 0.95 | ME | 3° |
Clupeidae | Brevoortia gunteri | Bg | 0.03 | 0.01 | 0.16 | 0.19 | ME | 1° |
Clupeidae | Harengula jaguana | Hj | 0.10 | 0.04 | 0.16 | 0.30 | ME | 1° |
Engraulidae | Cetengraulis edentulus | Ce | 0.30 | 0.03 | 0.32 | 0.65 | ME | 1° |
Engraulidae | Anchovia clupeoides | Ac | 0.13 | 0.02 | 0.48 | 0.63 | ME | 1° |
Engraulidae | Anchoa hepsetus | Ah | 0.15 | 0.03 | 0.32 | 0.49 | ME | 1° |
Engraulidae | Anchoa lamprotaenia | Ala | 0.28 | 0.06 | 0.32 | 0.66 | ME | 1° |
Engraulidae | Anchoa mitchilli | Am | 30.20 | 4.20 | 9.38 | 43.75 ♦ | ME | 1° |
Ariidae | Sciades felis | Sf | 0.08 | 2.30 | 0.64 | 2.98 | ME | 3° |
Ariidae | Cathorops aguadulce | Ca | 17.60 | 33.0 | 10.49 | 61.07 ♦ | E | 2° |
Batrachoididae | Opsanus beta | Ob | 0.05 | 0.21 | 0.48 | 0.74 | ME | 3° |
Mugilidae | Mugil curema | Mc | 2.20 | 9.10 | 4.93 | 16.17 | ME | 1° |
Atherinidae | Atherinomorus stipes | Ast | 0.02 | 0.01 | 0.16 | 0.18 | ME | 1° |
Hemirhamphidae | Hemirhamphus brasiliensis | Hb | 0.30 | 0.53 | 0.64 | 1.45 | ES | 1° |
Hemirhamphidae | Hyporhamphus meeki | Hm | 0.02 | 0.07 | 0.16 | 0.24 | ME | 1° |
Belonidae | Strongylura timucu | St | 0.03 | 0.08 | 0.32 | 0.43 | ES | 3° |
Belonidae | Strongylura notata | Sn | 2.80 | 9.20 | 2.23 | 14.20 | ME | 3° |
Synbranchidae | Ophisternon aenigmaticum | Oa | 0.02 | 0.02 | 0.16 | 0.19 | D | 3° |
Centropomidae | Centropomus undecimalis | Cun | 0.30 | 1.60 | 2.38 | 4.27 | ME | 3° |
Centropomidae | Centropomus pectinatus | Cpe | 0.08 | 0.40 | 0.48 | 0.96 | ME | 2° |
Centropomidae | Centropomus parallelus | Cp | 2.66 | 3.13 | 6.52 | 12.31 | ME | 3° |
Carangidae | Oligoplites saurus | Os | 0.20 | 0.04 | 1.43 | 1.66 | ME | 3° |
Carangidae | Trachinotus falcatus | Tf | 0.03 | 0.01 | 0.32 | 0.35 | ES | 3° |
Carangidae | Selene vomer | Sv | 0.08 | 0.27 | 0.79 | 1.14 | ES | 3° |
Carangidae | Hemicaranx amblyrhynchus | Ha | 1.41 | 0.70 | 2.23 | 4.30 | ES | 3° |
Lutjanidae | Lutjanus cyanopterus | Lc | 0.02 | 0.02 | 0.16 | 0.19 | ME | 3° |
Lutjanidae | Lutjanus griseus | Lg | 0.52 | 1.40 | 1.59 | 3.49 | ME | 3° |
Gerreidae | Gerres cinereus | Gc | 0.03 | 0.14 | 0.16 | 0.33 | ES | 1° |
Gerreidae | Eucinostomus melanopterus | Em | 3.98 | 0.64 | 3.97 | 8.60 | ME | 1° |
Gerreidae | Eucinostomus gula | Eg | 0.17 | 0.03 | 0.64 | 0.83 | ME | 1° |
Gerreidae | Eucinostomus argenteus | Ea | 1.54 | 0.02 | 0.79 | 2.36 | ME | 1° |
Gerreidae | Diapterus rhombeus | Dr | 18.50 | 6.22 | 10.17 | 34.88 ♦ | ME | 1° |
Gerreidae | Diapterus auratus | Da | 8.94 | 2.80 | 9.06 | 20.80 ♦ | ME | 1° |
Gerreidae | Eugerres plumieri | Ep | 1.24 | 4.31 | 4.29 | 9.84 | ME | 1° |
Haemulidae | Anisotremus surinamensis | As | 0.02 | 0.07 | 0.32 | 0.40 | ME | 2° |
Haemulidae | Pomadasys crocro | Pc | 0.02 | 0.02 | 0.16 | 0.19 | ME | 3° |
Sparidae | Archosargus probatocephalus | Ap | 0.08 | 1.32 | 0.32 | 1.72 | ME | 2° |
Sciaenidae | Micropogonias undulatus | Mu | 0.13 | 0.25 | 0.16 | 0.54 | ME | 3° |
Sciaenidae | Bairdiella chrysoura | Bc | 0.47 | 1.08 | 2.07 | 3.61 | ME | 2° |
Sciaenidae | Bairdiella ronchus | Br | 2.61 | 7.40 | 6.84 | 16.84 | ME | 2° |
Cichlidae | Mayaheros urophtalmus | Mur | 0.02 | 0.19 | 0.16 | 0.37 | D | 1° |
Eleotridae | Gobiomorus dormitor | Gd | 0.10 | 0.60 | 0.64 | 1.34 | E | 3° |
Gobiidae | Bathygobius soporator | Bs | 0.07 | 0.01 | 0.16 | 0.24 | ME | 3° |
Gobiidae | Gobionellus oceanicus | Go | 0.89 | 0.32 | 4.13 | 5.34 | ES | 2° |
Sphyraenidae | Sphyraena barracuda | Sb | 0.07 | 0.42 | 0.32 | 0.81 | ES | 3° |
Paralichthyidae | Citharichthys spilopterus | Cs | 1.21 | 0.62 | 4.93 | 6.76 | ME | 3° |
Paralichthyidae | Citharichthys uhleri | Cuh | 0.02 | 0.01 | 0.16 | 0.18 | ME | 3° |
Achiridae | Achirus lineatus | Ali | 0.30 | 0.04 | 2.07 | 2.41 | ME | 2° |
Abreviaturas (AB) utilizadas en la figura 8 (Análisis de Correspondencia Canónica). Índice de Valor de Importancia (IVI): ♦ especies dominantes. Componente comunitario (CC): ME = marino eurihalino, ES = marino estenohalino, E = estuarino, D = dulceacuícola. Categoría trófica (CT): 1° = primer orden, 2° = segundo orden, 3° = tercer orden. Abbreviations (AB) used in figure 8 are shown (Canonical Correspondence Analysis). Importance Value Index (IVI): ♦ dominant species. Community inhabitant (CC): ME = marine euryhaline, ES = marine estenohaline, E = estuarine, D = freshwater. Trophic category (CT): 1° = first order, 2° = second order, 3° = third order.
Parámetros ecológicos de la comunidad-escala espacial: Comparando los parámetros de diversidad (H’), riqueza (D) y equitatividad (J’) para el 2005, los valores promedio más altos de H´ (1.48 ± 0.22) se registraron en la estación 4, D (1.98 ± 0.20) en la estación 6 y J´(0.87 ± 0.09) en la estación 5, y los más bajos (H´= 0.43 ± 0.61; D = 0.56 ± 0.78 y J´= 0.31 ± 0.43) en la estación 7 (Figura 4A a C). Durante 2006 y 2007-2008, se observó un comportamiento similar de los índices, con los valores promedio más altos en la estación 3 (H´= 1.44 ± 0.3 y 1.38 ± 0.49; D = 1.84 ± 0.34 y 1.69 ± 0.67) y J´ (0.81 ± 0.65 y 0.78 ± 0.15) en la estación 4 (Figura 4D a I). El análisis de varianza mostró diferencias significativas entre sitios sólo para H´ en 2005 (Kruskal Wallis-H = 5.17, p < 0.0001).
Parámetros ecológicos de la comunidad-escala temporal: Entre los meses de muestreo, los valores promedio más altos de diversidad (H´= 1.52 ± 0.23) se presentaron en septiembre 2005 (Figura 5A), D (1.99 ± 0.46) en noviembre y J´ (0.87 ± 0.05) en marzo, ambos en 2007 (Figura 5H y 5I). Por otra parte, los promedios más bajos para los parámetros de diversidad y equitatividad (H´= 0.59 ± 0.67 y J´= 0.28 ± 0.29) se presentaron en febrero 2006 (Figura 5D y F) y para la riqueza (D= 0.77 ± 0.3) en julio 2007 (Figura 5H). Se encontraron diferencias significativas entre meses durante 2006 para H’ (ANOVA-F = 0.20, p < 0.05) y J’ (Kruskal Wallis-H= 14.21, p < 0.05); para 2007-2008 en diversidad (ANOVA-F = 3.0, p < 0.05) y riqueza de especies (ANOVA-F = 3.1, p < 0.05). El Cuadro 3 presenta los valores promedios entre años de estos parámetros, no se presentaron diferencias significativas (p > 0.05).
Parámetro | 2005 | 2006 | 2007-2008 | Estimador | p |
---|---|---|---|---|---|
Diversidad (H’) | 1.30±0.48 | 1.24±0.49 | 1.19±0.49 | F=0.32 | 0.73 |
Riqueza (D) | 1.61±0.59 | 1.51±0.62 | 1.43±0.65 | F=0.56 | 0.57 |
Equitatividad (J’) | 0.66±0.24 | 0.69±0.23 | 0.68±0.21 | F=0.09 | 0.91 |
Densidad (ind./m2) | 0.32±0.43 | 0.11±0.08 | 0.09±0.06 | H=5.63 | 0.06 |
Biomasa (g/m2) | 1.92±2.04 | 1.15±1.85 | 1.32±1.18 | F=1.54 | 0.22 |
Peso promedio (g/ind.) | 11.81±9.87 | 12.99±14.06 | 18.56±15.84 | F=2.11 | 0.13 |
Abundancia-escala espacial: Entre sitios de muestreo, los valores promedio más altos de densidad (0.62 ind./m2 ± 0.68) y biomasa (3.43 g/m2 ± 3.01) se presentaron en la estación 2 del 2005 (Figura 6A y B), y el peso promedio (29.6 g/ind. ± 22.10) en la estación 5 del 2007-2008 (Figura 6I). Los promedios más bajos de densidad (0.02 ind./m2 ± 0.04), biomasa (0.12 g/m2 ± 0.15) y peso promedio (5.10 g/ind. ± 1.13) se registraron en la estación 7 del 2005 (Figura 6A a C). Las especies con mayor densidad fueron A. mitchilli y D. rhombeus (estaciones 2 y 3), mientras que la mayor biomasa la aportó Hypanus sabinus, Strongylura notata, Mugil curema (estación 1 y 2), y Cathorops aguadulce (estación 6 y 7). El análisis de varianza para densidad (ANOVA-F = 0.003, p < 0.001), biomasa (Kruskal Wallis-H = 8.20, p < 0.05) y peso promedio (Kruskal Wallis-H = 3.09, p < 0.05) indicaron diferencias significativas entre sitios sólo en 2006.
Abundancia-escala temporal: Entre meses de muestreo se observó que los valores promedio más altos de densidad (1.17 ind./m2 ± 0.60) y biomasa (4.33g/m2 ± 2.58) se presentaron en abril del 2005 (Figura 7A y 7B), y el peso promedio (30.82 g/ind. ± 23.1) en julio del 2007 (Figura 7I). Los valores promedio más bajos de densidad (0.06 ind./m2 ± 0.05) se registraron en julio de 2007, de biomasa (0.68 g/m2 ± 0.30) en agosto de 2006 y peso promedio (3.85 g/ind. ± 1.13) en abril de 2005 (Figs. 7A a 7I). Durante abril (2005) las especies A. mitchilli y D. rhombeus fueron las más abundantes en número, mientras que C. aguadulce (junio 2006), H. sabinus (noviembre 2007) y S. notata (septiembre 2007) presentaron la mayor biomasa. Durante julio y enero se registraron las tallas más grandes de especies como B. ronchus (190-295 mm LT), C. aguadulce (140-272 mm LT), C. parallelus (120-255 mm LT), D. auratus (105-152 mm LT) y S. notata (192-450 mm LT) que corresponden a organismos preadultos y adultos, mientras las tallas jóvenes (< 150 mm LT) se registraron en abril-junio (secas). El análisis de varianza entre meses indicó diferencias significativas para densidad (Kruskal Wallis-H =12.35, p < 0.001) en 2005, para biomasa (Kruskal Wallis-H = 9.68, p < 0.05) en 2006, y para peso promedio (Kruskal Wallis-H = 12.83, p < 0.05) en 2007-2008. Los valores más altos de abundancia promedio anuales se muestran en el cuadro 3, el análisis de varianza no mostró diferencias significativas entre años (p > 0.05).
El índice del valor de importancia (IVI) indicó que cuatro especies fueron dominantes para el sistema: C. aguadulce, A. mitchilli, D. rhombeus y D. auratus, las cuales representaron el 75 % en número y el 46 % en peso de la captura total. La especie C. aguadulce presentó el valor más alto de IVI con los mayores porcentajes en peso (Cuadro 2). De acuerdo con la clasificación de componentes comunitarios, 37 (74 %) especies fueron marinas eurihalinas, 9 (18 %) fueron estenohalinas, 2 (4%) estuarinas y 2 (4 %) dulceacuícolas. Por su alimentación y hábitos alimentarios, 23 (46 %) especies son consumidoras de tercer orden, 19 (38 %) son de primer orden y 8 (16 %) de segundo orden. Se observó que especies carnívoras de tercer orden se distribuyeron principalmente en la zona norte-este del sistema (e.g. C. parallelus, H. amblyrhynchus), donde son abundantes los crustáceos en áreas de pastos marinos, mientras que la mayor abundancia de consumidores de primer orden (e.g. A. clupeoides, E. melanopterus) se registraron en la parte central y sur de la laguna.
Análisis peces-hábitat: El análisis de correspondencia canónica (ACC) mostró que la relación entre la abundancia de las especies, los sitios de muestreo y las variables ambientales fue explicada en un 57 % de la varianza total para los dos ejes de ordenación, con eigenvalores de 0.29 (eje 1) y 0.16 (eje 2). El ACC explicó el 36 % de la variación especies-ambiente para el eje 1, señalando que los factores que determinan en mayor medida la composición y distribución de las especies son la temperatura (r = 0.93), salinidad (r = 0.56) y oxígeno disuelto (r = 0.55) (Figura 8). Las especies que mostraron mayor correlación con el eje 1 fueron A. surinamensis, A. probatocephalus, Bathygobius soporator, C. edentulus, G. cinereus, H. brasiliensis, H. meeki, Lutjanus griseus, Pomadasys crocro, Selene vomer, S. notata, S. barracuda y T. falcatus de los componentes eurihalino y estenohalino (Cuadro 2). Estas especies presentaron abundancias de bajas a intermedias, asociadas con zonas someras (estaciones 1 y 2) de mayor salinidad y transparencia (r = -0.34). El eje 2 explicó el 21 % de la variación especies-ambiente, donde el oxígeno disuelto (r = 0.45), la salinidad (r = 0.40) y la profundidad (r= -0.83) son los factores más importantes en la distribución de las especies a lo largo de este eje. Las especies que mostraron mayor correlación con el eje 2 fueron A. hepsetus, A. clupeoides, A. stipes, B. ronchus, C. aguadulce, D. rhombeus, E. plumieri, E. saurus, Harengula jaguana, L. cyanopterus, M. undulatus, Citharichthys uhleri, M. atlanticus, M. urophtalmus, S. felis y O. aenigmaticum que pertenecen a los componentes eurihalino y dulceacuícola (Cuadro 2). Estas especies presentaron abundancias bajas y altas, asociadas a zonas de mayor profundidad con menor salinidad y transparencia (estaciones 3, 5, 6 y 7).
Discusión
En el sistema lagunar La Mancha, la variabilidad de los parámetros físicoquímicos se relacionó principalmente con la mayor descarga de los ríos en la porción sur, así como por la entrada de agua marina en el norte, ambos procesos determinados por los cambios estacionales en las distintas épocas climáticas de la región (Díaz-Ruiz et al., 2008). El sistema presentó un gradiente salino que varió desde condiciones marinas en el extremo norte de la laguna (estaciones 1 a 4), a condiciones mesohalinas con los menores valores en la parte sur (estaciones 5 a 7) debido al aporte del Río Caño Grande y Río Sábalo. La variación temporal de la salinidad dentro de la laguna estuvo relacionada con una alta precipitación y máxima descarga de los ríos durante los meses de la época de lluvias como ocurrió en julio del 2005 y 2007, respectivamente, donde se registraron las salinidades más bajas. Mientras que las salinidades más altas se registraron en abril 2006 y septiembre del 2007, como resultado de la poca o nula precipitación, y al incremento en la insolación y evaporación que ocurrió durante esos meses (Lara-Domínguez et al., 2006). Durante los meses de nortes la salinidad fue menor (febrero 2006) por el efecto de dilución causado por las lluvias de la temporada previa. La temperatura del agua mantuvo un comportamiento similar entre los sitios evaluados, y mostró poca variación desde la boca hacia el interior de la laguna, debido a que el sistema es somero y con buena mezcla vertical (Snedden et al., 2013). La temperatura anual se relacionó con las temporadas climáticas, con la mayor temperatura en los meses más cálidos (abril-mayo-julio) del periodo primavera-verano y la menor durante los meses más fríos (noviembre-diciembre) del periodo otoño-invierno. Los niveles más altos de oxígeno en la columna de agua se presentaron en la zona norte de la laguna, debido a los procesos de mezcla y corrientes que ocurren en las cercanías de la boca con el mar (Lara-Domínguez et al., 2006). El oxígeno disuelto también presentó un patrón relacionado a la estacionalidad climática, con mayor concentración durante los meses de noviembre y diciembre de 2005 y 2007, relacionado con las bajas temperaturas y acción del viento de la temporada de nortes y la menor durante la época de lluvias, atribuida al incremento de materia orgánica por la descarga de los ríos en la parte sur del sistema (De la Lanza Espino, Rodríguez, Baus, & Pulido, 2011; Aguirre-León et al., 2014). Respecto a los valores de pH, estos se comportaron dentro del intervalo normal de 6.5 a 8.5 los cuales están dentro de la variación encontrada en otras lagunas costeras similares (De la Lanza Espino et al., 2011). Los valores de pH fueron mayores en el periodo otoño-invierno (2005 y 2006), mientras que los más bajos en el periodo primavera-verano (2005-2007). Las variaciones de pH estuvieron asociadas al intercambio de agua de mar, a la lluvia, cantidad de escurrimientos, actividad biológica por organismos, procesos biogeoquímicos, degradación de materia orgánica o la remoción de sedimentos por corrientes causando remineralización (López-Ortega et al., 2010; Torres-Alvarado, Calva-Benítez, Álvarez-Hernández, & Trejo-Aguilar, 2016). La profundidad y transparencia mostraron un comportamiento inverso entre las localidades de muestreo, con la menor profundidad en la boca de conexión, debido a los procesos de azolve que existe en esta zona, mientras que la menor transparencia, se relacionó con las descargas de los ríos en la parte sur del sistema. Durante los meses de lluvias de 2005 y 2007, se registró la mayor profundidad debido al incremento del volumen de agua por los aportes fluviales y las precipitaciones en el sistema. Un comportamiento similar de estos parámetros fue observado en ecosistemas parecidos por Aguirre-León, Díaz-Ruiz y Ramírez-Huerta, (2010), Díaz-Ruiz et al. (2012), Sosa-López et al., (2012) quienes indicaron que gran parte de los procesos biológicos y ecológicos de las especies de peces, son afectados por las variaciones físicoquímicas del agua en estos sistemas.
En este trabajo se identificaron 50 especies de peces entre 2005 y 2008. En otros estudios realizados en este sistema, Reséndez-Medina y Kobelkowsky (1991) registraron 44 especies, mientras que Juárez, Rojas-Galaviz, Mora Pérez y Zárate-Lomelí (2006) reportaron 42 especies y Díaz-Ruiz et al. (2008) 38 especies, por lo tanto, el número total registrado más las identificadas en el presente estudio ascendió a 77 especies para este sistema. De acuerdo con los resultados, se presentaron especies exclusivas durante cada año, destacando C. edentulus (2005), A. lamprotaenia(2006), M. atlanticus (2007) y Lutjanus cyanopterus (2008). Al respecto, se considera que un esfuerzo adicional de pesca en esta laguna resultaría en la captura de otras especies en función de otros ambientes como las zonas cercanas a la desembocadura de los ríos y la boca de conexión con el mar (Díaz-Ruiz et al., 2008). Las especies que habitan en el sistema La Mancha están adaptadas a la variación de los factores ambientales a través del tiempo, principalmente al gradiente de salinidad, el cual controla los movimientos de migración de las especies dentro y fuera del sistema en las diferentes épocas climáticas del año. Esto se reflejó en el mayor porcentaje de especies (92 %) del componente eurihalino (e.g. A. mitchilli, B. ronchus, D. rhombeus) y estenohalino (e.g. Gobionellus oceanicus, H. amblyrhynchus, Sphyraena barracuda, S. vomer). Estos resultados sugieren que los peces han desarrollado estrategias y rutas de migración para el uso de los hábitats lagunares-estuarinos como áreas de crianza, alimentación, desove y protección en alguna etapa de sus ciclos de vida con relación a las condiciones ambientales en espacio y tiempo (Díaz-Ruiz, Cano-Quiroga, Aguirre-León, & Ortega-Bernal, 2004; Able, 2005; Cowan, Yáñez-Arancibia, Sánchez-Gil, & Deegan, 2013).
En lagunas costeras tropicales, se conoce la utilización de variados hábitats por las especies de peces, determinada por respuestas complejas a los factores ambientales y biológicos, lo cual se refleja en las variaciones estacionales de la estructura comunitaria (Able, 2005; Aguirre-León et al., 2014). En este estudio se observaron cambios en los valores de diversidad y riqueza de especies relacionados con las condiciones hidrológicas y los diferentes ambientes del sistema, donde estos parámetros fueron más altos en la zona norte de la laguna donde predomina un ambiente polihalino. Este comportamiento fue similar a lo reportado por Castillo-Rivera et al. (2002) en la Laguna de Pueblo Viejo (Veracruz), donde la mayor diversidad y riqueza se presentó en la boca de comunicación con el Golfo, y donde las variaciones de salinidad jugaron un papel preponderante en la composición de peces; adicionalmente, otros estudios han indicado que el oxígeno disuelto, profundidad y turbidez también influyen en la diversidad de peces (Deegan et al. 1986; Barletta et al., 2005; Greenwood, Matheson, McMichael, & MacDonald, 2007). Para el periodo 2005-2008 en La Mancha, los valores de diversidad y abundancia se relacionaron principalmente con las variaciones de salinidad, temperatura y profundidad durante cada época climática de la región. Autores como Thiel, Sepulveda, Kafemann y Nellen (1995), se refieren a estas variables como “controles” de los patrones temporales de la diversidad, distribución y composición de las especies en un sistema. En este sentido, se observó que la diversidad y riqueza fueron más altas al final de la temporada de lluvias e inicios de nortes (septiembre y noviembre), debido a la mayor abundancia de especies visitantes eurihalinas y estenohalinas del componente marino, que optimizan el uso del sistema, ligado a las variaciones estacionales de salinidad y al intercambio de especies con el mar durante las épocas climáticas de la región.
Las variaciones en los patrones de abundancia de las poblaciones también están influenciadas por los cambios en las condiciones ambientales del sistema y por los patrones de migración estacional de las especies durante sus ciclos de vida (Paperno, & Brodie, 2004; Miller, & Skilleter, 2006; Cowan, Grimes, & Shaw, 2008). Así, la mayor densidad y biomasa de las especies, se registró principalmente en la zona norte de la laguna durante los meses de abril y junio (secas), cuando la salinidad y la temperatura fueron elevadas, lo cual sugiere que especies estenohalinas y eurihalinas se desplazan del mar hacia la laguna como juveniles y adultos, para crecer y alimentarse durante esta época como A. mitchilli, C. edentulus, D. rhombeus, E. melanopterus, H. amblyrhynchus y S. notata. Mientras que en el mes de julio (lluvias) los valores promedio más altos de estos parámetros se presentaron en la zona sur del sistema, donde individuos preadultos y adultos del componente estuarino, eurihalino y dulceacuícola, utilizan el sistema para alimentarse, madurar y/o protegerse como Sciades felis, C. aguadulce, C. parallelus, D. auratus, M. curema y O. aenigmaticum. Variaciones similares en la abundancia de estas especies de peces, también han sido observadas en otras lagunas costeras del Golfo de México por Díaz-Ruiz, Pérez-Hernández y Aguirre-León (2003), Sosa-López, et al. (2012), Cowan et al. (2013), Aguirre-León et al. (2014).
En el sistema lagunar La Mancha existe una relación estrecha entre la abundancia de peces con los gradientes de salinidad, temperatura y profundidad. Esto se refleja en el evidente éxito ecológico de las especies dominantes, que, por su alta abundancia, frecuencia y distribución, han podido adaptar sus estrategias alimenticias, reproductivas, de crianza o protección, lo que ha favorecido su permanencia en el ecosistema (Lara-Domínguez et al., 2011; Aguirre-León et al., 2014). Tal es el caso de C. aguadulce, A. mitchilli, D. rhombeus y D. auratus que se consideraron como dominantes en el sistema lagunar en todos los años de muestreo. La especie C. aguadulce se distribuye como adulto preferentemente en la zona sur de la laguna, donde la salinidad y transparencia son menores por la descarga de los ríos Caño Grande y Sábalo, así, este bagre utiliza el sistema con fines tróficos y reproductivos durante lluvias. Cabe mencionar que la especie C. aguadulce ha sido citada anteriormente como C. melanopus (Marceniuk, & Betancur, 2008). De esta manera, los resultados son similares con las observaciones realizadas en C. melanopus (sensu lato C. aguadulce) como especie dominante en otros sistemas lagunares del Golfo de México, donde ha sido la especie más abundante y con amplia distribución en áreas de mayor turbidez, menor salinidad e influencia fluvial, para fines reproductivos en los meses de lluvias y nortes (Yáñez-Arancibia, & Lara-Domínguez, 1988; Ayala-Pérez, Ramos-Miranda, Flores-Hernández, Vega-Rodríguez, & Moreno-Medina, 2008; Aguirre-León et al., 2010). La especie pelágica A. mitchilli, por su comportamiento eurihalino, estuvo relacionada con una salinidad variable, y presentó amplia distribución en áreas cercanas a la boca de conexión. Esta anchoa se capturó con mayor frecuencia durante los meses de estiaje, y realiza migraciones cortas entre áreas someras y más profundas del sistema para crianza y alimentación. Cowan et al. (2013) mencionaron que esta especie es más abundante durante las temporadas de secas en Laguna de Términos. Sin embargo, Castillo-Rivera, Zárate-Hernández, Ortiz-Burgos y Zavala-Hurtado (2010) señalaron que la mayor abundancia de esta especie se registró en época de lluvias para el sistema lagunar Pueblo Viejo al norte de Veracruz. Las mojarras D. rhombeus y D. auratus son especies con dependencia estuarina y gran capacidad eurihalina; en el presente estudio, se distribuyeron principalmente en áreas someras de manglar y de mayor salinidad, cercanas a la boca de conexión de La Mancha; además, realiza migraciones entre los ambientes del sistema, para alimentación y protección. D. rhombeus fue más abundante durante los meses de la temporada seca, mientras que D. auratus en los meses de lluvias, lo que mostró un comportamiento similar al de otras lagunas costeras estudiadas en el Golfo de México por Aguirre-León y Yáñez-Arancibia (1986); Aguirre-León y Díaz-Ruiz (2006); Aguirre-León et al. (2010).
El análisis de las especies dominantes reflejó sus adaptaciones en el uso del sistema durante sus ciclos de vida, se observó una alternancia estacional de especies relacionadas a los cambios ambientales que ocurren en el sistema a través del tiempo. Estos cambios se deben tanto al régimen de lluvias y al aporte fluvial, como a la influencia marina de una época climática a otra, lo que modifica las variables físicoquímicas de la laguna. Al respecto, otros estudios realizados en sistemas costeros tropicales también han mencionado un comportamiento similar sobre la influencia de los factores ambientales, como la salinidad, temperatura y profundidad, entre otros, en la estructura de la comunidad de peces (Bouchereau, Chaves, & Monti, 2008; Lara-Domínguez et al., 2011; Díaz-Ruiz et al., 2012; Cowan et al., 2013).
La mayoría de los sitios de muestreo se ubicaron en áreas con vegetación sumergida de R. maritima y circundante de R. mangle. Estas áreas actúan como áreas de protección y alimentación de alevines y juveniles de especies marinas y estuarinas de importancia ecológica y económica (Yáñez-Arancibia, Lara-Domínguez, & Day, 1993; Bouchereau et al., 2008). Las zonas con vegetación tienen gran disponibilidad de alimento y mantienen poblaciones bénticas (Kanouse, La Peyre, & Nyman, 2006; Vaslet, Bouchon-Navaro, Charrier, Louis, & Bouchon, 2010), que a su vez son alimento para los peces y protegen a los juveniles de los depredadores. La mayoría de las especies de peces en estos sistemas tienen ciclos de vida típicamente estuarinos, con migraciones estacionales, que coinciden con los periodos de alta productividad y disponibilidad de alimento (Aguirre-León et al., 2014). Esto se reflejó en la composición trófica de los peces en La Mancha, ya que durante el periodo de estudio 2005-2008, predominaron especies consumidoras de primer y tercer orden en el sistema. Este comportamiento se relacionaría con la disponibilidad del alimento en el tiempo, las preferencias alimenticias de los peces o con la etapa del ciclo de vida de cada especie para utilizar el sistema (Rivera-Arriaga, Lara-Domínguez, Villalobos, & Yáñez-Arancibia, 2003; Pineda-Peralta, Ramos-Miranda, Castillo-Rivera, & Arreguín-Sánchez, 2016). Estudios realizados en otras lagunas del Golfo de México, reportan mayor abundancia de peces consumidores de segundo orden (Yáñez-Arancibia, Lara-Domínguez, Aguirre-León, & Díaz-Ruiz, 1986; Christensen & Pauly, 1993). Este comportamiento difiere de lo encontrado para la laguna La Mancha, sin embargo, es necesario realizar mayor investigación sobre los hábitos alimenticios y relaciones tróficas de las especies de este sistema, ya que actualmente es muy escasa.
En la relación peces-hábitat se mostró que la salinidad, temperatura, oxígeno disuelto y profundidad, fueron las variables más importantes en la variación de la composición y distribución espacio temporal de las especies de peces en La Mancha. Dicho análisis permitió identificar la presencia de conjuntos de peces relacionados con los cambios físicoquímicos, donde la abundancia de las especies estuvo condicionada por las migraciones estacionales y por especies residentes del sistema. De esta manera, las especies asociadas con mayor influencia marina cercana a la boca de conexión fueron componentes estenohalinos y eurihalinos estacionales tales como A. surinamensis, S. notata, H. meeki, T. falcatus, H. brasiliensis, E. argenteus, H. amblyrhynchus, C. edentulus y M. curema, distribuyéndose en áreas someras con vegetación sumergida. El conjunto de especies asociadas a menores salinidades estuvo caracterizado por aquellas que toleran amplias variaciones de este parámetro, como C. aguadulce, A. mitchilli, D. auratus, D. rhombeus, E. plumieri y B. ronchus, que utilizan el sistema principalmente para alimentación, maduración y protección. Este comportamiento de grupos de peces también ha sido analizado por Díaz-Ruiz et al. (2004), Greenwood et al. (2007), Bouchereau et al. (2008), Arceo-Carranza et al. (2010), Vaslet, Bouchon-Navaro, Charrier, Louis y Bouchon (2010).
El análisis de la información presentada permite establecer un evidente gradiente ambiental determinado por los procesos hidrológicos que ocurren en el sistema, y que están relacionados con los cambios estacionales de la región. De la misma manera, los valores de diversidad, abundancia y frecuencia de la ictiofauna responden a ese gradiente en espacio y tiempo. Esta información, también nos permite entender los procesos ecológicos que ocurren en ecosistemas similares e identificar hábitats prioritarios de conservación y hacer uso adecuado de los recursos naturales. Este estudio constituye un punto de partida para evaluaciones posteriores, y una herramienta de monitoreo del estado de salud ambiental de esta laguna, a fin de contribuir a la conservación de su biodiversidad y protección como área importante para la crianza, alimentación, protección y/o reproducción de recursos pesqueros en este sitio con categoría Ramsar.