Tesis I Borrador
Tesis I Borrador
Tesis I Borrador
diferentes estudios como fuente de datos. Es por esto que se necesita analizar su
han convertido en una valiosa herramienta debido a que su registro es continuo a diferencia de
los registros de las estaciones meteorolgicas que en la toma suelen cometerse errores, y que se
Uno de los sensores ms utilizados es el TRMM debido a que presenta mayor data histrica de
la estimacin de lluvia (desde 1997) pero el 2014 sali el sensor GPM, como reemplazo al
TRMM, este presenta mejores caractersticas de resolucin temporal y espacial. Por lo que en
este trabajo se utiliz ambos sensores para hacer la comparacin con los registros de
este trabajo es evaluar las estimaciones de precipitacin por satlite utilizando la serie de datos
de las 5 estaciones activas de la provincia de Huamanga para el periodo 2014 al 2016. Los
satisfactorios al 100% para la escala mensual. Esta diferencia es atribuida por la topografa, el
clima complejo y la densidad de los pluvimetros en la zona que son utilizados por el sensor
GPM y TRMM para la correccin de las imgenes satelitales. Se obtuvo sobreestimaciones con
el sensor TRMM para las estaciones: Inia Cannan, La Quinua, San Pedro de Cachi y
El clima es por naturaleza variable, debido a que en el influyen distintos factores como
el relieve, latitud, corrientes marinas, etc. (Zuniga y Cresto, 2010; citado por Hurtado et
al., 2014) y uno de los elementos principales del sistema climtico es la precipitacin.
aplicaciones y usos tales como: estudios climticos, gestin del agua para la agricultura
y eventos extremos (Pea, 2016). Debido a que los pluvimetros ofrecen informacin
del punto donde son colocados, y por consecuencia, cuando las redes de pluvimetros
Durante los ltimos aos surgieron diversas tcnicas para estimar la precipitacin
pluvimetros (Joyce et al., 2004; Huffman et al., 2007 citado por Jos, 2009) que
precipitacin.
Es por esto que los productos de satlite se han convertido en una valiosa herramienta
metodologas (Gandia & Melia 1991). Uno de los productos recientes que fue diseado
para la percepcin remota de la precipitacin son las imgenes satelitales de la misin
GPM (Global Precipitation Measurement) (Andrade, 2016), que es una red internacional
Evaluar las estimaciones de precipitacin por satlite, en comparacin a los datos de las
Realizar las medidas de comparacin entre los datos de las estaciones meteorolgicas y las
meteorolgicas se va incrementando en diversas partes del mundo, por lo que; en el 2013 Lia
Ramos Fernndez realizo la tesis sobre el estudio de Utilidad de la precipitacin obtenida por
satlite en la modelacin hidrolgica aplicada a la cuenca del Rio Jucar concluye que los
tienden a ser ms precisos en un punto pero su significado espacial decae con la distancia, en
cambio las estimaciones por satlite ofrecen una mejor representacin espacial aunque sesgado.
mensual) en Regin Costa, Sierra y oriente del Ecuador teniendo como autor a Edisson Santiago
los valores del algoritmo TRMM 3B43 son ms acertados en la regiones Costa y Oriente que en
GPM (Global Precipitation Measurement) a escala sub-diaria para la provincia del Azuay del
autor Omar Fabricio Andrade Pillaga concluye en la demostracin de que las imgenes
satelitales tienen mejor comportamiento a medida que se aumenta la escala temporal, es decir
Teniendo estos estudios como base se piensa evaluar las estimaciones de Precipitacin en
nuestra regin.
2.2. FUNDAMENTO TERICO
2.2.1. PRECIPITACIN
2.2.1.1. Concepto
Villn (2002), menciona que la precipitacin, es toda forma de humedad que originndose en las
nubes, llega hasta la superficie del suelo, de acuerdo a esta definicin la precipitacin puede ser
primaria del agua de la superficie terrestre, y sus mediciones y anlisis, forman el punto
2.2.1.2. Clasificacin
en:
a. Precipitacin convectiva
Es causada por el ascenso de aire clido, que es ms liviano que el aire frio de los
b. Precipitacin orogrfica
Se producen cuando el vapor de agua que se forma sobre la superficie es empujada por
el viento hacia las montaas, aqu las nubes siguen por las laderas de las montaas, y
consiguiente precipitacin.
c. Precipitacin ciclnica
Se produce cuando hay un encuentro de dos masas de aire, con diferente temperatura y
humedad, las nubes ms calientes son violentamente impulsadas a las partes ms altas,
importantes.
2.2.1.3.1. Pluvimetros
Los instrumentos de medicin de precipitacin emplean una variedad de tecnologas
mantenimiento de rutina, pero son ms precisos que los pluvimetros basculantes. Los
El pluvimetro tradicional, almacena la lluvia diaria para luego ser medida, todo los
das a las 7:00 am, por un observador con la ayuda de una probeta calibrada.
Adems se tiene:
Pluvimetro electrnico: son usados en estaciones automticas para medir la lluvia.
Utiliza un sistema de cazoletas, que generan pulsos elctricos cada vez que se llenan y
lluvia durante el da. Con este instrumento es posible estimar la duracin de cada
evento.
precisas y en tiempo real para facilitar el ajuste a los sesgos de las estimaciones de
precipitacin de los radares y satlites, aportar datos de lluvia para los modelos
hidrolgicos y de crecidas repentinas y apoyar los pronsticos generales del estado del
pluvimetros que sean tan precisas, confiables y oportunos como sea posible.
como un rayo en una direccin predeterminada por una antena mvil. El ancho y la
irradiada, que viaja a la velocidad de la luz, se refleja parcialmente por las nubes y por
las partculas de precipitacin y regresa al radar, donde es recibida por la misma antena.
UCAR (2012) menciona que las redes de radares meteorolgicas tienen la funcin
cuadricula sobre una regin de inters. Los radares meteorolgicos son herramientas
Linsley (1977), sugiere que la informacin obtenida por satlites meteorolgicos pueden
ser utilizados para estimar cantidades de lluvia para periodos de un mes o ms. UCAR
productos que son subidos a enlaces desde el servicio meteorolgico que controla el
satlite.
conjunto numeroso de datos (serie de datos), que se necesita procesarlos para su anlisis y
a. Anlisis de consistencia
errores sistmicos, que se ponen en manifiesto como saltos o tendencias, afectando las
1. Anlisis grafico
Se basa en el criterio de que los valores acumulados de la precipitacin total anual de una
estacin, graficados con los valores acumulados de la estacin promedio, para un periodo
considerado, debe ser una lnea recta, de pendiente constante (sin quiebres), dependiendo de
3. Anlisis estadstico
Luego de observar la variacin del comportamiento de las series mensuales de precipitacin en
significativa y en base al anlisis de doble masa, se detectan los puntos de cambio, por lo que
iguales estadsticamente.
Si Tc >Tt entonces X1 X2 : el salto es significativo y las medias muestrales son
la desviacin estndar.
Si Sc >St entonces S1 S2: el salto es significativo, existe diferencia estadstica y es
necesario su correccin.
4. Correccin de la informacin dudosa.
En los casos en que los parmetros media y desviacin estndar resulten estadsticamente
periodo dudoso, mediante una ecuacin que permite mantener los parmetros del periodo ms
largo a otra con registro incompleto y corto. Las tcnicas que se utilizan son las de
2.2.2. TELEDETECCIN
2.2.2.1. Concepto
Cualquier procedimiento o tcnica de adquisicin de informacin sin tener contacto directo con
ella. Adems es entendida como una tcnica que tiene el objetivo de captura, tratamiento y
Visin panormica
En la tabla 2.1. se describen las principales bandas espectrales descritas por los sensores
Azul (0.45 0.5 um) Penetracin del agua, uso de la tierra, caractersticas de la
vegetacin, sedientos.
Infrarrojo reflectivo (0.75 0.9) Biomasa, identificacin de cultivos, lmites de tierra y agua.
Infrarrojo medio (1.5 1.75) Cantidad vegetacin, zonas ridas, nubosidad, hielo, nieve.
Microonda de onda larga (5 Medidas del suelo, lmites entre tierra y agua, penetracin de la
24 cm) vegetacin.
Fuente: http://www.javeriana.edu.co/biblos/tesis/ingenieria/tesis123.pdf
2.2.2.2. Aplicaciones
medidas de aguas superficiales y humedales para evaluar la situacin del habitad para las
aves acuticas.
Hidrologa: balance del agua, agro hidrolgico, inventarios del agua superficial,
acumulacin nival.
Agricultura: cartografa de cultivos y usos del suelo, control de erosin, manejo del agua,
Entre otros.
Sensores pasivos:
Estos miden energa radiante reflejada o emitida por el sistema tierra atmosfera. La energa
Sensores activos
Estos sensores lanzan rayos de radiacin sobre el sistema tierra-atmosfera miden la radiacin
1. Fotografas areas
corresponder con la visin ocular normal. La fotografa, adems del espectro visible, puede
las ondas electromagnticas que se reflejan o que emiten los objetos, captando radiaciones
georeferenciados.
1. Resolucin espacial
Capacidad del sistema para distinguir objetos en funcin de su tamao, interpretada de otra
forma es el tamao del objeto ms pequeo que puede ser distinguido en una imagen: tamao
2. Resolucin espectral
Indica el nmero y anchura de las bandas espectrales que puede discriminar el sensor. Entre
mayor sea la resolucin espectral tendr informacin del comportamiento de una misma
3. Resolucin radiomtrica
formato binario (nmero de bits) para codificarse. A mayor numero de bits reservados para
4. Resolucin temporal
Indica el tiempo de paso del sateltes sobre la vertical de un punto. Con ello se determina la
5. Resolucin angular
Capacidad que tiene un sensor para tomar imgenes oblicuas. Esta propiedad permite no
solamente acortar los ciclos temporales de recubrimiento, sino que posibilita la generacin de
Cobertura global, pero solo una o dos o menos mediciones al da por sensor. Existen
temporal.
Menos de una imagen al da. Cobertura no global y existen lagunas orbitales. Mientras
Ramos (2013) menciona que los sensores acoplados en satlites meteorolgicos realizan
lecturas en 5 bandas del espectro electromagntico para estimar la lluvia con diversas
tcnicas que estn en constante avance hacia tcnicas ms directas basadas fsicamente,
que han ido evolucionando desde medidas de la radiancia en el espectro visible (Vis) e
infrarrojo (IR), tcnicas basadas en microondas activas y pasivas (MW) hasta tcnicas
tabla 2.2.
emitida por las nubes, es decir es estimada a partir de la radiacin microonda emitida o
precipitacin principalmente a travs de la energa emitida por las gotas de lluvia que es
Estos satlites son una misin conjunta entre la NASA y la Agencia Japonesa de
El satlite TRMM Tropical Rainfall Measuring Mission entra en rbita en 1997 para
monitorear y estudiar lluvias tropicales, lleva un radar que transmite en una longitud de
Radar (PR), TRMM Microware Imager (TMI) (Captador de imgenes microondas TRMM),
Visible and Infrared Scanner (VIRS) (Escaner visible e infrarrojo) y un sensor de imgenes de
relampagos (LIS)
Figura 2.13. Instrumentos a bordo del satlite TRMM: PR, VIRS, TMI y LIS y un sistema de energa radiante de ls superficie
terrestre que emplea una base de datos de imgenes de nubes de altsima resolucin.
Fuente: https://arset.gsfc.nasa.gov/sites/default/files/users/Spanish_GPM_Week2_March24.pdf
precipitacin)
Tambin conocida como TRMM 3B43, que es un algoritmo del TRMM, version 7 que combina
los resultados de varios sensores para obtener mejor cobertura espacio temporal:
La resolucin temporal del TRMM 3B43 es mensual y la resolucin de 0.25 x 0.25 (~27x27
km) cubre un cinturn global desde los 50latitud Sur a los 50 Latitud Norte.
2.2.3.4.2. GPM
extremos que causan riesgos naturales y desastres y extender las capacidades actuales de
usar informacin precisa y oportuna de precipitacin para beneficiar directamente
sociedad.
El satlite GPM Core tiene un sensor de lluvia activo: Dual-frequency Precipitation Radar
de la intensidad de la lluvia, y uno pasivo: GPM Microware Imager (GMI), capaz de generar
imgenes con intensidades de precipitacin; estos dos sensores mejoran sobre el producto de
TRMM PR y TMI.
Existen diferentes algoritmos que usa el observatorio Core de la misin GPM para el
Combina datos del GPM GMI/DPR con los satlites de la constelacin del GPM para rendir
de 0.1 x 0.1 (~10km x 10km) y una cobertura espacial global de 60S a 60N.
Coordenadas utm
Cuadricula: 27
Figura 3.2. Mapa de ubicacin satelital de las 5 estaciones meteorolgicas de la provincia de Huamanga
los registros de estimacin de precipitacin por imagen satelital, los registros obtenidos
Fuente: propia
confiabilidad de la evaluacin del producto satelital. Por ello es necesario contar con
una base de datos veraz para comprobar las estimaciones del satlite. Al registrar los
datos de lluvia los pluvimetros pueden sufrir errores, por problemas con el equipo, mal
manejo del instrumento por parte del operador, cambios y ubicacin del instrumento,
2014), el porcentaje de vacos de las estaciones se debe tener un valor mximo de 20%
Para la descarga de las imgenes satelitales de estimacin de precipitacin de los dos algoritmos
o KMZ y Geotiff, para este trabajo se descarg en NetCDF (nc), con esta plataforma uno puede
lluvia estimada por satlite frente a los registros de precipitacin por estacin meteorolgica.
Las medidas de sesgo, raz del error medio cuadrado, coeficiente de correlacin de
a. Porcentaje de sesgo
satlite.
( Ppxi )
i=1
Donde:
N: nmero de estaciones
Evala el error de las estimaciones por satlite en comparacin con las estaciones
n
1
RMSE= ( Ps xiPp xi)2
n i=1
Es una medida de concordancia lineal entre las estimaciones por satlite y los datos de
mejor concordancia.
d. ndice de concordancia
Es una medida estandarizada del grado de error de prediccin. Este ndice tiene un valor
Donde:
meteorolgicas utilizadas
Fuente: propia
Las 5 estaciones cumplen con el valor admitido de 20%. Por lo que se trabaj con las 5
estaciones.
Anlisis de consistencia
Object 13
Figura 4.1. Histograma del registro de lluvia de la estacin Inia Cannan del periodo 2011-2016
Object 15
Figura 4.2. Histograma del registro de lluvia de la estacin La Quinua del periodo 2010-2015
Object 17
Figura 4.3. Histograma del registro de lluvia de la estacin San Pedro de Cachi del periodo 2010-2015
Object 20
Figura 4.4. Histograma del registro de lluvia de la estacin Cuchoquesera del periodo 1991-2016
Object 23
Figura 4.5. Histograma del registro de lluvia de la estacin Wayllapampa del periodo 2010 - 2015
En todas las estaciones el Tc fue menor que el Tt, para la media, por lo que no se
requiri la correccin. Tambin el Fc fue menor que el Ft para cada estacin (saltos),
Las tendencias de las series de datos a lo largo de los aos de registro, se observan en la
curva de doble ms relacionando a las 5 estaciones con una estacin promedio, que nos
Tabla 4.1. Registro de PP anual PP cumulada para las 5 estaciones del periodo 2010
2016
Object 25
Tabla 4.2. Registro de PP mensual promedio para las 5 estaciones del periodo 2010
2016
Object 27
Al tener las imgenes descargadas en formato nc de la plataforma Giovanni, para cada uno de
los algoritmos, se extrajo el registro de precipitacin intersecando con cada pixel que coincida
complemento de QGis: Point Sampling Tool, extrayndolo en una nueva tabla de atributos.
Para poder extraer los registros se tiene que tener en cuenta que el raster (la imagen satelital)
este en las mismas coordenadas que los Shape (de cada estacin), por lo general los raster
Adems se tiene que tener en cuenta la resolucin espacial de la imgenes debido a que en la
imagen del algoritmo TRMM 3B43 que es de 0.25 x 0.25, al obtener el registro de las
estimaciones de las estaciones Wayllapampa, Inia Cannan y La Quinua presentan los mismos
valores en la serie de precipitacin, ya que por estar cercanas, coinciden en la misma grilla del
Una vez extrados los registros de lluvia de los dos algoritmos, se procedi a analizar la
estacionalidad con respecto a cada estacin pluviomtrica, para ello se realiz el histograma del
En el caso del algoritmo TRMM 3B43 se utiliz el periodo 2010-2016 debido a que este
algoritmo presenta registros desde 1998. Pero para el anlisis de la estacionalidad se utiliz del
2014 2016.
Object 30
Object 32
Object 34
Object 36
Object 38
Figura 4.11. Histogramas de los registros de precipitacin de las estaciones meteorolgicas frente a los registros obtenidos de las
imgenes satelitales del algoritmo TRMM 3B43,
En el caso del algoritmo GPM 3IMERG en el periodo 2014-2016, ya que este satlite fue
lanzado recientemente como un reemplazo al TRMM y tiene registros desde abril del 2014, este
Object 42
Object 44
Object 46
Object 48
Figura 4.12. Histogramas de los registros de precipitacin de las estaciones meteorolgicas frente a los registros obtenidos de las
imgenes satelitales del algoritmo GPM 3IMERG
Object 50
Object 52
Object 54
Object 56
Object 58
Figura 4.13. Precipitacin acumulada de las estaciones meteorolgicas frente a la precipitacin acumulada de las imgenes
satelitales del algoritmo GPM 3IMERG
Object 60
Object 62
Object 64
Object 66
Object 68
Figura 4.14. Precipitacin acumulada de las estaciones meteorolgicas frente a la precipitacin acumulada de las imgenes
satelitales del algoritmo GPM 3IMERG
precipitacin en forma cuantitativa se realiz los errores estadsticos: Porcentaje de sesgo, raz
del error medio al cuadrado (RMSE), coeficiente de correlacin de Pearson (y) y concordancia
de acuerdo (d). En las tablas 4.3. y 4.4. muestran los valores promedios de cada error estadstico
a. Porcentaje de sesgo
TRMM
El sesgo relativo muestra valores de sobreestimacin en las estaciones de: Inia Cannan
(20.426). En el caso de la estacin Inia Cannan nos muestra un porcentaje elevado del
sesgo a diferencia de las dems esto nos indica que se estara considerando como
subestimacin.
GPM
(11.842 %). Para el caso de las estaciones: La Quinua (-8.358), San Pedro de Cachi
Figura 4.15. Porcentaje de Sesgo calculado para cada estacin meteorolgica con respecto a los registros de estimacin de
precipitacin de los dos satlites TRMM y GPM
La raz del error medio cuadrado, presenta valores moderados debido a la cantidad de
datos con los que se cont para la evaluacin, estas se muestran en la figura 4.16., si se
contase con ms registro o a nivel diario se podra observar el RMSE menores alas
obtenidas.
Object 73
Figura 4.16. RMSE calculado para cada estacin meteorolgica con respecto a los registros de estimacin de precipitacin de los
dos satlites TRMM y GPM
Al observar la figura 4.16. se muestra que el error que presenta el algoritmo de GPM
Los clculos realizados para el coeficiente de Pearson nos dieron valores cercanos a la
Object 75
Figura 4.17. Coeficiente de Pearson calculado para cada estacin meteorolgica con respecto a los registros de estimacin de
precipitacin de los dos satlites TRMM y GPM
Vemos la diferencia entre los dos algoritmos de los satlites TRMM y GPM, el
algoritmo TMPA segn muestra la figura 4.17, excepto para la estacin Cuchoquesera,
d. ndice de concordancia
Figura 4.18. ndice de concordancia calculado para cada estacin meteorolgica con respecto a los registros de estimacin de
precipitacin de los dos satlites TRMM y GPM
con el algoritmo del satlite GPM a diferencia del satlite TRMM y en el resto de las
estaciones se da el caso contrario ya que el algoritmo del satlite TRMM tiene mayor
concordancia frente al GPM, pero la diferencia nos es demasiado por lo que ambas
meteorolgica.
Tabla 4.3. Resumen de los errores estadsticos realizados entre las estacione
Tabla 4.4. Resumen de los errores estadsticos realizados entre las estacione
meteorolgicas frente a los registros del satlite GPM
V. CONCLUSIONES Y RECOMENDACIONES
5.1. Conclusiones
registro utilizado fue del periodo 2010 al 2016, pero para la comparacin con las estimaciones
de precipitacin estimada por los dos algoritmos TRMM 3B43 y GPM 3IMERG solo se utiliz
el registro de 2014 al 2016, en el caso del algoritmo 3IMERG se utiliz a partir del mes de Abril
debido a que recin la data estuvo disponible en este mes pero fue lanzado 27 de febrero.
Para tener datos de lluvia de calidad de las 5 estaciones se analiz la homogeneidad mediante el
anlisis grfico, anlisis de los histogramas y curva de doble masa, y el anlisis estadstico,
pruebas estadsticas paramtricas de T de Student y F de Fisher, que dieron como reultado que
Se obtuvo los registros de precipitacin de las imgenes satelitales mediante la plataforma Nasa
cualidad, con un error menor con el sensor GPM a diferencia del sensor GPM.
Para la comparacin de los registro de lluvia estimado por satlite frente a los registro de
estacin meteorolgica se emplearon los errores estadsticos lo cual nos indica que en forma
una correlacin cercana a 1 y tienen una comportamiento similar y son adecuados, por lo que es
escasa, y es una alternativa favorable para la utilizacin en zonas donde no contamos con
estaciones para el registro de lluvia, para los diferentes estudios donde se requiere la obtencin
de registros de lluvia.
5.2. Recomendaciones
poder utilizarlo recurrimos a diversos mtodos para poder obtener para una zona (cuenca) donde
plataformas como: Mirador y Giovanni, que no es necesario tener una cuenta, para poder
descargar imgenes de estimacin de precipitacin por satlite donde que nos permite elegir la
resolucin temporal y espacial de la zona de estudio, por lo que es recomendable para su uso.
Al momento de elegir el satlite se tiene que tener en cuenta las caractersticas del satlite, es
decir la resolucin temporal, como es el caso de los algoritmos empleados en este trabajo que
que presentan como TMPA tiene registros desde noviembre de 1997 e IMERG desde abril de
2014.
requerimos licencias, pero ltimamente se fue viendo el ingreso de software libres como es el
caso del QGis y otros, la cual es recomendable su uso debido a que uno no requiere obtener la
licencia para su uso y adems podemos encontrar diversos tutoriales para su aprendizaje en el
YouTube.
Para poder utilizar un registro de datos de lluvia siempre es necesario hacer un anlisis de
bondad para que nuestros datos sean confiables y hacer el respectivo proceso para su utilizacin,
VI. BIBLIOGRAFA
Hurtado, P. B.; Sanches, I.; Velasquez,A. M.; Esquivel, G.; Palomo, M. (2014).
Propuestas metodolgicas para el anlisis de calidad de datos de precipitacin.
Recuperado de: http://www.agrofaz.mx/wp-
content/uploads/articulos/2014143IV_1.pdf
VII. ANEXO
7.1. Plano de ubicacin
7.2. Anlisis estadstico de los registros de precipitacin de las estaciones meteorolgicas