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.
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
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
Boccolari MA, Malmusi S (2013) Changes in temperature and precipitation extremes observed in Modena, Italy. Atmos Res 122:16–31
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
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
Buishand TA (1982) Some methods for testing the homogeneity of rainfall records. J Hydrol 58:11–27
Buishand TA (1984) Tests for detecting a shift in the mean of hydrological time series. J Hydrol 73:51–69
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
Cai Z, Liu Q, Cao S (2020) Real estate supports rapid development of China’s urbanization. Land Use Policy 95:104582
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
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
De Boor C (2001) A practical guide to splines. Springer
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
Fawcett L, Walshaw D (2012) Estimating return levels from serially dependent extremes. Environmetrics 23:272–283
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
Findley DF (1991) Counterexamples to parsimony and BIC. Ann Inst Stat Math 43(3):505–514. https://doi.org/10.1007/BF00053369
Friedman J (2001) Greedy function approximation: a gradient boosting machine. Ann Stat 29:1189–1232
Friedman J (2002) Stochastic gradient boosting. Comput Stat Data Anal 38:367–378
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
Gilbert RO (1987) Statistical methods for environmental pollution monitoring. Wiley, New York
Granger CWJ (1969) Investigating causal relations by econometric models and cross-spectral methods. Econometrica 37:424–438
Granger CWJ, Ghysels E, Swanson NR, Watson MW (2001) Essays in econometrics: collected papers of Clive W. J. Granger. Cambridge University Press
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
Hastie T, Tibshirani R (1986) Generalized additive models (with discussion). Stat Sci 1:297–318
Hastie T, Tibshirani R (1990) Generalized additive models. Chapman & Hall/CRC
Hawkins M (1977) Testing a sequence of observations for a shift in location. J Am Stat Assoc 72:180–186
Hipel KW, McLeod AI (1994) Time series modelling of water resources and environmental systems. Elsevier Science, New York
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
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
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
Kalnay E, Kay M (2003) Impact of urbanization and land-use change on climate. Nature 423:528–531. https://doi.org/10.1038/nature01675
Karl TR, Diaz HF, Kukla G (1988) Urbanization: its detection and effect in the United States climate record. J Clim 1(11):1099–1123
Kendal MG (1975) Rank correlation methods, 4th edn. Charles Griffin, London
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
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
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
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
Mann HB (1945) Non-parametric tests against trend. Econometrica 13:163–171
McQuarrie ADR, Tsai CL (1998) Regression and time series model selection. World Scientific
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
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
Morlini I, Franco-Villoria M, Orlandini S (2023) Modelling local climate change using site-based data. Environ Ecol Stat 30(2):205–232
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
Pettitt AN (1979) A non-parametric approach to the change point problem. J R Stat Soc: Ser C Appl Stat 28:126–135
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
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
Reinsch CH (1967) Smoothing by Spline Functions. Numer Math 10(3):177–183. https://doi.org/10.1007/BF02162161
Sen PK (1968) Estimates of the regression coefficient based on Kendall’s tau. J Am Stat Assoc 63:1379–1389
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
Schwarz GE (1978) Estimating the dimension of a model. Ann Stat 6(2):461–464. https://doi.org/10.1214/aos/1176344136
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
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
Von Neumann J (1941) Distribution of the ratio of the mean square successive difference to the variance. Ann Math Stat 12:367–395
Wahba G (1990) Spline models for observational data. SIAM, Philadelphia
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
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
Wilks SS (1932) Certain generalizations in the analysis of variance. Biometrika 24:471–494
Winjngaard JB, Klein Tank AMG, Konnen GP (2003) Homogeneity of 20th century European daily temperature and precipitation series. Int J Climatol 23:679–692
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
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
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
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
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.
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.
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.
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.
About this article
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
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00477-023-02622-x