13.

MAXIN 14/6/07 14:51 Pgina 137

MAXIN: Aplicacin SIG para la

estimacin de valores
Intensidad-Duracin-Frecuencia de
precipitaciones en la Espaa peninsular

L. DE SALAS REGALADO (*), L. CARRERO DEZ (*) y J. A. FERNNDEZ YUSTE (*)

RESUMEN MAXIM es una aplicacin informtica desarrollada a partir de un sistema de informacin geogrfica. MAXIN se apo-
ya en los recientes estudios desarrollados por de Salas (2005) para la actualizacin y adecuacin de las leyes Intensidad-Dura-
cin-Frecuencia (IDF) en la Espaa peninsular, y permite, mediante un mapa interactivo, calcular la intensidad mxima para
cada duracin y perodo de retorno en cualquier punto de la Espaa peninsular. Esta aplicacin, de libre difusin, est disponible
en: http://www.forestales.upm.es (abriendo la seccin, organizacin docente, unidades docentes, hidrulica e hidrologa, enlaces).

MAXIN: A GIS APPLICATION TO ESTIMATE INTENSITY-DURATION-FRECUENCY VALUES OF

PRECIPITATION IN THE SPANISH PENINSULAR TERRITORY
ABSTRACT MAXIN is a Geographic Information System application. It was developed as a result of de Salas investigations
to up-date and obtain more adequate Intensity-Duration-Frequency (IDF) relationships in the Spanish peninsular territory
(2005). MAXIN estimates maximum precipitation intensity for d-duration and T- return period everywhere in the Spanish
peninsular area when pointing at the place of study in an interactive map. Downloading is free and it is available in:
http://www.forestales.upm.es.

Palabras clave: SIG; Leyes IDF, Anlisis regional.

1. ANTECEDENTES entonces no estaban disponibles y hoy ofrecen series

continuas y amplias.
La ley IDF actualmente utilizada en Espaa, y recogida en
la Instruccin 5.2-IC (1990), se desarroll a finales de los b) Adems, la funcin patrn elegida para caracterizar el
aos setenta (Ministerio de Obras Pblicas, Transporte y comportamiento de la variable aleatoria, Gumbel, con
Medio Ambiente, 1978), modificndose ligeramente el mapa un coeficiente de asimetra fijo 1,14, subestima siste-
de isolneas para la generalizacin a todo el territorio en mticamente los cuantiles extremos al menos en el
1987 (Ministerio de Obras Pblicas, Transporte y Medio rea de influencia mediterrnea (Bachro and Chao-
Ambiente, 1987). La metodologa entonces propuesta, su- che, 2006).
pone el ajuste de cada una de las series de intensidades m- Con estos antecedentes, de Salas (2005) plante la nece-
ximas anuales, Id, (series correspondientes a duraciones de sidad de actualizar estas leyes, incorporando:
5, 10, 15, 20, 30, 60, 120, 180, 360, 720 minutos y 24 horas) a) Los nuevos registros recogidos en las estaciones ya
a la funcin Gumbel, en cada una de las 21 estaciones con utilizadas para desarrollar las leyes IDF en los aos
pluvigrafo entonces disponibles en Espaa con un nmero setenta, y las nuevas estaciones que hoy ya dispo-
significativo de registros. nen de un nmero adecuado de datos. Todas las es-
Del prrafo anterior se deduce que esa ley, que es la que taciones pertenecen a la red pluviogrfica del Insti-
actualmente se utiliza, presenta serias limitaciones: tuto Nacional de Meteorologa (INM) quien facilit
a) En primer lugar por el nmero de datos y estaciones los datos.
manejados, ya que en los aos transcurridos, las esta- b) Una funcin de distribucin, SQRT-ETmax, ms ade-
ciones utilizadas han incrementado sustancialmente cuada para la variable aleatoria de valores extremos
sus registros y, adems, pueden utilizarse otras que estudiada.
c) La aplicacin de tcnicas de regionalizacin intra-
estacin, que permiten mejorar la robustez de las
(*)E.U.I.T. Forestal. Universidad Politcnica de Madrid. Avda Ramiro de estimas correspondientes a perodos de retorno
Maeztu s.n. 28040 Madrid. leticia.salas.regalado@upm.es altos.

Ingeniera Civil 146/2007 137

13.MAXIN 14/6/07 15:45 Pgina 138

MAXIN: APLICACIN SIG PARA LA ESTIMACIN DE VALORES INTENSIDAD-DURACIN-FRECUENCIA DE PRECIPITACIONES EN LA ESPAA PENINSULAR

2. MATERIAL Y MTODOS Una vez establecidas dos regiones en cada estacin se

procedi conforme a la siguiente secuencia:
2.1. DATOS
a) Eleccin de una funcin de distribucin. La muestra
Se han utilizado las series de intensidades mximas anuales nica y adimensional formada por el conjunto de las
para duraciones de 5, 10, 15, 20, 30, 60, 120, 180, 360, 720 series correspondientes a distintas duraciones de la
minutos y 24 horas, facilitados por el INM en 63 estaciones regin homognea, se ajust a la funcin de distribu-
del territorio peninsular espaol. cin SQRT-ET mx porque es ms conservadora que
la tradicional ley de Gumbel, por tener buena capaci-
2.2. METODOLOGA dad descriptiva de la rama derecha de la distribucin,
Estimacin IDF en las estaciones y porque ha sido utilizada satisfactoriamente en Es-
paa para la estimacin de cuantiles de precipitacio-
Con el fin de mejorar la robustez de las estimas para pero- nes mximas diarias anuales (Ferrer y Ardiles, 1994).
dos de retorno altos, se opt por aplicar la tcnica de regio-
nalizacin. Bsicamente esta tcnica consiste en establecer b) Estimacin de los parmetros de la funcin en cada
regiones homogneas, entendiendo como tales aquellas regin. Se utilizaron los mtodos de momentos y de
en la que las estaciones pueden diferir una de otra en me- mxima verosimilitud, y se adopt en cada caso aqul
dia, pero no difieren significativamente ni en coeficiente de que proporcionaba unas estimas lo ms prximas a los
variacin ni en coeficiente de asimetra. As, si en las se- valores observados, pero siempre del lado de la seguri-
ries de una regin homognea se elimina el factor local di- dad para los valores usuales de perodo de retorno que
vidiendo por su media se obtienen muestras que se puede se manejan. Los criterios de decisin se basaron en el
aceptar corresponden a una nica poblacin subyacente de Test de Kolmogorov-Smirnov y el anlisis grfico.
media uno, pudindose entonces agruparse y conformar c) Validacin. Se generaron 100 muestras sintticas por
una muestra nica representativa de la regin. Una vez el mtodo de simulacin de Montecarlo, y se compara-
ajustada una funcin de distribucin con los registros de la ron los coeficientes de variacin y asimetra con los co-
serie nica de la regin y obtenidos los cuantiles adimen- rrespondientes a los valores observados.
sionales para los perodos de retorno habituales, la intensi- d) Finalmente se comprob que este anlisis regional
dad mxima anual para una duracin, d, y un perodo de proporcionaba estimas ms robustas que el anlisis
retorno, T, -I(d, T)- se puede obtener a partir de la si- local.
guiente expresin (1):
Los detalles tcnicos de la secuencia descrita pueden con-
_ sultarse en de Salas y Fernndez (en prensa).
I(d,T)estacin = CA(T)regin* I d estacin (1)

Donde: IDF generalizada

CA(T)regin, es el cuantil adimensional, para un perodo de El siguiente paso consisti en extender los resultados obte-
retorno,
_ T, obtenido en la regin. nidos en las estaciones que disponen de pluvigrafo al resto
I d estacin, es la media de las intensidades mximas anuales del territorio peninsular.
para una duracin, d, obtenida en la estacin. Para ello se estableci una relacin entre I(d,T) obtenida
El enfoque regional pretende resolver el problema de ca- tras el proceso de regionalizacin y la I(24,T) en la misma
rencia de datos, paliando la carencia en el tiempo con la estacin, dato disponible en todo el territorio peninsular a
abundancia en el espacio (Hosking y Wallis, 1997). partir de la publicacin Mximas lluvias diarias en la Es-
En este caso, aunque el nmero de estaciones y regis- paa peninsular del Ministerio de Fomento (1999). Este co-
tros es muy superior al utilizado para la estimacin de las ciente se puede expresar como la ecuacin (2):
leyes IDF actualmente en uso, siguen siendo claramente I (I(d; )
d; TT) CA PG (T )cortas/largas I PG d CA PG (T )cotas/largas I PG d
insuficientes para caracterizar el territorio, mxime si, = = (2)
como ocurre en Espaa, la variabilidad climtica es tan I(I(24; )PG
24; TT) PG CA PG (T )largas I PG 24 CA PG (T )Ferrer I PG 24
importante, desde un clima semirido en el sureste espa-
ol, hasta atlntico en el norte. Y esa insuficiencia se ma- Donde:
nifiesta tambin en la imposibilidad de establecer regio- PG hace referencia a una estacin con pluvigrafo.
nes homogneas con un nmero significativo de CAPG(T)cortas/largas es el cuantil adimensional deducido en la
estaciones. Por esa razn, se propuso una regionalizacin regin correspondiente (cortas o largas duraciones).
intraestacin, conformado en cada estacin como re- CAPG(T)Ferrer, es el cuantil adimensional obtenido tras el an-
gin por un lado las series correspondientes a duraciones lisis regional de las precipitaciones mximas diarias, dispo-
menores o iguales a una hora, que se denomin regin nible en la citada publicacin. El subndice Ferrer hace refe-
cortas duraciones, por otro lado las mayores de una rencia
_ al apellido de uno de los autores de dicho trabajo.
hora, largas duraciones. Diversos estudios ponen de ma- I PGd, es la intensidad mxima media para la duracin d
nifiesto diferencias significativas en cuanto al comporta- obtenida en una estacin con pluvigrafo.
miento de la intensidad mxima anual al considerar ven- En las estaciones sin pluvigrafo, los dos trminos del
tanas temporales menores o mayores de una hora numerador no son conocidos, por lo que, en un intento de ge-
(Hershfield, 1961; Bell, 1969; Ferreri y Ferro, 1990; Porras neralizar los resultados obtenidos en las estaciones con plu-
y Porras, 2001; Llasat, 1998), diferencias que se pudieron vigrafo, se analizaron dos funciones, una que relacionara
acreditar siguiendo el protocolo propuesto por Hosking y el cociente de cuantiles con el perodo de retorno, h(T), y otra
Wallis (1997). Esta metodologa est basada en el anlisis que relacionara el cociente de intensidades medias con la
de los L-Momentos y establece tres filtros cada vez ms duracin, g(d), ecuacin (3).
exigentes en cuanto a control de la homogeneidad regio-
nal: 1 Grfico; 2 Test Medida de la Discordancia; 3 Test I (d; T ) CA PG (T ) cortas / l arg as I PG d
= PG h(T ) g(d) (3)
de Heterogeneidad. I (24; T )PG CA PG (T ) Ferrer I 24

138 Ingeniera Civil 146/2007

13.MAXIN 14/6/07 14:51 Pgina 139

MAXIN: APLICACIN SIG PARA LA ESTIMACIN DE VALORES INTENSIDAD-DURACIN-FRECUENCIA DE PRECIPITACIONES EN LA ESPAA PENINSULAR

Zona 1

Zona 2 Zona 1

Zona 1

FIGURA 1. Mapa de zonas para la

funcin h(T) con cortas duraciones.

Para la funcin h(T) se hizo un anlisis previo del cociente cuando se consideraban las regiones de cortas duraciones,
de cuantiles y se comprob que su comportamiento era dife- y otras dos cuando se contemplaban las de largas duracio-
rente creciente o decreciente con el perodo de retorno- se- nes (figuras 1 y 2), y se ajust una funcin de tipo logart-
gn la localizacin de la estacin considerada y variando mico en cada zona: h(T) = a(LnT)2+b(Ln T)+c, donde Ln, es
tambin segn se tratase de la regin cortas o de la de el logaritmo neperiano, T es el perodo de retorno en aos, a,
largas. Por lo tanto, se definieron dos zonas geogrficas b y c son tres parmetros que adoptan los valores de la tabla
para caracterizar el comportamiento de la funcin h(T) 1 en cada zona, y R2 es el coeficiente de correlacin.

Zona 1 Zona 1

Zona 2 Zona 1

Zona 1

FIGURA 2. Mapa de zonas para la

funcin h(T) con largas duraciones.

Ingeniera Civil 146/2007 139

13.MAXIN 14/6/07 14:51 Pgina 140

MAXIN: APLICACIN SIG PARA LA ESTIMACIN DE VALORES INTENSIDAD-DURACIN-FRECUENCIA DE PRECIPITACIONES EN LA ESPAA PENINSULAR

Cortas duraciones Largas duraciones

Zona 1C Zona 2C Zona 1L Zona 2L

a 4*10 -4
7*10 -3
1,2*10 -3
3,7*10-3

b 9,2*10-3 106,6*10-3 13,6*10-3 55*10-3

TABLA 1. Valores de los c 1,004 0,909 1,022 0,954

parmetros a, b, y c de la
funcin h(T) en cada zona. R 2
0,991 0,998 0,991 0,959

Para la funcin g(d) se propuso una funcin biparam- mxima diaria anual para el punto considerado, que
trica, ecuacin (4), siendo K un parmetro con significacin puede obtenerse con MAXPLU, aplicacin informtica dis-
fsica e igual al cociente de las intensidades mximas me- ponible en la citada publicacin del Ministerio de Fomento
dias en 1 hora y la correspondiente a 24 horas. (1999) y F el factor de mayoracin que permite pasar de
precipitacin diaria a 24 horas (suelen utilizarse valores
24a d a
Id que se mueven entorno a 1,13).
g(d) = =K 24a 1 (4)
I 24
3. RESULTADOS
Los valores de los parmetros K y a se obtuvieron en cada
una de las 63 estaciones y su generalizacin a todo el territo- La aplicacin MAXIN permite estimar I(d,T) en cada punto
rio peninsular se hizo trazando sendos mapas de isolneas, fi- de la Espaa peninsular a partir de la expresin (5). Para
guras 3 y 4. En el caso de K se utiliz el mtodo de interpola- desarrollar la aplicacin MAXIN fue necesario georreferen-
cin SPLINE, con 6 estaciones de referencia y una resolucin ciar toda la informacin espacial, utilizando para ello las co-
de 1.000m, mientras que para a se realiz mediante el m- ordenadas UTM referidas al huso 30, y la resolucin espa-
todo de interpolacin del inverso del cuadrado de la distancia cial adoptada fue 1.000m x 1.000m. Las capas de
con 12 estaciones de referencia y resolucin 1.000 m. informacin utilizadas fueron:
Finalmente el valor de I(d,T) en cualquier punto puede
estimarse a partir de la ecuacin (5): PMDA Precipitacin mxima diaria anual. Mapa de
isolneas obtenido mediante mtodo Kriging,
P(24, T ) punto resolucin 1000 m, y malla de 2.500m.
I (d, T ) punto = I (24, T ) punto h(T ) g(d) = h(T ) g(d) (5) Cv: Coeficiente de variacin. El mapa de isolneas de este
24
parmetro fue obtenido mediante interpolacin espacial por
donde P(24,T) punto puede obtenerse como el mtodo del inverso del cuadrado de la distancia, con las
PMDA(T) punto*F, siendo PMDA(T) punto la precipitacin 12 estaciones ms prximas, y una resolucin de 1.000m.

FIGURA 3. Mapa de isolneas del

parmetro K.

140 Ingeniera Civil 146/2007

13.MAXIN 14/6/07 14:51 Pgina 141

MAXIN: APLICACIN SIG PARA LA ESTIMACIN DE VALORES INTENSIDAD-DURACIN-FRECUENCIA DE PRECIPITACIONES EN LA ESPAA PENINSULAR

FIGURA 4. Mapa de isolneas del

parmetro a.

Tanto esta capa como la anterior se obtuvieron a partir de Zonas 1L, 2L, 1C y 2C, zonas geogrficas en cortas y lar-
los datos ofrecidos por el software incluido en Mximas llu- gas duraciones, para la eleccin de la ecuacin h(T) ade-
vias diarias en la Espaa Peninsular (1999). cuada.
I1 MAXIN presenta al usuario un mapa interactivo (figura
K= y a, parmetros de la funcin g(d). 5) en el que se debe sealar con el cursor el punto en estu-
I24
dio. Se muestra el polgono donde se encuentra el punto, po-

FIGURA 5. Mapa interactivo de la pantalla principal de MAXIN.

Ingeniera Civil 146/2007 141

13.MAXIN 14/6/07 14:51 Pgina 142

MAXIN: APLICACIN SIG PARA LA ESTIMACIN DE VALORES INTENSIDAD-DURACIN-FRECUENCIA DE PRECIPITACIONES EN LA ESPAA PENINSULAR

FIGURA 6. Tabla de estimas IDF ofrecidas por MAXIN.

lgono en el que los valores de los distintos parmetros son Actualizacin de la informacin pluviogrfica, con los
constantes. Esos valores aparecen reflejados en una tabla a registros de intensidad-duracin obtenidos hasta el
la derecha de la pantalla. En el mapa se muestran los lmi- ao 2000 y facilitados por el INM.
tes provinciales y la cuadrcula de las coordenadas UTM Aplicacin de mtodos de regionalizacin intraesta-
para facilitar al usuario la localizacin del punto de su inte- cin para mejorar la robustez de las estimas.
rs. Adems se pueden realizar dos zoom.
Utilizacin de una funcin de distribucin, SQRT-ET-
Una vez localizado el punto en estudio y obtenidos los pa-
mx desarrollada especficamente para el anlisis de
rmetros necesarios, se presiona calcular los valores de
precipitaciones mximas.
I(d,T); aparece otra pantalla, y los parmetros debern in-
troducirse en las celdas correspondientes en la hoja de cl- Aplicacin de una nueva metodologa para la generali-
culo preparada para tal efecto. Automticamente se actua- zacin de resultados obtenidos en estaciones pluvio-
liza la tabla en la parte inferior de la pantalla que calcula la grficas al resto del territorio peninsular.
intensidad mxima para distintas duraciones y perodos de Utilizacin de un SIG que ha permitido el desarrollo
retorno en dicho punto (Figura n6). Se puede obtener ade- de una aplicacin sencilla y fcil de utilizar para el
ms la precipitacin mxima diaria y en 24 horas siguiendo usuario.
la metodologa propuesta por el Ministerio de Fomento
(1999). En la hoja de clculo aparece tambin una celda bajo
el ttulo F; se trata del factor de mayoracin que permite 5. REFERENCIAS
pasar de precipitaciones mximas diarias a mximas en 24 BACHRO, J-N. AND CHAOUCHE,A (2006): Incertitude Des-
horas, y ser el usuario el que elija el factor. A falta de infor- timation des Pluies Extrmes du Pourtour Mditerranen
macin, se puede aceptar un valor orientativo de 1,13 pro- Illustration par les Donnes de Marseille. Hydrological Scien-
puesto por la Organizacin Mundial de Meteorologa. ces-Journal-des Sciences Hydrologiques, 51(3), pp 389-403.
Tambin se ofrece una gua de usuario donde, adems de
indicar los pasos para el uso de la aplicacin, se describe y BELL, F. C., (1969): Generalized rainfall-duration frequency
justifica brevemente la metodologa desarrollada para obte- relationship. Journal Hydraulic Div. ASCE 95(1), pp 311-
ner las ecuaciones resultantes. 327.
DE SALAS REGALADO, L. (2005): Regionalizacin de leyes
IDF para el uso de modelos hidrometeorolgicos de estima-
4. CONCLUSIN cin de caudales. Tesis Doctoral. Universidad Politcnica de
La aplicacin MAXIN permite estimar en el territorio penin- Madrid.
sular espaol, la intensidad mxima anual para distintas DE SALAS REGALADO, L y FERNNDEZ YUSTE, J.A.:
duraciones y perodos de retorno. La estimacin se realiza Annual maximum rainfall intensity intra-station regiona-
aplicando una metodologa que aporta las siguientes modifi- lization to estimate an IDF law: a solution to scarce spatial
caciones respecto a anteriores trabajos: data in Spain. Hydrological Processes ( en prensa).

142 Ingeniera Civil 146/2007

13.MAXIN 14/6/07 14:51 Pgina 143

MAXIN: APLICACIN SIG PARA LA ESTIMACIN DE VALORES INTENSIDAD-DURACIN-FRECUENCIA DE PRECIPITACIONES EN LA ESPAA PENINSULAR

FERRER, F. J., Y L. ARDILES (1994): Anlisis estadstico PORRAS, P. J. Y J. PORRAS (2001): New Perspective on
de las series anuales de mximas lluvias diarias en Espaa. Rainfall Frequency Curves. Journal of Hydrologic Enginee-
Ingeniera Civil/95, pp 87-100. ring, January/February, pp82-85.
FERRERI, G.B. Y V. FERRO (1990): Short Duration Rain- MINISTERIO DE FOMENTO.DIRECCIN GENERAL DE
falls in Sicily. Journal of Hydraulic Engineering 116, 3, pp CARRETERAS (1999): Mximas lluvias diarias en la Es-
430-435. paa Peninsula, 28pp 2h., 25 h pleg de map. Centro de Pu-
HERSHFIELD, D.M. (1961): Rainfall frequency Atlas of the blicaciones. Madrid
United States for durations from 30 minutes to 24 hours MINISTERIO DE OBRAS PBLICAS, TRANSPORTE Y
and return periods from 1 to 100 years. Tech. Paper 40, U.S. MEDIO AMBIENTE. Direccin General de Carreteras
Weather Bureau, Washington, D.C. (1978): Clculo hidrometeorolgico de caudales mximos en
HOSKING, J.R.M. Y J.R.WALLIS (1997): Regional Fre- pequeas cuencas naturales. Secretara de Publicaciones.
quency Analysis. Cambridge University Press. 224 pp. Madrid, 123 pp.
LLASAT, M. C.(1998): Una clasificacin de los episodios plu- MINISTERIO DE OBRAS PBLICAS, TRANSPORTE Y
viomtricos para su utilizacin en hidrologa. Aplicacin a la MEDIO AMBIENTE. Direccin General de Carreteras
serie de intensidad de lluvia en Barcelona. Ingeniera Civil (1987): Clculo hidrometeorolgico de caudales mximos en
112, pp 35-44. pequeas cuencas naturales. Alanmer, Madrid, 123 pp.

Ingeniera Civil 146/2007 143