Is the comparative analysis of temperature trends at the Geophysical Observatory of Modena different from that at the Monte Simone Observatory, Italy? Which analysis is better and more appropriate? Can a favorable conclusion be drawn for the climate future?

Global warming has become a critical environmental, social and economic threat, with the increasing frequency and intensity of extreme weather events. This study analyses temperature trends and climate indices in the Po Valley, an important economic and agricultural region in Italy, by examining data from two historical stations: the Modena Urban Observatory and the Monte Simone Rural Observatory. The analysis extends previous studies up to 2018, assessing the magnitude of climate change since the 1950s and isolating the urban heat island (UHI) effect in Modena. Significant warming trends were confirmed in both locations, with maximum (TX) and minimum (TN) temperature trends almost doubling from 1981 to 2018 compared to 1951-2018. For example, the Texas temperature trend at Modena was 0.84 Cdecade−1 and at Mont Simon was 0.62 Cdecade−1, while the Northeast temperature trend was 0.77 and 0.80 Cdecade−1, respectively. Extreme weather indices showed a significant increase in warm days and nights (TX90p and TN90p, respectively). Specifically, we found TX90p at Modena to be 27.5 days decade−1... Specifically, frost days decreased by 1.88 days decade−1 (37% of urban participation, UC), tropical nights increased by 5.16 days decade−1 (57% UC), warm nights increased by 12.7 days decade−1 (65% UC), and cold nights decreased by 3.19 days decade−1 (39% UC). Overall, this study highlights the importance of considering both global and local factors in the analysis of regional climate trends. Keywords :Apennines, climate change, climate extremes, global warming, Po Valley, UHI, urban health Island. INTRODUCTION Global warming has become increasingly apparent since the last decades of the previous century, posing an environmental, social, and economic threat to the planet due to its widespread and intensifying effects. Several studies in scientific literature have found that global warming is increasing the frequency of extreme weather events, including heatwaves, intense precipitation, floods, and droughts (Planton et al., 2008; Stott, 2016; Xie et al., 2015). Regarding Europe, heat waves in recent decades, especially in 2003, have had unprecedented impacts, causing significant social, economic, and environmental damage (Della-Marta et al., 2007a, 2007b; Klein Tank et al., 2005; Luterbacher et al., 2004). In Italy, temperature trends analysed over various regions show a consistent upward trend, with notable increases in recent decades (Scorzini & Leopardi, 2019; Toreti et al., 2010; Ventura et al., 2002). Researchers have also used extreme indices like the Diurnal Temperature Range (DTR) to study climate patterns, highlighting a reduction in temperature ranges due to more significant increases in minimum temperatures (Bartolini et al., 2008; Brunetti et al., 2000). The Expert Team on Climate Change Detection and Indices (ETCCDI) (Klein Tank et al., 2009) has identified a set of 27 climate indices crucial for analysing temperature and precipitation, with a specific focus on extreme characteristics (Chervenkov & Slavov, 2021; Kang et al., 2014; Klein Tank et al., 2005; Meehl et al., 2000; Sillmann et al., 2013a, 2013b). The present study fits within this context, aiming to analyse temperature trends and ETCCDI indices to estimate local climate change. The area considered is the Po Valley, which hosts a third of Italy's population. This territory has been experiencing increased economic and agricultural vulnerability due to the rise in annual temperatures from the mid-20th century to the present with significant impacts on crop production and in the industrial/civil sector (i.e., increased cooling demands for systems) (Tomozeiu et al., 2006; Zullo et al., 2019). Thus, the paper aims to assess both the magnitude of the climate changes experienced by the territory since the 1950s and to isolate the signal of urbanization in the temperature series to understand its contribution to global warming, as done by Zhang et al. (2021). For this purpose, we consider two historical stations with a long-standing tradition of measuring and recording meteorological parameters located in the same geographic region (Emilia Romagna region, northern Italy) but with distinct characteristics. The first station is the Geophysical Observatory of the University of Modena and Reggio Emilia (hereafter Modena Observatory), established in 1826 within the city centre of Modena in the Po Valley, an area that has experienced considerable urban development. The second station is the Mount Cimone Observatory (hereafter Cimone Observatory), an historical station of the Italian Air Force Meteorological Service. This observatory, located at an elevated position of 2165 m a.s.l. on Mount Cimone, has been collecting data since 1947. It resides in the free troposphere, free from urban influences (Carbone et al., 2014; Ramponi, 2023). This study starts evaluating climate changes in the two locations extending the analysis of Boccolari and Malmusi (2013) to 2018 for Modena, and also includes an assessment of climate changes at Mount Cimone. It also incorporates insights from Cundari and Colombo (1992), who correlated temperature data from the Cimone Observatory with a large area in the Po Valley, demonstrating its utility for long-term temperature representation. Despite extensive use of the Cimone Observatory for CO2 and pollution studies (Alemanno et al., 2014; Bonasoni et al., 1995, 2002; Cristofanelli et al., 2018; Fratticioli et al., 2023), there is a gap in the literature on its meteorological data since Colombo's study Finally, we analysed the potential contribution of urbanization to the temperature trends at the Modena Observatory station. The expansion of built-up areas near the station over time could have caused temperature increases attributable not to global warming but to the Urban Heat Island (UHI) effect (Barbieri et al., 2018; Costanzini et al., 2022). To assess the impact of the UHI, auxiliary data from the ERACLITO databases were utilized (Antolini et al., 2016). This comprehensive approach ensures that the observed temperature trends can be accurately attributed, distinguishing between global climate change and local urban influences. 2 | STUDY AREA AND DATASETS This section outlines the data sources employed in the study, focusing on the two distinct weather stations of Modena Observatory and Cimone Observatory shown in Figure 1. Both stations play crucial roles in assessing climate changes and offer unique insights due to their disparate locations. 2.1 | Meteorological stations 2.1.1 | Modena Observatory The Modena Observatory is housed in the East Tower of the Ducal Palace in Modena's city centre (44380 52.5900N, 10550 47.2200E WGS84, at 34.6 m a.s.l.). The meteorological instruments are located on a historic balcony, positioned at a barometric height of 64.2 m a.s.l. Established in 1826 by Francesco IV, Duke of Modena, the observatory began collecting meteorological data in 1827, encompassing parameters such as temperature, humidity, precipitation, pressure, wind speed and direction, and cloud cover. Past temperature records are comparable with present ones since 1861, marking the construction of the first meteorological window, later replaced by the historical balcony. Recognized as a Centennial Observing Station by the World Meteorological Organization (WMO), Modena Observatory attests to its longstanding commitment to meteorological observations (Lombroso et al., 2020; World Meteorological Organization, 2022). The city of Modena resides in the Po Valley, marked by extensive urban areas, agricultural fields, intense breeding, and wide manufacturing districts. The region exhibits a warm temperate subcontinental climate (Köppen–Geiger classification Cfa) (Peel et al., 2007), characterized by hot summers and frequent temperature inversions, particularly in the cold period (Caserini et al., 2017). Prevailing winds, generally weak, follow two main directions: WNW and ESE, aligning with the longitudinal axis of the Po Valley. Anthropogenic emissions, combined with topographic and meteorological conditions, contribute to poor air quality (Costanzini et al., 2018; Johnstone & Dawson, 2010). The Northern Apennines act as a climatological divide between continental Europe to the north and the Mediterranean Basin to the south (Cristofanelli et al., 2018). Modena Observatory characterizes an urban context due to its location in the city centre. Modena features a hot summer climate with intense droughts in July and August, occasionally interrupted by severe thunderstorms, along with cold and wet winters. The transitional seasons are typically rainy, with peaks in fall (October) and spring (April and May). 2.1.2 | Mount Cimone Observatory The Cimone Observatory, situated at 2165 m a.s.l. on the highest peak in the Northern Apennines (44110 370 N, 10420 000 E, WGS84), serves as a strategic location for the Italian Air Force Meteorological Service. Operational since 1937, this observatory has played and currently plays a crucial role in telecommunications, meteorological forecasting, climatology, and air navigation assistance. Positioned above the Planetary Boundary Layer, in the lower free troposphere, and isolated from anthropogenic activities,the Cimone Observatory holds significance for studying atmospheric composition, especially in measuring background levels of greenhouse gases (Carbone et al., 2014). It is an active participant in the Global Atmosphere Watch program (Cristofanelli et al., 2019; Galli et al., 2019) of the WMO and is part of the Italian Meteorological Network managed by the Air Force Meteorological Service (SYNOP 16134-Metar LIVC). The climate of Mt. Cimone is classified as alpine, with minimum temperatures ranging from −22C in winter to 18C in summer (Colombo et al., 2000; Tositti et al., 2014). Prevailing wind directions include SW and NE, with speeds reaching intensities of 216 kmh−1 , resulting in a perceived wind chill temperature as low as −45C (Ciattaglia et al., 2010). Mt. Cimone experiences two rainfall peaks, with the maximum occurring in November and a secondary peak in April. 2.1.3 | Synoptic patterns of the two stations From a meteorological standpoint, the synoptic weather patterns in Modena and Mt. Cimone exhibit significant similarities, but with some notable exceptions. Instances of heavy rainfall in both locations can occur when a lowpressure system is situated above the Tyrrhenian Sea. Winter cold air flows lead to snowy precipitation in both locations. Strong winds, particularly from the southwest, affect Mt. Cimone. Under these conditions, the Tyrrhenian side of the Apennines experiences cloudiness due to the Stau effect, while Modena may experience the Föhn phenomenon. During intense southwest flows, heavy precipitation has been observed as a result of the spillover effect even at distances 10–20 km downwind from the ridge (Lombroso & Fazlagic, 2000). In both locations, the influence of large highpressure fields fosters stable weather conditions. Summer heat waves are driven by the presence of a strong African anticyclone, which causes very high temperatures in both locations, that is, in June, July, and August 2003, August 2017, and June 2019. However, summer in Mt. Cimone is also characterized by convective systems, in some cases with orographic storms (Frontero & Lombroso, 1988). Stable conditions lead to robust thermal inversions during winter in Modena, resulting in lower temperatures compared to Mt. Cimone. During winter, the phenomenon of thermal inversion is frequently observed, resulting in notable differences between the two stations. This aspect will be further explored in this paper to enhance our understanding of the temperature series and trends calculations. 2.2 | Dataset For this study, the primary datasets comprise daily historical temperature records, encompassing TX and TN values from 1951 to 2018. The DTR series was also derived to provide a comprehensive perspective. The Modena Observatory has daily temperature data available from 1861 to the present. To ensure data consistency, the homogenization process was applied to daily TX and TN, following a method by Boccolari and Malmusi (2013). This method underwent some revisions to minimize undocumented corrections made to the original series. In the original study covering 1861 to 2010, homogenization was performed using statistical tests (SNH, Pettitt, Buishand, and Von Neumann) on mDTR (annual average of DTR) and vDTR (annual average of absolute day-to-day differences of the DTR) series (Wijngaard et al., 2003). These tests assess the statistical significance of temperature series homogeneity and identify breakpoints (apart from the Von Neumann test). Corrections, whether documented or not, aim to adjust the series for homogeneity. Notably, more corrections were needed after 1900 compared to the initial 40 years. Evaluation of the first 40 years is challenging due to extensive documented interventions on measuring instruments. Although mathematically homogeneous, corrections in this period may not convincingly show consistent trends. Subsequent corrections for the 1901–2020 period, using the same methodology, yielded a satisfactory outcome. Only two corrections were needed over the 120-year period, passing all four statistical tests at a 0.01 signifi cance level, specifically for the mDTR series. The rejection of vDTR is not surprising, given its relation to the daily variability of temperature. The corrections made symmetrically on TN and TX were as follows: • Until July 31, 1967 [TX + 0.3C] and [TN − 0.3C]; a documented correction corresponding to the calibration of the maxima and minima thermometer. • Until December 31, 1953 [TX + 0.15C] and [TN − 0.15C]; an undocumented correction. Comparing such new TX and TN trends with those obtained with the previous homogenization by Boccolari and Malmusi (2013), differences of only a few hundredths of a Celsius degree per decade were recorded. Therefore, the results presented in that paper can still be considered valid. As for the Cimone Observatory, meteorological observations began in 1947 but became systematic only since 1951. The TX and TN series for this observatory were extracted from the European Climate Assessment and Dataset (ECA&D) platform (Klein Tank et al., 2002). This dataset includes homogenized series derived from data collected by climatological divisions of the National Meteorological and Hydrological Services, observatories, and research centres across Europe and the Mediterranean. The available TX and TN series include data from January 1, 1951 to December 31, 2020. Notably, there are missing data in the two series, accounting for 352 days for TX and 333 for TN. To address this, a control procedure was applied (Desiato et al., 2012; Dufek et al., 2008) to discard years with more than 30% missing data (i.e., 2019 and 2020 with 63% and 68%, respectively). Uniformly spread gaps over various years were filled through linear interpolation with the other available data. Additionally, the homogeneity of these series was checked using the same method mentioned above, and no corrections were found to be necessary. The datasets from Modena Observatory and Cimone Observatory, including the computation of seasonal and annual anomalies for TX, TN, and DTR, are presented in section 4.2. 3 | METHODOLOGY The methodology applied in this study starts from the examination of TX and TN series recorded by both stations spanning from 1951 to 2018. Initially, computations were made for seasonal and annual temperature anomalies, utilizing cross-correlation techniques to identify shared tendencies between the two stations. An analysis about thermal inversions frequency for Modena Observatory has been also included in this section to help the assessment of the differences achieved in winter between the two stations. Subsequently, an iterative procedure that systematically computed trends and their statistical significance across various periods was employed to partition the time series into meaningful subperiods. Then trends over TX and TN are evaluated, and extreme climate indices from ETCCDI are calculated. The tendencies of these indices and their significance were assessed to elucidate the nuanced changes observed in both stations. Finally, the UHI contribution has been evaluated using the Urban Minus Rural method (Bian et al., 2014; Manalo et al., 2022; Park et al., 2017; Wu et al., 2019; Zhong et al., 2023), to assess the contribution of urbanization on temperatures and ETCCDI indices trends. 3.1 | Occurrence of daytime temperature inversions in Modena The Po Valley, where the Modena Observatory is situated, experiences limited ventilation due to its morphology. During the winter, also a weak radiative forcing occurs, contributing to the persistence of stagnant events characterized by the phenomenon of thermal inversions (Caserini et al., 2017) in these situations, lower temperatures can be recorded in Modena compared to Mt. Cimone. To gain deeper insights into these conditions, an analysis of the seasonal frequency of thermal inversions in Modena was conducted, utilizing data from both observatories. For each day, the lapse rate inC/100 m was calculated, representing the difference between the temperatures at Mt. Cimone and Modena. This calculation included minimum, maximum, and mean temperatures for every day from 1951 to 2018. We started examining periods where strong stagnant events have been documented, as during January 1989 and 1990, February 1993, November 1994, January 2012, December 2013, 2015 and 2016, etc. Although thermal inversions typically cease at 1000–1500 m, during these events minimal temperature differences between Modena and Cimone have been recorded. In these periods, the presence of thermal inversions has also been assessed observing radiosonde data. Based on these findings, two thresholds have been selected to identify the presence of thermal inversion in Modena. Aggregating data from all available days across all 12 months spanning from 1951 to 2018, the frequencies of different lapse rate were computed, taking into account the seasonal distribution. Fre quency are reported grouped for meteorological seasons, as: • Winter (DJF): months of December, January, and February. • Spring (MAM): months of March, April, and May. • Summer (JJA): months of June, July, and August. • Fall (SON): months of September, October, and November. Note that winter season includes the months of December, January, and February. We consider the month of December of 1 year and the January–February of the following year, comprising one winter season. The results of this analysis, useful for the assessment of differences in temperature anomalies between the two sites (and subsequent trends), are reported in section 4.1. 3.2 | TX and TN anomalies and DTR Daily series of TX and TN temperatures from 1951 to 2018 for both stations were used to compute daily anomalies. This was done by subtracting the climatological mean for each specific date from the daily values of each series. The climatological mean was based on the Climatic Normal (CLINO) for the period from 1961 to 1990, as defined by the World Meteorological Organization (WMO) as the standard reference period for long-term climate change assessments. To smooth out large variations due to insufficient samples, a 5-day moving window centred on each specific date was used to average the values. The daily anomalies were subsequently grouped by meteorological seasons to better understand temporal variations and differences between the two stations. A cross-correlation was computed between the series of seasonal and annual TX and TN anomalies of Modena Observatory and Cimone Observatory. Furthermore, the DTR was computed as the difference between daily TX and the TN at both the Modena Observatory and Cimone Observatory. This calculation provides insights into the variation between the highest and lowest temperatures recorded in a single day, offering a more comprehensive perspective on the temperature dynamics at the two stations. The DTR is particularly valuable in discerning differences in temperature increases between minimum and maximum values. The results are presented in section 4.2. 3.3 | Temperature trends The analysis of trends in time series is essential in the context of climate change as it provides insights into the evolution of a phenomenon over a specific period. Trends TX, TN, and DTR were evaluated within the time series of both Modena and Mt. Cimone. This analysis firstly covered the entire period from 1951 to 2018, and then specific significant periods were selected based on the subsequent analysis. The goal is to pinpoint statistically significant trends, typically at a confidence level of 99%, by dividing the entire observation period in time intervals covering at least 10 years within this range. This exploration was carried out through an automated procedure that systematically computed trends and their statistical significance across various sequentially periods. These periods include, for example: 1951–2018, 1952– 2018, 1953–2018, …, 2009–2018, 1951–2017, 1952–2017, …, 2008–2017, …, 1951–1961, 1952–1961, 1951–1960. The procedure reports the results on a diagram called a “shifting trend,” indicating the subperiods to consider for trend calculations. To evaluate the presence and significance of trends on temperatures on the periods identified above, the modified Mann–Kendall test (Kendall, 1957; Pohlert, 2023; Wilks, 2019) was applied, while the Theil– Sen method (Sen, 1968) was used to estimate the slope (section 4.3). These methods are widely used in literature to pursue this aim (Alexander et al., 2006; Mallick et al., 2022; Patra & Satpati, 2022). Indices from the ETCCDI were selected and applied to the temperature series of both stations. These indices including fixed-value thresholds, absolute measures, percentiles, and duration parameters (Alexander et al., 2006; Klein Tank et al., 2009) were employed to analyse the shifts in extreme climate trends. Fixed threshold indices indicate the number of days per year when a certain condition occurred, whereas percentile indices tally the annual days surpassing or falling below a specified percentile threshold. To compute these indices MATLAB® scripts were developed and applied. For the Modified Mann–Kendall test MATLAB® pre-existing software was retrieved (Aalok, 2023), developed, and applied. Table 1 lists the selected indices. The extreme climate indices were calculated based on the TX and TN time series for both stations, spanning the entire observational period as well as other significant periods identified by the shifting trend analysis. The trends are expressed in days per decade, and their statistical significance at the 99% and 95% confidence levels was determined using the modified Mann– Kendall test. Comprehensive results are presented in section .

CONCLUSIONS:

In this study, the analysis focused on the maximum and minimum temperature series of two historical stations in the Po Valley. The first station is located in an urban area, in the centre of Modena, while the second is a station in the lower free atmosphere, the Monte Cimone Observatory. The purpose was to assess the magnitude of climate changes experienced by both stations since the 1950s and to isolate the urban effect due to the presence of extensive residential and industrial areas in the territory by evaluating the UHI effect in Modena. From the study, we found positive anomalies in both TX and TN series for both stations, with TX displaying slightly higher hotspots. Subsequently, the analysis explores into the trends of maximum and minimum temperature series over different time periods: 1951–2018 and 1981–2018. TX and TN trends in the period 1981– 2018 are almost double those computed on the period 1951–2018, confirming the anomalous rise in temperatures in recent years in both sites. We retrieved TX annual temperature trends (1981–2018) reaching 0.84Cdecade−1 for Modena and 0.62Cdecade−1 for Mt. Cimone while for TN values of 0.77Cdecade−1 for Modena and 0.80Cdecade−1 for Mt. Cimone. This general warming trend is reflected in the calculation of extreme climate indices. For instance, warm days (TX90p) in Modena increased by 27.5 daysdecade−1 in the period 1981–2018, almost twice the increase observed at the Cimone Observatory (15 daysdecade−1 ). Similarly, the rise in warm nights (TN90p) is more pronounced in Modena, with an increase of 29.5 daysdecade−1 in the period 1981–2018, compared to 22 daysdecade−1 at the Cimone Observatory. These results have been confirmed also by the comparison between the results of the present paper (in the period 1981–2018) with the findings of a past study focused on Modena in the period 1981–2010. Building on the findings of Cundari and Colombo (1992), which established that a single station (Cimone Observatory) can serve as a representative for studying long-term temperature changes in the Po Valley, it can be affirmed that this holds true in broad terms. However, when examining temperature trends in Modena, it becomes apparent that retrieving the urban influence reveals notable differences. This station, located in the city centre, has been subjected to increasing urbanization over the years. The attempt to quantify the UHI effect in terms of Urbanization Effect and Urbanization Contribution using the reference period 1961–2018, as earlier data were not available, reveals that the UHI effect is significant concerning the rise in minimum temperatures, with notably annual UE and UC, respectively: FD −1.88 no. daysdecade−1 and 37%, TR 5.16 no. daysdecade−1 and 57%, TN90 12.7 no. daysdecade−1 and 65%, and TN10 – 3.19 no. daysdecade−1 and 39%. These results align with studies found in the literature. In conclusion, this research has contributed to characterizing the Po Valley in terms of urban and free atmosphere stations, highlighting how the impact of climate change affects the two sites. During the last period, the increase in trends for minimum temperatures (particularly notable for Mt. Cimone), the negative diurnal temperature range (DTR) trend observed on Mt. Cimone, and the overall assessment of trends on extreme indices (although some are higher in absolute terms for Modena, they require adjustment for Urban Effects) might even indicate a more pronounced warming trend for Mt. Cimone. These findings underscore the importance of considering diverse local influences in the analysis of regional climate trends and provide a significant foundation for formulating targeted adaptation and mitigation strategies. Future outlooks should include analysing changes in precipitation trends, with a focus on northern Italy's climatic divide, to gain a nuanced understanding of the climatic changes the territory is experiencing.

DATA AVAILABILITY STATEMENT The Modena data were provided by the Geophysical Observatory of Modena, University of Modena and Reggio Emilia, Italy (www.ossgeo.unimore.it), and are available upon request from the corresponding author. The Cimone Observatory data are available at https://www. ecad.eu. ORCID Sofia Costanzini https://orcid.org/0000-0003-2227-622X Mauro Boccolari https://orcid.org/0000-0003-1837-5081 Stephanie Vega Parra https://orcid.org/0009-0009- 3850-5817 Francesca Despini https://orcid.org/0000-0002-6813- 131X Luca Lombroso https://orcid.org/0009-0006-3547-8508 Sergio Teggi https://orcid.org/0000-0001-7375-0599 REFERENCES Aalok, A. (2023) Modified Mann–Kendall test. GitHub. Available from: https://github.com/atharvaaalok/Modified-MannKendallTest/releases/tag/v2.0 [Accessed on 10th January 2024] Acquaotta, F., Fratianni, S. & Garzena, D. (2015) Temperature changes in the North-Western Italian Alps from 1961 to 2010. Theoretical and Applied Climatology, 122(3), 619–634. Available from: https://doi.org/10.1007/s00704-014-1316-7 Alemanno, M., Di Diodato, A., Lauria, L. & Santobuono, N. (2014) Environmental measurements at Monte Cimone GAW station. International Journal of Global Warming, 6(4), 424–454. Available from: https://doi.org/10.1504/IJGW.2014.066048 Alexander, L.V., Zhang, X., Peterson, T.C., Caesar, J., Gleason, B., Klein Tank, A.M.G. et al. (2006) Global observed changes in daily climate extremes of temperature and precipitation. Journal of Geophysical Research: Atmospheres, 111(D5), D05109. Available from: https://doi.org/10.1029/2005JD006290 Antolini, G., Auteri, L., Pavan, V., Tomei, F., Tomozeiu, R. & Marletto, V. (2016) A daily high-resolution gridded climatic data set for Emilia-Romagna, Italy, during 1961–2010. International Journal of Climatology, 36(4), 1970–1986. Available from: https://doi.org/10.1002/joc.4473 Arias, P.A., Bellouin, N., Coppola, E., Jones, R.G., Krinner, G., Marotzke, J. et al. (2021) Technical summary. In: MassonDelmotte, V., Zhai, P., Pirani, A., Connors, S.L., Pean, C., Berger, S. et al. (Eds.) Climate change 2021: the physical science basis. Contribution of Working Group I to the sixth assessment report of the Intergovernmental Panel on Climate Change. Cambridge and New York, NY: Cambridge University Press, pp. 33–144. Barbieri, T., Despini, F. & Teggi, S. (2018) A multi-temporal analyses of land surface temperature using Landsat-8 data and open source software: the case study of Modena, Italy. Sustainability, 10(5), 1678. Available from: https://doi.org/10.3390/ su10051678 Bartolini, G., Morabito, M., Crisci, A., Grifoni, D., Torrigiani, T., Petralli, M. et al. (2008) Recent trends in Tuscany (Italy) summer temperature and indices of extremes. International Journal of Climatology, 28(13), 1751–1760. Available from: https://doi. org/10.1002/joc.1673 Bian, T., Ren, G., Zhang, B. & Zhang, L. (2014) Urbanization effect on long-term trends of extreme temperature indices at Shijiazhuang station, North China. Theoretical and Applied Climatology, 115, 407–418. Available from: https://doi.org/10.1007/ s00704-014-1127-x Bigi, A. & Ghermandi, G. (2014) Long-term trend and variability of atmospheric PM10 concentration in the Po Valley. Atmospheric Chemistry and Physics, 14(10), 4895–4907. Available from: https://doi.org/10.5194/acp-14-4895-2014 Bigi, A. & Ghermandi, G. (2016) Trends and variability of atmospheric PM2.5 and PM10–2.5 concentration in the Po Valley, Italy. Atmospheric Chemistry and Physics, 16(24), 15777– 15788. Available from: https://doi.org/10.5194/acp-16-15777- 2016 Boccolari, M. & Malmusi, S. (2013) Changes in temperature and precipitation extremes observed in Modena, Italy. Atmospheric Research, 122, 16–31. Available from: https://doi.org/10.1016/j. atmosres.2012.10.022 Bonasoni, P., Cristofanelli, P., Calzolari, F., Evangelisti, F., Bonafe, U., Petritoli, A. et al. (2002) The Mt. Cimone research station. Paper presented at Workshop on atmospheric research at the Jungfraujoch and in the Alps. Bonasoni, P., Giovanelli, G., Evangelisti, F., Colombo, T., Santaguida, R., Corazza, E. et al. (1995) Study of ozone transport at Mt. Cimone. In: Air pollution and visibility measure ments. Munich, Germany: SPIE, pp. 460–467. Brugnara, Y., Auchmann, R., Brönnimann, S., Bozzo, A., Berro, D.C. & Mercalli, L. (2016) Trends of mean and extreme temperature indices since 1874 at low-elevation sites in the southern Alps. Journal of Geophysical Research: Atmospheres, 121(7), 3304–3325. Available from: https://doi.org/10.1002/ 2015JD024582 .Brunet, M., Sigro, J., Saladie, O., Aguilar, E., Jones, P., Moberg, A. et al. (2005) Long-term change in the mean and extreme state of surface air temperature over Spain (1850–2003). Geophysical Research Abstracts, 7, 04042. Brunetti, M., Buffoni, L., Maugeri, M. & Nanni, T. (2000) Trends of minimum and maximum daily temperatures in Italy from 1865 to 1996. Theoretical and Applied Climatology, 66(1), 49–60. Available from: https://doi.org/10.1007/s007040070032 Brunetti, M., Maugeri, M., Monti, F. & Nanni, T. (2006) Temperature and precipitation variability in Italy in the last two centuries from homogenised instrumental time series. International Journal of Climatology, 26(3), 345–381. Available from: https:// doi.org/10.1002/joc.1251 Büttner, G., Kosztra, B., Maucha, G., Pataki, R., Kleeschulte, S., Hazeu, G. et al. (2021) User manual: Copernicus land monitoring service_corine land cover. European Environment Agency (EEA). https://land.copernicus.eu/en/technical-library/clcproduct-user-manual/@@download/file Carbone, C., Decesari, S., Paglione, M., Giulianelli, L., Rinaldi, M., Marinoni, A. et al. (2014) 3-year chemical composition of free tropospheric PM1 at the Mt. Cimone GAW global station– South Europe–2165 m asl. Atmospheric Environment, 87, 218–227. Caserini, S., Giani, P., Cacciamani, C., Ozgen, S. & Lonati, G. (2017) Influence of climate change on the frequency of daytime temperature inversions and stagnation events in the Po Valley: historical trend and future projections. Atmospheric Research, 184, 15–23. Available from: https://doi.org/10.1016/j.atmosres. 2016.09.018 Chervenkov, H. & Slavov, K. (2021) In: Dimov, I. & Fidanova, S. (Eds.) ETCCDI climate indices for assessment of the recent climate over Southeast Europe. Cham: Springer, pp. 398–412. Chung, Y.-S., Yoon, M.-B. & Kim, H.-S. (2004) On climate variations and changes observed in South Korea. Climatic Change, 66(1), 151–161. Available from: https://doi.org/10.1023/B: CLIM.0000043141.54763.f8 Ciattaglia, L., Cundari, V. & Colombo, T. (2010) Atmospheric carbon dioxide measurements at Mt. Cimone, Italy: 1979–1983. Tellus B, 39, 13–20. Available from: https://doi.org/10.1111/j. 1600-0889.1987.tb00266.x Ciccarelli, N., von Hardenberg, J., Provenzale, A., Ronchi, C., Vargiu, A. & Pelosini, R. (2008) Climate variability in northwestern Italy during the second half of the 20th century. Global and Planetary Change, 63(2), 185–195. Available from: https:// doi.org/10.1016/j.gloplacha.2008.03.006 Colangelo, G., Sanesi, G., Mariani, L., Parisi, S.G. & Cola, G. (2022) A circulation weather type analysis of urban effects on daily thermal range for Milan (Italy). Atmosphere, 13(9), 1529. Available from: https://doi.org/10.3390/atmos13091529 Colombo, T., Santaguida, R., Capasso, A., Calzolari, F., Evangelisti, F. & Bonasoni, P. (2000) Biospheric influence on carbon dioxide measurements in Italy. Atmospheric Environment, 34(29), 4963–4969. Available from: https://doi.org/10. 1016/S1352-2310(00)00366-6 Costanzini, S., Despini, F., Beltrami, L., Fabbi, S., Muscio, A. & Teggi, S. (2022) Identification of SUHI in urban areas by remote sensing data and mitigation hypothesis through solar reflective materials. Atmosphere, 13(1), 70. Available from: https://doi. org/10.3390/atmos13010070 Costanzini, S., Teggi, S., Bigi, A., Ghermandi, G., Filippini, T., Malagoli, C. et al. (2018) Atmospheric dispersion modelling and spatial analysis to evaluate population exposure to pesticides from farming processes. Atmosphere, 9(2), 38. Available from: https://doi.org/10.3390/atmos9020038 Cristofanelli, P., Brattich, E., Decesari, S., Landi, T.C., Maione, M., Putero, D. et al. (2018) The “O. Vittori” Observatory at Mt. Cimone: a “lighthouse” for the Mediterranean troposphere. In: Cristofanelli, P., Brattich, E., Decesari, S., Landi, T.C., Maione, M., Putero, D. et al. (Eds.) High-mountain atmospheric research: the Italian Mt. Cimone WMO/GAW Global Station (2165 m a.s.l.). Cham: Springer, pp. 1–14. Cristofanelli, P., Putero, D., Naitza, L., Marinoni, A., Calzolari, F., Roccato, F. et al. (2019) The Italian high-altitude atmospheric observation network. NextData. Cundari, V. & Colombo, T. (1992) Atmospheric temperature time series analysis at Mt. Cimone Observatory, Italy, 1946–1991. Il Nuovo Cimento C, 15(4), 381–390. Available from: https://doi. org/10.1007/BF02511738 Curci, G., Guijarro, J.A., Di Antonio, L., Di Bacco, M., Di Lena, B. & Scorzini, A.R. (2021) Building a local climate reference dataset: application to the Abruzzo region (Central Italy), 1930–2019. International Journal of Climatology, 41(8), 4414– 4436. Available from: https://doi.org/10.1002/joc.7081 Della-Marta, P.M., Haylock, M.R., Luterbacher, J. & Wanner, H. (2007a) Doubled length of western European summer heat waves since 1880. Journal of Geophysical Research: Atmospheres, 112(D15), D15103. Available from: https://doi.org/10. 1029/2007JD008510 Della-Marta, P.M., Luterbacher, J., von Weissenfluh, H., Xoplaki, E., Brunet, M. & Wanner, H. (2007b) Summer heat waves over western Europe 1880–2003, their relationship to large-scale forcings and predictability. Climate Dynamics, 29(2), 251–275. Available from: https://doi.org/10.1007/s00382-007- 0233-1 Desiato, F., Fioravanti, G., Fraschetti, P., Perconti, W. &Piervitali, E. (2012) Elaborazione delle serie temporali per la stima delle tendenze climatiche. In: Rapporto ISPRA, Vol. 32. Stato dell'Ambiente. https://www.isprambiente.gov.it/it/ pubblicazioni/stato-dellambiente/elaborazioneserietemporali perlastimadelletendenzedelclima.pdf Dufek, A.S., Ambrizzi, T. & Da Rocha, R.P. (2008) Are reanalysis data useful for calculating climate indices over South America? Annals of the New York Academy of Sciences, 1146(1), 87–104. Available from: https://doi.org/10.1196/annals.1446.010 Field, C.B., Barros, V., Stocker, T.F. & Dahe, Q. (Eds.). (2012) Managing the risks of extreme events and disasters to advance climate change adaptation: special report of the intergovernmental panel on climate change. Cambridge: Cambridge University Press. Fioravanti, G., Piervitali, E. & Desiato, F. (2016) Recent changes of temperature extremes over Italy: an index-based analysis. Theoretical and Applied Climatology, 123(3), 473–486. Available from: https://doi.org/10.1007/s00704-014-1362-1 Fioravanti, G., Piervitali, E., Desiato, F., Perconti, W. & Fraschetti, P. (2013) Variazione e tendenze degli estremi di temperatura e precipitazione in Italia. Rome, Italy: Istituto Superiore per la Protezione e la Ricerca Ambientale-ISPRA. Fratticioli, C., Trisolino, P., Maione, M., Calzolari, F., Calidonna, C., Biron, D. et al. (2023) Continuous atmospheric

in-situ measurements of the CH4/CO ratio at the Mt. Cimone station (Italy, 2165 m a.s.l.) and their possible use for estimating regional CH4 emissions. Environmental Research, 232, 116343. Available from: https://doi.org/10. 1016/j.envres.2023.116343 Frontero, P. & Lombroso, L. (1988) Analisi sulle diverse frequenze temporalesche in Appennino e nelle Prealpi. Paper presented at Atti Della Conferenza Internazionale Di Meteorologia Alpina., 20 Congresso Internazionale di Meteorologia alpina (CIMA 88) Sestola (MO), September 18–25, 1988. Galli, M., Caracciolo di Torchiarolo, L., Proietti, A. & Lauria, L. (2019) Il C.A.M.M. di Monte Cimone Stazione - Stazione Globale del programma Global Atmosphere Watch dell'Organizzazione Meteorologica Mondiale. Revista di Meteorologia Aeronautica, 2, 8–28. Available from: https://www.aeronautica. difesa.it/wp-content/uploads/2023/02/RIVISTA_4_2019.pdf Hartmann, D.L., Klein Tank, A.M.G., Rusticucci, M., Alexander, L.V., Brönnimann, S., Charabi, Y.A.R. et al. (2013) Observations: atmosphere and surface. In: Climate change 2013: the physical science basis. Contribution of Working Group I to the fifth assessment report of the Intergovernmental Panel on Climate Change. Cambridge: Cambridge University Press, pp. 160–254. Johnstone, J.A. & Dawson, T.E. (2010) Climatic context and ecological implications of summer fog decline in the coast redwood region. Proceedings of the National Academy of Sciences of the United States of America, 107(10), 4533–4538. Available from: https://doi.org/10.1073/pnas.0915062107 Kang, K.K., Lee, D.S., Hwang, S.H. & Kim, B.S. (2014) Analysis of extreme weather characteristics change in the Gangwon Province using ETCCDI indices. Journal of Korea Water Resources Association, 47(12), 1107–1119. Available from: https://doi.org/ 10.3741/JKWRA.2014.47.12.1107 Karl, T., Nicholls, N. & Ghazi, A. (1999) CLIVAR/GCOS/WMO workshop on indices and indicators for climate extremes. Weather and Climate Extremes, 42, 3–7. Kendall, M.G. (1957) Review of rank correlation methods. Biometrika, 44(1–2), 298. Available from: https://doi.org/10.2307/ 2333282 Klein Tank, A.M.G., Können, G.P. & Selten, F.M. (2005) Signals of anthropogenic influence on European warming as seen in the trend patterns of daily temperature variance. International Journal of Climatology, 25(1), 1–16. Available from: https://doi. org/10.1002/joc.1087 Klein Tank, A.M.G., Wijngaard, J.B., Können, G.P., Böhm, R., Demarée, G., Gocheva, A. et al. (2002) Daily dataset of 20thcentury surface air temperature and precipitation series for the European climate assessment. International Journal of Climatology, 22(12), 1441–1453. Available from: https://doi.org/10. 1002/joc.773 Klein Tank, A.M.G., Zwiers, F.W. & Zhang, X. (2009) Guidelines on analysis of extremes in a changing climate in support of informed decisions for adaptation. Climate data and monitoring. Geneva: WMO. WMO/TD-No. 1500. Li, H., Meier, F., Lee, X., Chakraborty, T., Liu, J., Schaap, M. et al. (2018) Interaction between urban heat Island and urban pollution Island during summer in Berlin. Science of the Total Environment, 636, 818–828. Available from: https://doi.org/10.1016/ j.scitotenv.2018.04.254 Li, H., Sodoudi, S., Liu, J. & Tao, W. (2020) Temporal variation of urban aerosol pollution Island and its relationship with urban heat Island. Atmospheric Research, 241, 104957. Available from: https://doi.org/10.1016/j.atmosres.2020.104957 Liu, W., Ji, C., Zhong, J., Jiang, X. & Zheng, Z. (2007) Temporal characteristics of the Beijing urban heat Island. Theoretical and Applied Climatology, 87(1), 213–221. Available from: https:// doi.org/10.1007/s00704-005-0192-6 Lombroso, L., Costanzini, S., Despini, F. & Teggi, S. (2020) Annuario 2020 dell'Osservatorio Geofisico di Modena. Atti della Societa dei Naturalisti e Matematici di Modena, 152, 6–35. Lombroso, L. & Fazlagic, S. (2000) Analisi delle precipitazioni alluvionali del 6-7novembre 2000. Atti Della Società Dei Naturalisti e Matematici di Modena, 131, 5–7. Lombroso, L. & Quattrocchi, S. (2008) L' Osservatorio di Modena: 180 anni di misure meteoclimatiche. Available from: www. libreriauniversitaria.it Luterbacher, J., Dietrich, D., Xoplaki, E., Grosjean, M. & Wanner, H. (2004) European seasonal and annual temperature variability, trends, and extremes since 1500. Science, 303(5663), 1499–1503. Available from: https://doi.org/10.1126/science. 1093877 Mallick, J., Islam, A.R.M.T., Ghose, B., Islam, H.M.T., Rana, Y., Hu, Z. et al. (2022) Spatiotemporal trends of temperature extremes in Bangladesh under changing climate using multistatistical techniques. Theoretical and Applied Climatology, 147(1), 307–324. Available from: https://doi.org/10.1007/ s00704-021-03828-1 Manalo, J., Matsumoto, J., Takahashi, H., Villafuerte, M., II, Olaguera, L.M., Ren, G. et al. (2022) The effect of urbanization on temperature indices in The Philippines. International Journal of Climatology, 42, 850–867. Available from: https://doi.org/ 10.1002/joc.7276 Meehl, G.A., Karl, T., Easterling, D.R., Changnon, S., Pielke, R., Changnon, D. et al. (2000) An introduction to trends in extreme weather and climate events: observations, socioeconomic impacts, terrestrial ecological impacts, and model projections. Bulletin of the American Meteorological Society, 81(3), 413–416. Available from: https://doi.org/10.1175/1520-0477(2000) 0812.3.CO;2 Oke, T.R. (1995) The Heat Island of the urban boundary layer: characteristics, causes and effects. In: Cermak, J.E., Davenport, A.G., Plate, E.J. & Viegas, D.X. (Eds.) Wind climate in cities. Dordrecht: Springer, pp. 81–107. Otto, I., Reckien, D., Reyer, C., Marcus, R., Le Masson, V., Jones, L. et al. (2017) Social vulnerability to climate change: a review of concepts and evidence. Regional Environmental Change, 17, 1651–1662. Available from: https://doi.org/10.1007/s10113-017- 1105-9 Park, B., Kim, Y., Min, S., Kim, M., Choi, Y., Boo, K. et al. (2017) Long-term warming trends in Korea and contribution of urbanization: an updated assessment. Journal of Geophysical Research: Atmospheres, 122(20), 10637–10654. Available from: https://doi.org/10.1002/2017JD027167 Patra, P. & Satpati, L. (2022) Characteristics and trend analysis of absolute and relative temperature extremes indices and related indices of Kolkata. Theoretical and Applied Climatology, 148(3)943–954. Available from: https://doi.org/10.1007/s00704-022- 03975-z Peel, M.C., Finlayson, B.L. & McMahon, T.A. (2007) Updated world map of the Köppen-Geiger climate classification. Hydrology and Earth System Sciences, 11(5), 1633–1644. Available from: https://doi.org/10.5194/hess-11-1633-2007 Planton, S., Déqué, M., Chauvin, F. & Terray, L. (2008) Expected impacts of climate change on extreme climate events. Comptes Rendus Geoscience, 340(9), 564–574. Available from: https://doi. org/10.1016/j.crte.2008.07.009 Pohlert, T. (2023) Non-parametric trend tests and change-point detection. CC by ND, 2016(4), 1–18. Ramponi, V. (2023) Lo sguardo scientifico del monte Cimone. Clionet, 7, 187–194. Available from: https://doi.org/10.30682/ clionet2307y Rohli, R.V. & Vega, A.J. (2017) Climatology, 4th edition. Burlington, MA: Jones and Bartlett Learning Publishers. Romano, B., Fiorini, L., Marucci, A. & Zullo, F. (2020) The urbanization run-up in Italy: from a qualitative goal in the boom decades to the present and future unsustainability. Land, 9(9), 301. Scorzini, A.R. & Leopardi, M. (2019) Precipitation and temperature trends over central Italy (Abruzzo region): 1951–2012. Theoretical and Applied Climatology, 135(3), 959–977. Available from: https://doi.org/10.1007/s00704-018-2427-3 Sen, P.K. (1968) Estimates of the regression coefficient based on Kendall's tau. Journal of the American Statistical Association, 63(324), 1379–1389. Available from: https://doi.org/10.1080/ 01621459.1968.10480934 Seneviratne, S.I., Zhang, X., Adnan, M., Badi, W., Dereczynski, C., Di Luca, A. et al. (2021) Weather and climate extreme events in a changing climate. In: Masson-Delmotte, V., Zhai, P., Pirani, A., Connors, S.L., Pean, C., Berger, S. et al. (Eds.) Climate change 2021: the physical science basis. Contribution of Working Group I to the sixth assessment report of the Intergovernmental Panel on Climate Change. Cambridge and New York, NY: Cambridge University Press, pp. 1513–1766. Sillmann, J., Kharin, V.V., Zwiers, F.W., Zhang, X. & Bronaugh, D. (2013a) Climate extremes indices in the CMIP5 multimodel ensemble: part 1. Model evaluation in the present climate. Journal of Geophysical Research: Atmospheres, 118(4), 1716–1733. Available from: https://doi.org/10.1002/ jgrd.50203 Sillmann, J., Kharin, V.V., Zwiers, F.W., Zhang, X. & Bronaugh, D. (2013b) Climate extremes indices in the CMIP5 multimodel ensemble: part 2. Future climate projections. Journal of Geophysical Research: Atmospheres, 118(6), 2473–2493. Available from: https://doi.org/10.1002/jgrd.50188 Sippel, S., Barnes, C., Cadiou, C., Fischer, E., Kew, S., Kretschmer, M. et al. (2024). Could an extremely cold central European winter such as 1963 happen again despite climate change? Weather and Climate Dynamics, 5(3), 943–957. Available from: https:// doi.org/10.5194/wcd-5-943-2024 Stott, P. (2016) How climate change affects extreme weather events. Science, 352(6293), 1517–1518. Available from: https://doi.org/ 10.1126/science.aaf7271 Tomozeiu, R., Pavan, V., Cacciamani, C. & Amici, M. (2006) Observed temperature changes in Emilia-Romagna: mean values and extremes. Climate Research, 31(2–3), 217–225. Toreti, A. & Desiato, F. (2008) Temperature trend over Italy from 1961 to 2004. Theoretical and Applied Climatology, 91(1), 51–58. Available from: https://doi.org/10.1007/s00704-006-0289-6 Toreti, A., Desiato, F., Fioravanti, G. & Perconti, W. (2010) Seasonal temperatures over Italy and their relationship with low-frequency atmospheric circulation patterns. Climatic Change, 99(1), 211–227. Available from: https://doi.org/10. 1007/s10584-009-9640-0 Tositti, L., Brattich, E., Cinelli, G. & Baldacci, D. (2014) 12 years of 7Be and 210Pb in Mt. Cimone, and their correlation with meteorological parameters. Atmospheric Environment, 87, 108–122. Available from: https://doi.org/10.1016/j.atmosenv.2014.01.014 Türke¸s, M., Yozgatlıgil, C., Batmaz, _ I., _ Iyigün, C., Koç, E.K., Fahmi, F.M. et al. (2016) Has the climate been changing in Turkey? Regional climate change signals based on a comparative statistical analysis of two consecutive time periods, 1950– 1980 and 1981–2010. Climate Research, 70(1), 77–93. Available from: https://doi.org/10.3354/cr01410 Twardosz, R. & Kossowska-Cezak, U. (2016) Exceptionally cold and mild winters in Europe (1951–2010). Theoretical and Applied Climatology, 125(1), 399–411. Available from: https://doi.org/ 10.1007/s00704-015-1524-9 Ventura, F., Rossi Pisa, P. & Ardizzoni, E. (2002) Temperature and precipitation trends in Bologna (Italy) from 1952 to 1999. Atmospheric Research, 61(3), 203–214. Available from: https://doi. org/10.1016/S0169-8095(01)00135-1 Wijngaard, J.B., Klein Tank, A.M.G. & Können, G.P. (2003) Homogeneity of 20th century European daily temperature and precipitation series. International Journal of Climatology, 23, 679–692. Available from: https://doi.org/10.1002/joc.906 Wild, M. (2013) Relevance of decadal variations in surface radiative fluxes for climate change. AIP Conference Proceedings, 1531(1), 728–731. Available from: https://doi.org/10.1063/1.4804873 Wild, M. (2016) Decadal changes in radiative fluxes at land and ocean surfaces and their relevance for global warming. WIREs Climate Change7(1), 91–107. Available from: https://doi.org/10.1002/wcc.372 Wilks, D.S. (2019) Statistical methods in the atmospheric sciences. Amsterdam: Elsevier. World Meteorological Organization. (2022) Centennial observing stations: state of recognition report: 2021. Geneva: WMO. Wu, Z., Zhang, Q., Song, C., Zhang, F., Zhu, X., Sun, P. et al. (2019) Impacts of urbanization on spatio-temporal variations of temperature over the Pearl River Delta. Xie, S.-P., Deser, C., Vecchi, G.A., Collins, M., Delworth, T.L., Hall, A. et al. (2015) Towards predictive understanding of regional climate change. Nature Climate Change, 5(10), 921– 930. Available from: https://doi.org/10.1038/nclimate2689 Zhang, P., Ren, G., Qin, Y., Zhai, Y., Zhai, T., Tysa, S.K. et al. (2021). Urbanization effects on estimates of global trends in mean and extreme air temperature. Journal of Climate, 34(5), 1923–1945. Zhang, X., Hegerl, G., Zwiers, F.W. & Kenyon, J. (2005) Avoiding inhomogeneity in percentile-based indices of temperature extremes. Journal of Climate, 18(11), 1641–1651. Available from: https://doi.org/10.1175/JCLI3366.1 Zheng, Z., Ren, G., Wang, H., Dou, J., Gao, Z., Duan, C. et al. (2018) Relationship between fine-particle pollution and the Urban Heat Island in Beijing, China: observational evidence. Boundary-Layer Meteorology, 169, 93–113. Available from: https://doi.org/10.1007/s10546-018-0362-6 .Zhong, R., Song, S., Zhang, J. & Ye, Z. (2023) Spatial–temporal variation and temperature effect of urbanization in Guangdong Province from 1951 to 2018. Environment, Development and Sustainability, 26, 9661–9683. Available from: https://doi.org/ 10.1007/s10668-023-03113-3 Zullo, F., Fazio, G., Romano, B., Marucci, A. & Fiorini, L. (2019) Effects of urban growth spatial pattern (UGSP) on the land surface temperature (LST): a study in the Po Valley (Italy). Science of the Total Environment, 650, 1740–1751. Available from: https://doi.org/10.1016/j.scitotenv.2018.09.331

More Abbas Kashani's questions See All
Similar questions and discussions