Front shape similarity measure for shape-oriented sensitivity analysis and data assimilation for Eikonal equation

We present a shape-oriented data assimilation strategy suitable for front-tracking problems through the example of wildfire. The concept of " front " is used to model, at regional scales, the burning area delimitation that moves, undergoes shape and topological changes under heterogeneous orography, biomass fuel and micrometeorology. The simulation-observation discrepancies are represented using a front shape similarity measure deriving from image processing and based on the Chan-Vese contour fitting functional. We show that consistent corrections of the front location and uncertain physical parameters can be obtained using this measure applied on a level-set fire growth model solving for an eikonal equation. This study involves a Luenberger observer for state estimation, including a topological gradient term to track multiple fronts, and of a reduced-order Kalman filter for joint parameter estimation. We also highlight the need – prior to parameter estimation – for sensitivity analysis based on the same discrepancy measure, and for instance using polynomial chaos metamodels, to ensure a meaningful inverse solution is achieved. The performance of the shape-oriented data assimilation strategy is assessed on a synthetic configuration subject to uncertainties in front initial position, near-surface wind magnitude and direction. The use of a robust front shape similarity measure paves the way toward the direct assimilation of infrared images and is a valuable asset in the perspective of data-driven wildfire modeling.

[1]  P. Moireau,et al.  The mechanism of monomer transfer between two structurally distinct PrP oligomers , 2017, PloS one.

[2]  M. Finney FARSITE : Fire Area Simulator : model development and evaluation , 1998 .

[3]  Arnaud Trouvé,et al.  Towards predictive data-driven simulations of wildfire spread - Part II: Ensemble Kalman Filter for the state estimation of a front-tracking simulator of wildfire spread , 2014 .

[4]  Guillaume Charpiat,et al.  Converting Level Set Gradients to Shape Gradients , 2010, ECCV.

[5]  Dominique Chapelle,et al.  A Luenberger observer for reaction-diffusion models with front position data , 2015, J. Comput. Phys..

[6]  M. Yousuff Hussaini,et al.  Quantifying parametric uncertainty in the Rothermel model , 2008 .

[7]  Ronald G. Rehm Fire-Front Propagation Using the Level Set Method , 2019 .

[8]  Nicolas Papadakis,et al.  A Variational Technique for Time Consistent Tracking of Curves and Motion , 2008, Journal of Mathematical Imaging and Vision.

[9]  Didier Lucor,et al.  Interactive comment on “Towards predictive data-driven simulations of wildfire spread – Part I: Reduced-cost Ensemble Kalman Filter based on a Polynomial Chaos surrogate model for parameter estimation” , 2014 .

[10]  Dominique Chapelle,et al.  Fundamental principles of data assimilation underlying the Verdandi library: applications to biophysical model personalization within euHeart , 2013, Medical & Biological Engineering & Computing.

[11]  Sai Ravela,et al.  Amplitude-Position Formulation of Data Assimilation , 2006, International Conference on Computational Science.

[12]  P. Tallec,et al.  Joint state and parameter estimation for distributed mechanical systems , 2008 .

[13]  Bruno Sudret,et al.  Global sensitivity analysis using polynomial chaos expansions , 2008, Reliab. Eng. Syst. Saf..

[14]  Chad M. Hoffman,et al.  A comparison of level set and marker methods for the simulation of wildland fire front propagation , 2016 .

[15]  R. Reichle Data assimilation methods in the Earth sciences , 2008 .

[16]  Takemasa Miyoshi,et al.  Estimating Model Parameters with Ensemble-Based Data Assimilation: A Review , 2013 .

[17]  L. Vese,et al.  A level set algorithm for minimizing the Mumford-Shah functional in image processing , 2001, Proceedings IEEE Workshop on Variational and Level Set Methods in Computer Vision.

[18]  James A. Sethian,et al.  Theory, algorithms, and applications of level set methods for propagating interfaces , 1996, Acta Numerica.

[19]  Patricia L. Andrews,et al.  Examination of the wind speed limit function in the Rothermel surface fire spread model , 2013 .

[20]  D. Mumford,et al.  Optimal approximations by piecewise smooth functions and associated variational problems , 1989 .

[21]  B. Fisher,et al.  The product of distributions , 1971 .

[22]  Soroosh Sorooshian,et al.  Dual state-parameter estimation of hydrological models using ensemble Kalman filter , 2005 .

[23]  D. Xiu Numerical Methods for Stochastic Computations: A Spectral Method Approach , 2010 .

[24]  Isabelle Bloch,et al.  Fuzzy sets for image processing and understanding , 2015, Fuzzy Sets Syst..

[25]  James A. Sethian,et al.  The Fast Construction of Extension Velocities in Level Set Methods , 1999 .

[26]  Nicolas Papadakis,et al.  Active contours without level sets , 2012, 2012 19th IEEE International Conference on Image Processing.

[27]  Dominique Chapelle,et al.  Sequential State Estimation for Electrophysiology Models with Front Level-Set Data Using Topological Gradient Derivations , 2015, FIMH.

[28]  M. L. Martins,et al.  Multiscale models for the growth of avascular tumors , 2007 .

[29]  B. Engquist,et al.  Discretization of Dirac delta functions in level set methods , 2005 .

[30]  A. Sullivan,et al.  Wildland surface fire spread modelling, 1990–2007. 2: Empirical and quasi-empirical models , 2007, 0706.4128.

[31]  Lin He,et al.  Incorporating topological derivatives into shape derivatives based level set methods , 2007, J. Comput. Phys..

[32]  Yunmei Chen,et al.  Region Based Image Segmentation Using a Modified Mumford-Shah Algorithm , 2007, SSVM.

[33]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods , 1999 .

[34]  A. Jazwinski Stochastic Processes and Filtering Theory , 1970 .

[35]  Loic Le Gratiet,et al.  Metamodel-based sensitivity analysis: polynomial chaos expansions and Gaussian processes , 2016, 1606.04273.

[36]  Heinz Pitsch,et al.  A consistent level set formulation for large-eddy simulation of premixed turbulent combustion , 2005 .