Assessing Inhomogeneities in Extreme Annual Rainfall Data Series by Multifractal Approach

: Testing the homogeneity in extreme rainfall data series is an important step to be performed before applying the frequency analysis method to obtain quantile values. In this work, six homogeneity tests were applied in order to check the existence of break points in extreme annual 24-h rainfall data at eight stations located in the Umbria region (Central Italy). Two are parametric tests (the standard normal homogeneity test and Buishand test) whereas the other four are non-parametric (the Pettitt, Sequential Mann–Kendal, Mann–Whitney U, and Cumulative Sum tests). No break points were detected at four of the stations analyzed. Where inhomogeneities were found, the multifractal approach was applied in order to check if they were real or not by comparing the split and whole data series. The generalized fractal dimension functions D q and the multifractal spectra f( α ) were obtained, and their main parameters were used to decide whether or not a break point existed.


Introduction
Frequency analysis of extreme rainfall data (defined as yearly maxima of a certain duration) is a common application in hydrologic engineering. It should be applied to recorded series that follow the characteristics of independency, stationarity and homogeneity [1]. Inhomogeneities in station data records can be due to observational routines (station relocations or changes in measuring techniques) and in that case, statistical methods and metadata information are effective for identifying these [2]. However, some shifts in climate dynamics are a consequence of human activities [3][4][5], or are related to El-Niño Southern Oscillation (ENSO) events such as those that occurred in 1976-77 [6], or to volcanic eruptions and changes in atmospheric composition and circulation [7].
Regardless of the reason, it is important to interpret past variability and recognize the shifts in climatic data in order to make adequate future projections. Therefore, it is crucial to investigate the effects of climate change on hydrologic time series, and especially on extreme rainfall data in order to accurately estimate or even update the intensity-duration-frequency (IDF) curves in a certain place [8]. The IDF curves are important tools for hydrologic and hydraulic design, which nowadays faces the problems of increased urbanization and population and of climate change [9], among others.
Numerous statistical methods such as the standard normal homogeneity (SNH) test, Buishand test, Pettitt test or Sequential Mann-Kendall test, among others, can be applied to evaluate the homogeneity of climate time series and have been widely used [13,18,19]. Some of the shifts detected under these classical tests turn out to be the regular behavior of time series under the scaling assumption [20,21]. Therefore, authors such as [22] recommend the consideration of scaling behavior of annual maximum rainfall when analyzing the trends in a time series, especially considering that the scaling hypothesis can increase the probability of rejection of correct null hypothesis [23]. Thus, accounting for scaling may help to avoid contradictions in trend analysis studies [23].
Authors such as [24] or [25] have analyzed changes in temporal data series during scale invariance, considering that several works have demonstrated the multifractal properties of meteorological systems [26][27][28]. Those systems exhibit the same properties for different measuring scales, thus certain behavioral information can be transferred from one scale to another, which is very useful for modeling purposes.
Sometimes, different results for the year when a shift or break point is detected in a data series are obtained by different statistical homogeneity tests, and metadata or other local and climate information cannot help to determine whether a break point really exists or not. In this work, we propose a novel use of the scale invariance properties of data sets as a useful tool to confirm whether inhomogeneities previously found in extreme rainfall data series are real or not. Moreover, the methodology can also be used to confirm the existence of a shift when only one or a few statistical tests have detected it. If a break point exists in an extreme rainfall series, the two data series derived by splitting the whole one before and after that break point should have different characteristics and can be considered as different. Therefore, the length decreases and can influence rainfall data estimation when applying frequency analysis to the two data series; this element is of great importance, especially for low return period quantile estimation, and so for IDF estimation.
For these purposes, extreme annual 24-h rainfall data from eight rain gauge stations with sufficiently long series were selected in the Umbria region (Central Italy). Six statistical homogeneity tests were applied to the selected extreme rainfall series in order to highlight inhomegeneities. The multifractal properties of the data series were then used to check the existence or not of break points detected by one or a few or all tests.

Data Source
Rainfall data from the Umbria region (central Italy) are used in this work. The Umbria region has an area of 8456 Km 2 , and is located in the Tiber river basin, which crosses the region from the north to the south-west. Its landscape is partially mountainous due to the presence of the Apennine Mountains that reach up to 2000 m a.s.l., and partially hilly in the central and western areas with altitudes ranging from 100 to 800 m a.s.l.
The mean annual rainfall for the last century is about 900 mm, with values varying in space from 650 mm to 1450 mm. The highest monthly rainfall values generally occur during the autumn-winter period, together with floods caused by widespread rainfall. The highest and lowest rainfall depths typically take place in November and July, respectively.  [29].
The study area is currently monitored through a dense rain gauge network (about 1 rain gauge every 90 Km 2 ). For this study, only the stations with at least 50 continuous years of extreme 24hour duration rainfall data have been considered (Table 1): Bastia Umbra, Bevagna, Nocera Umbra, Petrelle, Ponte Nuovo di Torgiano, Spoleto, Terni and Todi. Figure 1 shows the location of the eight sites considered in the Umbria region.

Homogeneity Tests
There are many statistical tests to check the homogeneity of a data series y1, y2, …, yn, with n the number of data, with specific characteristics. Some tests are known as parametric because they assume that the analyzed variable is normally distributed. Those that do not make this assumption are called non-parametric tests.
Under the null hypothesis, the annual values yi (i = 1, …, n) of the testing variable y are independent and identically distributed and the series is considered as homogeneous. Under the alternative hypothesis, the series shows a break in the mean and is considered as inhomogeneous. Moreover, there are tests that are able to give information about the year the break occurred, whereas those that assume the series is not randomly distributed under alternative hypothesis cannot provide this information.
In the present work, six homogeneity tests were applied. All of them are able to detect the year of the break if it exists. Two parametric tests were selected: the Buishand range (BR) test and the

Homogeneity Tests
There are many statistical tests to check the homogeneity of a data series y 1 , y 2 , . . . , y n , with n the number of data, with specific characteristics. Some tests are known as parametric because they assume that the analyzed variable is normally distributed. Those that do not make this assumption are called non-parametric tests.
Under the null hypothesis, the annual values y i (i = 1, . . . , n) of the testing variable y are independent and identically distributed and the series is considered as homogeneous. Under the alternative hypothesis, the series shows a break in the mean and is considered as inhomogeneous. Moreover, there are tests that are able to give information about the year the break occurred, whereas those that assume the series is not randomly distributed under alternative hypothesis cannot provide this information.
In the present work, six homogeneity tests were applied. All of them are able to detect the year of the break if it exists. Two parametric tests were selected: the Buishand range (BR) test and the standard normal homogeneity (SNH) test. The non-parametric tests were: the Pettitt (PT) test, Tequential Mann-Kendall test (SQMK) test, Mann-Whitney U (MWU) test and Cumulative Sum (CUSUM) test.
(1) Buishand Range test This test is based on the adjusted partial sums or cumulative derivations from the mean of the n data points, y: For a homogeneous record, S * k values fluctuate around zero. The values of S * k can be rescaled by the sample standard deviation, D y : with The Q, which is sensitive to departures from homogeneity, can then be obtained as: For different significance levels, the critical values for the test statistics (Qc) depend on the number of data and can be found in [30]. (2) Standard Normal Homogeneity test This test [31] compares the mean of the first k years of data with the last (n − k) years by using the T k of the statistics: where The year k shows a break if the value of T k is the maximum. If the T 0 = max 1≤k≤n T k is greater than the critical values [32], the null hypothesis is rejected. (3) Pettitt test The Pettitt test [33] is based on the ranking r i of the y i values. The ranks r i are obtained by ordering the data in crescent order, so that the smallest one gets ranked 1 and the highest gets the n-rank. The U k is then obtained as: If a break occurs in the year K, then the U k is: The statistical significance of the break point is checked by comparing the value of U K (Equation (7)) to its theoretical value: with α the significance level. (4) equential Mann-Kendall test This test is usually applied to evaluate trends in data series and also to check the moment when they start being significant [34][35][36]. It calculates two series of statistic values: one (U(t)) for the progressive temporal data set (y 1 , y 2 , . . . , y n ) and the other (U'(t)) for the corresponding regressive data set (y n , y n−1 , . . . , y 1 ), both with an average value of zero and a standard deviation of one. Usually, the ranks of the y i values are preferred. The t k has to be obtained by: with n i being the number of cases for which y i > y j ( j = 1, 2, . . . , i − 1). Under the null hypothesis, it is assumed that no trend exists, and t k follows a normal distribution with average and variance values given by: The sequential values of U(t k ) are finally obtained as: A positive value of U(t k ) shows a positive trend, whereas a negative trend appears if the U(t k ) value is negative. The null hypothesis is rejected if is the critical value of the typified normal distribution with a probability higher than α 2 . For a significance level α = 5%, the critical value of Z (1− α 2 ) is 1.9604. The values of U'(t k ) are also obtained from regressive temporal data series, t k' , E[t k' ], Var[t k' ] by using Equations (9)-(11), and finally: By plotting the curves of the U(t) and U'(t) as a function of the year, a break point is detected where the curves cross and diverge. (5) Mann-Whitney U test This test can detect a step change in a time series by assessing the difference in the means of two sub-series arising from splitting the complete data series. The original time series (y i , i = 1, 2, . . . , n) is broken into two subseries, one from y 1 to y n1 and the other from y n1+1 to y n , of sizes n 1 and n 2 (n − n1 + 1), respectively. A new data series, z t (t = 1, . . . , n) is obtained by rearranging the original time series y i in increasing order of magnitude. The Mann-Whitney test statistic is obtained by: where R(y i ) is the rank of observation y i in the series z t . If |Z c | > Z (1− α 2 ) the hypothesis of equal means in both subseries is rejected, being Z (1− α 2 ) the (1 − α/2) quantile of the normal distribution. (6) Cumulative Sum test The distribution-free cumulative sum technique was proposed by [37] to detect if a significant step change occurred in a given time series at a significance level α = 0.05 [38]. The null hypothesis considers that no step change exists. For the given data series y i (i = 1, . . . , n), the test statistics V k is obtained as: where k = 1, 2, . . . , n; y median is the median of the time series, and the sign function is: A step-change point consists of any year when the maximal or minimal V k falls outside the 95% confidence limit of ±1.36

Multifractal Characterization
Processes that look the same regardless of the scale at which they are observed are considered to be fractal type. Natural fractal processes, such as rainfall [26] have a statistical nature and their scaling properties can be expressed by statistical relationships [39,40]. Usually, the fractal analysis of data sets is performed by means of fixed-size algorithms that rely on the subdivision of the system into smaller parts with equal size, changing that size iteratively [41].
The scaling exponents of the qth moments of the process under study can be analyzed thorough the Generalized Fractal Dimension (GFD) function, D q , also known as the Rénji spectrum [42] by: where Z(q, r) is the partition function, depending on the values of the mass distribution (µ) of the points that represent a set data composed of n data. Considering a box of length r to fill the set, the number of ni points of data included can be obtained and also the value of µ i = n i n . The partition function is then given as Z(q, r) = n i=1 µ q i . Considering the limit of Equation (16) when q→1, the value of D1 is obtained as: The value of D q for q = 0 is known as the box-counting dimension or the fractal dimension of the system. For q = 1 the information entropy is obtained. The so-called correlation dimension is obtained for q = 2.
The log-log plot of Z(q, r) versus r gives the value of the linear segments of the mass exponent τ q . The relation between the Rénji spectrum (D q ) and τ q is obtained as: The Lipschitz-Hölder exponent, α q = − dτ q dq , obtains the multifractal spectrum f (α) as: The multifractal spectrum shows an inverted parabolic shape and its parameters can be related to interesting characteristics of the data set [25]. The smoothest and most extreme events in the process are related to the α max and α min values, respectively. The structure of the process can be explained through the value of α 0 . Processes with low correlation and fine structure show high values of α 0 .
Low values of α 0 are expected for regular processes. The asymmetry parameter is related to the role that extreme events play in the temporal structure of the data set. The width of the spectrum measures the degree of multifractality of the processes.

Break Point Detection
The BR and SNH tests were applied to all the extreme annual daily rainfall data series for each station selected in the Umbria region. The results obtained for both tests are shown in Table 2, where the statistics and their critical values are given for a significance level of 0.05. The results of the BR test show that only at the Spoleto station, the value of Q is higher than the critical one Q c , and an inhomogeneity is detected in 1977. The SNH test results show that no inhomogeneity is present in any of the data series analyzed, considering that all the statistics T 0 are below their critical values T 0c .
The four non-parametric tests described in the Methodology section were also applied to the data series. The results obtained for the PT test are shown in Table 3. For a significance level of α = 0.05, the value of K is higher than the critical value K c only for the Spoleto extreme data series. A break point is then detected in Spoleto in 1969. For the rest of the data series analyzed in the Umbria region, no inhomogeneities were detected by the PT test. The SQMK test was applied to all the extreme annual data series. Inhomogeneities were found for the Spoleto and Terni stations, and the results are shown in Figure 2. For Spoleto, a decreasing trend is found after 1965 according to the behavior of the U(t k ) curve. The progressive U(t k ) and regressive U'(t k ) curves cross and diverge in 1980 where a break point is detected; the average extreme annual rainfall for the 1959 to 1980 period is 70.65 mm while for the period from 1981 to 2017 it is 60.21 mm.
After 1981, the trend is decreasing and is significant starting from 2003 when the U(t k ) values become lower than the confidence level limit considered (α = 0.05). For the Terni station, a break point is detected in 1957 according to the behavior of progressive and regressive SQMK test curves. A decreasing trend in the extreme annual rainfall data values starts in hat year, and is statistically significant for short periods in the early seventies (from 1969 to 1974) and for the years 1994, 1995 and 1996.
The step changes detected by the MWU test for the analyzed series are shown in Figure 3. Considering the differences in the means, the years with break points are 1968, 1994, 1969 and 1962, for the Nocera Umbra, Petrelle, Spoleto and Terni stations, respectively. The average values before the break point are higher than those after it for all sites except for Petrelle, where the average value of extreme rainfall data is higher after the inhomogeneity date. For the Terni station, a break point is detected in 1957 according to the behavior of progressive and regressive SQMK test curves. A decreasing trend in the extreme annual rainfall data values starts in hat year, and is statistically significant for short periods in the early seventies (from 1969 to 1974) and for the years 1994, 1995 and 1996.
The step changes detected by the MWU test for the analyzed series are shown in Figure 3. Considering the differences in the means, the years with break points are 1968, 1994, 1969 and 1962, for the Nocera Umbra, Petrelle, Spoleto and Terni stations, respectively. The average values before the break point are higher than those after it for all sites except for Petrelle, where the average value of extreme rainfall data is higher after the inhomogeneity date.
Finally, the application of CUSUM technique to the available data series detected no step-changes in any years. As an example, the results for the Bastia Umbra and Terni stations are shown in Figure 4. The V k values obtained for both stations data series are always below the confidence limit considered, being close for year 1969 in Bastia Umbra and for 1960 and 1962 in the Terni station.
In order to summarize the results obtained with all of the applied tests,  Finally, the application of CUSUM technique to the available data series detected no stepchanges in any years. As an example, the results for the Bastia Umbra and Terni stations are shown in Figure 4. The Vk values obtained for both stations data series are always below the confidence limit considered, being close for year 1969 in Bastia Umbra and for 1960 and 1962 in the Terni station.
In order to summarize the results obtained with all of the applied tests, Table 4 shows the year of the break points detected at each station. For Bastia Umbra, Bevagna, Ponte Nuovo and Todi, no inhomogeneities were detected by any of the tests applied. For the Nocera and Petrelle stations, only the MW test was able to detect a break point in 1968 and 1994, respectively. Two tests indicated inhomogeneities for the rainfall data at Terni station: the SQMK test in 1957 and MWU test in 1962. Three different years were detected as break points in the Spoleto rainfall data: 1969 by the PT test, 1977 by BR test, and 1980 by the SQMK test. Table 4. Summary of the years in which break points were detected by the applied tests.

Station Parametric Tests Non-Parametric Tests BR SNH PT SQMK MWU CUSUM
Step changes detected by the Mann-Whitney U test at the Nocera Umbra, Spoleto, Petrelle and Terni stations.  According to the results shown above, the multifractal behavior of the available data series will be used to analyze whether a break point detected can be considered as real or not, when only one test is able to detect it or when different years with inhomogenities arise.  Table 4. Summary of the years in which break points were detected by the applied tests.

Station Parametric Tests Non-Parametric Tests
According to the results shown above, the multifractal behavior of the available data series will be used to analyze whether a break point detected can be considered as real or not, when only one test is able to detect it or when different years with inhomogenities arise.

Multifractal Behavior
For all stations at which a break point was detected (Nocera Umbra, Petrelle, Spoleto and Terni), multifractal characterization of extreme annual 24-h rainfall was performed for the complete data sets, and also for the series obtained from splitting the original data sets for the year of the inhomogeneity. For all of the data series, the GFD function was obtained and the main parameters are shown in Table 5. For all situations, the value of the fractal dimension D 0 is equal to 1. This value is related to the number of boxes needed to cover the fractal object or system under analysis. For all the sites and series, the value of D 0 = 1 fills the entire 1D domain.
The value of the fractal dimension for q = 1, D 1 , describes the degree of heterogeneity of the measure and characterizes the distribution and intensity of singularities with respect to the mean [43,44]. For the complete extreme rainfall data sets, the most heterogeneous is the one at the Spoleto station, followed by the Petrelle, Nocera Umbra and Terni stations. This behavior can also be analyzed through the value of (D 0 -D 1 ). The greater the value, the less uniform is the data set. The most uniform data series are those at Nocera Umbra and Terni stations, where fewer differences in the data than those in Petrelle and Spoleto are expected.
The correlation dimension values, D2, are also shown in Table 5. They can be used to analyze how often recurrence can be expected in a system, so thus high values are characteristic of less recurrent and predictable data series [45]. The highest value is obtained for Spoleto, followed by Petrelle, Terni and Nocera Umbra, respectively. On the contrary, high values of (D0-D2) ( Table 5) are expected for predictable data series.
Even though the fractal dimension values are the main important to be considered in order to analyze the fractal behavior of a system, the trend of D q for several q values (Rénji spectrum) is shown in Figure 5. The multifractal spectra f(α) were also obtained for all the available data series ( Figure 6) and several multifractal parameters are shown in Table 6. Table 5. Main results obtained for the Rénji spectrum of the extreme rainfall data series for which break points were detected by the statistical tests.

Station
Complete Data Series Before Break Point Series After Break Point Series    . Multifractal spectra f(α) for the complete rainfall data series and also for the rainfall data series before and after the break points detected in the Nocera Umbra, Petrelle, Spoleto and Terni stations.

Discussion
The same information can be analyzed by focusing on the similarities and differences between the fractal dimension values for complete data series and the series before and after the break points. For the Nocera Umbra station, the values obtained for the data set before the break point are similar to those obtained for the complete data set. More differences appear for the data set after the inhomogeneity. According to (D0-D1), more uniform data are found before the break point than after it. Moreover, the former data set is also more predictable than the latter, according to the (D0-D2) values.
The similarities found in Nocera Umbra between the complete and one of the split data sets are also found for the Terni station, where the break point year detected is 1957. Almost equal values of D1, D2, (D0-D1) and (D0-D2) are found, showing similar uniformity and degree of predictability in the data sets.
Different values are obtained for the Petrelle data sets, with the complete data series being less uniform than the one before the break point according to the (D0-D1) values, and more predictable than the one after the break point according to the (D0-D2) values. Figure 6. Multifractal spectra f(α) for the complete rainfall data series and also for the rainfall data series before and after the break points detected in the Nocera Umbra, Petrelle, Spoleto and Terni stations.

Discussion
The same information can be analyzed by focusing on the similarities and differences between the fractal dimension values for complete data series and the series before and after the break points. For the Nocera Umbra station, the values obtained for the data set before the break point are similar to those obtained for the complete data set. More differences appear for the data set after the inhomogeneity. According to (D 0 -D 1 ), more uniform data are found before the break point than after it. Moreover, the former data set is also more predictable than the latter, according to the (D 0 -D 2 ) values.
The similarities found in Nocera Umbra between the complete and one of the split data sets are also found for the Terni station, where the break point year detected is 1957. Almost equal values of D 1 , D 2 , (D 0 -D 1 ) and (D 0 -D 2 ) are found, showing similar uniformity and degree of predictability in the data sets.
Different values are obtained for the Petrelle data sets, with the complete data series being less uniform than the one before the break point according to the (D 0 -D 1 ) values, and more predictable than the one after the break point according to the (D 0 -D 2 ) values.
Similar values for the fractal dimension values D 1 , D 2 , and their differences with D 0, are found for the complete data set of Spoleto, and those obtained with the extreme rainfall data set after 1980.
Completely different values are found for both data series obtained before and after the break point in 1969, and those obtained for the complete data set.
For the Nocera Umbra station, the GFD function D q for the complete data series is almost coincident with the one obtained for the rainfall data series before the break point. Different values were obtained for the rainfall data series after the break point. The similarity in the behavior of the complete data series and the one before the potential break point means that no real break point exists since the addition of the data after it does not change the general pattern.
The Rénji spectra obtained for the Petrelle station data series are different. In this case, the three data series should be considered as different and thus the break point is real.
No similarities arise between the data series before or after the break points with respect to the complete data series at Spoleto station if all the D q values are considered. Nevertheless, for negative q values, similar D q functions are obtained for the data sets before and after the break points in 1969 and 1980. These similitudes are maintained for positive q values for D q functions of data sets before both break points. Nevertheless, D q functions are different for positive q values for data sets after the break points, with the one after the inhomogeneity at 1980 being the closest to the D q function for the complete data series. This behavior means that the break point took place in 1969.
For the Terni station, similar trends of D q function are found for the complete data series and those before and after the break point year of 1957, respectively. For the former series, similar values are found for q positive values whereas for the latter the similarities arise for negative q values. Different D q functions are obtained for those series derived from the break point in1962. This behavior implies that the inhomogeneity is placed at 1962 in the Terni extreme rainfall data series.
Regarding the shape of the multifractal spectra, similar f(α) are obtained for the Nocera Umbra complete data set and for the one before the break point year. Even though left-side skewed multifractal spectra are found for all the data series, the one for the data set after the break point differs from those previously mentioned.
Different multifractal spectra are also obtained for the Petrelle data sets, where not even the symmetry is similar; the complete data set spectrum is slightly right-skewed, whereas those obtained before and after the break point are skewed in the opposite direction.
Multifractal spectra before the break points appear similar at Spoleto, where the one for the complete data set also differs from those obtained before and after the break points, showing some similarities only for the left side of f(α) after the inhomogeneity in 1980.
Similarities in the multifractal spectra arise for the Terni station for both the complete data set and those before and after the break point in 1957, whereas different results are found for f(α) before and after the break point in 1962 and the whole rainfall data set.
For Nocera Umbra station, the complete extreme rainfall data series shows a negative value for the asymmetry parameter, which means that extreme events are important in the data series, and also there is a high complexity according to the value of parameter w. Even though the asymmetry parameter values are also negative for data sets before and after the break point, the last data series seems to be less influenced by extreme events than the previous one (or even than the complete series). The richness of the process is also lower for the data set after the break point than for both complete and before break point data sets. Almost similar values for all the parameters (α min , α max , a s and w) were obtained when analyzing both the complete data set and the one before the break point.
For the Petrelle station, different values are obtained for all the parameters and data sets. The influence of extreme events is present in both data sets before and after the break point according to the negative value of the asymmetry parameter (a s ), whereas the complete data series shows a greater influence of less extreme events, with a positive value of a s . The width of the spectra (w) is also different with a loss of richness (lower multifractality) in the complete data set compared to the one before the break point, maybe due to the influence of the data after the break point, where the richness has the lowest value.
For the Terni station, all the spectra exhibit negative values for the asymmetry parameter (a s ), which implies the dominance of extreme events in all the data series. The values closer to the αmin and αmax parameters are found for the complete data series and for the series before and after 1957, and thus a similar influence by extreme and smooth events can be expected for all the series. Similar complex processes are found according to the values of w. This behavior is not present for data series before and after the break point of 1962, where a loss of multifractality is found after the break point compared to the data series before it. This last data series also shows the strong influence of extreme events according to both the values of a s and αmin.
The multifractal spectrum for the Spoleto complete data series shows a negative value for the asymmetry parameter, which is influenced by the presence of extreme events in the data set. The shape changes to a positive value of the asymmetry parameter for the data sets before the break points. The most similar values of all the parameters are found for the complete data set and the one after the break point of 1980.
Authors such us [25] used the variability in the values of α 0 and w from the multifractal spectrum to assess changes in climate data series. For all the data sets studied in this work, the values of α 0 are always equal to 1 and thus no possible changes can be explained based on these values. Nevertheless, there are differences in the values obtained for the parameter that describes the degree of multifractality of the data series, w. Thus, the differences in the values of this multifractal parameter will be of great importance in deciding if a break point can really be identified in a data series.
According to all the results described above, it can be stated that no inhomogeneity exists in the Nocera extreme annual rainfall data series, whereas the break point detected in the Petrelle data set by the statistical test can be confirmed. For the Spoleto station the behavior is not as clear as in the previous stations if all the parameters are considered, but the closest similarity between complete and split data series is found in 1980, specially for the width of the multifractal spectrum (w) and the fractal dimensions D 1 and D 2 . Thus, the most important inhomogeneity is the one detected by the statistical tests in 1969. For the Terni station, the greatest differences in multifractal properties are found between the complete data series and those obtained considering the inhomogeneity in 1962, especially for the value of w, D 1 and D 2 . Therefore, no break point should be considered in 1957.

Conclusions
The multifractal character of rainfall data has been widely studied for extreme annual data of different durations and it is useful when dealing with rainfall predictive modeling or IDF estimation. Thus, the possibility of confirming the existence of a break point previously detected by statistical methods, based on the multifractal behavior of rainfall can help to enhance the results of methods and models that use rainfall as a direct or indirect input, especially, those that depend on the length of the rainfall data series.
The results of this work show that the multifractal characterization of extreme rainfall annual data series can be a suitable tool to decide whether or not an inhomogeneity exists in the data set, especially when there is no agreement in the results from statistical tests. Parametric and non-parametric statistical tests commonly applied to detect shifts in time data series were used to check if any break points were present in extreme annual 24-h rainfall data series in the Umbria region (Italy). For some of the stations considered, no inhomogeneities were found with the statistical tests.
For those stations where only one of the tests detected a break point, or where different years for the break points were identified by the different tests, the maintenance of multifractal parameters' values was used as the criterion to reject the potential break point. Because multifractal analysis results are sensitive to climate shifts if the spectra parameters of time series divided into subsets differ, the values of the heterogeneity fractal dimension D 1 , the correlation fractal dimension D 2 , and the complexity of the data set given by w provide the bases for the decision criteria. If the complete data series maintained the values of D 1 , D 2 or w with respect to any of the split series before or after the break point, it was considered that no inhomogeneity existed. The multifractal behavior of the split series was also found in the complete series. On the contrary, if the values of the multifractal parameters were different among the data series, the break point under consideration was accepted.
According to the results obtained in this work, the multifractal behavior of extreme rainfall data can be considered as an indicator of the changes in the dynamics of the atmospheric processes before and after observed shifts. This is of great importance in models based on extreme rainfall frequency analysis where the existence of a break point in a data series changes the length of it and can influence the estimated values. To make it easier for practitioners and engineers, a software that combine homogeneity tests and multifractal algorithms is being developed.