Comunidades de fondos blandos submareales de la bahía de Navidad (Jalisco, México): patrones espaciales y temporales

Soft-sediment benthic assemblages of Bahía de Navidad (Jalisco, México): spatial and seasonal patterns

Antonio Corgos^*, Enrique Godínez-Domínguez, Juan Ramón Flores-Ortega y Jorge A. Rojo-Vázquez

Departamento de Estudios para el Desarrollo Sustentable de Zonas Costeras, Universidad de Guadalajara. Av. Valentín Gómez Farías 82, 48980 San Patricio Melaque, Jalisco, México. * acorgos@gmail.com

Recibido: 27 agosto 2013
Aceptado: 28 agosto 2014

Resumen

Se analizaron los patrones de variación espacial y temporal de las comunidades bentónicas submareales de una bahía arenosa semiexpuesta en relación con las condiciones ambientales. Se seleccionaron 4 estaciones de muestreo con diferentes condiciones como tamaño de sedimento, exposición al oleaje y proximidad de descargas de agua continental, y se realizó un muestreo mensual utilizando una bomba de succión durante un ciclo anual. Se tomaron muestras de sedimento para su análisis granulométrico y su contenido en materia orgánica. Para los análisis estadísticos los datos fueron agrupados en 2 periodos hidroclimáticos: verano e invierno. La parte norte de la bahía presentó sedimento más fino como limos y arcillas, mayor porcentaje de materia orgánica, es decir, más del doble que la parte sur durante el verano, mayores densidades de individuos, con dominancia de escafópodos y poliquetos (88% de la abundancia), así como bajos niveles de diversidad y equitatividad, sobre todo, en el verano. La parte sur presentó arena fina con menor contenido en materia orgánica, menores densidades y mucha mayor diversidad y equitatividad sin grandes diferencias entre invierno y verano. Se observó una correlación muy elevada entre la composición de las comunidades bentónicas y el tamaño de sedimento y la cantidad de materia orgánica (Spearman r_s= 0.94).

Palabras clave: comunidades bentónicas, fondos blandos, zona costera, tropical, océano Pacífico, materia orgánica, tamaño sedimento.

Abstract

Spatial and seasonal patterns of benthic assemblages in a semiexposed sandy bay were analyzed in relation to sedimentary variables and organic load. Four sampling stations with different environmental conditions; sediment size, wave exposure and proxmity of freshwater runoff; were selected along the bay. Monthly samplings were carried out during a year using a suction sampler. Sediment samples were also taken to analyze granulometry and organic load. For analysis purposes data were pooled according to climatic periods present in the area, dry season-winter and rainy season-summer. The north part of the bay showed fine sediment, silt and clay, with higher organic matter percentage, doubling the south percentages during summer, higher density of macrofauna, with a clear dominance of scaphopods and polychaetes; 88% of total abundance; and lower diversity and evenness. Scaphopods showed densities up to 9 788 ind·m^-2. The south part of the bay showed coarser sediment; mainly fine sand; with lower organic matter content, lower macrofaunal density and higher diversity and evenness. Scaphopod densities were much lower in this area. High significant correlation was found between macrofaunal assemblages and the percentage of silt and organic matter (Spearman r_s= 0.94).

Introducción

Las áreas costeras son consideradas hábitats de cría para numerosas especies de organismos marinos. Tradicionalmente, los hábitats de cría han sido definidos como áreas donde los ejemplares juveniles se encuentran en elevadas densidades, evitan la depredación con más éxito o crecen más rápido (Gillanders et al., 2003). El paradigma de Gunter (1967) establece que los ambientes estuarinos representan a las áreas de cría por definición debido a la gran dependencia de este hábitat por parte de numerosas especies de peces en el sureste de Estados Unidos de América. Sin embargo, con el conocimiento de las historias de vida de peces de otras regiones, se encontró que donde la distribución de ambientes estuarinos es muy reducida, las zonas costeras someras y protegidas constituían verdaderas áreas de cría, y los macroinvertebrados bentónicos, el grupo trófico más crítico, al que se denominó "especies forraje" (Blaber, 1997) y que son la fuente de alimento de numerosas especies de peces. Las zonas someras de la costa de Jalisco fungen como zonas de alimentación y cría para muchas especies de peces de interés comercial, siendo la macrofauna bentónica de gran importancia en su dieta (Saucedo-Lozano et al., 1999; Saucedo-Lozano y Chiappa-Carrara, 2000; Pérez-España et al., 2005; Raymundo-Huizar y Saucedo-Lozano, 2008; Flores-Ortega et al., 2010, 2011; Aguilar-Palomino y González-Sansón, 2010).

Las zonas costeras están sujetas a alteraciones por descargas de ríos y por desechos domésticos, que propician cambios súbitos en salinidad, sedimentos, concentración de materia orgánica, así como de otros contaminantes (Ambrogi et al., 1990; Pagola-Carte y Saiz-Salinas, 2002; Salen-Picard y Arlhac, 2002). La macrofauna bentónica responde a la alteración del hábitat y su diversidad se relaciona con la intensidad y frecuencia de estas alteraciones (Weston, 1990; Nickell et al., 2003; Austen y Widdicombe, 2006). Las instalaciones de maricultura también producen un enriquecimiento de materia orgánica que altera los fondos marinos y produce cambios en las comunidades bentónicas (Weston, 1990; Nickell et al., 2003; Sutherland et al., 2007; D'Amours et al., 2008; Tomassetti et al., 2009). En estos casos la fuente de materia orgánica proviene principalmente de restos de alimento suministrado y no ingerido, así como de las heces de los peces enjaulados, cuyo impacto es más evidente en el bentos, sobre todo, debajo de las jaulas y en sus inmediaciones (Sanz-Lázaro y Marín, 2011). Los patrones de diversidad y abundancia de la macrofauna bentónica en ambientes alterados fueron descritos con detalle en el modelo de Pearson y Rosenberg (1978), donde se describen los cambios y sucesiones de las comunidades macrofaunales a lo largo de un gradiente de alteración por enriquecimiento en materia orgánica. En este modelo se indica que cuando hay un aporte importante de materia orgánica se produce un incremento de la abundancia de especies oportunistas y un descenso de la biomasa, la riqueza específica y la biodiversidad, aunque a niveles intermedios de perturbación también se pueden producir incrementos en la riqueza específica y la diversidad (Connell, 1978; Pearson y Rosenberg, 1978).

En la bahía de Navidad se encuentran una jaula de engorde de peces, una descarga estacional (arroyo El Pedregal), una laguna costera abierta permanentemente y exposición al oleaje variable según la orientación de las marejadas. Por otra parte, las condiciones hidroclimáticas de la zona presentan un patrón estacional muy claro, con 2 épocas evidentemente diferenciadas. De junio a noviembre prevalecen condiciones tropicales debido a la influencia de la contracorriente norecuatorial (CCNE) y la época lluviosa (Kessler, 2006; Lavín et al., 2006). Durante la estación seca, de diciembre a mayo, prevalecen las condiciones subtropicales debido a la influencia de una rama de la corriente de California con flujo hacia el sur (Kessler, 2006; Lavín et al., 2006) y, sobre todo, a eventos de afloramiento costero entre enero y mayo (López-Sandoval et al., 2009), que producen una notable disminución de la temperatura superficial del agua del mar.

Teniendo en cuenta lo anterior, se plantea la hipótesis de que existen diferencias espaciales y temporales en las comunidades de macrofauna bentónica de la bahía y tienen correlación con las condiciones ambientales, en concreto con los aportes de sedimento y materia orgánica provenientes de las descargas de aguas continentales y la jaula de maricultura. El objetivo de este trabajo es describir y analizar los patrones de variación espacial y temporal de la macrofauna submareal en la bahía de Navidad y determinar su relación con el tamaño del sedimento y su contenido en materia orgánica.

Materiales y métodos

Área de estudio. El área de estudio se localizó en la bahía de Navidad, situada en el sur de la costa de Jalisco (situada entre los puntos 19°10'30" y 19°14'00" N, 104°45'00" y 104°41'30" O), en el Pacífico tropical mexicano y ha sido descrita con detalle por Flores-Ortega et al. (2010). Se realizó un muestreo mensual de marzo de 2007 a marzo de 2008, salvo en los meses de abril y octubre de 2007 y febrero de 2008, debido a condiciones meteorológicas adversas. Se seleccionaron 4 estaciones de muestreo, para tener representación de secciones de la bahía con diferentes condiciones ambientales (Fig. 1). La estación 1 se ubica en la parte norte, en las inmediaciones de una jaula de engorde de peces (principalmente huachinango, Lutjanus peru) con capacidad para 10 000 ejemplares, y de una descarga del arroyo El Pedregal, que aporta sedimentos y materia orgánica durante la época de lluvias. Con una superficie de 0.49 km², es la más protegida del oleaje y presenta el sedimento más fino. La estación 2 se encuentra inmediata hacia el sur, a 700-1 400 m de la descarga y la jaula, con 0.49 km², presenta mayor exposición al oleaje y sedimento fino. La estación 3 es la más extensa y expuesta al oleaje con 0.77 km², presenta el sedimento más grueso y se encuentra entre 1 400-2 800 m de la descarga y la jaula, y a 750-2 150 m de la desembocadura de la laguna de Barra de Navidad. La estación 4 se encuentra en las inmediaciones de la desembocadura de la laguna de Barra de Navidad, que se mantiene abierta permanentemente, a pesar de que el aporte de agua continental tiene un marcado carácter estacional; con 0.27 km², está más expuesta al oleaje y presenta sedimento grueso. En todas las estaciones de muestreo la profundidad osciló entre 8 y 15 m.

Muestreo y procesado de muestras. Cada mes se tomaron 3 réplicas por estación empleando un cuadrado de 50x50 cm y se aspiró todo el interior hasta una profundidad de 15 cm en el sedimento mediante una bomba de succión operada por buceadores. La bomba de succión consistió en un tubo de PVC (180 cm de largo y 7.5 cm de diámetro) unido a un difusor de acero inoxidable de 30 cm de largo y conectado a un compresor de aire. En el extremo del tubo se acopló una bolsa de 0.5 mm de luz de malla. Las muestras se tamizaron a bordo de la embarcación con una malla de 0.5 mm. Los organismos fueron separados en el laboratorio in vivo, siempre que fue posible. En caso contrario, las muestras se congelaron a -20° C hasta su separación. Una vez separados, los organismos fueron fijados en una solución de formol al 4% en agua de mar y, posteriormente conservados en etanol al 70%. Los individuos se identificaron hasta el nivel taxonómico más bajo posible. Se formaron 26 grupos taxonómicos funcionales, órdenes, infraórdenes, superfamilias y familias, de acuerdo con su movilidad, alimentación, tamaño, forma y posición en el sedimento (Hewitt et al., 2008), además del criterio taxonómico (Cuadro 3), para describir los patrones de la macrofauna relacionados con factores ambientales como enriquecimiento en materia orgánica o contaminantes (Olsgard et al., 1997, 1998).

Análisis granulométrico y contenido de materia orgánica. Para los análisis de granulometría y del contenido de materia orgánica en el sustrato, de noviembre de 2007 a marzo de 2008, se recolectaron en bolsa 2 muestras de 300 g de sedimento en cada estación de muestreo y, posteriormente, una vez descongeladas se realizaron análisis granulométricos de acuerdo con el método descrito por Bale y Kenny (2005), el cual consiste en secar la submuestra de 200 g de sedimento a 60° C y, después tamizarla en húmedo a través de una torre de tamices de distinta luz de malla, de 2 mm a 0.0625 mm, para separar la arena del lodo, esta última fracción fue considerada como la mezcla de limo y arcilla. Las fracciones retenidas en los tamices, grava y arena, fueron secadas en una estufa a 60° C y, posteriormente se determinó su peso en gramos. La fracción de lodo se estimó como la diferencia entre el peso seco total de la muestra y el peso seco de la suma de las fracciones de arena y grava. El tamaño de grano fue expresado como porcentajes de grava, de arena y de lodo según la escala de Udden/Wentworth (Wentworth, 1922). El contenido de materia orgánica en el sedimento fue calculado de acuerdo con el método de oxidación: se colocan 200 g de cada muestra en un vaso de precipitados, se agregaron 200 ml de peróxido de hidrógeno al 40% y se mantuvo en movimiento continuo hasta la degradación total de la materia orgánica, normalmente 3 días. Se secó la muestra en la estufa a 60° C y se calculó la diferencia de peso, que corresponde al contenido en gramos de materia orgánica. El peso se estandarizó para 100 g de sedimento y se expresó como porcentaje de materia ánica.

Análisis estadísticos. Para realizar los análisis estadísticos se agruparon los datos según los patrones principales hidroclimáticos de la región, que consisten en 2 periodos principales: 1) de junio a noviembre, periodo de influencia de la contracorriente norecuatorial y estación lluviosa, y que denominaremos "verano"; 2) de diciembre a mayo, periodo de influencia de la corriente de California y estación seca, y que denominaremos "invierno". En este caso, es de destacar la importancia de diferenciar bien los periodos de lluvias y secas, ya que la descarga de aguas continentales es una importante fuente de sedimentos y materia orgánica en el bentos, sobre todo, en zonas cercanas a la costa.

La densidad media mensual (individuos· m^-2) de cada grupo taxonómico fue estimada como el promedio de las 3 réplicas en cada estación. Con estas medias mensuales se estimaron las medias para cada periodo hidroclimático. Para comparar la densidad media entre estaciones de muestreo y periodos hidroclimáticos, se llevó a cabo un análisis de varianza (Anova) de 2 vías. Los datos de densidad fueron transformados a su raíz cuarta, debido a la gran cantidad de ceros y pocos valores elevados obtenidos (Quinn y Keough, 2002). Con esta transformación se consiguió una distribución normal de los datos (test de Kolmogorov-Sminrov, p> 0.2) y la homogeneidad de varianzas (test de Levene, p> 0.29). Posteriormente, se utilizó el test de Duncan para realizar las comparaciones múltiples.

Se llevaron a cabo análisis multivariantes para detectar los patrones en la distribución de las asociaciones de macrofauna. Los métodos de Clarke (1993) y el software PRIMER v. 6 fueron utilizados (Clarke y Gorley, 2006). La abundancia total de cada grupo taxonómico en cada mes de muestreo y estación, suma de las 3 réplicas de los cuadrados, fue utilizada en estos análisis. Se aplicó una transformación a raíz cuadrada a la suma de abundancias, ya que es la que mejores correlaciones obtiene entre matrices de similitud (tanto de variables ambientales como de grupos faunísticos) cuando trabajamos a nivel taxonómico de familia o niveles superiores (Olsgard et al., 1997; Karakassis y Hatziyanni, 2000; Heino, 2008). Para analizar gráficamente el patrón de asociación de grupos taxonómicos de acuerdo con la estación de muestreo y periodo hidroclimático, se utilizó el método de escalamiento multidimensional no métrico (nMDS) y un análisis de clasificación numérica jerárquica aglomerativa, promedio aritmético de grupos de pares no ponderados (UPGMA, por sus siglas en inglés), utilizando el índice de Bray-Curtis como medida de similitud. La significancia en la formación de los grupos de asociaciones en el dendrograma se determinó con la prueba Simprof, la cual realiza una prueba de permutaciones a cada nodo de un dendrograma para identificar patrones de asociación significativos en muestras a priori no estructuradas (Clarke et al., 2008). Se realizó un análisis de similitud de 2 vías (Anosim) para analizar diferencias en las comunidades entre estaciones de muestreo y periodos hidroclimáticos. Por último, se utilizó un análisis de porcentaje de similitud (Simper) para identificar las especies que más contribuyen a la similitud/disimilitud media entre las estaciones de muestreo y periodos hidroclimáticos que presentaban diferencias significativas en el Anosim.

Con fines comparativos entre las distintas estaciones de muestreo y épocas hidroclimáticas, y la comparación con el modelo de Pearson y Rosenberg y la "hipótesis de alteración intermedia" de Connell (1978), se estimaron el índice H' de Shannon-Wiener utilizando logaritmos en base e (Krebs, 1999) como medida de la diversidad, y el índice (J') de Pielou (1966) como medida de la equitatividad. Para cada mes y estación de muestreo se calcularon ambos índices con los grupos taxonómicos determinados previamente acumulados en las 3 réplicas y, posteriormente se estimaron los H' y J' medios para cada estación y época hidroclimática. Para comparar las diferencias de ambos parámetros entre estaciones de muestreo se utilizó un análisis no paramétrico de Kruskal-Wallis seguido de un test de comparaciones múltiples. Para comparar estas diferencias entre periodos hidroclimáticos en cada estación de muestreo se utilizó un análisis no paramétrico U de Mann-Whitney.

Las variables ambientales con granulometría y contenido de materia orgánica fueron analizadas también mediante un MDS, a partir de una matriz de similitud calculada por índices de la distancia euclídea, con datos sin transformar. Para establecer la relación entre estas variables y las comunidades bentónicas se utilizó el procedimiento Bio-Env, incluido en el software PRIMER (Clarke y Warwick, 2001; Clarke et al., 2008). Este procedimiento utiliza una correlación de rangos de Spearman para comparar las matrices de similitud de las variables ambientales y las comunidades de organismos. Para este análisis se agruparon los datos de invierno y verano de las comunidades macrofaunales y se estimó la media anual para la materia orgánica.

Condiciones ambientales. El análisis granulométrico confirmó el gradiente ascendente de tamaño de sedimento observado de norte a sur. Las 2 estaciones de muestreo del norte (1 y 2) presentaron un tamaño de grano fino, principalmente compuesto por limo y arcilla (en más de un 95% en la estación 1 y cercano al 90% en la estación 2) (Cuadro 1), mientras que las estaciones de muestreo del sur (3 y 4) presentaron tamaño de grano más grueso, principalmente arenas finas y medias, con un tamaño medio superior en la estación 4 (Cuadro 1).

El contenido en materia orgánica siguió un patrón similar. Las estaciones de muestreo del norte duplicaron el contenido en materia orgánica de las del sur, haciéndose más evidentes las diferencias durante el verano (Cuadro 1). En la parte sur la estación 4 presentó un porcentaje superior a la estación 3 (Cuadro 1).

Patrones espaciales y temporales de abundancia. Durante el periodo de muestreo se capturaron un total de 47 985 ejemplares pertenecientes a 61 taxones diferentes, que se agruparon posteriormente en 26. La densidad media total fue de 1 744.9 ind·m^-2, obteniéndose grandes variaciones de abundancia entre grupos taxonómicos y periodos hidroclimáticos, con oscilaciones entre 0 y 10 884 ind·m^-2. Se encontraron diferencias significativas en la densidad media entre las estaciones de muestreo y entre los 2 periodos hidroclimáticos del año (Anova, p< 0.005), pero no su interacción (Anova, p> 0.5), por lo que en todas las estaciones de muestreo la densidad media fue significativamente superior en verano (Fig. 2). Las estaciones 1 y 2 presentaron densidades medias significativamente mayores que las demás estaciones (test de Duncan, p< 0.005). Entre las estaciones 1 y 2 y las estaciones 3 y 4, no se encontraron diferencias significativas (Fig. 2).

Asociaciones de macrofauna bentónica. Las comunidades bentónicas de las 4 estaciones de muestreo presentaron una similitud bastante elevada, aun así el análisis de clasificación numérico y la prueba Simprof identificaron diferencias significativas entre las de la parte norte y sur de la bahía (Pi= 2.65, p< 0.01; Fig. 3). Las estaciones 1 y 2 presentaron diferencias significativas entre periodos hidroclimáticos (Pi= 1.82, p< 0.05), pero no dentro del mismo periodo hidroclimático (Fig. 3). El análisis de similitud (Anosim) reveló diferencias significativas entre estaciones de muestreo, pero no entre periodos hidroclimáticos (Cuadro 2). Las comunidades de las estaciones 1 y 2 no presentaron diferencias significativas, al igual que las estaciones 3 y 4, con lo que existe una clara diferenciación entre las comunidades del norte y sur de la bahía (Cuadro 2). El resto de combinaciones presentó diferencias significativas, especialmente entre las estaciones 1 y 2 con la 4 (Cuadro 2). De los 26 grupos taxonómicos considerados a lo largo de todo el periodo de muestreo para las estaciones del norte, 4 grupos acumularon más del 90% de la abundancia total: escafópodos, poliquetos, bivalvos y anfípodos (Cuadro 2), y tan sólo 2: escafópodos y poliquetos, acumularon el 87.8% en la estación 1, encontrándose elevadas densidades de los primeros (Cuadro 3). En las estaciones del sur no se produjo tanta dominancia y son 9 los grupos necesarios para acumular el 90% de abundancia: anfípodos, poliquetos, cumáceos, misidáceos, isópodos, bivalvos, escafópodos, picnogónidos y gasterópodos (Cuadro 3). Se observa que la abundancia de escafópodos es muy elevada en las estaciones del norte de la bahía, y en las del sur su abundancia muestra un brusco descenso, pasando de una media de 1 290. 1 ind·m^-2 en la estación 2 a 29 ind·m^-2 en la estación 3 (Cuadro 3).

El análisis Simper indica que los escafópodos son el grupo taxonómico que más contribuye a la disimilitud o diferenciación entre las estaciones de muestreo de las 2 partes de la bahía, entre las que el Anosim detectó diferencias significativas, contribuyendo entre el 30 y el 44% a la disimilitud (Cuadro 4). Los poliquetos son el segundo grupo que más contribuye a la disimilitud, con porcentajes entre el 8.7 y el 14.2, seguidos de los anfípodos, quetognatos, bivalvos, picnogónidos, gasterópodos e isópodos con menores porcentajes (entre el 4 y el 10%; Cuadro 4).

Diversidad. El índice de diversidad de Shannon-Wiener (H') y la equitatividad (J') presentaron prácticamente el mismo patrón espacial y temporal. La diversidad presentó valores bajos y similares entre las estaciones 1 y 2, y valores más elevados y similares entre las estaciones 3 y 4 (Cuadro 5). Considerando todo el periodo de muestreo, la estación 3 presentó un índice H' significativamente mayor a todas las demás (Kruskal-Wallis, p< 0.01; Cuadro 5), en invierno no se encontraron diferencias significativas entre las 4 estaciones y en verano la estación 4 presentó el H' más elevado (Kruskal-Wallis, p< 0.01). Salvo en la estación 4, en todas las demás la diversidad media durante el verano fue inferior a la del invierno, sobre todo, en la estación 1, donde la diferencia es mayor (1.0 vs. 1.6), aunque estas diferencias no fueron significativas (U Mann-Whitney, p> 0.05 en todos los casos).

Considerando todo el periodo de muestreo, los valores de equitatividad en las estaciones 3 y 4 fueron significativamente mayores que en las estaciones 1 y 2 (Kruskal-Wallis, p< 0.01; Cuadro 5). No se encontraron diferencias significativas entre estaciones en los valores de invierno (Kruskal-Wallis, p> 0.05), pero sí en verano, donde las estaciones 3 y 4 presentaron valores más elevados que las estaciones 1 y 2 (Kruskal-Wallis, p< 0.01). No hubo diferencias significativas entre los valores de equitatividad de invierno y verano en cada estación (U Mann-Whitney, p> 0.05), aunque en la estación 1 el valor medio de invierno es casi el doble que en verano y la diferencia está cerca de ser significativa (p= 0.08).

Relación entre variables ambientales y comunidades de macrofauna. Se encontró una relación significativa muy alta entre las comunidades de macrofauna y la granulometría y la cantidad de materia orgánica. Se obtuvo un coeficiente de correlación de Spearman r_s= 0.94 para el porcentaje de arena fina, al igual que para el porcentaje de limo grueso, ambos característicos de las estaciones del norte y del sur. Asimismo, se obtuvo el mismo nivel de correlación para las combinaciones % arena muy gruesa + % materia orgánica, o % arena muy gruesa + % arena fina, entre otras.

Discusión

La bahía de Navidad presentó 2 zonas bien diferenciadas en cuanto a las condiciones ambientales del fondo marino y a las comunidades de macrofauna que lo habitan. La mitad norte, más protegida del oleaje, se caracteriza por un sedimento fino y rico en materia orgánica; y la parte sur, más expuesta al oleaje, por un sedimento más grueso y menor contenido en materia orgánica. Aunque no se observa que los aportes de materia orgánica y sedimento fino afecten a la parte sur de la bahía, tanto el arroyo Pedregal como la laguna de Barra de Navidad, que está permanentemente comunicada con el mar, son fuentes cercanas de estos materiales. Estudios previos indican poca variación de salinidad y turbidez en las inmediaciones de la boca de la laguna en todo el ciclo anual (Sandoval-Rojo et al., 1988), por lo que el efecto del agua continental y la sedimentación afecta, sobre todo, al interior de la laguna, que junto con el manglar deben actuar como filtro de materia orgánica y sedimentos. La diferenciación sedimentaria entre las partes norte y sur de la bahía podría explicarse por el patrón de corrientes de la bahía. Un estudio de corrientes reciente (datos no publicados) indica que las corrientes dominantes en la estación 1 son de dirección predominante SE-NO, y las de la estación 2 son de dirección S-N, lo que impediría la dispersión de grandes cantidades de materia orgánica y sedimento hacia la parte sur de la bahía, que está mucho más expuesta al oleaje y lavaría gran parte de la materia orgánica y el sedimento más fino.

Existe la posibilidad de que las instalaciones de maricultura situadas en las inmediaciones de la estación 1 también fueran causantes de un enriquecimiento en materia orgánica, como se ha demostrado en otros estudios (Weston, 1990; Nickell et al., 2003; Sutherland et al., 2007; D'Amours et al., 2008; Tomassetti et al., 2009). La actividad máxima de la jaula de engorde se registró entre marzo y abril, cuando los ejemplares alcanzaron mayor tamaño y fueron cosechados. Sin embargo, la cantidad de materia orgánica fue muy superior durante la época de lluvias. También es posible que haya otros focos de contaminación no identificados en el área, pero en estos casos es muy difícil separar los efectos antrópicos de los naturales (Pires-Vanin et al., 2013). Las diferencias estacionales observadas entre periodos hidroclimáticos indican que la influencia de la descarga del río es más importante que otras posibles fuentes de enriquecimiento por materia orgánica.

El patrón estacional de abundancia encontrado en el presente trabajo ha sido observado en otras comunidades macrobentónicas en el Pacífico central mexicano (Ahumada-Sempoal et al., 2002) y en otras zonas del golfo de México (Hernández-Arana et al., 2003), donde la proliferación de ejemplares se debe a que el gran aporte de materia orgánica por parte de las descargas de los ríos en la época de lluvias produce un incremento en la biomasa (Escobar-Briones y Soto, 1997; Hernández-Arana et al., 2003; Rodríguez-Valencia, 2004), aunque en algunos casos la hipoxia causada por grandes descargas de los ríos (Gastón, 1985) o la acumulación masiva de sedimento tras una tormenta (Rodríguez-Valencia, 2004) pueden producir una reducción importante de la densidad del macrobentos.

Los patrones de variación espacial de la abundancia, la diversidad y el contenido en materia orgánica observados en el presente estudio concuerdan claramente con el gradiente de contaminación por materia orgánica descrita en el modelo de Pearson y Rosenberg (1978), y que ha sido documentado en numerosos casos de alteración por enriquecimiento en materia orgánica (Weston, 1990; Nickell et al., 2003; Kutti et al., 2007; D'Amours et al., 2008; Tomassetti et al., 2009) (Fig. 4). En este modelo se describen 4 estados de sucesión bentónica: normal, transitoria, contaminada y fuertemente contaminada, dependiendo del nivel de materia orgánica en el sedimento. En este gradiente la diversidad decrece cuando aumenta el contenido en materia orgánica, hasta que la fauna desaparece, cuando la contaminación es muy severa. La abundancia también se incrementa con la cantidad de materia orgánica (Fig. 4) y se produce una simplificación de la estructura trófica de la comunidad, que pasa a ser dominada por unas pocas especies oportunistas como los poliquetos o los anfípodos. Estos son de pequeño tamaño y detritívoros, presentan reproducción continua o semicontinua y están adaptados para colonizar rápidamente (Lenihan y Micheli, 2001; Rosa y Bemvenuti, 2006) y son competidores exitosos cuando la disponibilidad de materia orgánica es alta, pero pobres competidores cuando el alimento y el espacio es escaso, por lo que rápidamente son desplazados por otras especies cuando las condiciones ambientales se restablecen (Lenihan y Micheli, 2001).

Numerosos estudios acerca de la respuesta de las comunidades bentónicas ante un enriquecimiento de materia orgánica coindicen en destacar que los poliquetos, sobre todo, capitélidos y espiónidos son los oportunistas más abundantes en condiciones de mayor enriquecimiento (Pearson y Rosenberg, 1978; Bigot et al., 2006; Vivan et al., 2009; entre otros), aunque los anfípodos (Lenihan y Micheli, 2001; Brooks et al., 2003; Bigot et al., 2006; Vivan et al., 2009), isópodos y cumáceos (Lenihan y Micheli, 2001; Vivan et al., 2009) y algunos bivalvos (Brooks et al., 2003; Kutti et al., 2007; Vivan et al., 2009) también forman parte de esta fauna oportunista. En el presente estudio la parte de la bahía enriquecida con materia orgánica presentó elevada abundancia de poliquetos, pero una alta dominancia de escafópodos del género Dentalium. En la estación 1 representan casi el 80% de las capturas, con una densidad media de 1 831 ind·m^-2, y llegando a un máximo de 9 788 ind·m^-2. En ninguno de los estudios citados anteriormente se encontraron escafópodos entre la fauna de zonas enriquecidas con materia orgánica, ya sea debido a cultivos marinos o a actividad antrópica. Ríos-Jara et al. (2003) destacan la escasa literatura existente sobre escafópodos en las costas del Pacífico mexicano e indican que aparecen ocasionalmente en las muestras de sedimento marino. Los datos obtenidos en trabajos tanto del Pacífico como del Atlántico indican que son organismos poco abundantes (10-220 ind·m^-2) en la plataforma de Jalisco y Colima (Ríos-Jara et al., 2003) y todavía menos abundantes (4-34 ind·m^-2) en la plataforma de Yucatán (Hernández-Arana et al., 2003) y se encuentran principalmente en zonas profundas, entre 72-75 m de profundidad, aunque algunos se capturaron a 23 m) en Jalisco y Colima, y 80-140 m en Yucatán, por lo que las elevadas densidades obtenidas en el presente trabajo indican que se trata de un fenómeno excepcional. La biología de los escafópodos es aún hoy en día la menos entendida de los moluscos (Reynolds, 2002). Se sabe que son organismos excavadores semisésiles de pequeño tamaño, viven enterrados en el sedimento y son más abundantes en sedimentos finos, fangosos y fango-arcillosos. Son carnívoros que se alimentan principalmente de foraminíferos y otros microorganismos (Reynolds, 2002), aunque algunas especies no son especialistas y se alimentan de cualquier cosa que puedan capturar (Shimek, 1990). Las densidades obtenidas en el presente estudio son muy superiores a las registradas hasta el momento, por lo que suponemos que en la parte norte de la bahía de Navidad encuentran el tipo de sedimento y la cantidad de alimento óptimos para soportar estas densidades. Los foraminíferos son indicadores de contaminación (Eichler et al., 2012) o enriquecimiento por materia orgánica (Gómez-Noguera y Hendrickx, 1997), y se han encontrado elevadas densidades en lugares con sedimento fino y contaminados, por materia orgánica u otro tipo de contaminantes, tanto en zona intermareal (Hack et al., 2007) como submareal (Gómez-Noguera y Hendrickx, 1997; Eichler et al., 2012), por lo que la parte norte de la bahía de Navidad podría contener elevadas densidades de estos organismos. La densidad de los escafópodos fue muy inferior en la parte sur de la bahía y son el grupo taxonómico que marca las diferencias entre comunidades del norte y del sur. La parte sur de la bahía no está tan afectada por un gran aporte de materia orgánica y sedimento fino y presentó una fauna más diversa, aunque los grupos más abundantes fueron anfípodos, poliquetos, cumáceos, misidáceos, isópodos y bivalvos, la mayoría de ellos también considerados como oportunistas y encontrados en otras zonas tropicales afectadas moderadamente por enriquecimiento en materia orgánica (Bigot et al., 2006; Vivan et al., 2009; Pires-Vanin et al., 2013), por lo que podemos suponer que el enriquecimiento en materia orgánica afecta a toda la bahía en mayor o menor medida.

Recientes estudios sobre alimentación de peces de arrecife (Aguilar-Palomino y González-Sansón, 2010; Flores-Ortega et al., 2010) y demersales (Flores-Ortega et al., 2011) en el área de estudio indican que las zonas someras de la bahía de Navidad podrían ser una fuente importante de alimento para estas especies, especialmente la parte sur de la bahía (estaciones 3 y 4), ya que los crustáceos: anfípodos, decápodos e isópodos, de pequeño tamaño son una parte muy importante de la dieta de especies como la sarangola, Microlepidotus brevipinnis y el chivo, Mulloidichthys dentatus (Flores-Ortega et al., 2010), o las rayas Urobatis halleri, Urotrygon munda y Urotrygon rogersi, que se alimentan principalmente de decápodos y anfípodos (Flores-Ortega et al., 2011). El 62% de la abundancia de los grupos taxonómicos capturados en la estación 3, y el 69% en el caso de la estación 4, son presas potenciales de estas especies, lo que indicaría el valor de la bahía de Navidad, al menos en su parte sur, como fuente de alimento de peces importantes para la pesquería.

Agradecimientos

El presente trabajo fue financiado por la Universidad de Guadalajara. Los autores agradecen a Yael Morales-Jáuregui su apoyo en el trabajo de campo y la separación e identificación de organismos.

Aguilar-Palomino, B. y G. González-Sansón. 2010. Composición y superposición de dietas en cuatro especies de Diplectrum (Perciformes: Serranidae) en el Pacífico central mexicano. Revista de Biología Tropical 58:1189-1198.         [ Links ]

Ahumada-Sempoal, M. A., S. J. Serrano-Guzmán y N. Ruiz-García. 2002. Abundancia, estructura poblacional y crecimiento de Atrina maura (Bivalvia: Pinnidae) en una laguna costera tropical del Pacífico mexicano. Revista de Biología Tropical 50:1091-1100.         [ Links ]

Ambrogi, R., D. Bedulli y G. Zurijni. 1990. Spatial and temporal patterns in structure of macrobenthic assemblages. A three-year study in the Northern Adriatic Sea in front of the Po River Delta. Marine Ecology 11:25-41.         [ Links ]

Austen, M. C. y S. Widdicombe. 2006. Comparison of the response of meio- and macrobenthos to disturbance and organic enrichment. Journal of the Experimental Marine Biology and Ecology 330:96-104.         [ Links ]

Bale, A. J. y A. J. Kenny. 2005. Sediment analysis and characterization. In Methods for the study of marine benthos, A. Eleftheriou y A. McIntyre (eds.). Blackwell Scientific, Oxford. p. 43-86.         [ Links ]

Bigot, L., C. Conand, J. Amouroux, P. Frouin, H. Bruggemann y A. Gremare. 2006. Effects of industrial outfalls on tropical macrobenthic sediment communities in Reunion Island (Southwest Indian Ocean). Marine Pollution Bulletin 52:865-880.         [ Links ]

Blaber, S. J. M. 1997. Fish and fisheries of tropical estuaries. Chapman & Hall, Londres. 367 p.         [ Links ]

Brooks, K. M., A. R. Stierns, C. V. W. Mahnken y D. B. Blackburn. 2003. Chemical and biological remediation of the benthos near Atlantic salmon farms. Aquaculture 219:355-377.         [ Links ]

Clarke, K. R. 1993. Non-parametric multivariate analyses of changes in community structure. Australian Journal of Ecology 18:117-143.         [ Links ]

Clarke, K. R. y R. M. Warwick. 2001. Change in marine communities. An approach to statistical analysis and interpretation. PRIMER-E. Plymouth.         [ Links ]

Clarke, K. R. y R. N. Gorley. 2006. PRIMER v6. User manual and tutorial. PRIMER-E. Plymouth.         [ Links ]

Clarke, K. R., P. J. Somerfield y R. N. Gorley. 2008. Testing of null hypotheses in exploratory community analyses: similarity profiles and biota-environment linkage. Journal of the Experimental Marine Biology and Ecology 366:56-69.         [ Links ]

Connell, J. H. 1978. Diversity in tropical rain forests and coral reefs. Science 199:1302-1309.         [ Links ]

D'Amours, O., P. Archambault, C. W. McKindsey y L. E. Johnson. 2008. Local enhancement of epibenthic macrofauna by aquaculture activities. Marine Ecology Progress Series 371:73-84.         [ Links ]

Eichler, P. P, B. B. Eichler, B. S. Gupta y A. R. Rodrigues. 2012. Foraminifera as indicators of marine pollutant contamination on the inner continental shelf of Southern Brazil. Marine Pollution Bulletin 64:22-30.         [ Links ]

Escobar-Briones, E. G. y L. A. Soto. 1997. Continental shelf benthic biomass in the Western Gulf of Mexico. Continental Shelf Research 17:585-604.         [ Links ]

Flores-Ortega, J. R., E. Godínez-Domínguez, J. A. Rojo-Vázquez, A. Corgos, V. H. Galván-Piña y G. González-Sansón. 2010. Interacciones tróficas de las seis especies de peces más abundantes en la pesquería artesanal en dos bahías del Pacífico central mexicano. Revista de Biología Tropical 58:383-397.         [ Links ]

Flores-Ortega, J. R., J. A. Rojo-Vázquez, E. Godínez-Domínguez, G. González-Sansón, A. Corgos y Y. Morales-Jáuregui. 2011. Feeding habits of three round stingrays (Rajiformes: Urotrygonidae) in the central Mexican Pacific. Ciencias Marinas 37:279-292.         [ Links ]

Gaston, G. R. 1985. Effects of hypoxia on macrobenthos of the inner shelf off Cameron, Louisiana. Estuarine, Coastal and Shelf Science 20:603-613.         [ Links ]

Gillanders, B. M., K. W. Able, J. A. Brown, D. B. Eggleston y P. F. Sheridan. 2003. Evidence of connectivity between juvenile and adult habitats for mobile marine invertebrates. Marine Ecology Progress Series 247:281-295.         [ Links ]

Gómez-Noguera, S. E. y M. E. Hendrickx. 1997. Distribution and abundance of Meiofauna in a subtropical coastal lagoon in the Southeastern Gulf of California, Mexico. Marine Pollution Bulletin 34:582-587.         [ Links ]

Gunter, G. 1967. Some relationships of estuaries to the fisheries of the Gulf of Mexico. In Estuaries, G. H. Lauff (ed.). Publication 8, American Association for the Advancement of Science, Washington, D. C. p. 62-638.         [ Links ]

Hack, L. A., L. A. Tremblay, S. D. Wratten, A. Lister y V. Keesing. 2007. Benthic meiofauna community composition at polluted and non-polluted sites in New Zealand intertidal environments. Marine Pollution Bulletin 54:1801-1812.         [ Links ]

Heino, J. 2008. Influence of taxonomic resolution and data transformation on biotic matrix concordance and assemblage-environment relationships in stream macroinvertebrates. Boreal Environment Research 13:359-369.         [ Links ]

Hernández-Arana, H. A., A. A. Rowden, M. J. Attrill, R. M. Warwick, R. M. y G. Gold-Bouchot. 2003. Large-scale environmental influences on the benthic macroinfauna of the Southern Gulf of Mexico. Estuarine, Coastal and Shelf Science 58:825-841.         [ Links ]

Hewitt, J. E., S. F. Thrush y P. D. Dayton. 2008. Habitat variation, species diversity and ecological functioning in a marine system. Journal of the Experimental Marine Biology and Ecology 366:116-122.         [ Links ]

Karakassis I. y E. Hatziyanni. 2000. Benthic disturbance due to fish farming analyzed under different levels of taxonomic resolution. Marine Ecology Progress Series 203:247-253.         [ Links ]

Kessler, W. S. 2006. The circulation of the Eastern Tropical Pacific: a review. Progress in Oceanography 69:181-217        [ Links ]

Krebs, C. J. 1999. Ecological methodology. Second edition. Benjamin Cummings. Menlo Park, California. 620 p.         [ Links ]

Kutti, T, P. K. Hansen, A. Ervik, T. Høisæther y P. J. Johannessen. 2007. Effects of organic effluents from a salmon farm on a fjord system. II. Temporal and spatial patterns in infauna community composition. Aquaculture 262:355-366.         [ Links ]

Lavín, M. F., P. C. Fiedler, J. A. Amador, L. T. Ballance, J. Färber-Lorda y A. M. Mestas-Núñez. 2006. A review of Eastern Tropical Pacific oceanography: summary. Progress in Oceanography 69:391-398.         [ Links ]

Lenihan, H. S. y F. Micheli. 2001. Soft sediment communities. In Marine community ecology, M. Bertness, M. E. Hay y S. D. Gaines (eds.). Sinauer, Sunderland, Massachusetts. p. 253-287.         [ Links ]

López-Sandoval, D. C., J. R. Lara-Lara, M. F. Lavín, S. Álvarez-Borrego y G. Gaxiola-Castro. 2009. Primary productivity observations in the Eastern Tropical Pacific off Cabo Corrientes, Mexico. Ciencias Marinas 35:169-182.         [ Links ]

Nickell, L. A., K. D. Black, D. J. Hughes, J. Overnell, T. Brand, T. D. Nickell, E. Breuer y S. M. Harvey. 2003. Bioturbation, sediment fluxes and benthic community structure around a salmon cage farm in Loch Creran, Scotland. Journal of the Experimental Marine Biology and Ecology 285-286: 221-233.         [ Links ]

Olsgard, F., P. J. Somerfield y M. R. Carr. 1997. Relationships between taxonomic resolution and data transformations in analyses of macrobenthic community along an established pollution gradient. Marine Ecology Progress Series 149:173-181.         [ Links ]

Olsgard, F., P. J. Somerfield y M. R. Carr. 1998. Relationships between taxonomic resolution, macrobenthic community patterns and disturbance. Marine Ecology Progress Series 172:25-36.         [ Links ]

Pagola-Carte, S. J. y I. Saiz-Salinas. 2002. Changes in the sublittoral faunal biomass induced by the discharge of a polluted river along the adjacent rocky coast (N. Spain). Marine Ecology Progress Series 212:13-27.         [ Links ]

Pearson, T. H. y R. G. Rosenberg. 1978. Macrobenthic succession in relation to organic enrichment and pollution of the marine environment. Oceanography and Marine Biology: An Annual Review 16:229-311.         [ Links ]

Pielou, E. C. 1966. The measurement of diversity in different types of biological collections. Journal of Theoretical Biology 13:131-144.         [ Links ]

Pérez-España, H., M. Saucedo-Lozano y A. R. Raymundo-Huizar. 2005. Trophic ecology of demersal fishes from the Pacific shelf off central Mexico. Bulletin of Marine Science 77:19-31.         [ Links ]

Pires-Vanin, A. M. S., E. Arasaki y P. Muñiz. 2013. Spatial pattern of benthic macrofauna in a sub-tropical shelf, São Sebastião Channel, southeastern Brazil. Latin American Journal of Aquatic Research 41:42-56.         [ Links ]

Quinn, G. P. y M. J. Keough. 2002. Experimental design and data analysis for biologists. Cambridge University Press, Nueva York. 537 p.         [ Links ]

Raymundo-Huizar, A. R. y M. Saucedo-Lozano. 2008. Hábitos alimentarios del pez triglido Prionotus ruscarius (Gilbert y Starks, 1904) durante 1996, en las costas de Jalisco y Colima, México. Revista de Biología Marina y Oceanografía 43:7-15.         [ Links ]

Reynolds, P. D. 2002. The scaphopoda. Advances in Marine Biology 42:37-236.         [ Links ]

Ríos-Jara, E., M. Pérez-Peña, E. López-Uriarte y E. Juárez-Carrillo. 2003. Moluscos escafópodos de la campaña Atlas V (Plataforma continental de Jalisco y Colima, México) a bordo del B/O El Puma. Ciencias Marinas 29:67-76.         [ Links ]

Rodríguez-Valencia, J. A. 2004. Respuesta de los poliquetos bentónicos a la variabilidad ambiental y condiciones El Niño en Bahía Petacalco, Guerrero, México. Ciencias Marinas 30:515-526.         [ Links ]

Rosa, L. C. y C. E. Bemvenuti. 2006. Temporal variability of the estuarine macrofauna of the Patos Lagoon, Brazil. Revista de Biología Marina y Oceanografía 41:1-9.         [ Links ]

Salen-Picard, C. y D. Arlhac. 2002. Long-term changes in a Mediterranean benthic community: relationships between the Polychaete assemblages and hydrological variations of the Rhône River. Estuaries 25:1121-1130.         [ Links ]

Sandoval-Rojo, L., F. J. Flores-Verdugo, U. Zaragoza-Araujo, J. W. Day, Jr. y A. Estrada-Mercado. 1988. Phytoplankton productivity in Barra de Navidad coastal lagoon on the Pacific coast of Mexico. Revue d' Hydrobiologie Tropicale 21:101-108.         [ Links ]

Sanz-Lázaro, C y A. Marín. 2011. Diversity patterns of benthic macrofauna caused by marine fish farming. Diversity 3:176-199.         [ Links ]

Saucedo-Lozano, M., G. González-Sansón y X. Chiappa-Carrara. 1999. Alimentación natural de Lutjanus peru (Nichols y Murphy, 1922) (Lutjanidae: Perciformes), en las costas de Jalisco y Colima, México. Ciencias Marinas 25:381-400.         [ Links ]

Saucedo-Lozano, M. y X. Chiappa-Carrara. 2000. Alimentación natural de juveniles de Lutjanus gutattus (Pisces: Lutjanidae), en las costas de Jalisco y Colima, México. Boletín del Centro de Investigaciones Biológicas 34:159-180.         [ Links ]

Shimek, R. L. 1990. Diet and habitat utilization in a Northeastern Pacific Ocean scaphopod assemblage. American Malacological Bulletin 7:147-169.         [ Links ]

Sutherland, T. F., C. D. Levings, S. A. Petersen, P. Poon y B. Piercey. 2007. The use of meiofauna as an indicator of benthic organic enrichment associated with salmonid aquaculture. Marine Pollution Bulletin 54:1249-1261.         [ Links ]

Tomassetti, P., E. Persia, I. Mercatali, D. Vani, V. Marussso y S. Porrello. 2009. Effects of mariculture on macrobenthic assemblages in a Western Mediterranean site. Marine Pollution Bulletin 58:533-541.         [ Links ]

Vivan, J. M., M. Di Domenico y T. C. Marques de Almeida. 2009. Effects of dredged material disposal on benthic macrofauna near Itajaí Harbour (Santa Catarina, south Brazil). Ecological Engineering 35:1435-1443.         [ Links ]

Wentworth, C. K. 1922. A scale of grade and class terms for clastic sediments. The Journal of Geology 30:377-392.         [ Links ]

Weston, D. P. 1990. Quantitative examination of macrobenthic community changes along an organic enrichment gradient. Marine Ecology Progress Series 61:233-244.         [ Links ]