Data Assimilation of Satellite Fire Detection in Coupled Atmosphere-Fire Simulation by WRF-SFIRE

Currently available satellite active fire detection products from the VIIRS and MODIS instruments on polar-orbiting satellites produce detection squares in arbitrary locations. There is no global fire/no fire map, no detection under cloud cover, false negatives are common, and the detection squares are much coarser than the resolution of a fire behavior model. Consequently, current active fire satellite detection products should be used to improve fire modeling in a statistical sense only, rather than as a direct input. We describe a new data assimilation method for active fire detection, based on a modification of the fire arrival time to simultaneously minimize the difference from the forecast fire arrival time and maximize the likelihood of the fire detection data. This method is inspired by contour detection methods used in computer vision, and it can be cast as a Bayesian inverse problem technique, or a generalized Tikhonov regularization. After the new fire arrival time on the whole simulation domain is found, the model can be re-run from a time in the past using the new fire arrival time to generate the heat fluxes and to spin up the atmospheric model until the satellite overpass time, when the coupled simulation continues from the modified state.

Data assimilation of dead fuel moisture observations from remote automated weather stations

Fuel moisture has a major influence on the behavior of wildland fires and is an important underlying factor in fire risk assessment. We propose a method to assimilate dead fuel moisture content observations from remote automated weather stations (RAWS) into a time-lag fuel moisture model. RAWS are spatially sparse and a mechanism is needed to estimate fuel moisture content at locations potentially distant from observational stations. This is arranged using a trend surface model (TSM), which allows us to account for the effects of topography and atmospheric state on the spatial variability of fuel moisture content. At each location of interest, the TSM provides a pseudo-observation, which is assimilated via Kalman filtering. The method is tested with the time-lag fuel moisture model in the coupled weather-fire code WRF-SFIRE on 10-hr fuel moisture content observations from Colorado RAWS in 2013. We show using leave-one-out testing that the TSM compares favorably with inverse squared distance interpolation as used in the Wildland Fire Assessment System. Finally, we demonstrate that the data assimilation method is able to improve fuel moisture content estimates in unobserved fuel classes.

Toward an integrated system for fire, smoke, and air quality simulations

In this study, we describe how WRF-Sfire is coupled with WRF-Chem to construct WRFSC, an integrated forecast system for wildfire and smoke prediction. The integrated forecast system has the advantage of not requiring a simple plume-rise model and assumptions about the size and heat release from the fire in order to determine fire emissions into the atmosphere. With WRF-Sfire, wildfire spread, plume and plume-top heights are predicted directly, at every WRF timestep, providing comprehensive meteorology and fire emissions to the chemical transport model WRF-Chem. Evaluation of WRFSC was based on comparisons between available observations to the results of two WRFSC simulations. The study found overall good agreement between forecasted and observed fire spread and smoke transport for the Witch-Guejito fire. Also the simulated PM2.5 (fine particulate matter) peak concentrations matched the observations. However, the NO and ozone levels were underestimated in the simulations and the peak concentrations were mistimed. Determining the terminal or plume-top height is one of the most important aspects of simulating wildfire plume transport, and the study found overall good agreement between simulated and observed plume-top heights, with some (10% or less) underestimation by the simulations. One of the most promising results of the study was the agreement between passive-tracer modeled plume-top heights for the Barker Canyon fire simulation and observations. This simulation took only 13h, with the first 24h forecast ready in almost 3h, making it a possible operational tool for providing emission profiles for external chemical transport models.