Next Article in Journal
Estimation of Heavy Metal Concentrations in the Water of Urban Lakes in the Russian Arctic (Murmansk)
Next Article in Special Issue
Impacts of Rainfall Data Aggregation Time on Pluvial Flood Hazard in Urban Watersheds
Previous Article in Journal
Advanced Modeling Strategies for Hydraulic Engineering and River Research
Previous Article in Special Issue
A Review on Rainfall Data Resolution and Its Role in the Hydrological Practice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Time-Resolution of Rainfall Data on Trend Estimation for Annual Maximum Depths with a Duration of 24 Hours

1
Department of Civil and Environmental Engineering, University of Perugia, via G. Duranti 93, 06125 Perugia, Italy
2
National Research Council, Research Institute for Geo-Hydrological Protection, via Madonna Alta 126, 06128 Perugia, Italy
*
Author to whom correspondence should be addressed.
Water 2021, 13(22), 3264; https://doi.org/10.3390/w13223264
Submission received: 11 October 2021 / Revised: 9 November 2021 / Accepted: 16 November 2021 / Published: 17 November 2021

Abstract

:
The main challenge of this paper is to demonstrate that one of the most frequently conducted analyses in the climate change field could be affected by significant errors, due to the use of rainfall data characterized by coarse time-resolution. In fact, in the scientific literature, there are many studies to verify the possible impacts of climate change on extreme rainfall, and particularly on annual maximum rainfall depths, Hd, characterized by duration d equal to 24 h, due to the significant length of the corresponding series. Typically, these studies do not specify the temporal aggregation, ta, of the rainfall data on which maxima rely, although it is well known that the use of rainfall data with coarse ta can lead to significant underestimates of Hd. The effect of ta on the estimation of trends in annual maximum depths with d = 24 h, Hd=24 h, over the last 100 years is examined. We have used a published series of Hd=24 h derived by long-term historical rainfall observations with various temporal aggregations, due to the progress of recording systems through time, at 39 representative meteorological stations located in an inland region of Central Italy. Then, by using a recently developed mathematical relation between average underestimation error and the ratio ta/d, each Hd=24 h value has been corrected. Successively, commonly used climatic trend tests based on different approaches, including least-squares linear trend analysis, Mann–Kendall, and Sen’s method, have been applied to the “uncorrected” and “corrected” series. The results show that the underestimation of Hd=24 h values with coarse ta plays a significant role in the analysis of the effects of climatic change on extreme rainfalls. Specifically, the correction of the Hd=24 h values can change the sign of the trend from positive to negative. Furthermore, it has been observed that the innovative Sen’s method (based on a graphical approach) is less sensitive to corrections of the Hd values than the least-squares linear trend and the Mann–Kendall method. In any case, the analysis of Hd series containing potentially underestimated values, especially when d = 24 h, can lead to misleading results. Therefore, before conducting any trend analysis, Hd values determined from rainfall data characterized by coarse temporal resolution should always be corrected.

1. Introduction

It is well known that climate change is mainly due to greenhouse gas emissions from human activities [1]. One of the most important consequences is the modification of the hydrologic cycle with significant implications for water resources [2,3,4,5]. In the last century, mean global surface temperatures showed an increase of approximately 1.1 °C [6] and, based on the Clausius–Clapeyron relation, for each 1 °C increase in global temperature, the precipitable water increases by ~7% [7,8], even though relative humidity appears to decrease at high temperatures [1,8,9,10]. Moreover, it is expected that temperature will increase near to the surface and will decrease in the upper troposphere, favoring atmospheric instability [11]. Considering that atmospheric warming and water vapor trends also have local non-uniformity, the associated variation in average rainfall is typically mutable over the planet [12].
Extreme precipitation is very erratic, and its trends are less spatially coherent than those in average rainfall. On a global scale, there are significant unevenness and places where heavy rainfall is increasing, seemingly prevailing over areas where they are decreasing [13]. For example, [14] evaluated the long-term (1950–2018) trends in daily precipitation extremes over more than 7000 stations, and found that 66% of them showed increasing trends, and the remaining 34% showed decreasing trends. Furthermore, approximately 10% of stations (mainly in Europe, North America, and South Africa) showed a statistically significant increasing trend, while only 2.1% of stations (in the western United States, Canadian Prairies, and northern China) showed a significant decreasing trend. Reference [15] studied changes in daily precipitation over the Australian continent (period 1966–2013). They found an increase in daily precipitation during the second half of the observed period, as compared to the first half. Considering the period 1950–2015 and using daily gridded rainfall data (at 0.25-degree spatial resolution), [16] showed a significant rise in extreme precipitation events over central India. Increasing daily rainfall trends have been found also in west China [17], Bangladesh [18], South Korea [19], southern west Africa [20], and parts of southeast Asia [21].
Clear evidence that temporal variations in the occurrence of extreme rainfall events can be due to large-scale atmospheric and oceanographic oscillations, such as the North Atlantic Oscillation (NAO) or El Niño-Southern Oscillation (ENSO), is described in the scientific literature (e.g., [22,23,24,25,26,27,28,29,30]). Indicative analysis by [31] showed cyclic variations with a period of 30–40 years, and concluded that trends over time periods of <40 years could be ascribed to large-scale atmospheric changes.
Therefore, the analysis of trends in annual maximum rainfall depths, Hd, for a given duration, d, should be performed only for long-term rainfall data recorded, for example, from the earliest decades of the last century (see also [32]). Furthermore, analyses on the Hd series for d < 1 h are rarely available since, in the last century, all rainfall data have been recorded by adopting different temporal aggregations (or time resolutions), ta, dependent on changes in the recording systems through time. Currently, rainfall amounts are measured by tipping bucket sensors and recorded in a data-logger for each tip-time associated with a fixed rainfall depth, but until the last decades of the 21st century they were recorded only over paper rolls (pluviograph), generally with hourly ta [33,34]. In addition, for many years, especially before the Second World War, only daily rainfalls are available, recorded daily at a specific local time, and measuring the accumulated depth during the previous 24 h [35].
On this basis it can be deduced that, before the advent of data-loggers, rainfall data were always characterized by coarse temporal aggregation, with probable effects on analyses based on their use [36,37,38,39,40,41,42,43,44,45,46,47]. In fact, in some cases, the correct values of Hd can be significantly underestimated up to 50%, especially when d = 1 h and 24 h, due to the high probability of the presence of values with ta/d = 1 [48]. Moreover, long series of Hd values, together with a percentage of values obtained from continuous data (more recently recorded), always contain a percentage of elements derived from data characterized by coarse temporal aggregation, that are therefore potentially underestimated. This is problematic since, as well as the replacement of stations, the use of various rain gauge types with time and the change of surrounding near the equipment, could produce important effects on associated analyses, including the determination of rainfall depth-intensity-frequency curves [47] and trend evaluation of intense rainfalls [48].
By using some of the more common climatic trend tests characterized by very different approaches (least-squares linear trend analysis, Mann–Kendall test, and Sen’s method) the main objective of this paper is to evaluate the effect of time-resolution of rainfall data on trend estimation for annual maximum depths with duration 24 h. We focus our attention on the Hd series with d = 24 h as they are among the longest and most frequently available series for analyses of climate trends. However, similar analyses could be carried out for different d.
Furthermore, another significant challenge that we launch with this paper, together with the special issue that we promoted, consists in stimulating similar analyses conducted in different geographical areas of the world.

2. Study Area and Rainfall Data

The study area (Umbria region, with a surface area of 8456 km2) is located in an inland zone of central Italy, and is characterized by a complex orography along the eastern boundary, where the Apennine Mountains exceed 2000 m a.s.l. In the central and western areas, orography is mainly of hilly type, with elevations ranging from 100 to 800 m a.s.l.. A wide percentage of the study area is included in the basin of Tiber River that crosses the region from north to south-west receiving water from many tributaries, mainly located on the hydrographic left side.
On the basis of observations made by the rain gauge network shown in Figure 1 and specified in Table 1, annual rainfall depth through the region ranges from 650 mm to 1450 mm, with mean value of about 900 mm. Higher monthly rainfall values generally occur during the autumn-winter period, with floods caused by widespread rainfall.
For each selected station, we considered all the Hd=24 h values already obtained and validated from the Regional Hydrographic Service (RHS) by using the available rainfall data. It is important to note that recently, mainly since 1992, it has been possible to obtain rainfall data recorded in data-loggers for each tip time associated with a fixed rainfall depth (no more than 0.2 mm). In this case, each rainfall event was summarized by aggregating the number of tips over a ta equal to 1 min. Nevertheless, for long time series, a considerable amount of rainfall data was available with hourly recording system, as a consequence of the paper rolls adoption. Furthermore, in the absence of other possibilities, daily information derived from direct observation made each day at 9:00 a.m. was used. Two examples of the Hd=24 h series used in this paper are shown in Table 2 and Table 3 for the Gubbio and Todi rain gauge stations, respectively. For all analyses carried out in this paper, 39 Hd=24 h time series for the rain gauge stations reported in Table 1 were selected.

3. Methods

Each of the selected Hd=24 h series, hereinafter referred to as “uncorrected” because they can contain some underestimated values (see also [47]), is used to verify the existence of possible trend produced by climate change. Specifically, we considered the following tests, selected as they are characterized by very different approaches: (1) least-squares linear trend analysis; (2) non-parametrical Mann–Kendall test [49,50]; and (3) Sen’s method [51].
The least-squares method uses a straight line in order to fit the given points and it is known as the method of linear or ordinary least squares. This line is found as the best fit from which the sum of squares of the distances from the points is minimized. Equations with certain parameters usually represent the results in this method. The method of least squares actually defines the solution for the minimization of the sum of squares of deviations or the errors in the result of each equation.
The non-parametrical Mann–Kendall test is commonly used in detecting trends of variables in many fields. Statistic S can be calculated by the following:
S = k = 1 n 1 j = k + 1 n s g n ( x j x k )
with:
s g n ( x j x k ) = { + 1 ,   i f   ( x j x k ) > 0 0 ,   i f   ( x j x k ) = 0 1 ,   i f   ( x j x k ) < 0
and where n is the length of the sample, 𝑥k and 𝑥j are from 𝑘 = 1, 2, …, n − 1 and 𝑗 = 𝑘 + 1, …, n. If n is bigger than 8, statistic S approximates to normal distribution. The mean of S is 0 and the variance of S can be acquired as follows:
v a r ( S ) = n ( n 1 ) ( 2 n + 5 ) 18
Then, the test statistic Z is denoted by:
Z = { S 1 v a r ( S ) ,   i f   S > 0 0 ,   i f   S = 0 S + 1 v a r ( S ) ,   i f   S < 0
If Z > 0, it indicates an increasing trend, and vice versa. Given a confidence level α, the sequential data would be supposed to experience statistically significant trend if |Z| > Z(1 − α/2), where Z(1 − α/2) is the corresponding value of P = α/2 following the standard normal distribution. In this study, a 0.05 confidence level was used.
Finally, the innovative trend analysis proposed by Sen (2012) is based on a sub-section time series plot on a Cartesian coordinate system. In this method, recently criticized by [52], a time series is divided into two equal parts that are separately sorted in ascending order. Then, the first sub-series is located on the X-axis, and the second sub-series is located on the Y-axis. If the investigated data are collected on the 1:1 line, then there is no trend. If data fall above the 1:1 line or below the 1:1 line, then an upward trend or a downward trend in the time series, respectively, exists [51].
We note that these three very common tests could be representative to many others. Furthermore, since they are characterized by very different approaches, they could potentially offer different results.
Successively, the above-mentioned tests were repeated on a new version of the same series (hereinafter referred to as “corrected”), where the underestimation error due to the coarse time resolution of historical rainfall data was eliminated/minimized by using an average correction identically applied to all Hd values characterized by the same ratio ta/d. We note that an alternative correction approach, based on a sophisticated stochastic method, would produce very similar results with respect to the simple deterministic approach used here [48]. Note that this last result cannot be easily generalized to any other analysis [10,53,54]. With the aim of applying the correction, the average underestimation errors, Ea%, have been determined through the relations proposed earlier by [47], expressed by:
E a % = 6.14 ( t a d ) 2 + 5.69 t a d [ % ]                       d     30   min  
E a % = 6.7 ( t a d ) 2 + 4.72 t a d [ % ]       30   min   < d < 180   min
E a % = 5.2 ( t a d ) 2 + 5.57 t a d [ % ]       d     180   min    

4. Results and Discussion

As previously mentioned, initially the Hd=24 h series were corrected through the use of Equation (5c) proposed by [47]. For correction, we could use any of the relationships proposed by [36,37,38,39,40,41,42,43,44,45,46], but there is no reason to assume that this choice could have a significant effect on the results we present. As an example, Table 4 shows the Hd=24 h series of the Spoleto station, both in the original version (“uncorrected”) and in the “corrected” version.
Then, we fitted a least-squares linear trend to the Hd=24 h in both versions, “uncorrected” and “corrected”, considering all 39 selected stations. For the “uncorrected” Hd=24 h series, the number of positive least-square linear trends, equal to 23, outnumbers the negative ones, equal to 16 (see also Table 5). Opposite results were obtained for the “corrected” Hd=24 h series, as cases with negative least-square linear trends become equal to 27, whereas the positive ones become 12.
Figure 2 shows a comparison between “uncorrected” and “corrected” annual maximum rainfall depths for duration 24 h regarding Terni station. As can be seen, the “uncorrected” series shows a positive trend, while the “corrected” series linear trend is negative. This is mainly due to the fact that many old Hd=24 h values were underestimated, due to being derived by daily observations, while in the last decades the original Hd=24 h values have been determined without errors by using rainfall data recorded in data-loggers for each tip time associated with a fixed rainfall depth equal to 0.1 or 0.2 mm. Then, from a geometric point of view, if older values increase due to corrections and the most recent ones remain unchanged, the linear regression will be characterized by a decreasing slope.
Figure 3 shows the same comparison of Figure 2 between “uncorrected” and “corrected” Hd=24 h series for a restricted number of representative stations. It can be seen that, in some cases, such as that of Figure 3h regarding San Savino station, the least-square linear trend was also negative before correction with Equation (5c), with slope of the linear regression equal to −0.0285 mm/year. However, after correction, the negative trend becomes exacerbated (slope of the linear regression equal to −0.1172 mm/year).
With the main purpose to obtain useful and intuitive graphic representations of the linear trend slopes for both “uncorrected” and “corrected” series, Figure 4 shows the frequency of a reasonable number of slope classes, while in Figure 5 positive and negative slope values have been located in the geographic position of each rain gauge station.
Figure 4 highlights that a Gaussian probability function could be adequate to represent the distribution of linear trend slope values independently if “uncorrected” or “corrected” series are considered, with the average value changing from positive (0.0287) in the case of the “uncorrected” series to negative (−0.0485) for the “corrected” series. As it can also be deduced in Table 5, the two standard deviations are almost indistinguishable.
Figure 5 highlights that, within the Umbria region, there are no geographic areas most affected by specific trends, both when considering the “uncorrected” or “corrected” series.
The analysis of the 39-time series was successively performed by the non-parametric Mann–Kendall test. For the “uncorrected” series, by using a significance level equal to 0.05, the percentage without significant trends is approximately 90%, with one series (Spoleto) characterized by a significant negative trend, and three series (Castelluccio di Norcia, Gubbio and Montelovesco) by a significant positive trend (see Table 5). After the correction, the Mann–Kendall test evidenced six cases with a significant negative trend (Abeto, Massa Martana, Perugia, San Gemini, San Savino and Spoleto), and no positive cases. Therefore, for the Mann–Kendall test, the transition from “uncorrected” and “corrected” Hd=24 h series produces different results and conclusions.
Similarly to that which has already emerged using the most classic analyses, Sen’s method also produces different results when applied to the “uncorrected” and “corrected” annual maximum rainfall depths series. As an example, Figure 6 shows the results of the method for both the “uncorrected” and “corrected” Hd=24 h series observed at San Savino station. As can be seen in Figure 6a, the “uncorrected” series shows a no trend condition, because all values are close to 1:1 line and cross it several times, while in Figure 6b, the “correct” series evidences a clear monotonic decreasing trend (details on the results interpretation can be found in [51]).
It is interesting to note Sen’s method appears less sensitive to the Hd underestimation, due to the rainfall coarse temporal aggregation than classical linear regression and Mann–Kendall methods. This depends by the graphical approach on which the method is based. The example of the series regarding Norcia station is emblematic. Moving from the “uncorrected” series to the “corrected” series, both the linear trend slope and the Mann–Kendall test statistic Z change from positive to negative values (see Table 5). Instead, as shown in Figure 7, Sen’s method suggests that both “uncorrected” and “corrected” series are non-monotonically increasing.
All the results we have presented refer to a region of central Italy. We believe that they can be generalized to any other geographical area of the world, since it is widely demonstrated that the problem of the Hd underestimation is dependent on the technological evolution of rainfall recording systems, with all places showing similarity. In any case, we hope that this work can stimulate other case studies.
Even outside of the scope of this paper, it can be stated that, recently, in a study region located in central Italy, extreme rainfall events characterized by duration 24 h have not clearly and significantly deviated to trends observed 50–80 years ago. Therefore, despite in this region the effects of climate change having produced a continuous decrease in the annual rainfall amount and a significant increase of annual and monthly air temperature, these proclamations that try to attribute to climatic changes also the responsibility for the worrying hydrological and geological instability of the regional territory are not at all justified.

5. Conclusions

In the scientific literature, it is possible to find many analyses conducted to assess whether climate change has influenced the trends of the annual maximum rainfall depths with a duration of 24 h. Although it is well known that the Hd values obtained from rainfall data characterized by coarse temporal resolution can be significantly underestimated, the aforementioned analyses have been always conducted without taking into account the data origin.
The main objective of this paper was to evaluate the effect of time-resolution of rainfall data on Hd=24 h series trend estimation.
A representative number of rain gauge stations working perfectly for approximately a century with available Hd=24 h series, validated by the Regional Hydrographic Service, was carefully chosen.
All selected series, referred to as “uncorrected” since there were certainly some data containing some underestimated values [47], were used to verify the existence of possible trend due to climate change. For this purpose, we considered the very common least-squares linear trend analysis, the non-parametrical Mann–Kendall test, and the Sen’s method. Successively, the same series were modified through the procedure proposed by [47], with the purpose of eliminating/minimizing the underestimations, obtaining the so-called “corrected” series.
The underestimation errors due to coarse time-resolution of rainfall data produce significant effects on the least-squares linear trend analysis (based on a geometric approach) of the Hd=24 h values. Specifically, since a prevalence of increasing and decreasing trends for “uncorrected” and “corrected” series, respectively, was observed, the correction procedure can change the sign of the trend. For the non-parametric Mann–Kendall test (based on a statistical approach) with a significance level 0.05, using the selected 39 “uncorrected” time series, only one case was characterized by a significant negative trend, while three cases exhibited a significant positive trend. After the corrections, cases with negative trends became six, and there were no cases with a positive trend. Finally, the innovative Sen’s method (based on a graphical approach) has been noted to be less sensitive to corrections of the Hd values than the least-squares linear trend and the Mann–Kendall method. However, its results on the Hd=24 h series were also affected by the correction.
Overall, it can be concluded that the analysis of Hd series containing potentially underestimated values, especially when d = 24 h, can lead to misleading results. Therefore, before conducting any trend analysis, Hd values determined from rainfall data characterized by coarse temporal resolution should always be corrected with appropriate procedures suggested by the scientific literature. This correction could be neglected in the analysis performed by using the temporal scale of rainfall records [55].
The analysis conducted in this work may have the limit of considering rainfall data from only one geographical area of the world. There is no reason to assume that choosing only the case study of Umbria instead of many different case studies could have a significant effect on the results we presented. In any case, in the future we hope to find in the scientific literature similar analyses from many other territories.

Author Contributions

Investigation, writing—original draft preparation and writing—review and editing, R.M., C.S., J.D., A.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was mainly financed by University of Perugia (Fondo Ricerca di Base 2019).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Koutsoyiannis, D. Revisiting the global hydrological cycle: Is it intensifying? Hydrol. Earth Syst. Sci. 2020, 24, 3899–3932. [Google Scholar] [CrossRef]
  2. Barnett, T.P.; Pierce1, D.W.; Hidalgo, H.G.; Bonfils, C.; Santer, B.D.; Das, T.; Bala, G.; Wood, A.W.; Nozawa, T.; Mirin, A.A.; et al. Human-induced changes in the hydrology of the Western United States. Science 2008, 319, 1080–1083. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Sharafati, A.; Pezeshki, E. A strategy to assess the uncertainty of a climate change impact on extreme hydrological events in the semi-arid Dehbar catchment in Iran. Theor. Appl. Climatol. 2020, 139, 389–402. [Google Scholar] [CrossRef]
  4. Hajani, E. Climate change and its influence on design rainfall at-site in new south wales state, Australia. J. Water Clim. Chang. 2020, 11, 251–269. [Google Scholar] [CrossRef]
  5. Homsi, R.; Sanusi Shiru, M.; Shahid, S.; Ismail, T.; Bin Harun, S.; Al-Ansari, N.; Chau, K.-W.; Mundher Yaseen, Z. Precipitation projection using a CMIP5 GCM ensemble model: A regional investigation of Syria. Eng. Appl. Comp. Fluid Mech. 2020, 14, 90–106. [Google Scholar] [CrossRef]
  6. IPCC. Global Warming of 1.5 °C. An IPCC Special Report on the Impacts of Global Warming of 1.5 °C above Pre-Industrial Levels and Related Global Greenhouse Gas Emission Pathways, in the Context of Strengthening the Global Response to the Threat of Climate Change, Sustainable Development, and Efforts to Eradicate Poverty; Masson-Delmotte, V., Zhai, P., Pörtner, H.-O., Roberts, D., Skea, J., Shukla, P.R., Pirani, A., Moufouma-Okia, W., Péan, C., Pidcock, R., et al., Eds.; World Meteorological Organization: Geneva, Switzerland, 2018; p. 32. [Google Scholar]
  7. Lenderink, G.; Van Meijgarrd, E. Increase in hourly precipitation extremes beyond expectations from temperature changes. Nat. Geosci. 2008, 1, 511–514. [Google Scholar] [CrossRef]
  8. Hardwick-Jones, R.; Westra, S.; Sharma, A. Observed relationships between extreme sub-daily precipitation, surface temperature, and relative humidity. Geophys. Res. Lett. 2010, 37, L22805. [Google Scholar] [CrossRef] [Green Version]
  9. Wang, J.-W.; Wang, K.; Pielke Sr, R.A.; Lin, J.C.; Matsui, T. Towards a robust test on North America warming trend and precipitable water content increase. Geophys. Res. Lett. 2008, 35, L18804. [Google Scholar] [CrossRef] [Green Version]
  10. Dimitriadis, P.; Koutsoyiannis, D.; Iliopoulou, T.; Papanicolaou, P. A global-scale investigation of stochastic similarities in marginal distribution and dependence structure of key hydrological-cycle processes. Hydrology 2021, 8, 59. [Google Scholar] [CrossRef]
  11. Kunkel, K. North American trends in extreme precipitation. Nat. Hazards 2003, 29, 291–305. [Google Scholar] [CrossRef]
  12. Re, M.; Barros, V.R. Extreme rainfalls in SE South America. Clim. Chang. 2009, 96, 119–136. [Google Scholar] [CrossRef]
  13. Fowler, H.J.; Ali, H. Changes in Extreme Rainfall Events under the Warming Climate. In Rainfall. Modeling, Measurement and Applications; Morbidelli, R., Ed.; Elsevier: Amsterdam, The Netherlands, 2021. [Google Scholar] [CrossRef]
  14. Sun, Q.; Zhang, X.; Zwiers, F.; Westra, S.; Alexander, L.V. A global, continental, and regional analysis of changes in extreme precipitation. J. Clim. 2021, 34, 243–258. [Google Scholar] [CrossRef]
  15. Guerreiro, S.B.; Fowler, H.J.; Barbero, R.; Westra, S.; Lenderink, G.; Blenkinsop, S.; Li, X.F. Detection of continental-scale intensification of hourly rainfall extremes. Nat. Clim. Chang. 2018, 8, 803–807. [Google Scholar] [CrossRef]
  16. Roxy, M.K.; Ghosh, S.; Pathak, A.; Athulya, R.; Mujumdar, M.; Murtugudde, R.; Rajeevan, M. A threefold rise in widespread extreme rain events over central India. Nat. Commun. 2017, 8, 1–11. [Google Scholar] [CrossRef] [PubMed]
  17. Wang, Y.; Zhou, L. Observed trends in extreme precipitation events in China during 1961–2001 and the associated changes in large-scale circulation. Geophys. Res. Lett. 2005, 32, L09707. [Google Scholar] [CrossRef]
  18. Shahid, S. Trends in extreme rainfall events of Bangladesh. Theor. Appl. Climatol. 2011, 104, 489–499. [Google Scholar] [CrossRef]
  19. Park, J.S.; Kang, H.S.; Lee, Y.S.; Kim, M.K. Changes in the extreme daily rainfall in South Korea. Int. J. Climatol. 2011, 31, 2290–2299. [Google Scholar] [CrossRef]
  20. Nkrumah, F.; Vischel, T.; Panthou, G.; Klutse, N.A.B.; Adukpo, D.C.; Diedhiou, A. Recent trends in the daily rainfall regime in southern West Africa. Atmosphere 2019, 10, 741. [Google Scholar] [CrossRef] [Green Version]
  21. Donat, M.G.; Alexander, L.V.; Herold, N.; Dittus, A.J. Temperature and precipitation extremes in century-long gridded observations, reanalyses, and atmospheric model simulations. J. Geophys. Res. Atmos. 2016, 121, 11174–11189. [Google Scholar] [CrossRef] [Green Version]
  22. Gershunov, A.; Cayan, D.R. Heavy daily precipitation frequency over the contiguous United States: Sources of climatic variability and seasonal predictability. J. Clim. 2003, 16, 2752–2765. [Google Scholar] [CrossRef] [Green Version]
  23. Schreck, C.J.; Semazzi, F.H.M. Variability of the recent climate of Eastern Africa. Int. J. Climatol. 2004, 24, 681–701. [Google Scholar] [CrossRef]
  24. Grimm, A.M.; Tedeschi, R.G. ENSO and extreme rainfall events in South America. J. Clim. 2009, 22, 1589–1609. [Google Scholar] [CrossRef]
  25. Haylock, M.R.; Peterson, T.C.; Alves, L.M.; Ambrizzi, T.; Anunciação, Y.M.T.; Baez, J.; Barroz, V.R.; Berlato, M.A.; Bidegain, M.; Coronel, G.; et al. Trends in total and extreme South America rainfall in 1960–2000 and links with sea surface temperature. J. Clim. 2006, 19, 1490–1512. [Google Scholar] [CrossRef]
  26. Cayan, D.R.; Redmond, K.T.; Riddle, L.G. ENSO and hydrologic extremes in the western United States. J. Clim. 1999, 12, 2881–2893. [Google Scholar] [CrossRef] [Green Version]
  27. Gershunov, A. ENSO influence on intraseasonal extreme rainfall and temperature frequencies in the contiguous United States: Implications for long-range predictability. J. Clim. 1998, 11, 3192–3203. [Google Scholar] [CrossRef]
  28. Aryal, S.K.; Bates, B.C.; Campbell, E.P.; Li, Y.; Palmer, M.J.; Viney, N.R. Characterizing and modelling temporal and spatial trends in rainfall extremes. J. Hydromet. 2009, 10, 241–253. [Google Scholar] [CrossRef]
  29. Kamruzzaman, M.; Beecham, S.; Metcalfe, A.V. Non-stationarity in rainfall and temperature in the Murray Darling Basin. Hydrol. Proc. 2011, 25, 1659–1675. [Google Scholar] [CrossRef]
  30. Kendon, E.J.; Blenkinsop, S.; Fowler, H.J. Will changes in short duration precipitation extremes be detectable before daily extremes? J. Clim. 2018, 31. [Google Scholar] [CrossRef]
  31. Willems, P. Adjustment of extreme rainfall statistics accounting for multidecadal climate oscillations. J. Hydrol. 2013, 490, 126–133. [Google Scholar] [CrossRef]
  32. Teegavarapu, R.S.V.; Nayak, A. Evaluation of long-term trends in extreme precipitation: Implications of in-filled historical data use for analysis. J. Hydrol. 2017, 550, 616–634. [Google Scholar] [CrossRef]
  33. Bonaccorso, B.; Cancelliere, A.; Rossi, G. Detecting trends of extreme rainfall series in Sicily. Adv. Geosci. 2005, 2, 7–11. [Google Scholar] [CrossRef] [Green Version]
  34. Fatichi, S.; Caporali, E. A comprehensive analysis of changes in precipitation regime in Tuscany. Int. J. Climatol. 2009, 29, 1883–1893. [Google Scholar] [CrossRef]
  35. Morbidelli, R.; García-Marín, A.P.; Al Mamun, A.; Atiqur, R.M.; Ayuso-Muñoz, J.L.; Bachir Taouti, M.; Baranowski, P.; Bellocchi, G.; Sangüesa-Pool, C.; Bennett, B.; et al. The history of rainfall data time-resolution in different geographical areas of the world. J. Hydrol. 2020, 590, 125258. [Google Scholar] [CrossRef]
  36. Hershfield, D.M.; Wilson, W.T. Generalizing of Rainfall-intensity-frequency Data. IUGG/IAHS Publication No. 43; AIHS. Gen. Ass. Tor. 1958, 1, 499–506. [Google Scholar]
  37. Hershfield, D.M. Rainfall Frequency Atlas of the United States for Durations from 30 Minutes to 24 Hours and Return Periods from 1 to 100 Years; US Weather Bureau Technical Paper N. 40; U.S. Department of Commerce: Washington, DC, USA, 1961. [Google Scholar]
  38. Weiss, L.L. Ratio of true to fixed-interval maximum rainfall. J. Hydraul. Div. 1964, 90, 77–82. [Google Scholar] [CrossRef]
  39. Harihara, P.S.; Tripathi, N. Relationship of the clock-hour to 60-min and the observational day to 1440-min rainfall. Ind. J. Meteorol. Geophys. 1973, 24, 279–282. [Google Scholar]
  40. Natural Environment Research Council. Flood Studies Report; Natural Environmental Research Council: London, UK, 1975. [Google Scholar]
  41. Van Montfort, M.A.J. Sliding maxima. J. Hydrol. 1990, 118, 77–85. [Google Scholar] [CrossRef]
  42. Faiers, G.E.; Grymes, J.M.; Keim, B.D.; Muller, R.A. A re-examination of extreme 24 hour rainfall in Louisiana, USA. Clim. Res. 1994, 4, 25–31. [Google Scholar] [CrossRef]
  43. Van Montfort, M.A.J. Concomitants of the Hershfield factor. J. Hydrol. 1997, 194, 357–365. [Google Scholar] [CrossRef]
  44. Young, C.B.; McEnroe, B.M. Sampling adjustment factors for rainfall recorded at fixed time intervals. J. Hydrol. Eng. 2003, 8, 294–296. [Google Scholar] [CrossRef]
  45. Yoo, C.; Park, M.; Kim, H.J.; Choi, J.; Sin, J.; Jun, C. Classification and evaluation of the documentary-recorded storm events in the Annals of the Choson Dynasty (1392–1910), Korea. J. Hydrol. 2015, 520, 387–396. [Google Scholar] [CrossRef]
  46. Papalexiou, S.M.; Dialynas, Y.G.; Grimaldi, S. Hershfield factor revisited: Correcting annual maximum precipitation. J. Hydrol. 2016, 524, 884–895. [Google Scholar] [CrossRef]
  47. Morbidelli, R.; Saltalippi, C.; Flammini, A.; Cifrodelli, M.; Picciafuoco, T.; Corradini, C.; Casas-Castillo, M.C.; Fowler, H.J.; Wilkinson, S.M. Effect of temporal aggregation on the estimate of annual maximum rainfall depths for the design of hydraulic infrastructure systems. J. Hydrol. 2017, 554, 710–720. [Google Scholar] [CrossRef] [Green Version]
  48. Morbidelli, R.; Saltalippi, C.; Flammini, A.; Corradini, C.; Wilkinson, S.M.; Fowler, H.J. Influence of temporal data aggregation on trend estimation for intense rainfall. Adv. Water. Resour. 2018, 122, 304–316. [Google Scholar] [CrossRef]
  49. Mann, H.B. Nonparametric tests against trend. Econometrica 1945, 13, 46–59. [Google Scholar] [CrossRef]
  50. Kendall, M.G. Rank Correlation Methods; Griffin: London, UK, 1975. [Google Scholar]
  51. Sen, Z. Innovative trend analysis methodology. J. Hydrol. Eng. 2012, 17, 1042–1046. [Google Scholar] [CrossRef]
  52. Serinaldi, F.; Chebana, F.; Kilsby, C.G. Dissecting innovative trend analysis. Stoch. Environ. Res. Risk Assess. 2020, 34, 733–754. [Google Scholar] [CrossRef]
  53. Iliopoulou, T.; Koutsoyiannis, D. Revealing hidden persistence in maximum rainfall records. Hydrol. Sci. J. 2019, 64, 1673–1689. [Google Scholar] [CrossRef]
  54. Iliopoulou, T.; Koutsoyiannis, D. Projecting the future of rainfall extremes: Better classic than trendy. J. Hydrol. 2020, 588. [Google Scholar] [CrossRef]
  55. Koutsoyiannis, D. Stochastics of Hydroclimatic Extremes—A Cool Look at Risk; Kallipos: Athens, Greece, 2021; 333p, ISBN 978-618-85370-0-2. [Google Scholar]
Figure 1. Morphology of the study area and rain gauges used in the analysis with identification numbers listed in Table 1. The background Digital Elevation Model (DEM), whose spatial resolution is 25 m, is derived by the Copernicus EU-DEM v1.1.
Figure 1. Morphology of the study area and rain gauges used in the analysis with identification numbers listed in Table 1. The background Digital Elevation Model (DEM), whose spatial resolution is 25 m, is derived by the Copernicus EU-DEM v1.1.
Water 13 03264 g001
Figure 2. Time sequence of “uncorrected” and “corrected” annual maximum rainfall depths for duration 24 h, Hd=24 h, with the respective linear trends, for Terni station. The red circles highlight periods where the Hd=24 h values have been obtained from rainfall data with the indicated aggregation time, ta. During the period 1936–2001 Hd=24 h values have been obtained from rainfall data with variable ta. Available data within the period 1921–2019.
Figure 2. Time sequence of “uncorrected” and “corrected” annual maximum rainfall depths for duration 24 h, Hd=24 h, with the respective linear trends, for Terni station. The red circles highlight periods where the Hd=24 h values have been obtained from rainfall data with the indicated aggregation time, ta. During the period 1936–2001 Hd=24 h values have been obtained from rainfall data with variable ta. Available data within the period 1921–2019.
Water 13 03264 g002
Figure 3. Time sequence of “uncorrected” and “corrected” annual maximum rainfall depths for duration 24 h, Hd=24 h, with the respective linear trends, for the following stations: (a) Umbertide; (b) Calvi dell’Umbria; (c) Ficulle; (d) Monteleone di Spoleto; (e) Narni Scalo; (f) Norcia; (g) Ponte Nuovo; and (h) San Savino. Available data within the period 1921–2019.
Figure 3. Time sequence of “uncorrected” and “corrected” annual maximum rainfall depths for duration 24 h, Hd=24 h, with the respective linear trends, for the following stations: (a) Umbertide; (b) Calvi dell’Umbria; (c) Ficulle; (d) Monteleone di Spoleto; (e) Narni Scalo; (f) Norcia; (g) Ponte Nuovo; and (h) San Savino. Available data within the period 1921–2019.
Water 13 03264 g003
Figure 4. Frequency of the linear trend slope for (a) “uncorrected” and (b) “corrected” Hd=24 h series considering classes amplitude equal to 0.04 mm/year. All selected rain gauge stations were considered.
Figure 4. Frequency of the linear trend slope for (a) “uncorrected” and (b) “corrected” Hd=24 h series considering classes amplitude equal to 0.04 mm/year. All selected rain gauge stations were considered.
Water 13 03264 g004
Figure 5. Values of the linear trend slopes for both “uncorrected” (left) and “corrected” (right) Hd=24 h series located in the corresponding geographic position of the rain gauge station.
Figure 5. Values of the linear trend slopes for both “uncorrected” (left) and “corrected” (right) Hd=24 h series located in the corresponding geographic position of the rain gauge station.
Water 13 03264 g005
Figure 6. Trend conditions according to Sen’s method for the annual maximum rainfall depths for San Savino station and duration equal to 24 h: (a) “uncorrected” series; and (b) “corrected” series. Available data within the period 1921–2018.
Figure 6. Trend conditions according to Sen’s method for the annual maximum rainfall depths for San Savino station and duration equal to 24 h: (a) “uncorrected” series; and (b) “corrected” series. Available data within the period 1921–2018.
Water 13 03264 g006
Figure 7. Trend conditions according to Sen’s method for the annual maximum rainfall depths for Norcia station and duration equal to 24 h: (a) “uncorrected” series; and (b) “corrected” series. Available data within the period 1921–2019.
Figure 7. Trend conditions according to Sen’s method for the annual maximum rainfall depths for Norcia station and duration equal to 24 h: (a) “uncorrected” series; and (b) “corrected” series. Available data within the period 1921–2019.
Water 13 03264 g007
Table 1. Main characteristics of the selected rainfall stations. The geographic position is expressed in Universal Transvers Mercator (UTM) coordinates computed by using the WGS84 ellipsoid model. ARDta is the percentage of available rainfall data characterized by specific temporal aggregation, ta.
Table 1. Main characteristics of the selected rainfall stations. The geographic position is expressed in Universal Transvers Mercator (UTM) coordinates computed by using the WGS84 ellipsoid model. ARDta is the percentage of available rainfall data characterized by specific temporal aggregation, ta.
ID
Number
Rain Gauge
Station
Altitude (m a.s.l.)UTM33 X
(m)
UTM33 Y
(m)
Available Data PeriodARDta (%)
ta = 1
min
ta = 60 minta = 1440 min
1Abeto946341,8054,744,5711951–20149.40.090.6
2Amelia321287,9594,714,8291921–201719.220.560.3
3Arrone221316,2894,716,8601921–201711.00.089.0
4Assisi408305,7994,771,4421921–20010.045.254.8
5Attigliano64277,4954,711,0221921–201511.80.088.2
6Bastia203301,3774,769,7161922–201731.80.068.2
7Bevagna212307,3704,757,3201921–201721.76.072.3
8Calvi dell’Umbria305299,1644,698,5611951–201720.40.079.6
9Cascia604338,4774,731,5921922–201719.60.080.4
10Castelluccio di Norcia1349354,0314,743,4091921–201712.70.087.3
11Castiglione del Lago260259,7604,779,5791921–201915.727.157.1
12Città di Castello304277,6434,815,7381921–201936.429.933.8
13Compignano240278,3944,758,5931922–201734.60.065.4
14Corciano306280,8714,776,2041921–201916.40.083.6
15Ficulle440260,1444,747,4801921–201510.80.089.2
16Foligno220310,6784,758,2251916–201525.731.442.9
17Gualdo Tadino599319,8704,789,9531921–201918.447.134.5
18Gubbio471302,7894,802,3291921–201920.042.537.5
19Lago di Corbara128273,6404,731,0141963–201934.00.066.0
20Massa Martana328297,4574,738,7411921–201924.64.970.5
21Monte del Lago260270,6574,780,2521923–201612.536.351.3
22Monteleone di Spoleto933331,8824,723,6181953–201923.70.076.3
23Montelovesco634290,4844,798,1421921–201941.20.058.8
24Narni Scalo109298,3814,713,9161921–201933.70.066.3
25Nocera Umbra534320,2814,776,4051921–201932.60.067.4
26Norcia691345,0424,740,1891921–201924.40.075.6
27Orvieto311263,1784,733,5591921–201521.346.132.6
28Perugia440288,0874,775,3491921–20193.843.852.5
29Petrelle342269,8304,803,5531921–201930.40.069.6
30Pianello233302,0034,779,6691921–201922.60.077.4
31Ponte Nuovo174290,4914,765,1441921–201921.77.271.1
32Prodo431273,7524,738,7901921–201713.60.086.4
33San Gemini299298,2754,720,3011921–201919.56.973.6
34San Savino260271,1704,776,4681921–201832.00.068.0
35Sellano604330,3074,750,4801951–201732.40.067.6
36Spoleto353314,9524,736,1621921–201920.740.239.1
37Terni123307,1234,714,6031921–201918.432.249.4
38Todi329288,0894,740,3191921–201929.840.429.8
39Umbertide305284,8674,798,8361921–201922.528.848.8
Table 2. Annual maximum rainfall depths (in mm) for duration (d) equal to 24 h, Hd=24 h, derived from data characterized by different aggregation times, ta. Gubbio rain gauge station.
Table 2. Annual maximum rainfall depths (in mm) for duration (d) equal to 24 h, Hd=24 h, derived from data characterized by different aggregation times, ta. Gubbio rain gauge station.
YearHd=24 hYearHd=24 hYearHd=24 hYearHd=24 h
192165194458.2197343.4199776.7
192249194672197440.4199862.7
192355.6194755197589.8199970.4
192445194861197672.4200063
192548194958197740200151.5
192656.4195031.7197854.2200252.6
192734.2195161197953.42003106.8
192855.7195241198066.8200455.8
192974.2195330.8198159.8200591.8
193083.4195455198277200683.6
193179.1195568.81984101200745.4
193257.2195650.5198544200865.6
193364.6195761198658.6200951
193471.4195851198768.4201078
193567.8195951.5198862.6201137.2
193659.81960127.2198991.42012131.4
1937741961761990782013104.2
193833.4196270199147.8201472.4
193944.8196352199260.6201554.2
194048196458.4199364.4201676.6
194144.61965115.4199487201744.8
194257196645.4199551.1201851.2
194360.8196868.2199678.1201982.2
Legend:
ta = 1 minute
ta = 1 hour
ta = 1 day
Table 3. Annual maximum rainfall depths (in mm) for duration (d) equal to 24 h, Hd=24 h, derived from data characterized by different aggregation times, ta. Todi rain gauge station.
Table 3. Annual maximum rainfall depths (in mm) for duration (d) equal to 24 h, Hd=24 h, derived from data characterized by different aggregation times, ta. Todi rain gauge station.
YearHd=24 hYearHd=24 hYearHd=24 hYearHd=24 h
192169.3195032197456.51998113.2
192256195177.3197565.6199988.2
192363.3195251.4197651.4200042.2
192439.4195350.2197741.8200154.6
192552.2195448.2197874.2200283.1
192646.6195545.8197953200341.7
192750.8195652.6198093.6200447.4
192861195748.4198129200570.6
192965195845198249.8200641.1
193037.2195957198363200736.8
193138.41960124198467.2200855.6
193256196182.8198543.4200983.3
193348.2196255198698.6201065.1
193469.6196388.5198763.2201137.5
1935140.4196490.8198860201267.5
1936207196584198945.8201345.8
193780196642.2199065.2201467.9
193865.8196764.2199146.8201539.5
193953.2196897199250.8201655.2
194061.41969103.4199384.1201747.8
194158197038.2199438.8201844.2
194247197136.51995137.5201958
194841.5197246.6199647.2
1949511973581997118.2
Legend:
ta = 1 minute
ta = 1 hour
ta = 1 day
Table 4. Annual maximum rainfall depths (in mm) for duration (d) equal to 24 h, Hd=24 h, in the original version (“uncorr.”) and in the modified version through Equation (5c) (“corr.”). Spoleto rain gauge station.
Table 4. Annual maximum rainfall depths (in mm) for duration (d) equal to 24 h, Hd=24 h, in the original version (“uncorr.”) and in the modified version through Equation (5c) (“corr.”). Spoleto rain gauge station.
YearHd=24 hYearHd=24 hYearHd=24 h
“Uncorr.”
Series
“Corr.”
Series
“Uncorr.” Series“Corr.”
Series
“Uncorr.”
Series
“Corr.”
Series
19217078.4195970.879.319895252.1
19226977.3196087.387.5199099.499.6
1923150168.1196172.172.319916161.1
1924142159.1196279.379.5199285.685.8
19256269.51963112112.2199360.860.9
19266067.21964113.6113.8199456.763.5
19275865.01965104104.2199558.565.6
192846.552.1196652.953.0199667.375.4
19297078.4196782.883.019976168.4
19307482.9196865.665.719997373.2
19315359.419697676.2200048.648.7
19326067.2197038.738.8200135.840.1
19337078.4197161.661.7200254.654.6
19346876.219725151.1200333.433.4
1935140156.9197350.450.520044343
193695106.5197442.642.7200599.499.4
19377381.81975107107.2200642.742.7
19394651.619769898.2200739.139.1
194058.465.4197774.374.5200845.645.6
194158.365.319785454.1200973.673.6
194968.468.5197948.153.9201052.652.6
195062.462.5198068.268.3201144.844.8
19518585.2198161.661.7201264.664.6
195271.271.3198263.263.3201366.466.4
195343.448.6198358.458.5201478.278.2
195456.663.419847373.2201560.860.8
195556.563.3198538.242.8201654.454.4
19568181.2198678.478.6201782.882.8
195760.968.3198762.662.7201879.279.2
195865.265.3198836.636.720196262
Legend:
Value modified through Equation (5c) with ta = 1 hour
Value modified through Equation (5c) with ta = 1 day
Table 5. Slope (in mm/year) of the least-squares linear regressions and Mann–Kendall test statistic Z of annual maximum rainfall depths for the selected stations and for duration 24 h, for the “uncorrected” and “corrected” series. For the Mann–Kendall test statistic Z, in bold the cases with significant trend.
Table 5. Slope (in mm/year) of the least-squares linear regressions and Mann–Kendall test statistic Z of annual maximum rainfall depths for the selected stations and for duration 24 h, for the “uncorrected” and “corrected” series. For the Mann–Kendall test statistic Z, in bold the cases with significant trend.
Rain Gauge
Station
Linear Trend Slope (mm/year)Mann-Kendall Test Statistic Z
“Uncorrected” Series“Corrected” Series“Uncorrected” Series“Corrected” Series
Abeto0.08150.01960.79−1.97
Amelia−0.0766−0.1913−0.59−1.45
Arrone0.0115−0.03330.29−0.19
Assisi−0.0071−0.1097−0.41−1.15
Attigliano−0.0212−0.0716−0.30−0.82
Bastia Umbra−0.0211−0.11930.16−0.99
Bevagna−0.0348−0.1158−0.84−1.81
Calvi dell’Umbria0.0227−0.09770.68−0.12
Cascia0.07510.01151.560.69
Castelluccio di Norcia0.42930.41501.971.69
Castiglione del Lago0.09070.01291.790.03
Città di Castello0.0413−0.01540.750.16
Compignano0.14580.05271.660.81
Corciano−0.0419−0.1482−0.55−1.13
Ficulle0.0363−0.02710.31−0.04
Foligno−0.0604−0.00590.34−0.64
Gualdo Tadino0.0237−0.02770.28−0.38
Gubbio0.16800.11142.161.25
Lago di Corbara0.25000.10261.460.84
Massa Martana−0.1288−0.2280−1.54−2.48
Monte del Lago0.06390.01570.04−0.64
Monteleone di Spoleto0.0347−0.09581.580.63
Montelovesco0.26390.18573.111.95
Narni Scalo0.0082−0.1076−0.32−1.21
Nocera Umbra0.11030.02171.05−0.18
Norcia0.0490−0.01780.42−0.63
Orvieto0.09760.05721.430.74
Perugia−0.1769−0.2129−1.72−1.99
Petrelle−0.0551−0.16040.22−0.91
Pianello0.11810.04281.760.83
Ponte Nuovo0.0099−0.06971.03−0.47
Prodo−0.0373−0.10390.38−0.17
San Gemini−0.1220−0.2188−1.60−2.36
San Savino−0.0285−0.1172−0.97−1.97
Sellano−0.0315−0.1575−0.30−1.22
Spoleto−0.2171−0.3159−2.08−3.16
Terni0.0368−0.05220.63−0.13
Todi−0.0473−0.1102−0.15−1.09
Umbertide0.0576−0.00830.58−0.15
average0.0287−0.04850.39−0.50
standard deviation0.11900.12811.151.16
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Morbidelli, R.; Saltalippi, C.; Dari, J.; Flammini, A. Effect of Time-Resolution of Rainfall Data on Trend Estimation for Annual Maximum Depths with a Duration of 24 Hours. Water 2021, 13, 3264. https://doi.org/10.3390/w13223264

AMA Style

Morbidelli R, Saltalippi C, Dari J, Flammini A. Effect of Time-Resolution of Rainfall Data on Trend Estimation for Annual Maximum Depths with a Duration of 24 Hours. Water. 2021; 13(22):3264. https://doi.org/10.3390/w13223264

Chicago/Turabian Style

Morbidelli, Renato, Carla Saltalippi, Jacopo Dari, and Alessia Flammini. 2021. "Effect of Time-Resolution of Rainfall Data on Trend Estimation for Annual Maximum Depths with a Duration of 24 Hours" Water 13, no. 22: 3264. https://doi.org/10.3390/w13223264

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop