Adam K. Kochanski

Coupled atmosphere-wildland fire modeling with WRF 3.3 and SFIRE 2011

Abstract. We describe the physical model, numerical algorithms, and software structure of a model consisting of the Weather Research and Forecasting (WRF) model, coupled with the fire-spread model (SFIRE) module. In every time step, the fire model inputs the surface wind, which drives the fire, and outputs the heat flux from the fire into the atmosphere, which in turn influences the atmosphere. SFIRE is implemented by the level set method, which allows a submesh representation of the burning region and a flexible implementation of various kinds of ignition. The coupled model is capable of running on a cluster faster than real time even with fine resolution in dekameters. It is available as a part of the Open Wildland Fire Modeling (OpenWFM) environment at http://openwfm.org , which contains also utilities for visualization, diagnostics, and data processing, including an extended version of the WRF Preprocessing System (WPS). The SFIRE code with a subset of the features is distributed with WRF 3.3 as WRF-Fire.

Temperature and Water Vapor Variance Scaling in Global Models: Comparisons to Satellite and Aircraft Data

AbstractObservations of the scale dependence of height-resolved temperature T and water vapor q variability are valuable for improved subgrid-scale climate model parameterizations and model evaluation. Variance spectral benchmarks for T and q obtained from the Atmospheric Infrared Sounder (AIRS) are compared to those generated by state-of-the-art numerical weather prediction “analyses” and “free-running” climate model simulations with spatial resolution comparable to AIRS. The T and q spectra from both types of models are generally too steep, with small-scale variance up to several factors smaller than AIRS. However, the two model analyses more closely resemble AIRS than the two free-running model simulations. Scaling exponents obtained for AIRS column water vapor (CWV) and height-resolved layers of q are also compared to the superparameterized Community Atmospheric Model (SP-CAM), highlighting large differences in the magnitude of CWV variance and the relative flatness of height-resolved q scaling in SP-...

Real time simulation of 2007 Santa Ana fires

Abstract In this study we test the feasibility of using a coupled atmosphere–fire model for real time simulations of massive fires. A physics-based coupled atmosphere–fire model is used to resolve the large-scale and local weather as well as the atmosphere–fire interactions, while combustion is represented simply using an existing operational surface fire behavior model. This model combination strikes a balance between fidelity and speed of execution. The feasibility of this approach is examined based on an analysis of a numerical simulation of two very large Santa Ana fires using WRF–Sfire, a coupled atmosphere–fire model developed by the Open Wild Fire Modeling Community (OpenWFM.org). The study demonstrates that a wind and fire spread forecast of reasonable accuracy was obtained at an execution speed that would have made real-time wildfire forecasting of this event possible.

Assimilation of Perimeter Data and Coupling with Fuel Moisture in a Wildland Fire - Atmosphere DDDAS

We present a methodology to change the state of the Weather Research Forecasting (WRF) model coupled with the fire spread code SFIRE, based on Rothermel’s formula and the level set method, and with a fuel moisture model. The fire perimeter in the model changes in response to data while the model is running. However, the atmosphere state takes time to develop in response to the forcing by the heat flux from the fire. Therefore, an artificial fire history is created from an earlier fire perimeter to the new perimeter, and replayed with the proper heat fluxes to allow the atmosphere state to adjust. The method is an extension of an earlier method to start the coupled fire model from a developed fire perimeter rather than an ignition point. The level set method can be also used to identify parameters of the simulation, such as the fire spread rate. The coupled model is available from openwfm.org, and it extends the WRF-Fire code in WRF release.

Simulation of the 2009 harmanli fire (bulgaria)

We use a coupled atmosphere-fire model to simulate a fire that occurred on August 14---17, 2009, in the Harmanli region, Bulgaria. Data was obtained from GIS and satellites imagery, and from standard atmospheric data sources. Fuel data was classified in the 13 Anderson categories. For correct fire behavior, the spatial resolution of the models needed to be fine enough to resolve the essential micrometeorological effects. The simulation results are compared to available incident data. The code runs faster than real time on a cluster. The model is available from openwfm.org and it extends WRF-Fire from WRF 3.3 release.

A slab model of the Great Salt Lake for regional climate simulation

A slab lake model was developed for the Great Salt Lake (GSL) and coupled to a regional climate model to enable better evaluation of regional effects of projected climate change. The GSL is hypersaline with an area of approximately 4400 km2, and its notable shallowness (the deeper sections average 6.5–9 m at current lake levels) renders it highly sensitive to climate change. A time-independent (constant) effective mixing depth of approximately 5 m was determined for the GSL by numerically optimizing model-observation agreement, and improvement gained using a time-dependent effective mixing depth assumption was smaller than the uncertainty in the satellite-based observations. The slab model with constant effective mixing depth accounted for more than 97% of the variance in satellite-based observations of GSL surface temperature for years 2001 through 2003. Using a lake surface temperature climatology in place of the lake model resulted in annual mean near-surface air temperature differences that were small (∼10−2 K) away from the lake, but differences in annual precipitation downstream reached 3 cm (4.5%) mainly because of enhanced turbulent heat fluxes off the lake during spring. When subjected to a range of pseudo global warming scenarios, the annual mean lake surface temperature increased by 0.8°C per degree of air temperature increase.

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

In this study, WRF-Sfire is coupled with WRF-Chem to construct WRFSC, an integrated forecast system for wildfire behaviour and smoke prediction. WRF-Sfire directly predicts wildfire spread, plume and plume-top heights, providing comprehensive meteorology and fire emissions to chemical transport model WRF-Chem, eliminating the need for an external plume-rise model. Evaluation of WRFSC was based on comparisons between available observations of fire perimeter and fire intensity, smoke spread, PM2.5 (particulate matter less than 2.5 μm in diameter), NO and ozone concentrations, and plume-top heights with the results of two WRFSC simulations, a 48-h simulation of the 2007 Witch–Guejito Santa Ana fires and a 96-h WRF-Sfire simulation with passive tracers of the 2012 Barker Canyon fire. The study found overall good agreement between forecast and observed local- and long-range fire spread and smoke transport for the Witch–Guejito fire. However, ozone, PM2.5 and NO concentrations were generally underestimated and peaks mistimed in the simulations. This study found overall good agreement between simulated and observed plume-top heights, with slight underestimation by the simulations. Two promising results were the agreement between plume-top heights for the Barker Canyon fire and faster than real-time execution, making WRFSC a possible operational tool.

One-Way Coupling of the WRF–QUIC Urban Dispersion Modeling System

AbstractSimulations of local weather and air quality in urban areas must account for processes spanning from meso- to microscales, including turbulence and transport within the urban canopy layer. Here, the authors investigate the performance of the building-resolving Quick Urban Industrial Complex (QUIC) Dispersion Modeling System driven with mean wind profiles from the mesoscale Weather Research and Forecasting (WRF) Model. Dispersion simulations are performed for intensive observation periods 2 and 8 of the Joint Urban 2003 field experiment conducted in Oklahoma City, Oklahoma, using an ensemble of expert-derived wind profiles from observational data as well as profiles derived from WRF runs. The results suggest that WRF can be used successfully as a source of inflow boundary conditions for urban simulations, without the collection and processing of intensive field observations needed to produce expert-derived wind profiles. Detailed statistical analysis of tracer concentration fields suggests that, fo...

Data management and analysis with WRF and SFIRE

We introduce several useful utilities in development for the creation and analysis of real wildland fire simulations using WRF and SFIRE. These utilities exist as standalone programs and scripts as well as extensions to other well known software. Python web scrapers automate the process of downloading and preprocessing atmospheric and surface data from common sources. Other scripts simplify the domain setup by creating parameter files automatically. Integration with Google Earth allows users to explore the simulation in a 3D environment along with real surface imagery. Postprocessing scripts provide the user with a number of output data formats compatible with many commonly used visualization suites allowing for the creation of high quality 3D renderings. As a whole, these improvements build toward a unified web application that brings a sophisticated wildland fire modeling environment to scientists and users alike.

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

Fuel moisture has a major influence on the behaviour of wildland fires and is an important underlying factor in fire risk assessment. We propose a method to assimilate dead fuel moisture content (FMC) 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 FMC. 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-h FMC observations from Colorado RAWS in 2013. Using leave-one-out testing we show that the TSM compares favourably 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 on FMC estimates in unobserved fuel classes.