Next Article in Journal
Association between the Rail Breakage Frequency in Beijing–Tianjin–Hebei High-Speed Railway and the Eurasian Atmospheric Circulation Anomaly
Next Article in Special Issue
Visualization of Prediction Methods for Wildfire Modeling Using CiteSpace: A Bibliometric Analysis
Previous Article in Journal
URANS Simulations of Vehicle Exhaust Plumes with Insight on Remote Emission Sensing
Previous Article in Special Issue
A Study of Two High Intensity Fires across Corsican Shrubland
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Applying Bayesian Models to Reduce Computational Requirements of Wildfire Sensitivity Analyses

1
Agriculture and Food, Commonwealth Scientific and Industrial Research Organisation (CSIRO), St Lucia, Brisbane, QLD 4067, Australia
2
Department of Infrastructure Engineering, Faculty of Engineering and Information Technology, University of Melbourne, Melbourne, VIC 3053, Australia
3
School of Public Health, Faculty of Medicine and Health, University of Sydney, Sydney, NSW 2006, Australia
4
Data61, Commonwealth Scientific and Industrial Research Organisation (CSIRO), Clayton, VIC 3168, Australia
5
Cloud Computing and Distributed Systems (CLOUDS) Lab, School of Computing and Information Systems, The University of Melbourne, Melbourne, VIC 3053, Australia
*
Author to whom correspondence should be addressed.
Current address: Building 80, The University of Queensland, 306/306 Carmody Rd., St Lucia, Brisbane, QLD 4072, Australia.
Atmosphere 2023, 14(3), 559; https://doi.org/10.3390/atmos14030559
Submission received: 2 December 2022 / Revised: 7 February 2023 / Accepted: 13 February 2023 / Published: 15 March 2023

Abstract

:
Scenario analysis and improved decision-making for wildfires often require a large number of simulations to be run on state-of-the-art modeling systems, which can be both computationally expensive and time-consuming. In this paper, we propose using a Bayesian model for estimating the impacts of wildfires using observations and prior expert information. This approach allows us to benefit from rich datasets of observations and expert knowledge on fire impacts to investigate the influence of different priors to determine the best model. Additionally, we use the values predicted by the model to assess the sensitivity of each input factor, which can help identify conditions contributing to dangerous wildfires and enable fire scenario analysis in a timely manner. Our results demonstrate that using a Bayesian model can significantly reduce the resources and time required by current wildfire modeling systems by up to a factor of two while still providing a close approximation to true results.

1. Introduction

Each year natural disasters such as wildfires cause harm to people and significant destruction of physical infrastructure. To better understand these phenomena and reduce their impact, various natural hazard models have been developed over the years [1,2]. The effectiveness of dynamic systems based on these models often depends on how quickly they can predict the unfolding of events. Currently, several simulation models, such as Firemap [3], SiroFire [4], Prolif [5], Farsite [6], Pyrocart [7], Firemaster [8], FireStation [9], Prometheus [10], Spark [11], and Phoenix [12], are used to predict the spread of wildfires across a landscape based on pre-existing fire spread models. Deriving accurate risk metrics often requires a large number of fire simulations, and the inputs to these wildfire spread models have associated uncertainties that can influence the resulting rate of speed and, consequently, the area burned by the fire. Quantifying the sensitivity of these parameters to the resulting output is useful for worst-case scenario analysis in operational risk management [13]. However, such scenario analysis for effective wildfire management requires additional simulations to be run under various combinations of input values, which can be both computationally expensive and time-consuming [14].
Several methods and tools are available for conducting sensitivity analysis in the environmental and wildfire context. These include variance-based approaches such as Sobol’s method [15], Cukier’s method [16], and Saltelli’s method [17]; density-based methods such as Krzykacz’s method [18], Plischke’s method [19], and Pianosi’s method [20]; neural network methods [21,22]; Taylor series expansion [23]; and polynomial chaos expansion [24]. Additionally, several tools such as the Monte Carlo Analysis Toolbox (MCAT) [25], Eikos [26], MATLODE [27], SAFE (Sensitivity Analysis for Everybody) [28], SALib [29], and OpenTURNS (Open source Treatment of Uncertainty, Risk ’N Statistics) [30] have been developed for sensitivity analysis. However, these tools can only be used once the sets of input and output values are available. While these methods and tools are useful for conducting sensitivity analyses of fire simulations, they do not consider the computational requirements of running the fire simulations under a large number of input combinations, as required by these methods and tools. As a result, using a conventional wildfire management system to run any of these methods or tools for scenario analysis in emergency planning and management may be prohibitively time-consuming, hindering the ability to make better-informed decisions to minimize the extent of damage caused by the disaster. In this work, we investigate whether Bayesian models can reduce the computational requirements of such analyses while maintaining a close approximation to true results and assessing wildfire impacts.
The Bayesian approach is a widely used method in the literature for constructing models to explain various phenomena. It involves fitting a probabilistic model to a given set of data to summarize and predict new observations. In the wildfire domain, the Bayesian approach has been used in several applications; prediction of the likelihood of large fires [31,32], projection of wildfire activities [32,33], estimation of fire suppression costs and resource allocation [34,35], estimation of the size of extreme fires [36,37,38], wildfire risk assessment [39,40,41,42], and prediction/modeling wildfire behavior [43,44,45,46,47]. One advantage of using Bayesian models in wildfire applications is that it allows for the incorporation of prior knowledge to any observation data sizes to make approximate estimations. Despite requiring a large number of simulations, wildfire management practices, such as scenario analysis, can benefit from the Bayesian model to reduce computational costs while still providing realistic predictions.
The main objective of this paper is to investigate the application of a Bayesian model for combining prior knowledge and data to make predictions for wildfire management practices, such as scenario analysis. These predictions will help to significantly reduce the computational requirements of such analyses as only a fraction of the entire input combinations would now be required to be run using the fire simulation tool. The proposed Bayesian model is used to estimate the impacts of wildfires at a location based on three major meteorological inputs: temperature, relative humidity, and wind speed collected from weather stations. The choice of these parameters is based on the experimental setup of previous works [48,49], and the goal of this investigation is to determine the potential time and resource savings that can be achieved through the use of a Bayesian model in wildfire modeling systems.
We build two Bayesian models for comparison, one with and one without a latent effect. The use of a latent effect is investigated to determine whether it improves the accuracy of the model and under what circumstances it can be ignored for the simplicity of the model. If the addition of the latent effect does not improve the model’s accuracy, it can be omitted for the sake of simplicity. Additionally, we evaluate the performance of the models under different priors for the hyperparameters associated with the input parameters to the fire simulations. We then apply the best-performing Bayesian model to estimate the impact of a wildfire (fire size) based on the available data and prior information. The predicted values, along with the data, are used to estimate the sensitivity of the fire area to the input parameters, which can enable scenario analysis in wildfire modeling systems.
The remainder of the paper is organized as follows. In Section 2, we provide a detailed overview of the workflow for this study. In Section 3, we explain the experimental setup for the study, while in Section 4, we present the results and discuss our findings. Finally, in Section 5, we offer conclusions and suggestions for future work.

2. Workflow

Figure 1 illustrates the overall workflow of the study that investigates the reduction of computational requirements for sensitivity analysis in wildfire management through the use of a Bayesian model. The following section outlines the Bayesian model fitting process and then describes the experimental setup for its application.

2.1. Bayesian Model Fitting

For Bayesian model fitting, we first investigated the influence of latent effects in the model before identifying the best prior for the precision parameter in the model.

2.1.1. Influence of Latent Process

In our Bayesian model fitting, the joint posterior distribution of parameter θ is conditional on the observed data y , and hence we write the posterior distribution as follows.
p ( θ | y ) = p ( y | θ ) p ( θ ) p ( y )
where p ( y | θ ) is the likelihood of the data given all the parameters θ , p ( θ ) is the prior distribution of the parameters, and p ( y ) is the marginal likelihood. p ( y ) is a normalizing constant and is mathematically equal to:
p ( y | θ ) p ( θ ) d θ
The dimension of the posterior distribution p ( θ | y ) usually depends on the parameter dimension, i.e., on θ , and can be obtained using joint posteriors as follows.
p ( θ | y ) p ( y | θ ) p ( θ )
In other words, the posterior distribution of the parameters can be estimated by scaling the product of the likelihood and the prior.
We considered the Bayesian model shown in Equation (4).
y N μ , I n σ ϵ 2
where y is the fire area in hectares that follows a normal distribution with latent mean μ and variance σ ϵ 2 . Here, I n is an identity matrix of order n. For each fire area, we use the Spark input data, i.e., the meteorological information, to identify the latent process μ , as given in Equation (5).
μ N β 0 1 n + j = 1 3 β j x j , Σ
where β 0 is the intercept and β 1 , β 2 , β 3 are the coefficients associated with the meteorological data x , i.e., we write β j , j = 0 , 1 , 2 , 3 .
We define the term Σ as a correlated fire area variance-covariance matrix. Note that for simplicity, we use a fixed correlation structure S to define the fire areas; thus we write Σ = σ s 2 S . This leads us to use the prior distribution for σ s 2 instead of the whole correlation structure Σ .
To relate the model parameters in Equations (4) and (5) with the posterior distribution in Equation (3), we define θ = μ , σ ϵ 2 , Σ and the log of the joint posterior distribution of the model parameters is given by:
log P ( θ | y ) n 2 log σ ϵ 2 1 2 σ ϵ 2 y μ y μ log | Σ | 1 / 2 1 2 μ Δ Σ 1 μ Δ + log P ( θ )
where Δ = β 0 1 n j = 1 3 β j x j , and P ( θ ) is the joint prior distributions of the model parameters: β j , j = 0 , 1 , 2 , 3 , σ s 2 , σ ϵ 2 . Under the INLA structure, we approximate the joint posterior distribution in Equation (6) and consider the posterior marginals. Hence, we write the marginals for the latent process μ as follows.
P ( μ j | y ) = P ( μ j | σ ϵ 2 , y ) π ( σ ϵ 2 ) d σ ϵ 2
In addition to the Bayesian model that takes latent correlated effects of fire areas into account, we also consider a simpler version of the model that does not consider these effects. For Bayesian model updating, we use the Integrated Nested Laplace Approximation (INLA) [50] framework. INLA is a newer approach for computing Bayesian models that is less computationally expensive than popular Markov chain Monte Carlo (MCMC) methods, such as Gibbs sampling [51], which can provide similar solutions for posterior distributions of the parameters of interest.

2.1.2. Sensitivity to the Priors

To analyze the sensitivity of the models to the priors, we considered six different prior distributions for the parameters in the model. In this paper, we present the sensitivity to the priors only for the precision parameters, as the distribution of the fire area is more sensitive in the tails. The priors we consider include the half-Cauchy, half-t, log-gamma, half-normal, penalized complexity (PC log-gamma), and uniform improper, as listed in Table 1. The values of the parameters of the priors are so chosen to make the priors weakly informative and let the observation data set drive the posteriors. We included the penalized complexity prior in our sensitivity analysis to investigate whether priors based on probability statements about the parameters (PC priors [52]) result in better model performance compared to default priors. We then analyzed the sensitivity by examining the posterior marginals of the three weather input parameters: temperature, relative humidity, and wind speed.

2.1.3. Evaluation Metrics

To compare the models with and without latent effects and different priors, we use several metrics: the marginal likelihood (MLK) [53], Deviance Information Criterion (DIC) [54], Watanabe–Akaike Information Criterion (WAIC) [55], and Conditional Predictive Ordinates (CPO) [56]. The MLK is the probability of the observed data values in the fitted model and can be used to estimate posterior probabilities in the model. DIC and WAIC are model performance criteria that measure the complexity of the model by considering the goodness of fit and penalty term, along with the effective number of parameters. CPO is the posterior probability of observing a value when the model is fitted using all the data except the observation in question. In this paper, we use the CPOs of all the observations, transformed through a log transformation, as shown in Equation (8).
C P O * = i = 1 n l o g ( C P O i )
where n is the total number of observations.

3. Experimental Setup

3.1. Fire Simulation Tool—Spark

We used Spark [11] to simulate and predict the spread of wildfires under different conditions and fuel types. It provides a flexible platform for simulating wildfire behavior in various vegetation types by allowing the integration of different packages and models, such as wind field generation and topographic correction, fire ignition models, fire-line interactions, fireband transport, and fire transmission models. The simulations require input data on fire behavior, land classification, fuel load, topography, and weather to produce output metrics such as total burned area, fire intensity, and the number of urban cells burned. In addition to predicting fire progression, Spark can also predict fireband dynamics and risk metrics for fire severity and impact. It can model firebreaks, spot fires, and the coalescence of different parts of the fire over time and is able to run simulations for multiple fire perimeters simultaneously. The calculations in Spark are parallelized using the OpenCL framework. More information on Spark can be found in [57].

3.2. Weather Inputs

For this study, we selected three weather inputs: temperature, relative humidity, and wind speed, based on the experimental setup of previous research [14,48]. The ranges and distributions for these inputs, as provided in Table 2, follow the same experimental design. These values were chosen to cover operational weather conditions for wildfire modeling in the Australian context and can easily be modified as needed. Wildfires tend to grow more aggressively under conditions characterized by high wind speed, temperature, and low relative humidity. The other static inputs for fire simulations were taken from configurations and records maintained by the Tasmanian government and the TFS [58]. All simulations were run for five hours, and the cumulative burned area during this period was reported as an output.

3.3. Wildfire Management Practice Use Case—Scenario Analysis

We created a wildfire management use case of scenario analysis to identify the worst conditions for aggressive wildfires. The scenario analysis is enabled by the results of the sensitivity analysis, which quantifies the relative influence of each input parameter on fire simulations. The sensitivity analysis uses two sets of data: one with true observations and the other with Bayesian model-predicted values. The study area and evaluation metrics for the application of the Bayesian model in the use case are discussed further below.

3.3.1. Study Area

Tasmania was chosen as the study area for the use case due to the frequent occurrence of wildfires in the region, the availability of high-quality land datasets that can be used in operational wildfire simulation tools, and the systematic grid configuration of fire start locations, which has been well-studied [14,59]. During the 2018–2019 wildfire season, Tasmania experienced 841 wildfires that burned 310,311 hectares of forest [59]. The Tasmania Fire Service has established a grid of 68,048 potential fire start locations at a distance of 1 km, regardless of the type of land, with locations on water bodies shifted to the nearest land location. The model fitting is based on a dataset of fire simulations run at a single location, considering different combinations of the input parameters. Tasmania also has a detailed, high-resolution dataset of simulations, as maintained in one of our previous works [60].

3.3.2. Sensitivity Analysis

We estimated the sensitivity indices of the input parameters in the fire simulation for various sizes of true observations (sample size) used in the Bayesian models. The choice of the priors of precision in the model was based on our initial findings, and the choice was complimented with varied sample sizes to predict the values of fire area for combinations of input parameters. For example, for a sample size of 4000, 4000 random true observations were used to construct a fitted Bayesian model, which then predicted the fire areas for the next 4000 sampled combinations of input parameters. The means of the predicted values were considered for the estimation of the sensitivity indices.
For the sensitivity analysis of fire size to input weather parameters, 8000 samples for different weather input combinations within their ranges were generated using Saltelli’s sampling method [61]. The sensitivity indices were estimated with a variance-based sensitivity analysis (SA) method (Sobol Analysis [15]) using the python framework SALib [29]. The choice of the sampling method aligns with the method for estimating the sensitivity indices.

3.3.3. Evaluation Metrics

To further investigate the potential of the Bayesian model in reducing the computational requirements of scenario analysis in wildfire management, we used the measure of similarity and the reduction in computational requirements. The measure of similarity, calculated using Pearson’s correlation coefficient [62], assessed the closeness between the true observed data and the Bayesian model predicted values for input combinations. The value of the correlation coefficient ranges from 0 to 1, with higher values indicating better model performance. Additionally, we estimated the reduction in computational requirements by calculating the total number of data values predicted by the Bayesian model, which gives a theoretical upper limit of the possible reduction in computational requirements. For example, if the Bayesian model predicted half of the data values used for sensitivity analysis, the computational requirements would be reduced by up to a factor of 2.

4. Results and Discussion

In this section, we present and discuss the Bayesian model for model fitting, the role of the latent process, and the model sensitivity to different priors. We also address uncertainty quantification in the Bayesian model with an increased dataset and the estimation of sensitivity indices using the model-predicted values and the available dataset for scenario analysis.

4.1. Model Fitting

4.1.1. Latent Effects

Table 3 shows a comparison of the two fitted models. All the metrics favored the model with a latent process over the model without one. The model without a latent process had a higher value of marginal likelihood than the model with random noise, indicating that the observed values were more likely to occur in the latter model. Additionally, the lower values of DIC and WAIC for the model with the latent process indicated a better fit with better goodness of fit. The higher calculated values of CPO* also showed that both models had a good fit, with the model with a latent process performing slightly better. These metrics demonstrate that the latent process in the Bayesian model fitting contributed to slightly improved performance. Based on our analysis, the model with the latent process should be preferred. However, due to the insignificant differences between the values of the metrics obtained for both models, both models can be used for predictions with the appropriate priors.

4.1.2. Sensitivity of Bayesian Modeling to Priors

Figure 2 presents the posterior marginals of the temperature, relative humidity, and wind speed model parameters (i.e., β ’s) for various precision priors in the model. The selection of priors can significantly affect the posterior distribution of the fitted model parameters, as demonstrated by the posterior marginals of the input parameters. While most priors yielded similar posterior marginals, the half-Cauchy prior resulted in a significantly displaced posterior margin. In addition, we compared the performance of the models with different priors based on four evaluation metrics and established a preference order for the priors. Figure 3 displays the preference orders for the priors according to these evaluation metrics, with 1 representing the highest performance and 6 the lowest. The preference order for DIC overlapped with the preference order for CPO, and both are depicted in the same order in the figure. Overall, the PC loggamma prior yielded the best model performance, while the uniform and half-Cauchy priors performed relatively poorly. The PC loggamma prior did, however, have the worst value for the MLK metric. Therefore, the Bayesian model with the PC loggamma prior to precision was the most effective in our analysis and should be employed in any further applications of the model in wildfire management practices.

4.2. Wildfire Management Practice Use Case—Scenario Analysis

4.2.1. Similarity between True and Predicted Values

Before estimating the sensitivity indices of the input variables in the fire simulations as an application of the Bayesian model fitting, we calculated the similarity score of the prediction values given by the fitted models (considering the means) with the true data values. Figure 4 illustrates the correlation coefficients calculated for various fitted models. It is evident that the similarity to the true observations increases as the number of true observations considered in the model fitting increases. Although an additional 1000 observations were considered in each model, the improvement in the correlation coefficients was not substantial (only 0.02 from 4000 to 7000 observations).

4.2.2. Scenario Analysis through Sensitivity Analysis to Input Parameters

Figure 5 displays the sensitivity indices estimated for the input parameters in the fitted model, compared to those estimated using actual values obtained from the wildfire simulations. The indices based on the full data sets from the actual model runs indicate that relative humidity, wind speed, and temperature had contributions of 72%, 19%, and 9%, respectively. When the fitted model predicted 4000 values using 4000 true observations, the contributions of the parameters in order were 80%, 13%, and 7%, respectively, which are similar to the true values. As the fitted model considered more data points in the data set for model fitting and predicted fewer data points, the estimated indices became closer to the true values. These levels of influence of the input parameters align with the mean values of the posterior marginals of the three input parameters. Estimating the sensitivity of the fire area to input parameters through sensitivity analysis, in combination with posterior marginals analysis provides important insights into the factors contributing to destructive wildfires and facilitates scenario analysis. Wildfires tend to grow rapidly under high values of temperature and wind speed coupled with low values of relative humidity.

4.2.3. Reduced Computational Requirements

The model fitted with only 4000 data points in our demonstration produced results that were significantly closer to the true values. These findings have important implications for state-of-the-art wildfire management systems, as they suggest that we can trade off the time and computational resources required to run 4000 simulations for which the Bayesian model predicts the values with the time and resources needed to build the Bayesian model. Without considering the Bayesian model fitting, the computational requirements of the fire model simulations are significantly higher. Our application of a Bayesian model for worst-case scenario analysis through sensitivity analysis in wildfire management practices demonstrated promising results, as we were able to obtain close-to-true values with only 4000 fire simulations instead of 8000, saving computational and time resources by up to a factor of 2. Obtaining close-to-true predictions quickly during wildfire emergencies is crucial for effective wildfire management.

5. Conclusions and Future Work

Wildfire modeling systems are critical for understanding the spread of fires and making informed decisions during emergencies. Obtaining close-to-true predictions quickly is crucial in these situations. However, state-of-the-art wildfire management practices, such as scenario analysis, often require a large number of wildfire simulations to be run, which can be computationally expensive and time-consuming. In this study, we demonstrated how probabilistic models built using Bayesian models can be used to improve wildfire management practices. We also examined the impact of a latent effect on the performance of Bayesian models and the sensitivity of the model to different priors for precision. Our application of Bayesian models in estimating the sensitivity of fire size to input parameters showed that this approach can significantly reduce the computational cost and time of wildfire applications by a factor of up to two while still providing close approximations to the true values.
The study presents promising results but has several limitations. The Bayesian model used assumes a linear relationship between inputs and outputs, and only a few parameters and fire area were considered. To further improve the study, the authors plan to build more complex hierarchical Bayesian models with multiple inputs and outputs from fire simulations. Additionally, the study is currently limited to a specific location in Tasmania, and future research will expand this to consider the influence of spatial and fuel characteristics that may vary with fire start locations. Finally, the complexity of building the Bayesian model in terms of algorithms and computation is not considered and will be examined in future research.

Author Contributions

Conceptualization, U.K. and J.A.; Formal analysis, U.K.; Methodology, U.K., J.A. and K.S.B.; Software, J.H.; Supervision, J.A. and J.H.; Validation, U.K., J.A., K.S.B., J.H. and R.B.; Visualization, U.K.; Writing—original draft, U.K., K.S.B. and R.B.; Writing—review and editing, U.K., J.A., K.S.B., J.H. and R.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used for the Bayesian model fitting is a subset of the data set available at https://doi.org/10.25919/XGAM-BS33 [60], accessed on 1 December 2022.

Acknowledgments

The authors would like to thank our academic colleagues and reviewers who helped in improving the quality of the manuscript at various stages of the work. Additionally, the authors are grateful to the guest editor for their generous support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kaizer, J.S.; Heller, A.K.; Oberkampf, W.L. Scientific computer simulation review. Reliab. Eng. Syst. Saf. 2015, 138, 210–218. [Google Scholar] [CrossRef] [Green Version]
  2. KC, U.; Garg, S.; Hilton, J.; Aryal, J.; Forbes-Smith, N. Cloud Computing in natural hazard modeling systems: Current research trends and future directions. Int. J. Disaster Risk Reduct. 2019, 38, 101188. [Google Scholar] [CrossRef]
  3. Vasconcelos, M.; Guertin, D.; Zwolinski, M. FIREMAP: Simulation of Fire Behavior—A GIS Supported System; General Technical Report; United States Department of Agriculture Forest Service: Washington, DC, USA, 1990; pp. 217–221. [Google Scholar]
  4. Coleman, J.R.; Sullivan, A.L. A real-time computer application for the prediction of fire spread across the Australian landscape. Simulation 1996, 67, 230–240. [Google Scholar] [CrossRef]
  5. Plourde, F.; Doan-Kim, S.; Dumas, J.; Malet, J. A new model of wildland fire simulation. Fire Saf. J. 1997, 29, 283–299. [Google Scholar] [CrossRef]
  6. Finney, M.A. FARSITE, Fire Area Simulator–Model Development and Evaluation; Number 4; US Department of Agriculture, Forest Service, Rocky Mountain Research Station: Ogden, UT, USA, 1998. [Google Scholar]
  7. Perry, G.L.; Sparrow, A.D.; Owens, I.F. A GIS-supported model for the simulation of the spatial structure of wildland fire, Cass Basin, New Zealand. J. Appl. Ecol. 1999, 36, 502–518. [Google Scholar] [CrossRef]
  8. Eklund, P. A distributed spatial architecture for bush fire simulation. Int. J. Geogr. Inf. Sci. 2001, 15, 363–378. [Google Scholar] [CrossRef]
  9. Lopes, A.; Cruz, M.; Viegas, D. FIRESTATION–An integrated system for the simulation of wind flow and fire spread over complex topography. In Proceedings of the III International Conference on Forest Fire Research and 14th Conference on Fire and Forest Meteorology, Luso, Portugal, 16–20 November 1998; pp. 16–20. [Google Scholar]
  10. Committee, C.S. Prometheus User Manual, v. 3.0. 1; Canadian Forest Service. Available online: https://prometheus.io/docs/introduction/overview/ (accessed on 24 May 2021).
  11. Miller, C.; Hilton, J.; Sullivan, A.; Prakash, M. SPARK–A bushfire spread prediction tool. In Proceedings of the International Symposium on Environmental Software Systems, Melbourne, VIC, Australia, 25–27 March 2015; Springer: Berlin, Germany, 2015; pp. 262–271. [Google Scholar]
  12. Tolhurst, K.; Shields, B.; Chong, D. Phoenix: Development and application of a bushfire risk management tool. Aust. J. Emerg. Manag. 2008, 23, 47–54. [Google Scholar]
  13. Matthews, S.; Fox-Hughes, P.; Grootemaat, S.; Hollis, J.J.; Kenny, B.J.; Sauvage, S. Australian Fire Danger Rating System: Research Prototype; International Association of Wildland Fire: Missoula, MT, USA, 2019; p. 384. [Google Scholar]
  14. KC, U.; Garg, S.; Hilton, J. An efficient framework for ensemble of natural disaster simulations as a service. Geosci. Front. 2020, 11, 1859–1873. [Google Scholar] [CrossRef]
  15. Sobol’, I.M. On sensitivity estimation for nonlinear mathematical models. Mat. Model. 1990, 2, 112–118. [Google Scholar]
  16. Cukier, R.; Fortuin, C.; Shuler, K.E.; Petschek, A.; Schaibly, J. Study of the sensitivity of coupled reaction systems to uncertainties in rate coefficients. I Theory. J. Chem. Phys. 1973, 59, 3873–3878. [Google Scholar] [CrossRef]
  17. Saltelli, A.; Tarantola, S.; Chan, K.S. A quantitative model-independent method for global sensitivity analysis of model output. Technometrics 1999, 41, 39–56. [Google Scholar] [CrossRef]
  18. Krzykacz-Hausmann, B. Epistemic sensitivity analysis based on the concept of entropy. In Proceedings of the SAMO2001, Madrid, Spain, 18–20 June 2001; pp. 31–35. [Google Scholar]
  19. Plischke, E.; Borgonovo, E.; Smith, C.L. Global sensitivity measures from given data. Eur. J. Oper. Res. 2013, 226, 536–550. [Google Scholar] [CrossRef]
  20. Pianosi, F.; Wagener, T. A simple and efficient method for global sensitivity analysis based on cumulative distribution functions. Environ. Model. Softw. 2015, 67, 1–11. [Google Scholar] [CrossRef] [Green Version]
  21. Freissinet, C.; Vauclin, M.; Erlich, M. Comparison of first-order analysis and fuzzy set approach for the evaluation of imprecision in a pesticide groundwater pollution screening model. J. Contam. Hydrol. 1999, 37, 21–43. [Google Scholar] [CrossRef]
  22. Richardson, A.D.; Hollinger, D.Y. Statistical modeling of ecosystem respiration using eddy covariance data: Maximum likelihood parameter estimation, and Monte Carlo simulation of model and parameter uncertainty, applied to three simple models. Agric. For. Meteorol. 2005, 131, 191–208. [Google Scholar] [CrossRef]
  23. Bachmann, A.; Allgöwer, B. Uncertainty propagation in wildland fire behaviour modelling. Int. J. Geogr. Inf. Sci. 2002, 16, 115–127. [Google Scholar] [CrossRef]
  24. Hilton, J.E.; Stephenson, A.G.; Huston, C.; Swedosh, W. Polynomial Chaos for sensitivity analysis in wildfire modelling. In Proceedings of the International Congress on Modelling and Simulation, Hobart, Australia, 3–8 December 2017; pp. 3–8. [Google Scholar]
  25. Wagener, T.; Kollat, J. Numerical and visual evaluation of hydrological and environmental models using the Monte Carlo analysis toolbox. Environ. Model. Softw. 2007, 22, 1021–1033. [Google Scholar] [CrossRef]
  26. Ekstrom, P.A. Eikos: A Simulation Toolbox for Sensitivity Analysis in Matlab. FACILIA AB. 2005. Available online: https://baixardoc.com/preview/eikos-a-simulation-toolbox-for-sensitivity-analysis-5c9a8be15eec9 (accessed on 12 March 2021).
  27. D’Augustine, A.F. MATLODE: A MATLAB ODE Solver and Sensitivity Analysis Toolbox. Ph.D. Thesis, Virginia Tech, Blacksburg, VA, USA, 2018. [Google Scholar]
  28. Pianosi, F.; Sarrazin, F.; Wagener, T. A Matlab toolbox for global sensitivity analysis. Environ. Model. Softw. 2015, 70, 80–85. [Google Scholar] [CrossRef] [Green Version]
  29. Herman, J.; Usher, W. SALib: An open-source Python library for sensitivity analysis. J. Open Source Softw. 2017, 2, 97. [Google Scholar] [CrossRef]
  30. Dutfoy, A.; Dutka-Malen, I.; Pasanisi, A.; Lebrun, R.; Mangeant, F.; Gupta, J.S.; Pendola, M.; Yalamas, T. OpenTURNS, an Open Source initiative to Treat Uncertainties, Risks’ N Statistics in a structured industrial approach. In Proceedings of the 41èmes Journées de Statistique, SFdS, Bordeaux, France, 25–29 May 2009. [Google Scholar]
  31. Preisler, H.K.; Chen, S.C.; Fujioka, F.; Benoit, J.W.; Westerling, A.L. Wildland fire probabilities estimated from weather model-deduced monthly mean fire danger indices. Int. J. Wildland Fire 2008, 17, 305–316. [Google Scholar] [CrossRef] [Green Version]
  32. Pimont, F.; Fargeon, H.; Opitz, T.; Ruffault, J.; Barbero, R.; Martin-StPaul, N.; Rigolot, E.; Rivière, M.; Dupuy, J.L. Prediction of regional wildfire activity in the probabilistic Bayesian framework of Firelihood. Ecol. Appl. 2021, 31, e02316. [Google Scholar] [CrossRef] [PubMed]
  33. Ager, A.A.; Barros, A.M.; Day, M.A.; Preisler, H.K.; Spies, T.A.; Bolte, J. Analyzing fine-scale spatiotemporal drivers of wildfire in a forest landscape model. Ecol. Model. 2018, 384, 87–102. [Google Scholar] [CrossRef]
  34. Preisler, H.K.; Westerling, A.L.; Gebert, K.M.; Munoz-Arriola, F.; Holmes, T.P. Spatially explicit forecasts of large wildland fire probability and suppression costs for California. Int. J. Wildland Fire 2011, 20, 508–517. [Google Scholar] [CrossRef]
  35. Penman, T.D.; Cirulis, B.; Marcot, B.G. Bayesian decision network modeling for environmental risk management: A wildfire case study. J. Environ. Manag. 2020, 270, 110735. [Google Scholar] [CrossRef]
  36. Mendes, J.M.; de Zea Bermudez, P.C.; Pereira, J.; Turkman, K.; Vasconcelos, M. Spatial extremes of wildfire sizes: Bayesian hierarchical models for extremes. Environ. Ecol. Stat. 2010, 17, 1–28. [Google Scholar] [CrossRef]
  37. Joseph, M.B.; Rossi, M.W.; Mietkiewicz, N.P.; Mahood, A.L.; Cattau, M.E.; St. Denis, L.A.; Nagy, R.C.; Iglesias, V.; Abatzoglou, J.T.; Balch, J.K. Spatiotemporal prediction of wildfire size extremes with Bayesian finite sample maxima. Ecol. Appl. 2019, 29, e01898. [Google Scholar] [CrossRef] [Green Version]
  38. Cisneros, D.; Gong, Y.; Yadav, R.; Hazra, A.; Huser, R. A combined statistical and machine learning approach for spatial prediction of extreme wildfire frequencies and sizes. arXiv 2021, arXiv:2112.14920. [Google Scholar] [CrossRef]
  39. KC, U.; Hilton, J.; Garg, S.; Aryal, J. A probability-based risk metric for operational wildfire risk management. Environ. Model. Softw. 2022, 148, 105286. [Google Scholar] [CrossRef]
  40. KC, U.; Aryal, J. Leveraging a wildfire risk prediction metric with spatial clustering. Fire 2022, 5, 213. [Google Scholar] [CrossRef]
  41. KC, U.; Garg, S.; Hilton, J.; Aryal, J. An adaptive quadtree-based approach for efficient decision making in wildfire risk assessment. Environ. Model. Softw. 2022, 160, 105590. [Google Scholar] [CrossRef]
  42. Carriger, J.F.; Thompson, M.; Barron, M.G. Causal Bayesian networks in assessments of wildfire risks: Opportunities for ecological risk assessment and management. Integr. Environ. Assess. Manag. 2021, 17, 1168–1178. [Google Scholar] [CrossRef] [PubMed]
  43. Khakzad, N. Modeling wildfire spread in wildland-industrial interfaces using dynamic Bayesian network. Reliab. Eng. Syst. Saf. 2019, 189, 165–176. [Google Scholar] [CrossRef]
  44. Storey, M.A.; Bedward, M.; Price, O.F.; Bradstock, R.A.; Sharples, J.J. Derivation of a Bayesian fire spread model using large-scale wildfire observations. Environ. Model. Softw. 2021, 144, 105127. [Google Scholar] [CrossRef]
  45. Jaafari, A.; Gholami, D.M.; Zenner, E.K. A Bayesian modeling of wildfire probability in the Zagros Mountains, Iran. Ecol. Inform. 2017, 39, 32–44. [Google Scholar] [CrossRef]
  46. Silva, G.L.; Soares, P.; Marques, S.; Dias, M.I.; Oliveira, M.M.; Borges, J.G. A Bayesian modelling of wildfires in Portugal. In Dynamics, Games and Science; Springer: Berlin, Germany, 2015; pp. 723–733. [Google Scholar]
  47. Zwirglmaier, K.; Papakosta, P.; Straub, D. Learning a Bayesian network model for predicting wildfire behavior. In Proceedings of the ICOSSAR 2013, New York, NY, USA, 16–20 June 2013. [Google Scholar]
  48. KC, U.; Garg, S.; Hilton, J.; Aryal, J. A cloud-based framework for sensitivity analysis of natural hazard models. Environ. Model. Softw. 2020, 134, 104800. [Google Scholar] [CrossRef]
  49. KC, U.; Aryal, J.; Hilton, J.; Garg, S. A Surrogate Model for Rapidly Assessing the Size of a Wildfire over Time. Fire 2021, 4, 20. [Google Scholar] [CrossRef]
  50. Lindgren, F.; Rue, H. Bayesian spatial modelling with R-INLA. J. Stat. Softw. 2015, 63, 1–25. [Google Scholar] [CrossRef] [Green Version]
  51. Gelfand, A.E.; Smith, A.F. Sampling-based approaches to calculating marginal densities. J. Am. Stat. Assoc. 1990, 85, 398–409. [Google Scholar] [CrossRef]
  52. Simpson, D.; Rue, H.; Riebler, A.; Martins, T.G.; Sørbye, S.H. Penalising model component complexity: A principled, practical approach to constructing priors. Stat. Sci. 2017, 32, 1–28. [Google Scholar] [CrossRef]
  53. Chib, S. Marginal likelihood from the Gibbs output. J. Am. Stat. Assoc. 1995, 90, 1313–1321. [Google Scholar] [CrossRef]
  54. Spiegelhalter, D.J.; Best, N.G.; Carlin, B.P.; Van Der Linde, A. Bayesian measures of model complexity and fit. J. R. Stat. Soc. Ser. B (Stat. Methodol.) 2002, 64, 583–639. [Google Scholar] [CrossRef] [Green Version]
  55. Watanabe, S. A widely applicable Bayesian information criterion. J. Mach. Learn. Res. 2013, 14, 867–897. [Google Scholar]
  56. Pettit, L. The conditional predictive ordinate for the normal distribution. J. R. Stat. Soc. Ser. B (Methodol.) 1990, 52, 175–184. [Google Scholar] [CrossRef]
  57. Spark: Predicting Bushfire Spread. Available online: https://data61.csiro.au/en/Our-Research/Our-Work/Safety-and-Security/Disaster-Management/Spark (accessed on 12 May 2021).
  58. Tasmania. List Data. 2021. Available online: https://listdata.thelist.tas.gov.au/opendata/ (accessed on 12 March 2021).
  59. Tasmania Fire Service. State Fire Commission Annual Report; Tasmania Fire Service: Hobart, Australia, 2019. [Google Scholar]
  60. Ujjwal, K.C.; Garg, S.; Hilton, J.; Aryal, J. Fire Simulation Data Set for Tasmania; CSIRO: Canberra, Australia, 2021. [Google Scholar] [CrossRef]
  61. Saltelli, A. Making best use of model evaluations to compute sensitivity indices. Comput. Phys. Commun. 2002, 145, 280–297. [Google Scholar] [CrossRef]
  62. Benesty, J.; Chen, J.; Huang, Y.; Cohen, I. Pearson correlation coefficient. In Noise Reduction in Speech Processing; Springer: Berlin, Germany, 2009; pp. 1–4. [Google Scholar]
Figure 1. Workflow Diagram. Before investigating the effects of a latent process, the entire dataset is used for Bayesian model fitting. After evaluating the performance of the model for six different priors of the precision parameter, the best model is selected. For sensitivity analysis and scenario analysis in wildfire management, we use a combination of true observations and the best model’s predicted values.
Figure 1. Workflow Diagram. Before investigating the effects of a latent process, the entire dataset is used for Bayesian model fitting. After evaluating the performance of the model for six different priors of the precision parameter, the best model is selected. For sensitivity analysis and scenario analysis in wildfire management, we use a combination of true observations and the best model’s predicted values.
Atmosphere 14 00559 g001
Figure 2. Posterior marginals of three input parameters—temperature, relative humidity, and wind speed (left to right) with different priors of precision. The choice of priors in Bayesian models can affect the posterior marginals of model parameters. Most of the priors produced slightly different but similar posteriors, with the exception of the half-Cauchy prior, which resulted in significantly displaced posteriors. The negative values for relative humidity indicate that higher relative humidity values are associated with a decrease in fire size.
Figure 2. Posterior marginals of three input parameters—temperature, relative humidity, and wind speed (left to right) with different priors of precision. The choice of priors in Bayesian models can affect the posterior marginals of model parameters. Most of the priors produced slightly different but similar posteriors, with the exception of the half-Cauchy prior, which resulted in significantly displaced posteriors. The negative values for relative humidity indicate that higher relative humidity values are associated with a decrease in fire size.
Atmosphere 14 00559 g002
Figure 3. Preference order (1 is the highest and 6 is the lowest) of models with different priors based on several evaluation metrics—DIC, WAIC, MLK, and CPO. The values for DIC and CPO are the same, so they are represented by a single graph in the plot. Overall, the penalized complexity (PC) prior (default loggamma) had better model performance, even though it had the worst value for MLK. The Cauchy prior had the overall worst model performance.
Figure 3. Preference order (1 is the highest and 6 is the lowest) of models with different priors based on several evaluation metrics—DIC, WAIC, MLK, and CPO. The values for DIC and CPO are the same, so they are represented by a single graph in the plot. Overall, the penalized complexity (PC) prior (default loggamma) had better model performance, even though it had the worst value for MLK. The Cauchy prior had the overall worst model performance.
Atmosphere 14 00559 g003
Figure 4. Measure of the similarity between the true observations and the Bayesian model predicted values. The similarity increased as more true observed values were considered in the Bayesian model (sample size). However, the improvement in similarity was not significant even when the sample size increased from 4000 to 7000.
Figure 4. Measure of the similarity between the true observations and the Bayesian model predicted values. The similarity increased as more true observed values were considered in the Bayesian model (sample size). However, the improvement in similarity was not significant even when the sample size increased from 4000 to 7000.
Atmosphere 14 00559 g004
Figure 5. Sensitivity indices of fire size to input weather inputs estimated with true observations and Bayesian model estimated values. The sample size of 8000 represents the complete dataset with all true observations. Our application of Bayesian model-predicted values gave close-to-true estimations of the sensitivity indices using as few as 4000 true observations. As the number of true observations used in the estimation increased, the estimations improved in their closeness to the complete dataset of true observations.
Figure 5. Sensitivity indices of fire size to input weather inputs estimated with true observations and Bayesian model estimated values. The sample size of 8000 represents the complete dataset with all true observations. Our application of Bayesian model-predicted values gave close-to-true estimations of the sensitivity indices using as few as 4000 true observations. As the number of true observations used in the estimation increased, the estimations improved in their closeness to the complete dataset of true observations.
Atmosphere 14 00559 g005
Table 1. Different Priors for the precision parameter used for model fitting.
Table 1. Different Priors for the precision parameter used for model fitting.
PriorsParameters NameParameters Value
half-Cauchymean, scale(0, 25)
half-tmean, shape(0, 3)
log-gammashape, rate(1, 0.00005)
half-normalmean, precision(0, 0.001)
PC log-gammashape, rate(5, 0.01)
uniform improperstandard deviation ( σ ) π U N ( σ ) 1
Table 2. Range and Probability Density Function (PDF) of three input parameters—temperature, relative humidity, and wind speed. The ranges and the distributions of the input parameters are chosen to cover operational weather conditions for wildfire modeling in the Australian scenario as per the experimental setup of our previous research.
Table 2. Range and Probability Density Function (PDF) of three input parameters—temperature, relative humidity, and wind speed. The ranges and the distributions of the input parameters are chosen to cover operational weather conditions for wildfire modeling in the Australian scenario as per the experimental setup of our previous research.
ParametersUnitPdfRange
Temperature CUniform Distribution[10, 40]
Relative Humidity%Uniform Distribution[5, 90]
Wind Speed kmh 1 Uniform Distribution[10, 60]
Table 3. Comparative Analysis of Bayesian models. All the evaluation metrics favor the model with random noise as the model had marginally improved performance than the model without random noise. Due to insignificant differences between the values of the metrics, both models may be used with proper priors.
Table 3. Comparative Analysis of Bayesian models. All the evaluation metrics favor the model with random noise as the model had marginally improved performance than the model without random noise. Due to insignificant differences between the values of the metrics, both models may be used with proper priors.
ComponentsWith Latent ProcessWithout Latent Process
MLK64,329.6064,561.05
DIC128,586.56128,945.93
WAIC128,584.76129,008.34
C P O * 8.048.063
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

KC, U.; Aryal, J.; Bakar, K.S.; Hilton, J.; Buyya, R. Applying Bayesian Models to Reduce Computational Requirements of Wildfire Sensitivity Analyses. Atmosphere 2023, 14, 559. https://doi.org/10.3390/atmos14030559

AMA Style

KC U, Aryal J, Bakar KS, Hilton J, Buyya R. Applying Bayesian Models to Reduce Computational Requirements of Wildfire Sensitivity Analyses. Atmosphere. 2023; 14(3):559. https://doi.org/10.3390/atmos14030559

Chicago/Turabian Style

KC, Ujjwal, Jagannath Aryal, K. Shuvo Bakar, James Hilton, and Rajkumar Buyya. 2023. "Applying Bayesian Models to Reduce Computational Requirements of Wildfire Sensitivity Analyses" Atmosphere 14, no. 3: 559. https://doi.org/10.3390/atmos14030559

APA Style

KC, U., Aryal, J., Bakar, K. S., Hilton, J., & Buyya, R. (2023). Applying Bayesian Models to Reduce Computational Requirements of Wildfire Sensitivity Analyses. Atmosphere, 14(3), 559. https://doi.org/10.3390/atmos14030559

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop