This paper is the first part in a series of two articles and presents a data-driven wildfire simulator for forecasting wildfire spread scenarios, at a reduced computational cost that is consistent with operational systems. The prototype simulator features the following components: an Eulerian front propagation solver FIREFLY that adopts a regional-scale modeling viewpoint, treats wildfires as surface propagating fronts, and uses a description of the local rate of fire spread (ROS) as a function of environmental conditions based on Rothermel's model; a series of airborne-like observations of the fire front positions; and a data assimilation (DA) algorithm based on an ensemble Kalman filter (EnKF) for parameter estimation. This stochastic algorithm partly accounts for the nonlinearities between the input parameters of the semi-empirical ROS model and the fire front position, and is sequentially applied to provide a spatially uniform correction to wind and biomass fuel parameters as observations become available. A wildfire spread simulator combined with an ensemble-based DA algorithm is therefore a promising approach to reduce uncertainties in the forecast position of the fire front and to introduce a paradigm-shift in the wildfire emergency response. In order to reduce the computational cost of the EnKF algorithm, a surrogate model based on a polynomial chaos (PC) expansion is used in place of the forward model FIREFLY in the resulting hybrid PC-EnKF algorithm. The performance of EnKF and PC-EnKF is assessed on synthetically generated simple configurations of fire spread to provide valuable information and insight on the benefits of the PC-EnKF approach, as well as on a controlled grassland fire experiment. The results indicate that the proposed PC-EnKF algorithm features similar performance to the standard EnKF algorithm, but at a much reduced computational cost. In particular, the re-analysis and forecast skills of DA strongly relate to the spatial and temporal variability of the errors in the ROS model parameters.

This paper is the second part in a series of two articles, which aims at presenting a data-driven modeling strategy for forecasting wildfire spread scenarios based on the assimilation of the observed fire front location and on the sequential correction of model parameters or model state. This model relies on an estimation of the local rate of fire spread (ROS) as a function of environmental conditions based on Rothermel's semi-empirical formulation, in order to propagate the fire front with an Eulerian front-tracking simulator. In Part I, a data assimilation (DA) system based on an ensemble Kalman filter (EnKF) was implemented to provide a spatially uniform correction of biomass fuel and wind parameters and thereby, produce an improved forecast of the wildfire behavior (addressing uncertainties in the input parameters of the ROS model only). In Part II, the objective of the EnKF algorithm is to sequentially update the two-dimensional coordinates of the markers along the discretized fire front, in order to provide a spatially distributed correction of the fire front location and thereby, a more reliable initial condition for further model time-integration (addressing all sources of uncertainties in the ROS model). The resulting prototype data-driven wildfire spread simulator is first evaluated in a series of verification tests using synthetically generated observations; tests include representative cases with spatially varying biomass properties and temporally varying wind conditions. In order to properly account for uncertainties during the EnKF update step and to accurately represent error correlations along the fireline, it is shown that members of the EnKF ensemble must be generated through variations in estimates of the fire's initial location as well as through variations in the parameters of the ROS model. The performance of the prototype simulator based on state estimation (SE) or parameter estimation (PE) is then evaluated by comparison with data taken from a reduced-scale controlled grassland fire experiment. Results indicate that data-driven simulations are capable of correcting inaccurate predictions of the fire front location and of subsequently providing an optimized forecast of the wildfire behavior at future lead times. The complementary benefits of both PE and SE approaches, in terms of analysis and forecast performance, are also emphasized. In particular, it is found that the size of the assimilation window must be specified adequately with the persistence of the model initial condition and/or with the temporal and spatial variability of the environmental conditions in order to track sudden changes in wildfire behavior. The present prototype data-driven forecast system is still at an early stage of development. In this regard, this preliminary investigation provides valuable information on how to combine observations with a fire spread model in an efficient way, as well as guidelines to design the future system evolution in order to meet the operational requirements of wildfire spread monitoring.

Data assimilation is widely used to improve flood forecasting capability, especially through parameter inference requiring statistical information on the uncertain input parameters (upstream discharge, friction coefficient) as well as on the variability of the water level and its sensitivity with respect to the inputs. For particle filter or ensemble Kalman filter, stochastically estimating probability density function and covariance matrices from a Monte Carlo random sampling requires a large ensemble of model evaluations, limiting their use in real-time application. To tackle this issue, fast surrogate models based on polynomial chaos and Gaussian process can be used to represent the spatially distributed water level in place of solving the shallow water equations. This study investigates the use of these surrogates to estimate probability density functions and covariance matrices at a reduced computational cost and without the loss of accuracy, in the perspective of ensemble-based data assimilation. This study focuses on 1-D steady state flow simulated with MASCARET over the Garonne River (South-West France). Results show that both surrogates feature similar performance to the Monte-Carlo random sampling, but for a much smaller computational budget; a few MASCARET simulations (on the order of 10–100) are sufficient to accurately retrieve covariance matrices and probability density functions all along the river, even where the flow dynamic is more complex due to heterogeneous bathymetry. This paves the way for the design of surrogate strategies suitable for representing unsteady open-channel flows in data assimilation.