TMCnet News
An Alternative Index for the Contribution of Precipitation on Very Wet Days to the Total Precipitation [Journal of Climate](Journal of Climate Via Acquire Media NewsEdge) ABSTRACT Daily precipitation series from more than 1800 stations across Europe are analyzed for the fraction of the total precipitation due to very wet days: that is, days with precipitation amounts exceeding the 95th percentile. This fraction is calculated on a seasonal (3-monthly) basis for the period 1961-2010.Anew index S95pTOT is introduced as an alternative to the frequently used index R95pTOT. Contrary to R95pTOT, which uses a fixed climatological 95th percentile, the new index assumes a separate 95th percentile for each year. Based on a Weibull distribution fit to the wet-day precipitation amounts, an analytical expression for S95pTOT is derived. It is shown that R95pTOT is strongly influenced by changes in the mean wet-day precipitation, whereas S95pTOT is more representative of changes in the distributional shape. The results for S95pTOT do not support the conclusion for a disproportional increase of extreme precipitation over northern Europe as was concluded from the trend in R95pTOT in earlier studies. Also, the contrast between trends in northern and southern Europe in winter is less pronounced for S95pTOT than for R95pTOT. (ProQuest: ... denotes formulae omitted.) 1. Introduction Now that consensus builds up on how the rise of the global temperature affects the intensity of the hydrological cycle (see, e.g., Trenberth 2011), also the interest in the changes in seasonal mean precipitation and heavy precipitation is growing.Most information on the changes in precipitation refers to the seasonal mean precipitation but, when it comes to societal impacts, the changes in heavy precipitation are more relevant (Field et al. 2012). The question that arises is how the changes in heavy precipitation relate to those in the mean precipitation. Changes of heavy precipitation relative to the mean have been studied from observations as well as climate model simulations for different regions, using a variety of methods. To separate changes in extreme precipitation from those in the mean, many studies consider the portion of the annual or seasonal precipitation total contributed by days with precipitation amounts exceeding some high level. The level is usually adapted to the probability distribution of daily precipitation amounts for the location of interest. This yields a standardized measure, suitable to compare locations with different precipitation regimes. The index R95pTOT, recommended by the Expert Team on Climate Change Indices (ETCCDI; Nicholls and Murray 1999; Zhang et al. 2011), is an example of such a measure. This index is used in many regional (Peterson and Manton 2008) and global (Frich et al. 2002; Alexander et al. 2006) studies. In the definition of R95pTOT, the 95th percentile of wet-day precipitation amounts, calculated for a chosen reference period, is used to identify very wet days. The sum of daily precipitation amounts exceeding this level is divided by the total precipitation to obtain the fraction of precipitation on very wet days. Trends in R95pTOT have been analyzed by KleinTank andKeuroonnen (2003) for the period 1949-99 for 151 stations across Europe on an annual basis. They found a significant positive trend in R95pTOT for a substantial part of those stations for which also a significant positive trend in the annual totals was found. They concluded that this supports the notion of an amplified response of the extremes relative to the precipitation total. In more recent literature R95pTOT is also considered (Turco and Llasat 2011; Ducic et al. 2012; Sillmann et al. 2013; Donat et al. 2013) as an indicator of the amplified response of extreme precipitation events to climate change. A disadvantage of the index R95pTOT at a seasonal scale, reported by Zolina et al. (2009), is that the number of days falling into the highest category can vary strongly over time or become zero in seasons for which the highest precipitation amount does not exceed the 95th percentile of the reference period. They introduced an alternative index, resembling R95pTOT in its definition, but based on the gamma distribution for the wet-day amounts and the associated theoretical distribution of the fractional contribution of the wettest days to the seasonal or annual total. This new index shows less yearto- year variation, although it is essentially similar to the procedure of Klein Tank and Keuroonnen (2003) for R95pTOT. In this paper, we address an additional disadvantage of R95pTOT as an indicator of an amplified response of extreme precipitation: namely, that a change in the mean without a change in the shape of the distribution also affects R95pTOT. Hence, a trend in R95pTOT not necessarily represents a change in the distributional shape associated with an amplified response of extreme precipitation (section 3). We introduce an index S95pTOT, which does not have this disadvantage and also shows less variation over the years than R95pTOT. The fundamental difference between R95pTOT and S95pTOT is that in the latter no fixed level is used to define the very wet days. Instead, this level is allowed to vary by year and season. In addition, an analytical expression is derived for S95pTOT in terms of the shape parameter of a fitted Weibull distribution. The 95th percentile is no longer explicitly needed but is instead an implicit part of this expression. Because the analytical expression for S95pTOT is derived in terms of the Weibull shape parameter, the association between S95pTOT and the distributional shape follows naturally. The trends in the new index are compared with those of the traditional index for more than 1800 stations across Europe. The central question is whether the trends of S95pTOT also support the earlier findings of a disproportional change of the extremes over Europe. The remainder of the paper is structured as follows: Section 2 describes the data used in this study. In section 3 we argue why R95pTOT is suboptimal as an indicator of a disproportional change of the extremes relative to the mean and define S95pTOT. Results of tests for trend in the index time series are presented in section 4. The final section presents the conclusions and some final remarks concerning the new index. 2. Precipitation data The daily precipitation data used in this study was obtained from the European Climate Assessment and Dataset (ECA&D; Klein Tank et al. 2002). The so-called blended station series in ECA&D covering the study period 1961-2010 were initially selected. All these series were tested for homogeneity, using the methods of Wijngaard et al. (2003), and labeled as ''useful'' for the somewhat longer test period 1951-2011. Homogeneity tests for this period are part of the standard quality control of ECA&D.We do not expect the homogeneity of the series for the study period to be substantially different. From the selected series those series were rejected for which less than 80% of the data were available over the study period. The resulting set of stations varies in size through the year between 1815 stations in June-August (JJA) to 1856 stations in March-May (MAM). Data from ECA&D were also used by Zolina et al. (2009). That study was restricted to the original, nonblended series to rule out potential artifacts introduced by the blending process applied within ECA&D. In this process, time series are filled up with data from neighboring stations, where possible, and augmented with synoptical data, distributed through the Global Telecommunication System up to the current month. However, the quality control procedures and the homogeneity checks in ECA&D, to which only the ''blended'' station series are subjected, substantially reduces the risks of potential artifacts. 3. Method a. Background of the used indices The index R95pTOT is defined for a season within a specific year as the ratio of the sum of all daily precipitation amounts exceeding the climatological (1961- 90) 95th percentile Q of the wet-day precipitation amounts to the total precipitation (see also the full definition by Klein Tank and Keuroonnen 2003). Wet days are here defined as those days with d51mmof precipitation or more. Let g be the probability density of all daily precipitation amounts (including those on dry days), m be the corresponding mean, and N be the number of days in the season of interest. The sum of the precipitation amounts exceeding Q can then be approximated by ... and the total precipitation can be approximated by Nm, which leads to the following approximation of R95pTOT: ... (1) This can be rewritten in terms of the probability density gw and mean mw of the wet-day precipitation amounts as ... (2) (see appendix A). Because Q is kept constant, a change in R95pTOT may be due to a change in the mean rather than a change in the shape of the distribution. This is best illustrated by introducing the standardized precipitation amount x0 5x/mw as integration variable in Eq. (2), ... (3) where g0 w(x0)5mwgw(x) is the density of x0. Now assume that g0 w does not change over time, whereas mw does. The latter induces a change in R95pTOT (to a first order) of ... (4) The change in R95pTOT induced by the change in mw is thus to a first order proportional to the relative change in mw. Though Q refers to the climatological 95th percentile over some reference period, it may not correspond with the 95th percentile in an individual season, and therefore R95pTOT does not represent the contribution of the upper 5%. Therefore, we introduce a modified index RS95pTOT based on the time-varying 95th percentile q for the season of interest. Analogous to Eqs. (2) and (3), this index can be written as ... (5) In contrast to R95pTOT, a change in the mean does not alter the integral as long as g0 w does not change, because then q/mw remains unchanged. We now assume that the daily precipitation amounts X on wet days follow a two-parameter Weibull distribution with scale parameter a and shape parameter c, shifted by the wet-day threshold d. The corresponding distribution function is given by ... (6) and the probability density is given by ... (7) The expression for RS95pTOT then becomes (appendix A) ... (8) where P(, ) is the normalized lower incomplete gamma function [Abramowitz and Stegun 1965, Eq. (6.5.1)] and G() is the complete gamma function [Abramowitz and Stegun 1965, Eq. (6.1.1)]. Equation (8) can be simplified considerably by taking the limit as d / 0, yielding the new index proposed in this paper, ... (9) S95pTOT can be seen as an approximation of RS95pTOT for small d and has the advantage of depending only on the Weibull shape parameter. The factor in Eq. (8) preceding the square brackets is the ratio between the mean excess over the wet-day threshold and the mean wet-day precipitation, which is very close to one for d sufficiently small. The first term within the square brackets equals 0.05 times the wet-day threshold divided by the mean excess, which is very small. Assuming a shifted gamma distribution for the wet-day precipitation amounts, an expression similar to Eq. (9) can be derived, depending on the gamma shape parameter only. The preference for the Weibull distribution is motivated in appendix B, by means of an L-moment ratio diagram. Alternatively, RS95pTOT and S95pTOT can be estimated empirically, without assuming a distribution. The index RS95pTOT can be calculated as the sum of the nq largest wet-day precipitation amounts divided by the sum over all n wet-day amounts, where nq equals 0.05n, rounded to the nearest integer. The empirical estimate of S95pTOT can be obtained similarly, except that in this case all wet-day precipitation amounts should be reduced by d beforehand. b. Application to station data For each station (3-monthly) seasonal time series of the original index R95pTOT and the new index S95pTOT for the period 1961-2010 were analyzed. The series of RS95pTOT were also studied to assess the effect of neglecting d. The series of R95pTOT were obtained from ECA&D to ensure correspondence with earlier studies based on this index. The series of RS95pTOT and S95pTOT were calculated using the empirical as well as the parametric approach. For the parametric estimates, a two-parameter Weibull distribution was fitted to the shifted wet-day amounts (over d) using probability weighted moments (Boes et al. 1989; Hosking and Wallis 1997). Only seasons with 80 or more days with data and 10 or more wet days (d 5 1 mm) were taken into account, to prevent the empirical estimate from becoming zero. The estimation of S95pTOT was tested with synthetic data from a knownWeibull distribution, as described in appendix C. The estimates of S95pTOT were found to be nearly unbiased. Furthermore, estimating S95pTOT parametrically rather than empirically reduced the uncertainty. To detect trends in the time series of R95pTOT, RS95pTOT, and S95pTOT, theMann-Kendall statistic (MK; Hamed and Rao 1998; Yue et al. 2002a) was determined for each station and each season. In the theoretical variance of MK, the possibility of ties is taken into account. The potential influence of a nonzero lag-1 autocorrelation is dealt with by trend-free prewhitening as discussed by Yue et al. (2002b). Apart from testing the statistical significance of the MK statistic for each station individually, the significance of the percentage r1 of stations with a positive value of the MK statistic was tested. This field significance was determined using a block permutation procedure with a block size of 3 yr to preserve the dependence between successive years. The seasonal values of R95pTOT, RS95pTOT, and S95pTOT were permutated for all stations simultaneously to preserve the spatial dependence. The latter can have a considerable influence on the variance of r1 (cf. Douglas et al. 2000). This permutation procedure is similar to the moving block bootstrap in Kiktev et al. (2003), except that the blocks are sampled without replacement, so every sample consists of a random permutation of blocks. In contrast with the bootstrap, the permutation procedure preserves the sample size for stations with missing data. For each permutation replication, the percentage r1 * of stations was determined for whichMK.0, and this number was compared with the value r1 for the observations. Let f be the fraction of permutation replications for which r1 * .r1. A twosided significance probability was obtained as p 5 2min (f, 1 2 f ). 4. Results Figure 1 shows the climatological mean (1961-2010) of R95pTOT (top) and S95pTOT determined either empirically (middle) or based on the Weibull distribution (bottom) in the December-February (DJF) season. It is seen that the mean of S95pTOT is somewhat higher than that of R95pTOT and that this difference is generally larger than the difference between both versions of S95pTOT. The maps give no indication of a bias resulting from the Weibull assumption. Figure 2 shows the relative standard deviation of R95pTOT and both versions of S95pTOT, expressed as a percentage of the mean. The standard deviation was determined here using the mean square of successive differences in order to rule out the influence of trends. The standard deviation of R95pTOT is larger than that of S95pTOT. This reflects the strong interannual variation of R95pTOT, partly related to variations in themean wet-day precipitation [see Eq. (4)]. The values of S95pTOT based on the Weibull distribution are slightly less variable than those determined empirically. Figure 3 shows the sign and the significance of the trend in R95pTOT (top) and S95pTOT (bottom) for the DJF season. Trends in RS95pTOT (not shown here) reveal a pattern very similar to that of S95pTOT. While particularly the northern part of Europe is dominated by a (significant) increase of R95pTOT, this is far less pronounced for S95pTOT. At many sites the signs of the trend in S95pTOT and R95pTOT differ. Insignificant (positive or negative) trends become significantly negative and significant positive trends are reduced to insignificant trends. This is particularly striking in southern Scandinavia, the Netherlands, Germany, and the United Kingdom. In Spain and southern France, the number of stations with a negative trend is less for S95pTOT than for R95pTOT. As a whole, the contrast between the trends in northern and southern Europe is more diffuse for S95pTOT than for R95pTOT. Table 1 displays the percentages of the stations for which a positive trend was found in the indices S95pTOT, RS95pTOT, and R95pTOT; the fraction of wet days fw; the mean wet-day amount mw; and the coefficient of variation (CV) of wet-day precipitation CVw. A distinction is made between the stations north and south of the 48th parallel. For each percentage the field significance (based on the block permutation procedure described in section 3b, with 1000 permutation replications) is also listed. In DJF and MAM the percentages of positive trends in R95pTOT clearly indicate an increase for the majority of the northern stations. This majority is field significant. The same holds for the corresponding percentages for mw. The percentage of positive trends in S95pTOT and RS95pTOT for northern Europe in the DJF and MAM seasons is substantially lower than for R95pTOT, in many cases even below 50%. For the southern stations, the percentage of positive trends in R95pTOT is smaller than for the northern stations in those seasons and not field significant. Consistently, no field significant trend is found for mw. These results confirm the suggested influence of changes in mw on R95pTOT. The percentages of positive trends in\ R95pTOT for the southern stations are of the same order as those in S95pTOT and RS95pTOT in DJF andMAM. Though the trends in R95pTOT are not field significant for these stations, the trends in S95pTOTand RS95pTOT are in a number of cases. Finally, it is observed that in general the percentage of positive trends in S95pTOT andRS95pTOT corresponds betterwith that inCVw than that in mw. A large contrast is observed between northern and southern Europe regarding the changes in the fraction of wet days in the seasons DJF and JJA. The majority of the northern stations show an increase, whereas for the majority of the southern stations a decrease is found, in particular in DJF. The latter is field significant. In MAMan overall decrease is seen whereas in September- November (SON)most stations in northern and southern Europe show an increase. Differences between the results for S95pTOT and RS95pTOT in Table 1 reflect the influence of the finite wet-day threshold d. These differences are generally small compared to the differences with the results for R95pTOT. Also (not shown here) the spatial pattern of the MK statistics for RS95pTOT is very similar to that for S95pTOT, shown in the bottom panel of Fig. 3 for the DJF season. In particular for the northern part of Europe it is likely that the results of Table 1 are predominated by the trends in regions with a high station density in Fig. 3. It is interesting to zoom in on these regions. We defined two subregions, NorSwe between 558 and 658N and between 58 and 27.58E, containing southern Norway, Sweden, and a small part of Finland, and NetGer between 478 and 558N and between 38 and 158E, covering primarily the Benelux and Germany. Figure 4 shows the same as Fig. 3 but zooms in on the subregions NorSwe (Figs. 4a,c) and NetGer (Figs. 4b,d). For many stations within both regions, the trends in S95pTOT and R95pTOT have opposite signs (increasing R95pTOT in Figs. 4a,b versus decreasing S95pTOT in Figs. 4c,d), though fewer significant trends are counted for S95pTOT than for R95pTOT. For NorSwe, the most striking are the Norwegian stations near the coast, some of which have significant opposite trends for R95pTOT and S95pTOT. In subregion NetGer, the trends in the Netherlands stand out. The trend in R95pTOT is in general positive throughout the country and significant near the coast and near the southern border. For S95pTOT the trend is on the whole negative (even significantly so at some eastern stations), except for the coastal area. In analogy with Table 1, Table 2 summarizes the results for the two subregions. Given the dense clusters of stations in both subregions, the results in Table 2 may strongly influence those for northern Europe in Table 1. Large percentages of stations with an increase in R95pTOT are found in DJF andMAMfor the NetGer subregion, as well as in JJAfor the NorSwe subregion. These increases again appear to be related to the increases in mw. For DJF, the majority of the stations in the NorSwe subregion show a decrease in CVw and S95pTOT, which is field significant for the empirical estimate of S95pTOT. The effect of this decrease on R95pTOT is compensated by an increase in mw, found at 82.6%of the stations in this subregion. The increase of fw in DJF in northern Europe can, at least partly, be ascribed to region NorSwe, while the decrease inMAMand the slight increase in SON are accounted for by NetGer. None of the changes in fw are field significant, as was also concluded for the whole of northern Europe, where these changes are partly averaged out. Figure 5 gives an example of how the values of R95pTOT and S95pTOT compare for an individual station. The station Enonkoski Simanala (Finland) was chosen because of the high MK statistic found for the index R95pTOT in the DJF season. Transformed to the standard normal scale, this value was 4.08 (highly significant rise) and the nonparametric estimate of the slope (Sen 1968) amounted to 4.5%decade21. On the contrary, the index S95pTOT calculated from the same daily data using fitted Weibull distributions resulted in an MK statistic of 1.43 (not significant at the 5% level) and a slope of 0.6%decade21. Both index series are displayed in the top panel of Fig. 5. The first half of the series of R95pTOT contains a number of zeroes. These correspond with seasons in which the highest daily precipitation amounts fail to exceed the climatological 95th percentile. The zeroes do by definition not occur in the S95pTOT series.Apart fromthat effect, the year-to-year variations are smaller for the S95pTOT series than for R95pTOT (also seen in Fig. 2), enhancing the detection of a trend if present. The bottompanel of Fig. 5 shows the difference R95pTOT 2 S95pTOT (leftordinate) versus time together with the mean wet-day precipitation (right ordinate). The strong correspondence between both sequences confirms the conclusion drawn from the tables that the trend in mw can explain much of the difference between the trends in R95pTOT and S95pTOT. 5. Conclusions and discussion The objective of this study was to assess changes in the relative contribution of very wet days (i.e., days with precipitation amounts exceeding the 95th percentile) to the total precipitation amounts. Particular attention was given to the detection of a disproportional increase of extreme precipitation amounts, relative to the total precipitation. The index R95pTOT has often been used to monitor such changes. However, the use of this index has been questioned because of its strong year-to-year variations (Zolina et al. 2009). In this paper it is shown that R95pTOT is also influenced by changes in the mean wet-day precipitation. An alternative index S95pTOT was introduced, exhibiting less year-to-year variation and being better suited to characterizing a disproportional change of precipitation extremes. Assuming the wet-day precipitation amounts are Weibull distributed, an analytical expression for S95pTOT was derived in terms of the Weibull shape parameter. For the quality-checked daily precipitation series of more than 1800 stations across Europe covering the period 1961- 2010 the seasonal series of R95pTOT and S95pTOT were analyzed. For all four seasons we compared the number of positive trends in R95pTOT and S95pTOT, the fraction of wet days, and the mean and the CV of the wet-day precipitation amounts over the period 1961-2010. The reason that R95pTOT is sensitive to the changes in the mean wet-day precipitation lies in the fact that a fixed, climatological 95th percentile is used to define very wet days. We argued that a trend in the mean wetday precipitation then contributes in a positive sense to the trend in R95pTOT, which was not fully recognized in earlier interpretations of R95pTOT as a measure of ''disproportional change in the extremes.'' An additional problem in the definition of R95pTOT is that year-to-year variations in the number of wet days and mean wet-day precipitation induce strong variations and zeroes in this index. These can also affect the observed trend in R95pTOT. The proposed index S95pTOT is based on a 95th percentile that is not assumed constant over time. Evaluation of S95pTOT using the proposed analytical expression does not require an explicit estimate of the 95th percentile but only the Weibull shape parameter. The contrast between the trends in northern and southern Europe was far less pronounced for S95pTOT than for R95pTOT. For many stations the trend in both indices even has an opposite sign. The net difference with the trend in R95pTOT was a decreased number of significant positive trends and amoderate increase in the number of negative trends. From the percentages of positive trends it is concluded that in Europe the trend in R95pTOT is related to the trend in the mean, whereas the trend in S95pTOT seems to bemore influenced by the trend in the CV of the wet-day precipitation amounts. The fact that the climatological 95th percentile over the reference period is not representative of the 95th percentile of an individual season is also the real cause behind the zeroes in the R95pTOT series reported by Zolina et al. (2009). Their parametric approach actually avoids these zeroes by tackling the discrete nature of the calculation ofR95pTOT, so that in those casesR95pTOT becomes very small rather than zero. However, because they also obtain the 95th percentile from a reference period, their trends more or less resemble those of R95pTOT. In this study, rather than presenting an alternative calculation method, we have presented an alternative index with its own definition and interpretation. For seasons in which either the number of valid daily precipitation amounts drops below 80 or the number of wet days falls below 10, S95pTOT was not calculated. This has no serious effects on the S95pTOT series for the northern stations. For southern Europe, the lack of wet days reduces the number of stations for which S95pTOT is calculated in the JJA season by about 20%. This number does not show a distinctive trend, which suggests no systematic influence on the sign of trends in S95pTOT. However, missing values in the series of S95pTOT reduce the power of the tests for trend. For seasons having very few wet days, one might question the relevance and interpretation of such an index anyway. For seasons having sufficient wet days, S95pTOT clearly is meaningful. Though R95pTOT is also undeniably a meaningful index, one should be careful with its interpretation. The 95th percentile determined over a climatological period identifies precipitation events that are perceived as ''severe'' using the climatological period as a reference for a certain location. This makes R95pTOT very useful within a climate change impact context. However, to infer if extremes increase disproportionally and characterize changes in the distributional shape, othermeasures are more appropriate. The index S95pTOT introduced here is one of them. Acknowledgments. We acknowledge the data providers in theECA&D project (http://www.ecad.eu). The research leading to this paper has received funding from the European Union Seventh Framework Programme (FP7/2007-2013) under Grant Agreement 242093. The authors thank Dr. E. J. Klok, who initiated this study. REFERENCES Abramowitz, M., and I. A. Stegun, 1965: Handbook of Mathematical Functions. Dover, 1046 pp. Alexander, L. V., and Coauthors, 2006: Global observed changes in daily climate extremes of temperature and precipitation. J. Geophys. Res., 111, D05109, doi:10.1029/2005JD006290. Boes, D., J.-H. Heo, and J. D. Salas, 1989: Regional flood quantiles estimation for a Weibull model. Water Resour. Res., 25, 979-990. Donat, M. G., and Coauthors, 2013: Updated analysis of temperature and precipitation extreme indices since the beginning of the twentieth century: The HadEX2 dataset. J. Geophys. Res., 118, 2098-2118, doi:10.1002/jgrd.50150. Douglas, E. M., R. M. Vogel, and C. N. Kroll, 2000: Trends in floods and low flows in the United States: Impact of spatial correlation. J. Hydrol., 240, 90-105. Duan, J., J. Selker, and G. E. Grant, 1998: Evaluation of the probability density functions in precipitation models for the Pacific Northwest. J. Amer. Water Resour. Assoc., 34, 617-627, doi:10.1111/j.1752-1688.1998.tb00959.x. Ducic, V., J. Lukovic, D. Buric, G. Stanojevic, and S. Mustavic, 2012: Precipitation extremes in the wettest Mediterranean region (Krivosije) and associated atmospheric circulation types. Nat. Hazards Earth Syst. Sci., 12, 687-697, doi:10.5194/ nhess-12-687-2012. Field, C. B., and Coauthors, 2012: Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation. Cambridge University Press, 582 pp. Frich, P., L. V.Alexander, P.Della-Marta, B. Gleason,M.Haylock, A. M.G. Klein Tank, and T. Peterson, 2002:Observed coherent changes in climatic extremes during the second half of the twentieth century. Climate Res., 19, 193-212. Goda, Y., M. Kudaka, and H. Kawai, 2010: Incorporation of the Weibull distribution in L-moments for regional frequency analysis of peaks-over-threshold wave heights. Proc. 32nd Int. Conf. on Coastal Engineering, Paper 59, Shanghai, China, American Society of Civil Engineers, 11 pp. Groisman, P. Ya., and Coauthors, 1999: Changes in the probability of heavy precipitation: Important indicators of climate change. Climatic Change, 42, 243-283. Hamed, K. H., and A. R. Rao, 1998: A modifiedMann-Kendall trend test for autocorrelated data. J. Hydrol., 204, 182- 196. Hosking, J. R. M., 1990: L-moments: Analysis and estimation of distributions using linear combinations of order statistics. J. Roy. Stat. Soc., 52B, 105-124. _____, and J. R. Wallis, 1997: Regional Frequency Analysis: An Approach Based on L-Moments. Cambridge University Press, 224 pp. Kiktev, D., D. M. H. Sexton, L. Alexander, and C. K. Folland, 2003: Comparison of modeled and observed trends in indices of daily climate extremes. J. Climate, 16, 3560-3571. Klein Tank, A. M. G., and G. P. Keuroonnen, 2003: Trends in indices of daily temperature and precipitation extremes in Europe, 1946-99. J. Climate, 16, 3665-3680. _____, and Coauthors, 2002: Daily dataset of 20th century surface air temperature and precipitation series for the European Climate Assessment. Int. J. Climatol., 22, 1441- 1453. Nicholls, N., and W. Murray, 1999: Workshop on Indices and Indicators for Climate Extremes: Asheville, NC, USA, 3-6 June 1997 breakout group B: Precipitation. Climatic Change, 42, 23-29. Peterson, T. C., and M. J. Manton, 2008: Monitoring changes in climate extremes: A tale of international collaboration. Bull. Amer. Meteor. Soc., 89, 1266-1271. Sen, P. K., 1968: Estimates of the regression coefficient based on Kendall's tau. J. Amer. Stat. Assoc., 63, 1379-1389. Sillmann, J.,V.V.Kharin, X. Zhang, F.W.Zwiers, andD.Bronaugh, 2013: Climate extreme indices in the CMIP5 multi-model ensemble: Part 1: Model evaluation in the present climate. J. Geophys. Res. Atmos., 118, 1716-1733, doi:10.1002/ jgrd.50203. Trenberth, K. E., 2011: Changes in precipitation with climate change. Climate Res., 47, 123-138, doi:10.3354/cr00953. Turco, M., and M. C. Llasat, 2011: Trends in indices of daily precipitation extremes in Catalonia (NE Spain), 1951-2003. Nat. Hazards Earth Syst. Sci., 11, 3213-3226, doi:10.5194/ nhess-11-3213-2011. Vogel, R. M., and N. M. Fennessey, 1993: L moment diagrams should replace product moment diagrams. Water Resour. Res., 29, 1745-1752. Wijngaard, J. B., A. M. G. Klein Tank, and G. P. Keuroonnen, 2003: Homogeneity of 20th century European daily temperature and precipitation series. Int. J. Climatol., 23, 679-692. Wilby, R. L., and T. M. L. Wigley, 2002: Future changes in the distribution of daily precipitation totals across North America. Geophys. Res. Lett., 29, 1135, doi:10.1029/2001GL013048. Yue, S., P. Pilon, and G. Cavadias, 2002a: Power of the Mann- Kendall and Spearman's rho tests for detecting monotonic trends in hydrological series. J. Hydrol., 259, 254-271. _____, _____, B. Phinney, and G. Cavadias, 2002b: The influence of autocorrelation on the ability to detect trend in hydrological series. Hydrol. Processes, 16, 1807-1829. Zhang, X., L. V. Alexander, G. C. Hegerl, P. Jones, A. Klein Tank, T. C. Peterson, B. Trewin, and F. W. Zwiers, 2011: Indices for monitoring changes in extremes based on daily temperature and precipitation data. Wiley Interdiscip. Rev.: Climate Change, 2, 851-870, doi:10.1002/wcc.147. Zolina, O., C. Simmer, K. Belyaev, A. Kapala, and S. Gulev, 2009: Improving estimates of heavy and extreme precipitation using daily records from European rain gauges. J. Hydrometeor., 10, 701-716. R. LEANDER, T. A. BUISHAND, AND A. M. G. KLEIN TANK Royal Netherlands Meteorological Institute (KNMI), De Bilt, Netherlands (Manuscript received 28 February 2013, in final form 21 September 2013) Corresponding author address: A. M. G. Klein Tank, Royal Netherlands Meteorological Institute (KNMI), P.O. Box 201, 3730 AE, De Bilt, Netherlands. E-mail: [email protected] APPENDIX A Relative Contribution of the Upper 5% of the Distribution The density g is related to the density gw of the wetday precipitation amounts by g(x) 5 fwgw(x) for x $ d, where fw is the wet-day frequency. Equation (1) can then be rewritten as ... (A1) The factor fwmw/m gives the fraction of the precipitation total contributed by the wet days. Even with a relatively high wet-day threshold d of 1 mm, this factor is virtually equal to one. Equation (A1) then results in Eq. (2). Replacing x by y 1 d, Eq. (5) gives the following for RS95pTOT: ... (A2) The integral in the first term on the right-hand side represents Pr(X . q) 5 0.05, because q denotes the 95th percentile of the distribution of X. Substituting for gw the density of the shifted Weibull distribution [Eq. (7)] in the second term on the right-hand side, we have ... (A3) where ... denotes the upper incomplete gamma function [Abramowitz and Stegun 1965, Eq. (6.5.3)]. Since q is the 95th percentile of the shiftedWeibull distribution [Eq. (6)] with parameters a and c, ... (A4) Combining Eqs. (A2), (A3), and (A4) leads to ... (A5) Substitution of mw 5d1aG[(1/c)11] finally yields ... (A6) where ... refers to the normalized lower incomplete gamma function [Abramowitz and Stegun 1965, Eq. (6.5.1)] and ... refers to the complete gamma function. APPENDIX B Choice of the Weibull Distribution The two-parameter gamma distribution has been frequently used to model wet-day precipitation amounts (e.g., Groisman et al. 1999; Wilby and Wigley 2002;Zolina et al. 2009). For daily precipitation data from the U.S. Pacific Northwest, however, Duan et al. (1998) showed that the diagram of the sample L skewness versus the sample L coefficient of variation (L-CV) was more consistent with the theoretical relation between the L skewness and L-CV for theWeibull distribution than with that for the gamma distribution. The L-CVand the L skewness are alternatives to the conventional CV and coefficient of skewness, based on linear combinations of the ranked observations. These L-moment ratios are less biased and more robust to outliers than the conventional moment ratios. L-moment ratio diagrams have been used to select a probability distribution, in particular within the hydrological context. They were introduced by Hosking (1990) and their preference to conventional moment ratio diagrams was further demonstrated by Vogel and Fennessey (1993). The theoretical L-CV t2 and L skewness t3 for the two-parameterWeibull distribution can be obtained from Goda et al. (2010) as ... (B1) These expressions can be combined to obtain the relation ... (B2) For the two-parameter gamma distribution, as a special case of the Pearson type-III distribution, it follows from Hosking and Wallis (1997) that ... (B3) where c is the gamma shape parameter and Ix(, ) is the normalized incomplete beta function [Abramowitz and Stegun 1965, Eq. (6.6.2)]. For the latter, Hosking and Wallis (1997) also provide a rational-function approximation (not shown here), which was used to plot the relation between t3 and t2 for the gamma distribution. In Fig. B1 the sample L skewness of the wet-day precipitation is plotted versus the sample L-CV for all stations and all individual DJF seasons, together with the theoretical curves for the Weibull and gamma distributions. Note that the wet-day precipitation amounts were reduced by the wet-day threshold d to comply with the zero lower bound of the two theoretical distributions. It can be seen that for these data the Weibull distribution performs at least as well as the gamma distribution, though the slope is still somewhat underestimated. The L-moment ratio diagrams for the other seasons are similar. APPENDIX C Testing the Estimation of S95pTOT To test the performance of the two estimators for S95pTOT, a set of 50 000 simulations was conducted. In each simulation 100 values were generated from a Weibull distribution with shape parameter c 5 0.6 and scale parameter a 5 1.0. A sample was formed by accepting values at random with a probability of 60%. This incorporates the effect of a varying number of wet days in a season as a source of uncertainty. From each sample S95pTOT was approximated empirically, as well as through the shape parameter of a Weibull distribution fitted to the sample. Figure C1 compares the histograms of the two estimators. The dashed vertical line signifies the theoretical value associated with the chosen shape parameter. From the histograms we can conclude that both estimators are nearly unbiased. Furthermore, the use of the Weibull shape parameter seems to narrow the histogram, compared to that of the empirical estimates: for example, the width of the equitailed 90% interval reduces from approximately 28% for the empirical estimates to 19% for the estimates based on the Weibull fits. This is in line with the smaller interannual standard deviation for the Weibull-based estimates, displayed in Fig. 2. (c) 2014 American Meteorological Society |