INTRODUCCIÓN
El monitoreo de los eventos hidrológicos presenta una limitación en la data disponible para la modelación computacional de los eventos y comprender su complejidad (Muñoz et al., 2023). Los datos son proporcionados por las estaciones meteorológicas que son de representación puntual, pero en zonas con una baja densidad de estaciones meteorológicas es necesario buscar alternativa que brinden solución a la problemática de escasa data a nivel espacial, surgen métodos de interpolación, reanálisis y validación de datos satelitales (Essou et al., 2017; C.M et al., 2022; Vaheddoost et al., 2023).
Las cuencas, sub cuencas y quebradas sin aforo requieren de valores que aproximen a determinar la disponibilidad hídrica para prevenir la escasez de agua, la cual es predominante en regiones áridas y semiáridas, el planteamiento de soluciones para una correcta gestión evitaría conflictos (Marco-Dos Santos et al., 2019). También es necesario para una gestión ante el riesgo de áreas con sensibilidad a las inundaciones (Nieto-Chaupis, 2019).
En el Perú se han realizado estudios para construcción de datos de lluvias grillados, de alta resolución para las cuencas hidrográficas, con ello se evalúan balances hídricos, eventos extremos, inundaciones, modelos hidrológicos con el fin de respaldar las políticas y prácticas del agua (Zubieta et al., 2015; Aybar et al., 2019; Fernandez-Palomino et al., 2022).
En Lima se realizó asociación de la variación de almacenamiento de agua entre correlación de dos productos grillado, Tropical Rainfall Measuring Mission (TRMM) y Gravity Recovery and Climate Experiment (GRACE), se obtuvo 0,8677 como mejor resultado (Mendoza et al., 2022).
En Áncash se buscó mejorar la correlación entre los productos TRMM y Global Precipitation Measurement (GPM), los resultados son aceptables (Villavicencio et al., 2022). La cuenca del río Piura se realizó una validación de los datos diarios del satélite TRMM para obtener valores extremos y aplícalos a la ingeniería hidráulica (Cabrera et al., 2016).
Es necesario realizar validaciones de alternativas de productos grillados de precipitación para respaldar estudios posteriores, ya que en Tumbes no se ha realizado. Es por ello que se busca validar el Producto Pisco v2.1 mensual, por estación meteorológica, elevación de la cuenca y región climática. Cabe destacar que el producto Pisco v2.1 es importante porqué contempla un proceso especializado de métodos combinación de datos de estaciones terrenas, reanálisis de productos grillados, un estricto proceso de control de calidad y métodos geoestadísticos.
MATERIAL Y MÉTODOS
El tipo de investigación es retrospectivo y transversal para un análisis en un solo corte y de datos pasados (Dagnino, 2014; Calderón & Alzamora de los Godos, 2018). La metodología realizada se describe en la Figura 1.
1. Obtención de Datos
Los datos de precipitaciones estimadas por el producto grillado PISCO (Peruvian Interpolated data of the SENAMHI’s Climatological and hydrological Observations) por sus siglas en inglés (Aybar et al., 2019). La validación fue en el mismo periodo de tiempo (36 años), el producto se encuentra alojado en la siguiente dirección web: https://iridl.ldeo.columbia.edu/SOURCES/.SENAMHI/.HSR/.PISCO/.Prec/.v2p1/ donde se puede descargar los datos estimado de precipitaciones mensuales del producto PISCO precipitaciones mensuales (PISCOpm_v.2,1). Por ser un producto grillado se presenta en cuadriculas para su ello se utiliza software R y R Core Team (2022), en la extracción de los datos se realiza de manera puntual en la Tabla 1 se detalla las estaciones meteorológicas correspondientes.
N | Nombre | Latitud | Longitud | Alti-tud |
1 | Puerto Pizarro | 3°30'13.67'' | 80°23'42.14'' | 6 |
2 | Los Cedros | 3° 38' 1'' | 80° 32' 1'' | 74 |
3 | Campamento Sede | 3°33'4.14'' | 80°25'46.64'' | 17 |
4 | Quebrada Angostura | 3° 46' 1'' | 80° 24' 1'' | 89 |
5 | El Tigre | 3°48'8.4'' | 80°27'24'' | 45 |
6 | Rica Playa | 3°48'42.84'' | 80°30'19.65'' | 68 |
7 | Cabo Inga | 3°58'33.39'' | 80°24'6.54'' | 231 |
8 | Capitán Hoyle | 4° 4' 1'' | 80° 29' 1'' | 235 |
9 | Condor Flores | 3° 53' 24.09'' | 80° 8' 22.73'' | 384 |
Fuente: adaptado de SENAMHI.
2.Preparativos y control de calidad de datos de precipitación
2.1. Análisis exploratorio data
Para el análisis se realizó sobre data con ordenamiento multianual y mensual para cada una de las estaciones meteorológica, se realizó una exploración visual, tendencia y saltos (Kundzewicz & Robson, 2000).
En la identificación de datos atipicos se empleó el método de Turkey, para la completación de datos se empleó el método de software HEC- 4 Monthly Streamflow Simulation, el cual cuenta con un respaldado en variedad de estudios hidrológicos.
3.Validación
La validación de las datas de precipitación mensuales del producto Pisco v.2,1 y la data de las estaciones meteorológicas del Servicio Nacional de Metereológia e Hidrología del Perú (SENAMHI). se empleó métricas estadísticas y la correlación entre ambas data.
3.1. Validación por estación meteorológica
El producto PISCO tiene una resolución espacial de 0,1°x 0,1° (10km2x10km2), la información puntual de las estaciones meteorológicas y las observadas por SENAMHI en series de tiempo desde enero de 1981 hasta diciembre de 2016 (dependiendo de la información disponible). La validación se realiza entre las series de tiempos con las medidas de comparación estadísticas como: Coeficiente de determinación(r2), error absoluto medio (MAE), raíz error medio cuadrático (RMSE), índice de eficiencia de nash-sutcliffe (NSE), sesgo estadístico porcentual (PBIAS).
3.2. Validación por elevación
Se adecua las series de precipitación en Precipitación Promedio Multianual (PPM), esto para la data estimada y observada. También se considera su altitud de cada una de ella, para ello se debe ordenar en columnas con los siguientes titulados: Estación, Altitud y Precipitación Promedio Multianual Observada (mm) y Precipitación Promedio Multianual Estimada PISCO (mm).
La correlación se realiza entre la altitud de las estaciones y su Precipitación Promedio Multianual (mm) observada, de igual manera la Precipitación Promedio Multianual Estimada PISCO (mm) se emplea el coeficiente de determinación, esto nos indicará el comportamiento y tendencia de los datos de precipitación y su altitud, el ajuste entre observadas y estimadas.
3.2. Validación por estacionalidad
Las estaciones del año están dadas por el movimiento de traslación de nuestro planeta, tiene una duración de 3 meses, para el Perú el otoño comienza en marzo, primavera inicia en septiembre, invierno en junio y verano en diciembre (Gob.pe, 2022). Teniendo en consideración lo anterior se genera 4 estaciones de acuerdo con sus meses se agrupará los datos de precipitación mensuales así:
Otoño (marzo-abril-mayo), invierno (junio-Julio-agosto), primavera (septiembre-octubre-noviembre) y verano (diciembre-enero-febrero). En la serie de datos se verá reflejado en las precipitaciones totales mensuales (mm), iniciando desde marzo 1981, abril 1981, mayo 1981 que es otoño de igual manera se agrupará en función del mes y la estación que corresponde hasta llegar al año de 2016, esto se realizará tanto para precipitación estimada y observada. se realizará la validación con los estadísticos de comparación entre las series estacionales de precipitaciones estimadas y observadas.
3.3. Validación por región climática
Para la validación se debe tener el mapa de regiones climáticas del Perú que se encuentra en https://www.senamhi.gob.pe/?p=mapa-climatico-del-peru el cual fue elaborado por el servicio nacional de meteorología e hidrología del Perú (SENAMHI), con ello se identifica la regiones climáticas teniendo en cuenta la precipitación anual de cada una de las estaciones con la data observada o estimada de precipitación, para ello se realiza series de tiempos de precipitación anual promediadas desde enero de 1981 a diciembre de 2016. luego las series de tiempo según a la región climática que corresponda se realiza la validación entre la data observada y estimada de acuerdo con los estadísticos comparativos.
3.4. Medidas estadísticas de comparación.
Para la validación del comportamiento de los datos de precipitación observada y estimada, aplican estadísticos comparativos en la Tabla 2 se detalla cada una de ellas como: coeficiente de determinación(r2), error absoluto medio (MAE), raíz error medio cuadrático (RMSE), índice de eficiencia de nash-sutcliffe (NSE), sesgo estadístico porcentual (PBIAS). Los resultados de los estadísticos comparativos permiten clasificar y analizar su comportamiento de acuerdo con algunos rangos establecidos.
Métrica Estadística | Unidad | Ecuación | Valor óptimo |
---|---|---|---|
Error Absoluto Medio | m m/h | 1 n 𝒊=𝟏 𝒏 𝑬𝒊−??𝒊 | 0 |
Raíz del Error Medio Cuadrático | mm/h | 𝟏 𝒏 𝒊=𝟏 𝒏 𝑬𝒊−𝑶𝒊 𝟐 | 0 |
Coeficiente de Determinación | ------- | 𝑟 𝑠 2 =1− 6 𝑑 2 𝑛 𝑛 2 −1 | ±1 |
Índice de Eficiencia de Nash-Sutcliffe | ------- | 1− 𝒊=𝟏 𝒏 𝑬𝒊−𝑶𝒊 𝟐 𝒊=𝟏 𝒏 𝑶𝒊− 𝑶 𝟐 | 1 |
Sesgo Estadístico Porcentual | ------- | 𝟏𝟎𝟎∗ 𝒊=𝟏 𝒏 𝑬𝒊−𝑶𝒊 𝒊=𝟏 𝒏 𝑶𝒊 | < ±10 |
Nota: n número de muestras, Ei Datos Pisco, Oi Datos Observados.
RESULTADOS Y DISCUSIÓN
Análisis exploratorio data
Se realizó el análisis de datos atípicos con el método de box-plot, con ello se pudo determinar los datos que no se ajustan al comportamiento normal de los conjuntos de datos (Figura 2). Se encontraron en promedio 3 valores atípicos, los cuales se asocian a los eventos del Niño costero en los años 1983, 1992, 1997, 1998. Consultando y comparando los datos con estaciones vecinas se decidió conservar esos datos, puesto que eliminarlos sería eliminar los rasgos históricos y negar su existencia Para el caso de las variables de precitación del producto PISCO no necesita mayor análisis debido a que esta ya tiene una compleja validación y análisis a profundidad.
Para el análisis de consistencia de los datos se realizó con la curva de doble masa, para el análisis estadístico se realizó consistencia en la media y en la desviación estándar. Todo esto se llevó a cabo con el software Hydrognomon 4 y Trend.
La completación de datos se hizo con Hec- 4 Monthly Streamflow Simulation el cual está basado en un lenguaje fortran. Es sencillo y ofrece buenos resultados, además tienen estudios y experiencia que lo respaldan.
Validación por estación meteorológica
El análisis puntual de los datos considerando la variabilidad espacial del producto PISCO se realizó considerando las estaciones meteorológicas y su disponibilidad de data.
Para ello se consideró las métricas estadísticas, se comenzó realizando la lectura de los datos observados y simulados en Rstudio para realizar la comparación con el package “hidroGOF” (Mauricio Zambrano-Bigiarini, 2020), quien cuenta con variedades de métricas. Se observó que en la métrica de Coeficiente de determinación (R2) que las estaciones como Campamento Sede, Tigre y Puerto Pizarro presentaron los máximos valores 0,89, 0,84 y 0,85 respectivamente considerado una correlación Muy Alta según (Mamani, 2019). En contrario la estación Condor Flores tiene el mínimo valor de 0,12 y las estaciones Cedros, Quebrada Angostura y Rica Playa una correlación Alta con valores de 0,61, 0,62 y 0,73. Cabo inga y Capitán Hoyle presentan valores de 0,46 y 0,54 con una clasificación Media.
Entonces se puede inferir a partir de los valores obtenidos que el 66,67% de las estaciones presentan un ajuste bueno para tener en cuenta.
Otro análisis de importancia es el Índice de Eficiencia de Nash- Sutcliffe (NSE), nos indica que tan preciso es el modelo generado, se encontraron solo dos estaciones que presenta una precisión Insuficiente que son Condor Flores y Capitán Hoyle, luego la estación Cabo Inga con 0,32 tiene una precisión satisfactoria. Una excelente precisión presenta las estaciones como Campamento Sede y Puerto Pizarro con valores superiores 0,80.
En el análisis el Sesgo Estadístico Porcentual (PBIAS), representa la sobreestimación en caso de valores positivos y con valores contrarios una subestimación (Moriasi, Arnold, Van Liew, & Harmel, 2007). Se tiene un 55,56% de las estaciones totales con sobreestimación, las cuales van desde valores de 5,4 a 49,5 y son Quebrada Angostura, Capitán Hoyle respetivamente.
Con subestimación se tiene los valores extremos de -5,5 a -30,8 que corresponde a Cabo Inga y Rica Playa. Así como se resumen los valores en la Tabla 3.
Estación | R^2 | MAE | RMSE | NSE | PBIAS % |
Puerto Pizarro | 0,46 | 38,9 | 138,89 | 0,32 | -5,5 |
Los Cedros | 0,89 | 11,52 | 37,4 | 0,88 | -13,9 |
Campamento Sede | 0,54 | 41,79 | 126,51 | -0,19 | 49,5 |
Quebrada Angostura | 0,61 | 14,81 | 57,62 | 0,6 | 8,6 |
El Tigre | 0,12 | 78,95 | 256,36 | -0,09 | -26 |
Rica Playa | 0,85 | 12,35 | 39,42 | 0,84 | 9,8 |
Cabo Inga | 0,62 | 20,46 | 62,32 | 0,53 | 5,4 |
Capitán Hoyle | 0,73 | 22,77 | 68,34 | 0,7 | -30,8 |
Condor Flores | 0,84 | 14,35 | 51,79 | 0,75 | 6,4 |
Validación por elevación
Se realizó un clúster por elevación de las estaciones para la comprobación de los datos, como resultado se obtuvo tres clasificaciones, Cuenca Alta (Condor Flores, Cabo Inga, Capitán Hoyle), Cuenca Media (Tigre, Rica Playa, Quebrada Angostura) y Cuenca Baja (Campamento Sede, Cedro, Puerto Pizarro). Para ello se consideró los promedios de las estaciones por cada clasificación.
De acuerdo a la Tabla 4 se encontró los siguientes resultados:
Para el análisis de Coeficiente de Determinación (R2), se encontró los mejores resultados en la Cuenca Alta con un valor de 0,57. El menor valor obtuvo la Cuenca Media que es de 0,22, al igual con un valor cercano la Cuenca Baja con 0,27.
Con el Índice de Eficiencia de Nash- Sutcliffe (NSE), que determina que tan preciso es el modelo.
La Cuenca Alta, Cuenca Media y Cuenca Baja se obtuvo resultados insuficientes. Con valores de 0,06, -0,45 y -9,27 respectivamente.
La variación en la precisión de los valores del modelo se puede dar por sobreestimación o subestimación, para ello se empleó el Sesgo Estadístico Porcentual (PBIAS), las Cuenca Baja tiene una sobreestimación con un valor de 7,1. En contrario, las Cuenca Media y Alta tienen una subestimación con valores de -11,9 y -25.9 respectivamente.
Se puede inferir que, en función de la elevación y el agrupamiento de estas estaciones, en la parte baja presenta una baja relación lineal de los datos, insuficiente precisión y una sobreestimación. En la parte media y alta presenta una subestimación.
Validación por estacionalidad
Se consideró el agrupamiento de los datos de precipitación en arreglo de acuerdo los meses de las estaciones de año, para cada estación tanto observada y simulada.
En verano que corresponde a los meses de diciembre-enero-febrero. En la Tabla 5 se tiene el resumen de los resultados obtenidos, pero se detalla que en la métrica del Coeficiente de Determinación (R2), las estaciones como Rica Playa con 0,96 alcanzo el máximo valor de ajuste de los datos, siguiéndole Campamento Sede, Puerto Pizarro y Tigre. En contra parte, la estación Condor Flores con 0,26 tiene la peor varianza o ajuste de los datos.
En términos de Índice de Eficiencia de Nash- Sutcliffe (NSE), las estaciones Campamento Sede y Puerto Pizarro presentan el máximo ajuste con 0,9, seguida por Tigre, Rica Playa y Cedros. Las demás estaciones presentan un ajuste de satisfactorio a insuficiente, teniendo en cuenta a la estación Quebrada Angostura con el peor ajuste.
Para términos del Sesgo Estadístico Porcentual (PBIAS), se tiene más subestimación en las estaciones, la estación Rica Playa presenta -42,7 y Puerto Pizarro con -3,3. La sobreestimación se da en las estaciones Capitán Hoyle (39,6), Tigre (9,1) y Quebrada Angostura (3,5).
Estación | R^2 | MAE | RMSE | NSE | PBIAS % |
Cabo Inga | 0,32 | 76,68 | 223,49 | 0,24 | -10,9 |
Campamento Sede | 0,9 | 18 | 45,98 | 0,9 | -11,7 |
Capitán Hoyle | 0,43 | 64,1 | 156,56 | -0,08 | 39,6 |
Cedros | 0,67 | 26,84 | 83,95 | 0,67 | -4,9 |
Condor Flores | 0,26 | 95,74 | 156,33 | -0,16 | -21,8 |
Puerto Pizarro | 0,91 | 18,62 | 42,3 | 0,9 | -3,3 |
Quebrada Angostura | 0,37 | 44,15 | 99,98 | -0,24 | 3,5 |
Rica Playa | 0,96 | 46,01 | 92,07 | 0,74 | -42,7 |
Tigre | 0,89 | 23,79 | 56,35 | 0,82 | 9,1 |
El invierno que se desarrolla entre los meses de junio, julio y agosto. Como se muestra en la Tabla 6 se presenta un muy alto ajuste de los datos en 3 estaciones de 9 las cuales son Campamento Sede, Puerto Pizarro y Rica Playa con valores superiores al 0,90, La estación Tigre con 0,61 tiene un ajuste alto y el resto de estaciones que son Cabo Inga, Capitán Hoyle, Cedros, Quebrada Angostura y Condor flores tienen un valor de 0 con respecto a esto (Mamani,2019) lo clasifica como mínima.
Con respecto al Índice de Eficiencia de Nash- Sutcliffe (NSE), se presenta un desempeño Bueno en las estaciones Puerto Pizarro, Rica Playa y Tigre con un rango de 0,41 - 0,53. En las demás estaciones se evidenció un desempeño insuficiente en la clasificación de (Molnar, 2011).
En desempeño de Sesgo Estadístico Porcentual (PBIAS), las estaciones como Condor Flores, Cabo Inga presentan los menores valores, pero igual junto con las demás estaciones tienen un desempeño insatisfactorio. Pero en cuestiones de sobreestimación tenemos a 4 de 9 estaciones que son: Capitán Hoyle, Cedros, Quebrada Angostura y Condor Flores. Los otros 5 restantes presentan subestimación.
Estación | R^2 | MAE | RMSE | NSE | PBIAS % |
Cabo Inga | 0 | 7,23 | 34,38 | -1,36 | -32,9 |
Campamento Sede | 0,98 | 5,4 | 37,77 | 0,29 | -73,9 |
Capitán Hoyle | 0 | 5,1 | 30,54 | -1148,59 | 2755,5 |
Cedros | 0 | 1,43 | 5,35 | -10,75 | 74,1 |
Condor Flores | 0 | 5,66 | 18,33 | -2,48 | 10,5 |
Puerto Pizarro | 0,98 | 3,95 | 25,21 | 0,47 | -58,5 |
Quebrada Angostura | 0 | 4,19 | 21,46 | -431,16 | 923 |
Rica Playa | 0,94 | 3,78 | 21,28 | 0,51 | -71,3 |
Tigre | 0,61 | 5,34 | 30,12 | 0,53 | -58,1 |
Los meses de marzo, abril y mayo, periodo donde se da la estación de otoño, las precipitaciones al inicio se presentan en periodo lluvioso y luego se llega al periodo seco. Teniendo en cuenta a la Tabla 7 con respecto a las métricas se detalla que el Coeficiente de Determinación (R2), las estaciones Campamento Sede, Puerto Pizarro presenta el mejor ajuste teniendo una clasificación de 0,92 y 0,91 que según (Mamani, 2019) es muy alta. Luego la estación Condor Flores con 0,05 presenta un ajuste con clasificación mínima y Cabo Inga tiene una clasificación de media junto a las demás estaciones.
Con el Sesgo Estadístico Porcentual (PBIAS), se pudo constatar que la estación Campamento Sede tiene una subestimación con una clasificación de buena, como la estación Tigre con una sobreestimación, En subestimación las estaciones tienen un comportamiento insatisfactorio con valores promedio a 25 y cabe mencionar a la estación Capitán Hoyle que presenta el máximo valores de sobreestimación con desempeño insatisfactorio según (Moriasi, Arnold, Van Liew, & Harmel, 2007) con un valor de 62,5.
Estación | R^2 | MAE | RMSE | NSE | PBIAS % |
Cabo Inga | 0,43 | 72,93 | 187,11 | 0,9 | 4,3 |
Campamento Sede | 0,92 | 20,95 | 44,74 | 0,88 | -12,2 |
Capitán Hoyle | 0,57 | 91,41 | 195,05 | -0,64 | 62,5 |
Cedros | 0,51 | 29,53 | 78,72 | 0,32 | 25,8 |
Condor Flores | 0,05 | 178,73 | 472,11 | -0,27 | -24,9 |
Puerto Pizarro | 0,91 | 25,83 | 61,51 | 0,75 | 28,3 |
Quebrada Angostura | 0,74 | 44,32 | 88,77 | 0,71 | -29,8 |
Rica Playa | 0,54 | 39,55 | 98,55 | 0,46 | -14,7 |
Tigre | 0,78 | 25,65 | 80,87 | 0,61 | 11,5 |
La estimación con la métrica de Índice de Eficiencia de Nash- Sutcliffe (NSE), la variación del ajuste en clasificación excelente con valor de 0,9 y 0,88 son las estaciones Cabo Inga y Campamento Sede. Condor Flores presenta un ajuste insuficiente teniendo en cuenta que cuando los valores son negativos se entiende que la media de los valores observados tiene mejor ajuste que los valores estimados.
De acuerdo la Tabla 8 que corresponde a la estación primavera considerando los meses de setiembre, octubre y noviembre. Se tiene indicando que el Coeficiente de Determinación (R2), presenta un excelente ajuste en estaciones como Campamento Sede con 0,96. De igual manera Puerto Pizarro con 0,98. También se tiene las estaciones Cedros y Rica Playa. Existen una proporción de variación mayor que genera un ajuste mínimo en las estaciones de Condor Flores, Quebrada Angostura.
Para el Índice de Eficiencia de Nash- Sutcliffe (NSE), se evidencia un buen desempeño en las estaciones como Campamento Sede, Cedros, Puerto Pizarro con valores desde 0,88 a 0,97 con una clasificación excelente bajo la clasificación de (Molnar,2011) también se tiene estaciones que tienen un desempeño insuficiente o de pésimo ajuste como son las estaciones Capitán Hoyle, Condor Flores y Quebrada Angostura.
También en la Tabla 8 se encuentran valores para Sesgo Estadístico Porcentual (PBIAS), el desempeño en esta métrica para la estación Puerto Pizarro es el mejor teniendo una mínima sobre-estimación con el valor de 4,6. Para las demás estaciones se tiene subestimaciones y sobreestimaciones con un desempeño insatisfactorio.
Estación | R^2 | MAE | RMSE | NSE | PBIAS % |
Cabo Inga | 0,77 | 2,78 | 9,48 | 0,6 | -41,8 |
Campamento Sede | 0,96 | 1,74 | 7,3 | 0,88 | 27,7 |
Capitán Hoyle | 0,25 | 6,54 | 23,07 | 0,02 | -75,8 |
Cedros | 0,93 | 1,45 | 2,8 | 0,91 | -33 |
Condor Flores | 0,07 | 35,66 | 123,37 | -0,13 | -54,1 |
Puerto Pizarro | 0,98 | 1,01 | 2,62 | 0,97 | 4,6 |
Quebrada Angostura | 0,01 | 6,04 | 23,85 | -6,02 | 57,6 |
Rica Playa | 0,94 | 1,75 | 6,47 | 0,68 | -60,8 |
Tigre | 0,42 | 2,62 | 10,28 | 0,35 | -36,1 |
Validación por región climática
Para la zona de estudio se encontraron 3 regiones climáticas con estaciones que nos brindan los datos para realizar la comparación de ajuste entre los valores observados y simulados.
En la Tabla 9 se puede observar que la región climática C(i,p)A’ se define como una región de precipitación efectiva semiseco cálido y con invierno, primavera seco que se encuentra a altitud mayor a 600 m.s.n.m. En Tumbes se encuentra en áreas de los cerros de Amotape con temperaturas máximas de 21 °C a 31° C y con lluvias anuales de 700 y 1200 mm. Para el primer análisis de desempeño en la métrica de Coeficiente de Determinación (R2), se encuentra un valor de 0,41 considerado un ajuste bueno y un valor de 146,82 para Raíz del Error Medio Cuadrático (RMSE).
Según Índice de Eficiencia de Nash- Sutcliffe (NSE), se tiene un valor de 0,23 se tiene un ajuste satisfactorio. Con el Sesgo Estadístico Porcentual (PBIAS) se determinó que existe un muy buen ajusto con una mínima subestimación teniendo un valor de -4,2.
La región climática D(i)A’ se define como semiárido con invierno seco cálido (Tabla 9). Este se desarrolla en altitudes mayores a los 300 m.s.n.m. en los Cerros de Amotape de Tumbes, caracteriza por tener precipitaciones anuales entre 700 a 1200 mm y temperatura de 17 °C a 29 °C.
Se obtuvo Coeficiente de Determinación (R2) de 0,12 (muy bajo) y una Raíz del Error Medio Cuadrático (RMSE) de 256,36. Para el Índice de Eficiencia de Nash- Sutcliffe (NSE) el modelo presenta un ajuste insatisfactorio con un valor de -0,09 considerando que la media de los observados tiene mejor ajuste a los valores estimados y como final tenemos el Sesgo Estadístico Porcentual (PBIAS) con el cual se encontró una subestimación insatisfactoria con un valor de -26.
En la región climática D(i,p)A’ un clima semiárido, cálido con invierno y primavera seco es el clima que se encuentra en gran ocupación en el área de estudio presentando temperaturas de 19 °C a 33 °C y precipitaciones anuales de 200 a 500 mm. Además, se da en altitudes menores a los 300 m.s.n.m. Teniendo en cuenta la métrica de Coeficiente de Determinación (R2) se tiene un ajuste muy alto con un valor 0,92 que nos indica que la variación es la mínima, para la Raíz del Error Medio Cuadrático (RMSE) se tiene un valor de 28,94 un error bajo. Tomado en cuenta el valor de 0,91 para el Índice de Eficiencia de Nash- Sutcliffe (NSE) se tiene una clasificación de excelente Estadístico Porcentual (PBIAS) se una mínima subestimación con una clasificación de muy buena con un valor de -3,6, esta estación presenta un ajuste muy bueno en las métricas.
Región Climática | R^2 | MAE | RMSE | NSE | PBIAS % |
C(i,p)A' | 0,41 | 39,91 | 146,82 | 0,23 | -4,2 |
D(i)A' | 0,12 | 78,95 | 256,36 | -0,09 | -26 |
D(i,p)A' | 0,92 | 8,56 | 28,94 | 0,91 | -3,6 |
Para Mamani (2019) que realizó una regionalización de precipitaciones mensuales en la Cuenca de la Vertiente del Lago Titicaca - Lado peruano, se hizo validación de los datos del Producto PISCO teniendo como indicadores el Coeficiente de Determinación (R2), para seis estaciones tuvo cuatro estaciones que presentaba un muy buen ajuste con valores superiores a 0,90 y el resto con valores superiores a 0,70, En contraste, esta validación se asemeja a la comprobación en punto grilla que para nuestras nueve estaciones tenemos un desempeño variable, siendo las estaciones Campamento Sede, Puerto Pizarro y Tigre las que se acercan a valor de 0,90.
El Sesgo Estadístico Porcentual (PBIAS) tuvo un comportamiento de 3 estaciones sobreestimadas y 3 subestimadas con el mínimo valor de 5,236 y máximo de -27,36. En nuestro caso, en este indicador tenemos una un desempeño similar teniendo 4 estaciones subestimadas y 5 sobreestimadas con valores extremos de 5,4 y 49,5 en las estaciones Quebrada Angostura y Capitán Hoyle.
Mamani (2019) Obtuvo valores para Índice de Eficiencia de Nash- Sutcliffe (NSE) de bueno con valores superiores a 0,90 y una sola estación de 0,875. En nuestro estudio obtuvimos desempeños diversos con mayor ajuste de 0,88 y peor ajuste de -0,19 para las estaciones de Campamento Sede y Capitán Hoyle.
Zapana (2019) en su estudio intitulado Análisis de la sequía meteorológica de la Intercuenca Alto Apurímac y las cuencas Alto Urubamba y Yavero utilizando el método de L-Momentos, se realizó la validación del Producto PISCO para 24 estaciones en punto grillado, utilizo como métricas el Coeficiente de Determinación (R2) el cual, ha presentado resultado aceptable para su estudio, teniendo como valores extremos 0,543 y 0,946. Con respecto a nuestros resultados tienen mejor desempeño promedio sin considerar a las estaciones Condor Flores con 0,12 y Cabo Inga con 0,46.
Porto (2021) que realizó la validación de las estaciones de la Cuenca de la Vertiente del Lago Titicaca - Lado peruano, considero un total de 34 estaciones con el método del vector regional los agrupo en 5 regiones. Como indicadores utilizo el Coeficiente de Determinación (R2), obteniendo un desempeño muy bueno con valores superiores a 0,90, Para nuestra zona de estudio obtuvimos valores menores, pero con desempeño bueno.
Considerando el Índice de Eficiencia de Nash- Sutcliffe (NSE) se obtuvo un desempeño excelente con valores frecuentes de 0,90, En comparación con los resultados obtenidos en nuestro estudio que son de desempeño diverso teniendo un valore destacable de 0,89.
También considero el Sesgo Estadístico Porcentual (PBIAS) donde obtuvo un diverso desempeño en las estaciones, la región I y región V las estaciones en su totalidad es subestimada. Se encontró valores extremos de subestimación -27,094 y -0,399. Para el desempeño de sobreestimación se tiene 22,119 y 0,837. En contraste con los valores obtenidos, estos presentan mayor sobreestimación y subestimación con valores superiores.
También realizó estudio a nivel estacional obteniendo buenos resultado con algunas estaciones que tiene un desempeño regular, para los meses de marzo, abril y mayo se encontró valores de R2 = 0,392. En los meses de junio, julio y agosto un R2 = 0,293. En setiembre, octubre y noviembre un R2 = 0,356 y por último para diciembre, enero y febrero un R2 = 0,273. Esto en comparación para los valores obtenidos en nuestra zona de estudio en marzo, abril y junio (R2 = 0,05). Junio, julio y agosto (R2 = 0). En setiembre, octubre y noviembre (R2 = 0,01) y para diciembre, enero y febrero (R2 = 0,26). Lo cual desestima el uso de esas estaciones para dichos periodos.
También se realizó una validación según la región climática donde encontró un desempeño excelente en todas las métricas con una mínima subestimación, en nuestra área de estudio solo una región climática llega a tener un desempeño similar.
CONCLUSIONES
Se realizó el análisis exploratorio y consistencia de los datos de precipitación. A nivel mensual, en las 9 estaciones que se validó su ajuste en su correlación se obtuvo valores altos en un intervalo de 0,61 a 0,89 donde se puede optar por la utilización de los datos estimados, también se debe considerar el índice de Nash que se obtuvo en un intervalo de 0,60 a 0,88 y para la consideración de que valores de precipitación están sobreestimados o subestimados en un ajuste bueno se tiene un valor más o menos 15 hasta la convergencia a 0, De acuerdo a lo anterior se debe destacar a estaciones como Campamento Sede, Puerto Pizarro, Tigre, Rica Playa y Cedros.
Para la validación en arreglo estacional verano tiene la mayor correlación en las estaciones de acuerdo a las clasificaciones, teniendo para R= 4 estaciones con muy alta, para NSE = 3 estaciones con clasificación excelente y 4 estaciones con el mínimo PBIAS. La siguen las estaciones como otoño y primavera con similar ajuste en las estaciones meteorológicas. En cambio, invierno presenta el peor desempeño tanto en correlación, eficiencia y sesgo.
Se pudo identificar las regiones climáticas donde se ubicaban las estaciones. Se generó un promedio multianual por cada región.
La región climática D (i, p) A' presenta el mejor ajuste tanto en correlación con un valor de R=0,92. Un NSE= 0,91 y tiene una mínima subestimación en los datos con un -3,6 %. Por lo tanto, esa región climática evidencia un desempeño excelente con el Producto PISCO. No como las otras regiones climáticas como lo son C (i, p) A' con un desempeño bueno y D(i)A' con un desempeño mínimo.
En cuanto a elevación se obtuvo correlación media a mínima, teniendo un mayor ajuste a mayor elevación, lo cual se da en la Cuenca Alta, para el índice de eficiencia de Nash-Sutcliffe el desempeño es insatisfactorio, de igual manera para el Sesgo Estadístico Porcentual. En la Cuenca Baja se evidencia una sobreestimación y en Cuenca Medía se tiene un desempeño mínimo en general con una subestimación.
De acuerdo con los resultados obtenidos en la comprobación del Producto Pisco en nuestra área de estudio a nivel estación meteorológica, elevación, estacionalidad y región climática. para ello se puede considerar el uso de los datos de precipitación Para la realización de estudios hidrológicos como son análisis de oferta hídrica, fines de irrigación, balances hídricos, para zonas con escasa información o con información faltante, todo esto teniendo en cuenta su desempeño.