Recibido: 10 junio 2021
Aceptado: 11 noviembre 2021
Disponible: 17 diciembre 2021
Con el fin de representar la precipitación y evaporación total mensual en una cuenca hidrográfica del Orinoco colombiano, este trabajo evaluó la capacidad de los modelos climáticos regionales incluidos en el Experimento regional coordinado de reducción de escala (CORDEX-CORE). Para ello, complementariamente, se incluyeron datos de precipitación y evaporación total de fuentes como Climate Hazards Center InfraRed Precipitation with Station data (CHIRPS), el reanálisis atmosférico (ERA5), Global Precipitation Climatology Center (GPCC) y Global Land Evaporation Amsterdam Model (GLEAM). Las comparaciones entre los ensambles de los modelos y las observaciones se hicieron utilizando métodos gráficos y métodos cuantitativos, entre ellos: diagramas de cajas, porcentajes de sesgo, eficiencia de Nash-Sutcliffe, entre otros. Los resultados evidencian que los valores promedio de precipitación están adecuadamente representados, en términos de su temporalidad y magnitud, por el ensamble del modelo RegCM, mientras que los valores promedio de evaporación total están mejor representados por el ensamble del modelo REMO en términos de la temporalidad, más no en su magnitud. Por otra parte, las estimaciones de caudal de largo plazo evidencian que los valores de evaporación total proporcionados por los modelos permiten una adecuada estimación del caudal promedio de largo plazo, pero no la adecuada estimación del ciclo anual de caudales. Este trabajo es pionero en la evaluación de los datos de precipitación y evaporación total mensual suministrados por CORDEX-CORE en el Orinoco colombiano, sienta precedentes para la incorporación de datos de modelos regionales para fines hidrológicos en zonas poco instrumentadas del país, y es el primer paso hacia la evaluación de escenarios regionalizados de cambio climático.
Palabras clave: Precipitación mensual, evaporación mensual, balance hídrico, comparación de modelos, modelos climáticos regionales, evaluación de desempeño.
We evaluate the capacity of the Regional Climate Models, included in the regional scale reduction experiment (CORDEX-CORE), to represent the monthly precipitation and evaporation in a watershed of the Colombian Orinoco region. In addition, we use precipitation and evaporation data from sources such as Climate Hazards Center InfraRed Precipitation with Station data, ERA5 atmospheric reanalysis, Global Precipitation Climatology Center and Global Land Evaporation Amsterdam Model. The comparisons between the assemblies of the RCMs and the observations were done using graphical and quantitative methods, among them: box plots, percentages of bias, Nash-Sutcliffe efficiency, among others. Our results allow concluding that the average precipitation values are adequately represented, in terms of their temporality and magnitude, by the assembly of the RegCM model while the average values of total evaporation are better represented by the assembly of the REMO model in terms of temporality, but not in its magnitude. Furthermore, the long-term flow estimation shows that the total evaporation values provided by the RCMs allow an adequate estimation of the long-term average flow, but not the adequate estimation of the annual flow cycle. This work is a pioneer in the evaluation of the monthly precipitation and evaporation data provided by CORDEX-CORE in the Colombian Orinoco, sets precedents for the incorporation of regional models data for hydrological purposes in poorly instrumented areas of the country, and is the first step towards the evaluation of regional climate change scenarios.
Keywords: Monthly precipitation, Monthly evaporation, Water Balance, Model Comparison, Regional Climate Models, Performance assessment.
El territorio continental de Colombia está conformado por las regiones naturales: Caribe, Pacífico, Andina, Amazonía y Orinoquía. Ésta última recibe su nombre debido a que la mayor parte de sus ríos se encuentran en la cuenca del río Orinoco y está compuesta por los departamentos de Arauca, Casanare, Meta y Vichada, abarcando el 18 % del territorio continental nacional, limitando al norte y al este con Venezuela, al sur con la región Amazónica y al oeste con la región Andina [
Mientras el Orinoco experimenta su temporada seca (noviembre-marzo), la región de bajas presiones atmosféricas, conocida como zona de convergencia intertropical (ZCIT), migra hacia el sur, intensificando los vientos alisios sobre el noreste de Venezuela [
Los modelos climáticos globales (MCG) permiten representar e investigar la dinámica climática global y pretenden simular la respuesta a cambios de patrones en variables del sistema climático, tales como el incremento de los gases de efecto invernadero o la radiación [
Los MCG representan el clima utilizando una rejilla tridimensional sobre el globo terráqueo, normalmente con una resolución horizontal superior a los 100 km, mientras la atmósfera se representa a través de 10 a 20 capas verticales y 30 para el océano. Esta resolución es muy gruesa para observar detalladamente los procesos de la superficie terrestre y del clima de manera directa en una región o lugar específico [
Dentro de los esfuerzos de modelación regional, se encuentra el Experimento regional coordinado de reducción de escala (CORDEX), cuya misión principal es regionalizar escenarios de cambio climático en el ámbito global y, a su vez, generar un espacio para que la comunidad que esté interesada en esos escenarios pueda compartir información, tecnología, resultados y soporte técnico, generando así una mejor interacción entre los modeladores y las personas interesadas en realizar estudios sobre vulnerabilidad y adaptación al cambio climático [
Los MCR REMO y RegCM son modelos atmosféricos hidrostáticos, que emplean más de veinte niveles para representar la atmosfera y se diferencian, principalmente, en que REMO mantiene las mismas parametrizaciones para todo el planeta, mientras que RegCM emplea diferentes parametrizaciones por dominio [
La evaluación del desempeño de los MCR se desarrolla por medio de métricas de valor agregado con respecto a los MCG [
Dichos trabajos han mostrado avances en la modelación regional con respecto a la modelación global en la representación de la precipitación o evaporación, dependiendo de la región, temporada o cuenca [
La reducción de escala de los MCG mediante los MCR permite diagnosticar, evaluar y gestionar los impactos del cambio climático para la adecuada planificación de los recursos hídricos, la conservación de los recursos naturales, el riego de plantaciones de gran escala, entre otras actividades de los diferentes sectores públicos y privados [
Especialmente, en cuencas hidrográficas con poca instrumentación, procesos de gran interés como la precipitación, la evaporación, la infiltración y el flujo subterráneo son difíciles de observar de forma continua, y sin ambigüedades, debido a la heterogeneidad de las cuencas y los procesos hidroclimáticos. Además, dichos procesos requieren de un adecuado monitoreo, evaluación y diagnóstico para modelar adecuadamente sus posibles cambios e impactos en escalas espaciales y temporales más finas [
En este sentido, los MCR ofrecen información sobre variables hidrológicas (precipitación, caudales, entre otras) que pueden ser útiles para evaluar el impacto hidrológico en cuencas hidrográficas considerando escenarios de cambio climático [
Para el caso particular de Colombia, pocos estudios han utilizado información de MCG y MCR para la verificación de balances hídricos de largo plazo y modelación hidrológica considerando escenarios de cambio climático [
Este estudio tiene como objetivo evaluar si los MCR, incluidos en el corte de Suramérica para el experimento CORDEX-CORE, representan adecuadamente la precipitación y la evaporación total en el Orinoco colombiano. Adicionalmente busca hacer un análisis comparativo entre observaciones de caudal in-situ y las estimaciones de caudal de largo plazo utilizando información de los MCR. Se evaluará el desempeño de las simulaciones del experimento por primera vez en esta zona, con respecto a las observaciones obtenidas por el Instituto de Hidrología, Meteorología y Estudios Ambientales (IDEAM) y otras bases de datos interpoladas espacialmente para el período de referencia utilizado en el quinto informe de evaluación (AR5) del Grupo Intergubernamental de Expertos sobre el Cambio Climático (IPCC), 1986-2005. Por lo tanto, este trabajo será útil para avanzar en aspectos como: (1) La verificación del desempeño, las simulaciones históricas ofrecidas por los MCR para la precipitación y la evaporación total sobre Colombia; (2) La incorporación de la información suministrada por los MCR para análisis hidrológicos e históricos en zonas con poca instrumentación, mediante el estudio del desempeño de los MCR y la estimación del caudal de largo plazo en la cuenca de análisis; (3) Brindar una herramienta metodológica para la selección de información que sirva para el diagnóstico y modelación de los impactos hidrológicos en cuencas hidrográficas de Colombia considerando la información del experimento CORDEX-CORE.
2.1 Datos
Como zona de estudio se seleccionó la cuenca hidrográfica del río Meta, la cual tiene un área de 78 072 km2 (Figura 1). Esta zona está definida aguas arriba de la estación limnigráfica Aguaverde (código 35267080) del IDEAM, ubicada en el municipio La Primavera, en el departamento de Vichada. Se tomó información mensual de precipitación y evaporación total a una resolución espacial de 0.25° x 0.25° de los MCR REMO2015 [
Se tomaron 85 estaciones de precipitación, con resolución temporal mensual y un porcentaje de datos faltantes inferior al 15 % en el periodo de referencia 1986-2005. Se utilizaron los datos de caudal de la estación Aguaverde (código 35267080), para el mismo periodo de referencia. Además, considerando que las series de tiempo seleccionadas tienen un porcentaje de datos faltantes inferior a 15 %, dichos datos faltantes fueron llenados utilizando los valores mensuales multianuales [
Además, se usó información de precipitación a escala temporal mensual y resolución espacial 0.25° x 0.25° de las siguientes bases de datos: Climate Hazards Center Infra Red Precipitation with Station (CHIRPS) [
2.2 Métodos
Para el análisis comparativo entre las fuentes de información se promediaron los datos de precipitación del IDEAM para todas las estaciones en la cuenca de estudio. Por su parte, con la información obtenida de ERA5, GPCC, CHIRPS y GLEAM se hicieron promedios sobre la cuenca de estudio. Además, para cada uno de los MCR se realizó un ensamble o promedio de las tres simulaciones disponibles, ya que el ensamble permite disminuir las incertidumbres de las regionalizaciones de los MCR [
2.2.1 Métodos gráficos
El diagrama de caja es una representación gráfica de la distribución y variación de los datos, así como los datos atípicos. La altura de la caja representa el rango intercuartílico entre el primer cuartil (porción inferior) y el tercer cuartil (porción superior). La línea horizontal de la caja representa la mediana mientras que las líneas verticales (bigotes) que salen de la caja se extienden hasta el mínimo y máximo del conjunto de datos. Los valores representados por puntos debajo o por encima de los bigotes se consideran valores atípicos [
2.2.2 Métodos cuantitativos
Para una comparación inicial entre los datos de referencia y los resultados de los MCR se utilizaron el coeficiente de correlación de Pearson (R), el coeficiente de determinación (R2), el error típico y la diferencia porcentual. Estas métricas permiten hacer un diagnóstico preliminar sobre la asociación lineal y las diferencias entre los valores de referencia y simulados en términos de su magnitud. Posteriormente, se implementaron para el análisis comparativo las métricas de desempeño compiladas por [
En (2) RSR se define como la relación entre la raíz del error cuadrático medio (RMSE) y la desviación estándar de los datos medidos (STDEVobs). El valor óptimo es cero e indica que la simulación del modelo es perfecta, sin variación entre los datos observados y simulados. Entre menor sea el RSR, mejor es el rendimiento del modelo.
En (3) la NSE se define como una métrica normalizada que ajusta los valores entre los datos observados y los simulados para medir la dispersión entre estos. Los valores posibles del NSE están en el rango de menos infinito hasta uno; cuando es igual a uno, indica una coincidencia perfecta entre las simulaciones y las observaciones, entre más cerca a uno esté el resultado, más preciso es el modelo.
Para todas las métricas, n es el número total de datos de la serie de tiempo, 𝑌𝑖𝑜𝑏𝑠 es el valor en el tiempo i de la variable observada, 𝑌𝑖𝑠𝑖𝑚 es el valor en el tiempo i de la variable simulada y 𝑌𝑖𝑚𝑒𝑎𝑛 es el promedio de los datos observados. Para evaluar el desempeño de los MCR con base en las métricas seleccionadas se tomaron como referencia los rangos definidos por [
Métricas | |||
Desempeño | RSR | NSE | PBIAS (%) |
Muy bueno* | 0.00 ≤ RSR ≤ 0.50 | 0.75 < NSE ≤ 1.00 | PBIAS < ±10 |
Bueno** | 0.50 < RSR ≤ 0.60 | 0.65 < NSE ≤ 0.75 | ±10 ≤ PBIAS < ±15 |
Satisfactorio+ | 0.60 < RSR ≤ 0.70 | 0.50 < NSE ≤ 0.65 | ±15 ≤ PBIAS < ±25 |
Insatisfactorio° | RSR > 0.70 | NSE ≤ 0.50 | PBIAS ≥ ±25 |
En (4) se tomará la ecuación de balance hidrológico de largo plazo [
Donde P representa la precipitación promedio de largo plazo, E la evaporación total promedio de largo plazo y R la estimación indirecta de la escorrentía promedio de largo plazo. Se estimará el valor de R con base en las simulaciones de los MCR y la información de referencia para precipitación y evaporación total. Es importante aclarar que la información sobre la variable E que ofrecen los MCR del experimento CORDEX-CORE corresponde a la evapotranspiración real [
La Figura 2 muestra la distribución espacial de la precipitación promedio multianual según las estaciones de IDEAM en la cuenca del río Meta. Se evidenció que la precipitación promedio mensual varía en la zona de estudio en un rango de 67 mm/mes a aproximadamente 600 mm/mes, donde los valores por encima de 250 mm/mes se encuentran en la zona del Piedemonte de la cordillera Oriental y los valores más bajos (67-100 mm/mes) en la zona del altiplano en los departamentos de Boyacá y Cundinamarca.
Particularmente, sobre el Piedemonte de la cordillera Oriental de los Andes se concentran más del 80 % de las estaciones pluviométricas disponibles (69 estaciones), lo que resulta en una estación pluviométrica cada 412 km2, mientras menos del 20 % (16 estaciones) se encuentran en el resto del territorio de estudio, una estación pluviométrica por cada 3100 km2.
Por otra parte, en la Figura 3 se muestran los campos de precipitación promedio mensual de los MCR con las condiciones de borde de tres los MCG: NorESM1-M, MPI-ESM-LR y HadGEM2. En general, se evidencia heterogeneidad en los valores promedio mensuales que arroja cada MCR. Por ejemplo, RegCM-Had exhibe alta precipitación (350-400 mm/mes) en el sur de la cuenca de estudio, mientras RegCM-MPI muestra alta precipitación en el norte.
De igual manera, REMO-MPI muestra precipitación alta en el Piedemonte (200-400 mm/mes), mientras REMO-NOR sugiere una precipitación promedio mensual casi homogénea en toda la cuenca de interés (100-200 mm/mes). Por lo tanto, la heterogeneidad que presentan los resultados particulares de los MCR en la estimación de la precipitación media mensual pone de manifiesto la importancia de utilizar ensambles para minimizar la incertidumbre en la selección de una realización en particular [
La Figura 4 muestra los ciclos anuales de precipitación y evaporación total sobre la zona de estudio, para las fuentes de referencia y los ensambles de MCR en el periodo 1986-2005.
La Figura 4-a evidencia que la mayoría de las fuentes de información de precipitación exhiben un ciclo anual con distribución unimodal. Es decir, sólo se presenta un mínimo y un máximo de precipitación mensual durante el año. Se evidenció una temporada de lluvia entre abril y noviembre, con precipitación máxima de 364 mm en el mes de junio. También una temporada seca entre diciembre y marzo, con precipitación mínima de 34.28 mm en el mes de enero, lo cual coincide con la temporada de actividad del Chorro del Orinoco y la ausencia de la ZCIT [
De los modelos de CORDEX-CORE, el ensamble del modelo REMO (línea azul) presenta los valores de precipitación más alejados de las observaciones, mostrando subestimación para el periodo noviembre-agosto con relación a los valores registrados por las estaciones del IDEAM (línea negra continua). También se muestra que el ensamble del modelo RegCM (línea roja) presenta valores promedio mensuales multianuales muy cercanos al promedio de las estaciones del IDEAM durante los meses de noviembre hasta abril. Sin embargo, durante los meses de mayo hasta junio dicho modelo sobrestima las observaciones, mientras que durante julio-noviembre, las subestima. Es de resaltar que ni REMO ni las fuentes satelitales logran representar de manera adecuada el máximo de precipitación reportado en el mes de junio por los datos del IDEAM. No obstante, no hay gran dispersión entre los valores observados por el IDEAM, los valores del ensamble del modelo RegCM y las demás fuentes de referencia.
La Figura 4-b muestra el ciclo anual de evaporación total para el mismo periodo de análisis sobre la cuenca de estudio, los resultados probaron que el ensamble del modelo RegCM difiere significativamente en orden de magnitud de los valores reportados por ERA5, GLEAM y el ensamble de las corridas del modelo REMO. En la Figura 4-c se presenta un acercamiento a la comparación gráfica entre el ensamble REMO y las fuentes ERA5 y GLEAM. Mientras ERA5 y GLEAM presentan variabilidad bimodal y evaporación total mensual promedio entre 80 y 100 mm/mes a lo largo del año, REMO sugiere un ciclo anual predominantemente unimodal con valores máximos alrededor de 140 mm/mes durante septiembre-octubre, y valores mínimos cercanos a 70 mm/mes durante abril-junio.
Lo que sugiere que ninguno de los dos ensambles hace una buena representación del ciclo anual de la evaporación total.
La Figura 5-a muestra los diagramas de cajas para la distribución de los valores mensuales de precipitación y evaporación total obtenidos de los ensambles de los MCR y las fuentes de referencia. Para la precipitación, se evidenció que todas las referencias estaban alrededor de 200 mm/mes, además las fuentes de información exhibieron mayor variabilidad en los valores asociados al primer cuartil (precipitación baja), mientras que los valores asociados al máximo del tercer cuartil (precipitación alta) presentaron similitud en orden de magnitud. También se demostró que la fuente que presenta valores promedio mensuales más cercanos a las observaciones es GPCC, teniendo un rango de variabilidad y una media similares a los del IDEAM. En cuanto a la evaporación total, la Figura 5-b muestra que las referencias ERA5 y GLEAM tienen valores alrededor de los 100 mm/mes, al igual que el ensamble del modelo REMO. Si bien REMO representa bien la mediana, presenta una mayor variabilidad en la distribución de la evaporación total con respecto a los valores de referencia.
Además, se revisó la variabilidad de los valores mensuales de los ensambles con relación a las observaciones utilizando diagramas de cajas mes a mes para los ensambles que presentaron mayor similitud en el ciclo anual promedio. La Figura 5-c muestra que los valores mensuales del ensamble RegCM (rojo) presentan menor variabilidad que los de referencia (azul), dicha diferencia en la dispersión de los valores mensuales puede atribuirse a que los registros pluviométricos son puntuales y no tienen una distribución espacial uniforme, mientras que la información de RegCM proviene de datos interpolados espacialmente en una retícula uniforme. La Figura 5-d muestra que los valores mensuales del ensamble REMO (naranja) presentan diferencias en la estacionalidad de la evaporación total, pero ambas fuentes de datos coinciden en que el periodo septiembre-diciembre presenta valores altos de evaporación total mensual entre 100 mm/mes y 140 mm/mes.
La Tabla 2 muestra los resultados de las métricas cuantitativas de rendimiento para la precipitación con relación a los datos IDEAM de los ensambles de MCR y otras fuentes de referencia. En general, se muestra que el ensamble de MCR que mejor simuló las observaciones fue RegCM. En este sentido, se puede resaltar que RegCM tiene una alta correlación (𝑅 = 0.91) con las observaciones obtenidas por el IDEAM, reflejando una importante asociación lineal entre las bases de datos. El coeficiente de determinación (𝑅2 = 0.82) indica que el modelo RegCM representa más del 80 % de la varianza de las observaciones. Mientras que el error típico es de 22.12 mm con relación al valor promedio mensual observado en los datos del IDEAM. Asimismo, el ensamble RegCM presenta una baja diferencia porcentual con un 3 % de desviación con respecto al IDEAM. La clasificación del rendimiento para cada una de las métricas se basa en los rangos presentados en la Tabla 1.
Fuentes | (mm) | (mm) | Error típico (mm) | Dif. | R | R2 | RSR | NSE | PBIAS |
RegCM | 201.83 | 123.50 | 22.12 | 3 % | 0.91 | 0.82 | 0.46* | 0.79* | 6.99* |
REMO | 180.86 | 116.37 | 47.71 | 13 % | 0.83 | 0.69 | 0.61** | 0.63+ | 11.40** |
CHIRPS | 204.75 | 113.98 | 22.90 | 1 % | 0.86 | 0.74 | 0.52** | 0.73** | 1.55* |
ERA5 | 197.36 | 94.61 | 33.15 | 5 % | 0.84 | 0.73 | 0.55** | 0.70** | 0.93* |
GPCC | 208.17 | 119.28 | 23.36 | -1 % | 0.85 | 0.73 | 0.55** | 0.70** | -0.51* |
Fuentes | (mm) | (mm) | Error típico (mm) | Dif. | R | R2 | RSR | NSE | PBIAS |
RegCM | 276.20 | 65.06 | 35.02 | 381 % | 0.67 | 0.45 | 16.56° | -273.08° | -181.77° |
REMO | 99.66 | 29.37 | 23.14 | 201 % | 0.52 | 0.27 | 2.23° | -3.95° | -2.69* |
GLEAM | 91.47 | 7.40 | 3.14 | 193 % | 0.77 | 0.59 | 0.86° | 0.26° | 6.65* |
Se siguió evidenciando que el modelo REMO es el que presenta mayor diferencia con relación a los valores de referencia de precipitación, este reporta una diferencia porcentual de 13 % por encima del valor de las observaciones IDEAM. También, se observó un error típico de 47.71 mm con relación a la media de las observaciones. Por su parte, para el ensamble REMO, métricas como el coeficiente de correlación (𝑅 = 0.83) y el coeficiente de determinación (𝑅2 = 0.69) presentan una reducción significativa en la representación de las observaciones comparado con el ensamble RegCM. Además, el 𝑃𝐵𝐼𝐴𝑆 ≌ 12.26 indica un aumento en el sesgo de REMO con relación a las observaciones IDEAM, mientras las métricas RSR y NSE cercanas a 0.60 reflejan un rendimiento satisfactorio pero inferior al del ensamble RegCM.
Por otra parte, para la evaporación total, se tomaron los datos de ERA5 como referencia y se evaluó el desempeño de los ensambles de los MCR, así como de los datos GLEAM. La Tabla 1 evidencia que el ensamble REMO presenta valores de evaporación total menos alejados de los reportados por ERA5. Sin embargo, ni REMO ni RegCM exhiben métricas de rendimiento satisfactorias para la evaporación total con relación a los datos de ERA5.
Para la estimación de caudales se utilizó la ecuación de balance hídrico de largo plazo (4) considerando la información de precipitación del ensamble con el mejor desempeño (RegCM), evaporación total del ensamble con el mejor desempeño (REMO), la evaporación total reportada por fuentes de referencia (GLEAM y ERA5) y la precipitación según las estaciones de IDEAM. Mientras que el caudal promedio multianual observado en la estación Aguaverde del IDEAM (código 35267080) es de 3346.28 m3/s, las estimaciones desarrolladas utilizando las bases de datos seleccionadas llevan a los siguientes resultados en magnitud y en porcentaje con respecto a las observaciones: IDEAM-Ens_REMO: 3237.83 m3/s (97 %); IDEAM-GLEAM: 3484.52 m3/s (104 %); IDEAM-E_ERA5: 3276.85 m3/s (98 %); Ens_RegCM-Ens_REMO: 3078.62 m3/s (92 %). En general, puede concluirse que tanto la estimación de caudal empleando la evaporación total del ensamble REMO y la información del IDEAM, la estimación de caudal usando evaporación total de GLEAM y ERA5, y la estimación del caudal empleando los dos mejores ensambles de MCR exhiben errores porcentuales inferiores al 5.0 % en la estimación del caudal promedio multianual.
No obstante, múltiples aplicaciones hidrológicas (Ej.: planeación del recurso hídrico, actividades agrícolas, planeación energética, entre otras) requieren de la adecuada caracterización y estimación del ciclo anual. Con este fin, se verificó si la información de los MCR permite hacer una adecuada estimación del caudal promedio mensual multianual. Al respecto, encontramos que al emplear la información de precipitación y/o evaporación total de los MCR no es posible hacer una estimación adecuada ni de la temporalidad del ciclo anual de los caudales ni de su magnitud promedio mensual. La Figura 6 muestra que los registros de caudal de IDEAM y las estimaciones exhiben un ciclo anual de caudales unimodal, pero estos ciclos difieren en las temporadas de caudales máximos y mínimos. Además, todos los MCR exhiben valores de evaporación total mensual superiores a la precipitación en el periodo octubre-marzo sugiriendo un déficit hídrico en la cuenca que difiere significativamente de las observaciones. Por su parte, durante la temporada abril-julio los caudales estimados utilizando información de los MCR arrojan valores superiores a los registrados por IDEAM, lo que sugiere que durante dicho periodo los MCR subestiman los valores de evaporación total que se presentan en la zona.
Gráficamente se identificó que el MCR que más se aproxima a las observaciones de precipitación obtenidas por el IDEAM para el Orinoco colombiano es RegCM. Sin embargo, cabe resaltar que las fuentes interpoladas satelitales también muestran valores similares a los valores obtenidos por el IDEAM para la zona de estudio y la escala temporal utilizada.
El análisis comparativo a través de métricas de rendimiento permite concluir y verificar que RegCM es el MCR que mejor representa el ciclo unimodal de la precipitación en la cuenca de estudio. En relación con esto, aunque RegCM es el modelo que muestra valores más aproximados a los del IDEAM, en algunas épocas del año, especialmente en mayo, subestima en gran medida las observaciones. Por otra parte, también se verificó que la precipitación según REMO presenta un mayor sesgo y dispersión que RegCM, lo cual nos permite concluir que tanto el ensamble de RegCM como las CHIRPS, GPCC y ERA5 podrían ser empleados como fuentes para el estudio de la precipitación en la zona de estudio. No obstante, la precipitación mensual del MCR REMO presenta resultados satisfactorios para la representación del ciclo anual promedio y podría ser utilizado para describir el ciclo anual de la precipitación en cuencas sin instrumentación.
Gráficamente es posible hacer una primera inspección de los resultados de evaporación total y determinar que ninguno de los MCR representa adecuadamente dicha variable. Una vez calculadas las métricas de desempeño, se encuentra que ninguna de las fuentes comparadas con ERA5 logra representar de forma satisfactoria la evaporación total. Estos resultados confirman las dificultades para la consecución de información de buena calidad de esta variable, de la cual no se cuentan redes consolidadas de medición en el país, por lo cual se recomienda para aplicaciones hidrológicas hacer una estimación indirecta de dicha variable a partir de la información de temperatura disponible.
Los ensambles de los MCR son una buena aproximación para la adecuada representación de la precipitación. Los resultados aquí expuestos confirman que el ensamble de RegCM representa adecuadamente la precipitación en la zona de estudio, esto podría estar relacionado con las diferentes parametrizaciones por dominio del modelo RegCM [
El análisis comparativo utilizando la ecuación de balance hídrico de largo plazo permite concluir que los datos de precipitación del ensamble de RegCM y los datos de evaporación total del ensamble REMO permiten una estimación adecuada del caudal promedio anual de largo plazo. Sin embargo, ni los datos de evaporación total reportados por REMO y RegCM, ni los reportados por las fuentes ERA5 y GLEAM, permiten hacer una estimación de la temporalidad y magnitud del ciclo anual de caudales en la cuenca de estudio. La selección de cada MCR dependerá del enfoque del estudio a desarrollar. Para estudios apoyados en el análisis del comportamiento temporal de las variables, será conveniente seleccionar el MCR que representa mejor el ciclo anual en la zona de estudio.
En general, los MCRs tie nen múltiples fuentes de incertidumbre, tales como la conceptualización y parametrización de cada MCR, la incertidumbre en los escenarios de emisión para proyecciones futuras en el contexto de cambio climático, incertidumbre en los procesos de reducción de escala (downscaling), incertidumbre en la sensibilidad de los modelos y la incertidumbre asociada a la adecuada representación de la complejidad y no-linealidad de la dinámica del sistema climático terrestre [
En particular, en el caso de este estudio, se encontró que los MCRs no hacen una adecuada representación de la evaporación total. Según la literatura, los problemas en la representación de esta variable son frecuentes debido a la dificultad para la parametrización de los procesos que configuran la evaporación total y a las pocas mediciones que se tienen en tierra para poder comprender la variable y validar las representaciones de las interacciones suelo atmósfera [
Como trabajo futuro, se propone el estudio de la regionalización de variables hidrológicas a través de la caracterización de patrones espaciales y de escalamiento en cuencas no instrumentadas utilizando información de los MCR.
CFV agradece al proyecto “Evaluación del desempeño y proyecciones para la representación de la disponibilidad de agua superficial de los modelos regionales en el dominio SAM-CORDEX para las cuencas del Río Amazonas y Del Plata” por su financiación durante el desarrollo de este proyecto. HDS agradece a Colciencias por la beca para doctorados nacionales convocatoria 617-2. Los autores también agradecen a las instituciones que permiten el acceso libre a sus bases de datos: Las simulaciones de MCR de CORDEX se encuentran disponibles gracias al Earth System Grid Federation en su página web (https://esgf.llnl.gov/), ERA5 está disponible en la página web del European Copernicus Programme (https://www.copernicus.eu/), CHIRPS se encuentra disponible gracias a la University of California Santa Barbara Climate Hazards Center (https://www.chc.ucsb.edu/data), GPCC está disponible gracias al Physical Sciences Laboratory (https://psl.noaa.gov/data/gridded/data.gpcc.html), GLEAM está disponible en la página web del proyecto (https://www.gleam.eu/) y la información observacional de estaciones en Colombia se encuentra disponible gracias al IDEAM (www.ideam.gov.co).
Los autores declaran no tener conflictos de intereses.
Carolina Florian Vergara se encargó de la obtención y procesamiento de los datos, preparación del manuscrito inicial y análisis de la información. Carolina Florian Vergara, Hernán D. Salas y Alejandro Builes-Jaramillo participaron en la conceptualización del estudio, en el análisis y discusión de resultados, escritura y edición del documento.