Contextualización: actualmente, se reconoce a nivel internacional que los embalses, además de alterar el natural ciclo del carbono, son fuertes emisores de gases de efecto invernadero (GEI) y, en la cuantificación de sus emisiones, prevalecen grandes incertidumbres sobre las que, solo hasta el momento, se ha comenzado a generar conocimiento, mediante mediciones directas en campo sobre las verdaderas magnitudes de las emisiones que tienen lugar en estos hidrosistemas.

Vacío de conocimiento: los altos costos del desarrollo de mediciones in situ, análisis de muestras y cuantificación de emisiones, han sido un impedimento mayor para la producción de información suficiente que conlleve a una comprensión global del comportamiento de estos hidrosistemas en términos de emisiones.

Propósito del estudio: con base en los resultados del Proyecto GEITRO (2012-2014), cuyo objetivo principal fue cuantificar las emisiones de GEI generados en el embalse de Riogrande II (Antioquia, Colombia) a partir de muestreos exhaustivos in situ, se abrió la posibilidad de explorar métodos alternativos que permitan la cuantificación de emisiones mediante la aplicación de modelos fundamentados en el uso de variables medibles basados en sistemas remotos como la teledetección. El objetivo principal de este estudio fue evaluar cómo el uso conjugado de datos, medidos desde el espacio, y modelos que relacionen dichos datos con las concentraciones superficiales de un gas en el agua, así como su coeficiente de velocidad de transferencia, se podrían usar para estimar de manera acertada los flujos difusivos del gas en el espacio y en el tiempo, reduciendo los esfuerzos de monitoreo en campo.

Metodología: los métodos estuvieron fundamentados en el análisis de la información proveniente de los datos de campo para el desarrollo y validación de los modelos, así como la deducción de temperaturas superficiales del agua medidas a partir sensores remotos para la estimación final de emisiones a partir de los modelos propuestos.

Resultados y conclusiones: como resultado se han desarrollado dos modelos: uno para la estimación de concentraciones superficiales de CO2 y otro para la estimación del coeficiente de velocidad de transferencia normalizado sobre número de Schmidt de 600 (K600). El resultado de la aplicación de estos dos modelos condujo a la estimación de flujos difusivos de CO2 del embalse. Los procesos de validación permitieron concluir que los modelos desarrollados logran reproducir, de manera más que aceptable (error relativo 15 ± 0,07% en promedio), los flujos difusivos de CO2 medidos en el embalse de RGII entre los años 2012 y 2014. Finalmente, a partir de estos resultados, se emplearon las temperaturas superficiales, medidas por MODIS-Terra, con el fin de estimar los flujos difusivos de CO2 sobre el embalse.

Los embalses son considerados como una fuente de emisiones de gases de efecto invernadero de origen antrópico a causa de su naturaleza. Su comportamiento, como ecosistema acuático, responde a los mismos procesos fisicoquímicos y a la degradación de la materia orgánica (MO) que un ecosistema acuático de origen natural. A pesar de que existen en la actualidad grandes incertidumbres con relación a la magnitud precisa de emisiones de estos sistemas, diversos estudios (Guérin y Abril, 2007; Guérin et al., 2008) sugieren que los embalses pueden contribuir en un 4% y 12% de las emisiones globales de origen antrópico de CO2 y CH4 respectivamente.

Hasta el momento, no hay evidencia científica concluyente respecto a la magnitud real de emisiones de GEI que tienen lugar en embalse, y es probable que las aproximaciones difieran de manera importante según la región geográfica. El Informe del IPCC (2011), junto con otra serie de estudios más recientes, como el de Raymond et al. (2013) y Werhli (2013), sostienen que existen muy pocas mediciones de emisiones de GEI en regiones tropicales a nivel de embalses y de ecosistemas acuáticos en general y que, por el momento, la poca información disponible resulta insuficiente para dar una conclusión relacionada con los aportes de emisiones de GEI a partir de dichas regiones.

Existen importantes estudios que aportan conocimiento en este dominio de la cuantificación de emisiones en ecosistemas acuáticos, los cuales están representados en aquellos desarrollados por Guérin y Abril (2007) y Guérin et al., (2007) sobre el embalse de Petit-Saut en Guyana francesa, Kemenes et al., 2007; Kemenes et al.,2011 en la represa de Balbina en Brasil, Chanudet et al., (2011) en los embalses de Nam Ngum y Nam Leuk, Fearnside (2016); Fearnside (2009) y Demarty y Bastien (2011) en diferentes embalses tropicales. Sin embargo, estas iniciativas son insuficientes para determinar con precisión la contribución de emisiones de GEI de embalses en el trópico.

Los elevados costos relacionados con el desarrollo de trabajos de campo periódicos, compra y mantenimiento de equipos de medición, tomas y análisis de muestras, entre otros, dificultan la ejecución el desarrollo de esta labor. En este sentido, se hace inminente el desarrollo de modelos que permitan estimar de forma confiable el conjunto de emisiones de estos ecosistemas a partir de información que sea accesible por medio de sensores remotos (estaciones de medición en tiempo real, teledetección, etc.).

El uso de los sensores remotos para este fin ha sido altamente reportado para ecosistemas marinos (Stephens et al, 1995); sin embargo, su aplicación en ecosistemas acuáticos continentales está muy débilmente documentada. En esto se basa la necesidad de reportar los posibles resultados, positivos o negativos, de tal aplicación en un ecosistema acuático tropical de alta montaña.

Esta investigación tuvo como objetivo fundamental aplicar el uso de herramientas de observación, como sensores remotos, para alimentar modelos de emisiones de GEI actualmente reportados en la literatura, con el fin de reconstituir flujos difusivos de CO2 medidos en campo sobre el embalse de Riogrande II (RGII) y poder determinar la viabilidad del uso de este tipo de tecnologías en la estimación de tales emisiones de este biogás en ecosistemas acuático tropicales.

Entre el año 2012 y el 2014, el Grupo de Investigación en Ecología y Contaminación Acuática de la Universidad Nacional de Colombia, Sede Palmira, desarrolló el proyecto de investigación “Cuantificación de emisiones de GEI en embalses tropicales: Caso del Embalse de Riogrande II” – GEITRO (GEITRO 2012). Este proyecto tuvo como objetivo fundamental adquirir datos de campo durante 18 meses, para estimar las emisiones de CO2 provenientes del embalse de Riogrande II, así como el análisis del comportamiento biogeoquímico del sistema que permitiría comprender las relaciones que influyeron sobre tales emisiones. El conjunto de datos de campo y análisis realizados para el embalse fueron puestos a disposición para el desarrollo de este estudio.

El embalse de Riogrande II (Figura 1), está situado a 2.270 msnm en la región Andina de Colombia, al norte del departamento de Antioquia. Presenta un área inundada, a la cota de vertimiento, de aproximadamente 1.214 ha, y represa, a esa cota, un volumen de agua estimado en 236 mm3 con tiempos de retención hidráulicos de 77 días en promedio. Sobre la superficie, los anchos máximos del espejo de agua no superan los 900 m y pueden fluctuar en promedio entre los 300 m y 320 m, aproximadamente. Las profundidades máximas, que pueden alcanzar los 48 m, se encuentran en varias áreas de la presa y en el centro de confluencia de los dos principales tributarios que surten el sistema: Río Grande y Río Chico. Estos dos ríos, condicionados por regímenes de precipitación que varían entre 1900 y 2100 mm al año, pueden llegar a aportar, en promedio, alrededor de 15 a 25 m3 s-1, a partir de un área drenada de aproximadamente 1.050 km2. El sistema está globalmente condicionado por un régimen bimodal marcado por precipitaciones moderadas entre los meses de marzo, abril y mayo; y por precipitaciones fuertes durante los meses de septiembre y octubre. El período de sequía más importante tiene lugar hacia finales de diciembre e inicios de febrero.

Entremarzo 2012 y febrero 2014, se realizaron 14 trabajos de campo. Estas comisionescientíficas tuvieron como objetivo general tomar muestras en los diferentespuntos de control (Figura 1) para el seguimiento de la calidad del agua entérminos de parámetros fisicoquímicos superficiales como pH, oxígeno disuelto,temperatura y conductividad; concentraciones de gases en la columna de agua;medición de flujos difusivos desde la superficie del lago a la atmósfera; ymedición de flujos de ebullición. Con respecto a lo anterior, se usó en formapermanente una estación meteorológica portátil que registró durante todo esteperiodo, con intervalos de 30 minutos entre un registro y otro, las siguientesvariables: temperatura máxima y mínima, precipitación, evapotranspiración,humedad relativa, presión atmosférica, radiación solar, velocidad y direccióndel viento.

Enel marco de este estudio se consideran dos fuentes de información con relacióna los flujos difusivos de CO2 en el embalse de RiograndeII: i) Los flujos que ya se midieron en campo durante el proyecto GEITRO (2012)y, ii) Los flujos que se estimarán mediante lautilización conjunta de los modelos desarrollados en esta investigación eimágenes satelitales. Para el segundo caso, es necesario considerar que, eltotal de las emisiones de un gas (g) que tienen lugar desde el espejo de aguahacia la atmósfera son dependientes de: i) la variación de velocidades detransferencia del gas en dicha interface a una temperatura superficial del aguadada y, ii) la diferencia entre la concentración delgas en el agua y su concentración de equilibrio con la atmósfera según laecuación 1:

DondeF (CO2) es el flujo difusivo de CO2, kg,T es la velocidad de transferencia del gaspara un gas (g) particular a una temperatura (T) dada, y ΔC es el gradiente deconcentración de CO2 expresado como la diferencia entre laconcentración en el agua (Cw) y la concentración deequilibrio entre el agua y la atmósfera (Ceq).

Cuando kg,T es conocido para un gas a una temperatura T1 (Kg1,T1), es posible calcular un valor diferente de kg,T a una temperatura T2 (Kg2,T2) gracias a la relación entre los números de Schmidt (Sc) (Jähne, 1987):

Elexponente n de la Ecuación 2 varía de 1 a -0,67 según el proceso dominante delfenómeno de la difusión (Jähne, 1987). El número de Schmidt (Sc) es calculado parael CO2 con la Ecuación propuesta por Wanninkhof(1992):

Donde T es la temperatura superficial del agua en ºC.

Jähne (1987) proponen una solución para la Ecuación 2 en función a la Ecuación 3:

Donde k600 es la velocidad de transferencia del gas normalizada por el número de Schmidt de 600 (Sc = 600, para CO2 a 20ºC), y sobre la superficie de un lago n puede asumirse como 0,66 para velocidades de viento inferiores a 3,5 m.s-1, o 0,5 para velocidades superiores a ese valor (Guérin et al., 2007).

Ha quedado establecido que los flujos difusivos de un gas dependen tanto de la concentración superficial de dicho gas, como del coeficiente de velocidad de transferencia del mismo (ver Ecuación 1). Para el primer caso, si bien para este estudio se podrían usar las concentraciones superficiales de CO2 medidas en campo en el proyecto GEITRO (2012), la aproximación que se tiene en cuenta en esta investigación ha considerado la posibilidad de desarrollar un modelo capaz de establecer una relación entre dichas concentraciones y las temperaturas superficiales del agua. En efecto, las concentraciones de CO2 en el agua tienen ocurrencia en función a la temperatura, las presiones parciales del gas, la alcalinidad y la salinidad del medio. En cuerpos de agua continentales con pequeños cambios en alcalinidad y salinidad (lo cual ocurre en la mayoría de los casos), la variación de concentraciones de CO2 están, por ende, principalmente gobernadas en estos casos por los cambios en la temperatura y las presiones del gas en agua. Para el caso de los océanos, existen modelos ampliamente documentados en los que se correlacionan las presiones parciales (pCO2) o concentraciones superficiales del gas, con las temperaturas superficiales. Entre los más utilizados está el modelo propuesto por Stephens et al, (1995); sin embargo, se estima que no existe evidencia de modelos reportados para el caso de ecosistemas acuáticos continentales.

La aproximación propuesta en este estudio consiste entonces en el uso de los datos de campo relacionados con la medición de concentraciones superficiales y temperaturas superficiales (GEITRO, 2012). El objetivo es ajustar entre ambas variables una relación matemática que permita estimar dichas concentraciones a partir de valores conocidos de temperatura superficial del agua. Esta temperatura, para efectos de esta investigación, sería obtenida a partir de imágenes satelitales tal como se describirá más adelante.

Siguiendo la misma lógica para la solución de la Ecuación 1, se ha establecido que, para la estimación de los flujos difusivos de un gas se requiere también conocer la velocidad de transferencia del mismo (kg,T). Igualmente que el caso anterior, si bien podría usarse para este estudio las velocidades de transferencia CO2 estimadas a partir de datos de campo por el proyecto GEITRO (2012), esta investigación ha considerado la posibilidad de desarrollar un modelo adicional que conlleve a deducir, para el caso de Riogrande II, un modelo ajustado entre tales valores de velocidad de trasferencia y la velocidad del viento.

Para los lagos y embalses, el coeficiente de velocidad de transferencia de un gas (k en la Ecuación 1) ha sido altamente documentado y parametrizado como una función de la velocidad del viento y normalizado usando el número de Schmidt de 600 (K600) (Cole y Caraco, 1998; Crusius y Wanninkhof, 2003).

Li, (2018) sintetiza cinco de los modelos más confiables para la estimación de K600 (ecuaciones 5, 6, 7, 8, 9):

Esta relación es precisapara velocidades de viento tomadas a 10 metros de altura (U10) y cuya magnitudno supera los 3.7 m.s-1.

Esta relación está ajustadapara velocidades del viento U10 superiores a 2 m.s-1.

La cual no presentarestricciones para límites de velocidad del viento.

Igualmente, sinrestricciones para límites de velocidad del viento.

Sin restricción para velocidad del viento.

Para cada una de estas relaciones, K600 (cm.h-1) es la velocidad de transferencia ajustada al número de Schmidt de 600 para el CO2 a 20ºC. U10 (m.s-1) es la velocidad del viento ocurrida a 10 m por encima del espejo de agua.

Conociendo entonces los valores de K600 a una velocidad del viento (U10), por medio de cualquiera de los modelos antes descritos, es posible deducir los valores de velocidad de transferencia del CO2 (kg,T) para estimar los flujos difusivos del gas (F(CO2)) según lo dispuesto en la Ecuación 1.

Se escogieron los productos MODIS Terra, generados por la misión espacial MODIS, para el objeto de este estudio, teniendo en cuenta exclusivamente la disponibilidad de imágenes sobre el embalse de RGII para las mismas fechas en que existen las mediciones de campo descritas anteriormente. MODIS Terra observa la totalidad de la cobertura terrestre cada 1 a 2 días con adquisición de datos en 36 bandas espectrales diferentes. Sus órbitas están ubicadas a 705 km de altura y proveen una resolución espacial de 1000 m para las bandas que registran la Temperatura de la Superficie Terrestre (LST por sus siglas del inglés Land Surface Temperature) (Fisher, 2013).

Wan (2013) propone la estimación de la LST como la muestra la Ecuación 10:

Donde: T31 y T32 son las bandas 31 y 32 de MODIS; Ɛ son el promedio de la emisividad de la superficie de las bandas 31 y 32 de MODIS, C, A1, A2, A3, B1 y B3 son coeficientes de regresión (Wang, Liang y Meyers, 2008).

El procedimiento específico consistió en descargar las imágenes disponibles para las fechas de interés, es decir, las mismas fechas en las cuales se realizaron los trabajos de campo sobre el embalse de RGII dentro del marco del proyecto GEITRO (2012). La descarga se realizó del portal web de la NASA, llamado Centro de Archivos Activos Distribuidos de Nivel-1 y Atmósfera (LAADS-DAAC). A través de esta plataforma, se buscaron, seleccionaron y ubicaron las fechas necesarias para su posterior descarga y tratamiento (NASA, 2018), asegurándose de acceder a los productos de Temperatura y Emisividad de la Superficie Terrestre (Wan, 2013).

Después de la preselección y descarga de los productos, se proyectaron las imágenes al sistema de coordenadas UTM zona 18N con Datum WGS 1984, usando la herramienta MODIS Reproyection Tools (MTR) de acuerdo con Dwyer & Schmidt (2006).

Los valores de las temperaturas superficiales o LST de las imágenes seleccionadas, inicialmente en escala absoluta de Kelvin, fueron convertidas a escala relativa en grados Celsius empleando la siguiente relación:

Donde DN es el Número Digital estimado por los algoritmos a bordo del satélite según la Ecuación 1.12, y corresponde a los valores de LST día o noche respectivamente contenidos en las imágenes seleccionadas.

El procedimiento del cálculo se realizó mediante la calculadora de raster en el software ArcGIS 10.6. Estas temperaturas se validaron teniendo en cuenta las temperaturas medidas en campo para una misma fecha y posteriormente usadas para la estimación de concentraciones superficiales de CO2 en el embalse de RGII.

Del proyecto GEITRO (2012), se extrajo el conjunto de mediciones de flujos difusivos y concentraciones de CO2 realizadas durante los trabajos de campo (12 en total), así como los datos relacionados con las temperaturas superficiales del cuerpo de agua y velocidades del viento actuando sobre el hidrosistema registrados a cada instante de tiempo en que se medían concentraciones y flujos difusivos. El número de trabajos de campo son relacionados por fecha (Tabla 1).

A partir de la medición de los flujos difusivos de CO2 mediante el método de las cámaras flotantes, la temperatura superficial del agua y las concentraciones superficiales del biogás, el desarrollo del proyecto GEITRO (2012) estuvo en capacidad de deducir los valores de velocidad de transferencia (Kg, T) y el número de Schmidt (ScT) relacionados en la Ecuación 2 en cada uno de los puntos de control (Figura 1). En consecuencia, se pudo estimar el valor de velocidad de transferencia normalizado a K600, y se pudo determinar la relación con la velocidad local del viento, que se midió en la estación meteorológica instalada en el embalse.

Teniendo en cuenta el conjunto de medidas de campo antes mencionadas, particularmente concentraciones y temperaturas superficiales, fue posible proponer un modelo propio de concentraciones superficiales para el embalse de RGII que relaciona esas dos variables usando el software de Microsoft EXCEL (este fue usado en permanencia para todo el análisis de los datos que son presentados en adelante).

La Figura 2 presenta la mejor correlación encontrada (R2=0,77; p<0,05) entre el conjunto de datos disponibles. El número total de datos correlacionados (n) fue de 23, lo que corresponde a un 70% del total de datos disponibles dentro de la base de datos descrita en la sección anterior.

La función exponencial encontrada y que será, enconsecuencia, usada para modelar en adelante las concentraciones superficialesde CO2 a partir de temperaturas superficiales medidas en campo, opor cualquier otro método disponible, está dada por la Ecuación 12.

Donde, Cw es la concentración superficial de CO2 en el agua expresada en nmol.L-1, y T es la temperatura superficial del agua expresada en ˚C.

Esta relación fue utilizada entonces para modelar las concentraciones superficiales a partir de temperaturas estimadas por teledetección.

Si bien este estudio pudo haber utilizado directamente los valores de Kg,T calculados por el proyecto GEITRO (2012), la intención es la de emplear herramientas que conlleven a disminuir las intervenciones directas en campo. En este sentido, un resultado adicional de esta investigación es estimar los valores de K600 a partir de valores conocidos de velocidad del viento, teniendo en cuenta la relación bien documentada que existe entre ambas variables (Li, 2018).

Con el fin de definir el mejor modelo para estimar los valores de K600 para el caso de RGII; se ha aplicado cada una de las Ecuaciones 5 a 9, empleando los valores de velocidad de viento (U10) medidos en campo; y se han comparado los valores estimados por dichos modelos con los valores de K600 ya deducidos dentro del marco del proyecto GEITRO (2012).

Los resultados de la aplicación de estos 5 modelos son presentados en la Figura 3. Como se puede observar en esta figura, solo uno de los modelos, el propuesto por Guérin et al. (2007), presenta una correlación que podría considerarse aceptable para poder estimar con la mejor precisión posible los valores de K600 (R2 = 0,62; p<0,05). Por el contrario, los resultados obtenidos a partir de los otros modelos no son satisfactorios. Esto no puede interpretarse como que los modelos sean errados. Indica, simplemente, que ninguno de ellos se ajusta al comportamiento que se observó en torno a K600 dentro del embalse RGII durante los años 2012 a 2014 cuando se realizaron los muestreos de campo. En efecto, la mayoría de estos modelos han sido desarrollados para ecosistemas acuáticos en zonas boreales o tropicales de bajas altitudes. En el caso de RGII estamos hablando de un embalse tropical de montaña, y hasta donde se tiene conocimiento ningún modelo ha sido propuesto para este tipo de ecosistema.

Graciasa la disponibilidad de datos de campo, resultó posible deducir un modelo propiode K600 para el embalse de RGII, al correlacionar las velocidadesdel viento medidas en campo (y ahora corregidas a U10) con los valores de K600deducidos a partir de los datos mediante la aplicación de las Ecuaciones 1, 3 y4. La Ecuación 13 corresponde al modelo de deducido para RGII en el marco deeste estudio.

En la Figura 3f se observa la correlación (R2= 0,63; p<0,05) de los datos estimados por el modelo y los datos conocidos de K600.

Del contexto anterior se puede deducir que, si bien el modelo propuesto por Guérin et al. (2007), particularmente aplicable a embalses tropicales de zonas bajas, presenta una correlación aceptable con el modelo deducido en el marco de este estudio, aunque presenta una correlación un poco mejor que la del modelo antes mencionado (R2=0,64, p<0,05). Está claro que su aplicabilidad está adaptada a un ecosistema acuático tropical de montaña, específicamente al embalse de RGII. En este sentido, resulta más consecuente el uso de este último modelo para estimar los valores de K600 a partir de valores conocidos de U10 en el embalse de RGII.

Al final del proceso de selección de imágenes seencontraron 8 productos MODOS Terra de día paralas fechas en que hubo mediciones de campo del proyecto GEITRO (2012). Para cadauna de esas fechas se dedujeron las temperaturas superficiales sobre los puntosde control R2, R4, R7 y R8 (disponibles con datos de campo de flujos difusivossegún losreportados por el Proyecto GEITRO(2012)). La Figura 4 muestra el resultado de este proceso.

Los valores de temperatura deducidos a partir de imágenesMODIS se correlacionaron con los datos de temperatura superficial medidos en campo.La Figura 5 muestra el resultado de dicha correlación. Este resultado (R2=0,87;p<0,05), demuestra la buena precisión que presenta el sensor en MODIS Terrapara la adquisición de este tipo de información, sobre ecosistemas acuáticoscontinentales, particularmente en zonas tropicales de alta montaña, como es elcaso de RGII. Estos resultados permiten concluir que los posibles errores sobrela estimación de flujos de CO2 en el embalse RGII tienen una mayor relación conel uso de los modelos antes descritos, que con la pequeña diferencia entre losdatos de temperaturas medidos en campo y estimados por teledetección.

Con base en la información de campo y en el análisis de resultados disponibles para el embalse RGII, junto con la aplicación de los modelos ajustados para la estimación de concentraciones superficiales de CO2 a partir de temperaturas superficiales, y para la estimación de K600 (lo que lleva a encontrar valores ajustados del coeficiente de velocidad de transferencia de este biogás (kg,T)), se han estimado los flujos difusivos de CO2 (F(CO2)) sobre el embalse a dos niveles diferentes: i) estimación de flujos con el propósito de validar los modelos propuestos; y, ii) estimación de los flujos empleando los datos de temperaturas superficiales deducidos a partir de imágenes satelitales del MODIS Terra.

En el primer caso, se llevó a cabo la estimación de flujos difusivos de CO2 (F(CO2)) a partir del modelo deducido para concentraciones superficiales en el marco de esta investigación (Ecuación 12), y del modelo propuesto para la estimación de K600 (Ecuación 13), con el fin de establecer la correlación y consecuente grado de error que presentan estos modelos al comparar sus resultados de estimaciones de flujos con relación a los flujos calculados a partir de los datos de campo ampliamente descritos en este documento por el proyecto GEITRO (2012). La Figura 6 muestra la correlación entre ambas variables.

Al observar la Figura 6 se puede concluir que la aplicación de los modelos usados en el marco de esta investigación, han permitido reproducir de manera más que aceptable (R2 =0,78; p<0,05, y con un error relativo sobre el conjunto de datos de 15 ± 0,07% en promedio), los flujos difusivos de CO2 medidos en el embalse de RGII entre los años 2012 y 2014.

Para el segundo caso, los resultados antes presentados ofrecen la confiabilidad necesaria para emplear las temperaturas superficiales medidas por MODIS Terra para la estimación de los flujos difusivos sobre el embalse. En este caso se aplicó la misma metodología anterior: a partir de los modelos propuestos en las ecuaciones 11 y 12, se estimaron los valores de concentraciones superficiales de CO2 y K600. Estos resultados fueron posteriormente usados en las ecuaciones 4 y 1 para obtener los valores de Kg,T y, finalmente F(CO2), respectivamente. La Figura 7 muestra el resultado de la correlación obtenida después este proceso.

De la Figura 7 se puede concluir que, si bien la correlación encontrada (R2=0,71; p<0,05) entre los flujos estimados a partir de los modelos propuestos, los cuales hacen uso de las temperaturas superficiales obtenidas con datos del satélite MODIS Terra, esta es inferior a la correlación encontrada en el proceso de validación de dicho modelo (Figura 6); estos resultados siguen siendo bastante aceptables. El uso de las temperaturas superficiales observadas por MODIS Terra usadas en la Ecuación 3 para la estimación del número de Schmidt, posteriormente K600 y finalmente Kg,T; puede ser aceptado para la estimación de F(CO2) a partir de la diferencia de concentraciones superficiales de CO2 y la concentraciones del gas al equilibrio con la atmósfera, junto con el Kg,T según la ecuación 1.

En la actualidad, no se conocen en su totalidad los aportes de emisiones de CO2 en embalses tropicales; sin embargo, diversos estudios (Guérin y Abril, 2007; Guérin et al., 2008, entre otros) sugieren que los embalses pueden contribuir en un 4% de las emisiones globales de origen antrópico de CO2, magnitud que no es despreciable y que amerita una cuantificación más precisa para, así mismo, ajustar mejor los balances globales de carbono propuestos por el mismo IPCC.

El IPCC (2011), junto con otra serie de estudios, tales como Raymond et al. (2013) y Werhli (2013), sostienen que existen muy pocas mediciones en regiones tropicales, no solo a nivel de embalses, sino a nivel de ecosistemas acuáticos en general y que, por el momento, la poca información disponible resulta insuficiente para determinar alguna conclusión relacionada con los aportes de emisiones de GEI a partir de dichas regiones. Si bien hasta el momento han existido importantes esfuerzos con el fin de aportar conocimiento en dicho dominio, los cuales están representados en los estudios como los desarrollados por Guérin y Abril (2007) y Guérin et al., (2007) sobre el embalse de Petit-Saut en Guyana francesa, Kemenes et al., 2007; Kemenes et al.,2011 en la represa de Balbina en Brasil, Chanudet et al., (2011) en los embalses de Nam Ngum y Nam Leuk, Fearnside (2016); Fearnside (2009) y Demarty y Bastien (2011) en diferentes embalses tropicales, que son alternativas necesarias, aunque insuficientes, para determinar con precisión la contribución de emisiones de GEI de embalses en el trópico. Los elevados costos ligados al desarrollo de trabajos de campo periódicos, compra y mantenimiento de equipos de medición, tomas y análisis de muestras, entre otros, dificultan la ejecución el desarrollo de esta labor. En este sentido, se hace inminente el desarrollo de modelos que permitan estimar de forma acertada el conjunto de emisiones de estos ecosistemas a partir de información que sea accesible por medio de sensores remotos (estaciones de medición en tiempo real, teledetección, etc.).

En este sentido, este estudio tuvo como objetivo fundamental aplicar el uso de imágenes satelitales en diferentes modelos para reconstituir flujos difusivos de CO2 medidos en campo sobre el embalse de Riogrande II (RGII) y poder así concluir sobre la viabilidad del uso de este tipo de tecnologías en la estimación de tales emisiones de este biogás en ecosistemas acuático tropicales.

Como resultado de esta investigación, se ha propuesto un modelo para estimar concentraciones superficiales de CO2 aprovechando la relación que existe entre este parámetro y la variación de las temperaturas superficiales del agua. A nuestro conocimiento, esta sería la primera vez que se reporta un modelo de estas características para un ecosistema acuático tropical de alta montaña. Asimismo, se ha logrado evaluar diferentes modelos reportados en la literatura y proponer un modelo propio para la estimación del coeficiente de velocidad de transferencia normalizado sobre número de Schmidt de 600 (K600), lo que conlleva estimar, de manera consecuente, los coeficientes de velocidad de transferencia (kg,T) para el CO2 en el embalse de RGII. El hecho de proponer un modelo propio para estimar K600 a partir de su relación ampliamente documentada con valores de velocidad del viento (U10), surge de la baja precisión que tuvieron dichos modelos al intentar representar los valores de K600 medidos en campo. Dicho modelo, presentado en la Ecuación 12, presentó una correlación ligeramente superior a la propuesta por Guérin et al. (2007) (R2=0,64), por lo que se ha aceptado como una solución más ajustada para la modelación del comportamiento del sistema de RGII.

Estando en capacidad de estimar concentraciones superficiales de CO2, concentraciones al equilibrio con la atmósfera y el cuerpo de agua, así como los coeficiente de velocidad de transferencia del CO2 (kg,T), fue posible modelar los flujos difusivos del biogás (F(CO2)) según la Ecuación 1, para las fechas en que dichos flujos fueron medidos por el proyecto GEITRO (2012), lo que ha permitido establecer, en consecuencia, el dominio de validez de los dos modelos desarrollados en el marco de esta investigación y que vienen de ser discutidos. Los resultados presentados permiten concluir que los modelos propuestos logran reproducir de manera más que aceptable (R2 =0,78, con un error relativo sobre el conjunto de datos de 15 ± 0,07% en promedio), los flujos difusivos de CO2 medidos en el embalse de RGII entre los años 2012 y 2014.

Teniendo en cuenta los resultados antes expuestos, fue posible emplear las temperaturas superficiales medidas por MODIS Terra para la estimación de los flujos difusivos de CO2 sobre el embalse. Los resultados fueron satisfactorios en la medida en que la correlación obtenida después este proceso entre los flujos estimados a partir de los modelos propuestos, los cuales hacen uso de las temperaturas superficiales obtenidas a partir de satélite MODIS-Terra fue de R2=0,71. Si bien esta correlación resultó inferior a la correlación encontrada en el proceso de validación de los modelos, estos resultados siguen siendo aceptables.

Se ha logrado demostrar que, de forma global, el uso de las temperaturas superficiales observadas por MODIS-Terra para la estimación del número de Schmidt, a partir del cual es posible estimar los valores de K600 usando el modelo desarrollado en el marco de este estudio, para finalmente calcular valores de kg,T; así como el uso de esas mismas temperaturas para la estimación de concentraciones superficiales de CO2 mediante la aplicación del modelo propuesto en esta investigación, puede ser aceptado para la estimación de flujos difusivos de este biogás (F(CO2)). Queda entendido que los valores de velocidad del viento (U10) utilizados en el marco de esta investigación y necesarios para la estimación de K600, fueron aquellos medidos en campo por la estación meteorológica dispuesta por el proyecto GEITRO. Sin embargo, dependiendo del lugar de la zona de estudio, estos valores de velocidad del viento pueden ser medidos por otros tipos de medio remotos como AEOLUS (ESA,2018) o por estaciones de las redes meteorológicas nacionales o locales, con lo cual, se reduciría le necesidad de medir esta variable en campo. No fue objeto de esta investigación el evaluar otras fuentes de información que aporten velocidades de viento. El asunto queda entonces planteado como expectativa para otras investigaciones relacionadas.