Muños-Rojas Et Al., 2014 PDF
Muños-Rojas Et Al., 2014 PDF
Muños-Rojas Et Al., 2014 PDF
0002
Resumen
Uno de los componentes estructurales más abundante en la RBLP son los peces, los cuales
desempeñan diversos papeles ecológicos entre ellos el transporte de energía, por lo que el objetivo
de este trabajo fue realizar un análisis comparativo de la composición y abundancia de los peces
se recolectaron muestras biológicas con una frecuencia de dos horas en cada ciclo nictemeral.
Las pruebas de comparación de las variables ambientales mostraron que no existen diferencias
significativas entre horas de día y horas de noche, a excepción del pH. Las comparaciones entre
y horas de noche. Se identificaron dos conglomerados entre las especies dominantes, así como
una independencia de las variables ambientales. Los peces en la RBLP, aprovechan las diversas
condiciones de hábitat para desarrollar todo o parte de su ciclo de vida, y son responsables de
1
Licenciatura en Biología. Universidad Autónoma Metropolitana Unidad Xochimilco. Calzada del Hueso 1100. Col. Villaquietud,
Coyoacan, 04960, Cd. de México.
2
Departamento El Hombre y su Ambiente. Universidad Autónoma Metropolitana Unidad Xochimilco. * luayala@correo.xoc.uam.mx
3
Instituto de Ecología, Pesquerías y Oceanografía del Golfo de México. Universidad Autónoma de Campeche. Av. Héroe de Nacozari
480. C.P. 24029. San Francisco de Campeche, Campeche. México
Abundancia nictemeral en Los Petenes
Revista Digital de Investigación Científica
Variabilidad ambiental.
Abstract
Los Petenes Biosphere Reserve (RBLP, for its acronym in Spanish) integrates inland and marine
relevance. One of the most abundant structural components in the RBLP are fishes that play
some ecological roles, such as energy transportation. Thus, the objective of this study was to
analyze the composition and abundance of fishes comparatively in nycthemeral cycles depending
on environmental variability. Four sampling campaigns were performed in February and May
both in 2013 and 2016, when environmental variables were measured and biological samples
were collected in each nycthemeral cycle in a two-hour frequency. The comparison tests for
environmental variables showed no significant differences between day and night hours, except
for pH. Comparisons between nycthemeral cycles showed some differences considering dissolved
oxygen and pH. A total of 38 species were identified grouped in nine orders and 20 families, 11
of which were dominant. Comparison for density and biomass showed no significant differences
between day and night. Two associations of dominant species were observed, as well as a high
degree of independence for environmental variables. The fish in the RBLP take advantage of the
habitat conditions in developing all or part of their life cycle, and they are responsible for energy
abundance.
Introducción
La reserva de la biósfera Los Petenes (RBLP) en el estado de Campeche, México, integra diversos
ambientes, tanto continentales, como marinos, conformando una unidad biogeográfica única
derivada de su alta diversidad de flora y fauna, así como de ecosistemas, brindándole una gran
importancia biológica, ecológica, económica y científica (Ayala-Pérez et al. 2014; Corbalá et al.
2007). En su porción marina, la reserva cuenta con la mayor superficie y en mejor estado de
conservación de pastos marinos a nivel nacional, dominados por Thalassia testudinum, que son
hábitat propicio para alimentación, protección y crianza de diversas especies, incluidos recursos
pesqueros con manejo comercial, como el pulpo Octopus maya y algunos peces de escamas (lisas,
Los peces costeros utilizan los diversos hábitats de una manera programada en tiempo y
espacio, y una forma de estimar dicho uso es mediante el análisis de los cambios en la composición
ambiental, tal y como lo han planteado Arceo-Carranza et al. (2012), Arceo-Carranza y Vega-
Cendejas (2009) y Lugendo et al. (2007). Para el caso particular de la RBLP donde la diversidad de
hábitats está determinada por una compleja interacción entre manglar, pastos marinos y petenes,
el conocimiento es escaso y algunas aproximaciones las han planteado Ayala-Pérez et al. (2014),
Dadas las diversas actividades económicas que el hombre desarrolla en la RBLP, es necesario
continuar con la generación de información sobre las interacciones que los peces establecen entre
ellos y con sus hábitats y que les permiten desarrollar todo o parte de sus ciclos de vida y de esta
manera, entender la relevancia de conservar tanto sus hábitats como su diversidad y abundancia
consientes de los múltiples bienes y servicios que nos proporcionan. Es por ello que, el objetivo
de este estudio fue analizar los cambios en la composición y abundancia de los peces de manera
Materiales y métodos
La RBLP se localiza en la costa norte del estado de Campeche, en el sureste del Golfo de México,
limita al norte con la Reserva de la Biósfera Ría Celestún, al este colinda con los municipios de
Tenabo, Hecelchakán y Calkiní, al sur con la ciudad de San Francisco de Campeche y su límite
oeste se encuentra doce millas náuticas mar adentro que corresponde al mar territorial (CONANP
2006). El sitio de estudio se ubicó en la localidad conocida como Balantaunche cerca de la isla
| 35
Abundancia nictemeral en Los Petenes
Revista Digital de Investigación Científica
un canal que comunica con un peten y la amplia distribución de parches de pastos marinos con
manglar de borde (Ayala-Pérez et al. 2014).
De acuerdo con Torres-Castro et al. (2008) y Villalobos-Zapata (2004) en la región se identifican
dos temporadas climáticas: la temporada de secas que abarca de noviembre a mayo y la de
lluvias de junio a octubre. El clima predominante en la zona centro-sur de la RBLP es Aw (cálido
subhúmedo con lluvias en verano), mientras que en su extremo norte es de tipo BS’h’w (semiseco
y seco-cálido).
Se realizaron cuatro campañas de muestreo en febrero y mayo de 2013 y 2016. En cada campaña
se registraron las mediciones de variables ambientales: temperatura (ºC), salinidad (UPS), oxígeno
disuelto (mg/L) y pH del agua en dos niveles de profundidad (superficie y fondo), empleando
una sonda multiparamétrica ISY 556. También se recolectaron muestras biológicas del necton con
una red de arrastre de cinco metros de largo, 2.5 m de abertura de trabajo y 19 mm de luz de malla,
operada por doce minutos a bordo de una lancha con motor fuera de borda, a una velocidad
promedio de 2 nudos para cubrir un área de arrastre aproximada de 1,800 m2. La frecuencia de los
muestreos fue de dos horas durante un periodo de 24 horas.
El ciclo nictemeral se dividió en horas de día y horas de noche, considerando las 7 horas como
la primera hora de día y las 19 horas como la primera hora de noche. Los organismos capturados
se fijaron con formaldehido al 15% y se almacenaron en bolsas de plástico apropiadamente
etiquetadas que fueron trasladadas al laboratorio para su procesamiento.
Todos los organismos capturados fueron identificados taxonómicamente con ayuda de literatura
especializada (Castro-Aguirre et al. 1999; Cervigón et al. 1992; Fischer 1978) y se ordenaron
sistemáticamente de acuerdo con lo propuesto por Nelson (2016). Así mismo, cada individuo fue
medido en longitud total (mm), longitud estándar (mm) y peso total (g) utilizando un ictiómetro
convencional y una balanza digital con capacidad de 6.0 Kg y 0.1 g de precisión.
El comportamiento nictemeral de las variables ambientales se analizó mediante la representación
gráfica de líneas de tendencia y se realizaron pruebas estadísticas de normalidad (W de Shapiro-
Wilk) y de comparaciones paramétricas (ANOVA) y no paramétricas (Kolmogorov-Smirnov), así
como la prueba a posteriori de Tukey. Estos procedimientos se realizaron con ayuda del software
Past3 (Hammer et al. 2001).
Para determinar el dominio ecológico de las especies de peces se utilizó el índice de importancia
relativa (IIR) modificado por Koranteng (2001), que se define como:
IIR = (%N)*(%W)*(%F)
36 | Áreas Naturales Protegidas Scripta. 2017
Ciprés-Chávez, Ayala-Pérez, Rojas-Rubio y Sosa López
Revista Digital de Investigación Científica
total; %W= Porcentaje de contribución en peso de la especie i a la captura total; %F= Frecuencia
de aparición porcentual de la especie i con relación al número total de muestreos. De acuerdo con
el autor, las especies con valores IIR ≥ 20 se consideran especies dominantes, valores entre 20 y 1
Por otro lado, para el análisis de la abundancia se estimaron los valores de densidad (ind/m2),
biomasa (g/m2) y peso promedio (g/ind) y se utilizaron los valores promedio por horas de día
y de noche para realizar las comparaciones. Para la determinación de las asociaciones entre las
especies dominantes, se realizó un análisis clúster (método Ward, distancia Gamma), utilizando
las matrices de abundancia numérica por especie por ciclo de muestreo, con ayuda del software
estadístico Systat 13.0 (Systat Software, Inc. 2009). Adicionalmente, se utilizaron las matrices de
número de individuos por especie y hora y la de variables ambientales por hora para desarrollar un
Análisis de Correspondencia Canónica (ACC), con ayuda del programa estadístico multivariado
MVSP (Kovach Computing Services, 2016) con la intención de representar gráficamente las
de las especies utilizando el código alfa que integra las dos primeras letras del género y dos de
la especie. Finalmente se analizan y discuten las diferencias observadas entre día y noche y entre
meses, considerando la variabilidad ambiental, así como las características biológicas y ecológicas
de las especies.
Resultados
representan en la Figura 1. La temperatura promedio del agua para el periodo de muestreo fue de
26.7°C con un intervalo de variación de 20.3 a 32°C. En escala temporal la temperatura promedio
fue mayor en mayo de 2016 y la menor en febrero de 2016. En escala nictemeral se observa que
la mayor temperatura se presentó a las 17 horas y la mínima a las 11 horas. Para el caso de la
en mayo 2013 se registraron los valores más altos, mientras que en febrero 2016 los menores y
| 37
ABUNDANCIA NICTEMERAL EN LOS PETENES
Revista Digital de Investigación Científica
en general las variaciones son mínimas en los ciclos nictemerales. La concentración de oxígeno
disuelto mostró un promedio de 2.7 mg/L con un intervalo de variación de 1.2 a 6.3 mg/L, donde
los valores promedio más elevados se registraron en mayo de 2013 y los mínimos en febrero de
2016. En escala nictemeral no se aprecia un comportamiento distintivo entre horas de día y horas
de noche. Por último, el pH osciló con valores entre 7.3 y 8.9 con un promedio de 8. En febrero y
mayo de 2016 se observan los valores más bajos mientras que los mayores se observaron en 2013.
Una ligera tendencia de incremento del pH se distingue para las 17-19 horas.
Figura 1. Variación nictemeral de las variables ambientales en la porción central de la RBLP, Campeche, México.
por día y noche presentaron una distribución normal (p<0.05), por lo cual se decidió aplicar
una prueba paramétrica ANOVA de una vía para efectuar las comparaciones entre los ciclos
entre horas de día y horas de noche (p<0.05). Para la comparación entre las campañas de muestreo,
entre campañas (p<0.05). Sin embargo, en cuanto a oxígeno disuelto y pH los meses febrero de
Se capturaron un total de 9 726 individuos con un peso conjunto de 156 kg. Se clasificaron
sistemático de las especies. Para el conjunto de muestras recolectadas se identificaron once especies
con dominio ecológico que representan el 96% en número de individuos y 94% en peso del total
de la captura (Tabla 2). La especie más importante por su abundancia fue Lagodon rhomboides,
más altos durante la noche particularmente en febrero de 2013, así como un valor destacado a las
19 horas en febrero de 2016. En cuanto a peso promedio, no se aprecian diferencias entre horas
de día y horas de noche, sin embargo, durante la noche se aprecia una ligera diferencia entre las
campañas de mayo y febrero de los dos años de estudio. Es notorio que entre las horas 13 a 17 en
La matríz de número de individuos de las especies dominantes agrupadas por horas de día y
Wilk cuyos resultados fueron p= 0.0033 y p= 0.0017, para horas de día y p= 0.0022 y p= 0.0007 para
El análisis clúster de la matriz de abundancia numérica por especie agrupada por mes y hora,
permitió generar un dendrograma donde se distingue la conformación de dos grupos (Figura 3).
De acuerdo con el criterio de varianza mínima de Ward en cada paso el par de clúster con
distancia mínima entre ellos son mezclados y de esta forma se representa el dendrograma, por
otra parte, la distancia Gamma permite analizar en una escala ordinal los datos de abundancia,
con lo cual se distingue la mayor disimilaridad de S. testudineus y C. penna, que son especies
| 39
Abundancia nictemeral en Los Petenes
Revista Digital de Investigación Científica
Tabla 1. Listado sistemático de acuerdo con Nelson (2016) de las especies de peces recolectados en la
porción central de la RBLP. Se muestran los códigos
alfa de las especies dominantes.
del 31.6% y 9.08% y una correlación especie-ambiente de 0.83 y 0.716 para los ejes uno y dos
respectivamente. En el cuadrante uno se presentan las especies C. penna, L. analis y A. rhomboidalis
con una aparente independencia de las variables ambientales. En el cuadrante dos se observa
la presencia aislada de E. gula. En el cuadrante tres se distingue el dominio de H. aurolineatum y
L. rhomboides y finalmente en el cuadrante cuatro se agrupan L. synagris, O. beta, H. plumierii, S.
testudineus y O. chrysoptera con una marcada asociación con las variables de temperatura, salinidad
y oxígeno disuelto (Figura 4).
Discusión
La variabilidad ambiental en la región de “Los Petenes”, está condicionada principalmente por la
temporalidad climática, procesos oceanográficos como el sistema de mareas, procesos climático-
meteorológicos, por ejemplo, la precipitación, y por una condición geográfica como la extensión
de la plataforma continental con condiciones de baja profundidad, lo que propicia una amplia
cobertura de pastos marinos (Ayala-Pérez et al. 2014; Castillo-Rivera et al. 2005; Lara-Lara et al.
2008).
excepción del pH, de esta forma existe evidencia de estabilidad ambiental entre horas de día y
horas de noche. Para el caso de la temperatura los intervalos de variación oscilaron entre 1.41
| 41
Abundancia nictemeral en Los Petenes
Revista Digital de Investigación Científica
Figura 2. Comportamiento de la abundancia de la comunidad de peces (densidad, biomasa y peso promedio) en los
cuatro ciclos nictemerales realizados en la RBLP.
y 2.08ºC, lo cual es un claro reflejo de las condiciones tropicales de la región. En este sentido,
meteorológicas aledañas a la RBLP, las cuales son similares a lo que se reporta en el presente
la temperatura entre los meses de mayo y febrero, pero además en febrero de 2016 los registros
fueron anormalmente bajos, debido a que los muestreos se realizaron bajo los efectos de un frente
Figura 3. Dendrograma del análisis clúster de la abundancia de especies dominantes en ciclos nictemerales
en la RBLP. Los códigos alfa de las especies son descritos en la Tabla 2.
atribuyéndose a que los muestreos fueron en la misma temporada climática de secas. Sin embargo,
entre campañas de muestreo existieron diferencias estadísticas, destacando los valores más bajos
obtenidos en febrero de 2016, observando nuevamente los efectos del frente frío que trajo lluvia.
de cenotes, petenes, ojos de agua y ríos subterráneos mantienen el balance hídrico de la región,
| 43
Abundancia nictemeral en Los Petenes
Revista Digital de Investigación Científica
De igual forma, el oxígeno disuelto no presentó fluctuaciones importantes entre horas día y horas
noche. Las comparaciones entre algunas campañas de muestreo, mostraron diferencias estadísticas,
ya que el comportamiento del oxígeno disuelto se puede interpretar como una respuesta al efecto
movimiento del agua, los vientos y la temperatura, que promueven el intercambio gaseoso con
la materia orgánica, provocando esta variación entre campañas pero no así a escala nictemeral
Figura 4. Análisis de Correspondencia Canónica (ACC) entre la abundancia de las especies dominantes con las variables ambientales
en la porción central de la RBLP.
con lo reportado por Ayala-Pérez et al. (2014) para un ciclo nictemeral. Dicho comportamiento se
debe al efecto del lavado de sedimentos carbonatados principalmente por el agua continental que
recorre la zona costera de origen cárstico propia de la península de Yucatán (Muñoz-Rojas et al.
2013). Así mismo, los valores promedio de pH para praderas de pastos marinos dominados por
T. testudinum reportados por MacDonald (2011) y Torres-Castro et al. (2008), concuerdan con los
de este trabajo.
Los sistemas estuarinos como es el caso de la RBLP, se caracterizan por tener una amplia
variedad de hábitats, lo que se refleja en una compleja estructura ecológica con alta productividad
primaria como los pastos marinos y el manglar, la cual es aprovechada para el desarrollo de
todo o parte de los ciclos de vida de los peces (Lugendo et al. 2007). Villalobos-Zapata (2004)
reporta 68 especies de peces para toda la reserva, mientras que en este estudio se identificaron
38 especies sólo en un sitio, lo cual destaca la variedad de hábitat. En diversos reportes se señala
que las familias Sparidae, Haemulidae y Lutjanidae son las más importantes por su diversidad y
Entre las especies dominantes destacan por su abundancia L. rhomboides y H. plumierii, que
muestran un amplio intervalo de tolerancia a las variables ambientales, lo que concuerda con lo
et al. (2012).
son comparativamente diferentes en promedio ya que son cuatro veces superiores en horas de
noche respecto de las horas de día. Autores como Castillo-Rivera et al. (2005) y Lugendo et al.
(2007) han descrito las preferencias de muchas especies de peces por una actividad nocturna tanto
para la búsqueda de alimento como para la evasión de sus depredadores. En la RBLP el hábitat
de pastos marinos ofrece las condiciones óptimas para el desarrollo de funciones como la crianza
y, por tanto, se convierte en una zona de alimento para peces de mayor tamaño. En cuanto al
peso promedio no hay diferencias generales entre horas de día y horas de noche, sin embargo, se
distingue una tendencia particular de incremento entre las 17 y 21 hrs, lo que se interpreta como el
periodo de penumbra donde los peces de mayor talla demuestran su preferencia de horario para
alimentación.
y C. penna, que en general tienen hábitos alimenticios basados en el zoobentos y una importancia
pesquera local ya que alcanzarían tallas superiores a los 30 cm (Fishbase 2017). El segundo grupo
es más diverso, con especies que tienen hábitos alimenticios contrastantes que demuestran una
| 45
Abundancia nictemeral en Los Petenes
Revista Digital de Investigación Científica
disponibilidad de alimento. Tanto la cadena del detritus como la cadena del pastoreo, se observa
en la comunidad de peces de la RBLP. Para la cadena del detritus la aportación de materia orgánica
tanto de los manglares como de los pastos marinos en la porción costera litoral es abundante
y por el otro lado, las condiciones que favorecen la producción fitoplanctónica como base de
la cadena del pastoreo son buenas, lo que asegura un soporte abundante de alimento para los
consumidores de segundo orden. Estas consideraciones tienen como base las observaciones de
campo acompañadas de los reportes sobre las cadenas tróficas basadas en el detritus y en el
García et al. (2009), Villalobos-Zapata y Mendoza-Vega (2010) y Day et al. (2012), lo cual también
sirve para argumentar que los peces son responsables de flujos energéticos que contribuyen a la
muestran las especies dominantes en dos grupos y por otro lado la alta dispersión de las variables
se justifica por las inusuales condiciones que se presentaron en febrero de 2016 cuando el muestreo
se desarrolló bajo los efectos de un frente frío estacionario. Adicionalmente, es claro que las
del hábitat, por lo cual la dependencia de alguna variable ambiental es mínima al menos para el
ciclo día-noche.
Sin duda, los resultados de este estudio están limitados a dos meses del año que si bien pueden
ambiental en un ciclo anual. Autores como Hoeinghaus et al. (2003), Kobza et al. (2004) y Allen et al.
(2006) han descrito el efecto que las variables ambientales tienen sobre la composición, abundancia
de ecosistema con gran homogeneidad ambiental, alimento abundante y diverso, así como
diversidad de hábitats que la comunidad de peces aprovecha y mantiene a través de los flujos
energéticos que desarrolla para cubrir las diferentes etapas de sus ciclos de vida.
Conclusiones
demuestra una mayor actividad, reflejada por su abundancia, durante las horas de noche. Sin
que en conjunto representan más del 90% de la captura total. Los peces desarrollan todo o parte de su
ciclo de vida acoplados a la variabilidad ambiental, desarrollando nichos que les permite reducir la fricción
por competencia, pero además contribuir significativamente a la resiliencia del ecosistema, por lo tanto,
es necesario desarrollar más estudios considerando las escalas nictemeral, espacial y temporal para otras
Literatura citada
Allen, T., M. Jiménez y S. Villafranca. 2006. Estructura y categorías tróficas de peces asociados a praderas
Arceo-Carranza, D., E. Vega-Cendejas y M. Hernández de Santillana. 2012. Day and night trophic
variations of dominant fish species in a lagoon influenced by freshwater seep. Journal of Fish Biology
82: 54-68.
costa de Campeche, México. Latin American Journal of Aquatic Research 40 (1): 63-78.
| 47
Abundancia nictemeral en Los Petenes
Revista Digital de Investigación Científica
Fischer, W. 1978. FAO Species identification sheets for fishery purposes. Western Central Atlantic.
Hammer O., A.T.D. Harper y P.D. Ryan. 2001. PAST: Paleontological statistics software package for
education and data analysis. Palaeontología Electrónica 4(1): 4-9.
Hoeinghaus, D., C. Layman, D. Albrey y K. Winemiller. 2003. Spatiotemporal variation in fish
assemblage structure in tropical floodplain creeks. Environmental Biology of Fishes 67: 379-387.
Kobza, R.M., J.C. Trexler, W.F. Loftus y S.A. Perry. 2004. Community structure of fishes inhabiting
aquatic refuges in a threatened karst wetland and its implications for ecosystem management.
Biological Conservation 116: 153-165.
Koranteng, K. 2001. Structure and dynamics of a demersal assemblage on the continental shelf an upper
slope off Ghana, West Africa. Marine Environmental Research 220: 1-12.
Kovach Computing Services. 2013. Multi Variate Statistical Package 3.22. Anglesey, Wales.
Consultado: 10-01-2017. En: https://www.kovcomp.co.uk/mvsp/mvspwbro.html
Landínez-García, R.M., S.P. Ospina-Guerrero, D.J. Rodríguez-Castro, R. Arango y E. Márquez.
2009. Análisis genético de Lutjanus synagris en poblaciones del Caribe Colombiano. Ciencias
Marinas 4 (35): 321-331.
Lara-Domínguez, A. L. y A. Yáñez-Arancibia. 1999. Productividad secundaria, utilización de hábitat
y estructura trófica. pp. 153-166. En: Yáñez-Arancibia A. y A.L. Lara-Domínguez (Eds).
Ecosistemas de manglar en América Tropical. Inst. Ecol. AC. México, UICN/HORMA Costa
Rica, NOAA/NMFS Silver Spring MD USA. 380 pp.
Lara-Lara, J.R., V. Arenas, C. Bazán, V. Díaz, E. Escobar, M. García, G. Gaxiola, G. Robles, R. Sosa,
L.A. Soto, M. Tapia y J.E. Valdez. 2008. Los ecosistemas marinos. pp. 135-159. En: Soberón, J.,
G. Halffter y J. Llorente-Bousquets (eds.). Capital natural de México, vol. I: Conocimiento
actual de la biodiversidad. CONABIO, México. 621 pp.
Lugendo, B.R., A. De Groene, I. Cornelissen, A. Pronker, I. Nagelkerken, G. Van-der-Velde and Y.D.
Mgaya. 2007. Spatial and temporal variation in fish community structure of a marine embayment in
Zanzibar, Tanzania. Hidrobiología 58:1-16.
MacDonald, B. 2011. Análisis de la diversidad de las praderas de pastos marinos en la laguna Marino-
costera de la Bahía La Graciosa, Izabal, Guatemala. Tesis de Licenciatura. Universidad de San
Carlos de Guatemala. Guatemala. 93 pp.
Muñoz-Rojas, S., L.A. Ayala-Pérez, A. Sosa-López y G.J. Villalobos-Zapata. 2013. Distribución y
abundancia de la comunidad de peces en la porción litoral de la Reserva de la Biósfera Los Petenes,
| 49
Abundancia nictemeral en Los Petenes
Revista Digital de Investigación Científica
Nelson, J.S., 2016. Fishes of the World. John Wiley and Sons 5th Ed. Nueva York, E.U. 752 pp.
Ramos, O.R., R. Sepúlveda y F. Villalobos. 2003. El agua en el medio ambiente, muestreo y análisis.
Systat Software, Inc. 2009. Systat 10.2. Consultado: 21-01-2017 En: https://systatsoftware.com/
contact/get-a-quote-on-systat-software-products/
Hernández. 2008. Ictiofauna de sistemas cárstico palustres con impacto antrópico: Los Petenes de
Villalobos-Zapata G.J. 2004. Reservas de la Biósfera: Los Petenes y Ría Celestún. 397-412 pp. En:
Cita
https://doi.org/10.18242/anpscripta.2017.03.03.02.0002