Printer Friendly

El colapso ambiental en el rio Patia, Colombia: variaciones morfologicas y alteraciones en los ecosistemas de manglar.

The environmental collapse in the Patia River, Colombia: morphological variations and impacts on mangrove ecosystems

INTRODUCCION

Las intervenciones de origen humano sobre los flujos terrestres de agua y sedimentos son en la actualidad un factor relevante dentro de la evolucion de los sistemas deltaicos y generan cambios significativos en cortos periodos de tiempo. La alteracion de los cursos fluviales, asi como las variaciones en los aportes fluviales de agua, sedimentos y nutrientes, son algunos de los cambios introducidos por las actividades humanas que tienen repercusiones sobre las zonas costeras (Vorosmarty & Sahagian, 2000; Wolanski & Spagnol, 2000; Syvitski et al., 2005). El deterioro de la dinamica natural de los deltas implica, ademas de graves impactos ambientales, serias problematicas socio-economicas ligadas al detrimento de los ecosistemas costeros y a la disminucion del recurso pequero, entre otros servicios directos e indirectos de que hace uso la poblacion (Overeem & Syvitski, 2009).

Los efectos de las presiones antropicas sobre la dinamica natural de los sistemas deltaicos se observan en las zonas costeras de todo el mundo. En deltas como el Mississippi, Ebro o Nilo, la progresiva regulacion de la descarga de los rios mediante el establecimiento de represas se ha traducido en procesos destructivos de la linea de costa y la consecuente perdida de terrenos altamente productivos debido al deficit de sedimentos fluviales (Coleman et al., 1998; Sanchez-Arcilla et al., 1998; Stanley & Warne, 1998). En el caso contrario, el aumento en la produccion de sedimentos al interior de las cuencas debido a la deforestacion, adecuacion de terrenos para actividades agropecuarias, extraccion minera e industrializacion, se reflejan en procesos acelerados de progradacion de la zona litoral y en el deterioro de los ecosistemas costeros, debido al incremento de nutrientes y turbidez de aguas estuarinas y costeras. Ejemplos de deltas a nivel mundial con estas caracteristicas son los rios Amarillo, Fly y Magdalena (Yang et al., 1999; Wolanski & Spagnol, 2000; Restrepo et al., 2006).

En el delta del rio Patia, costa del Pacifico colombiano (Fig. 1), las principales problematicas ambientales se derivan de la captura inducida del caudal del rio Patia por el rio Sanquianga mediante el canal Naranjo. La construccion de este canal artificial en 1973 por la sociedad Maderas Naranjo, permitio facilitar el transporte de madera desde el rio Patia hasta el rio Satinga. El pequeno canal de 1,5 m de ancho por 1,2 m de profundidad y 1,3 km de largo (Castillo, 1995) comenzo a ampliarse debido a las continuas crecientes del rio Patia hasta alcanzar un ancho de mas de 400 m en la actualidad. Los cambios hidraulicos del rio Patia y la formacion del nuevo brazo Patianga han originado graves problematicas ambientales como procesos acelerados de erosion y sedimentacion en la linea de costa y en los cauces, cambios en la salinidad de los distributarios del delta, y alteraciones en los ecosistemas estuarinos y de manglar.

Los impactos que el cambio de curso del rio Patia ha generado sobre los ecosistemas de manglar, son de especial interes para la zona norte del delta debido a la presencia del Parque Nacional Natural Sanquianga (PNNS), la mayor reserva de manglares del Pacifico norte sudamericano (Fig. 1). En esta region litoral, con un area de 800 [km.sup.2], se encuentra aproximadamente el 20% de los manglares del Pacifico colombiano y habitan ca. 8000 personas que dependen de los recursos que ofrecen estos ecosistemas (UAESPNN, 2005). Los efectos del transvase se evidencian principalmente en cambios en la composicion y distribucion de las coberturas de manglar y en la colonizacion de especies de agua dulce debido a la disminucion de la salinidad en la zona norte del delta (Restrepo & Cantera, 2011).

A los problemas descritos anteriormente, se suma la falta de documentacion completa y de mediciones cuantitativas de los cambios que ha presentado el sistema deltaico, y tambien, la carencia de informacion acerca de los factores fisicos y sociales que interactuan en el delta. De los estudios y trabajos que se han realizado sobre este sistema (Gomez, 1986; Moreno, 2003; FUNDEPUBLICO, 1994; IDEAM, 1996, 1997; MINAMBIENTE, 1996), no todos estan disponibles para el publico en general y los que se encuentran disponibles no documentan de manera clara la magnitud de los impactos ambientales derivados del canal Naranjo (Defensoria del Pueblo, 2009). Esto dificulta establecer el alcance real de los efectos que se presentan en la zona e impide desarrollar lineamientos y politicas de gestion adecuadas.

En este trabajo se presenta un diagnostico de la situacion ambiental del delta como producto del desvio del rio Patia y una discusion preliminar de las repercusiones socio-economicas. El diagnostico incluye: (1) la estimacion de los cambios morfologicos de la linea de costa del delta mediante el calculo de las tasas de erosion y acrecion para el periodo 1986-2008; (2) el analisis de la distribucion longitudinal de la salinidad en los estuarios del lobulo norte del delta, region de la reserva de manglares PNNS; y (3) el calculo de los cambios en la cobertura vegetal del PNNS con enfasis en las alteraciones del ecosistema de manglar, y en los procesos de sedimentacion y cierre de canales.

[FIGURA 1 OMITIR]

El delta del rio Patia

Ubicado al suroccidente de Colombia, el delta del rio Patia (Fig.1) abarca una franja litoral de mas de 120 km y se extiende en una amplia llanura aluvial con un area emergente de aproximadamente 1700 [km.sup.2] (Restrepo & Cantera, 2011). El rio Patia tiene una longitud de 415 km y cuenta con un area de drenaje de 23700 [km.sup.2] (Restrepo & Lopez, 2008b), que contiene mas de 1000 cursos de agua que drenan la cordillera occidental de los Andes y parte de la central. En su desembocadura, el rio Patia descarga un caudal de 1320 [m.sup.3] [s.sup.-1] y una carga de sedimentos de 27 Mton [ano.sup.-1] para una produccion de sedimentos de 1500 ton [km.sup.-2] [ano.sup.-1], la tasa de erosion mas alta de cualquier rio drenando el margen occidental de sudamerica (Restrepo & Cantera, 2011).

Al igual que en el resto de la region del Pacifico colombiano, el clima en el delta del Patia se caracteriza por altos valores de temperatura, precipitacion y humedad. El promedio anual de precipitacion oscila entre 3500 y 6500 mm, siendo las temporadas lluviosas entre abril-junio y septiembre-noviembre y las menos lluviosas en diciembre-marzo y julio-agosto. La temperatura promedio en la zona es de 26[grados]C y la humedad relativa promedio anual entre 80 y 90% (Eslava, 1993; Del Valle-Arango, 1994).

El regimen mareal en el delta es de tipo semidiurno puro y rango mesomareal con amplitudes entre 2,6 y 3 m. El oleaje presenta condiciones intermedias de energia, con olas dominantes del tipo swell provenientes en su mayoria del suroeste, con alturas promedio de 1,4 m y periodos de 6 s (Restrepo & Lopez, 2008a). Los vientos dominantes son de direccion suroeste y noroeste con una velocidad promedio de 2,5 m [s.sup.-1] (UAESPNN, 2005).

Localizado en una zona de alta actividad tectonica, debido a la triple colision de las placas Nazca, Cocos y Suramerica, este sistema deltaico esta expuesto a la ocurrencia de movimientos tectonicos y fenomenos relacionados como tsunamis, licuacion de suelos y subsidencia (INGEOMINAS e IGAC, 2008). El ultimo terremoto de gran magnitud registrado en el litoral pacifico colombiano, ocurrido en 1979 al suroeste de la costa de Tumaco (Fig. 1), genero un tsunami que destruyo seis caserios y afecto todo el frente deltaico del Patia. Tasas de subsidencia de hasta 1,6 m se registraron en el litoral del delta despues del tsunami (Herd et al., 1981). Ademas de los eventos sismicos, el area del delta tambien se encuentra afectada por varios sistemas de fallas activas que ejercen control estructural sobre el litoral y modelan su evolucion. El principal sistema de fallas es Remolino-El Charco, con direccion NE-SW paralelo a la linea de costa (Fig.1) (Gomez, 1986).

Evolucion del delta del Patia

La evolucion del delta del rio Patia ha estado marcada por un alto dinamismo producto de diversos procesos tectonicos, marinos y fluviales. En la formacion del delta se destacan dos periodos donde se generaron gran parte de sus caracteristicas geomorfologicas actuales. En el Pleistoceno tardio (Fig. 2a), el rio Patia drenaba hacia el norte, formando un lobulo paralelo a la cordillera occidental, siendo sus distributarios los actuales rios Tapaje, Satinga, Sanquianga y La Tola. De igual forma, el rio Telembi y los rios Iscuande y Guapi construian lobulos al sur y al norte del rio Patia, respectivamente. El lobulo formado por el rio Patia en este periodo es denominado delta antiguo (Gomez, 1986) y se caracteriza por canales de marea amplios y varios remanentes de cordones litorales sucesivos que demuestran un alto nivel de progradacion (Gomez, 1986; Moreno, 2003).

En el Holoceno temprano, los levantamientos finales de la cordillera Occidental reactivaron la tectonica de la zona a traves de la falla de Remolino-El Charco. Los movimientos estructurales obligaron al rio Patia a migrar progresivamente hacia el sur dejando en el trayecto cauces abandonados y lagunas. Durante el mismo periodo, el levantamiento de la serrania de Guandipa ejercio un control estructural en el suroriente del rio Telembi, obligandolo a fluir hacia el norte y a labrar su cauce a traves de las rocas sedimentarias de la serrania. En el Holoceno tardio (Fig. 2b), el rio Patia es finalmente capturado por el rio Telembi en el sector de Fatima, y ambos rios fluyen hacia el oeste a traves de la serrania y forman el lobulo de Salahonda (Fig. 2a). Este lobulo deltaico, denominado delta reciente (Gomez, 1986), se caracteriza por tener canales estrechos y un solo frente de cordones litorales (Gomez, 1986; Moreno, 2003).

La evolucion del delta del Patia en el Holoceno tardio ha estado caracterizada por la actividad tectonica en una zona litoral inestable, con procesos co-sismicos de subsidencia y levantamiento de la franja costera. Un analisis sobre la reconstruccion del desarrollo de la cuenca del Patia al norte del delta, junto con los cambios vegetacionales, estudios basados en registros de polen, diatomeas y caracteristicas sedimentarias (Velez et al., 2001), sugieren que: (1) el lobulo norte del delta esta influenciado por la actividad fluvial entre 3850 y 2700 anos AP; (2) despues del anterior periodo, el delta esta mas aislado de los aportes fluviales y solo algunos pulsos fluviales interrumpidos se presentan entre 2010 y 1430 anos AP; (3) en esta epoca, los manglares se desarrollan mas desde la franja litoral hacia las partes lagunares aguas arriba; (4) entre los anos 1430 y 810 AP se registra un avance de la franja litoral hacia el mar, debido posiblemente a un levantamiento estructural de origen tectonico; (5) entre 810 y 580 anos AP, el polen de manglar disminuye considerablemente y la costa evidencia un retroceso que sugiere eventos de subsidencia tectonica; y (6) para el periodo entre los anos 580 y 300 AP, no se registran mas periodos fluviales y la franja de manglar avanza de nuevo hacia el litoral, sugiriendo procesos de levantamiento tectonico; evidencia de polen de arboles como Cecropia indican el inicio de actividades humanas en el delta desde hace 300 anos aproximadamente. Esta evaluacion del desarrollo del Patia en el Holoceno reciente muestra la alternancia de procesos constructivos y destructivos de origen tectonico en la configuracion del delta, los cuales se han visto interrumpidos y modificados por el desvio de origen antropico del rio Patia hace 39 anos.

Actualmente, el delta se encuentra en un nuevo periodo de desarrollo debido a la construccion del canal Naranjo y al transvase del caudal del rio Patia por su antiguo cauce, el rio Sanquianga (Fig.2d). El cambio en la direccion del flujo se vio favorecido por el desnivel existente entre los dos rios, el cual permitia, en epocas de creciente, que el rio Patia se desbordara y parte de su caudal fluyera por el rio Sanquianga. Esta nueva configuracion del rio Patia ha permitido la reactivacion del delta antiguo, el cual se denomina en la actualidad lobulo norte o activo. Este se extiende desde la bocana Amarales hasta la bocana Pasacaballos (Fig.1b). Por otra parte, el delta reciente, llamado actualmente lobulo sur o inactivo, comprende el area entre la bocana Pasacaballos y la bocana Salahonda (Fig.1b). El cambio en las condiciones hidraulicas del rio Patia ha generado problematicas de tipo ambiental, economico y social con efectos opuestos en los lobulos del delta. Los principales impactos se derivan de procesos acelerados de erosion y sedimentacion en la linea de costa y en los cauces.

[FIGURA 2 OMITIR]

MATERIALES Y METODOS

Se emplearon diversas fuentes de informacion para construir la base de datos del delta, la cual abarca aspectos geofisicos, ecologicos y socio-economicos del area de estudio. Los datos espaciales y alfanumericos recopilados se manejaron mediante Sistemas de Informacion Geografica (SIG) para facilitar su almacenamiento, integracion y posterior analisis.

La cartografia base y tematica utilizada se obtuvo de distintas instituciones y comprende: mapas del Instituto Geografico Agustin Codazzi a escala 1:100000 del ano 1990 para toda el area del delta y a escala 1:25000 de 1988 correspondientes a la zona norte del delta; mapas topograficos base de los servicios en linea ESRI (ArcGis Online Map Services) (ESRI, 2009); cartografia base y tematica del Parque Nacional Natural Sanquianga de la Unidad Administrativa Especial del Sistema de Parques Nacionales Naturales; y cartografia base y tematica del Sistema de Informacion Geografica para la Planeacion y el Ordenamiento Territorial del Instituto Geografico Agustin Codazzi.

Para la evaluacion multitemporal de las variaciones morfologicas en el delta del rio Patia se emplearon, junto con la cartografia, imagenes satelitales Landsat TM y ETM+ de 30 m de resolucion espacial del Global Land Cover Facility de la Universidad de Maryland y del Servicio Geologico de Estados Unidos de distintas fechas entre 1986 y 2001, imagenes ASTER de la coleccion TerraLook de 15 m de resolucion para los anos 2007 y 2008 del Servicio Geologico de Estados Unidos y fotografias aereas de 1962 (escala 1:60000) y 1983 (escala 1:40000) del Instituto Geografico Agustin Codazzi. Adicionalmente, se realizaron inspecciones en campo de levantamiento de puntos GPS para zonas de interes y de registro fotografico de las diferentes coberturas vegetales encontradas en el PNNS y del estado de la linea de costa en algunas zonas del delta.

Toda la cartografia disponible fue transformada a un mismo sistema de coordenadas geograficas (World Geodetic System 1984), para ser empleada como punto de partida para la georreferenciacion de los demas datos espaciales. El pre-procesamiento de las imagenes satelitales incluyo: correccion geometrica, mejoramiento, evaluacion de estadisticos basicos, composiciones de distintas bandas, generacion de mascaras y delimitacion de regiones de interes.

Calculo de tasas de erosion y acrecion litoral

Para calcular las tasas de erosion y acrecion costeras se generaron mosaicos de las imagenes satelitales y se digitalizo la linea de costa. Los montajes se elaboraron teniendo en cuenta el porcentaje de nubosidad, el estado de la marea (bajamar) y las fechas de las imagenes. De esta manera se obtuvo cuatro mosaicos correspondientes a las fechas de 1986, 1997-1998, 2001 y 2007-2008. Las imagenes empleadas para la digitalizacion de la linea de costa se georreferenciaron mediante la cartografia IGAC 1:100000 la cual cubre la totalidad del area de estudio. El registro de las imagenes Landsat se hizo mediante la rectificacion imagen a mapa de las escenas del 10 de mayo de 1986 por ser las que mejor coincidian con la cartografia 1:100000. A partir de estas imagenes se registraron, de imagen a imagen, las demas fechas procurando siempre que el error de la raiz cuadrada media (RMS) no sobrepasara medio pixel (15 m). Por otra parte, las imagenes ASTER se ajustaron adecuadamente a la cartografia y por lo tanto, no se georreferenciaron para evitar distorsiones.

Para delimitar la linea de costa se empleo como indicador de la posicion real del borde costero la linea de vegetacion, ya que este indicador biologico se ve menos afectado por las variaciones a corto plazo ocasionadas por las condiciones de marea y del nivel relativo del mar de la costa pacifica (Morton & Speed, 1998). A pesar que pueden existir desigualdades entre la posicion del borde de vegetacion y la linea de costa, debido a que la vegetacion responde a diversas condiciones ambientales (Morton & Speed, 1998), este indicador es practico, considerando la informacion disponible para el area de estudio, asi como las restricciones en el trabajo de campo, ya que el borde de vegetacion se registra consistentemente a lo largo del delta y es posible identificarlo a partir de las imagenes satelitales (Pajak & Leatherman, 2002).

Para facilitar la identificacion de la linea de vegetacion y disminuir los errores durante la digitalizacion, se emplearon composiciones en falso color (RGB) 752 y se calculo el indice de vegetacion de diferencia normalizada (Normalized Difference Vegetation Index--NDVI) para los mosaicos de imagenes Landsat TM y ETM+. Las imagenes Aster del 2007 y 2008 son distribuidas por el USGS en formato JPG y solo cuentan con las bandas 1, 2 y 3 por lo cual se debio hacer la identificacion de la linea de vegetacion solo con la composicion en color verdadero (RGB) 321.

Para evaluar el cambio en la linea de costa se empleo el programa Digital Shoreline Analysis System (DSAS) (Thieler et al., 2009), que calcula las tasas de erosion y acrecion (End Point Rate-EPR) al medir la distancia, en sentido perpendicular, entre lineas de distintas fechas y luego dividirla entre los anos transcurridos. En este estudio las lineas base se trazaron siguiendo la forma de la linea de costa tierra adentro y la distancia entre transectos se establecio en 100 m.

La evaluacion de exactitud de los resultados se hizo teniendo en cuenta dos fuentes principales de error que generan incertidumbre en la posicion real de las lineas de costa: la georreferenciacion de las imagenes y la digitalizacion de la linea de vegetacion. El error de georreferenciacion ([E.sub.g]) se deriva del error de la raiz cuadrada media (RMS) obtenido durante la correccion geometrica de las imagenes, mientras que el error de digitalizacion ([E.sub.d]) se estimo a partir de valores de otros estudios referenciados por Moore (2000). Siguiendo la metodologia descrita por Hapke et al. (2006), el error total en la posicion de las lineas de costa ([E.sub.lc]) se calculo mediante la raiz cuadrada de la suma de los cuadrados de los errores (Ecuacion 1).

[E.sub.lc] = [raiz cuadrado de [E.sub.g.sup.2 + [E.sub.d.sup.2]] (1)

Distribucion longitudinal de la salinidad en los estuarios del lobulo activo del delta

Las observaciones de variacion longitudinal de la salinidad para establecer el limite de incursion salina a lo largo de los estuarios del lobulo deltaico activo en la region norte, se llevaron a cabo en los tres principales canales: Sanquianga, Guascama y Amarales. Para comparar esta distribucion salina con los estuarios del lobulo sur inactivo, tambien se realizaron perfiles en el canal del estuario Pasacaballos (Fig. 1b).

Se efectuaron muestreos longitudinales a lo largo del eje principal de cada canal distributario, desde la boca estuarina hacia el apice del delta, en estaciones secuenciales para la toma de perfiles verticales de salinidad. El proposito fue lograr una caracterizacion cuasi-sinoptica de la incursion de salinidad a traves de cada sistema distributario. Estas medidas se efectuaron usando un perfilador CTD YSI-Cast Away en tres comisiones de campo durante noviembre 2010 y febrero-mayo 2011. Las mediciones se iniciaron el primer y segundo dia de mareas "Spring Tides". Las amplitudes estuvieron entre 2,9 y 3,0 m, con valores dentro del promedio del rango mesomareal del delta del Patia.

El protocolo fue el mismo para todos los muestreos y consistio en la realizacion de perfiles longitudinales estuario arriba, los cuales se iniciaron en la boca de cada sistema una hora antes de marea alta y consecuentemente terminaron unas 1,5 h en la parte superior de cada distributario. Con el fin de comprobar los limites maximos de intrusion salina en cada estuario, se realizaron cuatro muestreos adicionales en terreno de un proyecto alterno entre WWF y el Ministerio de Ambiente y Desarrollo Sostenible de Colombia, durante junio-julio y septiembre-noviembre de 2012.

Clasificacion de coberturas vegetales

El analisis de cambio de la cobertura vegetal del PNNS se realizo a partir de imagenes satelitales del 10 de mayo y 23 de marzo de 1986 y del 10 de abril de 2001. Para el ano 1986 se emplearon dos imagenes debido a que no se encontro una sola con toda el area del parque despejada. Las imagenes satelitales fueron recortadas teniendo en cuenta el lindero del PNN Sanquianga y posteriormente, se georrefenciaron usando la cartografia IGAC 1:25000. Este proceso se realizo mediante la rectificacion imagen a mapa de la escena del 10 de mayo de 1986 por ser de una fecha cercana a la cartografia y por tener mayor numero de puntos de control. A partir de esta imagen se registraron, de imagen a imagen, las demas fechas disponibles para la zona del parque, procurando siempre que el error de la raiz cuadrada media (RMS) no sobrepasara medio pixel (15 m).

El metodo usado para el analisis de cambio de las coberturas fue el de comparacion post-clasificacion, mediante el cual se evaluaron los cambios en las coberturas luego de que cada una de las imagenes fue clasificada de manera independiente. De esta manera, es posible trabajar con datos de dos sensores diferentes (Landsat TM y ETM+), y no es necesario realizar la normalizacion radiometrica de las imagenes (Armenakis et al., 2003; Jensen, 2005). El esquema de clasificacion empleado fue el Corine Land Cover (CLC), adaptado para Colombia, para la produccion de cartografia escala 1:100000 de cobertura y uso del suelo (IDEAM, 2010).

El proceso de clasificacion fue el mismo para las tres imagenes y los datos de entrada fueron las bandas espectrales a excepcion de las termicas y la pancromatica. Inicialmente se realizo una clasificacion no supervisada empleando el metodo ISODATA (Iterative Self-Organizing Data Analysis Technique), y para facilitar el proceso de agrupacion del algoritmo, se aplicaron a cada una de las imagenes mascaras para el agua y nubes, las cuales se establecieron a partir de la banda 5 (rango de datos de 0 a 18 ND) y la banda 2 (rango de datos de 0 a 60 ND) respectivamente. El resultado fue la clasificacion de seis tipos de coberturas que se identificaron empleando fotografias aereas de 1983 para las escenas de 1986 y con verificacion en campo para la imagen de 2001. Para las seis coberturas establecidas, se delimitaron areas de entrenamiento a partir de las cuales se realizo una clasificacion supervisada empleando la metodologia de maxima similitud. Para evaluar cualitativamente la exactitud de los resultados, estos fueron nuevamente confrontados con fotografias aereas y revision en campo. A partir de esta revision se establecio de manera general, la correspondencia entre las coberturas a partir de la clasificacion y su distribucion real.

RESULTADOS

Variaciones morfologicas en el delta del rio Patia

Para el litoral Pacifico, los principales factores que generan variaciones en la linea de costa son los niveles de energia mareal y oleaje, subsidencia cosismica de los suelos deltaicos, aumentos temporales del nivel del mar ocasionados por los eventos El Nino y la Oscilacion del Sur (ENOS) y las alteraciones antropicas de la descarga de sedimentos en la zona costera (Gonzalez et al., 2002; Restrepo, 2008). En el delta del Patia, la desviacion de los sedimentos ocasionada por el canal Naranjo se ha convertido en uno de los principales agentes de cambio en la evolucion de la linea de costa. Al incrementarse el aporte de sedimentos en el lobulo norte, la respuesta esperada es la acrecion de la linea de costa, mientras que en el caso contrario, lobulo sur, se espera retroceso de la linea frontal del delta por la disminucion de los sedimentos fluviales.

Las tasas de erosion y acrecion calculadas para toda la linea de costa del delta del Patia ejemplifican los principales cambios morfologicos en los sectores norte y sur para los tres periodos de analisis 1986-1998, 1998-2001 y 2001-2008. A partir de los promedios de las tasas de cambio y del porcentaje de transectos en erosion, acrecion o equilibrio (Tabla 1, Fig. 3), es posible establecer a grandes rasgos el comportamiento de los lobulos deltaicos. En la zona norte se observa que para el periodo 1986-1998, la tasa de cambio es negativa con la mayoria de los transectos presentando erosion. Esto concuerda con el proceso de retroceso documentado por Martinez et al. (1995), quienes registran un aporte insuficiente de arenas, lo que indica que el efecto del transvase del Patia al Sanquianga aun no era evidente en la zona costera. Para el periodo de 1998-2001, la tasa de retroceso de la linea de costa disminuye y se registra un mayor numero de transectos en equilibrio. Este resultado se puede deber al aumento del aporte de sedimentos por el brazo Patianga (union del rio Patia con el Sanquianga a traves del canal Naranjo, Fig. 1b) y al inicio de un proceso de progradacion por procesos de depositacion de sedimentos, como lo demuestra el aumento de los bancos de arena en las bocanas y la colmatacion y cierre de esteros. Por ultimo, entre 2001 y 2008, la linea de costa norte presenta una tasa positiva con la mayoria de los transectos en acrecion, lo que se puede interpretar como la estabilizacion del lobulo norte, donde nuevas zonas de bajos internareales y playas han empezado a ser colonizados por la vegetacion, indicando condiciones de estabilidad del terreno.

El lobulo sur presenta, entre 1986-1998, tasas de cambio negativas con la mayoria de los transectos en erosion, pero con valores mucho menores comparados con los de la zona norte. Este resultado concuerda con lo descrito por Martinez et al. (1995), quienes encuentran evidencia de erosion para la zona sur pero tambien formacion activa de dunas debido al aporte de sedimentos litorales. Para el periodo 1998-2001 se registran tasas positivas con la mayoria de los transectos en equilibrio. Esta aparente progradacion no se debe al avance de la linea de costa sino al alargamiento de las islas barreras y al proceso de cierre de las bocanas debido a la disminucion del caudal por el brazo Patia. Para el periodo 2001-2008 se evidencia el retroceso de la linea de costa, con la mayoria de las transectas en estado de erosion y un mayor cierre de las bocanas con formacion de bajos intermareales al frente de los canales de ingreso (Fig. 4).

En general, el frente deltaico de todo el sistema Patia-Sanquianga esta en retroceso para el periodo 1986-1998. Este comportamiento erosivo se pudo deber al efecto post-tsunami 1979, cuando el sistema litoral del delta evidencio tasas de subisidencia de hasta 1,6 m a lo largo de una franja costera de 200 km (Herd et al., 1981).

Los sectores de la linea de costa que evidenciaron los cambios mas criticos muestran, las tendencias descritas anteriormente para cada lobulo (Fig. 3). En la zona norte el retroceso de la costa alcanzo tasas maximas de hasta 40 m [ano.sup.-1] (Fig. 3c), las cuales concuerdan con las observaciones realizadas en campo de turbas de manglar expuestas, restos de troncos y raices de manglares en la zona intermareal y tramos desprovistos de cobertura por remocion de la vegetacion (Fig. 5). La erosion de la linea de costa en este sector seria compensada en un futuro por aportes de sedimentos fluviales siempre que se mantengan las tasas de acrecion registradas para el periodo 2001-2008, las cuales alcanzan, en algunos tramos, valores maximos de hasta 30 m [ano.sup.-1] (Fig. 3d).

En el lobulo sur, como se menciono anteriormente, se presenta un retroceso del frente litoral y procesos de acrecion cerca de las bocanas debido al cierre de los canales de acceso por disminucion del caudal del brazo Patia. Las tasas maximas de retroceso en esta zona alcanzan 70 m [ano.sup.-1] (Fig. 3g), siendo el evento mas destacado el rompimiento de la isla barrera de San Juan. Esta isla barrera se vio afectada por los terremotos y posteriores tsunamis de 1906 y 1979, lo que obligo a la poblacion a reubicarse al interior del estero detras de la barra de arena donde se encontraba el caserio inicialmente. La ruptura definitiva de la isla ocurrio, al parecer en 1995, luego de lo cual parte de la poblacion se desplazo hacia otros caserios. Las tasas de erosion, en el sitio de la ruptura, alcanzaron a 78,67 m [ano.sup.-1] en el periodo 2001-2008 y el canal de apertura se amplio a mas de 2 km (Fig. 4). Durante el mismo periodo, al norte de la isla barrera, se aprecio el cierre de la bocana Guandipa debido a la formacion y posterior vegetalizacion de barras de arena al frente de los canales de ingreso. Este mismo proceso se presento en las demas bocanas del sector sur del delta, con colmatacion de los canales y avance de la cuna salina. El proceso de abandono del lobulo sur ha derivado en perdida de terrenos de playa, cultivos y ecosistemas de manglar debido al retroceso de la linea de costa y a la salinizacion de los cursos de agua (Defensoria del Pueblo, 2009).

Evaluacion de exactitud de las tasas de cambio

Al confrontar las tasas de cambio calculadas con las observaciones realizadas en campo, con los testimonios de los habitantes de la region y con estudios anteriores, se encuentra que los valores obtenidos reflejan adecuadamente los patrones de variacion de la linea de costa para los periodos estudiados. Adicionalmente, a las apreciaciones cualitativas realizadas, se evaluo la exactitud de las tasas de cambio para establecer de manera cuantitativa, el nivel de error derivado de la incertidumbre en la posicion de la linea de costa digitalizada para cada periodo.

Los valores de error obtenidos, a partir de la evaluacion de exactitud, se presentan en la Tabla 2. El error promedio, para todos los periodos fue de [+ o -]25 m, indicando que la distancia horizontal entre un punto sobre una linea de costa digitalizada y su posicion real esta en un rango maximo de 25 m. Teniendo en cuenta que el error de la raiz cuadrada media para cualquier cartografia no puede ser mayor de 0,3 mm a la escala del mapa (resolucion 64 de 1994 del Instituto Geografico Agustin Codazzi), el error maximo permitido para los datos derivados de imagenes Landsat, cuya escala aproximada es de 1:100000 (Lopez-Vizoso, 1986), seria de 30 m. Dado que los errores obtenidos no superan este valor y que las distancias de retroceso o acrecion calculadas son, en general, mucho mayores que los errores totales encontrados, es posible concluir que las lineas de costa estimadas son aceptables y por lo tanto, las tasas de cambio establecidas a partir de estas son una buena aproximacion a la realidad.

Impacto en los ecosistemas de manglar del Parque Nacional Sanquianga

El incremento del aporte fluvial de agua y sedimentos en el rio Sanquianga ha generado alteraciones en las condiciones fisicas y quimicas del sistema estuarino del lobulo norte del delta. Los procesos de sedimentacion y cierre de canales, aumento de la turbidez por solidos suspendidos, erosion de las bancas de los distributarios activos, disminucion de la salinidad y variacion en la concentracion de nutrientes son algunos ejemplos de las modificaciones que se han presentado a raiz del canal Naranjo.

El efecto del transvase sobre los ecosistemas de manglar del PNNS se evidencia principalmente en los cambios en la composicion y distribucion de las coberturas vegetales ocurridos durante las ultimas cuatro decadas (Restrepo & Cantera, 2011). Mediante la evaluacion de las variaciones de la cobertura de manglar del parque, entre 1986 y 2001, se realiza un acercamiento preliminar a la problematica ocasionada por la alteracion del sistema estuarino debido a la reactivacion del lobulo norte.

[FIGURA 3 OMITIR]

[FIGURA 4 OMITIR]

[FIGURA 5 OMITIR]

Del proceso de clasificacion de las imagenes de 1986 y 2001 se establecieron seis coberturas en el area del PNN Sanquianga (Figs. 6, 7):

i) Sistema de playas y bajos de arena intermareales.

ii) Pantanos de lodo intermareal en margenes de lagunas estuarinas y canales distributarios.

iii) Asociacion de plantas herbaceas, vegetacion de firmes, bosques de natal y claros expuestos por deforestacion.

iv) Bosques de manglar de baja densidad y/o reducida altura.

v) Bosques de manglar denso y/o altura superior a 20 m.

vi) Bosques de condiciones fluviales y suelos firmes con asociaciones caracteristicas de especies de natales y guandales.

El analisis de cambio de las coberturas vegetales se enfoco en los ecosistemas de manglar y en las asociaciones de plantas herbaceas y vegetacion de agua dulce. Esta ultima cobertura se incluyo en el analisis debido a que presento una gran variacion entre las dos fechas por la alteracion en las condiciones de salinidad de los distributarios activos (Figs. 6b, 6i), logrando desplazar en algunas zonas al manglar.

Para ejemplificar los cambios encontrados en el PNNS, se definieron tres zonas representativas (Fig. 7): la bocana del distributario activo Sanquianga, que moviliza la mayor cantidad de caudal y sedimentos y, por lo tanto, presenta las mayores alteraciones hidricas (Fig. 7a); la bocana Amarales, desembocadura del rio Aguacatal, donde las condiciones de caudal y salinidad no se han visto afectadas considerablemente por el transvase, por lo cual los factores ambientales se han mantenido relativamente constantes (Fig. 7b); y por ultimo, el apice del delta actual aguas abajo de Bocas de Satinga, donde el rio se divide en los distributarios y se presentan procesos acelerados de sedimentacion (Fig. 7c).

Como se menciono anteriormente, las alteraciones en el sistema estuarino del lobulo norte son mas evidentes en el principal distributario activo, el Sanquianga. Antes de la afectacion por el canal Naranjo, este cauce se caracterizaba por drenar lagos internos con bajos aportes fluviales de agua y sedimentos y por una alta influencia mareal con intrusion salina de hasta 35 km al interior de la planicie costera. El incremento en el aporte de sedimentos y la disminucion de la salinidad en este cauce han resultado en el cierre de canales interdistributarios y en el aumento de las asociaciones de vegetacion de agua dulce. Del analisis de cambio en la zona de la bocana Sanquianga (Fig. 7a), se destaca que en un periodo de 15 anos la cobertura de manglar alto-denso se ha incrementado por encima del manglar bajo-abierto, que puede ser el resultado de procesos de colonizacion de los pantanos de lodo intermareales que se desarrollan en las margenes de los canales y la posterior repoblacion de los esteros una vez que se cierran (Fig. 5). Durante el mismo periodo, en la zona de la bocana se observa la expansion de la vegetacion de agua dulce en casi el doble de hectareas (Tabla 3), proceso que se puede deber al reemplazo de la cobertura de manglar por ecosistemas de natal (Mora megitosperma) y naidizal (Euterpe cuatrecasana) junto con sus respectivas especies asociadas (UAESPNN, 2005).

No existen datos de mediciones de salinidad en los distributarios deltaicos del Patia. Para obtener los "paleo" limites de intrusion salina, especialmente para las epocas anteriores a la diversion del rio Patia en el rio Sanquianga, estos fueron comprobados en cada distributario con expertos de la region, principalmente los guias de Parques Nacionales. Ademas, por medio de analisis geomorfologicos y de asociaciones vegetales, visualizadas con imagenes de satelite y fotos aereas anteriores a 1990, los "paleo" limites de intrusion salina fueron marcados aproximadamente.

En la bocana Amarales (Fig. 7b), las condiciones estuarinas se han conservado estables debido al escaso caudal proveniente del rio Patianga que es evacuado por el rio Aguacatal. La influencia mareal y la intrusion salina se mantienen en este cauce, tal como se comprobo en campo mediante las mediciones longitudinales de salinidad (Fig. 6i). El patron de distribucion y las proporciones de hectareas de las coberturas de manglar estudiadas no muestran modificaciones drasticas entre las dos fechas. Al igual que en la bocana Sanquianga, esta zona presenta un aumento de la vegetacion de agua dulce en casi el triple de hectareas (Tabla 3), que en este caso se asociaria principalmente a procesos de deforestacion o establecimientos de cultivos en zonas de firmes o de manglares.

Al evaluar las transformaciones ocurridas en el apice del delta entre 1986 y 2001 (Fig. 7c), es evidente el cambio en la configuracion del cauce del rio Sanquianga ocasionado por la erosion de los margenes y los procesos acelerados de sedimentacion. La drastica disminucion de la salinidad y el crecimiento de barras y bajos de lodos debido al aporte de sedimentos, ha permitido la colonizacion y predominio de vegetacion de agua dulce en lo que antes eran ecosistemas de manglar. La variacion en la extension de las coberturas es notable para esta zona, con reducciones de 29% del area del manglar alto, 33% del manglar bajo y un aumento de 42% en el area de la vegetacion de agua dulce (Tabla 3). A partir de la valoracion en campo se pudo constatar la transformacion de esta zona, donde se destaca el predominio del helecho Acrostrichium aureum junto con otras especies arbustivas oportunistas, que sirven como bioindicadores de la magnitud del cambio observado en el area, ya que son especies que colonizan de manera rapida los espacios dejados por el manglar (Tavera, 2009; Restrepo & Cantera, 2011).

[FIGURA 6 OMITIR]

[FIGURA 7 OMITIR]

De acuerdo con expertos locales y funcionarios del PNNS (Wilfrido Ibarbo, com. pers.), el efecto del incremento de aportes fluviales en la ecologia de los manglares del Parque se ha incrementado en magnitud durante la ultima decada. Eventos no antes vistos como defoliacion masiva de cientos de hectareas debido a la plaga estacional de un gusano de la familia Saturnidae, como tambien, la aparicion de tumores "cancerigenos" en las raices y ramas del manglar, se comenzaron a observar hace ca. 15 anos. Ademas, en el apice del delta actual, una region extensa que estaba dominada por manglares Rhizophora, ha sido reemplazada por pastos caracteristicos de ambientes fluviales como Paspalum spp. y Panicum spp., conformando la parte exterior y frontal de los canales, seguidos por comunidades de juncos de agua Eichornia spp. Dicho cinturon es a su vez seguido de comunidades de la cana brava Gynerum sagittatum. Este patron vegetal es caracteristico de distributarios activos y de alta descarga fluvial, y ademas, su presencia ha sido reportado en deltas de Tabasco, Mexico (Thom, 1967), Patia y Mira, en la costa Pacifica de Colombia (West, 1957), delta del Atrato, sobre la costa Caribe de Colombia (Vann, 1959), y es similar al reportado en deltas tropicales como el Purari, Papua-Nueva Guinea (Conn, 1983).

Incursion longitudinal de la salinidad en el lobulo norte del delta, reserva natural PNNS

La distribucion longitudinal de la salinidad varia en principio segun la descarga fluvial y controla los procesos de inundacion, composicion de los suelos y distribucion de la vegetacion en los margenes de los canales distributarios del delta del Patia. Previo al cambio de curso del rio Patia hacia el rio Sanquianga, el marco vegetal en el lobulo norte del delta se caracterizaba por manglares densos y zonas estuarinas de amplia longitud, con limites de incursion salina hasta 35 km aguas arriba. Posterior al aumento de la descarga fluvial del rio Sanquianga a traves del canal Naranjo, que se acentuo mas a inicios de los anos 90', las condiciones de aportes fluviales al PNNS se incrementaron exponencialmente. Unos estuarios con formaciones de manglar de gran desarrollo se vieron expuestos al aumento de aportes de agua y sedimentos del Patia. En otras palabras, los antes estuarios Sanquianga y Guascama se convirtieron en clasicos canales distributarios de sistemas deltaicos. Las medidas longitudinales de salinidad a lo largo de estos distributarios activos indica la clara dependencia del patron de descarga fluvial en el plano deltaico norte (Fig. 6).

Los analisis de imagenes aereas de 1962 y cartografia tematica de formaciones vegetales del PNNS en 1989, muestran que la interfase de salinidad de 1 incursionaba aproximadamente 35 km aguas arriba por el distributario Sanquianga (Fig. 6a). Actualmente, la interfase de salinidad de 0 cerca al fondo de los canales se encontro a 13, 15 y 32 km en Sanquianga (Figs. 6f, 6h), Guascama (Figs. 6c, 6f) y Amarales (Fig. 10i), respectivamente. Estas variaciones de distribucion longitudinal de salinidad evidencian el efecto del incremento de los aportes fluviales en el rio Sanquianga despues del trasvase del Patia, cambios que han dado como resultado el movimiento aguas abajo de las fronteras de incursion de los frentes salinos. Por el contrario, el distributario Pasacaballos en el lobulo sur e inactivo del delta, que presenta ausencia de flujos de agua dulce durante ocho meses al ano, posee aguas mas salinas con amplia distribucion longitudinal de salinidad (Fig. 6b).

En el estuario Amarales (Fig. 6a), que no tiene influencia de las descargas del Sanquianga, las condiciones salinas son tipicas de estuarios del Pacifico colombiano con rango mesomareal (Restrepo & Kjerfve, 2002). Este canal esta caracterizado por: (1) salinidades neta-promedio en el perfil de profundidad de mas de 20; (2) salinidades superiores a 21 en condiciones de marea baja; (3) no presentar areas de erosion y alto hidrodinamismo en los margenes de los canales; (4) colonizacion e intrusion de bosques de manglar hasta 30 km aguas arriba; y (5) presencia de cinturones de manglar de gran desarrollo (Rhizophora) a lo largo del limite de las lagunas estuarinas y canales mareales, donde los suelos son lodos finos distribuidos por corrientes mareales.

DISCUSION

Variaciones morfologicas en el delta del rio Patia

Los procesos destructivos y constructivos que se presentan en los lobulos del delta del rio Patia evidencian la importancia de las presiones humanas dentro de la evolucion de estos sistemas. La expansion demografica y el desarrollo economico intensivo tanto al interior de las cuencas de drenaje como en las zonas costeras alteran el transporte de agua, sedimentos y nutrientes, ocasionan que muchos sistemas deltaicos a nivel mundial se encuentren alejados de su estado natural (Syvitski et al., 2005; Ericsson et al., 2006; Syvitski & Saito, 2007).

En el caso del lobulo sur, los procesos erosivos de la linea de costa son similares a las fases destructivas que experimentan los deltas que tienen deficit de sedimentos debido a la construccion de represas al interior de sus cuencas fluviales o a la desviacion de los cauces para usos agricolas o industriales. Dos importantes deltas a nivel mundial afectados por este tipo de procesos destructivos son el Mississippi y el Nilo. En el delta del Mississippi, la construccion de 65.000 represas en su cuenca hidrografica ha derivado en altas tasas de erosion en la franja litoral, las cuales son intensificadas por la subsidencia de los suelos deltaicos debido a la extraccion de petroleo en la zona costera (Coleman et al., 1998; Syvitski & Saito, 2007). En el rio Nilo, la construccion de represas para controlar las inunda-ciones anuales y el desarrollo de sistemas de irrigacion para agricultura ha generado erosion acelerada del frente litoral con enderezamiento de la linea de costa, aumento de la intrusion salina y perdida de areas de pantano y lagunas (Stanley & Warne, 1998).

Los procesos constructivos del lobulo norte se asemejan a las fases de progradacion de ciertos deltas, donde la cantidad de sedimentos que llega a la costa se incrementa debido a actividades como agricultura, ganaderia, mineria y explotacion forestal. Deltas a nivel mundial con estas caracteristicas son el Amarillo (Huanghe) y el Fly. El delta del rio Amarillo, ubicado en la costa oriental de China, es formado por el segundo rio mas largo de ese pais, el cual transporta una carga de sedimentos seis veces mayor que la del rio Mississippi debido a procesos de deforestacion y erosion del suelo, dando como resultado tasas de progradacion de la linea de costa de hasta 4 km [ano.sup.-1] (Yang et al., 1999; Syvitski & Saito, 2007). De manera similar, las actividades humanas en la cuenca del rio Fly han incrementado las tasas de erosion del suelo y en consecuencia ha aumentado la cantidad de sedimentos que llega a aguas estuarinas y costeras. La alteracion de las propiedades fisicas y biologicas del agua debido al exceso de sedimentos degrada los ecosistemas marinos debido a las alteraciones en las cadenas benticas y al sofocamiento del plancton y de corales (Wolanski & Spagnol, 2000).

Impacto en los ecosistemas de manglar del Parque Nacional Natural Sanquianga (PNNS)

Las problematicas ambientales encontradas en el PNNS siguen el patron predominante de degradacion que presentan los ecosistemas de manglar a nivel mundial. El desarrollo urbanistico e industrial en las zonas costeras, la reclamacion de tierras para agricultura y acuicultura, la sobre-explotacion del manglar, contaminacion y alteracion de los ciclos terrestres del agua, junto con otros impactos de origen antropico, han causado que el manglar sea uno de los ecosistemas mas amenazados en el mundo con tasas de reduccion anual de aproximadamente 2,1% global y hasta 3,6% en America (Valiela et al., 2001). En algunos casos, los efectos de las actividades humanas no se traducen en perdida de extension de la cobertura de manglar, sino en degradacion de la composicion floristica y biodiversidad del ecosistema, dando como resultado la dominancia de especies menos vulnerables que se adaptan facilmente a diferentes condiciones ambientales (Dahdouh-Guebas et al., 2005).

La importancia del manglar en la costa del Pacifico colombiano va mas alla de los servicios ecologicos y economicos que estos ofrecen, sino que tambien constituyen una importante defensa costera contra la erosion, tormentas e impacto de inundacion por tsunamis. Por esto el estudio de las implicaciones ambientales a largo plazo que tendra el transvase del rio Patia sobre los ecosistemas de manglar, no solo del PNNS, sino en toda la linea de costa del delta, debe ser primordial dentro de los planes de manejo y las politicas de desarrollo sostenible de la zona.

Por ultimo, es importante tener en cuenta que, desde el punto de vista cualitativo, la clasificacion de coberturas obtenidas para el PNNS se adapta adecuadamente a lo registrado durante los reconocimientos en campo y en las aerofotografias, por lo cual, su uso es valido como un acercamiento preliminar a los cambios de la vegetacion. Sin embargo, ya que no se cuenta con una evaluacion de exactitud cuantitativa, no es recomendable emplear los valores de hectareas obtenidos como datos de entrada para otro tipo de calculos debido a que se desconoce el nivel de error que estos contienen.

Implicaciones socio-economicas

A pesar de la riqueza natural de la region del Pacifico Colombiano y de su gran potencial economico favorecido por su posicion geoestrategica, las condiciones de competitividad y desarrollo de la zona presentan un gran atraso en comparacion con las demas regiones del pais, como se evidencia en las graves deficiencias en la prestacion de los servicios de salud, educacion y la baja cobertura de servicios publicos.

Al evaluar los indicadores de calidad de vida y necesidades basicas insatisfechas de los municipios del area del delta del rio Patia, se observa que para el ano 2005, el promedio del indice de calidad de vida fue de 43,7% frente al promedio nacional de 78,8%, una diferencia de mas de 30 puntos porcentuales, mientras que el indice de necesidades basicas insatisfechas para ese ano fue de 80,9%, en comparacion con el 27,8% a nivel nacional, lo que indica que la mayoria de la poblacion es pobre y vive en condiciones precarias (DNP, 2010; SIGOT, 2011).

A las problematicas sociales de la region se suma la economia netamente extractiva que se basa en actividades del sector primario como la pesca artesanal, agricultura y explotacion forestal, las cuales se caracterizan por una baja productividad con fines de subsistencia de la poblacion o de intercambio de productos. Los establecimientos industriales o de comercio presentes en la zona son en su mayoria de capital foraneo y se dedican principalmente a la explotacion pesquera o maderera y al aprovechamiento de la palma africana y el palmito (Bravo-Pazmino, 1998; DNP, 2010).

Otros factores socio-economicos que han cobrado importancia en la region durante los ultimos anos, son el conflicto armado y el narcotrafico. Las caracteristicas geograficas de la zona y su posicion estrategica como frontera maritima, la han convertido en un area ideal para el establecimiento de cultivos ilicitos y en un corredor de trafico de drogas que es ampliamente disputado por los distintos grupos ilegales (UNODC y DNE, 2011).

El desarrollo de cultivos ilicitos en las cuencas de los rios Patia y Telembi se ha convertido en una fuente de ingresos para la poblacion debido a la falta de oportunidades laborales que ofrecen las actividades economicas tradicionales. Segun la Oficina de las Naciones Unidas Contra la Droga y el Delito (UNODC, 2011), el departamento de Narino se mantiene, desde el 2003, como uno de los departamentos con mayor area sembrada de coca, llegando a ocupar el primer lugar en el 2010 con el 26% del total del area sembrada del pais. De la extension total de cultivos de coca en el departamento, el 90% se ubica en la zona del Pacifico siendo Tumaco el municipio con la mayor cantidad de area sembrada tanto a nivel departamental como nacional (8% del total nacional) y la mayor produccion de cocaina del pais, con un 21% de la produccion nacional (UNODC, 2011; UNODC y DNE, 2011).

El incremento acelerado de los cultivos ilicitos en la region Pacifica ha tenido repercusiones ambientales como la deforestacion de areas boscosas o de manglar, la contaminacion de la tierra y los cuerpos de agua superficiales por el uso intensivo de quimicos, el aumento de la erosion del suelo por la perdida de la cobertura vegetal, y en consecuencia, el incremento en la produccion de sedimentos al interior de las cuencas hidrograficas. A esta actividad ilegal tambien se asocian ciertas problematicas sociales como el desplazamiento forzado, la corrupcion administrativa, violencia, perdida de la gobernabilidad, abandono de los sistemas tradicionales de produccion, reclutamiento forzado, homicidios selectivos y masacres, entre otros (UAESPNN, 2005; INVEMAR et al., 2006).

Las distintas problematicas, mencionadas anteriormente, que configuran la dinamica socio-economica y demografica de la poblacion del delta del rio Patia, se ven agravadas por las transformaciones ambientales y ecologicas que se presentan en el rio Patia a raiz del canal Naranjo. La alta dependencia de los habitantes hacia los recursos naturales para subsistir, su fragilidad social y deficit de desarrollo economico, los hace vulnerables a los cambios del entorno, ya que, al no contar con suficientes herramientas politicas y de organizacion comunitaria, son incapaces de plantear alternativas de manejo para las nuevas condiciones que se presentan en el territorio.

En el sector norte del delta, la erosion de las margenes del rio Sanquianga, debido al aumento de caudal, ha derivado en la reduccion de la superficie de tierra cultivable, asi como en danos en propiedades de poblaciones riberenas. En Bocas de Satinga, cabecera municipal de Olaya Herrera y actual apice del delta, se han perdido mas de 2000 casas en los ultimos 15 anos, alcanzando en el 2008 una tasa de hasta tres viviendas destruidas por dia. A lo anterior se suma el deterioro de la infraestructura del municipio como las calles, el acueducto, centros de salud y locales comerciales, entre otros (Alcaldia Municipal Olaya Herrera, 2012).

De igual forma, las alteraciones en las condiciones de salinidad, el aumento de los sedimentos suspendidos en los ecosistemas estuarinos y la perdida de areas de manglar, han tenido como consecuencia la disminucion de recursos pesqueros importantes como la piangua (Anadara tuberculosa y Anadara similis), el cual es un molusco con gran demanda en Ecuador y que constituye una importante fuente de ingreso para cerca de 9.040 personas en el departamento de Narino (UAESPNN, 2005; WWF, 2005 en Delgado et al., 2010).

Estas circunstancias han llevado a que se generen dinamicas migratorias en la zona, con la reubicacion de poblados riberenos y el desplazamiento de la poblacion en busca de mejores condiciones economicas. Los movimientos de la poblacion se evidencian en el aumento acelerado del numero de habitantes en el area del PNN Sanquianga, que paso de 3.000 habitantes en 1977 a 8.075 en el 2005, y en la generacion de disputas entre las comunidades del sector por el aprovechamiento de los recursos naturales (Bravo-Pazmino, 1998; UAESPNN, 2005).

Por otro lado, el rio Sanquianga ha cobrado importancia en la region como via de transporte fluvial, lo que ha derivado en el establecimiento de aserrios en el sector norte del delta. El crecimiento de la actividad maderera, en detrimento de los bosques de guandal, ha favorecido los procesos migratorios de la poblacion hacia esta zona, que en la actualidad concentra la mayor parte de la infraestructura para el procesamiento de la madera, siendo Bocas de Santinga, con 24 aserraderos, el principal productor del pais (UAESPNN, 2005; Defensoria del Pueblo 2009; Alcaldia Municipal Olaya-Herrera, 2012).

En el sector sur del delta, el abandono progresivo del Brazo Largo del rio Patia ha causado un notable deterioro de la economia local. La fuerte erosion costera, debida a la disminucion del caudal del rio, ha destruido areas de manglar con importantes recursos forestales y pesqueros, al igual que tierras de cultivo y poblados costeros, como es el caso de Salahondita y Hojas Blancas, dos pequenos poblados del municipio Francisco Pizarro, que fueron destruidos completamente por la accion del oleaje y su poblacion se vio obligada a reubicarse (INVEMaR et al., 2006).

Ademas, el aumento de la salinidad del agua ha afectado la productividad del suelo y el recurso pesquero y por ende el rendimiento de las cosechas ha disminuido y los artes de pesca de la poblacion han cambiado. En los municipios de Mosquera y Francisco Pizarro, se reporta una reduccion progresiva de los volumenes de captura de la pesca artesanal, que se traduce en el detrimento economico de la poblacion que depende de esta actividad. Estos eventos han llevado a los habitantes de la zona a desplazarse hacia otros lugares, en busca de distintas formas de sustento, o a explotar areas alejadas de sus poblados, lo que deriva en conflictos territoriales con otras comunidades (Bravo-Pazmino, 1998; INVEMAR et al., 2006; Alcaldia Municipal Mosquera, 2008; Alcaldia Municipal Francisco Pizarro, 2012).

La perdida de navegabilidad del rio Patia en el sector sur, ha dejado a la poblacion practicamente aislada e imposibilitada para transportar sus productos, lo que ha obligado a muchos establecimientos comerciales a cerrar o a reubicarse en otros poblados. Para el caso del municipio Francisco Pizarro, el 40% de los poblados en zona rural tienen problemas de acceso con mareas bajas y en epocas de fuerte verano, asi mismo, el acceso maritimo a la cabecera municipal, Salahonda, solo puede realizarse con marea alta (Alcaldia Municipal Francisco Pizarro, 2012). Las dificultades en el transporte fluvial han contribuido a que gran parte de las actividades economicas de los poblados de este sector, como Salahonda y Mosquera, se trasladen a la zona norte, como es el caso de la actividad maderera en Bocas de Satinga (Defensoria del Pueblo, 2009).

CONCLUSIONES

La reactivacion del lobulo norte y la perdida de caudal en el sector sur del delta del rio Patia, han derivado en importantes afectaciones sobre los complejos sistemas estuarinos del area, asi como en cambios morfologicos en la zona costera y la red fluvial del delta. Los procesos erosivos y de sedimentacion, sumados a la alteracion de la productividad ecosistemica y los cambios en la biodiversidad de fauna y flora, se ven agravados por las condiciones de desarrollo socioeconomico de la poblacion, dando como resultado una problematica ambiental y social de grandes dimensiones.

Las tasas de acrecion y erosion presentadas en este articulo, evidencian la acelerada transformacion de la linea de costa del delta y la extension y magnitud de los procesos constructivos y destructivos que se registran en los lobulos norte y sur, respectivamente. Estos resultados demuestran el peso que las actividades humanas pueden tener dentro de los procesos evolutivos de las zonas costeras.

Por otro lado, los cambios ecologicos en los ecosistemas de manglar abarcan no solo la perdida de extension de la cobertura vegetal, sino que tambien involucran la perdida de la biodiversidad floristica debido a la colonizacion de vegetacion de agua dulce y la dominancia de especies de manglar menos vulnerables a las nuevas condiciones de salinidad que se presentan en el lobulo norte del delta. Los impactos en los ecosistemas estuarinos y de manglar repercuten sobre los servicios ambientales que estos prestan a la poblacion, lo que deriva en una mayor explotacion de los recursos naturales por parte de las comunidades costeras debido a la perdida de productividad de estos ecosistemas.

DOI: 103856/vol42-issue1-fulltext-4

AGRADECIMIENTOS

El desarrollo de este estudio fue posible gracias al apoyo del Instituto Nacional para el Desarrollo de la Ciencia y la Tecnologia "Francisco Jose de Caldas", contrato 1216-452-21267 (Morfodinamica del delta del rio Patia, Costa Pacifica Colombiana), el Departamento de Geologia y el Laboratorio de Sistemas de Informacion Geografica de la Universidad EAFIT, la Direccion General Maritima DIMAR, a traves de su Centro de Investigaciones Oceanograficas e Hidrograficas del Pacifico (CCCP) y la Unidad Administrativa Especial del Sistema de Parques Nacionales Naturales (UAESPNN).

REFERENCIAS

Alcaldia Municipal Francisco Pizarro. 2012. Plan de Desarrollo Municipal 2012-2015 "unidos por el bienestar de todos". Municipio de Francisco Pizarro, Salahonda, 222 pp.

Alcaldia Municipal Olaya Herrera. 2012. Proyecto de acuerdo Plan de Desarrollo 2012-2015 "De la Mano con el Pueblo". Municipio de Olaya Herrera, Bocas de Satinga, 22 pp.

Alcaldia Municipal Mosquera. 2008. Plan de Desarrollo Municipal 2008-2011 "Mosquera Positiva". Municipio de Mosquera, Mosquera, 61 pp.

Armenakis, C., F. Leduc, I. Cyr, F. Savapol & F. Cavayas. 2003. A comparative analysis of scanned maps and imagery for mapping applications. ISPRS J. Photogramm., 57(5-6): 304-314.

Bravo-Pazmino, H. 1998. Diversidad cultural y manglares del Pacifico colombiano. Ministerio del Medio Ambiente, Bogota D.C., 360 pp.

Castillo, R. 1995. El canal Naranjo: historia de una tragedia socio-ambiental en la cuenca baja del rio Patia. Universidad Nacional Sede Medellin, Medellin, 86 pp.

Coleman, J., H. Roberts & G. Stone. 1998. Mississippi River Delta: an overview. J. Coastal Res., 14(3): 698-716.

Conn, B.J. 1983. Aquatic and semi-aquatic flora of the Purari River system. In: T. Petr (ed.). The Purari-Tropical environment of a high rainfall river Basin. Dr. W. Junk Publishers, The Hague, pp. 283-293.

Dahdouh-Guebas, F., S. Hettiarachchi, D. Lo Seen, O. Batelaan, S. Sooriyarachchi, L.P. Jayatissa & N. Koedam. 2005. Transitions in ancient inland freshwater resource management in Sri Lanka affect biota and human populations in and around coastal lagoons. Curr. Biol., 15(6): 579-586.

Defensoria del Pueblo. 2009. Informe Defensorial Canal Naranjo: impactos y situacion actual. Defensoria Delegada para los Derechos Colectivos y del Ambiente, Bogota D.C., 39 pp.

Del Valle-Arango, J.I. 1994. Anotaciones sobre el clima de los bosques de guandal del delta del rio Patia. Rev. Fac. Nal. Agr., 47(1-2): 145-159.

Delgado, M.F., W. Gualteros, S. Espinosa, C. Lucero, A.M. Roldan, L.A. Zapata, J.R. Cantera, C. Candelo, C. Palacio, O. Munoz, G. Mayor & D.L. Gil-Agudelo. 2010. Pianguando-estrategias para el manejo de la piangua (Cartilla). Serie de Publicaciones Generales INVEMAR 45, Cali, 20 pp.

Departamento Nacional de Planeacion (DNP). 2010. Documento CONPES 3660. Politica para promover la igualdad de oportunidades para la poblacion negra, afro colombiana, palenquera y raizal. Departamento Nacional de Planeacion, Bogota D.C., 148 pp.

ESRI. 2009. World shaded relief, [http://www.esri.com/ software/arcgis/arcgisonline/standard-maps.html]. Reviewed: 6 October 2010.

Ericsson, J.P., C.J. Vorosmarty, L. Dingman, L.G. Ward & M. Meybeck. 2006. Effective sea-level rise and deltas: Causes of change and human dimension implications. Global Planet. Change, 50(1-2): 63-82.

Eslava, J. 1993. Climatologia. In: P. Leyva (ed.). Colombia Pacifico Vol. 1. Fondo para la proteccion del medio ambiente "Jose Celestino Mutis" & Financiera Electrica Nacional, Bogota D.C., pp. 136-147.

Fondo Mundial para la Naturaleza (WWF). 2005. Diagnostico socioeconomico del sector pianguero y estado actual del recurso hidrobiologico piangua Anadara tuberculosa en la costa Pacifica de Narino. Fondo Mundial para la Naturaleza Colombia, Cali, 35 pp.

Fundepublico. 1994. Diagnostico preliminar de evaluacion de los impactos socioeconomicos por la elaboracion del canal Naranjo. Fundacion para la Defensa del Interes Publico, Botoga, 39 pp.

Gonzalez, J., I.D. Correa & O. Aristizabal. 2002. Evidencias de subsidencia cosismica en el delta del San Juan. In: I.D. Correa & J.D. Restrepo (eds.). Geologia y oceanografia del delta del rio San Juan. Litoral Pacifico colombiano. Fondo Editorial Universidad EAFIT, Medellin, pp. 91-110.

Gomez, H. 1986. Algunos aspectos neotectonicos hacia el suroeste del litoral Pacifico Surcolombiano. Rev. CIAF, 11(1-3): 281-298.

Hapke, C.J., D. Reid, B.M. Richmond, P. Ruggiero & J. List. 2006. National assessment of shoreline change. Part 3D: historical shoreline change and associated coastal land loss along sandy shorelines of the California Coast National. U.S. Geological Survey, 79 pp.

Herd, D.G., T.L. Youd, H. Meyer, J.L. Arango, W.J. Person & C. Mendoza. 1981. The great Tumaco, Colombia earthquake of December 1979. Science, 211: 441-445.

Instituto de Hidrologia, Meteorologia y Estudios Ambientales (IDEAM). 1996. Morfodinamica del delta del rio Patia. Instituto de Hidrologia, Meteorologia y Estudios Ambientales, Universidad Nacional de Colombia, Bogota D.C., 145.

Instituto de Hidrologia, Meteorologia y Estudios Ambientales (IDEAM). 1997. Morfodinamica, poblacion y amenazas naturales en la Costa Pacifica de Colombia. Instituto de Hidrologia, Meteorologia y Estudios Ambientales, Universidad Nacional de Colombia, Bogota D.C., 127 pp.

Instituto de Hidrologia, Meteorologia y Estudios Ambientales (IDEAM). 2010. Leyenda nacional de coberturas de la tierra. Metodologia CORINE Land Cover adaptada para Colombia Escala 1:100.000. Instituto de Hidrologia, Meteorologia y Estudios Ambientales, Bogota D.C., 72 pp.

Instituto Colombiano de Geologia y Mineria (INGEOMINAS) & Instituto Geografico Agustin Codazzi (IGAC). 2008. Investigacion Integral del Anden Pacifico Colombiano. Tomo 1. Geologia. Instituto Colombiano de Geologia y Mineria, Instituto Geografico Agustin Codazzi, Bogota D.C., 165 pp.

Instituto de Investigaciones Marinas y Costeras Jose Benito Vives de Andreis (INVEMAR), Corporacion Autonoma Regional del Cauca (CRC) & Corporacion Autonoma Regional de Narino (CORPONARINO). 2006. Unidad ambiental costera de la llanura aluvial del sur: caracterizacion, diagnostico y zonificacion ambiental. Instituto de Investigaciones Marinas y Costeras Jose Benito Vives de Andreis, Corporacion Autonoma Regional de Narino, Corporacion Autonoma Regional del Cauca, Santa Marta, 383 pp.

Jensen, J.R. 2005. Introductory digital image processing: a remote sensing perspective. Pearson Prentice Hall, New Jersey, 526 pp.

Lopez-Vizoso, J.M. 1986. La observacion de la tierra desde el espacio: el mapa de ocupacion del suelo de la Comunidad Economica Europea. Estud. Geogr., 50(196): 409-434.

Martinez, J.O., J.L. Gonzalez, O.H. Pilkey & W.J. Neal. 1995. Tropical Barrier Islands of Colombia's Pacific coast. J. Coastal Res., 11(2): 432-453.

Ministerio del Medio Ambiente (MINAMBIENTE). 1996. El ultimo gran cambio del rio Patia. Causas, efectos y soluciones. Ministerio del Medio Ambiente, Bogota D.C., 51 pp.

Moore, L. 2000. Shoreline mapping techniques. J. Coastal Res., 16(1): 111-124.

Moreno, J.M. 2003. Avulsion de canales en la evolucion morfodinamica del delta del rio Patia, Colombia. In: N.C. Castillo-Murillejo & D.N. Alvis-Palma (eds.). El mundo marino de Colombia: investigacion y desarrollo de territorios olvidados. Universidad Nacional de Colombia, Bogota D.C., pp. 121-137.

Morton, R.A. & F.M. Speed. 1998. Evaluation of shorelines and legal boundaries controlled by water levels on sandy beaches. J. Coastal Res., 14(4): 1373-1384.

Oficina de las Naciones Unidas Contra la Droga y el Delito (UNODC). 2011. Censo de cultivos de coca 2010. Oficina de las Naciones Unidas Contra la Droga y el Delito, Sistema Integrado de Monitoreo de Cultivos Ilicitos, Bogota D.C., 111 pp.

Oficina de las Naciones Unidas Contra la Droga y el Delito (UNODC) y Direccion Nacional de

Estupefacientes (DNE). 2011. Persistenacia y productividad de la coca en la Region Pacifico 2009-2010. Oficina de las Naciones Unidas Contra la Droga y el Delito, Direccion Nacional de Estupefacientes, Bogota D.C., 97 pp.

Overeem, I. & J.P.M. Syvitski. 2009. Dynamics and vulnerability of Delta Systems. LOICZ Report and Studies No. 35. GKSS Research Center, Geesthacht, 54 pp.

Pajak, M.J. & S. Leatherman. 2002. The high water line as shoreline indicator. J. Coastal Res., 18(2): 329-337.

Restrepo, J.D. 2008. Vulnerabilidad de los deltas ante el cambio global: contexto mundial y Colombia. In: J.D. Restrepo (ed.). Deltas de Colombia: morfodinamica y vulnerabilidad ante el cambio global. Fondo Editorial Universidad EAFIT, Medellin, pp. 241-276.

Restrepo, J.D. & B. Kjerfve. 2002. San Juan River delta, Colombia: tides, circulation, and salt dispersion. Cont. Shelf. Res., 22: 1249-1267.

Restrepo, J.D. & S. Lopez. 2008a. Morfodinamica de los deltas de Colombia. In: J.D. Restrepo (ed.). Deltas de Colombia: morfodinamica y vulnerabilidad ante el cambio global. Fondo Editorial Universidad EAFIT, Medellin, pp. 75-100.

Restrepo, J.D. & S. Lopez. 2008b. Morphodynamics of the Pacific and Caribbean deltas of Colombia, South America. J. S. Am. Earth Sci., 25(1): 1-21.

Restrepo, J.D. & J. Cantera. 2011. Discharge diversion in the Patia River delta, the Colombian Pacific: geomorphic and ecological consequences for mangrove ecosystems. J. S. Am. Earth Sci., 46: 183-198.

Restrepo, J.D., P. Zapata, J.M. Diaz, J. Garzon-Ferreira & C.B. Garcia. 2006. Fluvial fluxes into the Caribbean Sea and their impact on coastal ecosystems: the Magdalena River, Colombia. Global Planet. Change, 50: 33-49.

Sanchez-Arcilla, A., J. Jimenez & H. Valdemoro. 1998. The Ebro Delta: morphodynamics and vulnerability. J. Coastal Res., 14(3): 745-772.

Stanley, D.J. & A. Warne. 1998. Nile delta in its destruction phase. J. Coastal Res., 14(3): 794-825.

Sistema de informacion geografica para el ordenamiento territorial nacional (SIGOT). 2011. Sistema de informacion geografica para el ordenamiento territorial nacional. Instituto Geografico Agustin Codazzi (IGAC). Bogota D.C. [http://sigotn.igac. gov.co]. Reviewed: 20 February 2011.

Syvitski, J.P.M. & Y. Saito. 2007. Morphodynamics of deltas under the influence of humans. Global Planet. Change, 57(3-4): 261-282.

Syvitski, J.P.M., C.J.S. Vorosmarty, A.J. Kettner & P. Green. 2005. Impact of humans on the flux of terrestrial sediment to the global coastal ocean. Science, 308: 376-380.

Tavera, H. 2009. Proyecto caracterizacion, diagnostico y zonificacion de los manglares en el Departamento de Narino, Documentos sintesis. Ministerio de Ambiente, Vivienda y Desarrollo Territorial, Corporacion Autonoma Regional de Narino, Fondo Mundial para la Naturaleza Colombia, Bogota D.C., 50 pp.

Thieler, E.R., E.A. Himmelstoss, J.L. Zichichi & A. Ergul. 2009. Digital shoreline analysis system (DSAS) version 4.0 An ArcGIS extension for calculating shoreline change. U.S. Geological Survey Open-File Report 2008-1278. [http://pubs.usgs.gov/ of/2008/1278/]. Reviewed: 10 January 2010.

Thom, B.G. 1967. Mangrove ecology and deltaic geomorphology, Tabasco, Mexico. J. Ecol., 55: 301-343.

Unidad Administrativa Especial del Sistema de Parques Nacionales Naturales (UAESPNN). 2005. Plan de manejo del Parque Nacional Natural Sanquianga. 2005-2009. Unidad Administrativa Especial del Sistema de Parques Nacionales Naturales, Cali, 294 pp.

Valiela, I., J. Bowen, & J. York. 2001. Mangrove forests: on of the world's threatened major tropical environments. BioScience, 50(10): 807-815.

Vann, J. 1959. Landform-vegetation relationships in the Atrato Delta. Ann. Ass. Am. Geogr., 49: 345-360.

Velez, M.I., M. Wille, H. Hooghiemstra, S. Metcalfe, J. Vanderberghe, & K. Van der Borg. 2001. Late Holocene environmental history of the southern Choco region, Pacific Colombia; sediment, diatom and pollen analysis of core El Caimito. Paleogeogr. Paleoclimatol. Paleoecol., 173: 197-214.

Vorosmarty, C.J.S. & D. Sahagian. 2000. Anthropogenic disturbance of the terrestrial water cycle. BioScience, 50(9): 753-765.

West, R.C. 1957. The Pacific lowlands of Colombia--a Negroid area of the American tropics. Louisiana University Press, Baton Rouge, 278 pp.

Wolanski, E. & S. Spagnol. 2000. Environmental degradation by mud in tropical estuaries. Reg. Environ. Change, 1(3-4): 152-162.

Yang, X., M.C.J. Damen & R.A. Zuidam. 1999. Satellite remote sensing and GIS for the analysis of channel migration changes in the active Yellow River Delta, China. Int. J. Appl. Earth Obs., 1(2): 146-157.

Received: 10 September 2012; Accepted: 27 November 2013

Adriana S. Parra (1) & Juan D. Restrepo Angel (2)

(1) Facultad de Ciencias del Mar, Universidad Autonoma de Sinaloa Programa de Maestria en Ciencias de la Tierra, Universidad Escuela de Administracion Finanzas e Instituto Tecnologico (EAFIT), A.A. 3300, Medellin, Colombia

(2) Departamento de Geologia, Escuela de Administracion, Finanzas e Instituto Tecnologico (EAFIT) A.A. 3300, Medellin, Colombia

Corresponding author: Adriana S. Parra (aparrar@eafit.edu.co)
Tabla 1. Tasas de cambio promedio de erosion y acrecion
litorales para los lobulos norte y sur del delta del rio
Patia, durante los periodos analizados.

                             1986-1998    1998-2001    2001-2008

Lobulo norte

Tasa de cambio                   -4,76        -3,02        +1,78
  promedio (m/ano)
Distancia de cambio             -57,59        -7,17       +13,39
  promedio (m)
Transectos erosion (%)          66,94%       30,90%       22,43%
Transectos acrecion (%)         26,91%       10,47%       43,52%
Transectos equilibrio (%)        6,15%       58,64%       34,05%
Numero total de transectos         602

Lobulo sur

Tasa de cambio                   -2.21        +4,81        -2,31
  promedio (m/ano)
Distancia de cambio             -26,75       +13,59       -16,18
  promedio (m)
Transectos erosion (%)          60,15%       10,83%       52,22%
Transectos acrecion (%)         35,20%       26,31%        6,38%
Transectos equilibrio (%)        4,64%       62,86%       41,39%
Numero total de transectos         517

Tabla 2. Medidas de error estimadas para las lineas de
costa del delta del rio Patia durante los anos analizados.
[E.sub.g]: Error de georreferenciacion, [E.sub.d]: Error de
digitalizacion, [E.sub.lc]: Error total linea de costa.

Fecha                    Medidas de error (m)

            [E.sub.g]    [E.sub.g]   [E.sub.d]   Error total
            zona norte   zona sur                ([E.sub.lc])

1986            10          8,8                     26,16
1997/1998     12.31         9,0        22.5         27,18
2001           7.5          8,6                     25,21
2007/2008       --          --                      22,50

Tabla 3. Comparacion de areas de cobertura vegetal en
1986 y 2001 en tres zonas del Parque Nacional Natural
Sanquianga (PNNS), delta del rio Patia.

Ano                         1986   2001

Recuadro A

Manglar denso/alto (ha)     2820   3892
Manglar abierto/bajo (ha)   3973   1891
Vegetacion mixta (ha)        727   1500
Total                       7520   7283

Recuadro B

Manglar denso/alto (ha)     2137   1842
Manglar abierto/bajo (ha)   3270   3296
Vegetacion mixta (ha)        219    624
Total                       5626   5761

Recuadro C

Manglar denso/alto (ha)     1077    761
Manglar abierto/bajo (ha)   1496    993
Vegetacion mixta (ha)       2299   3269
Total                       4871   5022
COPYRIGHT 2014 Pontificia Universidad Catolica de Valparaiso, Escuela de Ciencias del Mar
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2014 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Author:Parra, Adriana S.; Angel, Juan D. Restrepo
Publication:Latin American Journal of Aquatic Research
Date:Mar 1, 2014
Words:12900
Previous Article:Hematologia de juveniles de pacu, Piaractus mesopotamicus (Holmberg, 1887) alimentado con raciones suplidas con manan oligosacaridos (MOS).
Next Article:Distribucion y fuentes de fitoesteroles en sedimentos costeros y de rios del centro-sur de Chile.
Topics:

Terms of use | Privacy policy | Copyright © 2021 Farlex, Inc. | Feedback | For webmasters |