Geostatistical simulation of spatial variability of convective storms in Mexico City Valley
Javier Méndez-Venegas,1* Martín A. Díaz-Viera,2 Graciela S. Herrera1 and Arturo Valdés-Manzanilla3
1 Instituto de Geofísica, Universidad Nacional Autónoma de México, Ciudad Universitaria, Delegación Coyoacán, 04510 México D.F., México. *Corresponding author: lemendez84@yahoo.com.mx
2 Programa de Recuperación de Yacimientos, Instituto Mexicano del Petróleo, Del. Gustavo A. Madero.
3 División Académica de Ciencias Biológicas, Universidad Juárez Autónoma de Tabasco, Villahermosa, Tabasco.
]]>Received: November 11, 2011
Accepted: January 15, 2013
Published on line: March 22, 2013
Resumen
La precipitación es uno de los factores principales del ciclo hidrológico y el conocimiento de su distribución espacial es fundamental para la predicción de otras variables ambientales íntimamente relacionadas como son: el escurrimiento, las inundaciones, la recarga de los acuíferos. La mayor parte de la precipitación en la Ciudad de México es producida por tormentas convectivas, caracterizadas por una alta variabilidad espacial, lo cual implica que la modelación de su comportamiento sea muy compleja. En el presente estudio se aplicaron técnicas de simulación estocástica con enfoque geoestadístico para modelar la variabilidad espacial de la precipitación de tres tormentas convectivas. El análisis de los resultados muestra que usando la metodología propuesta se obtienen distribuciones espaciales de lluvia que reproducen las características estadísticas presentadas en la información disponible.
Palabras clave: geoestadística, variabilidad espacial de la precipitación, simulación secuencial Gaussiana, cosimulación, tormentas convectivas, radar meteorológico.
Abstract
]]> Precipitation is one of the main components of the hydrological cycle and knowledge of its spatial distribution is fundamental for the prediction of other closely related environmental variables, for example, runoff, flooding and aquifer recharge. Most of the precipitation in Mexico City is due to convective storms characterized by a high spatial variability, implying that modeling its behavior is very complex. In this work stochastic simulation techniques with a geostatistical approach were applied to model the spatial variability of the rainfall of three convective storms. The analysis of the results shows that using the proposed methodology spatial distributions of rain are obtained that reproduce the statistical characteristics present in the available information.Key words: geostatistics, rainfall spatial variability, sequential Gaussian simulation, cosimulation, convective storms, meteorological radar.
Introduction
One of the most modern instruments to estimate rainfall is meteorological radar. It cover a large area (about 200 km in radius); although the estimates are not precise, due to inherent errors of the instrument itself: anomalous propagation, attenuation, etc.; to its surroundings: beam blocking due to mountains, false echoes, evaporation, etc. (Zawadzki, 1984); and to the estimation algorithms (Seo and Krajewski, 2011). The rain gauge has been the traditional instrument for rainfall estimation due to its good precision, though representativeness of its measurements is of a few meters around the instrument. Many countries of the world, to take advantage of both instruments, have systems that estimate rainfall based on a combination of meteorological radar and rain gauge estimates. However, rainfall estimation becomes very complicated when the spatial distribution is very variable, which is the case of convective or electrical storms. Various geostatistical techniques to estimate rainfall using Kalman filters (Anhert et al., 1986) or geostatistical estimation methods such as kriging (Krajewski, 1987) have been developed.
Since 1995 a network of 13 C-band Doppler meteorological radar equipment exists in Mexico (Valdés-Manzanilla and Aparicio, 1997). Their main objectives are to monitor the tropical cyclones in or near the Mexican national territory and to estimate rainfall with hydrological purposes. One of these radar stations is near the metropolitan area of Mexico City, at the top of Cerro de la Catedral. Also, a network of digital rain gauges with telemetry, owned by the government of Mexico City, covers much of the city area.
Because of that, Valdés-Manzanillaand Herrera (2002) designed arainfall estimation method using both sources of meteorological information. A Kalman filter was used to calculate optimally, in real time, the mean error between rainfall estimated by radar and the one estimated by rain gauges. After applying this technique to two convective storms, the root mean square error was reduced by 1.3 and 1.9 mm during the entire storm.
Díaz-Viera, et al. (2009) explored different variants of kriging to estimate rainfall in the Mexico City metropolitan areausing radarand rain gauge data. Their estimates obtained by cokriging with a model of linear corregionalization and collocated cokriging generated better estimates of the rainfall than obtained by ordinary kriging.
Becerra-Soriano (2009), in her master's thesis, continued these two investigations. Her objective was to evaluate the cokriging method for estimating rainfall combining radar and rain gauges measurements and using all radar images from two storms in the Mexico City area. As part of the assessment, a calculation of runoff volume was included, in order to estimate the water volume that would go into the Mexico City drainage system.
Geostatistical estimation techniques like kriging-thebest unbiased linearestimators (Chilès and Delfiner, 1999)-may be optimal in the sense of minimizing the estimation error variance, but are strongly dependent on data quantity, spatial position and, the worst, they do not reproduce the spatial correlation. These techniques can generate unrealistic rainfall spatial distributions (Young, 2008; and Curtis and Clyde, 1999).
]]> An alternative method for spatial estimation is a simulation approach, which, by definition, reproduces the statistical behavior of the phenomenon. Specifically, geostatistical simulation methods can generate multiple realizations that are statistically equivalent in terms of the first and second-order moments (Chilès and Delfiner, 1999). Here, the application of geostatistical simulation methods to model rainfall spatial variability is considered.An antecedent to the present work is the master's thesis of Méndez-Venegas (2008), where he performed a simulation using only rain gauge data and a cosimulation using both rain gauge data and radar images for a single storm. The applied simulation method was sequential Gaussian (Alabert and Massonat, 1990). This paper is an extension of his work to a set of three convective storms in Mexico City.
Here, two simulations for each storm: a univariate simulation (ZS) using only rain gauge data and a cosimulation (ZCS) using rain gauge data and radar images are performed. To evaluate the results, their statistics were satisfactorily compared with those of the data.
Rain gage and radar image data
The radar data was obtained at the C-band Doppler meteorological radar station of the National Meteorological Service on Cerro de la Catedral, overlooking the metropolitan area of Mexico City (Figure 1). The radar images used are 8 bits images of 240 x 240 km with a resolution of 1 km2. A pixel covers an area of 1 km x 1 km in a pseudo-CAPPI presentation at 4 km above sealevel everyfifteenminutes(Valdés-Manzanilla and Aparicio, 1997). The precipitation datais from 61 rain gauges of the Water System of Mexico City and radio reporting, every minute, of the accumulated rainfall during the storm. These rain gauges are of tipping-bucket kind with telemetry and a density of one rain gauge for every 30 km2 (Díaz-Viera, et al., 2009).
The accumulated rainfall per hour for each type of measurement is calculated. The radar records an image with values of reflectivity (Z) every 15 minutes, these images to values of rain intensity (R) using a Z-R relationship are converted, subsequently four consecutive radar images are averaged for rainfall intensity to obtain effective cumulative rainfall in one hour. The relation Z = 300R1.4 recommended by the manufacturer is used (Valdés-Manzanilla and Herrera- Zamarrón, 2000).
For rain gauge data digitalfiles for the date and time of the storm are used. Each rain gauge has a counter that is incremented by one each time it registers a shower of 1/4 mm (Rosengaus, 2000). The cumulative rainfall per hour is calculated.
Rain gauges Zg (Table 1) and radar images Zr (Figure 2, Figure 3, Figure 4) were recorded in Mexico City during 13, 15 and 16 July 1997 (referred thereafter as storm 1, 2 and 3). Storm 2 has the largest number of gauge data (50), while storm 1 has only23gaugemeasurements.On Figure 2, Figure 3 and Figure 4 the gray scale images correspond to one hour accumulated precipitation given in millimeters (mm), while the cross symbols represent the locations of gauge data for this storm.
]]>
Geostatistical simulation methodology
Geostatistical methodology basically consists of three phases: exploratory data analysis, variographic analysis and estimation and/or simulation. The geostatistical simulationis applied in this work. Standard procedures are followed (Díaz-Viera, et al., 2009; and Méndez-Venegas, 2008).
Sequential Gaussian Simulation
Since the early 1990's, sequential Gaussian simulation has gained in popularity (Deutsch, 2002). A new simulated valueis obtained fromthe estimated conditional probability distribution function using observational and previously simulated values in a neighborhood of a given location applying a kriging method (Chilès and Delfiner, 1999).
]]> The theory behind sequential Gaussian simulation is based on using previously simulated value and input data throughout the simulation process. In practice, only the closest conditioning data are used.
Geostatical simulation of three storms
The sequential Gaussian method was applied using the rain gauge data and radar data. During the exploratory data analysis, several statistical parameters were computed (Table 1) and histogram graphics were generated. It was found that the data do not have normality; consequently, an anamorphosis transformation was applied to them, which ensured normality in the transformed data (Chilès and Delfiner, 1999).
Variograms were calculated and a model was adjusted to each using weighted least squares. The model with the lowest sum of squares errors was chosen and validated using cross validation. The leave-one-out method (Journel and Huijbregts, 1978) was used for cross-validation which involves removing each sample and estimating the value at that point using the kriging equations and the variogram model obtained. As a result, a map of the differences between actual and estimated values is obtained. The adjusted models were spherical.
Results and discussion
For the simulation of storm 1 (Figure 5) rain gauge data and the model shown in the first line of Table 4 were used. For the cosimulation of this storm (Figure 6) rain gauge and radar data of the corresponding storm and the model showed in Table 2 were used. The simulation of storm 2 (Figure 8) was done using rain gauge data and the model in the third line of Table 4; for the cosimulation (Figure 9) the rain gauges and radar data of storm 2 and the model shown in Table 7 were used. For storm 3, as in the previous two cases, the univariate simulation (Figure 11) only uses rain gauge data and a model (fifth line, Table 4) and the cosimulation (Figure 12) was done with all the information available and amodel (Table 8).
Results for each storm were compared with the corresponding sample information. The Figure 6. Cosimulation of one hour accumulated rainfall in millimeters for storm 1. simulations reproduce adequately the statistical values (Table 3, Table 6 and Table 9), the histograms and box plots (Figure 7, Figure 10 and Figure 13), as the variogram model of the data (Table 4 and Table 5).
For all storms, simulations estimate adequately precipitation. Comparison of radar images with simulations is complex because radar images are not taken at ground level. In radar images the precipitation falls in a small region and is greater than those recorded by rain gauges and from simulations, this may be due to evaporation of rainfall before it reaches the ground, because of high elevation of the radar beam over the valley of Mexico (Zawadzki,1984).
Conclusions
Spatial stochastic simulations using the sequential Gaussian method reproduce adequately data statistics (minimum, maximum, mean value, median, variance, histogram, variogram model, etc.) in both univariate and bivariate cases. Simulationcanbeanideal tool tomodel thespatial distribution of rainfall.
When there is enough information to accurately estimate the variogram (Storm 2), simulations using only rain gauge data generated consistent estimations with the variability and the spatial distribution of the rainfall, but cosimulations with rain gauge data and radar images generated more precise and detailed estimations of the spatial distribution.
Using the simulation approach, rainfall distributions in storms could be generated from their statistical properties. Simulation is a powerful tool for studying the phenomena involved in precipitation.
For optimal performance of the simulation procedures, it is necessary to follow a methodology consistent with hypothesis, as those described here.
]]> Bibliography
Anhert P., Krajewski W., Johnson E., 1986, Kalman Filter estimation of radar-rainfall fild bias. XXIII Conferencia en meteorología de radar. Americ. Meteor. Soc., Snowmass, 33-37. [ Links ]
Becerra-Soriano L., 2009, Estimación de lluvia en el Distrito Federal utilizando datos de pluviógrafos y de radar meteorológico, tesis de Maestría, UNAM. [ Links ]
Chilès P.J., Delfiner P., 1999, Geoestatistics: Modeling Spatial Uncertainty. Wiley. New Cork. 695 pp. [ Links ]
Collier C.G., 1983, A weather radar procedure for real-time procedure estimation of surface rainfall. Q.J.R.M.S., 109, 589-608. [ Links ]
Curtis D.C., Clyde B.S., 1999, Comparing Spatial Distributions of Rainfall Derived from Rain Gages and Radar. NEXRAIN Corporation Folsom. [ Links ]
Deutsch C., 2002, Geostatistical Reservoir Modeling. Oxford University Press, New York. [ Links ]
Díaz-Viera M., Herrera-Zamarrón G.S., Valdés-Manzanilla A., 2009, A linear coregionalization model for spatial rainfall estimation in the Mexico City valley combining rain gages data and meteorological radar images. Revista Ingeniería Hidráulica en México, vol. XXIV, No. 3. pp. 63-90. Julio-septiembre 2009. [ Links ]
Fitzwilliams P., Rios T., Curtis D., Thornhill R., 2006, Use of Radar-Rainfall in GIS-Based Sewer Modeling. Government Engineering. [ Links ]
Journel A.G., Huijbregts Ch.J., 1978, Mining Geostatistics. Academic Press. New York, 590 pp. [ Links ]
Krajewski W.F., 1987, Cokriging radar-rainfall and rain gage data, J. Geophys,, 92. 9571-9580. [ Links ]
Matheron G. 1963. Principles of Geostatistics. Economic Geology, 58, 1246-1266. [ Links ]
Méndez-Venegas J., 2008, Modelación de la Distribución Espacial de la Precipitación en el Valle de la Ciudad de México Usando Técnicas Geoestadísticas, tesis de Maestría en Estadística, Colegio de Postgraduados, Campus Montecillo, Chapingo. [ Links ]
Rosengaus M., 2000, Manejo de Emergencias Hidrometeorológicas en la Ciudad de México. Primer Simposio Internacional Sobre Riesgos Geológicos y Ambientales de la Ciudad de México. México D.F. [ Links ]
Seo B.C., Krajewski W.F., Investigation of the scale-dependent variability of radar-rainfall and rain gauge error covariance. Advances in Water Resources, 34, 1, January 2011, 152-163. [ Links ]
Seo D.J., Krajewski W.F., Bowles D.S., 1990a, Stochastic interpolation of rainfall data from rain gages and radar using Cokriging. 1. Design of experiments. Water Resources Research, 26, 3, 469-477. [ Links ]
Seo D.J., Krajewski W.F., Bowles D.S., 1990b, Stochastic interpolation of rainfall data from rain gages and radar using Cokriging. 2. Results. Water Resources Research, 26, 5, 915-924. [ Links ]
Valdés-Manzanilla A., Aparicio F.J., 1997, The Mexican Doppler radar network. XVIII Conference of radar meteorology, Austin Tx, Amer. Meteor. Soc. 35-36. [ Links ]
Valdés-Manzanilla A., Herrera-Zamarrón G., 2000, Informe final del proyecto: diseño de un sistema de estimación de lluvia usando radar meteorológico. Jiutepec, México: Coordinación de Tecnología Hidrológica, Subordinación de Hidrometeorología, IMTA. [ Links ]
Valdés-Manzanilla A., Herrera G.S., 2002, Design of a rain estimation system using a meteorological radar. Developments in water science. Computational methods in water resources, 47, 2, 1765-1772. [ Links ]
Young H.B., Byung K.S., 2008, Radar Rainfall Adjustment by Kalman-Filter Method and Flood Simulation using Two Distributed Models. The fifth European conference on radar in meteorology and hydrology. [ Links ]
Zawadzki I., 1984, Factors affecting the precision of radar measurements of rain. Preprints of the 22nd. Conference on radar meteorology. AMS, Boston, Mass. 251-256. [ Links ]
Zhanga Z., Switzerb P., 2007, Stochastic space-time regional rainfall modeling adapted to historical rain gauge data. American Geophysical Union. [ Links ]
]]>