Were the 2010 Pakistan floods predictable?
P. J. Webster,
- School of Earth and Atmospheric Science, Georgia Institute of Technology, Atlanta, Georgia, USA
V. E. Toma,
- School of Earth and Atmospheric Science, Georgia Institute of Technology, Atlanta, Georgia, USA
- School of Earth and Atmospheric Science, Georgia Institute of Technology, Atlanta, Georgia, USA
During July 2010, a series of monsoonal deluges over northern Pakistan resulted in catastrophic flooding, loss of life and property and an agricultural crisis that may last for years. Was the rainfall abnormal compared to previous years? Furthermore, could a high probability of flooding have been predicted? To address these questions, regional precipitation is analyzed using three dataset sets covering the 1981–2010 time period. It is concluded that the 2010 average May to August (MJJA) rainfall for year 2010 is somewhat greater in magnitude than previous years. However, the rainfall rate of the July deluges, especially in North Pakistan was exceptionally rare as deduced from limited data. The location of the deluges over the mountainous northern part of the country lead to the devastating floods. The European Centre for Medium Range Weather Forecasts (ECMWF) 15-day Ensemble Prediction System (EPS) is used to assess whether the rainfall over the flood affected region was predictable. A multi-year analysis shows that Pakistan rainfall is highly predictable out to 6–8 days including rainfall in the summer of 2010. We conclude that if these extended quantitative precipitation forecasts had been available in Pakistan, the high risk of flooding could have been foreseen. If these rainfall forecasts had been coupled to a hydrological model then the high risk of extensive and prolonged flooding could have anticipated and actions taken to mitigate their impact.
Two main factors control South Asian rainfall. On 2–5 year time scales, the El Niño-Southern Oscillation (ENSO) phenomena is associated with above average summer precipitation during a La Niña and deficits during an El Niño [Shukla and Paolina, 1983; Kumar et al., 2006]. Far more dramatic and higher amplitude modulations occur on subseasonal time scales. Over much of Asia the summer monsoon is divided into a series of “active” (rainy) and “break” (dry) periods following a roughly 20–40 days cycle [Lawrence and Webster, 2001; Webster and Hoyos, 2004; Hoyos and Webster, 2007] associated with the boreal summer Madden-Julian Oscillation [Madden and Julian, 1972] that produce a northeasterly excursion of large-scale convective anomalies under the action of a strong cross-equatorial pressure gradient [Stephens et al., 2004; Wang et al., 2005, 2006]. The arrival of convection over the Indian subcontinent heralds an active pluvial period. Summer rainfall in Pakistan is also monsoonal and, as such, has active and break periods. However, the total summer rainfall is far less than in the east (Figure 1a) decreasing from the Bay of Bengal (16 mm/day) across the plains of northern India (8–10 mm/day) to values of about 6–8 mm/day in northern Pakistan. Pakistan is at the western edge of the pluvial region of the monsoon.
During the late boreal spring of 2010, the tropical Pacific Ocean entered a La Niña phase and during July 2010 the monsoon over the northern part of the Indian subcontinent was “active” with rainfall extending across the Gangetic Plains between the Bay of Bengal in the east to northern Pakistan in the west (Figure S1 of the auxiliary material). Embedded in this active period were the deluges that caused the devastating floods in Pakistan. In late July, some Pakistan stations recorded rainfall amounts exceeding 300 mm over a four-day period (http://www.pakmet.com.pk/FFD/index_files/rainfalljuly10.htm). During the following days and weeks, flooding extended through the entire Indus Valley eventually reaching the Arabian Sea leaving behind a wake of devastation and destruction. In the end, the death toll was close to 2000 and over 20 million people were affected. An estimated 20,000 cattle were drowned. Power stations and transmission towers were destroyed along with other major infrastructure such as barrages, bridges and roads. Irrigation systems were destroyed and planting of subsequent crops delayed or abandoned with agricultural costs exceeding $US500M. Overall, estimates of damage exceed $US40B. In general, it was the poor that suffered the most and many of these will face the prospect of intergenerational poverty as a result of the floods [Webster and Jian, 2011]. Most assessments of the 2010 Pakistan floods have appeared on the internet and in relief organization reports (http://www.pakistanfloods.pk/; http://en.wikipedia.org/wiki/2010_Pakistan_floods). Eventually, scholarly articles on the flooding will be forthcoming discussing, in more detail, the climate and meteorological conditions that produced the flooding (e.g., R. A. Houze Jr. et al., Anomalous atmospheric events leading to the summer 2010 floods in Pakistan, submitted to Bulletin of the American Meteorological Society, 2010). However, to date there has been an absence of any comment about the predictability of the deluges or the associated risk of floods. Eventually, skill in predicting floods reduces to the predictability of precipitation and the use of an adequately sophisticated hydrological model. Thus, an immediate and critical question is the degree to which rainfall at the western edge of the South Asian monsoon system is predictable on time scales of 1–2 weeks. Is the predictability of precipitation in the western edge of the monsoon comparable to that seen over the Ganges and Brahmaputra basins [Hopson and Webster, 2010; Webster et al., 2010]?
In this study we focus on the predictability of 1–15-day ECMWF EPS forecasts [Buizza et al., 2007] over Pakistan. In the next section details of the observation and numerical model data are introduced. Section 3 discusses the uniqueness of the July-August flooding events and examines the prediction skill of 15-days rainfall forecast followed by conclusions related to the predictability of floods in Pakistan.
2. Data and Analysis
Three precipitation data sets are used to assess the variability of the precipitation over the Pakistan region: the Global Precipitation Climatology Project (GPCP) data [Adler et al., 2003] for the 1981–2009 period, the Tropical Rainfall Measuring Mission (TRMM) [Huffman et al., 2005, 2007] TRMM_3B42 product for 1998–2010, and the NOAA CPC Morphing Technique (CMORPH) Precipitation Product for the 2003–2010 period [Joyce et al., 2004]. GPCP (a merging of rain gauge data with satellite geostationary and low-orbit infrared and passive microwave information) and TRMM data sets (specifically the TRMM_3B42 set) were chosen for their temporal extension (29 and 13 years, respectively). All of these precipitation products had a 0.25° × 0.25° horizontal resolution facilitating a comparison with model output. Figure S4 shows time series of monthly rainfall anomalies for each of the data sets.
A comparison of the CMORPH and TRMM data sets (Figures S3) reveals considerable differences in the magnitude of estimates of precipitation during the third precipitation pulse of July 2010 that occurred over the higher terrain of northern Pakistan (Figure S2d). The TRMM rainfall estimate was considerably higher than CMORPH by about a factor of two consistent with the discussion of Gopalan et al.  who suggested that TRMM may overestimate precipitation rates over substantial terrain. Comparisons during earlier periods, when the precipitation maxima occurred over the plains of southern Pakistan and northwestern India are more comparable (Figure S2). Consequently, we use CMORPH as the principal data set for determining the sequence of events during 2010 and also as the principal agent for the statistical rendering of the quantified precipitation forecasts.
The ECMWF EPS forecasts consist of 51 ensemble members initialized twice per day (00 and 12 UTC), each ensemble member having a 15-day forecast horizon. The horizontal resolution of the model is 50 × 50 km from 0 to 10 days and then 80 km × 80 km from day 10–15 [Buizza et al., 2007]. For this initial study, model forecast precipitation for the months of July and August from 2007 to 2010 was converted into daily cumulative amounts. To minimize systematic model bias differences between the distributions of the ECMWF forecasts and the observed rainfall, a quantile-to-quantile (q-to-q) mapping technique was implemented following Hopson and Webster  and Webster et al.  (see method description in the auxiliary material). All rainfall forecasts presented here are adjusted using the q-to-q technique.
Beginning in early July 2010, there were six major pulses of torrential rainfall occurring over Pakistan, each separated by about a week (Figure 1b). One of the most intense periods occurred between July 27–30 over the mountainous regions of the north. Figure S2 shows the distribution of rainfall for the major pulses of monsoon rain. The earlier rainfall events caused flooding in Balochistan in central Pakistan. Flooding followed across northern Pakistan in the Khyber Pakhtunkhwa province with the later July rains extending to the Punjab in late July/early August (http://www.unitar.org/unosat/node/44/1469). Here we address the uniqueness and predictability of the floods.
There have been 67 reported flooding events in Pakistan occurring since 1900 with a clustering of 52 events of various severity in the last 30–40 years (International Disaster Data Base, http://www.emdat.be). Some of these events (e.g., 1950, 1973, 1976, 1977, 1992, 2001, 2007 and 2008) were also accompanied with large loss of life and property. This recent increase is consistent with the increase in intensity of the global monsoon accompanying the last three decades of general global warming (B. Wang et al., Recent intensification of global monsoon and precipitation, submitted to Nature, 2011) or perhaps changes in water management strategies, increases of damage due to a rapidly growing population or improved reporting through advances in communication.
Figure S4 shows the temporal variability of seasonal (MJJA) precipitation averaged in Pakistan (62°–74°E, 24°–36°N, blue rectangle in Figure 1a) and northern Pakistan (70°E–74°E, 30°N–36°N, red rectangle in Figure 1a) relative to the seasonal climatology for each of the data sets: GPCP and CMORPH. While there are amplitude differences between datasets, each shows substantial variability, with seasons of excessive rainfall and drought occurring irregularly over the past 30 years (Figures S3 and S4).
Were the rainfall events of 2010 worse than previous extreme events? Using a 13-year TRMM precipitation record, extreme events can be counted. An extreme event is defined here to occur when the two-days accumulated rainfall exceeds over 10 mm over all Pakistan and 20 mm over the northern Pakistan (Figures 1c and 1d). Note that the chosen thresholds for this analysis are much smaller than maximum daily rainfall measurements at specific stations (see http://www.pakmet.com.pk/FFD/index_files/rainfalljuly10.htm) due to a broader averaging area. Although there is considerable interannual variability, the number of extreme events over entire Pakistan, so defined, is larger in 2010 than in previous years, greater, for example than in 2008. In summary, 2010 stands out as a period of above average rainfall events over northern Pakistan. The number of extreme events over northern Pakistan is far more unique which, based on the very limited TRMM data set would have return periods of >30 years. Long-term variability for extreme events is calculated with GPCP pentad data set from 1981 (Figure S5) to 2007 overlapped with CMORPH pentad from 2003 to 2010. Although, there are differences between data sets, the high occurrence of Northern Pakistan extreme events in 2010 is relatively rare. Rainfall data is not sufficiently reliable prior to 1987 when GPCP data was generated on a daily basis. However, we do have CMORPH and TRMM data for 2008. As shown in Figure S6, the cumulative July-August rainfall for northern Pakistan is larger in 2010 than 2008, with values larger than 0.5 m in several areas.
The next step is to examine the predictability of the rainfall events depicted in Figure 1b. Figure 2a shows the total average precipitation [mm/day] for July 28–29, based on the CMORPH observational dataset and the ECMWF forecast ensemble mean initialized 4 days before the event (Figure 2b). The q-to-q correction was applied to the precipitation forecast data. The forecasts compare well with the observed rainfall with ECMWF slightly underestimating the rainfall intensity in the northern part of the region. The ECMWF forecast showed average precipitation larger than 40 mm/day in some areas which is over 3 times larger than the CMORPH climatological average for the region.
Figures 2c and 2d show the temporal evolution of the ECMWF forecast commencing on 22nd and 24th July, 2010 through August 9, 2010 for the Khyber-Pakhtunkhwa province, located in the north west of the country (red rectangle in Figure 1a). The diagram shows the probability distribution of precipitation based on the 51 ensemble members with the ensemble mean plotted as the black dotted line. The blue line represents the CMORPH observed rainfall. Good predictive skill of the July 28–29 event is found up to 6 days in advance. The same analysis done for various other monsoon pulses have resulted in similar conclusions (Figure 3).
Figure S7 shows an assessment of precipitation predictability in northern Pakistan using all available hindcast data. Predictability is shown as correlations between predicted and observed CMORPH rainfall values as a function of lead time for July based on 2007–2010 period. Note that for 2007, the model prediction extends only up to 10 days but up to 15 days for the 2008–2010 period. Correlations ≥0.7 were found for predictions 5 days in advance indicating useful predictive skill. Thus, the quantitative rainfall forecasts could be used as a robust variable in a flood forecasting scheme for Pakistan region.
In order evaluate whether the model can provide useful information with regards to the actual severity of the major rainfall events of July-early August 2010, all ECMWF forecasts made during the period were extracted and bias corrected. Then, the probability that the predicted rainfall would exceed the observed climatological average plus 1 standard deviation was computed. In other words, for each forecast, at each lead time, the percentage of ensemble members exceeding the threshold was computed. The exceedance threshold is calculated using 2003–2010 CMORPH data, with mean and standard deviation based on July-August daily average data. Results are shown in Figure 3 as shaded contours. The blue line represents the observed CMORPH rainfall averaged for the same region and the same time period. For example, the July 28 event was predicted almost 8 days in advance with a probability >60% over the climatological average plus 1 standard deviation (Figure 3). All the other events appear to have similar skill at the 8 to 10 day horizon.
From a climatological perspective, July and August precipitation rates were above average in Pakistan although not exceptionally so. However, in terms of rainfall rate, the monsoon pulses were extreme events compared to other years in the period 1998–2010. The devastating flooding occurred from a conspiracy of events. The summer of 2009 was a severe drought period with rainfall well below average (http://www.pakmet.com.pk/monsoon2009ver.pdf) so that vegetation may have been sparser during 2010. The region is mountainous with steep valleys and ridges. Furthermore, deforestation in northern Pakistan has been severe [e.g., Ali et al., 2006]. Deforestation and sparse undergrowth would exacerbate runoff through the steep valleys of the heavy rains that occurred during the month of July and early August.
The major result of the study is that the heavy rainfall pulses throughout July and early August were predictable with a high probability 6–8 days in advance. If these forecasts had been available to the regions of northern Pakistan, government institutions and water resource managers could have anticipated rapid filling of dams, releasing water ahead of the deluges. A high probability of flooding could have been anticipated.
Finally, it appears that Pakistan would benefit from a hydrological forecasting scheme similar to that developed for Bangladesh [Hopson and Webster, 2010; Webster et al., 2010]. The Bangladesh system incorporates the same form of statistically rendered ensemble precipitation forecasts as discussed above but coupled to a hybrid hydrological model. Working with Government of Bangladesh authorities, these 10-day river forecasts were communicated to the union (county) and village level allowing time to prepare for the floods for three major Brahmaputra floods during 2007/8 allowing the saving of household and agricultural effects and the successful evacuation of those in peril [Webster et al., 2010; Webster and Jian, 2010] (ADPC, Flood forecasts application for disaster preparedness: Post flood forecasts assessment 2008: Community response to CFAN forecasts, 2009, available at http://www.adpc.net/v2007/ and http://pacific.eas.gatech.edu/∼pjw/FLOODS).
This research has been supported by the Climate Dynamics Division of the National Sciences Foundation under Award NSF-ATM 0965610. Once again, we are indebted to ECMWF for providing data to make this analysis possible. We wish to thank J. A. Curry for interesting discussions as summarized at http://judithcurry.com/2010/09/20/pakistan-on-my-mind/and comments on the paper.
Noah Diffenbaugh thanks the two anonymous reviewers.
Format AvailableFull text: HTML | PDF
Copyright 2011 by the American Geophysical Union.
- Issue online:
- Version of record online:
- Manuscript Accepted:
- Manuscript Revised:
- Manuscript Received:
- flood forecast;
- 2010 Pakistan flooding
Auxiliary material for this article contains a text file and seven figures.
Auxiliary material files may require downloading to a local drive depending on platform, browser, configuration, and size. To open auxiliary materials in a browser, click on the label. To download, Right-click and select “Save Target As…” (PC) or CTRL-click and select “Download Link to Disk” (Mac).
Additional file information is provided in the readme.txt.
|grl27810-sup-0001-readme.txtplain text document, 3K||Readme.txt|
|grl27810-sup-0002-txts01.pdfPDF document, 588K||Text S1. Method description.|
|grl27810-sup-0003-fs01.epsPS document, 804K||Figure S1. Latitude-time longitude-time cross section for observed outgoing longwave radiation during June-July 2010.|
|grl27810-sup-0004-fs02.epsPS document, 2276K||Figure S2. Observed CMORPH precipitation for the 6 monsoon pulses during July-August 2010.|
|grl27810-sup-0005-fs03.epsPS document, 91K||Figure S3. Scatter diagram of monthly precipitation anomaly between CMORPH and TRMM over the period from 2003 to 2010.|
|grl27810-sup-0006-fs04.epsPS document, 217K||Figure S4. Seasonal mean precipitation for GPCP and CMORPH averaged over Pakistan and northern Pakistan.|
|grl27810-sup-0007-fs05.epsPS document, 187K||Figure S5. Number of heavy rainfall events over the summer in GPCP and CMORPH.|
|grl27810-sup-0008-fs06.epsPS document, 637K||Figure S6. Observed July-August cumulative CMORPH precipitation for years 2008 and 2010.|
|grl27810-sup-0009-fs07.epsPS document, 63K||Figure S7. Overall estimates of the predictability of precipitation in the Pakistan region versus lead-time for July based on 15-day forecasts from 2007–2010.|
Please note: Wiley-Blackwell is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.
- Adler, R. F., et al. (2003), The version 2 Global Precipitation Climatology Project (GPCP) monthly precipitation analysis (1979–present), J. Hydrometeorol., 4, 1147–1167
Unusually heavy monsoon rains in July and August 2010 left large swaths of Pakistan underwater. At least 18 million people were affected by the flood, and it is estimated that, more than six months later, several hundred thousand remain without even temporary shelter. As a result of lost crops and livelihoods from the flood and inadequate relief supplies, malnutrition continues to kill people. Like most floods, the Pakistani poor have suffered far more than those with resources to avoid the flood, or at least its aftermath.
Remains of a school destroyed by flooding, near Jacobabad by UK Department for International Development, on Flickr. Used under a Creative Commons license.
A paper in press in Geophysical Research Letters shows that the 2010 floods were extraordinary. Monsoonal rains tend to occur in pulses, with multi-day wet periods followed by multi-day dry periods, and while the total rainfall over Pakistan during the 2010 monsoon season was not unprecedented, the number and intensity of extremely heavy rains over northern Pakistan was very unusual. The authors are working with very limited historical and satellite data, but they estimate that the number of intense rain bursts that occurred in 2010 had a probability of less than 3% in any given year.
Using data from the European Centre for Medium Range Weather Forecasts collection of meteorological models, the authors of the new paper show that the timing and intensity of northern Pakistan’s monsoon rain bursts are predictable up to 6 to 8 days in advance – including the rains that caused the flooding in 2010.
Lead author, Peter Webster, and his coauthors from the Georgia Institute of Technology, draw the following conclusion from their analysis:
We conclude that if these extended quantitative precipitation forecasts had been available in Pakistan, the high risk of flooding could have been foreseen. If these rainfall forecasts had been coupled to a hydrological model then the high risk of extensive and prolonged flooding could have anticipated and actions taken to mitigate their impact.
The floods really kicked off with a burst of rain on 28-29 July 2010, and according to Webster’s reanalysis, that rainfall was predictable with good skill 7 days in advance (21 July). Webster and colleagues argue that if that forecast was available in Pakistan, lives would have been saved and the immensity of the disaster reduced. But, C. Christine Fair, writing on the Foreign Policy magazine website suggests that the flood was forecast in Pakistan.
In the middle of July, the PMD began tracking a storm brewing in the Bay of Bengal. This eastern weather system developed interactively with a western weather system to produce the massive rains and the subsequent super flood of 2010. On July 24, the PMD issued a flood warning to the provincial government of Khyber-Pakhtunkhwa (KPK). Despite these increasingly severe warnings, KPK’s citizenry did not believe them. … The PMD kept issuing warnings to KPK as the rains began to fall. However, as fate would have it, on July 28, … a passenger jet coming to Islamabad from Karachi crashed …With the media beset upon this tragic spectacle, the PMD’s warnings went unheeded as the rain began to fall.
So the Pakistani government did forecast the flood – at least four days out – in plenty of time to get people in northern Pakistan’s valleys out of the way. The problem was not with the meteorological and hydrologic science either internationally or in Pakistan. Instead, disaster was ensured when flood warnings were not taken sufficiently seriously by regional authorities, media, and residents.
Why wouldn’t flood warnings be heeded? Perhaps more could have been done to communicate to Pakistanis through channels whose authority they respected. Webster cites an example of flood warnings in Bangladesh being disseminated by imams at local mosques. The Foreign Policy article quoted above places some blame on media distractedness.
But there was also a more insidious reason the forecasted flood was ignored. It was a rare event, but it was also part of a new climatic pattern for Pakistan. As the Foreign Policy article describes it:
in recent years there has been a slow but steady change in the location where Pakistan’s major rainfalls concentrate. In the past, monsoon rains fell most intensely over the Punjab. Slowly and steadily, the concentration of rainfall has moved north and west to KPK. This redistribution of concentrated rainfall away from the Punjab and towards KPK explains why no one in KPK had any reason to believe the predicted weather.
Flooding frequency and intensity have increased in Pakistan in the last 30-40 years compared to earlier in the 20th century. Webster and coauthors state, “This recent increase is consistent with the increase in intensity of the global monsoon accompanying the last three decades of general global warming.” The flood warnings were ignored, in part, because the statistics of monsoon rain patterns are changing. Human memory and historical records are not good guidance if the weather system is changing. In situations like this one, the past is not the key to the present.
There are lots of things that should have been improved to lessen the magnitude of the Pakistani flood disaster – reservoir management should have been altered; emergency relief supplies should have been distributed more equitably, broadly, and consistently; international assistance should have been much more generous – but the two big lessons for hazard mitigation coming out of the Pakistan floods seem to be: “find a system for making sure that warnings are issued and that they actually make it to people in harm’s way” and “don’t assume the climate of living memory is a very good indicator of the weather of the present and future.”
Webster, P. J., Toma, V.E., & Kim, H.-M. (2011). Were the 2010 Pakistan floods predictable? Geophysical Research Letters : 10.1029/2010GL046346
Categories: by Anne, climate science, geohazards, hydrology, paper reviews, society
Tags: climate change, disasters, floods, hazard mitigation, hydrology, Pakistan