4.1. ARDL Analysis
This section presents the findings of the empirical calculations derived from the first proposed method. To begin with, we performed a cross-sectional dependence test as shown in
Table 2, which revealed that all three models employed—namely Pesaran, Frees, and Friedman—strongly rejected the hypothesis of cross-sectional independence, with statistical significance at the 1% level. This outcome indicates that the residual correlation value of the proposed model demonstrates a clear cross-sectional dependence under the fixed effect specification.
In order to further evaluate the validity of our proposed model, the current study conducted a stationary test, utilizing both the panel root test developed by Choi [
50] and the IPS test proposed by Im, et al. [
51]. As shown in
Table 3, the results of the stationary test reveal that all variables employed in our analysis strongly reject the null hypothesis at a significance level of 1%.
The subsequent analysis involved the cointegration test panel aimed at establishing the long-term association among variables. To achieve this objective, the study employed a set of models developed by Kao [
47], Westerlund [
48], Westerlund [
49], Pedroni [
45], and Pedroni [
46]. As presented in
Table 4, the findings indicate that all the statistical tests used achieved a significant level of 1%. In this regard, it is evident that all the models rejected the null hypothesis, which postulates that the data series were not related in the long term. To determine the optimal lag lengths, this test used three different criteria: Bayesian, Akaike, and Quinn. By applying these widely accepted criteria, we were able to select the most appropriate lag lengths for our model, providing a more accurate and comprehensive analysis of the underlying patterns and dynamics in our data. Hence, the results provide empirical evidence that the data in the current research exhibit a cointegrated behavior. The outcomes of this study, therefore, indicate that the proposed models are suitable for testing the cointegration of panel data.
After conducting thorough examinations, the findings suggest that the data used in this research are appropriate for analysis utilizing the panel autoregressive distributed lag (ARDL) model. The estimation results of the DTR effect are presented in
Table 5. The estimation results in
Table 5 shows two distinguished estimators used in this study, MG and PMG.
In this research, the pooled mean group (PMG) and mean group (MG) estimators were compared to obtain the best estimation results. The Hausman test was used to determine the best and most efficient model, which was found to be the MG estimator. According to the MG calculations, a 1 °C decrease in the DTR over the long term was associated with a 17.1% increase in the risk of an RTA. Interestingly, the effect of DTR on RTAs was observed to be smaller in the short term, as the MG model’s estimation indicated a 4% increase in RTAs with a decrease in DTR. In comparison to the effect of air pollution, specifically particulate matter, the magnitude of the DTR effect on RTAs was found to be smaller. However, it should be noted that estimating results using air pollution variables such as PM
10 can lead to bias if calculated directly [
3]. Therefore, when compared to estimates from other studies utilizing air pollution variables, such as the study by Sager [
3], the effect of DTR on RTAs was observed to be larger, particularly in Taiwan. These findings suggest that DTR may play a significant role in the incidence of RTAs in certain regions and should be taken into consideration when implementing measures to reduce the incidence of these accidents.
The present study’s findings corroborate previous research that has reported a positive correlation between air pollution and traffic accidents [
3,
4,
52]. Specifically, the estimation outcomes derived from the ARDL model reveal that PM
10 exert a positive influence on road traffic accidents. This phenomenon is observed across both short-term and long-term estimates. In fact, long-term estimates indicate that PM
10 possess a significant effect on road traffic accidents. Evidently, over the long term, the deleterious effects of air pollution are set to be exacerbated. It is worth noting that previous studies have uncovered a close association between air pollution and human cognitive performance [
53].
4.2. DLNM Analysis
The subsequent phase of this investigation delves into the non-linear association between DTR-lag-road traffic accident via the utilization of distributed lag non-linear models (DLNMs). The examination encompasses six major cities in Taiwan, specifically Taipei, New Taipei, Taoyuan, Taichung, Tainan, and Kaohsiung.
The lag time plays a crucial role in understanding the temporal dynamics of the relationship between DTR and road traffic accidents. In previous research, lag effects of DTR on health outcomes have often been observed in the order of days [
11,
12]. These shorter lag periods are appropriate when studying the immediate health impacts of temperature fluctuations on individuals. However, in this study, the lag effects in the order of months is driven by the recognition that certain consequences of DTR variations may have a delayed and cumulative nature, akin to the lag effects observed in relation to broader climate phenomena. We extended the lag time to encompass several months to capture the potential lagged and cumulative effects of DTR on road traffic accidents. This approach is deemed necessary because while DTR and driving are indeed experienced on a shorter timescale, the influence of DTR on road safety may manifest over more extended periods due to factors such as gradual changes in road conditions, driver behavior adaptations, and the cumulative impact of temperature variations on road infrastructure.
Figure 2 and
Figure 3 illustrate the non-linear and lagged relationship between the diurnal temperature range (DTR) and road traffic accidents in six major metropolitan areas in Taiwan. The purpose of presenting two types of figures is to help readers better understand the findings in cases where one figure may not be clear. These figures reveal how the urban areas in Taiwan respond differently to changes in DTR with regards to road traffic accidents. The results indicate that a lower DTR effect tends to increase road traffic accidents in New Taipei, particularly at lags of 0–10 months. However, at higher DTRs (≥4.6 °C), the relationship between DTR and road traffic accidents is relatively lower at lags of 0–25 months. This trend becomes more prominent at higher DTR levels, where it reinforces the association between DTR and road traffic accidents. The maximum lag of 30 months in these figures provides a comprehensive understanding of the temporal relationship between DTR and road traffic accidents.
The results for Taipei are comparable to those obtained in New Taipei. A low diurnal temperature range has a higher association with road traffic accidents in Taipei, as revealed in
Figure 2 and more comprehensively in
Figure 3. The relative risk (RR) value is ≥1 at the DTR level of 2–2.5 °C and lag of 0–15 months, indicating an elevated risk of road traffic accidents. As the DTR value increases, the condition tends to improve, as evidenced by the lower RR value. The pattern observed in Taoyuan is quite similar to that observed in Taichung. At low DTR levels, road traffic accidents are moderate, and there is a positive correlation with an increase in DTR, particularly at lag 0–5 months. However, the relationship between road traffic accidents and DTR tends to improve at lag ≥ 10 month, as demonstrated by lower RR values.
Additionally, the study found that Tainan exhibits a unique pattern regarding the relationship between diurnal temperature range (DTR) and road traffic accidents (RTA). At a DTR level of 3.5 °C, the RTA response tends to be high at lag 0, indicating an immediate effect. However, conditions improve at lags of 1–30 months, indicating that the low level of DTR is not accompanied by an increase in road traffic accidents at lag ≥ 1 month. When the DTR level reaches 4.0–5.5 °C, the response to RTAs is high. On the other hand, the situation in Tainan is nearly identical to that in Kaohsiung, except for the fact that in Kaohsiung, a low DTR level is accompanied by a high RR rate of road traffic accidents, particularly at lags of 0–25 months with a DTR level <3.5 °C. This implies that the relationship between DTR and road traffic accidents in these two cities may be influenced by similar factors. However, the magnitude and timing of the effect vary, indicating that local factors may also play a role in shaping this relationship.
The present analysis investigates the response of road traffic accidents to specific DTR and lag values. As highlighted by Gasparrini, et al. [
35] and Gasparrini [
54], explaining the RTA response without specific DTR and lag values can be weak and is considered a limitation in the analysis shown in
Figure 2 and
Figure 3. The inadequacy of the figures to explicate the impact of diurnal temperature range (DTR) on particular quantitative parameters, coupled with their constraining nature with regards to inferential objectives, is the root cause of this phenomenon. The insufficiency of these aforementioned visual representations to elucidate the nuanced effects of DTR on discrete variables and their incapacity to cater to inferences and conjectures concerning said variables are what underlie the observed trend.
In an attempt to circumvent this aforementioned constraint,
Figure 4 and
Figure 5 have been introduced to portray the response of road traffic accidents at specific DTR and lag values, proffering an unparalleled advantage in comprehending the impact of DTR on RTAs, particularly where it is most pronounced. This critical analysis presents a unique opportunity to either validate or invalidate the findings of prior analyses, including those employing panel autoregressive distributed lag (ARDL) in all cities and counties in Taiwan. As such, the present inquiry endeavors to provide a more comprehensive and precise comprehension of the interplay between DTR, lag, and fatal traffic accidents, the outcome of which can have far-reaching implications for devising efficacious road safety strategies in Taiwan.
Upon evaluating the computational outcomes demonstrated in the graphical representation exhibited in
Figure 4, it is discernible that divergent urban regions across Taiwan display unique responses. It is important to note that the DTR values depicted in
Figure 4 have been adjusted to correspond to the DTR values for each respective city. As evinced by the aforementioned figure, a significant surge in the response to DTR levels from road traffic accidents occurred at the lag of 0–5 months in three metropolises, namely New Taipei, Taipei, and Taoyuan, while the escalation in Taichung, Tainan, and Kaohsiung was comparatively less pronounced. It is worth highlighting that with regard to road traffic accidents, the response exhibited to low and high DTR values is distinct. In the case of Taoyuan, high DTR rates were accompanied by an elevated number of road traffic accidents. On the contrary, in the cities of New Taipei, Taipei, and Taichung, a low DTR rate was found to be associated with a considerable surge in road traffic accidents. The responses in Tainan and Kaohsiung were similar, with moderate impacts of both low and high DTR levels in both cities. It is noteworthy that this study posits that with the exception of Taoyuan, low DTR rates were linked to a high response from rates of road traffic accidents across all the cities.
The estimation continues by scrutinizing the DTR-RTA response at a particular lag, as evidenced in
Figure 5. Based on the evidence presented in this figure, it can be deduced that employing a specified lag of 5 month means that the response elicited from road traffic accidents at a low DTR level appears to be substantially amplified in virtually all cities, except for New Taipei. In addition,
Figure 5 illustrates that in the context of Tainan, when the DTR level is low, the response observed from road traffic accidents is relatively subdued and tends to augment as the DTR level rises. Nevertheless, it is noteworthy that this trend undergoes a shift when the DTR level reaches 4 °C. At this threshold, the aforementioned condition is noted as ameliorating.
The variation in the effect of DTR on road traffic accidents between different cities with specific DTR values can be attributed to several factors. First, each city has its own unique local weather patterns influenced by its geographical location and proximity to bodies of water. Second, cities situated at higher elevations often encounter more drastic temperature swings between day and night. This can lead to increased expansion and contraction of road surfaces, contributing to DTR-related road hazards. Third, local precipitation patterns can interact with DTR to affect road conditions. In summary, differences in local weather, geographical characteristics, and climate-related factors can lead to varying effects of diurnal temperature range (DTR) on road traffic accidents (RTAs) between different cities. These factors influence the intensity and frequency of DTR, contributing to the diversity of road safety outcomes across urban areas.