Colombia presenta una gran red hidrográfica en la región de los llanos orientales donde el Instituto de Hidrología, Meteorología y Estudios Ambientales IDEAM cuenta con estaciones hidrometeorológicas capaces de llevar registros de lluvia en el tiempo, siendo estas necesarias para la realización de estudios que permitan efectuar seguimiento a esta red y así fundamentar la toma de decisiones en sectores que se ven afectados por las variaciones climáticas que se presentan. De acuerdo con la situación mencionada, se pretende pronosticar el caudal circundante al Municipio de Yopal y obtener elementos que sirvan de base para prever y minimizar riesgos especialmente al sector del agro colombiano el cual, es característico de esta región. Se plantea un estudio de caso del caudal en la cuenca media del Río Cravo Sur en una serie temporal de periodicidad mensual desde el año 1975 hasta el 2014, luego se realiza un análisis exploratorio de datos el cual, permite visualizar los estadísticos descriptivos que sirven para evaluar el comportamiento general de la serie, posteriormente se analiza la misma con el fin de determinar el modelo más apropiado para el pronóstico.

Como resultado, se determinó que, a lo largo del tiempo, los caudales tienden levemente a aumentar, pero sin generar riesgos superiores a los que en la actualidad está expuesta la región especialmente en épocas de constantes lluvias; esta estimación se ajusta a las características propias de la estación y depende única y exclusivamente de su historia.

En la región de la Orinoquía o Llanos Orientales de Colombia como es conocida comúnmente, se presenta una variación climática en diferentes tiempos del año siendo que para el primer trimestre del mismo se cuenta con un tiempo de sequía como lo demuestran los datos históricos de los últimos 40 años los cuales, fueron suministrados por el Instituto de Hidrología, Meteorología y Estudios Ambientales IDEAM "institución pública de apoyo técnico y científico al Sistema Nacional Ambiental, que genera conocimiento, produce información confiable, consistente y oportuna, sobre el estado y las dinámicas de los recursos naturales y del medio ambiente" (IDEAM, 2014), estos datos, muestran que existen problemas de déficit hídrico en tiempos de sequía, acarreando situaciones fuera de lo normal en el suministro de agua a las poblaciones y el abastecimiento a las actividades del agro incluido la ganadería donde se ven limitadas las mismas por la falta de continuidad y oferta de los productos; por tal razón, se hizo necesario realizar un análisis exploratorio de datos en un periodo comprendido entre el año 1975 y 2014 a la estación hidrometeorológica Puente El Yopal, ubicada en la ciudad de El Yopal-Casanare sector del MORRO - PUENTE LA CABUYA específicamente en la cuenca media del río debido a que en esta se encuentra situada la Ciudad de Yopal Casanare, siendo que es una de las más beneficiadas por este recurso de acuerdo a los datos recopilados por la empresa Enviromental Ingenieros Consultores en la ejecución de contrato de consultoría suscrito con la Corporación Autónoma Regional de la Orinoquia CORPORINOQUIA en el año 2007 con el fin de desarrollar el Plan de Ordenación Y Manejo de la Cuenca del Río Cravo Sur POMCH, allí se especifica que la cuenca es un territorio de estructura rural el cual, del 100% de su área, solo el 0,26% alberga los cascos urbanos de los municipios de Labranzagrande, Paya, Pisba, Nunchía y Yopal, donde este último, corresponde al 94.3% del territorio urbano de la cuenca.

Se pretende pronosticar el caudal del río 15 años y obtener elementos que sirvan de base para la toma de decisiones con el fin de prever y minimizar los riesgos derivados de variaciones extremas en los caudales. En el orden de las ideas anteriores, este documento servirá como herramienta de planeación del territorio.

La situación existente asociada al problema o necesidad, es la presencia de un alto nivel de incertidumbre referente a las condiciones del Recurso Hídrico de las cuencas, específicamente la del rio Cravo Sur, información que viene siendo recopilada por entidades gubernamentales como: Ideam, Ingeominas, IGAC, Gobernaciones, Alcaldías, entre otras.

Algunas características generales de la Cuenca son: el Río Cravo Sur, desde su nacimiento en la cordillera Oriental en territorio del Departamento de Boyacá, específicamente el municipio de Tasco a 3.800 m.s.n.m. entre la Serranía de Peña Negra y el Páramo de Cadilla el cual, se encuentra protegido dentro de los límites del Parque Nacional Natural de Pisba, hasta su desembocadura en el Río Meta, municipio de Orocué Casanare a los 150 m.s.n.m. presenta una importante corriente que baña los municipios de Tasco, Socotá, Gámeza, Mongua, Aquitania, Labranzagrande, Pisba y Paya por el Departamento de Boyacá perteneciendo a este un 35.47% de su extensión total y Támara, Nunchía, Yopal, San Luis de palenque y Orocué por el Departamento de Casanare con el 64.53% restante; su área total es de 565.113 hectáreas y su forma es oval oblonda a rectangular oblonda, alberga una corriente principal de 205 kilómetros de longitud que transcurre sobre un amplio gradiente altitudinal (CORPORINIQUIA, CORPOBOYACA, & UAESPNN, 2007).

Con el propósito principal de generar un documento que presenten elementos que sirvan de base para la toma de decisiones con el fin de prever y minimizar los riesgos derivados de variaciones extremas en los caudales y a su vez sirva como herramienta de planeación del territorio, se desarrolla una investigación con enfoque cuantitativo debido a que se basa en el análisis exploratorio de una serie de datos históricos con un alcance descriptivo ya que, se analiza información de manera independiente y teniendo en cuenta el siguiente diseño metodológico:

Localización de la zona de estudio: la cuenca se localiza geográficamente entre los 4°41'13" hasta los 5°56'37" de Latitud Norte y entre los 71°34'09" y 72°46'28" de Longitud al Oeste de Greenwich, en conjunto con la influencia que ejerce la Zona de Convergencia Intertropical (ZCIT), condicionan las propiedades físicas de la atmósfera que definen el tiempo y el clima. Actúan también sobre las características climáticas de la cuenca del río Cravo Sur, todas las condiciones físicas que le son propias (forma, orientación, rango de altitudes, topografía, vegetación, etc.) y que producen o modifican los elementos que crean el clima (CORPORINIQUIA, CORPOBOYACA, & UAESPNN, 2007).

La investigación se desarrolla con los registros históricos de caudales desde el año 1975 hasta el año 2014 provenientes de la estación hidrometeorológica Puente el Yopal operada por el Instituto de Hidrología, Meteorología y Estudios Ambientales (IDEAM) ubicada en el sector comprendido entre el corregimiento El Morro y el Puente La Cabuya el cual, une al municipio de Yopal con el Oriente del llano Casanareño, datos que fueron facilitados por este ente gubernamental para fines académicos.

Fase I. Revisión y control de calidad de la serie de datos históricos y análisis exploratorio de datos: el IDEAM, como institución estatal encargada de "Fortalecer la capacidad tecnológica, científica, administrativa y financiera para producir la información hidrológica, meteorológica y ambiental de manera oportuna y con la calidad que requieran la ciudadanía, los sectores públicos y privados del país" (IDEAM, 2014,) fue la entidad facilitadora de los datos históricos de caudales existentes de los últimos 40 años recopilados mediante dos estaciones hidrometeorológicas activas actualmente, ubicadas en el municipio de Yopal, Departamento de Casanare y descritas a continuación:

Recibidos y tabulados los datos históricos existentes de caudales de las dos estaciones más cercanas a la ciudad de El Yopal descritas anteriormente, logra evidenciarse que no se cuenta con series completas por lo cual, se realiza proceso de imputación de datos para cada una de ellas, optando por completar por el método de los promedios y obteniendo como resultado dos series completas con 480 datos independientes por cada estación.

Se realiza el análisis exploratorio de datos, donde se obtuvo los estadísticos descriptivos básicos, que permitieron realizar una evaluación preliminar de las dos series, teniendo en cuenta la media aritmética, media recortada al 5%, mediana de cada estación, además de la dispersión y asimetría. Dicho análisis permitió también revisar la existencia o indicios de normalidad.

Fase II. Selección de estación a pronosticar: Para determinar si lo conveniente para la investigación era trabajar con la serie de datos de una de las dos estaciones o si por lo contrario se podría promediarlas y pronosticar los caudales con una sola serie se procede a evaluar si existen diferencias significativas en dichas series; previamente se evalúa el cumplimiento del supuesto de normalidad para determinar si se emplearían pruebas paramétricas o no paramétricas; para ello, se emplea la prueba de Kolmogorov - Smirnov rechazándose la hipótesis de que existe presencia de normalidad en ambas series debido a que su significancia se encuentra por debajo del 5% motivo por el cual, se aplican pruebas no paramétricas para comparar las series y evaluar las diferencias significativas entre ellas.

Al no contar con presencia de normalidad en las series, se procede a comparar visualmente entre ellas las medianas dentro de los gráficos de caja, los cuales, presentan pocos valores extremos y se evidencia una mediana mayor en la estación Puente El Yopal.

Basados en los resultados, se aplica la prueba no paramétrica de Mann Whitney evidenciándose diferencias significativas entre las medianas de las series debido a que su significación se encuentra por debajo del 5%; por lo anterior y de acuerdo al criterio del investigador, se toma la decisión de realizar la pronostico final con la serie de la estación hidrometeorológica Puente El Yopal debido a que esta cuenta con caudales mayores y con el fin de no subestimar el pronóstico.

Este gráfico presenta de manera preliminar que no existe tendencia en la serie por lo cual de manera intuitiva se afirmar que existe estacionariedad en media y en varianza. Por otra parte, se observan picos que evidencian estacionalidad la cual, es una característica propia de la serie y se corrobora mediante los factores estacionales relacionados en la tabla No 6.

Con el fin de corroborar la existencia de estacionariedad, se realiza análisis al grafico No 5 el cual, muestra que, al no decrecer lentamente, la estacionariedad está presente en la serie.

Utilizando el módulo para pronósticos del Software IBM SPSS Statistics el cual, "es el principal software estadístico que ofrece técnicas de recolección de datos y analítica predictiva para solucionar múltiples problemas empresariales y de investigación. Brinda varias técnicas, que incluyen pruebas de hipótesis lo que facilita la gestión de los datos, la selección y la ejecución de análisis y el intercambio de resultados, así como análisis multivariados y de sobrevida" (Statistics, IBM COLOMBIA S.A., 2018) y de acuerdo a las características de la serie, se obtiene un Modelo Estacional Autoregresivo Integrado de Medias Móviles ARIMA (p,d,q)(P,D,Q).

De acuerdo a la tabla No 8, se evidencia que, al evaluar la independencia de los residuales, estos no se encuentran correlacionados dado que son independientes indicándo que el modelo de pronostico seleccionado es robusto.

Realizado el pronóstico, se realiza prueba de Mann Whitney con el fin de determinar si existen diferencias significativas entre los datos iniciales y los pronosticados, evidenciándose que debido a que su significación se encuentra por encima del 5%, no existen diferencias significativas entre las series, por lo tanto, se puede afirmar que el modelo utilizado es eficiente para el pronóstico.

Con el pronostico final se procede a realizar análisis exploratorio de datos para cada mes, donde se obtuvo los estadísticos descriptivos básicos, que permitieron evaluar preliminarmente la serie pronosticada identificando el promedio del caudal el cual, tiende a aumentar anualmente, mientras que, promediando el caudal de cada mes durante los quince años pronosticados, los meses con mayores caudales tienden a ser los del tercer trimestre de cada año. Tabla 10 (Anexo).

Los datos aquí presentados servirán como fundamento teórico para la comunidad en general y a entes gubernamentales quienes, de acuerdo a su criterio, puedan tomar decisiones en pro del desarrollo de la región teniendo en cuenta, que son datos de caudales que posiblemente puedan ocurrir en la región y que podrían o no afectar los sectores afines a esta área como lo son el agro y ganadero además de la calidad de vida de la población del municipio El Yopal principalmente.

Se observa que el comportamiento del caudal en la cuenca media del Río Cravo Sur es uniforme mostrando ciclos anuales de estacionalidad presentados en las estaciones hidrometeorológicas ubicadas en esta zona de estudio.

En el análisis exploratorio realizado a los datos iniciales se obtuvo información descriptiva suficientemente útil sobre el comportamiento general de los datos, los cuales fueron corroborados con los resultados obtenidos en los análisis confirmatorios; cabe destacar también que este es un paso necesario para realizar cualquier tipo de investigación ya que, permite tener una visión descriptiva de los datos con ayuda de los estadísticos básicos, los gráficos y pruebas de normalidad.

Realizado el análisis del comportamiento de las series históricas, logra evidenciarse que el caudal en algunos sectores de la zona de estudio se encuentra por encima de los caudales aguas arriba y que este aumenta de manera importante principalmente en el sector circundante al municipio de El Yopal en el Departamento de Casanare.

De acuerdo con el pronóstico realizado, el promedio del caudal en la zona media del río Cravo Sur comprendida entre EL CORREGIMIENTO DEL MORRO - PUENTE LA CABUYA tiende a aumentar anualmente, mientras que, promediando el caudal de cada mes durante los quince años proyectados, los meses con mayores caudales tienden a ser los del tercer trimestre de cada año.

Siempre debe buscarse la posibilidad de realizar análisis con series de datos de tiempo completas con el fin de corroborar la veracidad en el comportamiento de cada variable.

Es claro que se presenta un pronóstico estimados de caudales derivados de datos históricos hasta el año 2014 pero ello no implica que los datos obtenidos sean cien por ciento reales por lo cual, podría ocurrir un evento anormal de caudal en la zona llegando a niveles no vistos antes debido a la constante variación del clima en los últimos años, es recomendable que la comunidad en general solicite a las entidades locales (alcaldías, corporaciones regionales ambientales entre otras), estudios más especializados con el fin de prevenir eventos no deseados.

Es importante continuar realizando investigaciones de este tipo, incorporando al análisis otras variables como lo son: temperatura, humedad, precipitación, velocidad del viento, evaporación, evapotranspiración entre otros, con el fin de evaluar el efecto de las mismas sobre los niveles del caudal.