SPI Trend Analysis of New Zealand Applying the ITA Technique

A natural temporary imbalance of water availability, consisting of persistent lower-than-average or higher-than-average precipitation, can cause extreme dry and wet conditions that adversely impact agricultural yields, water resources, infrastructure, and human systems. In this study, dry and wet periods in New Zealand were expressed using the Standardized Precipitation Index (SPI). First, both the short term (3 and 6 months) and the long term (12 and 24 months) SPI were estimated, and then, possible trends in the SPI values were detected by means of a new graphical technique, the Innovative Trend Analysis (ITA), which allows the trend identification of the low, medium, and high values of a series. Results show that, in every area currently subject to drought, an increase in this phenomenon can be expected. Specifically, the results of this paper highlight that agricultural regions on the eastern side of the South Island, as well as the north-eastern regions of the North Island, are the most consistently vulnerable areas. In fact, in these regions, the trend analysis mainly showed a general reduction in all the values of the SPI: that is, a tendency toward heavier droughts and weaker wet periods.


Introduction
Recently, the adverse impacts of climate change have become the focus of considerable international attention due to the increase in phenomena such as flood, heat waves, forest fires, and droughts [1,2].Among these damaging climate events, drought phenomena play a significant role in socio-economic and health terms, even though their impact on populations depends on the vulnerable elements [3].Moreover, understanding drought phenomena is paramount for the appropriate planning and management of water resources [4].For example, different drought events have been detected during the last decades [5][6][7], and drought is expected to become more frequent in the 21st century in some seasons and areas [8] following precipitation and/or evapotranspiration variability [9].
In recent years, several researchers have analyzed drought events in several parts of the world [10][11][12][13][14][15][16][17], even though drought phenomena are difficult to detect and to monitor due to their complex nature.Usually, drought severity is evaluated by means of drought indices since they facilitate communication of climate anomalies to diverse user audiences; they also allow scientists to assess quantitatively climate anomalies in terms of their intensity, duration, frequency, recurrence probability, and spatial extent [3,18].In the past few decades, numerous indices were proposed for identifying and monitoring drought events.Some of these indices refer to meteorological drought (scarcity of precipitation) and are based on the analysis of the rainfall information only.Thus, different categories of drought can be investigated by choosing appropriate temporal scales addressing different categories of users.Other indices, however, are more suitable to describe hydrological drought (scarcity in surface and subsurface water supplies), agricultural drought (water shortage compared to the typical needs crops irrigation), and socio-economic drought (referred to the global water consumption).Specifically, meteorological drought consists of temporary lower-than-average precipitation and results in diminished water resources availability [19], which impact on economic activities, human lives, and the environment [20].The most well-known index for analyzing the meteorological drought is undoubtedly the Standardized Precipitation Index (SPI) proposed by McKee et al. [21], which has been extensively applied in different countries [22][23][24][25][26][27][28].This drought index can be considered one of the most robust and effective drought indices, as it can be evaluated for different time scales and allows the analysis of different drought categories [29].Moreover, the evaluation of the SPI requires only precipitation data, making it easier to calculate more than complex indices, and allows for the comparison of drought conditions in different regions and for different time periods [30][31][32][33].Due to its intrinsic probabilistic nature, the SPI is the ideal candidate for carrying out drought risk analysis [34,35].With this aim, several authors focused on the SPI trend [36][37][38].These studies are mainly based on non-parametric tests, which are better suited to deal with non-normally distributed hydrometeorology data than the parametric methods.Recently, Şen [39] proposed the Innovative Trend Analysis (ITA) technique, which allows a graphical trend evaluation of the low, medium, and high values in the data.The ITA technique was widely applied to the trend detection of several hydrological variables.Haktanir and Citakoglu [40] analyzed the annual maximum rainfall series by means of the ITA method.Kisi and Ay [41] studied some water quality parameters registered at five Turkish stations by means of the ITA and the MK.Şen [42] and Ay and Kisi [43] applied the ITA to Turkish temperature data.The ITA technique was also used to analyze the trends of heat waves [44], monthly pan evaporations [45], and streamflow data [46].
Since agriculture is one of the largest sectors of the tradable economy, a period of drought in New Zealand can have significant ecological, social, and economic impacts [47].In fact, New Zealand experiences rainfall deficits and short duration of dry spells that are not as unusual as isolated drought events at regional level.For example, the widespread drought event that affected New Zealand from late 2007 to the end of autumn 2008 caused damages of about 2.8 billion New Zealand dollars [48].The 2013 drought in New Zealand was estimated to have caused GDP (Gross Domestic Product) to fall by 0.6% [49].Regional scenarios of drought in New Zealand evidenced an increase in drought trends during this century in all the areas presently subject to drought [50].Furthermore, based on the latest climate and impact modelling, more droughts can be expected in the future in some locations such as the agricultural regions on the Eastern coast and particularly the Canterbury Plains, as well as Northland [51].
In this article, drought events in several regions of New Zealand have been studied by applying the SPI at various time scales (3,6,12, and 24 months) starting from a database of 294 monthly rainfall series in the period 1951-2010.In particular, this work aims to identify the most drought-prone regions of New Zealand by analyzing its evolution through the identification of the SPI trend at different timescales by means of the Innovative Trend Analysis (ITA), which allows the trend identification of the low, medium, and high values of a series.

Standardized Precipitation Index
In this study, dry and wet periods were evaluated using the SPI at different time scale (3,6,12 and 24 months).In fact, while the 3-and 6-month SPI describe droughts that affect plant life and farming, the 12-and 24-month SPI influence the way how water supplies/reserves are managed [52,53].Angelidis et al. [54] offered a meticulous description of the method to compute the SPI.
In order to calculate the index, for each time scale, an appropriate probability density function (PDF) must be fitted to the frequency distribution of the cumulated precipitation.In particular, a gamma function is considered.The shape and the scale parameters must be estimated for each month of the year and for each time aggregation, for example by using the approximation of Thom [55].
Since the gamma distribution is undefined for a rainfall amount x = 0, in order to take into account the zero values that occur in a sample set, a modified cumulative distribution function (CDF) must be considered.
with G(x) the CDF and q the probability of zero precipitation, given by the ratio between the number of zero in the rainfall series (m) and the number of observations (n).Finally, the CDF is changed into the standard normal distribution by using, for example, the approximate conversion provided by Abramowitz and Stegun [56].
with c 0 , c 1 , c 2 , d 1 , d 2 , and d 3 as mathematical constants.Table 1 reports the climatic classification according to the SPI provided by the National Drought Mitigation Center.This index is now habitually used in the classification of wet periods, even though the original classification provided by McKee et al. [21] was limited to drought periods only.

Innovative Trend Analysis (ITA)
The ITA method has been first proposed by Şen [39]: unlike the MK test or other methods, the ITA's greatest advantage is the fact that it does not require any assumptions (serial correlation, non-normality, sample number, and so on).First, the time series is divided into two equal parts, which are separately sorted in ascending order.Then, the first and the second half of the time series are located on the X-axis and on the Y-axis, respectively.In Figure 1, a graphical representation of the innovative method on a Cartesian coordinate system is shown.If the data are collected on the 1:1 ideal line (45 • line), there is no trend in the time series.If data are located on the upper triangular area of the ideal line, an increasing trend in the time series exists, while if data are accumulated in the lower triangular area of the 1:1 line, there is a decreasing trend in the time series [39,42].Thus, trends of low, medium, and high values of any hydro-meteorological or hydro-climatic time series can be clearly identified through this method.
For example, the results shown in Figure 1 identify a positive trend of the lowest values and a negative trend of the highest ones, while no trends can be detect for the medium values that lies near the 1:1 ideal line.
In this work, in order to easily and better identify the possible trend of the severe dry and wet conditions, two vertical bands have been added in Figure 1: a red band corresponding to the severe drought limit (SPI = −1.5)and a blue band corresponding to the severe wet limit (SPI = 1.5).

Study Area
New Zealand, with an elongated shape of 1930 km and a maximum width of 400 km and a surface area of about 270,000 km 2 , is located in the southern hemisphere 2500 km east from Australia.Its altitude distribution is mainly mountainous; in fact, 76% of its territory presents an elevation higher than 600 m a.s.l. with a maximum altitude of 3724 m (Figure 2).

Study Area
New Zealand, with an elongated shape of 1930 km and a maximum width of 400 km and a surface area of about 270,000 km 2 , is located in the southern hemisphere 2500 km east from Australia.Its altitude distribution is mainly mountainous; in fact, 76% of its territory presents an elevation higher than 600 m a.s.l. with a maximum altitude of 3724 m (Figure 2).

Study Area
New Zealand, with an elongated shape of 1930 km and a maximum width of 400 km and a surface area of about 270,000 km 2 , is located in the southern hemisphere 2500 km east from Australia.Its altitude distribution is mainly mountainous; in fact, 76% of its territory presents an elevation higher than 600 m a.s.l. with a maximum altitude of 3724 m (Figure 2).New Zealand's climate is mainly influenced by the following physical factors [57]: a.The water masses surrounding the country, which results in cool summers and moderately cold winters, and the cold air masses from Antarctica, which cause snow and frost in some areas of the country; b.The mountains crossing both the islands, which constitute a barrier against air flows coming in from south-south west, cause significant differences in the rainfall amounts, also within short distances; c.The proximity of the Australia's eastern land/sea boundary, which is characterized by a low-pressure region of cyclonic circulation toward the Tasman Sea.
For a more thorough description of the New Zealand's climate, interested readers can refer to the National Institute of Water and Atmosphere Research (NIWA) of New Zealand, which published a detailed analysis of New Zealand's climate [58].
In order to carry out an SPI analysis at different time scales, the New Zealand National Climate database of the National Institute of Water and Atmospheric Research (NIWA) has been selected.Several studies on New Zealand's climate [59][60][61][62][63][64][65] made extensive use of this database, which presents high-quality data and complete, or near complete, records for the period 1951-2010.Until 2010, the New Zealand National Climate database consisted of measurements collected at 3011 stations, with a density of one station per 89 km 2 .In particular, for the present analysis, the monthly rainfall data were extracted and used after performing record error checks and metadata analyses for inhomogeneities detection.Following these procedures, a number of station series, which presented either low quality records or few years of observation for statistical purposes (<50 available years of data), were discarded.The series ending before 2010 and those showing more than 5% of lacking data were also discarded.Thus, the final selection included 294 series longer than 50 years with a density of one station per 913 km 2 (Figure 2).

Results and Discussion
Following the NIWA, which provides climate maps at regional scale, for every region of New Zealand shown in Figure 2, an average SPI series has been evaluated for each time scale, and a trend analysis has been performed through the application of the ITA approach.In particular, the SPI has been evaluated for each rain gauge and for each time scale and then a simple arithmetic average of the obtained SPI values has been computed for each region.
Before applying the ITA, the number of months showing severe or extreme dry and wet conditions was evaluated.As a result, different conditions have been detected considering the different time scales.
As to what concerns the 3-month SPI (Figure 3a) in nine out of 14 regions the number of months showing severe or extreme wet conditions is higher than the ones showing severe or extreme dry conditions.In particular, in the South Island, in the Canterbury and Southland areas, only in 15 and 18 months dry conditions have been detected while, instead, 32 and 27 months evidenced wet conditions, respectively.Differently from these regions, in the Hawke's Bay area (North Island), the months that showed dry conditions (35 months) clearly outperform the months (24) in which wet conditions have been detected.
The 6-month SPI (Figure 3b) showed a different behavior with respect to the 3-month SPI.In fact, in eight regions, the number of months showing severe or extreme dry conditions is higher than the one presenting wet conditions.Dry conditions have been detected especially in the North Island and, specifically, in the Auckland (27 months against 16), Bay of Plenty (36 months against 29), East Cape (37 months against 31), Wanganui-Manawatu (28 months against 13), and Hawke's Bay (31 months against 19) areas.In the South Island, only the Nelson-Marlborough region showed dry conditions in 26 months (17 wet months), while in the other regions, wet conditions have been evidenced.
Concerning the 12-month SPI (Figure 3c) an equal number of regions showed prevailing dry or wet conditions.In particular, in the North Island, in six out of nine regions, the number of months showing dry conditions is higher than the ones showing wet conditions.Conversely, in the South Island, four regions showed wet conditions, and only one region (Nelson-Marlborough) evidenced dry conditions.Specifically, in the North Island, relevant results have been obtained in the Bay of Plenty and in the East Cape regions, with 38 and 40 months showing severe or extreme dry conditions, respectively, while, in the South Island, the Westland (45 months against 24) and the Southland (38 months against 22) regions marked dry conditions have been detected.
Finally, the 24-month SPI (Figure 3d) showed a clear difference between the two islands.In fact, in the North Island, severe or extreme dry conditions have been detected in all the regions, while severe or extreme wet conditions have been detected in all the regions of the South Island.In particular, the driest conditions in the North Island have been detected in the Taranaki (73 months against 34) and in the Bay of Plenty (61 months against 37) regions, while, in the South Island, the Southland (73 months against 34) and the Canterbury (61 months against 37) areas showed the highest number of months with wet conditions.
Geosciences 2018, 8, x FOR PEER REVIEW 6 of 14 showing dry conditions is higher than the ones showing wet conditions.Conversely, in the South Island, four regions showed wet conditions, and only one region (Nelson-Marlborough) evidenced dry conditions.Specifically, in the North Island, relevant results have been obtained in the Bay of Plenty and in the East Cape regions, with 38 and 40 months showing severe or extreme dry conditions, respectively, while, in the South Island, the Westland (45 months against 24) and the Southland (38 months against 22) regions marked dry conditions have been detected.Finally, the 24-month SPI (Figure 3d) showed a clear difference between the two islands.In fact, in the North Island, severe or extreme dry conditions have been detected in all the regions, while severe or extreme wet conditions have been detected in all the regions of the South Island.In particular, the driest conditions in the North Island have been detected in the Taranaki (73 months against 34) and in the Bay of Plenty (61 months against 37) regions, while, in the South Island, the Southland (73 months against 34) and the Canterbury (61 months against 37) areas showed the highest number of months with wet conditions.With the aim to detect possible trends in the 3-, 6-, 12-, and 24-month SPI values, for each region the ITA method was applied to the monthly series of the index.The ITA method allowed to evidence the tendencies of both low and high SPI values, thus including values referred to wet conditions.As a result of the ITA approach, Figures 4-7 show the results obtained at regional level for the 3-, 6-, 12- With the aim to detect possible trends in the 3-, 6-, 12-, and 24-month SPI values, for each region the ITA method was applied to the monthly series of the index.The ITA method allowed to evidence the tendencies of both low and high SPI values, thus including values referred to wet conditions.As a result of the ITA approach, Figures 4-7 show the results obtained at regional level for the 3-, 6-, 12-and 24-month SPI, respectively.All the SPI series were divided into two 30-year sub-series: from 1951 to 1980, and from 1981 to 2010.
Generally, the main result obtained for the 3-month SPI values was a negative trend of the highest values of the index, which is related to weaker wet periods (Figure 4).This tendency has been detected in nine out of 14 regions but with a different behavior of the lowest SPI values.In fact, in four regions of the North Island (Northland, Auckland, Bay of Planty and East Cape) and in the Canterbury region in the South Island, a negative trend of both the lowest and the highest values of the index has been detected, thus evidencing heavier droughts and weaker wet periods.At the same time, in the Waikato, Wellington, Southland, and Westland regions, a positive trend of the lowest values and negative trend of the highest ones has been evidenced, both indicating weaker droughts and wet periods.Differently from the previous regions, the Otago, Wanganui-Manawatu, and Taranaki areas evidenced a tendency through weaker droughts and heavier wet periods given by a positive trend of both the lowest and the highest SPI values.Generally, the main result obtained for the 3-month SPI values was a negative trend of the highest values of the index, which is related to weaker wet periods (Figure 4).This tendency has been detected in nine out of 14 regions but with a different behavior of the lowest SPI values.In fact, in four regions of the North Island (Northland, Auckland, Bay of Planty and East Cape) and in the Canterbury region in the South Island, a negative trend of both the lowest and the highest values of the index has been detected, thus evidencing heavier droughts and weaker wet periods.At the same time, in the Waikato, Wellington, Southland, and Westland regions, a positive trend of the lowest values and negative trend of the highest ones has been evidenced, both indicating weaker droughts and wet periods.Differently from the previous regions, the Otago, Wanganui-Manawatu, and Taranaki areas evidenced a tendency through weaker droughts and heavier wet periods given by a positive trend of both the lowest and the highest SPI values.Finally, in the Nelson-Marlborough region, a negative trend of the lowest values (heavier droughts) and a positive trend of the highest values (heavier wet periods) have been detected.The results of the ITA methods on the Hawke's Bay region did not show a clear tendency, with the lowest and highest values falling close to the no trend line.As regards the 6-month SPI, results confirm the ones obtained for the 3-month SPI, with a spreading negative trend of the highest values of the index (Figure 5).In fact, similar results to the 3-month SPI have been obtained in Northland, Auckland, Bay of Planty, Wanganui-Manawatu, Taranaki, Nelson-Marlborough, Westland, and Southland.Different from the results obtained for the 3-month SPI, in the Waikato and in the Hawke's Bay regions, a negative trend of both the lowest and the highest values of the SPI index has been detected.Moreover, in the Canterbury and Otago region, a positive trend of the lowest values (weaker droughts) and a negative trend of the highest values (weaker wet periods) have been evidenced.Finally, the Wellington and East Cape regions did not show a clear tendency.
As regards the 6-month SPI, results confirm the ones obtained for the 3-month SPI, with a spreading negative trend of the highest values of the index (Figure 5).In fact, similar results to the 3month SPI have been obtained in Northland, Auckland, Bay of Planty, Wanganui-Manawatu, Taranaki, Nelson-Marlborough, Westland, and Southland.Different from the results obtained for the 3-month SPI, in the Waikato and in the Hawke's Bay regions, a negative trend of both the lowest and the highest values of the SPI index has been detected.Moreover, in the Canterbury and Otago region, a positive trend of the lowest values (weaker droughts) and a negative trend of the highest values (weaker wet periods) have been evidenced.Finally, the Wellington and East Cape regions did not show a clear tendency.Studies on the 3-and 6-month SPI trends, which impact vegetation and agricultural practices, are paramount for New Zealand because agriculture is one of the largest sectors of the economy.As a summary of the results of the trend analysis on these time scales, in the South Island, which is generally characterized by wet conditions, a negative trend of the highest SPI value has been detected in the majority of the regions, thus indicating a tendency toward weaker wet periods.In particular, in the agricultural area of the Canterbury, the 3-month SPI evidenced a decreasing trend of both the lowest and the highest values of the index, which lead to heavier droughts and weaker wet periods.In the North Island, where the majority of the regions showed a higher percentage of dry months Studies on the 3-and 6-month SPI trends, which impact vegetation and agricultural practices, are paramount for New Zealand because agriculture is one of the largest sectors of the economy.As a summary of the results of the trend analysis on these time scales, in the South Island, which is generally characterized by wet conditions, a negative trend of the highest SPI value has been detected in the majority of the regions, thus indicating a tendency toward weaker wet periods.In particular, in the agricultural area of the Canterbury, the 3-month SPI evidenced a decreasing trend of both the lowest and the highest values of the index, which lead to heavier droughts and weaker wet periods.In the North Island, where the majority of the regions showed a higher percentage of dry months than wet months, a general negative trend of both the lowest and the highest SPI values has been identified, with the exception of some region such as Wanganui-Manawatu and Taranaki.
Considering the long time scales, the 12-month SPI trend showed similar results than the 6-month SPI but with an increase in the negative trend (Figure 6).In fact, in the Northland, Waikato, Auckland, Bay of Planty, East Cape, Hawke's Bay, and Canterbury regions, a reduction of all the SPI values has been detected, thus confirming a clear tendency toward heavier droughts and weaker wet periods.On the contrary, the Wanganui-Manawatu, Nelson-Marlborough, and Otago regions showed trend results that indicate weaker droughts and heavier wet periods.An increase of the lowest values and a decrease in the highest values (weaker wet and dry periods) has been evidenced in the Wellington, Westland, and Southland regions.In the Taranaki area, no clear tendencies have been detected in the severe and extreme SPI values.
than wet months, a general negative trend of both the lowest and the highest SPI values has been identified, with the exception of some region such as Wanganui-Manawatu and Taranaki.
Considering the long time scales, the 12-month SPI trend showed similar results than the 6month SPI but with an increase in the negative trend (Figure 6).In fact, in the Northland, Waikato, Auckland, Bay of Planty, East Cape, Hawke's Bay, and Canterbury regions, a reduction of all the SPI values has been detected, thus confirming a clear tendency toward heavier droughts and weaker wet periods.On the contrary, the Wanganui-Manawatu, Nelson-Marlborough, and Otago regions showed trend results that indicate weaker droughts and heavier wet periods.An increase of the lowest values and a decrease in the highest values (weaker wet and dry periods) has been evidenced in the Wellington, Westland, and Southland regions.In the Taranaki area, no clear tendencies have been detected in the severe and extreme SPI values.As regards the 24-month SPI, in eight regions (Northland, Waikato, Auckland, Bay of Planty, East Cape, Wanganui-Manawatu, Hawke's Bay, and Canterbury), the results of the trend analysis with the ITA method evidenced a clear reduction of the lowest SPI values.Among these regions, the same result has been obtained for the highest values, with the exception of the Waikato and Auckland, in which the highest values showed a positive trend.By contrast, the Taranaki, Nelson-Marlborough, Otago Westland, and Southland regions showed an increase in the lowest values, with concomitant decrease in the highest values in the Taranaki and Otago areas and increase in the As regards the 24-month SPI, in eight regions (Northland, Waikato, Auckland, Bay of Planty, East Cape, Wanganui-Manawatu, Hawke's Bay, and Canterbury), the results of the trend analysis with the ITA method evidenced a clear reduction of the lowest SPI values.Among these regions, the same result has been obtained for the highest values, with the exception of the Waikato and Auckland, in which the highest values showed a positive trend.By contrast, the Taranaki, Nelson-Marlborough, Otago Westland, and Southland regions showed an increase in the lowest values, with concomitant decrease in the highest values in the Taranaki and Otago areas and increase in the highest values in the Nelson-Marlborough, Westland, and Southland regions.As also evidenced for the 3-month SPI, the Wellington region did not show a clear tendency.
The 12-and 24-month SPI are a broad proxy for water resource management.Results of the trend analysis on these time scales confirm the ones obtained on a short time scale but with an increase in the number of regions where a negative trend has been detected.In fact, current global level assessments suggest that droughts are expected to both increase and decrease following future climate change depending upon geographic location [66].Based on the latest climate and impact modelling, New Zealand can expect more droughts in the future in some locations [51].
As a result, this work has evidenced increase in drought trend in all the areas that are presently subject to drought, supporting what has been evidenced in past studies [50,64].In fact, the results of this paper confirm the geographic pattern of change found by Mullan et al. [50] and Caloiero [64], which mainly detected a drought increase in the future projections on the East Coast and no change in drought projections for the West Coast of the South Island.Specifically, as also evidenced by Clark et al. [51] and Caloiero [64], the results of this paper highlight that key agricultural regions on the Eastern side such as the Canterbury Plains are the most consistently vulnerable areas in the South Island, together with other regions in the North Island, including a key primary industry region like Waikato.highest values in the Nelson-Marlborough, Westland, and Southland regions.As also evidenced for the 3-month SPI, the Wellington region did not show a clear tendency.The 12-and 24-month SPI are a broad proxy for water resource management.Results of the trend analysis on these time scales confirm the ones obtained on a short time scale but with an increase in the number of regions where a negative trend has been detected.In fact, current global level assessments suggest that droughts are expected to both increase and decrease following future climate change depending upon geographic location [66].Based on the latest climate and impact modelling, New Zealand can expect more droughts in the future in some locations [51].
As a result, this work has evidenced an increase in drought trend in all the areas that are presently subject to drought, supporting what has been evidenced in past studies [50,64].In fact, the results of this paper confirm the geographic pattern of change found by Mullan et al. [50] and Caloiero [64], which mainly detected a drought increase in the future projections on the East Coast and no change in drought projections for the West Coast of the South Island.Specifically, as also evidenced by Clark et al. [51] and Caloiero [64], the results of this paper highlight that key agricultural regions on the Eastern side such as the Canterbury Plains are the most consistently vulnerable areas in the South Island, together with other regions in the North Island, including a key primary industry region like Waikato.

Conclusions
Generally, the methods used for trend detection allow researchers to analyze general tendencies of climatic variables, but they do not make it possible to identify sub-trends.In this paper, for every region of New Zealand, an average SPI series has first been evaluated at various time scales (3,6,12, and 24 months).Then, for each region, the number of months showing severe extreme dry and wet conditions has been detected.Finally, a graphical technique, based on the Innovative Trend Analysis (ITA) proposed by Şen [39], was applied to the SPI series.As a result, a different behavior emerged between the two islands.In fact, for the several time scale, in the majority of the regions of the North Island, the number of months showing dry conditions is higher than the ones showing wet conditions, while in the South Island, the number of months showing wet conditions is higher than the ones showing dry conditions.In particular, for the 24-month time scale, dry conditions have been detected in all the regions of the North Island, while wet conditions have been detected in the South Island ones.As to what concerns the trend analysis, on a short time scale, the results showed a negative trend of the highest SPI value in the South Island, thus indicating a tendency toward weaker wet periods.In particular, in the agricultural area of the Canterbury, the 3-month SPI evidenced a decreasing trend in both the lowest and the highest values of the index, which lead to heavier droughts and weaker wet periods.In the North Island, a general negative trend of both the lowest and the highest SPI values has been identified, thus evidencing heavier droughts and weaker wet periods.Results of the trend analysis on long time scales confirm the ones obtained on a short time scale but with an increase in the number of regions where a negative trend has been detected, highlighting that agricultural regions on the eastern side of the South Island, as well as the north-eastern regions of the North Island, are the most consistently vulnerable areas.

Figure 2 .
Figure 2. Location of the 14 regions of New Zealand and of the selected 294 rain gauge stations on a Digital Elevation Model (DEM) of New Zealand.

Figure 2 .
Figure 2. Location of the 14 regions of New Zealand and of the selected 294 rain gauge stations on a Digital Elevation Model (DEM) of New Zealand.

Figure 2 .
Figure 2. Location of the 14 regions of New Zealand and of the selected 294 rain gauge stations on a Digital Elevation Model (DEM) of New Zealand.
Finally, in the Nelson-Marlborough region, a negative trend of the lowest values (heavier droughts) and a positive trend of the highest values (heavier wet periods) have been detected.The results of the ITA methods on the Hawke's Bay region did not show a clear tendency, with the lowest and highest values falling close to the no trend line.Geosciences 2018, 8, x FOR PEER REVIEW 7 of 14 and 24-month SPI, respectively.All the SPI series were divided into two 30-year sub-series: from 1951 to 1980, and from 1981 to 2010.

Figure 4 .
Figure 4. Regional results of the Innovative Trend Analysis (ITA) method applied to the 3-month SPI.

Figure 4 .
Figure 4. Regional results of the Innovative Trend Analysis (ITA) method applied to the 3-month SPI.

Figure 5 .
Figure 5. Regional results of the ITA method applied to the 6-month SPI.

Figure 5 .
Figure 5. Regional results of the ITA method applied to the 6-month SPI.

Figure 6 .
Figure 6.Regional results of the ITA method applied to the 12-month SPI.

Figure 6 .
Figure 6.Regional results of the ITA method applied to the 12-month SPI.

Figure 7 .
Figure 7. Regional results of the ITA method applied to the 24-month SPI. Figure 7. Regional results of the ITA method applied to the 24-month SPI.

Figure 7 .
Figure 7. Regional results of the ITA method applied to the 24-month SPI. Figure 7. Regional results of the ITA method applied to the 24-month SPI.

Table 1 .
Climate classification according to the Standardized Precipitation Index (SPI) values.