OFERTA HÍDRICA SUPERFICIAL DEL HUMEDAL EL GALLINAZO, UBICADO EN AGUACHICA-COLOMBIA

SUPPLY OF SURFACE WATER IN “EL GALLINAZO” WETLAND, LOCATED ON AGUACHICA-COLOMBIA

Revista de Investigación Agraria y Ambiental

Universidad Nacional Abierta y a Distancia, Colombia

ISSN: 2145-6097

ISSN-e: 2145-6453

Periodicidad: Semestral

vol. 14, núm. 1, 2023

riaa@unad.edu.co

Recepción: 04 Julio 2022

Aprobación: 18 Septiembre 2022

Publicación: 18 Diciembre 2022

URL: http://portal.amelica.org/ameli/journal/130/1303815011/

DOI: https://doi.org/10.22490/21456453.5959

Financiamiento
Fuente: Universidad Popular del Cesar [seccional Aguachica]
Nº de contrato: Evaluación del uso de los recursos ecosistémicos del humedal el Gallinazo en Aguachica–Cesar”, con el código de aval “CIPSA2019-08” otorgado por la resolución institucional 3397 del 2019
Beneficiario: OFERTA HÍDRICA SUPERFICIAL DEL HUMEDAL EL GALLINAZO, UBICADO EN AGUACHICA-COLOMBIA

Esta obra está bajo una Licencia Creative Commons Atribución-NoComercial-CompartirIgual 4.0 Internacional.

CÓMO CITAR: Saldaña-Escorcia, R., Otalvarez Herrera, R. y Herrera Galviz, J.D. (2023). Oferta hídrica superficial del humedal El Gallinazo, ubicado en Aguachica - Colombia. Revista de Investigación Agraria y Ambiental, 14(1), 221 - 249. https://doi.org/10.22490/21456453.5959

Los humedales son ecosistemas con suelos saturados, ya sea de agua superficial o subterránea, con vegetaciones adaptadas a ese tipo de suelos. Estos ecosistemas se distinguen por ser uno de los más beneficiosos para las personas, puesto que proporcionan múltiples servicios como mitigación de inundaciones, amortiguamiento contra tormentas, tratamiento de aguas, regulación de la erosión, retención de nutrientes y contaminantes, así como la recarga de acuíferos (Chatanga et al., 2020; Rebelo et al., 2019). Estos ecosistemas son objeto de presiones [agricultura, deforestación, ganadería, crecimiento urbano, entre otras] que a menudo son intensivas (Mehvar et al., 2019) y requieren de grandes cantidades de agua, lo que ocasiona una competencia por la obtención del recurso hídrico. Asimismo, el cambio climático afecta aún más la disponibilidad del agua (Fay et al., 2016; Zhu et al., 2017) ya que hay una evapotranspiración de alto potencial (Přibáň & Ondok, 1985; Watson et al., 2018), lo que exacerba dicha competencia.

El recurso hídrico es fundamental para el equilibrio de los ecosistemas, al igual que para el desarrollo socioeconómico de las poblaciones (Bal et al., 2021) Sin embargo, en los últimos años, la expansión urbana y la polución hídrica, al igual que el cambio climático, han llevado a una disminución drástica de la disponibilidad del agua dulce a nivel mundial (Meraj et al., 2018). La crisis del agua sigue aumentando a niveles extremos, lo que ha causado que las zonas áridas y semiáridas, las cuales presentan dificultades por el acceso al recurso hídrico (Kumar et al., 2021), se estén volviendo comunes. Por ende, la cuantificación y preservación de los suministros de agua superficial son necesarias para consolidar la supervivencia de las poblaciones (Soulis, 2021).

Estas necesidades requieren estar respaldadas por un conocimiento adecuado de las escalas temporales y espaciales de variabilidad de los recursos hídricos para permitir la aplicación de estrategias apropiadas de gestión y desarrollo (Trigg & Tshimanga, 2020). Por ende, el conocimiento sobre los flujos de agua superficiales es fundamental para la predicción de inundaciones, transporte de sedimentos y contaminantes, así como para el diseño de sistemas de riego agrícolas e infiltración del agua en el suelo. No obstante, los estudios con descripciones detalladas sobre los flujos de ríos son escasos, debido al límite temporal y recursos disponibles. Por lo cual, los modelos matemáticos relacionados con los flujos de agua presentan una limitada aplicabilidad a causa de la disponibilidad adecuada de datos requeridos. Esto exige formular y/o establecer nuevos métodos simples para predecir dichos movimientos (Govindaraju et al., 1992; Negi et al., 2021).

Los modelos hidrológicos empleados [comúnmente] para abordar las necesidades de cuantificación y preservación del agua emplean datos y generan información confiable que ayuda a respaldar los dictámenes relacionados con la gestión integral del recurso (Hughes et al., 2020; Tshimanga et al., 2011); en el contexto del cambio climático, los modelos también proporcionan predicciones para la planificación futura (Hughes & Farinosi, 2020; Meng et al., 2016; Nonki et al., 2019). En la actualidad, existen varios modelos hidrológicos de lluvia-escorrentía que se diferencian por los detalles en los procesos, es decir, la conceptualización de los procesos, los datos requeridos de entrada y la resolución espacial-temporal (Nonki et al., 2021). Generalmente se debe elegir entre diversos modelos: continuos, agrupados, distribuidos, deterministas o estocásticos, dependiendo de los datos de entradas que requieren y se tengan para ejecutarlos Sin embargo, los modelos lluvia-escorrentía tienen un mayor uso, debido principalmente a los pocos datos de entrada que requieren, lo cual los hace útiles en zonas con escasez de información (Tegegne et al., 2017). Estos modelos también son simples y fáciles de implementar y Autores como Skaugen et al. (2015), Muñoz et al. (2016) y Parra et al. (2018) demostraron que estos modelos conceptuales son alternativas efectivas y/o confiables en la predicción hidrológica con información limitada.

Los modelos lluvia-escorrentía son un campo clave de las investigaciones hidrológicas, ya q que proporcionan información fundamental para la gestión del recurso hídrico (Tajiki et al., 2020; Xie et al., 2020; Xie, Liu, Zhang, Han, et al., 2021), principalmente en cuanto a desastres por sequías (Kang & Sridhar, 2017; Pan et al., 2020) e inundaciones (Gao et al., 2019; Hostache et al., 2018; X. Zhang et al., 2018). Del mismo modo, los modelos de lluvia-escorrentía convencionales pueden estimar los caudales a través de procesos hidrofísicos (Xie, Liu, Zhang, Wang, et al., 2021). Por otra parte, los subproductos que arroja el modelo están basados en mecanismos físicos no lineales, entre ellos, la evapotranspiración, la infiltración, la humedad del suelo, el movimiento de las aguas subterráneas y las actividades antrópicas como el uso de los suelos (Xie, Liu, Zhang, Han, et al., 2021). No obstante, estimar estos procesos en los modelos hidrológicos es complejo debido a los insumos y a la información disponible.

La relación lluvia-escorrentía ha sido una motivación fundamental para el desarrollo de procesos hidrológicos que utilicen cuencas como fuente de agua, puesto que la lluvia termina depositándose en cuencas. Esta relación entre la lluvia y la escorrentía es influenciada, principalmente, por variables hidrológicas y fisiográficas del área, tales como la pendiente, cobertura vegetal, entre otras (Gabriels et al., 2021). No obstante, las actividades humanas han alterado estas variables; es decir, los ciclos de precipitaciones están siendo erráticos debido al cambio climático antropogénico (Pachauri et al., 2014). Así mismo, el uso de los suelos circundantes ha cambiado, debido al crecimiento demográfico (urbanización), afectando las características fisiográficas de la zona (Braud et al., 2013; Yan et al., 2013) y, por ende, su capacidad de retener el agua.

Teniendo en cuenta lo anterior, para la valoración de la influencia de las actividades humanas en las variables que afectan la relación lluvia-escorrentía, se emplean modelos basados en el número de curvas [SCS-CN]. Este método, establecido por Servicio Nacional de Conservación de Recursos [NRSC] y el Departamento de Agricultura de los Estados Unidos [USDA] en 1969, es un modelo simple y conceptual empleado para estimar la escorrentía directa en una zona determinada y se fundamenta en la profundidad de las precipitaciones (Shi & Wang, 2020). Este método es ampliamente utilizado, debido a su facilidad de implementación, porque el número de curvas (CN) es un parámetro accesible mediante tablas ya establecida por otros autores, que describe los efectos adversos de las actividades humanas en los usos de los suelos y por ende en la escorrentía (American Society of Civil Engineers, 2009).

Ahora bien, esta investigación busca estimar la oferta hídrica en el humedal El Gallinazo, usando un modelo lluvia-escorrentía. Este humedal se encuentra en la periferia sureste de Aguachica, Cesar, la cual es un área proyectada para la expansión del municipio, y esto hace que el ecosistema se encuentre sometido a presiones antrópicas por el incremento en áreas deforestadas y el desarrollo de actividades agropecuarias, lo que genera contaminación de los cuerpos de aguas y compactación y desnivelación de los suelos, teniendo efecto directo en los servicios ecosistémicos, que, en este caso, sería la oferta hídrica. Por eso, esta investigación posee como objetivo principal estimar la oferta hídrica neta disponible en el humedal, puesto que no existe información básica para una gobernanza del agua eficaz e incluyente en la zona. Por el contrario, tener información sobre la oferta hídrica del humedal contribuiría en la formulación de políticas públicas sobre la gestión integral del recurso hídrico, logrando maximizar el bienestar social y económico en el municipio de Aguachica.

2.1 Área de estudio

Esta pesquisa se llevó a cabo en el humedal El Gallinazo [Figura 1], el cual se encuentra ubicado en la periferia sureste a 260 metros del casco urbano del municipio de Aguachica-Cesar (Saldaña-Escorcia & Otalvarez-Herrera, 2021). El municipio se encuentra en las siguientes coordenadas: latitud 8°17’58.124” Norte con longitud 73°36’24.825” Oeste, con un área de 58823 m² (5.88 hectáreas) y un perímetro de 1322 m. Dicho humedal pertenece a la microcuenca el Gallinazo y esta, a su vez, a la cuenca de la quebrada Buturama y quebrada Guaduas, que desemboca en el río Magdalena.

Figura 1.
Zona de estudio.
Autores.

2.2 Método SCS-CN

Este método presenta un enfoque simple, empleado para estimar la escorrentía superficial directa dentro de una cuenca hidrográfica y/o fluvial (Kumar & Jhariya, 2017; Shi & Wang, 2020; Singh et al., 2015; Verma et al., 2020). El primer paso de este modelo es el cálculo del área del drenaje de la cuenca del afluente (Kumar et al., 2021) y la ecuación del modelo tiene una relación marcada entre el volumen de la escorrentía superficial directa y la precipitación, como se expresa a continuación:

(1)

Q = profundidad de la escorrentía (mm)

P = profundidad de la lluvia (mm)

I_a = abstracción inicial (mm)

S = retención máxima potencial (mm)

Teniendo en cuenta que el volumen de la precipitación que cae antes de que dé comienzo la escorrentía se establece como I_a y es asuminda como 0,2*S (Soulis & Valiantzas, 2021), según Ling et al. (2019), la ecuación generalmente se escribe como:

(2)

Condicionando que la S está relacionada con el CN y que de esta se deriva el I_a a través de la multiplicación con el parámetro λ (-) se obtiene que:

(3)

El método SCS-CN asume que la variable λ toma valores de 0.2 y 0.8 (Shi et al., 2009) y emplea valores de CN en los que este es un parámetro adimensional [Rango de 0-100] que se determina mediante tablas basadas en el uso que se les da a las tierras y/o la cobertura vegetal que estas posean. En este caso se empleó, para definir las categorías de cobertura vegetal, la metodología Corine Land Cover adaptada para Colombia, la cual consistió en el uso de imágenes Landsat 8 [2019-2020] ortorrectificadas, combinadas y ajustadas bajo las firmas espectrales de las bandas con el programa ArcGIS 10.8; también se configuraron de acuerdo con la leyenda nacional de coberturas a escala 1: 100.000 adaptada por el Instituto de IDEAM.

Por otro lado, las tablas que se utilizaron están relacionadas con los Grupos Hidrológicos de Suelos [GHS] y estos, a su vez, con la tasa de infiltración del suelo dentro de la cuenca. Los Grupos Hidrológicos de Suelos se determinaron por medio de 57 muestras de suelo y se estableció la textura de cada una empleando la metodología establecida por el Centro Internacional de Agricultura Tropical [CIAT], en conjunto con la Corporación Biotec y la Universidad Nacional de Colombia (CIAT et al., 2006). Generalmente, los GHS se dividen en cuatro grupos [A, B, C y D] en los que el grupo de suelo hidrológico “A” posee una tasa de infiltración alta, el “B” una tasa moderada, el “C” una tasa moderada-baja y el grupo “D” tiene suelo con baja tasa de infiltración (Lian et al., 2020).

2.3 Precipitación media anual

2.3.1 Datos hidrometeorológicos

El estudio se desarrolló a partir de los datos obtenidos de 7 estaciones pluviométricas pertenecientes al IDEAM. La distribución espacial de estas estaciones de medición sobre y alrededor del área de estudio se muestra en la Figura 2.

Figura 2.
Posicionamiento espacial de las estaciones de medición a lo largo y alrededor del humedal El Gallinazo.
Autores.

2.4 Método de Kriging

Los métodos de interpolación espacial [basados en la geografía y en la estadística] se implementan hace décadas, utilizando datos de puntos específicos y áreas. Estas técnicas de interpolación espacial permiten obtener y/o predecir datos de otros sitios sobre la zona de estudio (Bechler et al., 2015; Camera et al., 2014; Manz et al., 2016). El método de Kriging y sus variantes se convirtieron en modelos de interpolación de referencia en el campo de la hidrología, puesto que autores como Das & Islam (2021), así como Wagner et al. (2012), indican que estos ofrecen resultados ajustados y con mayor precisión con relación a otros métodos como el polígono IDW y Thiessen.

El método de interpolación de Kriging se puede desempeñar de forma no lineal, debido a las coordenadas espaciales (Biau et al., 1999; Xu et al., 2018). Además, el método consta de dos partes: un análisis de variaciones espaciales a través de variogramas y un cálculo de pesos ponderados para las distancias interpoladas, con lo cual se crea un red adaptativa sistemáticamente que busca la mejor estimación de los resultados (Peng et al., 2016; Rousis et al., 2017; Wang & Lin, 2015). De manera general, la fórmula del modelo (1) se emplea para apreciar los valores de las ubicaciones no muestreadas $\hat{p}$($u_0$), en la que los pesos $w_i$($u_0$) están asociados con datos de estaciones calibradas cercanas a la zona de estudio w(u_i) y están determinadas así:

(4)

Los datos espaciales incorporados en la ecuación de Kriging son utilizados para estimar de forma robusta los pesos de los puntos desconocidos, calibrándose de tal forma que la interpolación es insesgada, mientras que la varianza del modelo disminuye (Das & Islam, 2021; Webster & Oliver, 2007). La ecuación (5) es adoptada para establecer el valor de los pesos en los puntos desconocidos de manera ordinaria:

(5)

$\text{γ}$(u_i,u_j) y $\mathrm{\gamma }$(u_j,u_o) correspondientemente, son las semivarianzas entre los puntos muestreados (u_i y u_j); los puntos muestreados y los no muestreados es (u_j y u_o); μ es la medida de Lagrange.

Por otra parte, en el método de Kriging con información externa, la media del parámetro primario depende de los coparámetros. Este método se utiliza en función a los residuos obtenidos de la estimación del valor medio y se evalúan las siguientes ecuaciones con el fin de establecer una ponderación para este procedimiento (Das & Islam, 2021; Webster & Oliver, 2007):

(6)

$\mathrm{\gamma }$(u_i,u_j) es las semivarianza residual entre los puntos muestreados; $\mathrm{\gamma }$(u_i,u_o) es la semivarianza entre los puntos muestreados y no muestreados; $m_k$ para k = 0,1,....,K, son los multiplicadores de Lagrange; $y_k$ para k = 1,2,....,K son las covariables. La valoración empírica del semivariograma fundada en agrupaciones de datos observados tiene la siguiente expresión:

(7)

$\hat{y}$ (h) = Semivarianza

h = Distancia de separación entre dos ubicaciones

N(h) = Número de pares de datos

(z(u_i) - z(u_i + h)) = Diferencia de valores entre puntos (u_i + h)

Para efectos del presente trabajo se empleó el software ArcGIS 10.8 que tiene incorporados algoritmos programados con diferentes métodos de interpolación, de los cuales Kriging es uno de ellos. Se realizó la ejecución del método utilizando la configuración de Kriging ordinario con un tipo de salida de superficie: predicción. Además, se configuró el tipo de transformación logarítmica con un orden de eliminación de tendencias secundario. El tipo de modelo fue esférico con la activación de la opción de anisotropía. Al ejecutar la herramienta se obtuvo un semivariograma con rango entre 0.118 y 0.237, una meseta de 1,0 y una pepita de 0,01035.

2.5 Precipitación efectiva

El concepto de precipitación efectiva [Pe] puede tener una interpretación diferente para cada campo científico o para diferentes objetivos de investigación. Por ejemplo, Pe puede considerarse, respectivamente, como la cantidad total de la precipitación que entra en un embalse, el porcentaje de la precipitación que contribuye al aumento en los niveles de las aguas subterráneas, el volumen de agua que utilizan los sistemas radiculares de las plantas, etc. (Tigkas et al., 2016). En este trabajo nos centramos en el concepto de Pe como factor que contribuye a la recarga de las aguas subterráneas, y se determina a través de la ecuación planteada por Magaña-Hernández et al. (2013):

(8)

P_e = Precipitación efectiva (mm)

P = Precipitación media anual (mm)

CN = Número de curvas (adimensional)

2.6 Oferta hídrica superficial neta

Para la obtención final de la oferta hídrica superficial [neta] disponible en la zona de estudio se aplicaron dos de los factores establecidos en la Resolución 865 de 2004 (Ministerio de Ambiente Vivienda y Desarrollo Territorial, 2004), por la cual se fija la metodología para hallar el índice de escasez para aguas superficiales:

· Reducción por calidad de agua: Corresponde a la probable variación de la calidad del agua debido a la manifestación de materia orgánica y sintética; de acuerdo con la resolución, la oferta hídrica total se debe multiplicar por el 25 %.

· Reducción por caudal ecológico: Este factor constituye el caudal mínimo demandado para el equilibrio del ecosistema, adoptando un 25 % sobre la oferta hídrica total.

2.7 Análisis Estadísticos

Este tipo de análisis únicamente se ejecutó para el cálculo de isoyetas, toda vez que la determinación de coberturas vegetales requirió la estimación de puntos de control en campo con los cuales se desarrolló un proceso de delimitación bastante eficiente, por tratarse de un área pequeña. Para la realización del mapa de isoyetas se empleó la herramienta de análisis geoestadísticos, en la cual se ejecutaron los análisis de tendencias, histograma, gráfico cuantil–cuantil [QQPlots], el gráfico de semivariograma y el de covarianza con los cuales se determinó el comportamiento de los datos. De acuerdo con ello, se establecieron las configuraciones establecidas para la interpolación mediante Kriging, del cual se obtuvieron las respectivas isoyetas, y que se describió anteriormente.

3.1 Coberturas del Suelo

Las coberturas vegetales se determinaron a través de reconocimiento mediante salidas de campo en el humedal y teniendo en cuenta la metodología Corine Land Cover para Colombia, con lo cual se determinaron dos tipos de coberturas vegetales mediante puntos de control: La primera, pastizales enmalezados, en su mayoría bajo suelos inundados o con niveles freáticos altos, ocupando aproximadamente el 27.29 %. La segunda, vegetación secundaria, se encontró principalmente alrededor del cauce del humedal El Gallinazo y en zonas secas, ocupando el 72.71 % del área del humedal [Figura 3].

Figura 3.
Coberturas del suelo del humedal El Gallinazo.
Autores.

La cobertura vegetal y los cambios en el uso de los suelos influyen directamente en la infiltración y en la escorrentía superficial, siendo estos procesos fundamentales para la generación de caudales, la retención de humedad en los suelos y la regulación hídrica dentro de la zona (Mozo et al., 2020) . Con relación a esto, las coberturas de vegetación secundaria baja, presentes en el humedal El Gallinazo, han disminuido durante el lapso de esta investigación debido a las actividades de extracción de material forestal ilegal. La disminución de la cobertura genera inconvenientes dentro del ciclo hidrológico regional, pues se alteran la variación y regulación del balance hídrico superficial y la evapotranspiración como lo mencionan Bao et al. (2015) ; Bégué et al. (2011) y Chen et al. (2014) . Igualmente, los autores Revueltas et al. (2020) y Zhang et al. (2019) indican que la cobertura vegetal es un factor fundamental en la estabilización de la temperatura, intercepción y redistribución de las lluvias, retención hídrica y escorrentía superficial.

3.2 Grupos Hidrológicos de Suelo

En la figura 4 se muestran los diferentes GHS pertenecientes al humedal El Gallinazo. El grupo A posee un 15.1 % de suelos arenosos; el grupo B, un 56.6 % de texturas moderadamente gruesas a moderadamente finas; mientras el grupo C, un 28.3 % exhibiendo texturas moderadamente finas.

Figura 4.
GHS del humedal El Gallinazo.
Autores.

La presencia de los diversos GHS dentro del humedal demuestra que los suelos poseen un potencial de recarga hídrica de alta a moderada, una capacidad de infiltración moderada y que pueden presentar un potencial de escorrentía moderadamente alto cuando presentan humedad, así como cuando se encuentran desprovistos de coberturas. Específicamente, el 28,3 % de los suelos del grupo C, con texturas finas e índices altos de humedad, poseen potencial moderadamente alto de escorrentía, y desprovistos de cobertura vegetal cambian los flujos de las corrientes superficiales, lo que aumenta la erosión (Duan et al., 2016; Liu et al., 2019) .

3.3 Número de Curva

En la tabla 1 se muestran los conjuntos de las coberturas vegetales y GHS obtenidos para el humedal El Gallinazo; con esto, a cada combinación establecida le corresponde un Número de Curva (CN), el cual representa el potencial de escorrentía superficial. Así los datos más bajos se hallaron en los suelos del grupo A con coberturas de vegetación secundaria y pastos enmalezados, al igual que en los suelos del grupo B con coberturas de vegetación secundaria baja.

Tabla 1.
CN para cada combinación suelo-cobertura.

Autores.

Kumar et al. (2021) afirman que la evaluación y/o determinación de la escorrentía superficial es fundamental para el establecimiento del de la oferta y la eficiencia hídrica de la zona, así como para la determinación de los niveles de riesgos en relación con las inundaciones aguas abajo; el método SCS-CN permitió establecer el parámetro de los niveles de riesgo de inundación con éxito, requiriendo de un menor empleo del tiempo con un mayor manejo de datos (Caletka et al., 2020) . Estos insumos servirán como base o línea base para la formulación de estrategias, planes o políticas de uso adecuado de los recursos hídricos, que en opinión de Zhou et al. (2021) es uno de los desafíos más importantes.

3.4 Precipitación Media Anual

Con los datos de las estaciones se obtuvo la precipitación media anual de la zona de estudio a través de las isoyetas mediante el método de Kriging [Figura 5], se analizaron los intervalos de precipitación y la influencia de estos en las áreas de la zona de estudio dando como resultado un valor aproximado de 1687.13 mm.

Figura 5.
Isoyetas.
Autores.

3.5 Determinación de la oferta hídrica

Con la determinación de las coberturas [Tabla 1] se estableció el número de curvas promedio para el área de estudio (CN) en 59,29. Con este valor se procedió al cálculo de la relación lluvia-escorrentía para determinar el almacenamiento máximo mensual de la zona de estudio:

El valor S representa la disimilitud entre la precipitación efectiva y la escorrentía superficial directa del periodo de estudio. Para esta ocasión, el valor de S alcanzado fue de 174,40 mm. Antes de hallar la escorrentía directa se estableció la precipitación efectiva, teniendo en cuenta la precipitación media anual y el número de curvas de la zona de estudio, de la siguiente manera:

Obtenidos dichos valores, se procede a determinar la escorrentía directa por medio de la ecuación 2, como se muestra a continuación:

Para el humedal El Gallinazo, el valor 𝑄 fue de 1467,52 mm al año. Este dato de escorrentía directa se convirtió a unidades de caudal:

Obteniendo este valor, se procede a multiplicarlo por el área total del humedal (58822 m²) para obtener un volumen, dando como resultado 0,00271 m³ por segundo. Aplicando los factores de corrección de la resolución descrita en la metodología, la oferta hídrica superficial [neta] disponible en el humedal El Gallinazo es de 0,00135 m³ por segundo o 1,35 litros por segundo, con un volumen anual de 42573,6 m³.

La disponibilidad del agua en el humedal El Gallinazo está influenciada por dos factores: la precipitación relacionada directamente con las temporadas de lluvia, lo que permite el flujo constante del agua a través de los drenajes, y las actividades ganaderas [bovina y equina] que generan afectaciones en los niveles del terreno, causando compactación del suelo y disminución de la profundidad de los drenajes, reduciendo el caudal efluente.

Al pasar los años, el crecimiento demográfico y la expansión de la frontera agropecuaria e industrial han aumentado la demanda hídrica, que en muchas ocasiones es mayor a la oferta hídrica neta de las cuencas (Brown et al., 2005) . Por lo tanto, si dicha presión antrópica persiste o aumenta en el humedal, puede tener efectos directos en la disponibilidad hídrica por los cambios en el uso de los suelos, como lo indican autores como Daza et al. (2014) y Ruales y Dávila (2019), quienes afirman que la labranza, la expansión de la frontera agropecuaria [eliminación de cobertura vegetal nativa por cultivos agrícolas y ganaderos] y la expansión demográfica disminuye entre un 40–60 % la disponibilidad hídrica.

Lo anterior sucede puesto que dichas actividades humanas alteran propiedades físicas de los suelos tales como la porosidad [macroporos y microporos], la densidad aparente, la retención de humedad y el contenido de carbono orgánico, lo que afecta la capacidad de almacenamiento de agua y, por ende, aumenta los problemas de erosión debido a la escorrentía que el suelo no puede absorber; conllevando, así, a problemas de degradación de los suelos, así como de disponibilidad y eficiencia hídrica. Asimismo, la calidad del agua se encuentra en riesgo y/o vulnerabilidad debido a la expansión urbana y agropecuaria hacia la zona periférica donde se ubica el humedal, que tiene como consecuencia la disminución del recurso hídrico con base en la disponibilidad [cantidad y calidad]. De acuerdo con Chávez-Jiménez y González-Zeas (2015) , esta escasez de agua genera conflictos socioeconómicos, al igual que altera los procesos hidrológicos desequilibrando los ecosistemas.

Por otra parte, al comparar los datos obtenidos en esta investigación sobre la oferta hídrica superficial [neta] disponible en el humedal El gallinazo [0,00135 m³ por segundo] con estudios similares, los cuales emplearon la metodología establecida en la Resolución 865 de 2004 (Ministerio de Ambiente Vivienda y Desarrollo Territorial, 2004), encontramos autores como Gutiérrez-Forero et al. (2015), que obtuvieron una oferta hídrica disponible entre un rango de 10,3 a 2,7 m³ por segundo, cuyos valores dependen de la estación climática de la zona. Salazar (2016) obtuvo en su investigación una oferta hídrica neta de 0,891 m³ por segundo, mientras que Caicedo (2019) registró una oferta neta de 1,13 m³ por segundo y, por último, Córdoba (2021) determinó una oferta 48,437 m³ por segundo.

Finalmente, se debe tener en cuenta que dichos estudios se realizaron a niveles de escala espacial diferentes, es decir, son investigaciones cuyas áreas de estudio fueron cuencas hidrográficas en comparación al área del humedal. Además, aunque la oferta hídrica superficial neta disponible saliente del humedal “El Gallinazo” no se encuentre entre el rango de estos estudios debido a su área, sigue siendo un ecosistema fundamental para el municipio y la región en relación con la lucha contra el cambio climático, la escasez del recurso hídrico y el desarrollo socioeconómico, al igual que para el logro de las metas pactadas en los Objetivos de Desarrollo Sostenible [ODS].

El cálculo de la oferta hídrica superficial, a través de la metodología lluvia-escorrentía, permitió el diagnóstico del recurso hídrico del humedal El Gallinazo, el cual indica que el humedal se enfrentará con problemas de escasez hídrica en un futuro no muy lejano, debido a las tendencias de expansión de la frontera agropecuaria y el crecimiento demográfico en la zona periférica. Esto aumentará las presiones ya generadas a este ecosistema por los cambios en los usos de la tierra, la demanda de agua para actividades socioeconómicas y la nula legislación de protección ambiental regional de estos ecosistemas, lo que afecta, de manera directa, el equilibrio de este.

Sin importar que Colombia posea una oferta hídrica generosa, dicho recurso no estará presente de forma infinita ni distribuido homogéneamente en los ecosistemas a nivel nacional. Por esta razón, la gestión adecuada del recurso hídrico superficial y/o subterráneo es esencial. Además, se deben tener en cuenta las tendencias de consumo, generación y escenarios futuros que, mediante métodos como la relación lluvia-escorrentía, permitan determinar de la oferta, demanda y los usos del agua; esto haría posible avanzar hacia una gestión integral del recurso hídrico [mediante políticas públicas ambientales] cumpliendo lo pactado en los ODS. Asimismo, en este contexto de futura escasez de agua se puede contribuir mediante el fortalecimiento de las estrategias de protección de los ecosistemas y, por ende, de su oferta hídrica existente; de igual forma, la reducción de la demanda de agua dentro de la zona de influencia, por medio de programas de uso eficiente y sostenible del recurso, es necesaria para garantizar la protección integral del recurso hídrico.

Por otra parte, la inexistencia de estudios o informes de investigación recientes sobre la calidad y demanda del recurso en la zona del humedal El Gallinazo, aparte de esta investigación, no permite concluir el estado del recurso hídrico en esa zona totalmente, por lo que se recomienda establecer y ejecutar otros estudios de investigación que permitan complementar la evaluación y/o determinación de la oferta hídrica superficial.

CÓMO CITAR: Saldaña-Escorcia, R., Otalvarez Herrera, R. y Herrera Galviz, J.D. (2023). Oferta hídrica superficial del humedal El Gallinazo, ubicado en Aguachica - Colombia. Revista de Investigación Agraria y Ambiental, 14(1), 221 - 249. https://doi.org/10.22490/21456453.5959

CONTRIBUCIÓN DE LA AUTORÍA: Primer autor: metodología, investigación, análisis de datos, conceptualización, escritura [borrador original], revisión y edición. Segundo autor: metodología, investigación, conceptualización, análisis de datos, escritura [borrador original]. Tercer autor: metodología, conceptualización, revisión y edición.

CONFLICTO DE INTERESES: Los autores declaran no tener ningún conflicto de intereses.

https://hemeroteca.unad.edu.co/index.php/riaa/article/view/5959 (html)

https://hemeroteca.unad.edu.co/index.php/riaa/article/view/5959/5962 (pdf)

Los autores agradecen al Grupo de Investigación en Gestión Ambiental y Territorios Sostenibles [GE&TES] de la Universidad Popular del Cesar [seccional Aguachica] por permitirnos ejecutar el subproyecto “Evaluación de la recarga hídrica en relación a la capacidad de retención y almacenamiento de agua del humedal El Gallinazo en Aguachica–Cesar” [Trabajo de grado]. Este subproyecto se ejecutó dentro del marco del macroproyecto de investigación “Evaluación del uso de los recursos ecosistémicos del humedal el Gallinazo en Aguachica–Cesar”, con el código de aval “CIPSA2019-08” otorgado por la resolución institucional 3397 del 2019.