Prediction model of moisture content of dead fine fuel in forest plantations on Maoer Mountain, Northeast China

Preventing and suppressing forest fires is one of the main tasks of forestry agencies to reduce resource loss and requires a thorough understanding of the importance of factors affecting their occurrence. This study was carried out in forest plantations on Maoer Mountain in order to develop models for predicting the moisture content of dead fine fuel using meteorological and soil variables. Models by Nelson (Can J For Res 14:597–600, 1984) and Van Wagner and Pickett (Can For Service 33, 1985) describing the equilibrium moisture content as a function of relative humidity and temperature were evaluated. A random forest and generalized additive models were built to select the most important meteorological variables affecting fuel moisture content. Nelson's (Can J For Res 14:597–600, 1984) model was accurate for Pinus koraiensis, Pinus sylvestris, Larix gmelinii and mixed Larix gmelinii—Ulmus propinqua fuels. The random forest model showed that temperature and relative humidity were the most important factors affecting fuel moisture content. The generalized additive regression model showed that temperature, relative humidity and rain were the main drivers affecting fuel moisture content. In addition to the combined effects of temperature, rainfall and relative humidity, solar radiation or wind speed were also significant on some sites. In P. koraiensis and P. sylvestris plantations, where soil parameters were measured, rain, soil moisture and temperature were the main factors of fuel moisture content. The accuracies of the random forest model and generalized additive model were similar, however, the random forest model was more accurate but underestimated the effect of rain on fuel moisture.


Introduction
Forest plantations provide timber and wood products, thereby contributing to the economic development of many countries (Carnus et al. 2006). Recent statistics of forest cover show that China has 24% of the world's forest plantations (Shen et al. 2019). In temperate and boreal regions, forests are threatened by fires whose frequency and intensity tend to increase due to human activities and climate change (Beach et al. 2010). However, fire is also one of the main factors contributing to forest renewal and succession, and plays an essential role in maintaining biodiversity and landscape heterogeneity (Chang et al. 2013). Despite their natural functions, fires of high intensity constitute a danger to forest ecosystems, human life, and infrastructure, and are a source of greenhouse gas emissions. Therefore, preventing and suppressing forest fires is one of the main tasks of forest agencies to reduce losses. This requires a good understanding of factors affecting fire occurrence.
According to the Global Fire Monitoring Centre (https :// gfmc.onlin e/inven tory/stati stic.html) of Germany's Freiburg University, fires occur on several hundred million hectares of forest and other vegetation around the world each year. In China, forest fires are reported in different provinces of the country, however up to now, forest fire statistics are not reported to the Global Fire Monitoring Centre. In Heilongjiang Province in northeast of China, the mean annual area of burned forest varies between 500 and 2500 km 2 (Li et al. 2015). This average represents 0.25− 1.2% of the forest cover of Heilongjiang Province, estimated at 205,328 km 2 according to the China Forest Administration website (http:// fores t.ckces t.cn/, 12 December 2019).
To fight forest fires, the China Forestry Administration carries out a national program to define forest fire prevention and firefighting measures, including management of dead forest fuels. Research is being carried out in the Maoer Mountain forest ecosystem on live and dead fuel loads, and on the prediction model of fuel moisture content for fire risk management. In this regard, the present study addresses the second aspect of this research project to establish models for predicting moisture content of dead fuels in plantations on Maoer Mountain forest ecosystem. Dead fuels are defined as "recently and partially decomposed tree leaves, needles, twigs, grass, bark and branches, with no living tissue" (Simard 1968). They are divided into different time classes (1, 10, 100, and 1000 h) as a function of their response time lag according to meteorological conditions, in particular relative humidity and temperature (Bradshaw et al. 1984).
Fuel load and moisture content, topography and wind speed are the main factors that control forest fires (Finney et al. 2010;Syphard et al. 2019), and fuel moisture is reported to be the most critical factor affecting fuel ignition (Possell and Bell 2013). Thus, identifying factors that affect fuel moisture content is crucial for fuel management, and fire prevention and suppression (Liu et al. 2013;Slijepcevic et al. 2015). As shown by Byram and Jemison (1943), fuel moisture content is a result of cumulative effects of past and current weather events. Moreover, Flannigan et al. (2016) showed that the moisture content of fuel is governed by a composite interaction between fuel characteristics, time, and environmental conditions such as rain and temperature. In general, factors that govern the dead fine fuel moisture content are relative humidity, temperature, precipitation, radiation, wind, as well as dew, soil humidity, soil temperature, slope, topography, aspect, and canopy density (Simard 1968;Nelson 2001;Chuvieco et al. 2004). For example, in rainless conditions, dead fine fuels have a moisture content less than 30%, which corresponds to the saturation point of wood fibres when the relative humidity is 100%. On the other hand, when fuels are exposed to rains, they store water in cell cavities and on their surface. In this case, their moisture content can reach 300% (Slijepcevic et al. 2013).
Models of fuel moisture levels based on weather conditions have been developed using empirical, semi-physical or physical methods (Simard 1968;Catchpole et al. 2001;Nelson 2001;Wittich 2005;Matthews 2006;Matthews et al. 2010). Empirical models provide a useful approach for estimating the danger of wildfires, however they are limited in their ability to represent complex processes such as precipitation or variation in solar radiation, or to be extrapolated to other types of fuel (Matthews et al. 2010). Approximately twenty years ago, physical models were developed (Weise et al. 2005;Matthews et al. 2010) based on previously built models (Fosberg 1975;Viney and Hatton 1990) which included some relevant processes but not all. The model by Matthews (2006) is process-based to representing the flux of energy and water in a litter bed composed of three elements: litter, air, and water. The flux of heat and water of each of these three elements is calculated using equations based on six variables: litter moisture content (kg of water per kg of dry litter), temperature of the free water on the litter surface, amount of water on litter surface (kg of water per m of litter), litter temperature, air temperature, and relative humidity. In the physical model, the flux of heat, water and radiation between the litter and the atmosphere or soil is predicted by replacing variables of the flux of heat and water in the process-model by specified boundary conditions: relative humidity, temperature, precipitation, solar radiation, thermal radiation, wind speed, soil moisture, and soil temperature (Slijepcevic et al. 2013).
In this study, it is hypothesized that relative humidity, precipitation and temperature govern the moisture content of dead fine fuels in forest plantations of the Maoer Mountain forest ecosystem. The objective was to evaluate existing equilibrium moisture content models and to predict fuel moisture content using meteorological variables. The specific objectives were: (1) to choose the best equilibrium moisture content function that could be used during prescribed burning in plantations; (2) to develop a prediction model of fuel moisture content with meteorological variables; and, (3) to determine the relationship between fuel moisture content with meteorological and edaphic variables.

Study area
A forest plantation in the Maoer Mountain forest ecosystem, 45° 43 ′N and 126° 37′ E in Harbin, Heilongjiang province at mean altitude of 255 m was selected (Fig. 1). The climate is cold and temperate with rainy summers. According to the Köppen and Geiger climate classification system, the area is a Dwb climate (humid continental). Yearly precipitation and average temperature are 649 mm and 2.7 °C, respectively. July is the hottest month with average temperatures of 21.8 °C, while January is the coldest with average temperatures of − 19.9 °C. The annual thermal amplitude is 41.7 °C, and the rain amplitude is 171 mm. The parent material is granite bedrock and the soil is dark brown forest soil (Wang 2006). The primary forest was dominated by Korean pine mixed with deciduous species such as Betula spp., Larix spp., Populus spp., and Quercus spp. From the beginning of the twentieth century, the original forest has been largely destroyed by large-scale industrial logging by Russian and Japanese companies as well as by the Chinese government and replaced by secondary forest and plantations (Chen et al. 1982). At present, there are three main types of secondary forest distributed over various site conditions: a Mongolian oak forest on steep upper arid and infertile slopes, a mixed deciduous forest on well-drained fertile, gentle mid-slopes and a deciduous forest on gentle, moist, and fertile slopes. There are two dominant plantations of Pinus spp. and Larix spp. (Wang 2006). This study was carried out in forest plantations because models on fuel moisture content in these ecosystems have not as yet been established in the area, and both natural forests and plantations are subject to wildfires. Plantations of Pinus koraiensis Siebold and Zucc., Pinus sylvestris L., Larix gmelinii Rupr., and mixed Larix gmelinii-Ulmus propinqua Koidz were selected because they produce large fuel loads that can easily ignite and spread fire. The fuel load and thickness estimates in Maoer Mountain forest ecosystem amount to 9.2 ± 1.9 t ha −1 and 2.9 ± 1.4 cm for P. koraiensis, 12.9 ± 3.9 t ha −1 and 4.3 ± 0.7 cm for L. gmelinii, and 8.5 ± 0.5 t ha −1 and 2.3 ± 1.4 cm for P. sylvestris (results not published). Forest plantations in Heilongjiang Province cover 40,420 km 2 , of which 87.4% are young and mid-aged. Coniferous plantations of L. gmelinii, P. koraiensis, P. sylvestris, Abies sibirica Ledeb. and A. nephrolepis (Trautv. ex Maxim.) Maxim. cover 24,560 km 2 , representing more than 60% of Heilongjiang's forest plantation area (Li Wenwu, n.d.). The fuel moisture content was measured in four different stands located between 127° 39′ 39.5″ E and 45° 24′ 13.75″ N. More details about sampling plots are presented in Table 1.

Fuel moisture content meter
A fuel moisture content (FMC) meter consists of an automatic balance and a mini weather station operating with electrical energy produced by solar panel. It was designed by the Forest Fire Prevention Laboratory of the Northeast Forestry University and automatically measures fuel mass, air temperature, relative humidity, wind speed, solar radiation, and rainfall reaching the fuels on the ground as well as soil moisture and temperature. To check the accuracy of each FMC meter balance, a sample of dead fuel was weighed and its mass was compared to that obtained using the FMC meter in the field (0.01 g of error). Components of the FMC meter are shown in Fig. 2.

Field sampling
Data collection was carried out from September 23 to November 15, 2018. For each plantation, altitude, canopy density (shade), slope, aspect, height and diameter at breast height (DBH) were recorded as they affect solar radiation, temperature and relative humidity on the soil surface under the canopy. In each stand, a 50 × 50 m sampling plot  was established and five sampling points were randomly installed where fuels were evenly distributed. Samples of dead leaves, needles and small twigs < 3 mm (Burgan and Rothermel 1984) were collected at random points (Schunk et al. 2017) and placed in the FMC 30 cm × 30 cm × 6 cm basket at each sampling point, and oven dried for 24 h at 105 °C. Samples were then weighed and placed back in the plot on the ground (duff layer) in a basket attached to the FMC meter scale (Fig. 2). To avoid disturbance of samples, a protective 1-cm mesh was put around the FMC meter. The initially sampled and dried fuels were kept on the FMC meter throughout the study period (i.e., there was no repeated manual sampling).
To analyse the variations in fuel moisture content, halfhourly meteorological and soil data provided by the FMC meter in different sampling plots were recorded. The FMC meter automatically recorded fuel weight, relative humidity, wind speed, solar radiation and temperature. Data were shared through hotspot from the FMC meter to an Android cell phone on which the fmcmeter software was installed (see supplementary Appendix 1). As a rain gauge and soil temperature and moisture sensors were not operational throughout the data collection, data were provided by the Global Ecology Laboratory of the Northeast Forestry University. This laboratory measures soil temperatures and humidity at a 5-cm depth by CS107L sensor (Campbell Company) connected to an automatic acquisition system (CR23X, Campbell).

Data analysis
The gravimetric method was used as it is regarded as most reliable, inexpensive, and accurate. It involves taking fuel samples of a predetermined size (fine, mid-size, large), and oven drying these samples (Pippen 2008). The fuel moisture content (FMC) determines the amount of water in the sample by calculating the quotient of difference between the wet mass and the dry mass (Slijepcevic et al. 2015) using the following function:FMC = where W w is the fresh weight and W d the drying weight. Various methods have been used for estimating fuel moisture content (Schroeder and Buck 1970). Among these, two kinds of models were developed in this study: the equilibrium moisture content model and the meteorological regression model. The random forest model was developed using the train function in randomForest and e1071 packages, and a generalized additive model using the gam function in mgcv package in RStudio version 1.1.453.0 (RStudio Team 2018) with generalized cross-validation (GCV) or restricted maximum likelihood (REML) processes. Graphs of fuel moisture content variation were drawn with OriginPro2018 (Origin Lab, Northampton, MA, USA).
Test and validation of selected equilibrium moisture content models Nelson (1984) model was used because it is widely accepted and incorporated into many physical models (Catchpole et al. 2001;Wittich 2005;Matthews 2006). Van Wagner and Pickett (1985)'s model was tested because it had been built using coniferous fuels and its use is still operational. Simard (1968) and Nelson and Hiers (2008) models have not been tested because they are weak predictors of fine fuel moisture content (Pippen 2008). In addition, models of Catchpole et al. (2001), Matthews (2006), and Wittich (2005) have not been tested because they are no longer operational. Other models not mentioned here were not tested because their operational use is considered unclear (Pippen 2008). Thus, among all of these functions only the ones of Nelson (1984) and Van Wagner and Pickett's (1985) were adapted using temperature and relative humidity. The following functions mathematically represent these models: (a) Nelson (1984) equilibrium moisture content response to environmental factors is defined as: where R is the universal gas constant (8.314 J K −1 mol −1 ); T is the ambient air temperature (°K); H is the air relative humidity (%); m is the H 2 O relative molecular mass (18 g mol −1 ); α, β are the parameters to be estimated. (b) Van Wagner and Pickett (1985) model is: where H is the relative humidity and T the temperature (°C).
To evaluate the equilibrium moisture content models, the root-mean square error (RSME) and mean absolute error (MAE) were used to choose the best model.

Meteorological regression method
Studies on fuel moisture content prediction have been conducted in natural forests in northeast China (Zhang et al. 2014;Sun et al. 2018). However, these studies have used an empirical method whose model's extrapolation are often limited. The present study used a physical method for developing predicting models that could be extrapolated to other sites or other types of fuels, taking into account the relevant processes of material, energy and heat flux between litter and atmosphere or soil. The prediction of fuel moisture content was developed using a random forest model and a generalized additive model.
Random forest model was selected as it is highly accurate, robust to outliers and noise, faster than bagging or boosting, and most interpretable (Breiman 2001). Variables of the best model were defined by combining the variable importance score to the number of variables (see "mtry") in the best model. For each model fitted, the optimal number of significant predictors was identified by "mtry" parameter with the 10-folds cross-validation method. The best model was determined by the optimal number of significant predictors, the lowest root-mean square error and mean absolute error, and the highest value of the adjusted R 2 .
Generalized additive model (GAM) was fitted to predict FMC models as it describes the relationship between dependent and independent variables without prejudging the issue by selecting any particular parametric method (De Leeuw 2009). GAMs cover an extensive range of models, such as multiple linear, generalised linear, and mixed effect models (Wood 2018). They extend linear and generalized linear models to include smooth functions of independent variables with smoothness determined by either a parameter that directly controls the curve smoothness or the estimated predictive accuracy. They are determined by the following function: where μ i is the mean of the dependent variable, s 0 the parametric coefficient, s i the smooth functions defining the additive component, and X i the independent variables; si(X i ) is the smooth terms, s i.j (X i ·X j ) the tensor product terms, and Ɛ the error. The smooth terms are text series giving the names of the covariates whose term is a function. The tensor product smoothers result from smoothing interaction of any or all of covariates. They are a natural way to represent smooth terms interaction in regression models as they are invariant compared to units in which covariates are measured, thus avoiding the need for arbitrary decisions regarding relative scaling of variables (Wood et al. 2013). They are suitable for models of variables with different units.
Generalized additive models were defined using a Gaussian family and an identity link function. For each FMC model, we first ran the model including all covariates, then consecutively removed covariates with the highest concurvity until all covariate concurvity estimates were < 0.4 (Lydersen et al. 2019). This restricted the effect of correlation between covariates in the model. Model prediction accuracy was determined by the generalized cross-validation (GCV) and restricted maximum likelihood (REML) methods. With the GCV, each sample was randomly divided into k (unfixed parameter, 10 in this case) equal-sized subsamples. One of ten subsamples was retained as validation data to test the model, and the nine remaining subsamples were used as training data. The cross-validation process was repeated ten times and resulted in a single estimation. The benefit of the GCV is that all observations are used for both training and validation (Golub et al. 1979). Selection of best models was made with the lowest GCV/REML and AIC, and highest adjusted R 2 .

Results
Validation test of selected equilibrium moisture content models Nelson (1984) and Van Wagner and Pickett (1985) models were tested with field data (temperature and relative humidity) to ensure their validity during controlled burning. After model fitting, only the Nelson's model showed good accuracy, with RMSE and MAE less than 1%, and was validated as suitable for prescribed burning in Maoer Mountain forest ecosystems. The Van Wagner and Pickett (1985) model was less accurate, with RMSE and MAE ≥ 5% (Table 2).

Ranking importance of variables using the random forest method
Plots of variable importance were built with the random forest method (Fig. 3). In the P. koraiensis plantation, two graphs, (a) one showing the variable importance of the model built with both meteorological and edaphic variables, and (b) the other showing the variable importance in the model with meteorological variables. Observations revealed that the soil moisture content significantly affected the FMC (Fig. 3a). In all plantations (without considering soil variables), relative humidity and temperature were the most important variables that affected the FMC.

Fuel moisture content variation
Temporal changes in FMC were due to different factors. After placing the fuels back into the FMC meter, a sudden increase had been observed with P. koraiensis fuel (Fig. 4) and mixed P. sylvestris-P. koraiensis fuel (Fig. 5a) as it had rained in the evening of that day. In the L. gmelinii-U. propinqua planation (Fig. 5b) and L. gmelinii plantation (Fig. 6), samples were put back one day after the rainfall beyond October 20 would be related to snowmelt, in addition to rainfall, temperature, relative humidity, solar radiation, wind speed, soil moisture, and temperature measured in this study (Figs. 4,5,6). Dahurian larch fuels had a higher water holding capacity than pine fuels. In the same way that these fuels absorbed rain water during rainy periods, they gradually lose it compared to pine fuels which were sensitive to small changes in temperature or sunshine. The solar radiation contributed considerably to the FMC in pine plantations, and wind speed significantly impacted the FMC  maximum values of parameters are presented in Table 3. Meteorological records at the time of FMC sampling varied from − 10.6 ºC to 25.5 ºC for temperature, 23.1-96.8% for relative humidity, zero to 2.8 mm for rainfall, zero to 2.95 m s −1 for wind speed, and zero to 1062 W m −2 for solar radiation. Soil variables varied from 1.7 ºC to 12.7 ºC for temperature and 21.5-28.4% for l humidity. The max value of FMC (38.7%) was observed in the L. gmelinii plantation and was related to higher values of relative humidity, low wind speeds (Table 3), canopy density and litter thickness (Table 1). Comparison between Fig. 4a, b which present the graphs of observed and predicted values of fuel moisture content in the P. koraiensis plantation shows that integration of soil humidity and temperature increases the models (RF and GAM) accuracy.

Model validation
To evaluate the accuracy of the random forest model, the variable importance score and cross-validation method were used. Random forest models showed that temperature and relative humidity were the most important factors affecting fuel moisture content. For all selected models, the RMSE and MAE values were below 5%. The model fitted with meteorological and soil variables (in the P. koraiensis plantation) was reasonably accurate (lower error values and a high value of R 2 -adj = 0.87) compared to the model fitted with meteorological variables only for which adjusted R 2 equalled 48%. In all plantations, models fitted with meteorological variables had RMSE and MAE vales lower than 5% and adjusted R 2 between 39 and 63% (Table 4).
GAM prediction showed that temperature, relative humidity and rainfall were the most important drivers of fuel moisture content. The estimated degrees of freedom and model accuracy parameters in all plantations are presented in Table 5 (for generalized cross-validation method) and in Table 6 (for restricted maximized likelihood method). Beyond the combined effect of temperature, rainfall and relative humidity, solar radiation or wind speed effect was also significant in some sites. Models fitted with soil variables (without soil variables) had explanatory power that varied from 32.8% (31.7%) to 50.0% (48.6%). However, in the P. koraiensis plantation where soil variables were measured, GAM explained 82.8% of FMC variation in the model fitted with both meteorological and soil variables; however, it explained only 39.1% (36.3%) of the FMC variation in the model fitted with only meteorological variables. Fuel moisture content varied mainly with the combined effect of soil moisture, rainfall, soil temperature, and secondly with the combined effect of relative humidity and temperature as well as solar radiation. The best model in the P. sylvestris and P. koraiensis plantation included both the combined effect of relative humidity, temperature and rainfall, and the combined effect of temperature and radiation and explained 38.1% (36.8%) of the deviance. The best model in the L. gmelinii and U. propinqua plantation included both the combined effect of relative humidity, temperature and rainfall and the combined effect of temperature and wind, and explained 50.1% (48.6%) of the deviance. The best model in the L. gmelinii plantation had particularly low explanatory power, accounting for 32.8% (31.7%) of the deviance. Comparing both methods using AIC, GCV accuracy was better than REML. Significant differences between models fitted with both meteorological and soil variables, and models fitted with only meteorological variables explained sufficiently the contribution of soil moisture and temperature on fuel moisture content.

Discussion
An overall outlook was made on the main weather factors affecting dead fine fuel moisture levels in Maoer Mountain forest plantations. Selected equilibrium moisture content (EMC) models were tested, importance score of variables measured to identify the best attributes influencing fuel moisture content, and prediction models of FMC developed and adopted. Nelson's (1984) model predicted the equilibrium moisture content better than the Van Wagner and Pickett (1985)   model. However, the EMC model is useful during predominantly dry conditions. In all sampling sites, the RMSE and MAE values of Nelson's (1984) model were lower than 1%, whereas for Van Wagner and Picket's (1985) model, these values were above 5%, and even reached 21%. Our results are consistent with previous findings of Sun et al. (2015). The relatively good accuracy of Nelson's (1984) model may be explained by the fact that it is a semi-physical model, whereas Wagner and Picket's (1985) model is an empirical one (Yu et al. 2018). The Nelson (1984) model accuracy was consistent with what was suggested by Pippen (2008), stating that a better model of fuel moisture content must have precision in the range of 1-2% accuracy to ensure its effectiveness in fire risk management. This concurs with Hu et al. (2016) who found the EMC method suitable for predicting FMC on small forest scales. The generalized additive model showed that the combined effects of temperature, relative humidity and precipitation on FMC was significant. This confirms our hypothesis that rainfall, relative humidity and temperature are the main drivers affecting fuel moisture contents in P. koraiensis, P. sylvestris, L. gmelinii, and mixed L. gmelinii-U. propinqua plantations on Maoer Mountain. Solar radiation or wind speed were additional predictors in the model. Their effect would be influenced by aspect, canopy density, litter thickness, and altitude. The random forest and generalized additive models responded differently to rainfall. The RF model appeared to provide more precision than GAM however, it underestimated the rainfall effect on fuel moisture content. In addition, because our models were developed with data taken on the surface where fuels are located provides greater accuracy compared to models developed with data from weather stations far from sampling sites. These results are consistent with Sen and Wan-long (2014) showing that precipitation strongly affects FMC and share a number of similarities with Zhang et al. (2014) and Chae et al. (2017) that showed that fuel moisture content is firmly related to rainfall and relative humidity, and with results of Bilgili et al. (2019) that showed that FMC correlates well with relative humidity and temperature for all types of fuels. The effect of wind speed on FMC was significant only in the L. gmelinii-U. propinqua plantation due to its proximity to an open area. However, its effect on FMC in other plantations was not significant, in agreement with Larsen (2015) who showed that the influence of wind speed on the moisture content of fuel is not apparent, mainly because the canopy of trees reduces wind speed or prevents it from reaching fuels on the soil surface. Although soil moisture and temperature were measured in one plantation, results show that their effect on fuel moisture was significant and are consistent with Krivtsov et al. (2008) and Matthews (2014) . In addition to the variables considered in this study, slope, aspect, and altitude should be incorporated in future models for a better understanding of environmental conditions on fuel moisture contents in Heilongjiang forest ecosystems.
In contrast to Byram and Jemison's (1943) method that developed a fuel moisture content model with current and past weather data, this study used current data because past data are documented on a daily average, while this study developed models with half-hour data. Although the models established in this study appear to be weak using current data, they are preliminary tools for monitoring moisture content of fuel in a shorter time interval. The use of Byram and Jemison's (1943) method also has drawbacks because its results are based on historical data while its application should relate to future forecast analysis. However, a comparative study on daily averages of current and past data will be carried out to evaluate the trend, especially since climate change is hypothesized to impact meteorological variables. Our results may have limitations for operational use. The first is that data in this study were tested in one season, thus giving a partial view of FMC variations. The second is that fire managers will use data of weather stations that are often in open areas far from local sites compared to measurements taken with an FMC meter in the field under tree canopies. In the next step in data collection, outdoor meteorological data will be included to make the study more applicable. The third limitation is that both random forest and generalized additive models do not have a mathematical formula that can be used without statistical software. But a general linear model could be fitted with significant variables, although its precision is lower than that of both.

Conclusion
In this study, the equilibrium moisture content function and meteorological regression method were used to predict the moisture content of dead fine surface fuels in plantations in the Maoer Mountain forest ecosystem. To achieve this, plantations of pine and Dahurian larch were selected as they produce large fuel loads that constitute a significant factor for forest fire management. Variable importance score was determined, and models were established. Nelson's (1984) function was validated for prescribed burning as it provided the best accuracy for predicting fuel moisture levels in sampled plantations. As expected from our hypothesis, the effect of temperature, relative humidity and rainfall on FMC was significant. Prediction model accuracy was substantially upgraded by including soil variables. The accuracy of the random forest model and generalized additive model was similar, however the accuracy of the random forest model was slightly higher but underestimated the effect of rainfall on FMC. The data used in this work were tested in one season, thus giving a partial view of fuel moisture content variations. Further investigations on dead and live fuels are needed to estimate fuel moisture contents of typical forest plantations in the region.