Aplicación de Algunos Métodos de Relleno A Series Anuales de Lluvia ...
Aplicación de Algunos Métodos de Relleno A Series Anuales de Lluvia ...
Aplicación de Algunos Métodos de Relleno A Series Anuales de Lluvia ...
, 7(1):1-20,2000
Instituto Meteorolgico Nacional,Ministerio del Ambiente y Energa, San Jos, Costa Rica
ABSTRACT
This paper shows the results obtained when different methods for filling gaps in annual precipitation series were
applied to records of meteorological stations, located in different regions of Costa Rica. Simple and multiple
regression were used. The ratio, the normal-ratio and the adjusted ratio methods were also utilized. The first, third
and fifth methods required just an additional station to fill the gaps, the others required two. The methods were
applied to 12 stations, from which 5 are located in the Caribbean Coast and 7 in the Pacifics. In total, the
information of 36 stations was used in this study. The Wilcoxon Sum Rank Test showed that all estimated series
have identical probability distributions to the observed series. The lower differences between observed and
estimated data are obtained when multiple regression is used. This is showed in the values obtained for MAE
(mean absolute error) and RMSE (root mean square error). Despite the former results, the maximum differences
between the observed and estimated values were very high in some cases (95%). A recent study, which uses an
autoregressive model of order 45 to estimate annual precipitation values shows that the highest differences with
this method are around 35%, which suggest the use of more complex methods if more accurate estimations are
required.
1. Introduccin
1
Direccin para correspondencia: Rosario. Alfaro, Instituto Meteorologico
Nacional, Apartado 7-3350-1000 San Jose, Costa Rica.
E-mail: [email protected]
42 ALFARO&PACHECO: Aplicacin de algunos mtodos de relleno se series anuales de lluvia de diferentes regiones de Costa Rica
Con el fin de comparar los resultados obtenidos, se Se determinaron las diferencias mximas
calcularon los errores MAE y RMSE, segn lo sugiere observadas entre los valores reales y los estimados.
Willmott (1982), ya que estos errores bsicamente miden el Adems, se compararon los resultados con aquellos
error promedio entre los valores estimados y los observados obtenidos en un estudio reciente (Leite and Peixoto. 1996)
mediante el clculo de las diferencias existentes entre dichos que utiliza un mtodo ms complejo para reproducir series
valores. Tambin se us el test de Wilcoxon (McClave and de datos de precipitacin anual.
Dietrich. 1991) para determinar si las poblaciones de las dos
muestras (la de valores observados y la de estimados) tenan
distribuciones de probabilidad idnticas.
CUADRO 2. Distancias y diferencias en altura existentes entre cada una de las estaciones analizadas y sus estaciones vecinas.
TOPICOS METEOROLOGICOS 43
Estacin en estudio D1 (km) D2 (km) Dif. alt. 1 (m) Dif. alt. 2 (m)
Puntarenas 8 7 137 447
Coto 47 5 3 -2 -1
CATIE 4 5 78 18
San Jos 13 22 -472 -10
Limn 33 52 35 10
Fabio 8 10 92 -380
Nicoya 22 63 -100 -105
Damas 10 4 -1 0
L. Vista 27 28 40 430
C.Quesada 8 24 -567 800
Pindeco 40 48 21 303
Upala 36 44 -10 5
Promedios 17.8 25.8 129.6 209.1
Rica, ya que existen regiones del pas que no cuentan con tom en cuenta el hecho de que las estaciones fueran
dos series de datos de estaciones vecinas cuyo perodo de cercanas a la estacin en estudio y se consideraron valores
registro sea igual al de la estacin en estudio. En las promedios de lluvia en diferentes perodos en vez de
regiones con una densidad alta de estaciones, el mtodo de precipitaciones normales.
regresin mltiple podra utilizarse en un futuro con ms Para obtener la precipitacin anual de una estacin
estaciones, con el fin de determinar si los valores estimados X (PX) para un ao determinado, utilizando los datos de dos
mejoran significativamente con respecto a aqullos estaciones A y B, se utiliz la siguiente ecuacin:
calculados con base en la informacin de dos estaciones
vecinas. PX = 1 ( NX PA + NX PB ) (3)
2 NA NB
2.2. Mtodo de la razn
donde NX es la lluvia promedio de la estacin X para el
El mtodo de la razn q (Barger, 1960, WMO 1966 mismo perodo que se obtiene la lluvia promedio de la
y WMO 1983) se basa en el hecho de que para pares de estacin A (NA) y B (NB). PA y PB son los valores
estaciones, la razn entre sus valores mensuales, anuales o correspondientes a PX (precipitacin anual) de las estaciones
medios, tiende a ser constante. Segn este mtodo, si se A y B.
tienen dos estaciones (A y B), se requiere determinar q de la
siguiente manera: 2.4. Mtodo de la razn ajustada
q = bi/ ai (1)
desde i = 1 hasta N Este mtodo es bsicamente el mtodo de la razn
donde N= nmero total de datos de la serie (seccin 2.2), nicamente que utiliza los promedios de lluvia
bi = dato i de la estacin B de aos lluviosos, secos y normales en vez del promedio
ai = dato i de la estacin A anual normal.
Seguidamente, si el valor bj falta en la estacin B, ste se Este mtodo, utilizado en forma ms meticulosa por
obtiene de la siguiente manera: Villalobos y Retana (1998) en otro estudio, para otros fines,
bj = q x aj (2) consiste en obtener una clasificacin de la precipitacin
Este mtodo es utilizado en la bibliografa anual de la estacin en estudio (B) y una estacin cercana
consultada para rellenar medias de diferentes perodos, pero (A), con base en el valor medio de sus series (Xp) y sus
se menciona que puede ser utilizado para valores mensuales desviaciones estndar (S), clasificando los totales anuales de
y anuales. En este estudio se utiliz particularmente para dichas estaciones segn se indica a continuacin:
obtener datos anuales, pero su aplicabilidad a datos Lluvioso si bi ( o ai ) > Xp +S,
mensuales ser investigada posteriormente. Normal si Xp-S bi ( o ai ) Xp+S
Seco si bi ( o ai ) < Xp-S
2.3. Mtodo de la razn-normal Una vez que se clasifican los totales anuales de
precipitacin de cada una de las estaciones, para calcular un
Este procedimiento estima cantidades de lluvia con dato faltante bi de la estacin en estudio, se utiliza la
base en tres estaciones cercanas y uniformemente espaciadas informacin obtenida anteriormente de la estacin cercana
con respecto a la estacin en estudio (Paulhus and Kkohler, para determinar a qu categora corresponde el ao que se
1952). Estos requerimientos son difciles de cumplir en desea estimar y de acuerdo a su clasificacin, se obtiene el
algunas regiones de nuestro pas, por lo que se utiliz el valor de q (q = bi/ ai) haciendo la sumatoria de los aos
mtodo con dos estaciones cercanas nicamente. correspondientes a esa categora. Posteriormente, si el total
Este mtodo se sugiere para cuando las diferencias anual faltante es bj , ste se calcula con base en la ecuacin
en las precipitaciones anuales normales de las estaciones (2), en donde aj representa el total anual correspondiente de
consideradas son mayores que un 10%. En el caso de las la estacin cercana.
estaciones del I. M. N, son pocas las estaciones que cuentan
con perodos de registro largos, por lo tanto, nicamente se 3. Datos utilizados
TOPICOS METEOROLOGICOS 45
4000
3000
2000
1972
1974
1976
1978
1980
1982
1984
1986
1988
1990
1992
1994
1996
Aos
Limn La Lola Fca el Carmen
Fig. 1. Series de datos anuales de la estacin Limn y estaciones cercanas (La Lola y Finca El Carmen).
.
3500
2500
1500
500
1963
1965
1967
1969
1971
1973
1975
1977
1979
1981
1983
Aos
Puntarenas Nagatac San M.Barranca
Fig. 2. Series de datos anuales de la estacin Puntarenas y estaciones cercanas (Nagatac y San Miguel de Barranca).
TOPICOS METEOROLOGICOS 47
3500
Lluvia anual en mm
2500
1500
500
1963
1965
1967
1969
1971
1973
1975
1977
1979
1981
1983
Aos
Puntarenas Q con Nag. Q con Barr.
Fig. 3. Datos reales de la estacin Puntarenas y valores estimados con el mtodo de la razn (Q) utilizando la informacin de las estaciones vecinas de Nagatac
(Nag) y San Miguel de Barranca (Barr).
2500
Lluvia anual en mm
1500
500
1963
1965
1967
1969
1971
1973
1975
1977
1979
1981
1983
Aos
Puntarenas Q con Nag. QA con Nag.
Fig. 4. Serie de datos observados de la estacin Puntarenas y valores estimados segn el mtodo de la razn (Q) y el de la razn ajustada (QA), usando la
informacin de la estacin vecina de Nagatac.
48 ALFARO&PACHECO: Aplicacin de algunos mtodos de relleno se series anuales de lluvia de diferentes regiones de Costa Rica
2500
Lluvia anual en mm
2000
1500
1000
500
1963
1965
1967
1969
1971
1973
1975
1977
1979
1981
1983
Aos
Puntarenas RG Nagatac RG Barranca
Fig. 5. Valores reales de la estacin de Puntarenas y datos estimados utilizando regresin simple (RG) con las estaciones de Nagatac y San Miguel de Barranca.
3000
Lluvia anual en mm
2500
2000
1500
1000
500
1963
1965
1967
1969
1971
1973
1975
1977
1979
1981
1983
Aos
Puntarenas RM QN
Fig. 6. Aplicacin del mtodo de regresin mltiple (RM) y de la razn normal (QN) para obtener datos estimados de lluvia anual para la estacin de
Puntarenas. Las estaciones cercanas utilizadas son Nagatac y San Miguel de Barranca.
TOPICOS METEOROLOGICOS 49
El estadstico TA se obtiene sacando el rango de las de aceptacin de TA para muestras con n=9 tiene como
observaciones de la muestra de valores observados y la lmite inferior y superior 63 y 108, respectivamente
muestra de valores estimados, como si las observaciones (McClave & Dietrich. 1991).
pertenecieran a la misma poblacin. El nmero de datos de Segn los resultados anteriores, en todos los casos,
la primera muestra se denota como n1 y los de la segunda se debe aceptar la hiptesis nula de que las poblaciones de
como n2 y corresponden a valores menores o iguales a 10. El las muestras tienen distribuciones de probabilidad idnticas.
test estadstico TA es la suma de rangos de la muestra ms Sin embargo, vale la pena destacar las diferencias mximas
pequea o en el caso en que n1=n2, se puede usar cualquier en % obtenidas entre los valores estimados y los observados
suma de rangos. segn cada mtodo, las cuales se muestran a continuacin en
Cuando el tamao de las dos muestras es igual a 10 los cuadro 7 y 8.
o ms grande, se usa el test estadstico Z, el cual est dado Obsrvese que de las 12 estaciones analizadas, 9 de
por la siguiente ecuacin: ellas tienen las diferencias mximas ms pequeas cuando se
utiliza el mtodo de la regresin mltiple para obtener los
TA n1(n1+ n2 + 1) valores estimados.
Z= ________2_____ (4) El mtodo de la razn ajustada muestra las
n1n2(n1+ n2+1) diferencias ms grandes en la mitad de las estaciones; sin
12 embargo, un detalle muy importante que vale la pena
destacar es el hecho de que la estacin de Limn, presenta la
La hiptesis nula utilizada en este test (Ho) mayor coincidencia de clasificacin de sus totales anuales de
establece que la poblacin de las 2 muestras tienen lluvia (84%) con respecto a los de la estacin vecina
distribuciones de probabilidad idnticas. La hiptesis utilizada en este mtodo (La Lola), la cual se encuentra a 33
alternativa (Ha) afirma que la distribucin de probabilidad kms de distancia y a 35 metros ms de altura. En este caso
para la muestra A est desplazada hacia la izquierda o la en especial, la diferencia mxima ms pequea se obtiene
derecha de la distribucin de probabilidad de la muestra B. utilizando este mtodo.
Para un nivel de significancia del 5%, la regin de
rechazo de Z corresponde a Z < -1.96 o Z > 1.96. La regin
CUADRO 5. Valores mximos y mnimos correspondientes a los errores MAE y RMSE obtenidos, segn cada mtodo de relleno.
CUADRO 6. Resultados (estadstico Z o TA) de la aplicacin de los mtodos de relleno para estimar totales anuales de precipitacin. Los valores de
50 ALFARO&PACHECO: Aplicacin de algunos mtodos de relleno se series anuales de lluvia de diferentes regiones de Costa Rica
CUADRO 7. Valores de las diferencias mximas en % obtenidas con cada uno de los mtodos, para cada una de las estaciones analizadas.
CUADRO 8. Diferencias mximas absolutas en % encontradas entre los mtodos de relleno utilizados.
valores reales de las series y los valores estimados con cada uno de los
TOPICOS METEOROLOGICOS 51