Temporal Variability of Temperature Extremes in the Sardinia Region (Italy)

: In this paper, the temporal tendencies of temperature data from the island of Sardinia (Italy) were analyzed by considering 48 data series in the period 1982–2011. In particular, monthly temperatures (maximum and minimum), and some indices of daily extremes were evaluated and tested to detect trends using the Mann-Kendall non-parametric test. Results showed a positive trend in the spring months and a marked negative trend in the autumn-winter months for minimum temperatures. As regards maximum temperatures, almost all months showed positive trends, although an opposite behavior was detected in September and in the winter months. With respect to the extreme indices, a general increasing trend of the series was detected for the diurnal temperature range (DTR), frost days (FD), summer days (SU25), warm (WSDI) and cold (CSDI) spells. As regards tropical nights (TR20), an equal distribution of positive and negative trends has emerged. Results of the spatial analysis performed on the trend marks suggested that Sardinia’s topography could inﬂuence temperature variability. selected 48 stations on a study-area digital elevation model. Author Contributions: Conceptualization, T.C. and I.G.; methodology, T.C.; software, T.C.; validation, T.C. and I.G.; formal analysis, T.C.; investigation, I.G.; resources, T.C. and I.G.; data curation, I.G.; writing—original draft preparation, T.C. and I.G.; writing—review and editing, T.C. and I.G.; visualization, T.C. and I.G.; supervision, T.C.


Introduction
In its latest special report, the IPCC evidenced a temperature rise of approximately 1.0 • C above pre-industrial levels, caused by human activities, and forecast that it will likely reach 1.5 • C between 2030 and 2052 [1]. Due to global warming, over the last decades, temperature behavior has been the focus of several studies, which confirmed the mean temperature increase at different spatial scales [2][3][4][5]. Even though trend assessment of mean temperature is significant, analysis of the extreme temperature tendencies, particularly at small spatial scales is paramount, given its impacts on several natural processes and human activities, such as hydrology, food security assessment, crops, agro-ecological zoning, as well as human health and comfort [6,7]. In fact, temperature increase can affect hydrological processes, speeding up the circulation of water vapor and influencing the spatiotemporal distribution and intensity of precipitation [8]. It also directly influences hydrological features such as evaporation, runoff and soil water, which could easily raise the number of extreme climate events. Therefore, the occurrence frequency and strength of droughts and floods will increase, and the contradictory relationship between water supply and demand will become more serious [9]. Furthermore, when high temperature values rise above critical levels, a rise in human death rates can be observed [10]. For example, a severe heat wave affected Europe during the summer of 2003, with considerable health, ecological, societal, and economic impacts [11]. Among these, a high death rate among the elderly, documented in several countries was undoubtedly the most significant [12]. The ecological impacts involved wildfires, increased contamination, a decrease in livestock, low quality crops, and a reduction of forest areas, flora and fauna [13][14][15][16]. The economic damages in 2003 have been evaluated to be more than 10 billion of US dollars [11]. Several climate follows Sicily (about 26,000 km 2 ) and precedes Cyprus (10,000 km 2 ) in the list of the largest islands in the Mediterranean Sea. The orography of the island is mainly characterized by plains, especially in the northwestern and in the southern areas, but also by alluvial valleys and mountains [54], with an average and a maximum altitude of 337 m a.s.l. and 1834 m a.s.l., respectively (Figure 1). Due to its position, the Köppen-Geiger classification [55] identifies the climate of the island as a hot-summer Mediterranean climate, thus presenting relatively mild winters (with rain) and very hot summers (often very dry).
Hydrology 2020, 7, x FOR PEER REVIEW 3 of 18 islands in the Mediterranean Sea. The orography of the island is mainly characterized by plains, especially in the northwestern and in the southern areas, but also by alluvial valleys and mountains [54], with an average and a maximum altitude of 337 m a.s.l. and 1834 m a.s.l., respectively ( Figure  1). Due to its position, the Köppen-Geiger classification [55] identifies the climate of the island as a hot-summer Mediterranean climate, thus presenting relatively mild winters (with rain) and very hot summers (often very dry). In the present study, complete or near-complete records of daily temperature data have been downloaded from the website of the Sardinia Region (http://www.regione.sardegna.it) for the period 1922-2011 [56]. The database consists of high-quality rainfall and temperature data which have been extensively used in past climatological studies performed on the region, e.g., [57][58][59]. In particular, daily temperature data collected at 191 stations are available online for the period 1982-2011, with an average density of one station per 125 km 2 . These data have been collected following the same methodology across the different stations, and the temperature values used in the paper are the absolute maximum and minimum values registered in a day from automatic stations. In order to In the present study, complete or near-complete records of daily temperature data have been downloaded from the website of the Sardinia Region (http://www.regione.sardegna.it) for the period 1922-2011 [56]. The database consists of high-quality rainfall and temperature data which have been extensively used in past climatological studies performed on the region, e.g., [57][58][59]. In particular, daily temperature data collected at 191 stations are available online for the period 1982-2011, with an average density of one station per 125 km 2 . These data have been collected following the same methodology across the different stations, and the temperature values used in the paper are the absolute maximum and minimum values registered in a day from automatic stations. In order to perform a reliable statistical analysis in this work, temperature series with more than 10% of missing data were discarded. Thus, the analyses have been performed on 48 temperature series (Figure 1), with an average density of one station per 490 km 2 . Table 1; Table 2 show the main statistics of the minimum and maximum temperature series, respectively.

Methodology
The aim of this study is to analyze the temporal characteristics of extreme temperatures in Sardinia. To achieve this, a set of well-known indices, summarizing frequency, intensity and persistence of daily temperatures were extracted from the list proposed by the Expert Team of the World Meteorological Organization and Climate Variability and Predictability (ETCCDMI), e.g., [60,61]. In particular, 14 indices have been selected (Table 3) representing the annual daily maximum and minimum temperatures extremes (TXx, TNx, TXn, TNn, all abbreviations defined in Table 3), their mean difference (DTR), the annual number of days in which temperatures exceeded a fixed threshold (FD, TR20, SU25) or percentiles (TN10p, TX10p, TN90p, TX90p, WSDI, CSDI). As a result, these indices can be grouped into three categories: warm, cold and variability indices [62]. Specifically, due to the different sampled temperature distribution, which can be detected over large areas by using indices based on fixed thresholds, the percentile thresholds-based indices are particularly important and can be considered more appropriate than the former for a spatial comparison of the extremes [63]. For example, at mid-latitude climates, a useful indicator for cases of extreme cold could be the number of days evidencing minimum temperature below 0 • C (FD) which, at extreme latitudes, represents normality during winter nights. Additionally, in mild climate settings, where the average summer maxima are around 18 • C, the number of summer days above 25 • C (SU25) could be a useful index to indicate abnormally warm conditions. In order to detect the existence of temporal tendencies in the monthly temperatures and in the temperature extremes, the non-parametric test of Mann-Kendall (MK) [52,53] was applied considering the statistical significance of the trend at the 90%, 95% and 99% confidence intervals. For a series with dimension n, the MK statistic is obtained as: in which x i and x j are the values of the variable at times i and j, respectively, with i < j.
In particular, if the x i values are independent and randomly ordered, and for a dimension n > 10, the statistic S can be represented by a normal distribution presenting zero mean and variance given by: with t i ties number with dimension i.
Finally, the test statistic Z MK is standardized as: Using a two-tailed test for a specified significance level α, the null hypothesis is accepted if |Z MK | is lower than Z 1-α/2 , otherwise the null hypothesis is rejected and the trend can be considered significant. Specifically, in April, 58.3%, 39.6% and 20.8% of the series showed a positive trend with a significance level (SL) of 90%, 95% and 99% respectively. On the contrary, a negative trend has been detected only in 4.2% (SL = 90%) and 2.1% (SL = 95%) of the series. No significant negative trends have been identified for an SL = 99%. In May, positive (negative) trends have been detected for 40.4% (8.5%), 34.0% (2.1%) and 14.9% (no trend) of the series with a significance level (SL) of 90%, 95% and 99% respectively. In June, the percentage of series showing a positive trend is similar to the one identified in May (41.7%, 33.3% and 14.6% for the ascending SL) while a negative trend has been detected in 12.5% (SL = 90%), 8.3% (SL = 95%) and 2.1% (SL = 90%) of the series.

Analysis of the Monthly Minimum Temperature
The spatial distribution of the trends (SL = 95%) shows that in April, May and June the temperature series presenting increasing values are distributed across the region, while the negative trends are observed in the southern stations, especially in June, and in the eastern mountain of the region ( Figure 3).
An opposite trend behavior than in April, May and June has been detected in the other months, and in particular in the autumn-winter period ( Figure 2). In fact, in February, in October and particularly in September, a marked negative trend has been identified respectively in 34.0%, 39.6% and 42.6% of the series, with a significance level of 95%, and respectively in 14.9%, 18.8% and 23.4% of the series, with a significance level of 99%. In the same months a positive trend has been detected for an SL = 90% only in 6.4%, 6.3% and 4.3% of the gauges. From the spatial analysis, the results obtained in these months (SL = 95%) evidenced a positive trend distributed across the region and in particular in the south-eastern side of the island, and scattered negative trend in the central area ( Figure 3). In addition to autumn and winter, a negative trend, albeit less marked, has also been detected in March, July and August. In effect, in these months, there are few differences among the percentage of series showing a positive or a negative trend ( Figure 2). In particular, in March the percentage of series showing a negative trend with an SL = 90% (28.3%) and an SL = 95% (19.6%) is slightly higher than the one presenting a positive trend, equal to 19.0% and 13.0%, respectively. On the contrary, for an SL = 99%, a higher percentage of positive trend series has been detected, 6.5% against 4.3% of negative trend series. The spatial distribution of the trend (SL = 95%) in these months did not show any peculiarity, although the positive trend seems to be mainly localized in the central area of the region (Figure 3). Hydrology 2020, 7, x FOR PEER REVIEW 7 of 18 The spatial distribution of the trends (SL = 95%) shows that in April, May and June the temperature series presenting increasing values are distributed across the region, while the negative trends are observed in the southern stations, especially in June, and in the eastern mountain of the region (Figure 3).
An opposite trend behavior than in April, May and June has been detected in the other months, and in particular in the autumn-winter period ( Figure 2). In fact, in February, in October and particularly in September, a marked negative trend has been identified respectively in 34.0%, 39.6% and 42.6% of the series, with a significance level of 95%, and respectively in 14.9%, 18.8% and 23.4% of the series, with a significance level of 99%. In the same months a positive trend has been detected for an SL = 90% only in 6.4%, 6.3% and 4.3% of the gauges. From the spatial analysis, the results

Analysis of the Monthly Maximum Temperature
As regards maximum temperatures, a marked positive trend has been identified for all the monthly values, with the exception of January, September, December and, partially, October ( Figure  4).

Analysis of the Monthly Maximum Temperature
As regards maximum temperatures, a marked positive trend has been identified for all the monthly values, with the exception of January, September, December and, partially, October (Figure 4). The highest numbers of stations with positive trends have been identified in the spring and summer months (from March to August), with a maximum percentage in June. During this month 62.5% 56.3% and 35.4% of the stations evidenced this tendency with an SL = 90%, 95% and 99%, respectively, and a negative trend has been detected in 6.3% (SL = 90%), 6.3% (SL = 95%) and 4.2% (SL = 90%) of the series. Relevant positive trends have also been detected in April and May, when significant trends have been identified respectively in 56.3%, and 50.0% of the series, with a significance level of 95%, and respectively in 33.3% and 16.7% of the series, with a significance level of 99%. In the same months, only in 6.3% of the gauges has a positive trend has been detected for a SL = 90% (Figure 4). The spatial distribution of the trend (SL = 95%) in these months is similar to the one obtained for the minimum temperature. In fact, the temperature series presenting increasing values are distributed across the region, while the negative trends are observed around the eastern mountains of the region ( Figure 5). In February, October and November, there are few differences among the percentage of series showing positive or negative trends, although the positive trend prevails, especially for lower SLs (Figure 4). In fact, the percentage of series showing a negative trend with a SL = 99% in February (6.3%), October (6.3%) and November (6.7%) is equal or slightly higher than the one presenting a positive trend, equal to 2.1%, 4.2% and 6.7%, respectively. From the spatial analysis, the results obtained in these months (SL = 95%) evidenced a positive trend distributed in the southern side of the region, and scattered negative trends in the central area ( Figure 5). Finally, in January, September and December a prevailing negative trend has been identified respectively in 21.3%, 25.5% and 21.3% of the series, with a significance level of 95%, and in 10.6%, 8.5% and 17.0% of the series respectively, with a significance level of 99%. In the same months, only in 6.4%, 12.8% and 10.6% of the gauges a positive trend has been detected for an SL = 90% (Figure 4). The spatial analysis of the results obtained in these months (SL = 95%) evidenced relevant results only in January and December, with the positive trend affecting the eastern side of the island ( Figure 5).   The highest numbers of stations with positive trends have been identified in the spring and summer months (from March to August), with a maximum percentage in June. During this month 62.5% 56.3% and 35.4% of the stations evidenced this tendency with an SL = 90%, 95% and 99%, respectively, and a negative trend has been detected in 6.3% (SL = 90%), 6.3% (SL = 95%) and 4.2% (SL = 90%) of the series. Relevant positive trends have also been detected in April and May, when significant trends have been identified respectively in 56.3%, and 50.0% of the series, with a significance level of 95%, and respectively in 33.3% and 16.7% of the series, with a significance level of 99%. In the same months, only in 6.3% of the gauges has a positive trend has been detected for a SL = 90% (Figure 4). The spatial distribution of the trend (SL = 95%) in these months is similar to the one obtained for the minimum temperature. In fact, the temperature series presenting increasing values are distributed across the region, while the negative trends are observed around the eastern mountains of the region ( Figure 5). In February, October and November, there are few differences among the percentage of series showing positive or negative trends, although the positive trend prevails, especially for lower SLs (Figure 4). In fact, the percentage of series showing a negative trend with a SL = 99% in February (6.3%), October (6.3%) and November (6.7%) is equal or slightly higher than the one presenting a positive trend, equal to 2.1%, 4.2% and 6.7%, respectively. From the spatial analysis, the results obtained in these months (SL = 95%) evidenced a positive trend distributed in the southern side of the region, and scattered negative trends in the central area ( Figure 5). Finally, in January, September and December a prevailing negative trend has been identified respectively in 21.3%, 25.5% and 21.3% of the series, with a significance level of 95%, and in 10.6%, 8.5% and 17.0% of the series respectively, with a significance level of 99%. In the same months, only in 6.4%, 12.8% and 10.6% of the gauges a positive trend has been detected for an SL = 90% (Figure 4). The spatial analysis of the results obtained in these months (SL = 95%) evidenced relevant results only in January and December, with the positive trend affecting the eastern side of the island ( Figure 5).
Hydrology 2020, 7, x FOR PEER REVIEW 10 of 18 Figure 5. Spatial distribution of monthly maximum temperature trend.

Analysis of the Extreme Temperature Indices
Given the minimum and the maximum daily temperature values recorded in Sardinia in the period 1922-2011, some extreme temperature indices have been analyzed for trend detection. With respect to TXx and TNx, which are the maximum extreme values of daily temperature, more than

Analysis of the Extreme Temperature Indices
Given the minimum and the maximum daily temperature values recorded in Sardinia in the period 1922-2011, some extreme temperature indices have been analyzed for trend detection. With respect to TXx and TNx, which are the maximum extreme values of daily temperature, more than 50% of the series evidenced significant trends (SL = 90%). In fact, considering the TXx index, 41.7% (20.8%), 31.3% (12.5%) and 12.5% (8.3%) of the gauges evidenced a positive (negative) trend, with an SL = 90%, 95% and 99%, respectively ( Figure 6).
An opposite trend behavior has been observed for TNx with 31.9% (SL = 90%), 23.4% (SL = 95%) and 4.3% (SL = 99%) of the series presenting a negative trend, and 19.1%, 17.0% and 10.6% of the series evidencing a positive trend for an SL equal to 90%, 95% and 99%, respectively ( Figure 6). Spatially, the temperature series of TXx and TNx presenting increasing values are distributed across the region, while the negative trends are mainly concentrated around the eastern mountains and in the north-eastern side of the region (Figure 7). As regards the annual extremes minimum temperatures (TXn and TNn), an opposite trend behavior has been detected, with the temperature series showing a prevailing positive trend of TXn and a negative trend of TNn, although for both TXn and TNn the majority of the series, 70.2% and 66.7%, respectively, showed an absence of trend. In particular, for TXn, 23.4%, 12.8% and 2.1% of the series showed a positive trend with SLs of 90%, 95% and 99%, respectively. On the contrary, a negative trend has been detected only in 6.4% (SL = 90% and SL = 95%) of the series ( Figure 6). For TNn 22.9% (SL = 90%), 16.7% (SL = 95%) and 2.1% (SL = 99%) of the series showed a negative trend and only 10.4% showed a positive trend for an SL = 90% ( Figure 6). The spatial analysis of the trend results, for an SL = 95%, evidenced that for both TXn and TNn the significant trends are localized in the southern part of the region (Figure 7). The trend analysis performed on the extreme percentile-based indices (TN10p, TN90p, TX10p and TX90p) evidenced prevailing positive trends for TX10p, TN10p and TX90p ( Figure 6). In fact, for the TX10p index, 27.1%, 18.8% and 10.4% of the series showed a positive trend with SLs of 90%, 95% and 99%, respectively, mainly localized around the eastern mountain and in the northeastern side of the region (Figure 8).
Hydrology 2020, 7, x FOR PEER REVIEW 12 of 18 99%) of the series showed a negative trend and only 10.4% showed a positive trend for an SL = 90% ( Figure 6). The spatial analysis of the trend results, for an SL = 95%, evidenced that for both TXn and TNn the significant trends are localized in the southern part of the region (Figure 7). The trend analysis performed on the extreme percentile-based indices (TN10p, TN90p, TX10p and TX90p) evidenced prevailing positive trends for TX10p, TN10p and TX90p ( Figure 6). In fact, for the TX10p index, 27.1%, 18.8% and 10.4% of the series showed a positive trend with SLs of 90%, 95% and 99%, respectively, mainly localized around the eastern mountain and in the northeastern side of the region ( Figure 8).  Table 3 for abbreviation definitions).
The tropical nights (TR20) show an almost equal number of positive and negative trends ( Figure  6), with the latter mainly located in the southern part of the region (Figure 9). In particular, the percentage of negative trends ranges from 16.7% (SL = 99%) to 29.2% (SL = 90%) while the percentage of positive trends ranges from 12.5% (SL = 99%) to 22.9% (SL = 90%).
Finally, both the warm (WSDI) and the cold (CSDI) spell duration indicators denote an increasing trend involving 22.9% of the stations for both the WSDI and the CSDI with an SL = 95% ( Figure 6). The majority of these trends are localized in the southern part of the region (Figure 9).

Discussion
Recently, a number of extreme climate indices, whether using observed, reanalysis, or climate projection data, have become available in the literature; however, small-scale analysis based on highdensity networks are paramount to perform reliable analysis. As an example, in Italy, Zollo et al. [64] detected a general underestimation of both minimum and maximum temperature extremes for simulations driven by the CMCC-CM global model. They also found that for indices depending on maximum temperature, the spatial pattern is less skillfully represented in the case of regions  Table 3 for abbreviation definitions).
On the contrary, a negative trend was detected in 25.0% (SL = 90%), 18.8% (SL = 95%) and 8.3% (SL = 99%) of the series especially in the southern part of the region (Figure 8). A more marked difference between positive and negative trends has been observed for TN10p and TX90p. In particular, for TN10p, 35.4% (18.8%), 29.2% (12.5%) and 10.4% (4.2%) of the series showed a positive (negative) trend with SLs of 90%, 95% and 99%, respectively ( Figure 6). As regards TX90p, 58.3% (SL = 90%), 41.7% (SL = 95%) and 22.9% (SL = 99%) of the series showed a negative trend while only a maximum of 16.7% of the series showed a positive trend for an SL = 90% ( Figure 6). The spatial analysis of the trend results for both TN10p and TX90p evidenced positive trends distributed across the region with few negative values in the central areas of the region (Figure 8). On the other hand, for the TN90p, a decreasing trend has been observed in more than 40% of the series, for SL = 90% and 95%, and in about 23% of the series for an SL = 99% (Figure 6), but also in this case no clear spatial patterns are visible in the region (Figure 8).
The analysis of the DTR index evidenced a marked increasing trend with 54.2% (SL = 90%), 47.9% (SL = 95%) and 41.7% (SL = 99%) of the series presenting positive values, and 22.9%, 10.4% and 6.3% of the series evidencing a negative trend for an SL equal to 90%, 95% and 99%, respectively ( Figure 6). The spatial analysis of the trend results showed an almost uniform distribution of the positive trend across the region with few negative values in the central areas of the region ( Figure 9). As regards the frost days (FD), 26.1%, 19.6% and 6.5% of the series showed a positive trend with an SL of 90%, 95% and 99%, respectively (Figure 6), mainly localized around the eastern mountain and in the southwestern side of the region (Figure 9). Conversely, a negative trend has been detected in 13.0% (SL = 90%) and 6.5% (SL = 95%). With respect to the summer days (SU25), a positive trend in 38.3% (SL = 90%), 27.7% (SL = 95%) and 17.0% (SL = 99%) and a negative trend in 23.4% (SL = 90%), 14.9% (SL = 95%) and 8.5% (SL = 99%) of the stations have been observed ( Figure 6). Spatially, this evidence appears to be distributed all over the region without a clear difference between areas with positive or negative values (Figure 9). minimum temperatures at a planetary scale evidenced an increase in the minimum values rather than the maximum ones [66], results of this study evidenced an opposite behavior, thus confirming the findings achieved by previous research performed in Italy [41]. Similarly, these results only partially conform to the rising tendency identified since 1976 in warm days and nights in some European areas, and disagree with the weakness of the decreasing trends in cold days and nights [67]. Relating to diurnal temperature range (DTR), as verified in previous surveys, an important growing trend was detected due to a faster warming in maximum temperatures than in minimum ones, although some other DTR research works at a global scale proved an opposite behavior, showing a large decreasing trend [68][69][70][71]. In Sardinia, increasing frequency and intensity of warmer extremes, together with a weak increasing trend of the colder extremes, occurs. These outcomes agree with previous studies on daily temperature extremes in Mediterranean countries [42,48,72,73]. Figure 9. Spatial distribution of extreme temperature indices trend (DTR, FD, SU25, TR20, WSDI and CSDI; see Table 3 for abbreviation definitions).
In accordance with the warming in global mean temperatures [1], a decrease in cold extremes and an increase in warm ones has occurred since the middle of the twentieth century in a vast majority of land regions. In Europe, this trend has been already experienced since 1960 [74], when a statistically significant increase in the number of warm days and nights, and an opposite trend in the number of cold days and nights, have been detected. However, there is no agreement in the literature on trend behavior, and differing tendencies for temperature extremes have been recognized over different time-periods during the last century [14,75].

Conclusions
The temporal variability assessment of monthly and extreme temperatures in the Sardinia region (Italy) was carried out to identify the trends, and to analyze the behaviors of the extreme Figure 9. Spatial distribution of extreme temperature indices trend (DTR, FD, SU25, TR20, WSDI and CSDI; see Table 3 for abbreviation definitions).
The tropical nights (TR20) show an almost equal number of positive and negative trends ( Figure 6), with the latter mainly located in the southern part of the region (Figure 9). In particular, the percentage of negative trends ranges from 16.7% (SL = 99%) to 29.2% (SL = 90%) while the percentage of positive trends ranges from 12.5% (SL = 99%) to 22.9% (SL = 90%).
Finally, both the warm (WSDI) and the cold (CSDI) spell duration indicators denote an increasing trend involving 22.9% of the stations for both the WSDI and the CSDI with an SL = 95% ( Figure 6). The majority of these trends are localized in the southern part of the region (Figure 9).

Discussion
Recently, a number of extreme climate indices, whether using observed, reanalysis, or climate projection data, have become available in the literature; however, small-scale analysis based on high-density networks are paramount to perform reliable analysis. As an example, in Italy, Zollo et al. [64] detected a general underestimation of both minimum and maximum temperature extremes for simulations driven by the CMCC-CM global model. They also found that for indices depending on maximum temperature, the spatial pattern is less skillfully represented in the case of regions characterized by a marked orography such as Calabria and Sardinia. Moreover, as evidenced by Kyselý and Plavcová [65] the E-OBS v2.0 gridded dataset shows large biases when compared to a high-density network. Within this context, although previous investigations on maximum and minimum temperatures at a planetary scale evidenced an increase in the minimum values rather than the maximum ones [66], results of this study evidenced an opposite behavior, thus confirming the findings achieved by previous research performed in Italy [41]. Similarly, these results only partially conform to the rising tendency identified since 1976 in warm days and nights in some European areas, and disagree with the weakness of the decreasing trends in cold days and nights [67]. Relating to diurnal temperature range (DTR), as verified in previous surveys, an important growing trend was detected due to a faster warming in maximum temperatures than in minimum ones, although some other DTR research works at a global scale proved an opposite behavior, showing a large decreasing trend [68][69][70][71]. In Sardinia, increasing frequency and intensity of warmer extremes, together with a weak increasing trend of the colder extremes, occurs. These outcomes agree with previous studies on daily temperature extremes in Mediterranean countries [42,48,72,73].
In accordance with the warming in global mean temperatures [1], a decrease in cold extremes and an increase in warm ones has occurred since the middle of the twentieth century in a vast majority of land regions. In Europe, this trend has been already experienced since 1960 [74], when a statistically significant increase in the number of warm days and nights, and an opposite trend in the number of cold days and nights, have been detected. However, there is no agreement in the literature on trend behavior, and differing tendencies for temperature extremes have been recognized over different time-periods during the last century [14,75].

Conclusions
The temporal variability assessment of monthly and extreme temperatures in the Sardinia region (Italy) was carried out to identify the trends, and to analyze the behaviors of the extreme temperatures, in the period 1922-2011, by means of 14 commonly used daily temperature indices. From this analysis we draw some important concluding remarks: (a) For the minimum temperature series, diffuse negative tendencies have been identified, with a positive trend observed only in the spring months; (b) For the maximum temperature series, a general increasing trend has been identified in almost all the months, with the exclusion of January, September, December and, partially, October. In particular, a consistent positive trend in the spring and in the summer months has been detected, especially in June; (c) Among the evaluation of the indices of extreme temperatures, the DTR results evidenced a marked increasing trend distributed across the region. The frost days (FD) showed a positive trend of the series mainly localized around the eastern mountain and in the southwestern side of the region. The summer days (SU25) displayed more positive trends than negative ones, spatially distributed all over the region. The tropical nights (TR20) showed an equal distribution of positive and negative trends. Finally, both the warm (WSDI) and the cold (CSDI) spell indices displayed an increasing trend localized in the southern part of the region. (d) Results indicate that Sardinia's topography may cause temperature variability. In fact, the mountains crossing the island constitute a barrier between the tropical air masses coming from the African coasts on eastern side, and air masses carried by western winds originating from the Atlantic Ocean on the western side. The changes in the temperature regime in this region could have regionally specific impacts on its ecosystems. (e) These results can be very useful for several stakeholders. In particular: (i) planners, disaster management agencies and policy makers who need to be alert in a region with a high-density population; (ii) hydrologists and climatologists who need to investigate the potential causes of the observed temporal trends to delineate the variability of extremes due to natural and anthropogenic effects.