{"id":20830,"date":"2017-11-21T11:53:11","date_gmt":"2017-11-21T16:53:11","guid":{"rendered":"http:\/\/www.realclimate.org\/?p=20830"},"modified":"2017-11-21T11:53:11","modified_gmt":"2017-11-21T16:53:11","slug":"a-brief-review-of-rainfall-statistics","status":"publish","type":"post","link":"https:\/\/www.realclimate.org\/index.php\/archives\/2017\/11\/a-brief-review-of-rainfall-statistics\/","title":{"rendered":"A brief review of rainfall statistics"},"content":{"rendered":"
\n

There have been a number of studies which show that we can expect more extreme rainfall with a global warming (e.g. Donat et al., 2016)<\/a><\/span>. Hence, there is a need to increase our resilience to more rainfall in the future. <\/p>\n

We can say something about how the rainfall statistics will be affected by a global warming, even when the weather itself is unpredictable beyond a few days. <\/p>\n

Statistics is remarkably predictable for a large number of events where each of them is completely random (welcome to thermodynamics and quantum physics). <\/p>\n

The normal distribution<\/a> has often been used to describe the statistical character of daily temperature, but it is completely unsuitable for 24-hr precipitation. Instead, the gamma distribution<\/a> has been a popular choice for describing rainfall. <\/p>\n

I wonder, however, if there is an even better way to quantify rainfall statistics.<\/p>\n

<\/p>\n

I have played around with the gamma distribution in an attempt to model daily rainfall statistics and its dependency on a set of physical factors. Without much success. <\/p>\n

However, then I noticed that most daily rain gauge appeared to be almost exponentially distributed<\/a> if I only included the rainy days (e.g. setting the threshold for a wet day at 1 mm). <\/p>\n

When I plotted the histogram for rainfall on wet days with a log-y axis, I would mostly get a straight line of dots (see a typical example below). <\/p>\n

\"\"<\/a>
Historgam of 24-hr precipitation measured at Bj\u00f8rnholt in a forest near Oslo. There will always be some clutter at the upper end of plots like these because there are so few data points representing these extreme values.<\/figcaption><\/figure>\n

The nice thing with the exponential distribution (which is a particular case of the gamma function) is that it only requires one<\/em> parameter to specify the mathematical curve: it’s the inverse of the mean value \"\mu\". <\/p>\n

I then used Bayes’ theorem<\/a> to account for dry and wet days, where the probability for rainfall was taken to be the wet-day frequency \"f_w\". <\/p>\n

The advantage of this approach is that I now had two<\/em> parameters which were easy to estimate: the wet-day mean precipitation (or mean rainfall intensity) \"\mu\" and the wet-day frequency \"f_w\". <\/p>\n

Furthermore, it turned out that\"f_w\" is often closely connected to the wind direction, and can easily be predicted based on circulation patterns or sea-level pressure anomalies. <\/p>\n

It was harder to find a systematic influence on \"\mu\", as it is likely affected by several factors, including the air moisture (which depends on temperature) and cloud top heights. <\/p>\n

The total precipitation is the product of \"n f_w \mu\", where \"n\" is the number of days. <\/p>\n

In other words, \"f_w\"\u00a0and \"\mu\" tell me many things I needed to know about the rainfall statistics (there are other aspects too, such as the mean duration of dry\/wet spells, the spatial extent, and whether it comes as rain, sleet, snow or hail).<\/p>\n

The equation for estimating the probability for a rain event with amounts exceeding \"x\" can be written as (using 1-CDF for the exponential distribution<\/a>):<\/p>\n

(1) <\/span>   <\/span>\"\begin{equation*}<\/p>\n

I have called it the “rain equation”, both because the name has not been taken and because it can provide many answers concerning rainfall. <\/p>\n

It can address questions about the likelihood of heavy rainfall and whether it is due to an increase in the number of rainy days (e.g. due to changes in circulation) or because the rains have become more intense. <\/p>\n

It is also on par with the normal distribution – in both cases, they are not meant to provide accurate probabilities for extreme events far out in the tails. <\/p>\n

However, they are both capable of quantifying the probability of more moderate values, which can be illustrated in the figure below:<\/p>\n

\"\"<\/a>
Figure 1. A comparison between probabilities estimated with the rain equation and the observed fraction of events with more than 30 mm rain in Groningen in the Netherlands. Here \"H(X - x)\" refers to the Heaviside function<\/a>, which is a mathematical way of expressing that I only counted the number of events with more than 30 mm\/day each year in the observervations (the plot was made with the R-package esd<\/a><\/code> and the command test.rainequation(loc='GRONINGEN-1',threshold=20)<\/code>).<\/figcaption><\/figure>\n

The rain equation captures long-term changes as well as inter-annual variations. In this example, I used the annual wet-day mean precipitation \"\mu\" and frequency \"f_w\" estimated from the observations themselves to show its potential. <\/p>\n

It can also be assessed against observations in a more systematic way, as in Figure 2:<\/p>\n

\"\"<\/a>
Figure 2. A scatter plot of probabilities and corresponding fractions of events from long rain gauge records in Europe, based on the wet-day mean precipitation and frequency from the observations (the plot was made with the R-package esd<\/a><\/code> and the command scatterplot.rainequation()<\/code>).<\/figcaption><\/figure>\n

A correlation of 0.98 is quite impressive, however, the rainfall is not perfectly exponentially distributed (Benestad et al., 2012)<\/a><\/span>. It nevertheless provides a means to address climate change<\/a> connected to a change in either \"f_w\" or \"\mu\". <\/p>\n

We have used the rain equation in an attempt to downscale seasonal and decadal forecasts for precipitation (Benestad and Mezghani, 2015)<\/a><\/span>. <\/p>\n

One thing that puzzles me, however, is that I cannot see this equation being used very much, despite the fact that it is so simple, seems so obvious, and can demonstrate impressive capabilities. <\/p>\n

I would have thought it is an old formula. Perhaps one that has gotten out of fashion, but is documented in old papers that are not yet digitized and easy to google. Perhaps with a different name. Or have I missed something?<\/p>\n

References<\/h2>\n
    \n
  1. <\/a>\nM.G. Donat, A.L. Lowry, L.V. Alexander, P.A. O\u2019Gorman, and N. Maher, \"More extreme precipitation in the world\u2019s dry and wet regions\", Nature Climate Change<\/i>, vol. 6, pp. 508-513, 2016. http:\/\/dx.doi.org\/10.1038\/nclimate2941<\/a>\n\n\n<\/li>\n
  2. <\/a>\nR.E. Benestad, D. Nychka, and L.O. Mearns, \"Spatially and temporally consistent prediction of heavy precipitation from mean values\", Nature Climate Change<\/i>, vol. 2, pp. 544-547, 2012. http:\/\/dx.doi.org\/10.1038\/nclimate1497<\/a>\n\n\n<\/li>\n
  3. <\/a>\nR.E. Benestad, and A. Mezghani, \"On downscaling probabilities for heavy 24-hour precipitation events at\n seasonal-to-decadal scales\", Tellus A: Dynamic Meteorology and Oceanography<\/i>, vol. 67, pp. 25954, 2015. http:\/\/dx.doi.org\/10.3402\/tellusa.v67.25954<\/a>\n\n\n<\/li>\n<\/ol>\n\n<\/div> ","protected":false},"excerpt":{"rendered":"

    There have been a number of studies which show that we can expect more extreme rainfall with a global warming . Hence, there is a need to increase our resilience to more rainfall in the future. We can say something about how the rainfall statistics will be affected by a global warming, even when the […]<\/p>\n","protected":false},"author":11,"featured_media":0,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"_exactmetrics_skip_tracking":false,"_exactmetrics_sitenote_active":false,"_exactmetrics_sitenote_note":"","_exactmetrics_sitenote_category":0,"_genesis_hide_title":false,"_genesis_hide_breadcrumbs":false,"_genesis_hide_singular_image":false,"_genesis_hide_footer_widgets":false,"_genesis_custom_body_class":"","_genesis_custom_post_class":"","_genesis_layout":"","footnotes":""},"categories":[1],"tags":[],"class_list":{"0":"post-20830","1":"post","2":"type-post","3":"status-publish","4":"format-standard","6":"category-climate-science","7":"entry"},"aioseo_notices":[],"post_mailing_queue_ids":[],"_links":{"self":[{"href":"https:\/\/www.realclimate.org\/index.php\/wp-json\/wp\/v2\/posts\/20830","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.realclimate.org\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.realclimate.org\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.realclimate.org\/index.php\/wp-json\/wp\/v2\/users\/11"}],"replies":[{"embeddable":true,"href":"https:\/\/www.realclimate.org\/index.php\/wp-json\/wp\/v2\/comments?post=20830"}],"version-history":[{"count":73,"href":"https:\/\/www.realclimate.org\/index.php\/wp-json\/wp\/v2\/posts\/20830\/revisions"}],"predecessor-version":[{"id":20913,"href":"https:\/\/www.realclimate.org\/index.php\/wp-json\/wp\/v2\/posts\/20830\/revisions\/20913"}],"wp:attachment":[{"href":"https:\/\/www.realclimate.org\/index.php\/wp-json\/wp\/v2\/media?parent=20830"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.realclimate.org\/index.php\/wp-json\/wp\/v2\/categories?post=20830"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.realclimate.org\/index.php\/wp-json\/wp\/v2\/tags?post=20830"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}