Spectral diagonal ensemble Kalman filters

A new type of ensemble Kalman filter is developed, which is based on replacing the sample covariance in the analysis step by its diagonal in a spectral basis. It is proved that this technique improves the aproximation of the covariance when the covariance itself is diagonal in the spectral basis, as is the case, e.g., for a second-order stationary random field and the Fourier basis. The method is extended by wavelets to the case when the state variables are random fields, which are not spatially homogeneous. Efficient implementations by the fast Fourier transform (FFT) and discrete wavelet transform (DWT) are presented for several types of observations, including high-dimensional data given on a part of the domain, such as radar and satellite images. Computational experiments confirm that the method performs well on the Lorenz 96 problem and the shallow water equations with very small ensembles and over multiple analysis cycles.

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.

Convergence of the Square Root Ensemble Kalman Filter in the Large Ensemble Limit

Ensemble filters implement sequential Bayesian estimation by representing the probability distribution by an ensemble mean and covariance. Unbiased square root ensemble filters use deterministic algorithms to produce an analysis (posterior) ensemble with prescribed mean and covariance, consistent with the Kalman update. This includes several filters used in practice, such as the Ensemble Transform Kalman Filter (ETKF), the Ensemble Adjustment Kalman Filter (EAKF), and a filter by Whitaker and Hamill. We show that at every time index, as the number of ensemble members increases to infinity, the mean and covariance of an unbiased ensemble square root filter converge to those of the Kalman filter, in the case a linear model and an initial distribution of which all moments exist. The convergence is in $L^{p}$ and the convergence rate does not depend on the model dimension. The result holds in the infinitely dimensional Hilbert space as well.

BDDC by a frontal solver and the stress computation in a hip joint replacement

A parallel implementation of the BDDC method using the frontal solver is employed to solve systems of linear equations from finite element analysis, and incorporated into a standard finite element system for engineering analysis by linear elasticity. Results of computation of stress in a hip replacement are presented. The part is made of titanium and loaded by the weight of human body. The performance of BDDC with added constraints by averages and with added corners is compared.