“…Recently, an increasing number of studies have used non-linear time series analysis methods, such as GAM [ 5 , 23 , [38] , [39] , [40] , [41] , [42] , [43] , [44] ] and DLNM [ 28 , [44] , [45] , [46] ]. However, these studies encompass analyses at the country [ 43 , 44 ], state [ 5 , 26 , 40 ] or city [ 23 , 28 , 38 , 39 , 41 , 45 , 46 ] level. In studies with a country or state as the unit of analysis, the regional scope is relatively large, and acquiring data on meteorological indicators and their representativeness is prone to significant errors.…”