GEOPHYSICAL RESEARCH LETTERS, VOL.???, XXXX, DOI:10.1029/, Trends in 20th Century Drought over the Continental United States Konstantinos M. Andreadis Civil and Environmental Engineering, University of Washington, Seattle, WA Dennis P. Lettenmaier Civil and Environmental Engineering, University of Washington, Seattle, WA K.M. Andreadis, Department of Civil and Environmental Engineering, University of Washington, Wilson Ceramic Lab, Box 352700, Seattle, WA 98195-2700. (kostas@hydro.washington.edu)
X - 2 ANDREADIS AND LETTENMAIER: U.S. DROUGHT TRENDS We used a simulated dataset of hydro-climatological variables to examine for 20th century trends in soil moisture, runoff, and drought characteristics over the conterminous United States (U.S.). A predominant wetting trend is apparent in both model soil moisture and runoff over much of the U.S., with a few drying trends in parts of the Southwest. The trend patterns were qualitatively similar to those found in streamflow records observed at a station network minimally affected by anthropogenic activities. This wetting trend is consistent with the general increase in precipitation in the latter half of the 20th century. Droughts have, for the most part, become shorter, less frequent, and cover a smaller portion of the country over the last century. The main exception is the Southwest and parts of the interior of the West, where, notwithstanding increased precipitation (and in some cases increased soil moisture and runoff), increased temperature has led to trends in drought characteristics that are mostly opposite to those for the rest of the country especially in the case of drought duration and severity, which have increased.
1. Introduction ANDREADIS AND LETTENMAIER: U.S. DROUGHT TRENDS X - 3 Droughts are one of the most costly natural disasters, with estimated annual U.S. losses between $6-8B [FEMA, 1995]. According to NOAA s National Climate Data Center (NCDC), the recent western U.S. drought cost over $10B for 2002 alone [NCDC, 2003]. Potential changes in the characteristics of drought, such as severity and frequency, would have adverse effects on water management and aquatic ecosystems. Changes in the severity of drought as climate changes associated with mid-continental drying predicted by some climate models over the next century [Wetherald and Manabe, 1995] are therefore of great concern. Furthermore, although droughts are generally associated with decreased precipitation (and on average, global precipitation is projected to increase in the future by most climate models); increased rainfall does not necessarily mean less intense droughts. For example, less frequent but heavier precipitation events could lead to longer dry periods, even though rainfall amounts might have increased. There is a general consensus that temperature and precipitation increased in the latter half of the 20th century over the continental U.S. [IPCC, 2001]. Nonetheless, the nature of evaporation changes is unclear (due largely to the paucity of records that measure actual, as opposed to potential, evaporation), and so the nature of changes in the characteristics of drought is not immediately apparent from the observed precipitation and temperature (or streamflow, for that matter) changes. Furthermore, long-term records of a key variable related to drought, soil moisture, are essentially nonexistent. For this reason, the need for an objective method for investigating changes in the frequency and severity of droughts becomes apparent [Trenberth et al., 2004].
X - 4 ANDREADIS AND LETTENMAIER: U.S. DROUGHT TRENDS We suggest an approach to the problem that is based on reconstructions, using a wellvalidated land surface hydrology model, of key hydro-climatic variables that have not been observed directly over the long term. Andreadis et al. [2005] successfully used a long- term simulated dataset to reconstruct 20th century drought conditions over the continental U.S., and compared the extent, severity and durations of the largest events for that period. While that study did not evaluate changes in drought characteristics, the underlying methods nonetheless form the basis for the analyses we present here. The goal of this paper, therefore, is to determine whether changes in drought characteristics have occurred over the continental U.S. in the past 80 years, and if so, to explore the nature of the changes that have occurred. 2. Hydrologic Model and Datasets We used the Variable Infiltration Capacity (VIC) model [Liang et al., 1994] to simulate soil moisture and runoff over the continental U.S., at a spatial resolution of one-half degree and a daily timestep. VIC essentially solves an energy and water balance over each model grid cell, and accounts for the effects of sub-grid scale variability in soil, vegetation, precipitation, and topography by partitioning each grid cell into different tiles. The soil column is represented as three layers, while vegetation and topography data as well as calibration parameters were identical to those used in [Maurer et al., 2002]. The forcing data, which include precipitation, air temperature, and wind speed, were obtained from the NCDC using the methods outlined in [Maurer et al., 2002], and extended from 1915 to 2003 (also used as the simulation period). In addition, the data were corrected for temporal heterogeneities using the procedure described in [Hamlet and Lettenmaier,
ANDREADIS AND LETTENMAIER: U.S. DROUGHT TRENDS X - 5 2005] that essentially adjusts the decadal variability to match that of the high-quality U.S. Historical Climatology Network (HCN). In previous studies (e.g. [Maurer et al., 2002]) the VIC model has been applied at higher resolutions (typically one-eighth degree), however droughts tend to be spatially persistent over large areas, and for this purpose the coarser one-half degree resolution is adequate. VIC has been successfully used to simulate streamflow over a number of large river basins in North America and elsewhere, and at continental to global scales as well [Nijssen et al., 1997; Lohmann et al., 1998]. In previous applications to the continental U.S., Maurer et al. [2002] and others have demonstrated good agreement between simulated soil moisture and point measurements, between observed and simulated snow water equivalent, and between observed and simulated streamflow. 3. Drought Characteristics A generally accepted definition of drought is elusive, although many definitions appear in the literature. Here, we use the threshold approach described in [Dracup et al., 1980], where a drought is said to have occurred when the variable of interest (e.g. streamflow) is below a pre-defined threshold. In general four different types of drought can be considered: meteorological, agricultural, hydrological, and socio-economic. In this study we examine agricultural and hydrological droughts, using model-derived soil moisture and runoff as drought indicators. Model soil moisture and runoff were aggregated to monthly values, with the first 10 years of simulated data truncated to minimize model initialization effects. Soil moisture and runoff values were expressed as percentiles relative to all simulated values for a given
X - 6 ANDREADIS AND LETTENMAIER: U.S. DROUGHT TRENDS grid cell and month. This facilitates droughts across the domain, because anomalies in absolute terms reflect different severities in different parts of the domain due to climatic, vegetation, soils, and other differences. Furthermore, percentiles are convenient in that they provide an ordinal range from zero to one, over which the variables vary uniformly by construct, and seasonal variations are removed as well. 4. Trend Analysis One of the most widely used statistical tests for trend in hydro-climatic time series, is the Mann-Kendall (MK) test [Mann, 1945; Helsel and Hirsch, 1992]. The test is nonparametric, so it does not make any assumptions about normality or linearity. Hirsch et al. [1982] proposed a modified version of the MK test more appropriate for seasonal data, termed the Seasonal MK test. An advantage of the seasonal MK test for hydrologic data is that trends in seasons with small values are not dominated by larger values in other seasons [Lettenmaier et al., 1994]. The seasonal MK test also accounts for serial dependence, e.g. correlation between seasons (months), by inflating the variance of the test statistic [Hirsch and Slack, 1984]. While the seasonal test does not account for serial correlation over years (e.g. correlation from one January to the next), this issue can be addressed using an effective independent sample size as proposed by Hamed and Rao [Hamed and Rao, 1998]. 5. Results 5.1. Trends in soil moisture and runoff We applied the seasonal MK test to the monthly simulated soil moisture and runoff data, and the standard MK test to the series of individual months. A value of 0.05 was
ANDREADIS AND LETTENMAIER: U.S. DROUGHT TRENDS X - 7 chosen as the local significance level for a 2-sided test. The statistically significant trends for model soil moisture can be seen in Figure 1, which shows that 1450 cells (43.6% of the domain) exhibit an upward trend, while many fewer trended downward (95 cells, or 2.9%). Figure 1 shows that the wetting trends cover the majority of the country, with most of the drying trends being in the Southwest. Annual trends in runoff are very similar to those found for soil moisture. In particular, upward trends were found for 41.5% of the U.S. (1377 cells), while the number for downward trends was 2.3% (77 cells). These results agree with previous studies that have suggested a general increase in streamflow over the conterminous U.S. [Lins and Slack, 1999; McCabe and Wolock, 2002; Groisman et al., 2004]. In order to provide a rough assessment of the inferred trends, we tested stations in the Hydro-Climatic Data Network (HCDN) for trends. This dataset contains streamflow daily data from a set of 256 stations that have minimal upstream regulation making them suitable for identifying long-term trends. The stations used here are a subset of the stations utilized by [Lettenmaier et al., 1994] for which we were able to obtain data from 1925-2003 so as to be consistent with our model results. Results from the application of the seasonal MK test to the HCDN dataset are qualitatively similar to the model results (Figure 2), with 72 stations (28.1%) having an upward trend and only 6 stations (2.3%) having downward trends. It is important to note that we would not expect an exact correspondence because, among other reasons, the HCDN stations observe streamflow, which is effectively an integrator of runoff over an area, whereas the model produces grid cell runoff. In addition, spatial coverage of the HCDN stations that provided a long-term
X - 8 ANDREADIS AND LETTENMAIER: U.S. DROUGHT TRENDS record is much smaller. Nonetheless, the general consistency in the number of trends is encouraging. When examining the spatial patterns of statistically significant trends for each month (not shown), a set of persistent clusters of wetting trends exist over the Great Lakes region, and over parts of the northeast U.S. which is most evident during the winter months. Another upward trend cluster appears in the Lower Mississippi and Gulf regions during October-December, however during the summer months there is a group of drying trends in the Southeast. In addition, a few drying trends occur over the Colorado River basin and California during the summer months. One of the features of the annual trends map is the spatial clustering of uptrends and downtrends. This is somewhat expected, because climatological data are highly correlated spatially. The question arises as to how this spatial correlation affects the results of the trend analysis. This issue was first addressed in the context of climatic trends by Livezey and Chen [1983] who proposed a construct for local (e.g., grid cell) and field (e.g., spatial domain) significances. They proposed to resolve the issue of spatial correlation via a Monte Carlo procedure that replicates the spatial structure of the input fields, and performs the statistical test on the generated fields (where the null hypothesis is true by construct). From the number of rejections of the null hypothesis one can build an empirical distribution and determine the threshold fraction for field significance. We applied this procedure and evaluated field significance by resampling the model-derived data by creating 500 permutations of the years of data and shuffling the order of the time series, thus retaining the spatial correlations. Using a level of 0.05 for field significance,
ANDREADIS AND LETTENMAIER: U.S. DROUGHT TRENDS X - 9 the fraction thresholds were exceeded (20.4% for soil moisture, and 14.1% for runoff) indicating that the trends observed are unlikely to have occurred due to chance given the observed spatial correlation. 5.2. Trends in drought duration, frequency, severity, and spatial extent Time series of drought durations were constructed for each model grid cell by creating sequences of the number of consecutive months in which the soil moisture/runoff percentiles were below specified thresholds (from 10 to 50%). Figure 3 shows a map of the locations where significant (0.05 level) trends were found using the MK test for a drought threshold P = 20%. The number of significant trends is not large (47 upward and 102 downward for runoff, and 62 upward and 161 downward for soil moisture). There are clear spatial patterns, however. Figure 3 shows that the Great Lakes, Mississippi basin, and Pacific Northwest have decreasing drought durations, while parts of the west including California and the interior have increased drought durations. A noteworthy feature is that although the Pacific Northwest had a predominantly wetting trend, extreme droughts (P = 10%) have become longer. Drought severities were computed for each grid cell as the cumulative departure from the 20th percentile threshold, averaged sequentially for the duration of each event identified. Figure 4 shows the location of statistically significant trends in drought severity for model runoff. Decreasing trends (198 cells) are mostly found in the Northeast, Great Lakes, Lower Mississippi, and the Pacific Northwest. Upward trends are fewer (121 cells) and are located mostly in Texas, the Southwest and intermountain West. The results were similar for model soil moisture, with 73 grid cells having upward and 180 cells having
X - 10 ANDREADIS AND LETTENMAIER: U.S. DROUGHT TRENDS downward trends. Field significance cannot be evaluated properly in the case of drought characteristics, since the time series do not have the same length for different model grid cells. The MK test was also applied to time series of drought frequency for different thresholds (not shown). There is a predominant reduction in drought frequency for the eastern U.S. and Midwest, with a few significant upward trends in parts of the Southwest. As the drought threshold was reduced so did the area covered by the downward trends. In contrast, the number of upward trends in the Southwest was almost invariant with respect to the drought threshold. Finally we examined trends in the spatial extent of drought, that is the percent area of the U.S. that was experiencing drought (using different thresholds) at each month. Soil moisture (and runoff) drought spatial extent showed a downward trend which however was insignificant for all thresholds (10-50%). 6. Discussion An obvious question arises as to the causes of the observed trends. Soil moisture results from a complicated balance of precipitation (the partitioning of which into infiltration and runoff is a nonlinear function of soil moisture), and evaporation, which depends on surface atmospheric conditions as well as soil moisture. Over much of the U.S., there has been a general upward trend in precipitation, but temperature has increased as well. Aside from the Southwest and parts of the interior west, model soil moisture has been increasing, as have most soil moisture-related drought indicators. While no long observation records of soil moisture exist, our results for runoff are similar to those for soil moisture, and
ANDREADIS AND LETTENMAIER: U.S. DROUGHT TRENDS X - 11 our runoff trends are comparable to those in observed streamflow records (most trends are upward, aside from parts of the Southwest and West). The area covered by decreasing soil moisture and runoff (and increases in most drought indicators) is interesting, as precipitation has increased in this area as well. While a detailed analysis is beyond the scope of this paper, we can examine the spatial patterns of trends in precipitation and temperature and compare them qualitatively with the soil moisture and runoff trends. Over the period of our analysis, precipitation has either increased or remained (statistically) unchanged at the vast majority of locations (499 uptrends and only 19 downward trends). In addition, the spatial patterns of the precipitation trends (not shown) agree fairly well with those for runoff. On the other hand, the effects of temperature trends on the spatial patterns of soil moisture/runoff are not so clear. Over the period of analysis, temperature increased for most of the western U.S. and the Northeast (750 cells), while there were significant decreases for the Gulf region and part of the Southeast (345 cells). It appears that in those locations of the Southwest and West where downward trends in soil moisture were observed, the effect of increased evaporative demand associated with the higher temperatures more than canceled the effects of increased precipitation. 7. Conclusions We used a constructed time series of soil moisture and runoff over the continental U.S. to examine trends in soil moisture and runoff, and drought characteristics related to these variables for the period 1925-2003. Over much of the country, there has been a wetting trend which is reflected in a predominance of upward trends in both model-derived soil moisture and runoff. These trends are generally consistent with increases in precipitation
X - 12 ANDREADIS AND LETTENMAIER: U.S. DROUGHT TRENDS during the latter half of the 20th century observed over most of the U.S. [Groisman et al., 2004]. Furthermore, trends in the simulated runoff are similar to those in observed records of streamflow at a set of index stations that have been minimally affected by anthropogenic activities. Trends in most drought characteristics (duration, frequency, severity, and extent) are similar to those in soil moisture and runoff, i.e., droughts have, for the most part, become shorter, less frequent, less severe, and cover a smaller portion of the country over the last century. The main exception is the Southwest and parts of the interior of the West, where, notwithstanding increased precipitation (and in some cases increased soil moisture and runoff), increased temperature has led to trends in drought characteristics that are mostly opposite to those for the rest of the country especially in the case of drought duration, severity and frequency, which have increased. References Andreadis, K. M., E. A. Clark, A. W. Wood, A. F. Hamlet, and D. P. Lettenmaier (2005), 20th century drought in the conterminous united states, J. Hydrometeor., 6, 885-1001. Dracup, J. A., K. S. Lee, and E. G. Paulson (1980), On the definition of droughts, Water Resour. Res., 16, 297 302. FEMA (1995), National Mitigation Strategy: Partnerships for Building Safer Communities, 40 pp., Federal Emergency Management Agency, Washington D.C. Groisman, P. Y., R. W. Knight, T. R. Karl, D. E. Easterling, B. Sun, and J. H. Lawrimore (2004), Contemporary changes of the hydrological cycle over the contiguous united states: trends derived from in situ observations, J. Hydrometeor., 5, 64 85.
ANDREADIS AND LETTENMAIER: U.S. DROUGHT TRENDS X - 13 Hamed, K. H., and A. R. Rao (1998), A modified mann-kendall trend test for autocorrelated data, J. Hydrology, 204, 182 196. Hamlet, A. F., and D. P. Lettenmaier (2005), Production of temporally consistent gridded precipitation and temperature fields for the continental united states, J. Hydrometeor., 6, 330 336. Helsel, D. R., and R. M. Hirsch (1992), Statistical methods in water resources, 522 pp., Elsevier, New York. Hirsch, R. M., and J. R. Slack (1984), A nonparametric trend test for seasonal data with serial dependence, Water Resour. Res., 20, 727 732. Hirsch, R. M., J. R. Slack, and R. A. Smith (1982), Techniques of trend analysis for monthly water quality data, Water Resour. Res., 18, 107 121. IPCC (2001), Climate Change 2001: The Scientific Basis, Contribution of Working Group I to the Third Assessment Report of the Intergovermental Panel on Climate Change, 881 pp., Cambridge University Press. Lettenmaier, D. P., E. F. Wood, and J. R. Wallis (1994), Hydro-climatological trends in the continental united states, 1948-88, J. Climate, 7, 586 607. Liang, X., D. P. Lettenmaier, E. F. Wood, and S. J. Burges (1994), A simple hydrologically based model of land surface water and energy fluxes for gcms, J. Geophys. Res., 99, 14,415 14,428. Lins, H. F., and J. R. Slack (1999), Streamflow trends in the united states, Geophys. Res. Let., 26, 227 230. Livezey, R. E., and W. Y. Chen (1983), Statistical field significance and its determination by monte carlo techniques, Mon. Wea. Rev., 111, 46 59.
X - 14 ANDREADIS AND LETTENMAIER: U.S. DROUGHT TRENDS Lohmann, D., E. Raschke, B. Nijssen, and D. P. Lettenmaier (1998), Regional scale hydrology. ii. application of the vic-2l model to the weser river, germany, Hydrol. Sci. J., 43, 143 158. Mann, H. B. (1945), Nonparametric tests against trend, Econometrica, 13, 245 259. Maurer, E. P., A. W. Wood, J. C. Adam, and D. P. Lettenmaier (2002), A long-term hydrologically based dataset of land surface fluxes and states for the conterminous united states, J. Climate, 15, 3237 3251. McCabe, G. J., and D. M. Wolock (2002), A step increase in streamflow in the conterminous united states, Geophys. Res. Let., 29, 2185. NCDC (2003), Climate of 2002 - Annual Review U.S. Drought, National Climatic Data Center, http://www.ncdc.noaa.gov/oa/climate/research/2002/ann/drought-summary.html. Nijssen, B., D. P. Lettenmaier, X. Liang, S. W. Wetzel, and E. F. Wood (1997), Streamflow simulation for continental-scale river basins, Water Resour. Res., 33, 711 724. Trenberth, K. E., J. T. Overpeck, and S. Solomon (2004), Exploring drought and its implications for the future, Eos, Trans. Amer. Geophys. Union, 85, 27. Wetherald, R. T., and S. Manabe (1995), The mechanisms of summer dryness induced by greenhouse warming, J. Climate, 8, 3096 3108.
ANDREADIS AND LETTENMAIER: U.S. DROUGHT TRENDS X - 15 Figure 1. Annual trends in model soil moisture. Blue triangles show upward trends, while downward trends are shown as red inverted triangles.
X - 16 ANDREADIS AND LETTENMAIER: U.S. DROUGHT TRENDS Figure 2. Annual trends in streamflow observed from the HCD Network. Upward (blue triangles) and downward (red inverted triangles) trends are shown in the larger map. The smaller map shows the spatial coverage of the HCDN stations.
ANDREADIS AND LETTENMAIER: U.S. DROUGHT TRENDS X - 17 Figure 3. Trends in drought duration. Blue inverted triangles show downward trends, and red triangles show upward trends.
X - 18 ANDREADIS AND LETTENMAIER: U.S. DROUGHT TRENDS Figure 4. Trends in drought severity. Upward trends are shown as red triangles, while downward trend as blue inverted triangles.