Statistical Prediction of Wet and Dry Periods in theComahue Region (Argentina)
Marcela H. González, Diana Dominguez
DOI: 10.4236/acs.2012.21004   PDF    HTML     5,530 Downloads   9,234 Views   Citations


General features of rainy season with excess or deficits are analyzed using standardized precipitation index (SPI) in Limay and Neuquen River basins. Results indicate that most of dry and wet periods persist less than three months in both basins. Furthermore, an increase of rainfall variability over time is observed in the Limay river basin but it is not detected in the Neuquen river basin. There is a tendency for wet (dry) periods to take place in El Ni?o (La Ni?a) years in both basins. Rainfall in both basins, have an important annual cycle with its maximum in winter. In addition, possible causes of extreme rainy seasons over the Limay River Basin are detailed. The main result is that the behavior of low level precipitation systems displacing over the Pacific Ocean in April influences the general hydric situation during the whole rainy season. In order to establish the existence of previous circulation patterns associated with interannual SPI variability, the composite fields of wet and dry years are compared. The result is that rainfall is related to El Ni?o- Southern Oscillation (ENSO) phenomenon and circulation over the Pacific Ocean. The prediction scheme, using multiple linear regressions, showed that 46% of the SPI variance can be explained by this model. The scheme was validated by using a cross-validation method, and significant correlations are detected between observed and forecast SPI. A polynomial model is used and it little improved the linear one, explaining the 49% of the SPI variance. The analysis shows that circulation indicators are useful to predict winter rainfall behavior.

Share and Cite:

M. González and D. Dominguez, "Statistical Prediction of Wet and Dry Periods in theComahue Region (Argentina)," Atmospheric and Climate Sciences, Vol. 2 No. 1, 2012, pp. 23-31. doi: 10.4236/acs.2012.21004.

1. Introduction

The Andes Mountain range lies all along the western part of South America; and the Comahue region is located in that area, between 38˚S and 43˚S (Figure 1). Two important rivers—Limay and Neuquen-run in this area. The Alicura, Piedra del Aguila, Pichi Picun Leufu and El Chocón hydroelectric dams are on these rivers. The hydropower system is approximately 5000 Mw with an annual energy generation of 14,500 Gwh, 20% of the Argentinean budget. While operating the dams, there are several aspects to be taken into account regarding water level, such as emergency, flood smoothing, normal and extraordinary levels, probability of flooding and high drainage level from irrigated valleys. This conventional work operation level is not enough because the prospective demand of water depends not only on the meteorological situation but also on daily electricity demand. Therefore, it is necessary a better understanding of rainfall behaviour over the basin through the knowledge of atmospheric predictors which allows to anticipate seasonal precipitation and the building of statistical models. The scientific basis of seasonal climate predictability lies in the fact that slow variations in the earth boundary conditions (i.e. sea surface temperature or soil wetness) can influence global atmospheric circulation, and thus precipitation. As the skill of seasonal numerical prediction models is still limited, the statistical study of the probable relationships between some local or remote forcing and rainfall is essential. Some authors have analyzed these relationships in the southern hemisphere-Gissila et al. [1] in Ethiopia; Reason [2] in South Africa and Zheng and Frederiksen [3] in Australia. In Argentina, Gonzalez and Vera [4] detected rainfall patterns analyzing interannual rainfall variability in the Comahue region using a principal component analysis. Gonzalez et al. [5] derived rainfall prediction schemes, using multiple linear regressions which explained the 51% of the winter rainfall variance in the Limay River basin and the 44% in the Neuquen River one; and Gonzalez and Cariaga [6] use predictors for the application of the Climate Prediction Tool (CPT) from IRI. CPT software is based on canonical correlation analysis, and the authors showed that correlation between observed and forecast winter rainfall is significant all over the area of study and increases towards the west and the northwest. Gonzalez and Murgida

Figure 1. Stations used in the study.

[7] studied the characteristics of rainfall in the Chaco region of Argentina and detected previous circulation patterns associated with wet summer seasons. Some studies focused on studying the relation with El Niño-Southern Oscillation (ENSO) [8-14]. The aim of this paper is to detect the possible relationships between the accumulated rainfall during rainy season in the Comahue region and both, the ocean and atmospheric circulation patterns, previously observed. The paper is organized as follows: Section 2 describes the dataset and the methodology; Section 3 presents the SPI features in Limay (LB) and Neuquen (NB) river basins. Section 4 shows the relation between SPI calculated over the period April to September and atmospheric circulation. The association with SST was detailed in Section 5. Section 6 presents the building of regression models to estimate wet and dry events and Section 7 presents the main conclusions.

2. Data and Methodology

The area under study includes two river basins: the Limay river basin (LB) in the south and the Neuquen river basin (NB) in the north (Figure 1). To carry out this study, there were used monthly rainfall data derived from 20 stations at different sources—the National Meteorological Service, the Secretary of Hydrology of Argentina and the Territory Authority of the Limay, Neuquén and Negro river basins and the 1980-2007 record thereof. The period was selected because all the stations have less than 20% of missing monthly rainfall data and their quality has been carefully proved.

Two mean rainfall series were calculated to obtain the average of monthly precipitation of twelve stations in the NB and eight stations in the LB, in order to be representative of the precipitation over each one of the basins (Figure 1). The standardized precipitation index (SPI) was used to quantify the conditions of deficit or excess of precipitation for a six-month interval, which is the period that better works for hydrological matters [15,16], and with the purpose of improving the detection of rainfall excesses and deficits which could have significant consequences for the operation of dams. Computation of the SPI involves fitting a gamma probability density function to a given frequency distribution of six-month precipitation totals for LB and NB [15]. The SPI greater (lower) than zero indicates water excess (deficit). A wet (dry) period is defined as the period during which the SPI is continuously positive (negative). The magnitude of the index allows classifying the six-month accumulated rainfall in categories that go from extreme drought to extreme excess (extremely dry, severe dry, moderate dry, normal, moderate wet, severe wet and extremely wet). SPI series low frequency variability was analyzed using a linear trend method of minimum squares, and statistics significance was tested using a T-Student test.

There are also used monthly sea surface temperatures (SST), 500 HPa (G500), 1000 Hpa (G1000) and 200 Hpa (G200) geopotential heights, zonal (U) and meridional (V) winds at 850 Hpa and precipitable water (PW) from National Center of Environmental Prediction (NCEP) reanalysis [17]. Monthly anomalies were determined removing the climatologically monthly means from the original values. Composite fields of these variables for wet and dry years and the difference between them are plotted. The statistical significance is checked using a Student’s t-test (95% confidence level) on the difference between the wet and dry means. These significant areas are used to find the existing relation between SPI calculated with accumulated rainfall from April until September (SPI9) in LB and SST, G1000, G500, G200, U and V during the same period and only in April. The results permit to define some predictors, which are used to develop a statistical forecast model using the forward stepwise regression method [18] which retained only the variables correlated with a 95% significance level. Forward stepwise regression is a model-building technique that finds subsets of predictor variables that most adequately predict responses on a dependent variable by linear regression, given the specified criteria for adequacy of model fit [19]. Predictors available to carry out the regression scheme are carefully selected, based on statistical significance and physical reasoning. A polynomial model is also derived using a standard regression method.

The skill of the schemes was proved using contingency tables and a Chi-square test was applied to detect if they differ significantly from random ones. Some measures of accuracy [18] were calculated for wet and dry cases. The hit rate (H) or the right proportion is the fraction of all the cases when the categorical forecast correctly anticipated the subsequent event. The probability of detection (POD) is defined as the fraction of those occasions when the forecast event occurred on which it was also forecast. The false alarm relation (FAR) is the proportion of forecast events that fail to happen. Additionally, empirical estimated and observed SPI9 probability functions were calculated using frequency distributions and a chi-square test was used to prove that they do not differ significantly.

3. SPI Features in LB and NB

Figure 2 shows the six-months SPI for LB and NB series and their respective linear trend adjustment. It is notable that there is no trend in both cases. The duration of each period is defined as the number of months that the SPI retains the same sign and the intensity is the mean SPI over each period. Thirty wet and thirty dry periods were detected in SPI series in LB. Table 1 shows that 53% of the wet periods persists less than 3 months and 23% more than 10 months, meanwhile 60% of dry periods lasts less than 3 months and only 3% lasts more than 10 months, indicating that most of the dry periods are short. Table 2 shows the number, mean duration and mean intensity of wet and dry periods in the sub-periods 1980- 1988; 1989-1997 and 1998-2006. Wet and dry periods during 2007 are disregarded (only in this table) in order for the 3 periods to have the same length. It can be seen that the number of wet periods in LB slightly increases and their duration is shorter and the number of dry periods increases but their intensity and duration remain similar. This result indicates an increase of rainfall variability over time.

Twenty-four wet and twenty-four dry periods are detected in SPI series in NB. Table 1 shows that 37.5% of wet periods lasts more than 10 months and 29% less than 3 months meanwhile most of dry periods are shorter (less

Figure 2. Twelve months SPI for Limay River Basin (upper pannel) and Neuquen River Basin (bottom pannel).

Table 1. Percentage of wet and dry periods with different duration for LB (upper panel) and NB (bottom panel).

Table 2. Number, mean intensity and mean duration of wet and dry periods in 1980-1988, 1989-1997 and 1998-2006 in LB (upper panel) and NB (bottom panel).

than 3 months). Table 2 reveals that there is no specific trend when different time periods are considered.

It is important to notice that there is a tendency for wet (dry) periods to take place in El Niño (La Niña) years in both basins. In fact, in four wet periods in LB—out of the nine wettest ones—the warm phase on ENSO was in its mature stage, meanwhile in two of them La Niña event was developing. But in three dry periods—out of the eight driest ones—La Niña event was present and in two dry periods El Niño events developed. The signal is clearer in NB where, in five—out of the six wettest periods—El Niño was present and La Niña developed in none of them and in three—out of the six driest periods—La Niña was present and El Niño developed in none of them. Therefore, the correlation between sixmonth SPI and SST in EN34 region in the last month of the period is 0.38 in LB and 0.37 in NB, both significant at the 95% confidence level. This particular relation will be described in detail in Section 5.

Figure 3 shows the mean (1980-2007) monthly rainfall in both basins. It reveals that precipitation has an important annual cycle with a maximum in winter. As the most important quantity rainfall amount occurs from April to September, the analysis of SPI9 (SPI for the precipitation accumulated in this period) will be detailed hereinafter. Results in LB and NB were similar but only LB ones will be presented in this paper.

The percentage frequency of SPI9, a summary of the definitions and the years when the different categories took place are presented in Table 3.

4. Relation between SPI9 and the Atmospheric Circulation

Figure 4 shows the G1000 anomaly composite fields for wet and dry cases during the period April-September and the difference between wet and dry composites, and Figure 5 shows the same only for April. The most important feature is the dipole between high and middle latitudes (Figure 4). In dry composites the subtropical heights and the sub-polar lows are intensified and the subtropical heights extend towards Argentina. Meanwhile, in wet cases both belts are weakened. This pattern remains when composites in April are considered (Figure 5). This is a relevant result because the April geopotential height field could be used for predicting purposes. The same pattern is observed in the composites for each one of individual months in the period May-September (figures not shown). It is important to notice that precipitation systems usually displace over the Pacific Ocean, arrive in the South American coast and can go through the southern portion of the Andes mountains (south of 38˚S), because the latter are lower than in the north. Therefore

Conflicts of Interest

The authors declare no conflicts of interest.


[1] T. Gissila, E. Black, D. I. F. Grime and J. M. Slingo, “Seasonal Forecasting of the Ethiopian Summer Rains,” International Journal of Climatology, Vol. 24, No. 11, 2004, pp. 1345-1358. doi:10.1002/joc.1078
[2] C. Reason, “Subtropical Indian Ocean SST Dipole Events and Southern Africa Rainfall,” Geophysical Research Letters, Vol. 28, No. 11, 2001, pp. 2225-2227. doi:10.1029/2000GL012735
[3] X. Zheng and C. Frederiksen, “A Study of Predictable Patterns for Seasonal Forecasting of New Zealand Rainfall,” Journal of Climate, Vol. 19, No. 13, 2006, pp. 3320-3333. doi:10.1175/JCLI3798.1
[4] M. H. González and C. S. Vera, “On the Interannual Winter Rainfall Variability in Southern Andes,” International Journal of Climatology, Vol. 30, No. 5, 2010, pp. 643-657
[5] M. H. González, M. M. Skansi and F. Losano, “A Statistical Study of Seasonal Winter Rainfall Prediction in the Comahue Region (Argentine),” Atmosfera, Vol. 23, No. 3, 2010, pp. 277-294.
[6] M. L. Cariaga and M. H. González, “Estimating Winter and Spring Rainfall in the Comahue Region (Argentine) Using Statistical Techniques,” Advances in Environmental Research 11, Justin A. Daniels, Nova Science Publishers Inc., New York, 2011, pp. 103-128.
[7] M. H. González and A. M. Murgida, “Seasonal Summer Rainfall Prediction in Bermejo River Basin in Argentina,” Climate Variability, Abdel Hannachi, Intech, Reading, 2011, in Press.
[8] C. Schneider and D. Gies, “Effects of El Ni?o-Southern Oscillation on Southernmost South America Precipitation at 53oS Revealed from NCEP-NCAR Reanalysis and Weather Station Data,” International Journal of Climatology, Vol. 24, No. 9, 2004, pp. 1057-1076. doi:10.1002/joc.1057
[9] P. Aceituno, “On the Functioning of the Southern Oscillation in the South American Sector. Part I: Surface Climate,” Monthly Weather Review, Vol. 116, No. 3, 1988, pp. 505-524. doi:10.1175/1520-0493(1988)116<0505:OTFOTS>2.0.CO;2
[10] P. Aceituno and R. Garreaud, “Impactos de Los Fenó- menos El Ni?o y La Ni?a Sobre Regímenes Fluvio- métricos Andinos,” Revista de la Sociedad Chilena de Ingeniería Hidráulica, Vol. 10, No. 2, 1995, pp. 33-43.
[11] J. Rutland and H. Fuenzalida, “Synoptic Aspects of the Central Chile Rainfall Variability Associated with the Southern Oscillation,” International Journal of Clima- tology, Vol. 11, No. 1, 1991, pp. 63-76. doi:10.1002/joc.3370110105
[12] P. Aceituno, “Seasonality of the ENSO Related Rainfall Variability in Central Chile and Associated Circulation Anomalies,” Journal of Climate, Vol. 16, No. 2, 2003, pp. 281-296. doi:10.1175/1520-0442(2003)016<0281:SOTERR>2.0.CO;2
[13] R. Compagnucci and W. Vargas, “Inter-Annual Variability of the Cuyo River Streamflowin the Argentinean Andean Mountains and ENSO Events,” International Journal of Climatology, Vol. 18, No. 14, 1998, pp. 1593-1609. doi:10.1002/(SICI)1097-0088(19981130)18:14<1593::AID-JOC327>3.0.CO;2-U
[14] R. Compagnucci and D. Araneo, “Identificación de áreas de Homogeneidad Estadística Para los Caudales de los Ríos Andinos Argentinos y su Relación Con la Circu- lación Atmosférica y la Temperatura Superficial del Mar,” Meteorologica, Vol. 30, No. 1, 2005, pp. 41-54.
[15] T. B. Mckee, N. Doesken and J. Kleist, “The Relationship of Drought Frequency and Duration to Times Scales,” Proceeding 8th Conference of Applied Climatology, Anaheim, 17-22 January 1993.
[16] T. B. McKee, N. Doesken and J. Kleist, “Drought Moni- toring with Multiple Time Scales,” Proceeding 9th Con- ference on Applied Climatology, Dallas, 15-20 January 1995.
[17] E. Kalnay, M. Kanamitsu, R. Kistler, W. Collins, D. Deaven, L. Gandin, M. Iredell, S. Saha, G. White, J. Woollen, I. Zhu, M. Chelliah, W. Ebisuzaki, W. Higgings, J. Janowiak, K. C. Mo, C. Ropelewski, J. Wang, A. Leetmaa, R. Reynolds, R. Jenne and D. Joseph, “The NCEP/NCAR Reanalysis 40 Years-Project,” Bulletin of the American Meteorological Society, Vol. 77, 1996, pp. 437-471. doi:10.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2
[18] D. S. Wilks, “Statistical Methods in the Atmospheric Sciences (An Introduction),” International Geophysics Series Academic Press, San Diego, 1995.
[19] R. B. Darlington, “Regression and Linear Models,” McGraw- Hill, New York, 1990.
[20] E. Kalnay, K. C. Mo and J. Paegle, “Large-Amplitude, Short Scale Stationary Rossby Waves in the Southern Hemisphere: Observations and Mechanistics Experiments on Determine Their Origin,” Journal of the Atmospheric Sciences, Vol. 43, No. 3, 1986, pp. 252-275. doi:10.1175/1520-0469(1986)043<0252:LASSSR>2.0.CO;2
[21] NOAA/ESRL Physical Sciences Division, Boulder Colorado, 2011.
[22] H. Panofsky and G. Brier, “Some Applications of Sta- tistics to Meteorology,” The Pennsylvania State University, Pennsylvania, 1965.

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.