A Time-Series Analysis of Climate Variability in Urban and Agricultural Sites (Rome, Italy)

Identifying early signals of climate change and latent patterns of meteorological variability requires tools analyzing time series data and multidimensional measures. By focusing on air temperature and precipitation, the present study compares local-scale climate regimes at two sites in Central Italy (urban Rome and a peri-urban cropland 10 km west of Rome), using descriptive and inferential statistics on both variables and a drought index (the Standardized Precipitation Index, hereafter SPI) recorded over the last 60 years (1958–2017). The present work assumes the importance of urban-rural gradients shaping local-scale climate regimes and spatial variability, with differential impacts on individual variables depending on territorial background and intrinsic biophysical characteristics. Considering together precipitations and minimum/maximum air temperature at month and year scale, the analysis developed here illustrates two coexisting climatic trends at distinctive spatial scales: A general trend toward warming—specifically influencing temperature regimes—and a more specific pattern evidencing changes in local-scale climate regime along the urban gradient, with a more subtle impact on both precipitations and temperatures. Empirical results indicate that climate variability increased over the study period, outlining the low predictability of dry spells typical of Mediterranean climate especially in the drier season (spring/summer). On average, absolute annual differences between the two sites amounted to 70 mm (more rainfall in the peri-urban site) and 0.9 ◦C (higher temperature in the urban site). A similar trend toward warming was observed for air temperature in both sites. No significant trends were observed for annual and seasonal rainfalls. SPI long-term trends indicate high variability in dry spells, with more frequent (and severe) drought episodes in urban Rome. Considering together trends in temperature and precipitation, the ‘urban heat’ effect was more evident, indicating a clearer trend toward climate aridity in urban Rome. These findings support the adoption of integrated strategies for climate change adaptation and mitigation in both agricultural systems and relict natural ecosystems surrounding urban areas.


Introduction
Climate regimes showed sequential 'cold' and 'hot' waves in last centuries impacting significantly ecosystems and societies [1]. The average surface temperature grew by more than 0.7 • C in the last century and the time period between 1995 and 2007 was recognized as the warmest since 1850 [2][3][4]. Incorporating time series data available up to now, future climate projections indicate a different scenario from the current averages, foreshadowing a global warming [5] associated with a moderate decrease in rainfalls and locally-severe changes in wind direction and intensity, air pressure and soil humidity [6]. A higher variability of pertinent weather variables was also observed, exacerbating land sensitivity to both droughts and floods [7][8][9][10][11]. Being situated in a transitional region between two climatic systems (tropical and polar), the Mediterranean basin is one of the most sensitive areas to climate change worldwide. Geographical position and territorial factors (e.g., elevation and distance from the sea coast), play a crucial role in shaping climate regimes in the area, determining a particularly high rainfall and temperature variability along the elevation gradient [12]. Climate variability was associated with (i) an increase in average air temperature and (ii) a rise in vegetation stress caused by a mix of progressively higher summer temperatures, lower air humidity, lower rainfall rates and improved extreme events [6,13,14].
Together with climate aridity [15], land-use change and increased human pressure, environmental processes emerged with global warming [16][17][18][19][20][21][22]. For instance, Land Degradation (LD) was increasingly associated with climate change and its effects negatively impact soil quality determining a decrease in land productivity, biodiversity and ecosystem services [23][24][25][26][27]. Therefore, climate variations are related with life-threatening ecological conditions [6,[28][29][30] and require appropriate mitigation and adaptation policies [31]. Considering urbanization as a powerful factor of landscape transformations in Mediterranean Europe, a comparison of changes over time in local-scale climate regimes (e.g., between urban and peri-urban areas) is particularly stimulating in the light of global vs local climate change, drought and land degradation [32]. While climate change is a recognized phenomenon worldwide, an additional factor related to local changes in climate regimes is urban expansion, which generate a well-known effect known as the Urban Heat Island, hereafter UHI [33]. The UHI usually leads to local climate change, with rising temperatures and extreme rainfall events [33][34][35][36]. The UHI consequences involve planning challenges that should increasingly consider the negative environmental impact of climate change and urbanization together [37], mitigating less effectively the rapid increase in temperatures within urban areas in respect with the surrounding areas [38,39]. Earlier studies have focused on inner cities, since UHI stress resulted more pervasive in those places than in peri-urban sites [40,41]. However, a specific investigation of long-term trends in rainfall and temperature comparing urban and neighboring peri-urban sites is still lacking in Mediterranean large cities experiencing rapid urban expansion.
In such contexts, strategies of mitigation and adaptation incorporated in town planning are specifically requested to consider the expected increase in urban heat stress [38,42,43]. However, identifying signals of climate change and increased variability in rainfall and temperature regimes needs improved tools investigating large data sets, possibly based on multidimensional procedures [44]. Defining a homogeneous temporal unit is a prerequisite in the analysis of linear and non-linear long-term trends in temperatures and rainfall, time series break-points and latent relations among meteorological variables. This knowledge forms the necessary base to policy implementation in the field of mitigation and adaptation of climate change in large urban regions.
Based on these premises, our study investigates variability in precipitation and air temperature over a sufficiently long time period  in Rome (Central Italy), comparing trends at an urban site (Collegio Romano, downtown Rome) and at a peri-urban site (Casalotti, a suburb west of Rome). Using descriptive and inferential statistics on both original variables and a drought index (the Standardized Precipitation Index, hereafter SPI), the quantitative analysis developed in this work is aimed at identifying two coexisting trends in regional-scale climate regimes: (i) A general trend toward warming and (ii) a more specific pattern evidencing changes in local-scale climate regime along the urban gradient, and likely involving both temperature and precipitations. Similarly, to other Mediterranean sites [5,13,45], the study area experienced increased climate variations in the last decades. Climate regimes in central Italy, and especially in Rome's metropolitan region, modified relatively rapidly over time and space due to a different influence of environmental conditions and socioeconomic forces at local and regional scales [46]. Results of the present study are particularly relevant also for the intrinsic impact of local-scale climate change on specific economic activities. For instance, agriculture is a traditional activity in peri-urban Rome, recognized and practiced since millennia. Impact of global and local warming on peri-urban agriculture should be better investigated in Mediterranean Europe with the final objective to inform local policies increasing resilience of fringe cropping systems to climate change [47]. In the light of global warming and reinforced urban heat effects, results of this study may contribute to identify sustainable management options for these contexts, improving adaptation to climate change.

Methodology
Two sites were considered in this study: Collegio Romano (downtown Rome, 16 m at the sea level) and Casalotti (western Rome, 80 m at the sea level) ( Figure 1). These contexts were taken as representative of a strictly urban area and a peri-urban area, respectively. Meteorological data were acquired from professional, international standard weather stations operating continuously and producing high-quality data (with no missing values) over the entire study period. The two databases are homogeneous, without missing data, and this has made gap filling procedures unnecessary. Collegio Romano station was one of the historical weather stations in Italy, measuring meteorological variables since the end of the 18th century. Data from both gauging stations were extensively used in official statistics and national/regional reporting of climate conditions prepared by public authorities in charge of public works, environmental management, agriculture and forestry. Casalotti station is placed in a peri-urban context rich in green spaces, with sparse forest patches embedded in the 'Agro Romano' landscape, forming a typical mosaic of farmland, pastures, woods and wetlands. This landscape was increasingly exposed to urban expansion, soil degradation, wildfires and indirect human pressure [48][49][50][51]. The municipal territory of Rome extends 1285 km 2 , being the inner more urbanized part, encircled by a ring-shaped motorway called Grande Raccordo Anulare, with a total length of 345 km 2 .
The urban system extended from about 7% of the municipal area in 1960 to nearly 27% in 2015 [52]. During the first 20 years of the study period, urban expansion was mainly restricted to the inner part of the municipality, followed by peripheral, less compact urbanization. Changes in land-use between 1980 and 2001 show, beside the increment of urban areas, a decrease of extensive agricultural land-use (not irrigated arable land, complex cultivation patterns, sparsely vegetated areas), and the increase of intensive agriculture and forest systems (permanently irrigated arable land, reforestations), as well as semi-natural mixed systems (forests and woodland shrubs) [53]. At the end of the 1990s, about 50% of the municipal area remained cultivated, mainly in peri-urban areas. Figure 1 illustrates aerial photographs of the two sites derived from Google Earth imagery (Figure 1a,b) and is enriched with specific maps of early 1950s (Figure 1c,d), other more recent urbanization derived from Corine Land Cover maps (Figure 1e,f) and a comprehensive land imperviousness assessment carried out by European Environment Agency in the framework of Land Copernicus initiative (Figure 1g     Descriptive statistics were used to evaluate variability over time in air temperatures and precipitations at a year scale. The Standardized Precipitation Index (SPI) is defined as and can be interpreted as the number of standard deviations by which the observed x (in this case the precipitation value) anomaly deviates from the long-term mean. Being used at three different time scales (1 month, 3-months, 6-months), the SPI was further calculated to provide a general overview of drought regimes, comparing the two locations. By computation on cumulated rainfalls over a given time horizon, the SPI is a convenient tool for monitoring meteorological drought, whose assessment is traditionally based on precipitation data. It gives a flexible method to assess drought at a different scale [52,53] extending from near normal (−0.99) to extreme drought conditions (≤−2.0) [54][55][56][57][58][59][60]. SPI values greater than 2 indicate extremely wet conditions, while those less than −2 indicate extreme drought conditions. Using different time windows with the aim to investigate seasonal regimes, total precipitations were recorded. Finally, a composite index proposed by Perini et al. [61] was run to evaluate together precipitation and temperature trends over time. By cumulating the local effect of both variables, this index ranges from 1 (indicating a wet climate regime, with lower temperatures and higher precipitations) to 0 (indicating a drier climate regime, with higher temperatures and lower precipitations). A progressive decrease of the index indicates a trend over time toward climate aridity.

Results
Time series data comparing climate trends in Rome and Casalotti between 1958 and 2017 were analyzed considering together average air temperatures and precipitations at both month and year time scale.

Rainfall
At both sites, precipitation was particularly heterogeneous over time as expected in a Mediterranean dry zone. Precipitations in Casalotti showed high variability over time with a minimum of 443 mm in 2007 and a maximum of 1294 mm in 1976 (Table 1 and Figure 2a). The annual average precipitation was 802 mm with a standard deviation of 188 mm and a coefficient of variation amounting to 23.4% (Table 1) A one-way ANOVA test on the two annual rainfall time series allows to reject the null hypothesis that the means of the two sites are equal (F = 8.25, p = 0.0048) and the observed differences are due to statistical fluctuations. Moreover, the necessary condition for ANOVA test that the two variances are similar is guaranteed (Levene' test, F = 2.62, p = 0.108); this endorses that the observed greater average precipitation in Casalotti is attributable to random fluctuations. A normality test (Kolmogorov-Smirnov) on the two-time series finally allows to accept (p > 0.8 for both sites) the assumption of Gaussian distribution of time series data. Figure 3 exposes the time series of the annual SPI (broken lines) and 5-order polynomial fitting curves (red curves) in Casalotti and Collegio Romano annually (from 1958 to 2017). From the same time period, Figure 4 reports the temporal trend concerning air temperature in Casalotti (a) and Collegio Romano (b). Furthermore, absolute difference between Casalotti and Collegio Romano by year was portrayed ( Figure 4c). The blue line outlines the observed values, defining the trends of air temperature; while, the red line permits to draw the trend line in time, offering the corresponding equation.      By computation on total precipitation over a given time horizon, a Standardized Precipitation Index (SPI) provided a general overview of drought regimes in the area ( Figure 5). Considering a 3months' time window with the aim to investigate seasonal drought regimes, a total precipitation of  By computation on total precipitation over a given time horizon, a Standardized Precipitation Index (SPI) provided a general overview of drought regimes in the area ( Figure 5). Considering a 3months' time window with the aim to investigate seasonal drought regimes, a total precipitation of precipitation observed in Casalotti was 672 mm in 1962 autumn; the minimum precipitation was observed in summer 2000 (0 mm); in Collegio Romano, the maximum cumulated rainfall was observed in autumn 1991 (533 mm) and the minimum in summer 1993 (1.4 mm). Although few severe drought spells were observed in the study period, lower SPI scores were observed in such episodes in Collegio Romano than in Casalotti, especially in winter and autumn, corresponding to climatically wet periods in a typical Mediterranean climate, when the largest amount of rainfall is expected. Two severe drought episodes were observed in autumn (SPI < −2) for Collegio Romano.  Considering a time period of 6 months, a marked gap in cumulated rainfalls was observed between Collegio Romano and Casalotti. Average cumulated precipitations during wet (autumn and winter) and dry seasons (spring and summer) were respectively 546 mm and 255 mm in Casalotti; in Collegio Romano, average cumulated rainfall in autumn and winter was 476 mm and only 235 mm in spring and summer (Figures 6 and 7). Although particularly severe droughts were not observed at this time window, at least 5 and 7 episodes of moderately-severe droughts were observed in the area (SPI < −1.5) respectively during dry and wet seasons (Figure 6a,b). Under time periods with negative SPI, more negative values were systematically observed in Collegio Romano than in Casalotti.
Agriculture 2019, 10, x FOR PEER REVIEW 9 of 18 Considering a time period of 6 months, a marked gap in cumulated rainfalls was observed between Collegio Romano and Casalotti. Average cumulated precipitations during wet (autumn and winter) and dry seasons (spring and summer) were respectively 546 mm and 255 mm in Casalotti; in Collegio Romano, average cumulated rainfall in autumn and winter was 476 mm and only 235 mm in spring and summer (Figure 6 and 7). Although particularly severe droughts were not observed at this time window, at least 5 and 7 episodes of moderately-severe droughts were observed in the area (SPI < −1.5) respectively during dry and wet seasons (Figure 6a, b). Under time periods with negative SPI, more negative values were systematically observed in Collegio Romano than in Casalotti.   Considering a time period of 6 months, a marked gap in cumulated rainfalls was observed between Collegio Romano and Casalotti. Average cumulated precipitations during wet (autumn and winter) and dry seasons (spring and summer) were respectively 546 mm and 255 mm in Casalotti; in Collegio Romano, average cumulated rainfall in autumn and winter was 476 mm and only 235 mm in spring and summer (Figure 6 and 7). Although particularly severe droughts were not observed at this time window, at least 5 and 7 episodes of moderately-severe droughts were observed in the area (SPI < −1.5) respectively during dry and wet seasons (Figure 6a, b). Under time periods with negative SPI, more negative values were systematically observed in Collegio Romano than in Casalotti.   Assuming a time window of 12 months, average cumulated rainfall was higher in Casalotti (802 mm) than in Rome (711 mm) (Figure 7). Moderately severe episodes of drought were concentrated in the last two decades and show a relatively similar intensity in both Casalotti and Collegio Romano.
Analysis of climate variability based on ordinary least square regressions (Figure 2a,b), has been supplemented with a Mann-Kendall test-a non-parametric inference mainly used with non-gaussian data distributions-assessing the statistical relevance of monotonic upward or downward trends over time. Results of this test (Table 2) indicate the absence of any significant trend for both annual and seasonal precipitations along the analyzed time-series. These last results do not exclude a-priori changes in specific climatic variables, especially for extreme events that represent the most attractive topic under the climate change literature [62,63]. Extreme drought and precipitation can cause considerable damage to agriculture, human economies, people's livelihood and property, and even influence the social stability of entire regions [63][64][65]. Earlier studies suggest that the mean of a distribution has very little relationship to extremes and it has been shown that extremes may change to a greater extent than changes in mean [66][67][68]. A 5-order polynomial fitting method applied to the average annual FED (Frequency of Extreme Drought) has been used in order to study differences as well as commonality among different subareas from 1961 to 2015 in a Chinese site [69]. A similar approach is applied here to analyze possible variations in time for the SPI annual values of the two sites. To some extent, the SPI is a sensitive variable reflecting meteorological droughts, soil moisture variations and extreme precipitation events [70,71]. In  Table 3, indicate that, on average, temperature regime in Rome was 0.9 • C-1.0 • C higher than in Casalotti, while temperature increases over time were substantially similar in both sites. Table 3. Summary statistics of annual average temperature at Casalotti and Collegio Romano sites, 1958-2017.

Site
Min A One-way ANOVA test on the two average temperature time series allows to reject the null hypothesis, with a 0.05 level of significance, that the means of the two sites are equal, and that the observed differences are due to statistical fluctuations (F = 53.0, p < 0.001). Moreover, the necessary condition for ANOVA test that the two variances are similar has been confirmed (Levene test, F = 0.36, p = 0.55). A normality test (Kolmogorov-Smirnov) on the two-time series allows to accept (p > 0.57 for both sites) the assumption of Gaussian distribution of data for both locations.
For both sites, average temperature time series exhibited a slight upward trend, as revealed by a significant regression trend (R 2 > 0.4) with a similar positive slope ( 0.026) and intercept ( 15 • C). An overall F-test allows to accept the hypothesis of a statistically significant (at a 0.05 level of significance) positive liner trend for both sites (F > 40, p < 0.001). A composite index, proposed by Perini et al. [72] and evaluating together standardized precipitations and temperatures (1: wet climatic regime; 0: dry climate regime), indicates a trend toward climate aridity over the study period. This trend was significant and particularly intense in Rome (R 2 = 0.34) and less evident in Casalotti (R 2 = 0.16), as clarified with the use of regression coefficients (Figure 8). For both sites, F-test of overall significance provides sufficient evidence (p < 0.001) that the regression model fits the data better than the model with no independent variables. A steep regression slope in Collegio Romano indicates increasing aridity and a milder increase was observed in Casalotti. Taken together, these results indicate a clear trend toward local warming in both sites, with climate aridity being more evident in the urban site in the last decade.
Agriculture 2019, 10, x FOR PEER REVIEW 11 of 18 A One-way ANOVA test on the two average temperature time series allows to reject the null hypothesis, with a 0.05 level of significance, that the means of the two sites are equal, and that the observed differences are due to statistical fluctuations (F = 53.0, p < 0.001). Moreover, the necessary condition for ANOVA test that the two variances are similar has been confirmed (Levene test, F = 0.36, p = 0.55). A normality test (Kolmogorov-Smirnov) on the two-time series allows to accept (p > 0.57 for both sites) the assumption of Gaussian distribution of data for both locations.
For both sites, average temperature time series exhibited a slight upward trend, as revealed by a significant regression trend (R² > 0.4) with a similar positive slope (≅ 0.026) and intercept (≅ 15 °C). An overall F-test allows to accept the hypothesis of a statistically significant (at a 0.05 level of significance) positive liner trend for both sites (F > 40, p < 0.001). A composite index, proposed by Perini et al. [72] and evaluating together standardized precipitations and temperatures (1: wet climatic regime; 0: dry climate regime), indicates a trend toward climate aridity over the study period. This trend was significant and particularly intense in Rome (R 2 = 0.34) and less evident in Casalotti (R 2 = 0.16), as clarified with the use of regression coefficients (Figure 8). For both sites, F-test of overall significance provides sufficient evidence (p < 0.001) that the regression model fits the data better than the model with no independent variables. A steep regression slope in Collegio Romano indicates increasing aridity and a milder increase was observed in Casalotti. Taken together, these results indicate a clear trend toward local warming in both sites, with climate aridity being more evident in the urban site in the last decade.

Discussion
Earlier studies focused on climate change and local-scale variability in last decades, examining mostly thermometric records [61,[73][74][75][76][77][78][79][80], under the assumption that surface air temperature raised by 0.5-0.7 °C since the beginning of the 20th century. Climate analysis over a relatively long-time interval provides a relevant knowledge base to infer the impact of meteorological variables on soil moisture, vegetation growth, crop productivity and ecosystem stability [81], particularly in areas under widespread human pressure [82]. Comparison of local-scale climate regimes between neighboring sites with different land-use and landscape configurations (e.g., along the urban-rural gradient in large metropolitan regions) may provide useful indications to a better understanding of climate variability patterns, distinguishing regional patterns of change (e.g., global warming trends) from modifications in local climate mostly due to rapid urban expansion. The latter process contributed to determine 'urban heat' effects on both temperature and precipitation regimes.
Climate variability resulted to increase over the study period, delineating the low predictability of dry spells typical of Mediterranean climate particularly in the drier season (both spring and summer). The empirical results of this study evidenced (i) local warming (0.7 °C over 60 years) with

Discussion
Earlier studies focused on climate change and local-scale variability in last decades, examining mostly thermometric records [61,[73][74][75][76][77][78][79][80], under the assumption that surface air temperature raised by 0.5-0.7 • C since the beginning of the 20th century. Climate analysis over a relatively long-time interval provides a relevant knowledge base to infer the impact of meteorological variables on soil moisture, vegetation growth, crop productivity and ecosystem stability [81], particularly in areas under widespread human pressure [82]. Comparison of local-scale climate regimes between neighboring sites with different land-use and landscape configurations (e.g., along the urban-rural gradient in large metropolitan regions) may provide useful indications to a better understanding of climate variability patterns, distinguishing regional patterns of change (e.g., global warming trends) from modifications in local climate mostly due to rapid urban expansion. The latter process contributed to determine 'urban heat' effects on both temperature and precipitation regimes.
Climate variability resulted to increase over the study period, delineating the low predictability of dry spells typical of Mediterranean climate particularly in the drier season (both spring and summer).
The empirical results of this study evidenced (i) local warming (0.7 • C over 60 years) with a comparable trend between sites; (ii) site-specific differences in the aridity regime (considering together the joint effect of temperature and precipitations), with drier conditions in downtown Rome than in suburban Rome; (iii) this last statement is supported by a non-trivial, stable trend of the 12-months SPI in downtown Rome, turning toward a drier regime in the last few years. On average, absolute annual differences among the two sites amounted to 70 mm (more rainfall in the peri-urban site) and 0.9 • C (higher temperature in the urban site). A comparable trend toward warming was detected for air temperature in both sites. The persistent gap in average temperatures between sites (0.9 • C higher in Rome than in Casalotti, on average) can be attributable to specific site characteristics, including massive urbanization of inner Rome before the late 1950s, corresponding with the beginning of the observation period. With urban settlements expanding rapidly in Rome between 1960 and 2015 [52], the assumption of a direct influence of urban growth on local warming seems reasonable. Results also exposed that no significant trends were observed for annual and seasonal rainfalls. SPI long-term trends designate high variability in drought spells, with more frequent (and severe) episodes in the urban city of Rome than the countryside's area. Bearing in mind together trends in temperature and precipitation, the 'urban heat' effect was more apparent in Rome, demonstrating a clearer tendency toward climate aridity in the urban city than in peri-urban area. On the contrary, the specific characteristics of the suburban site, which is substantially far away from downtown Rome (approximately 10 km) and surrounded by large green areas mixed with low-density residential settlements, may explain the inherent difference in aridity regimes. A distinctive effect of the UHI on air temperatures could be identified when comparing closer urban and peri-urban sites experiencing more rapid settlement expansion. Conversely, a visible effect of the urban heat island was likely associated with climate aridity rather than with individual rainfall or temperature variables. A significant trend toward climate aridity was observed at both sites, but the annual increase was higher in Collegio Romano than in Casalotti. In this case, climate aridity (more than local warming) reflects (i) a regional trend toward warming and (ii) the indirect impact of urban expansion in Rome. More generally, climatic trends incorporate a marked variability over time that is typical of the Agro Romano district with evident seasonality (contrasting a dry season including spring and summer and a wet season encompassing autumn and winter). Prolonged and anomalous weather conditions were observed in 2003, with average temperature above average (17.2 • C) and particularly low rainfall (569 mm); the 2006-2007 period was particularly hot with average temperatures above 16.5 • C and rainfall below 600 mm.
Profiling urban climate may benefit from the comparative analysis of weather variables in strictly urban and neighboring rural sites [83][84][85]. Unfortunately, data obtained from several meteorological stations are limited to rainfall and thermometric measures [6,57]. At the same time, urban climate can be regarded as the product of a 'cultural mediation' among the variety of approaches in design and planning, as far as materials and structures are concerned. In these regards, compact urban forms may recover definite micro-climate conditions [86][87][88]. A detailed analysis of the urban climate would help an organization of different kinds of settlement along the urban-rural gradient [33,34], e.g., a standard classification based on Local Climatic Zones [61]. Urban climatic regimes have features that explain a spatially-explicit approach for enduring monitoring and adaptation strategies. For instance, new technologies can be applied to the investigation of the effects of urban areas on the formation of precipitation, clouds and the storms [55][56][57], under the assumption that urbanization affects precipitation regimes by increasing the hygroscopic nuclei of condensation of atmospheric moisture. This was due to the increased air turbulence produced by the improved 'roughness' of the land surface and to convection triggered by the properties and different thermal states of the materials [79]. Furthermore, air conditioning of buildings, vehicular traffic and the quality of material for the covering of land surfaces contribute to heating, although the scarcity of green areas and low ventilation diminish the effectiveness of the natural systems for mitigation during these extreme events [37,42]. This hints at negative effects of climate change that can be worsened in urban zones and reduced along the urban gradient [89]. Focusing on natural resources, short and medium-term climatic changes and especially climate aridity may impact on water availability, agriculture water requirements and production of agro-forest landscapes in the Mediterranean region [89][90][91]. 'Dry' and 'hot' years are the most serious time periods since they feature harsh situations possibly harmful for vegetation and crop growth, particularly during the growing season. This information may cause overexploitation of water resources [46], e.g., from groundwater, triggering critical state due to rising withdrawals by civil use and agriculture [20] and other problems, e.g., salinization of groundwater [19].
Multi-scalar mitigation strategies at both national, regional and local scale are essential to manage with climate changes in metropolitan regions through explicit measures to adapt to climate change [92,93]. Policies for mitigation and adaptation to climate changes should improve sustainable strategies aimed at the most effective use of natural resources close to areas with high human pressure [37,42,46,[94][95][96]. Approaches based on long-term climate forecasts and integrated with local-scale analysis of spatial variability in climate trends provide meaningful information to scenario investigation for future stability of peri-urban natural and agricultural environments in the Mediterranean basin [31,[97][98][99][100][101][102][103][104]. Rome's expansion was assessed in the last years by focusing on the spatial demographic dynamics and the decline of the agricultural sector. Implications of global/local warming on peri-urban agriculture should be further investigated with the final objective to inform municipalities' policies supporting a sustainable land management [47,96,[104][105][106][107][108][109].

Conclusions
This study provides insights in local-scale climate regimes, distinguishing exogenous impact of global warming and regional influence of metropolitan expansion in a large urban region in Italy. Empirical evidence indicates a substantial impact of urban expansion in shaping local-scale climate variability, mainly temperature. Fringe agricultural land is recognized as effective green spaces mitigating the 'urban heat' effect. Design of more effective policies linking green infrastructures with more sustainable urban design are particularly relevant in mitigation and adaptation of peri-urban communities exposed to rapid urban expansion and accelerated climate variability.