Skip to main content
Log in

Characterization of annual urban air temperature changes with special reference to the city of Modena: a comparison between regression models and a proposal for a new index to evaluate relationships between environmental variables

  • ORIGINAL PAPER
  • Published:
Stochastic Environmental Research and Risk Assessment Aims and scope Submit manuscript

Abstract

It is commonly recognized that the observed increase in global mean annual air temperature is strongly related to the increase in global carbon dioxide concentration C, and that both these variables are related to global development. It remains, however, unclear the degree to which local mean annual urban air temperature T is affected by local variables such as annual precipitation depth P and urban area extent A. This study assumes that A is a proxy of local development and C is a proxy of global development and investigates the commingled effects of A, P, and C on T by using long-term annual data observed over the years 1881–2019 from the Modena Observatory in Italy. Linear relationships between T, C and A are found to be spurious since all these series have a monotonic increasing trend with time. Parametric analytic models like logistic functions are found to lack flexibility. Smoothing splines can only give insights into the strength of the relationships but not on their shape defining the functional relationship between variables. Advanced nonlinear models like generalized additive models, instead, are found to combine flexibility in a parametric form, and appear therefore to be suitable models for explaining the complex relationships between A, P, and C on T. The different models are evaluated using traditional goodness of fit statistics like R2, AIC, BIC, and a new index of relation IR which is introduced to jointly evaluate the goodness-of-fit of relationships between variables that may either be dependent or independent.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7

Similar content being viewed by others

Data availability

Data sets used in this study are available at https://doi.org/10.1594/PANGAEA.938740.

References

  • Addinsoft (2023) XLSTAT Statistical and data analysis solutions. New York, USA. https://www.xlstat.com

  • Alexandersson H (1986) A homogeneity test applied to precipitation data. J Climatol 6:661–675

    Article  Google Scholar 

  • Anderson TW (2003) An introduction to multivariate statistical analysis. Wiley Series in Probability and Statistics

  • Barnston AG, Livezey RE (1987) Classification, seasonality and persistence of low-frequency atmospheric circulation patterns. Mon Wea Rev 115:1083–1126

    Article  ADS  Google Scholar 

  • Boccolari MA, Malmusi S (2013) Changes in temperature and precipitation extremes observed in Modena, Italy. Atmos Res 122:16–31

    Article  Google Scholar 

  • Buritica G, Naveau P (2022) Stable sums to infer high return levels of multivariate rainfall time series. Environmetrics. https://doi.org/10.1002/env.2782

    Article  Google Scholar 

  • Burnham KP, Anderson DR, Huyvaert KP (2011) AIC model selection and multimodel inference in behavioural ecology. Behav Ecol Sociobiol 65:23–35. https://doi.org/10.1007/s00265-010-1029-6

    Article  Google Scholar 

  • Buishand TA (1982) Some methods for testing the homogeneity of rainfall records. J Hydrol 58:11–27

    Article  Google Scholar 

  • Buishand TA (1984) Tests for detecting a shift in the mean of hydrological time series. J Hydrol 73:51–69

    Article  Google Scholar 

  • Burt SD (2023) A twice-daily barometric pressure record from Durham observatory in north-east England, 1843–1960. Geosci Data J 10(1):3–17. https://doi.org/10.1002/gdj3.135

    Article  ADS  Google Scholar 

  • Cai Z, Liu Q, Cao S (2020) Real estate supports rapid development of China’s urbanization. Land Use Policy 95:104582

    Article  Google Scholar 

  • Caporali E, Lompi M, Pacetti T, Chiarello V, Fatichi S (2021) A review of studies on observed precipitation trends in Italy. Int J Climatol 41:E1–E25

    Article  Google Scholar 

  • Corradini E (2014) From the geophysical-meteorological Observatory of Modena to the Italian network of observatories. In: Proceedings of the 12th conference of the international committee of ICOM for University Museums and Collections (UMAC), Singapore, 10–12 Oct, 2012. https://doi.org/10.18452/8741

  • Daramola MT, Xu M (2022) Recent changes in global dryland temperature and precipitation. Int J Climatol 42(2):1267–1282. https://doi.org/10.1002/joc.7301

    Article  Google Scholar 

  • De Boor C (2001) A practical guide to splines. Springer

    Google Scholar 

  • Etheridge DM, Steele LP, Langenfelds RL, Francey RJ (1998) Historical CO2 record derived from the Law Dome DE08 and DE08-2 ice cores. https://cdiac.ess-dive.lbl.gov/trends/co2/lawdome-data.html

  • Fawcett L, Walshaw D (2007) Improved estimation for temporally clustered extremes. Environmetrics 18:173–188

    Article  MathSciNet  Google Scholar 

  • Fawcett L, Walshaw D (2012) Estimating return levels from serially dependent extremes. Environmetrics 23:272–283

    Article  MathSciNet  Google Scholar 

  • Ferguson CR, Villarini G (2021) Detecting inhomogeneities in the Twentieth Century Reanalysis over the central United States. J Geophys Res 117: D05123. https://doi.org/10.1029/2011JD016988

    Article  Google Scholar 

  • Findley DF (1991) Counterexamples to parsimony and BIC. Ann Inst Stat Math 43(3):505–514. https://doi.org/10.1007/BF00053369

    Article  MathSciNet  Google Scholar 

  • Friedman J (2001) Greedy function approximation: a gradient boosting machine. Ann Stat 29:1189–1232

    Article  MathSciNet  Google Scholar 

  • Friedman J (2002) Stochastic gradient boosting. Comput Stat Data Anal 38:367–378

    Article  MathSciNet  Google Scholar 

  • Gao Z, Hou Y, Chen W (2019) Enhanced sensitivity of the urban heat island effect to summer temperatures induced by urban expansion. Environ Res Lett 14(9):094005. https://doi.org/10.1088/1748-9326/ab2740

    Article  ADS  CAS  Google Scholar 

  • Gilbert RO (1987) Statistical methods for environmental pollution monitoring. Wiley, New York

    Google Scholar 

  • Granger CWJ (1969) Investigating causal relations by econometric models and cross-spectral methods. Econometrica 37:424–438

    Article  Google Scholar 

  • Granger CWJ, Ghysels E, Swanson NR, Watson MW (2001) Essays in econometrics: collected papers of Clive W. J. Granger. Cambridge University Press

    Book  Google Scholar 

  • Guan Y, Lu H, Jiang Y, Tian P, Qiu L, Pellikka P, Heiskanen J (2021) Changes in global climate heterogeneity under the 21st century global warming. Ecol Indic 130:108075

    Article  Google Scholar 

  • Hastie T, Tibshirani R (1986) Generalized additive models (with discussion). Stat Sci 1:297–318

    Google Scholar 

  • Hastie T, Tibshirani R (1990) Generalized additive models. Chapman & Hall/CRC

    Google Scholar 

  • Hawkins M (1977) Testing a sequence of observations for a shift in location. J Am Stat Assoc 72:180–186

    Article  MathSciNet  Google Scholar 

  • Hipel KW, McLeod AI (1994) Time series modelling of water resources and environmental systems. Elsevier Science, New York

    Google Scholar 

  • Huang Y, Matsumoto KI (2021) Drivers of the change in carbon dioxide emissions under the progress of urbanization in 30 provinces in China: a decomposition analysis. J Clean Prod 322:129000

    Article  CAS  Google Scholar 

  • Jones PD, New M, Parker DE, MartinS RIG (1999) Surface air temperature and its changes over the past 150 years. Rev Geophys 37(2):173–199. https://doi.org/10.1029/1999RG900002

    Article  ADS  Google Scholar 

  • Jouzel J, Masson-Delmotte V, Cattani O, Dreyfus G, Falourd S, Hoffmann G, Minster B, Nouet J, Barnola JM, Chappellaz J, Fischer H, Gallet JC, Johnsen S, Leuenberger M, Loulergue L, Luethi D, Oerter H, Parrenin F, Raisbeck G, Raynaud D, Schilt A, Schwander J, Selmo E, Souchez R, Spahni R, Stauffer B, Steffensen JP, Stenni B, Stocker TF, Tison JL, Wernerm M, Wolff E (2007) Orbital and millennial Antarctic climate variability over the past 800,000 years. Science 317(5839):793–797. https://doi.org/10.1126/science.1141038

    Article  ADS  CAS  PubMed  Google Scholar 

  • Kalnay E, Kay M (2003) Impact of urbanization and land-use change on climate. Nature 423:528–531. https://doi.org/10.1038/nature01675

    Article  ADS  CAS  PubMed  Google Scholar 

  • Karl TR, Diaz HF, Kukla G (1988) Urbanization: its detection and effect in the United States climate record. J Clim 1(11):1099–1123

    Article  ADS  Google Scholar 

  • Kendal MG (1975) Rank correlation methods, 4th edn. Charles Griffin, London

    Google Scholar 

  • Krayenhoff ES, Voogt JA (2010) Impacts of urban albedo increase on local air temperature at daily–annual time scales: model results and synthesis of previous works. J Appl Meteorol Climatol 49:1634–1648

    Article  ADS  Google Scholar 

  • Kwiatkowski D, Phillips PCB, Schmidt P, Shin Y (1992) Testing the null hypothesis of stationarity against the alternative of a unit root. J Econometrics 54(1–3):159–178

    Article  Google Scholar 

  • Li J, Xie H, Li J, Yang G, Xie Y, Wang J, Zhou C, Zou S (2023) Influences of anthropogenic acids on carbonate weathering and CO2 sink in an agricultural karst wetland (South China). Ecol Indic 150:110192

    Article  CAS  Google Scholar 

  • Lou Y, Caruana R, Gehrke J (2012) Intelligible models for classification and regression. In: Proceedings of the 18th ACM SIGKDD international conference on knowledge discovery and data mining (KDD ’12), Beijing, China. ACM Press, pp 150–158

  • Lou Y, Caruana R, Gehrke J, Hooker G (2013) Accurate Intelligible Models with Pairwise Interactions. Proceedings of the 19th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining (KDD ’13) Chicago, Illinois, USA, ACM Press: 623–631

  • Luthi D, Floch ML, Bereiter B, Blunier T, Barnola JM, Siegenthaler U, Raynaud D, Jouzel J, Fischer H, Kawamura K, Stocker TF (2008) High resolution carbon dioxide concentration record 650,000–800,000 years before present. Nature 453:379–382. https://doi.org/10.1038/nature06949

    Article  ADS  CAS  PubMed  Google Scholar 

  • Mann HB (1945) Non-parametric tests against trend. Econometrica 13:163–171

    Article  Google Scholar 

  • McQuarrie ADR, Tsai CL (1998) Regression and time series model selection. World Scientific

    Book  Google Scholar 

  • Mehmood U, Tariq S, Haq ZU (2021) Effects of population structure on CO2 emissions in South Asian countries: evidence from panel estimation. Environ Sci Pollut Res 28:66858–66863

    Article  CAS  Google Scholar 

  • Milly PCD, Betancourt J, Falkenmark M, Hirsch RM, Kundzewicz ZW, Lettenmaier DP, Stouffer RJ (2008) Stationary is dead: whither water management? Science 319:573–574. https://doi.org/10.1126/science.1151915

    Article  CAS  PubMed  Google Scholar 

  • Morlini I, Franco-Villoria M, Orlandini S (2023) Modelling local climate change using site-based data. Environ Ecol Stat 30(2):205–232

    Article  Google Scholar 

  • Oceanic and Atmospheric Administration (2021) Trends in atmospheric carbon dioxide. https://gml.noaa.gov/ccgg/trends/

  • Pearl J (2000) Causality: models, reasoning, and inference. Cambridge University Press

    Google Scholar 

  • Pettitt AN (1979) A non-parametric approach to the change point problem. J R Stat Soc: Ser C Appl Stat 28:126–135

    Google Scholar 

  • Regoto P, Dereczynski C, Chou SC, Bazzanella AC (2021) Observed changes in air temperature and precipitation extremes over Brazil. Int J Climatol 41(11):5125–5142

    Article  Google Scholar 

  • Rehman A, Ma H, Ahmad M, Irfan M, Traore O, Chandio AA (2021) Toward environmental sustainability: devolving the influence of carbon dioxide emission to population growth, climate change, forestry, livestock and crops production in Pakistan. Ecol Indic 125:107460

    Article  CAS  Google Scholar 

  • Reinsch CH (1967) Smoothing by Spline Functions. Numer Math 10(3):177–183. https://doi.org/10.1007/BF02162161

    Article  MathSciNet  Google Scholar 

  • Sen PK (1968) Estimates of the regression coefficient based on Kendall’s tau. J Am Stat Assoc 63:1379–1389

    Article  MathSciNet  Google Scholar 

  • Stoica P, Selen Y (2004) Model-order selection: a review of information criterion rules. IEEE Signal Process Mag. https://doi.org/10.1109/MSP.2004.1311138

    Article  Google Scholar 

  • Schwarz GE (1978) Estimating the dimension of a model. Ann Stat 6(2):461–464. https://doi.org/10.1214/aos/1176344136

    Article  MathSciNet  Google Scholar 

  • Skytt T, Nielsen SN, Jonsson B (2021) Global warming potential and absolute global temperature change potential from carbon dioxide and methane fluxes as indicators of regional sustainability—A case study of Jämtland, Sweden. Ecol Indic 130:108075

    Google Scholar 

  • Verstraeten G, Poesen J, Demaree G, Salles C (2006) Long-term (105 years) variability in rain erosivity as derived from 10-min rainfall depth data for Ukkel (Brussels, Belgium): implications for assessing soil erosion rates. J Geophys Res 111:D22109

    ADS  Google Scholar 

  • Von Neumann J (1941) Distribution of the ratio of the mean square successive difference to the variance. Ann Math Stat 12:367–395

    Article  MathSciNet  Google Scholar 

  • Wahba G (1990) Spline models for observational data. SIAM, Philadelphia

    Book  Google Scholar 

  • Wang Q, Wang L (2021) The nonlinear effects of population aging, industrial structure, and urbanization on carbon emissions: a panel threshold regression analysis of 137 countries. J Clean Prod 287:125381

    Article  Google Scholar 

  • Wang J, Wu Y, Zhao Y, He S, Dong Z, Bo W (2019) The population structural transition effect on rising per capita CO2 emissions: evidence from China. Clim Policy 19(10):1250–1269

    Article  Google Scholar 

  • Wilks SS (1932) Certain generalizations in the analysis of variance. Biometrika 24:471–494

    Article  Google Scholar 

  • Winjngaard JB, Klein Tank AMG, Konnen GP (2003) Homogeneity of 20th century European daily temperature and precipitation series. Int J Climatol 23:679–692

    Article  Google Scholar 

  • Yule GU (1926) Why do we sometimes get nonsense-correlations between time-series? A study in sampling and the nature of time-series. J R Stat Soc 89(1):1–63. https://doi.org/10.2307/2341482

    Article  MathSciNet  Google Scholar 

  • Zou F, Li H, Hu Q (2020) Responses of vegetation greening and land surface temperature variations to global warming on the Qinghai-Tibetan Plateau, 2001–2016. Ecol Indic 119:106867

    Article  Google Scholar 

Download references

Acknowledgements

Air temperature and precipitation time series of Modena were provided by the Geophysical Observatory of the University of Modena and Reggio Emilia, Modena, Italy. Luca Lombroso is greatly acknowledged for providing specific comments and feedback on these data. The research reported in the present paper was supported by Fondazione Cassa di Risparmio di Modena through the grant 2018-0093, by the University of Modena and Reggio Emilia through the grant FAR 2020 Mission Oriented, and by the European Union NextGenerationEU/NRRP, Mission 4 Component 2 Investment 1.5, Call 3277 (12/30/2021), Award 0001052 (06/23/2022), under the project ECS00000033 “Ecosystem for Sustainable Transition in Emilia-Romagna,” Spoke 6 “Ecological Transition Based on HPC and Data Technology.” The authors thank the anonymous reviewers for comments that led to improvement in the manuscript.

Funding

The authors have no relevant financial or non-financial interests to disclose.

Author information

Authors and Affiliations

Authors

Contributions

All authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by Isabella Morlini, Sean Albertson and Stefano Orlandini. The first draft of the manuscript was written by Isabella Morlini and Stefano Orlandini and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Isabella Morlini.

Ethics declarations

Conflict of interest

The authors declare no conflicts of interest relevant to this study.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Appendix 1

Appendix 1

Considering annual values of the NAO teleconnection index, from 1950 to 2019, the statistical assessment of the presence or absence of a monotonic trend, inhomogeneity and change points, has been made by performing several nonparametric and parametric tests, using XLSTAT (Addinsoft 2023). P values of the different tests are reported in Table 6.

Table 6 P values of different homogeneity tests

At the 1% level, all tests agree in accepting the null hypothesis that the annual values are independent, and the series is homogeneous. Correlations with the other climatic series of P and T are then likely to be nonspurious. Correlation coefficients and p values are reported in Table 7.

Table 7 Correlation coefficients and related p values (in brackets)

The correlation between the NAO index and precipitation P, and the correlation between the NAO index and temperature T are clearly not significant. To further investigate the influence of the NAO index on the local temperature T and compare the nonlinear relationship between T and NAO with the nonlinear relationship between T and C and T and A, different smoothing splines and the 4-parameter logistic function are estimated. In Fig. 8, cubic smoothing splines with S = 1.0 and S = 0.2 and the 4-parameter logistic function LF are reported. Table 8 reports the NSE error (2) for smoothing splines. Clearly, smoothing splines with S = 1.0 overfit the data, while all the other estimated smoothing splines have similar AIC and BIC values (not reported) and can be selected as models with a proper degree of smoothness.

Fig. 8
figure 8

Time series of T and A (plot a), T and C (plot b), T and NAO (plot c), estimated values of T considering A (plot d) and C (plot e) and NAO (plot f) as dependent variable using smoothing splines (SS) with parameters 1.00 and 0.20 and the 4-parameter logistic function (LF)

Table 8 NSE errors

Both the graphical representation of the estimated functions and the values of the NSE indicate that local urbanization A and the global CO2 concentration C are more related to local temperature T than the NAO index and are better predictors of temperature changes in the last 70 years. In conclusion, results suggest that neither local precipitation P nor the NAO index seem to be related to the local temperature T. This may be because local geographical features may modify the impact of the NAO index on precipitation and temperature in specific areas.

Rights and permissions

Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Morlini, I., Albertson, S. & Orlandini, S. Characterization of annual urban air temperature changes with special reference to the city of Modena: a comparison between regression models and a proposal for a new index to evaluate relationships between environmental variables. Stoch Environ Res Risk Assess 38, 1163–1178 (2024). https://doi.org/10.1007/s00477-023-02622-x

Download citation

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s00477-023-02622-x

Keywords

Navigation