Investigating Drought Events and Their Consequences in Wildﬁres: An Application in China

: Understanding the impact of drought on ﬁre dynamics is crucial for assessing the potential effects of climate change on wildﬁre activity in China. In this study, we present a series of multiple linear regression (MLR) models linking burned area (BA) during mainland China’s ﬁre season from 2001 to 2019, across seven regions, to concurrent drought, antecedent drought, and time trend. We estimated burned area using Collection 6 Moderate Resolution Imaging Spectradiometer (MODIS) and drought indicators using either the Standardized Precipitation Evapotranspiration Index (SPEI) or the self-calibrated Palmer Drought Severity Index (sc-PDSI). Our ﬁndings indicate that the wildﬁre season displays a spatial variation pattern that increases with latitude, with the Northeast China (NEC), North China (NC), and Central China (CC) regions identiﬁed as the primary areas of wildﬁre occurrence. Concurrent and antecedent drought conditions were found to have varying effects across regions, with concurrent drought as the dominant predictor for NEC and Southeast China (SEC) regions and antecedent drought as the key predictor for most regions. We also found that the Northwest China (NWC) and CC regions exhibit a gradual decrease in burned area over time, while the NEC region showed a slight increase. Our multiple linear regression models exhibited a notable level of predictive power, as evidenced by the average correlation coefﬁcient of 0.63 between the leave-one-out cross-validation predictions and observed values. In particular, the NEC, NWC, and CC regions demonstrated strong correlations of 0.88, 0.80, and 0.76, respectively. This indicates the potential of our models to contribute to the prediction of future wildﬁre occurrences and the development of effective wildﬁre management and prevention strategies. Nevertheless, the intricate relationship among ﬁre, climate change, human activities, and vegetation distribution may limit the generalizability of these ﬁndings to other conditions. Consequently, future research should consider a broad range of factors to develop more comprehensive models.


Introduction
Wildfire plays a significant part in the terrestrial carbon cycle and ecosystem functioning as a natural disturbance process that is prevalent across the majority of the earth's surface [1,2]. The frequency and severity of wildfires have grown globally in recent decades as result of climate change [3,4], resulting in significant economic and human losses as well as substantial environmental pollution [5,6]. Understanding how wildfires respond to climate change is important for monitoring, predicting, and developing wildfire prevention and control programs.

Study Area
The Chinese mainland is geographically divided into seven distinct regions: Northeast China (NEC), Northwest China (NWC), North China (NC), Southwest China (SWC), Central China (CC), South China (SC), and Southeast China (SEC) [23,24]. These regions were chosen as the spatial units of analysis for several reasons. Adopting a region-based division approach offers the potential to effectively disseminate the findings of this study to national-level policymakers and personnel involved in forest-fire management. In terms of elevation, NEC, CC, and SEC generally have elevations below 500 m ( Figure 1a). SWC, western SC, and northern NWC comprise mountainous and plateau areas, with elevations reaching up to 3000 m (excluding the lower Sichuan Basin in northeast SC). The major vegetation types in the study area include needle and broadleaf forests in NEC, SEC, SC, and southern CC, while northern NC is predominantly covered by grasslands.
policymakers and civil protection agencies in improving early warning systems and implementing more effective fire management strategies. This study aims to investigate the suitability of a simple drought-fire model for modeling annual changes in burned area (BA) in China and to identify the drought index time window that exerts a more significant influence on BA.

Study Area
The Chinese mainland is geographically divided into seven distinct regions: Northeast China (NEC), Northwest China (NWC), North China (NC), Southwest China (SWC), Central China (CC), South China (SC), and Southeast China (SEC) [23,24]. These regions were chosen as the spatial units of analysis for several reasons. Adopting a region-based division approach offers the potential to effectively disseminate the findings of this study to national-level policymakers and personnel involved in forest-fire management. In terms of elevation, NEC, CC, and SEC generally have elevations below 500 m ( Figure 1a). SWC, western SC, and northern NWC comprise mountainous and plateau areas, with elevations reaching up to 3000 m (excluding the lower Sichuan Basin in northeast SC). The major vegetation types in the study area include needle and broadleaf forests in NEC, SEC, SC, and southern CC, while northern NC is predominantly covered by grasslands. The precipitation distribution in China follows a general pattern of decreasing from southeast to northwest, with the highest concentration of precipitation occurring between April and October. The regions of North China (NC) and Northwest China (NWC) receive the least amount of precipitation, with the highest monthly precipitation being less than 100 mm. Conversely, Southeast China (SEC), Southwest China (SWC), and Central China (CC) experience abundant precipitation, primarily influenced by the subtropical monsoon climate ( Figure 2). In terms of temperature distribution, China exhibits high temperatures during summer and notable variations between the northern and southern regions in winter. During January-March and November-December, the temperature in Northeast China (NEC), Northwest China (NWC), North China (NC), and Southwest China (SWC) The precipitation distribution in China follows a general pattern of decreasing from southeast to northwest, with the highest concentration of precipitation occurring between April and October. The regions of North China (NC) and Northwest China (NWC) receive the least amount of precipitation, with the highest monthly precipitation being less than 100 mm. Conversely, Southeast China (SEC), Southwest China (SWC), and Central China (CC) experience abundant precipitation, primarily influenced by the subtropical monsoon climate ( Figure 2). In terms of temperature distribution, China exhibits high temperatures during summer and notable variations between the northern and southern regions in winter. During January-March and November-December, the temperature in Northeast China (NEC), Northwest China (NWC), North China (NC), and Southwest China (SWC) remains below 0 • C. On the other hand, South China (SC), Central China (CC), and Southeast China (SEC) maintain temperatures above 0 • C throughout the year ( Figure 2).

Data
Burned area (BA) data were derived from the collection 6 Moderate Resolution Imaging Spectradiometer (MODIS) Burned Area (MCD64MCQ). The data are composed of monthly BA from 2001 to 2019 since the numerical values were lacking in the summer of 2020 due to a satellite sensor issue. The MCD64MCQ data product is convenient for regional modelling under continuous roughly spatiotemporal scale with spatial resolution of 0.25. The accuracy of the data was validated by numerous independent investigations using the high-resolution Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) image. Detailed information about this issue can be found elsewhere [25][26][27][28].

Data
Burned area (BA) data were derived from the collection 6 Moderate Resolution Imaging Spectradiometer (MODIS) Burned Area (MCD64MCQ). The data are composed of monthly BA from 2001 to 2019 since the numerical values were lacking in the summer of 2020 due to a satellite sensor issue. The MCD64MCQ data product is convenient for regional modelling under continuous roughly spatiotemporal scale with spatial resolution of 0.25. The accuracy of the data was validated by numerous independent investigations using the high-resolution Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) image. Detailed information about this issue can be found elsewhere [25][26][27][28].
Elevation, land cover, temperature, and precipitation data play a crucial role in providing the hydrological and topographical context for the study area. The elevation and land-cover data were sourced from the International Steering Committee for Global Mapping (ISCGM) through their website at https://globalmaps.github.io/ (accessed on 28 October 2022). Temperature and precipitation data were derived from the gridded Climatic Research Unit (CRU) Time-series (TS) dataset version 4.05 [29], which is also used for calculating the SPEI and sc-PDSI datasets by the University of East Anglia Climatic Research Unit. For assessing drought conditions, we utilized the Standardized Precipitation Evapotranspiration Index (SPEI) [19] and the self-calibrated Palmer Drought Severity Index (sc-PDSI) [17]. The sc-PDSI, similar to the PDSI, incorporates autoregression, where the sc-PDSI value for a specific month is determined by the weighted sum of the previous month's sc-PDSI and the moisture anomaly index. The moisture anomaly index is computed using the moisture departure (the difference between actual precipitation and climatically appropriate precipitation) and climatic characteristics. Palmer [18] derived empirical constants for climatic characteristics and duration factors based on data from several U.S. stations, while Wells [17] introduced a method to automatically adjust these empirical constants for a specific location using historical climate data. Elevation, land cover, temperature, and precipitation data play a crucial role in providing the hydrological and topographical context for the study area. The elevation and land-cover data were sourced from the International Steering Committee for Global Mapping (ISCGM) through their website at https://globalmaps.github.io/ (accessed on 28 October 2022). Temperature and precipitation data were derived from the gridded Climatic Research Unit (CRU) Time-series (TS) dataset version 4.05 [29], which is also used for calculating the SPEI and sc-PDSI datasets by the University of East Anglia Climatic Research Unit. For assessing drought conditions, we utilized the Standardized Precipitation Evapotranspiration Index (SPEI) [19] and the self-calibrated Palmer Drought Severity Index (sc-PDSI) [17]. The sc-PDSI, similar to the PDSI, incorporates autoregression, where the sc-PDSI value for a specific month is determined by the weighted sum of the previous month's sc-PDSI and the moisture anomaly index. The moisture anomaly index is computed using the moisture departure (the difference between actual precipitation and climatically appropriate precipitation) and climatic characteristics. Palmer [18] derived empirical constants for climatic characteristics and duration factors based on data from several U.S. stations, while Wells [17] introduced a method to automatically adjust these empirical constants for a specific location using historical climate data.
The SPEI, on the other hand, combines precipitation and temperature data, offering the advantage of capturing temperature variability's influence on drought assessment across multiple time scales [17]. The calculation procedure involves a climatic water balance, accumulation of deficit/surplus at different time scales, and adjustment to a log-logistic probability distribution [17]. We obtained the SPEI (version 2.7) and sc-PDSI datasets from the University of East Anglia Climatic Research Unit (UEACRU), with a spatial resolution of 0.5 degrees [30,31]. Overall, the combination of these datasets enables a comprehensive analysis of the hydrological and climatic conditions, which is essential for understanding drought dynamics in the study area.

Methodology
Many previous studies have attempted to investigate the effects of climate on wildfires at the original spatial scales of the dataset [32][33][34]. However, the process by which climate affects wildfire is very complex, and the smaller spatial scale is very sensitive to human factors, vegetation factors, and other arbitrary factors related to ignition. Therefore, the Chinese mainland is divided into seven distinct regions.
We used the SPEI and sc-PDSI indices, respectively, as drought indicators, and we processed the data as follows: We normalized the positively biased BA variables by applying a log transformation. The drought indicator is quantified by computing the mean of the indices over a duration of several months and expressed as Drought c,m,b . Here, m denotes the month used for drought indicator computation (fireseason start and fireseason end ), c represents the number of months, and b signifies the buffer period prior to the onset of the fire season. For concurrent drought, c is multiple months in the fire season (fireseason length ), m is the month when the fire season ends, and b is 0. For antecedent drought, c ranges from 1 to 12 months, m is the month when the fire season starts, and b ranges from 2 to 14 months.
For example, Drought 3,6,2 represents the average value from January to March.
Multiple linear regression (MLR) is commonly employed for establishing interannual climate-fire models owing to its ease of interpretation and the avoidance of overfitting associated with the use of non-linear models. For each region in China, we used the following drought-BA model to express the potential relationship between the annual variation of fire season BA and drought indices (SPEI, sc-PDSI), where the number of predictors is limited to a maximum of three (concurrent drought indicator, antecedent drought indicator, and time term): where BA(i, t) is the predicted BA in the i region and year t. β 0 is the intercept; β 1 represents the sensitivity of BA in each region to concurrent dry conditions as summarized by the drought index; β 2 represents the sensitivity of BA in each region to antecedent dry conditions as summarized by the drought index; β 3 is the coefficient of the time term T (in years) that characterizes the temporal trends of BA, thus taking into account the possible influence of slowly changing factors over the study period; ε is a stochastic noise term that captures all other factors that influence BA except for drought index. For each region, we considered these combinations of predictors: concurrent drought only, antecedent drought only, concurrent drought plus time, antecedent drought plus time, and concurrent and antecedent drought plus time. We employed the Akaike information criterion (AIC) [39] to select the models with the lowest AIC values. The AIC allows us to assess the quality of each model relative to the others in a given collection of models. The calculation method for the AIC value of each multiple linear regression (MLR) model is as follows: Here, L represents the maximum value of the likelihood function for the model, and k represents the number of parameters to be estimated in the model. The AIC combines the goodness of fit, as measured by the likelihood function, with a penalty that increases as the number of estimated parameters increases. This ensures that the AIC considers both model fit and complexity. To further evaluate the selected model, we performed leave-one-out cross-validation. This validation technique involves systematically leaving out one data point at a time, fitting the model to the remaining data, and then assessing the model's predictive performance using the left-out data point. This process is repeated for each data point, and the overall performance is evaluated. For a step-by-step description of the model development, selection, and evaluation process, please refer to Figure 3, which provides a visual representation of the procedure. models. The calculation method for the AIC value of each multiple linear regression (MLR) model is as follows:

AIC 2 2ln
Here, L represents the maximum value of the likelihood function for the model, and k represents the number of parameters to be estimated in the model. The AIC combines the goodness of fit, as measured by the likelihood function, with a penalty that increases as the number of estimated parameters increases. This ensures that the AIC considers both model fit and complexity. To further evaluate the selected model, we performed leaveone-out cross-validation. This validation technique involves systematically leaving out one data point at a time, fitting the model to the remaining data, and then assessing the model's predictive performance using the left-out data point. This process is repeated for each data point, and the overall performance is evaluated. For a step-by-step description of the model development, selection, and evaluation process, please refer to Figure 3, which provides a visual representation of the procedure. The statistical analysis and data processing mentioned above were carried out using Python. Specifically, the fitting and prediction of the MLR model were accomplished using the Scipy version 1.10.0, while the visualization was achieved using the Geopandas version 0.12.2 and related dependencies.

Results
During the period of 2001-2019, NEC, NC, and CC were the main regions where wildfires occurred in mainland China, accounting for approximately 28.49%, 19.81%, and 34.55% of the total burned area (BA), respectively. SC, SEC, NWC, and SWC accounted for about 7.16%, 6.35%, 3.23%, and 4.12%, respectively. The variations of the fire season The statistical analysis and data processing mentioned above were carried out using Python. Specifically, the fitting and prediction of the MLR model were accomplished using the Scipy version 1.10.0, while the visualization was achieved using the Geopandas version 0.12.2 and related dependencies.

Results
During the period of 2001-2019, NEC, NC, and CC were the main regions where wildfires occurred in mainland China, accounting for approximately 28.49%, 19.81%, and 34.55% of the total burned area (BA), respectively. SC, SEC, NWC, and SWC accounted for about 7.16%, 6.35%, 3.23%, and 4.12%, respectively. The variations of the fire season showed certain spatial distribution patterns. Most southern regions (SWC, SC, and CC) had their fire seasons concentrated from January to June (1-6), with a shorter average duration of the fire season (3 months; see Figure 4. In contrast, northern regions (NWC, NC, and NEC) had their fire season concentrated from February to October (2-10), with a longer average duration of the fire season (7.33 months; see Figure 4. SEC is the region with the lowest latitude, and its fire season lasts almost all year round (from January to November; see Figure 4. To capture the impacts of concurrent and antecedent droughts, we calculated the average drought index (SPEI/sc-PDSI) during the fire season. Additionally, we set a buffer of 2-14 months before the start of the fire season to differentiate between concurrent and antecedent droughts. For each buffer, we calculated the average for 1-12 months. Therefore, there are 156 antecedent drought indicators and 1 concurrent drought indicator and 1 time variable for each region. Finally, we considered these combinations of predictors: concurrent drought only, antecedent drought only, concurrent drought plus time, antecedent Fire 2023, 6, 223 7 of 12 drought plus time, and concurrent and antecedent drought plus time. Overall, we fitted 470 models (1 + 156 + 1 + 156 + 156) for each region and selected the model based on the AIC criterion. The identified key variables can be considered as potential predictors for each region. The concurrent drought is a potential predictor for NEC and SEC (Figure 5a), while the antecedent drought is a potential predictor for all regions except SEC (Figure 5b). In addition, the time variable is a potential predictor for NWC, CC, and NEC (Figure 5c). Figure 2a-c show the coefficients for concurrent drought, antecedent drought, and time, respectively. These coefficients indicate the sensitivity of BA to the predictors. Except for the antecedent drought coefficient for NEC, all drought coefficients are negative. The average coefficient for concurrent drought is −0.65, and the average coefficient for antecedent drought is −2.73 (for all regions, Table 1). The models show that, except for NEC, as the drought indicators decreases (i.e., drier climate conditions), BA increases. The model for NEC shows that if the concurrent drought indicator decreases (i.e., drier climate conditions), BA increases accordingly. However, if the previous winter was drier (antecedent drought indicator for NEC from September to December of the previous year), BA decreases. The time coefficient indicates that from 2001 to 2019, there was a slow declining trend in BA in the NWC and CC regions (Table 1), while the BA in the NEC region showed a slow increasing trend (Figure 5c).  To capture the impacts of concurrent and antecedent droughts, we calculated the average drought index (SPEI/sc-PDSI) during the fire season. Additionally, we set a buffer of 2-14 months before the start of the fire season to differentiate between concurrent and antecedent droughts. For each buffer, we calculated the average for 1-12 months. Therefore, there are 156 antecedent drought indicators and 1 concurrent drought indicator and 1 time variable for each region. Finally, we considered these combinations of predictors: concurrent drought only, antecedent drought only, concurrent drought plus time, antecedent drought plus time, and concurrent and antecedent drought plus time. Overall, we fitted 470 models (1 + 156 + 1 + 156 + 156) for each region and selected the model based on the AIC criterion. The identified key variables can be considered as potential predictors for each region. The concurrent drought is a potential predictor for NEC and SEC ( Figure  5a), while the antecedent drought is a potential predictor for all regions except SEC (Figure 5b). In addition, the time variable is a potential predictor for NWC, CC, and NEC (Figure 5c). Figure 2a-c show the coefficients for concurrent drought, antecedent drought, and time, respectively. These coefficients indicate the sensitivity of BA to the predictors. Except for the antecedent drought coefficient for NEC, all drought coefficients are negative. The average coefficient for concurrent drought is −0.65, and the average coefficient for antecedent drought is −2.73 (for all regions, Table 1). The models show that, except for   We used parsimonious models, as mentioned before, which are models that explain phenomena or predict results using as few parameters or variables as possible. Figure 6a shows the correlation coefficient (RhoIN) between the in-sample model fit values and the observed values. (The precise definition of each regional model is shown in Table 1.) The average of the correlations for all models is 0.73. Additionally, an important test for these models is to verify their ability to predict out-of-sample data (obtained through leave-oneout cross-validation). Figure 6b shows the correlation coefficient (RhoOUT) between the out-of-sample model fit values and the observed values. The average of the correlations for all models is 0.63. This indicates that the models have relatively good predictive power. We used parsimonious models, as mentioned before, which are models that explain phenomena or predict results using as few parameters or variables as possible. Figure 6a shows the correlation coefficient (RhoIN) between the in-sample model fit values and the observed values. (The precise definition of each regional model is shown in Table 1.) The average of the correlations for all models is 0.73. Additionally, an important test for these models is to verify their ability to predict out-of-sample data (obtained through leave-oneout cross-validation). Figure 6b shows the correlation coefficient (RhoOUT) between the out-of-sample model fit values and the observed values. The average of the correlations for all models is 0.63. This indicates that the models have relatively good predictive power.

Discussion
Climate plays a significant role in influencing the likelihood, intensity, and frequency of wildfires. Changes in climate patterns, such as temperature and precipitation, can create conditions that increase the risk of wildfires [40]. In addition, human activities such as land-use change and fire-suppression policies can exacerbate the effects of climate change on wildfire [41,42]. Due to China's enormous population, diverse climatic conditions, and vegetation types, there are considerable regional variations in the spatial and temporal distribution characteristics of wildfires. Nevertheless, we found that the fire season shows a spatial variation pattern that increases with latitude. In addition, NEC, NC, and CC were the main regions where wildfires occurred in mainland China. The spatial distribution of burned area (BA) is consistent with the results of many independent studies [24,43,44] but not consistent with the spatial distribution of fire frequency [45]. These results indicate the

Discussion
Climate plays a significant role in influencing the likelihood, intensity, and frequency of wildfires. Changes in climate patterns, such as temperature and precipitation, can create conditions that increase the risk of wildfires [40]. In addition, human activities such as land-use change and fire-suppression policies can exacerbate the effects of climate change on wildfire [41,42]. Due to China's enormous population, diverse climatic conditions, and vegetation types, there are considerable regional variations in the spatial and temporal distribution characteristics of wildfires. Nevertheless, we found that the fire season shows a spatial variation pattern that increases with latitude. In addition, NEC, NC, and CC were the main regions where wildfires occurred in mainland China. The spatial distribution of burned area (BA) is consistent with the results of many independent studies [24,43,44] but not consistent with the spatial distribution of fire frequency [45]. These results indicate the BA in the north of China may be larger than that in the south, but the frequency is lower.
In this study, we present a series of multiple linear regression (MLR) models linking burned area (BA) during mainland China's fire season from 2001 to 2019 and across seven regions to concurrent drought, antecedent drought, and time trend. Specifically, we investigated the dependency of wildfires during the fire season on concurrent drought conditions, which affect the dryness of fuel through climatic factors, as well as antecedent drought conditions, which impact the fuel load and the fine fuel structure through climate factors. It is worth noting that most regions included antecedent drought in their models, while only two regions included concurrent drought. This result differs significantly from that of the Mediterranean region [13,46]. Possible reasons for this discrepancy include the following: (1) The measurement of concurrent drought was solely based on the average values during the fire season, while the time window for antecedent drought was more extensive, and (2) concurrent drought may be highly correlated with early-stage drought and was thus excluded during model selection. (Studies have shown a clear drought trend throughout China [47].) In the case of China, Zhao [16] conducted a comparison of the fitting capabilities of different climate factors using a univariate linear regression approach. The study categorized the observed series into three groups: above normal, normal, and below normal. The assessment of the model's fitting capability was based on the count of observations and fitted values that fell into the above-normal category simultaneously. In our research, we rigorously evaluated the model's fitting capability and achieved higher rates of fit in most regions compared to the predictive fitting rates reported in Zhao's study.
The specific reasons warrant further investigation. The study is a necessary step in addressing the question of whether climate change will increase wildfire activity in China. The impact of climate change on wildfire activity is complex, as an increase in fire weather can be exacerbated or offset by changes in vegetation productivity [9].
The statistical analysis reported here indicates that concurrent drought conditions play a dominant role in the occurrence of fires in the NEC (β 1 > β 2 ) and SEC regions. The NEC region experiences high summer temperatures (averaging approximately 30 • C) and cold, dry winters (averaging approximately −30 • C). The majority of the annual rainfall occurs during the summer (around 500 mm) [48]. Temperatures gradually increase beginning in March, leading to melting snow and ice and a moistening of the land. Starting in October each year, temperatures significantly drop [49]. The main forest type is deciduous coniferous forest, which is difficult to dry and therefore has a relatively long fire season. However, the wood contains a large amount of resin, making it prone to large fires. A particularly interesting aspect is that the wetter climate conditions last year also led to an increase in BA in the NEC region. This may be because the previous wetter climate helped plant growth and ensured the continuity of fuels. The SEC region experiences consistently high temperatures throughout the year [50], leading to a prolonged fire season that spans from January to November. However, the SEC region has an annual average rainfall of up to 1780 mm [51], which results in not much burned area (BA). There is a moist humus layer under the forest, and the continuous high temperature in the fire season will make it dry, causing frequent fire. However, rainfall and a humid climate generally prevent large-scale fires in this area [45].
Antecedent drought is the key predictor in most regions, which may be because the fuel load and factors related to ignition are controlled by previous climatic conditions. For instance, drying of fine fuels during the early months of the fire season may encourage wildfire activity [38,52]. Fire occurrence is determined by a delicate balance between vegetation productivity (fuel production) and the frequency of dry conditions (fuel moisture) [9]. Most eastern regions (NC, CC, and SC) are more sensitive to antecedent drought than western regions (SWC and NWC). This is maybe due to the large terrain, low population density, limited vegetation coverage, and small BA in the western region (SWC and NWC). Its fuel load and continuity are considerably weaker when confronted with a variety of ecological issues including desertification, local forest destruction, grassland degradation, etc. [53][54][55][56][57].
In NWC and CC regions, the results of the model show that BA decreases slowly with time. This may be the result of a combination of vigilance and a set of preventive measures coordinated by local authorities, firefighters, and forest managers. The NEC region shows a trend of slow growth over time, which may be the impact of climate warming.

Conclusions
In general, SPEI and sc-PDSI demonstrate good applicability in the seven regions of China, although there are variations in the time window across different regions. The linear regression models we developed demonstrate relatively high ability in both the within-sample and out-sample. These relatively simple regression models link drought index with fire activity, which can be used to estimate the response of fires in the Chinese mainland to different climate-change scenarios, assuming that the climate-vegetationhuman-fire interaction and its feedback do not change significantly. However, the complex relationship between fire, climate change, human activities, and vegetation distribution may limit the applicability of these findings under very different conditions from the current situation. Future research lines should take into account a wide range of vegetation, human, and meteorological factors. The models showed good predictive power, indicating their potential usefulness in predicting future wildfire occurrences. Overall, these findings provide valuable insights into the spatial distribution and potential predictors of wildfires in mainland China, which can inform wildfire management and prevention strategies.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.