Next Article in Journal
Point Cloud Based Mapping of Understory Shrub Fuel Distribution, Estimation of Fuel Consumption and Relationship to Pyrolysis Gas Emissions on Experimental Prescribed Burns
Next Article in Special Issue
Using Structure Location Data to Map the Wildland–Urban Interface in Montana, USA
Previous Article in Journal
The Behaviour of Water-Mists in Hot Air Induced by a Room Fire: Effect of the Initial Size of Droplets
Previous Article in Special Issue
Human Fire Use and Management: A Global Database of Anthropogenic Fire Impacts for Modelling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Construction of Probabilistic Wildfire Risk Estimates for Individual Real Estate Parcels for the Contiguous United States

by
Edward J. Kearns
1,*,
David Saah
2,
Carrie R. Levine
2,
Chris Lautenberger
3,
Owen M. Doherty
4,
Jeremy R. Porter
1,
Michael Amodeo
1,
Carl Rudeen
2,
Kyle D. Woodward
2,
Gary W. Johnson
2,
Kel Markert
2,
Evelyn Shu
1,
Neil Freeman
1,
Mark Bauer
1,
Kelvin Lai
1,
Ho Hsieh
1,
Bradley Wilson
1,
Beth McClenny
4,
Andrea McMahon
2 and
Farrukh Chishtie
2
1
First Street Foundation, Brooklyn, NY 11201, USA
2
Spatial Informatics Group, Pleasanton, CA 94566, USA
3
Reax Engineering, Berkeley, CA 94704, USA
4
Eagle Rock Analytics, Sacramento, CA 95820, USA
*
Author to whom correspondence should be addressed.
Fire 2022, 5(4), 117; https://doi.org/10.3390/fire5040117
Submission received: 13 June 2022 / Revised: 26 July 2022 / Accepted: 7 August 2022 / Published: 15 August 2022

Abstract

:
The methodology used by the First Street Foundation Wildfire Model (FSF-WFM) to compute estimates of the 30-year, climate-adjusted aggregate wildfire hazard for the contiguous United States at 30 m horizontal resolution is presented. The FSF-WFM integrates several existing methods from the wildfire science community and implements computationally efficient and scalable modeling techniques to allow for new high-resolution, CONUS-wide hazard generation. Burn probability, flame length, and ember spread for the years 2022 and 2052 are computed from two ten-year representative Monte Carlo simulations of wildfire behavior, utilizing augmented LANDFIRE fuel estimates updated with all the available disturbance information. FSF-WFM utilizes ELMFIRE, an open-source, Rothermel-based wildfire behavior model, and multiple US Federal Government open data sources to drive the simulations. LANDFIRE non-burnable fuel classes within the wildland–urban interface (WUI) are replaced with fuel estimates from machine-learning models, trained on data from historical fires, to allow the propagation of wildfire through the WUI in the model. Historical wildfire ignition locations and NOAA’s hourly time series of surface weather at 2.5 km resolution are used to drive ELMFIRE to produce wildfire hazards representative of the 2022 and 2052 conditions at 30 m resolution, with the future weather conditions scaled to the IPCC CMIP5 RCP4.5 model ensemble predictions. Winds and vegetation were held constant between the 2022 and 2052 simulations, and climate change’s impacts on the future fuel conditions are the main contributors to the changes observed in the 2052 results. Non-zero wildfire exposure is estimated for 71.8 million out of 140 million properties across CONUS. Climate change impacts add another 11% properties to this non-zero exposure class over the next 30 years, with much of this change observed in the forested areas east of the Mississippi River. “Major” aggregate wildfire exposure of greater than 6% over the 30-year analysis period from 2022 to 2052 is estimated for 10.2 million properties. The FSF-WFM represents a notable contribution to the ability to produce property-specific, climate-adjusted wildfire risk assessments in the US.

1. Introduction

The threat of increasing wildfire risk across the United States has been described by a number of studies that discuss both the increasing incidence of wildfire and the increasing threat to forests and communities [1,2,3]. The implications of this growing risk threaten the economic stability, natural resources, and quality of life for the affected communities and local residents, and there are a number of resources (e.g., https://wildfireresearchcenter.org/, accessed on 13 June 2022; https://wildfirerisk.org/, accessed on 13 June 2022) now available to assist communities in meeting those growing risks. Westerling et al. [2] report that land management costs already exceeded USD 1 billion in costs nearly 10 years ago; however, a report from the Bureau of Land Management (BLM) and the Western Forestry Leadership Coalition (WFLC) highlighted the fact that this direct cost is simply a fraction of the larger economic costs of wildfires [4]. The WFLC report highlighted the fact that, beyond the direct dollars spent on land management and suppression, there are additional direct costs (such as firefighting crews), indirect costs (extensive and long-term implications of lost tax revenue, land recovery, and dips in property value), rehabilitation costs (watershed restoration, short term emergency loans, etc.), and additional uncharacterized costs (including human costs). The report estimates that the cost of wildfires reported through direct costs may only account for about 3% of all costs incurred from wildfires. In fact, NOAA reports over USD 79.8 billion in costs associated with the occurrence of wildfires between the most recent 5-year period of recorded events (2018 and 2021), not accounting for much of the cost associated with land management or long-term indirect and additional costs [5]. While the costs of wildfires have been exceedingly high in recent years, it is also growing at a rate that indicates its increasing impact on communities in the US, with the cost of the preceding 5 years of economic damages totaling only USD 8.5 billion (2012–2016) [5]. This increase in damages is nearly 10-fold and represents the growing risk to communities, and residents in those communities. A number of commercial fire risk products have been developed and are in wide use in the insurance industry (e.g., Verisk’s “FireLine” https://www.verisk.com/siteassets/media/downloads/underwriting/location/location-fireline.pdf, accessed on 1 June 2022), but these are statistically based solely upon past fires and related damages.
The growing risk has been linked to a series of different drivers in the literature. Some explanations have drawn on anthropogenic changes in industry-associated latent consequences, such as forest regrowth following a decline in logging in the late 19th century, which allowed for structural changes to the biomass (fuels) in those areas driven by the lack of the natural regulation from regularly occurring fires [2]. Competing explanations focus on the impact of variability in climate conditions associated with the increasing risk of wildfires, including increasing variability in moisture conditions, increasing drought frequency, and warming temperatures [6]. Finally, these explanations are further compounded by the fact that the areas most at risk of wildfires in direct relation to residential land uses have grown extensively in recent years [7]. This interface, referred to as the wildlands–urban interface (WUI), has shown significant growth in the last 20 years, with Radeloff and colleagues [7] reporting about an 8% growth in WUI area and a nearly 35% growth in population and housing units. In total, the research reports that half of all homes built in the 1990s, and about 40% in the 2000s, were built in the WUI. Recent statistical analyses at the property level have shown that 97% of home losses are found in the WUI [8]. Such rapid growth in high-risk areas means that even more properties are at risk of wildfire. Beyond the impact on magnitude, the larger WUI populations simply mean there is more opportunity for fire as the vast majority are ignited by human cases [9].
In response to the need to respond to this growing nationwide risk at the community level, the U.S. Federal Government supported the creation and publication of the publicly available Wildfire Risk to Communities (hereafter WRC; see WildfireRisk.org) [10], which conveys the relative risk for communities based on a 270 m horizontal resolution analysis. The tool is primarily intended to provide insight for community level wildfire solutions in a way that allows for communities to understand their relative risk comparatively with other areas, so that resources can be allocated in a measured and efficient way, with the goal of combating economic and human loss from wildfires. Wildfire Risk to Communities’ estimates are based on fire simulations that incorporate the US Forest Service’s 2014 Landscape Fire and Resource Management Planning Tools database v2.0.0 [11], with some modifications (Smail, personal comm. 2021), which provides open data describing the composition and state of fuels across the contiguous United States (CONUS). However, WRC’s focus is on community risk and actions to reduce those risks, and the metrics computed are not focused on individual properties and homes, nor does WRC include the impacts of climate change on future risk.
The development of the WRC tool served as a milestone in giving communities the ability to assess risk in their area and plan for resource allocation in relation to that risk. However, the developers of the tool acknowledge that it is a community level tool and should be used for community level purposes. This research aims to build upon the wildfire community’s considerable research on wildfire risk modeling [12], and to complement the WRC community level tool with a high-resolution model developed specifically for the property level at a national scale, the First Street Foundation-Wildfire Model (FSF-WFM). Given the increase in wildfire occurrence and the subsequent economic consequences [13], there remains a need to quantify the probable changes in wildfire exposure for US property owners and residents to provide to them with an improved awareness of their specific, property-level wildfire risk now and their expected risk in the future. The use of wildfire hazard estimates to provide property-level vulnerability estimates has been demonstrated in numerous studies, e.g., [14,15]. As the number of communities in the built environment suffering extensive losses grows (e.g., losses in the WUI exemplified by Gatlinburg, TN 2016; Paradise, CA 2018; Grand County, CO 2020; Boulder County, CO 2021), there is also a recognized need to describe the spread and risk of wildfire specifically within the WUI [16]. The development of such a model is based on the unique risk each individual property faces, based on property-level characteristics, and can be scaled nation-wide to provide homeowners with mitigation solutions, such as those included in the “resilience pathways” described in [17].
Building upon the WRC approach, the LANDFIRE database, climate projections, and existing open-source fire behavior models, the remainder of this document is designed to provide a transparent understanding of the framework and methodology that went into the development of the property-level wildfire model, taking an open science approach (https://earthdata.nasa.gov/esds/open-science, accessed on 13 June 2022). This study does not attempt to provide quantitative comparisons between the outputs of the FSF-WFM and the WRC approaches. While comparisons may be useful in understanding the nuances of the fuels used and model implementations, any direct quantitative estimates of the differences are difficult to interpret, not just due to those differences, but also because the models were developed with different purposes in mind. Direct quantitative comparisons with the aforementioned property-level statistical models typically used in insurance applications may be useful, since they are more similar in purpose, but due to the proprietary nature of and costs associated with those models’ outputs, the authors do not currently have access to those outputs at a sufficiently large scale to conduct such a comparison. Any such comparisons of results may be the subject of a future study, but would specifically be a comparison of methodological differences of scale and purpose versus a comparison of accuracy of the models. To that point, the model described in this paper is specifically designed to measure property risk and should be thought of as complementary to the larger community risk products.

2. Model Development

The FSF-WFM approach is based on the application of a fire behavior model to explore the incidence, severity, and probability of wildfires that occur at a property-level resolution across CONUS. This general approach has been shown to be useful at large scales in the aforementioned WRC using FSim [18], and on regional scales, such as the use of WyoFire [19]. Here, we use an open-source wildfire behavior model, ELMFIRE (Eulerian Level Set Model of Fire Spread), which has likewise been shown to produce useful results in this type of application [20], but also extends its use to estimate future wildfire hazards based on climate predictions.
The development of the FSF-WFM includes a series of steps associated with the integration of fuels, fire weather, and ignition locations into ELMFIRE. While each of these components will be explained in detail below, a definition/purpose of each component as they relate to the wildfire model is provided here for context.
  • Fuels: estimation of the fuels that support wildfires across the US at 30 m horizontal resolution, including assembly of new fuel estimates updated with disturbance descriptions for the previous 10 years and the conversion of buildings within the WUI into a burnable fuel type that allows the appropriate progression of wildfire throughout the WUI in the fire behavior model.
  • Fire weather: assembly of the weather data to drive the fire behavior model under a representative range of fire weather conditions for 2022 and 2052. Fire weather was derived from the National Oceanic and Atmospheric Administration’s (NOAA’s) surface weather reanalysis for 2011–2020 to create the 2022 hazard layers, and was driven by the same time series in 2052 with air temperature, precipitation, and humidity scaled to 2052 conditions, as represented by downscaled International Panel on Climate Change (IPCC) climate model ensemble results.
  • Ignition locations: identification of the likely ignition locations, temporal fire occurrence patterns, and conditions most likely for fire spread for future wildfires.
  • Fire behavior model: application of a fire incidence and landscape behavior model across the contiguous United States in a Monte Carlo simulation to build probabilistic estimates of 2022 and 2052 wildfire hazards in terms of burn likelihood, fire intensity, and spread of embers at 30 m horizontal resolution.
The resulting wildfire hazards product is based on the data sources listed, which were used to update the data to May 2021 (see Appendix A).

2.1. Fuels

The wildfire hazard estimate is heavily dependent upon estimates of the type, quantity, age, and condition of the combustible fuels across the US. Version 2.0.0 of the canonical U.S. Forest Service (USFS) LANDFIRE [11] fuels dataset at 30 m horizontal resolution is utilized as a baseline for provision of this fuel information, and is updated to characterize the risks in the present through the inclusion of all known disturbances from May 2021 to create a current fuels layer that is useful for assessing wildfire risk for the year 2022. One must note that not all disturbances were able to be adequately documented or described, and different US states exhibit different levels and styles of reporting. States with the highest fire risk in the Western and Southeastern US (e.g., California, Oregon, Arizona, Colorado, Washington, Idaho, and New Mexico) were prioritized to ensure their adequate inclusion in this study. These disturbances were incorporated as changes to surface and canopy fuels by modifying the geographically referenced LANDFIRE classifications, and include recent wildfires, prescribed burns, harvests, and other forest management practices, as reported by the data sources listed in Appendix B. Modification of the fuel descriptions was carried out in accordance with the LANDFIRE fuel classes and methodologies, and is congruent with the LANDFIRE disturbance code schema, which consists of thematic three-digit code values corresponding to disturbance type, severity, and time since disturbance, respectively, per the LANDFIRE Fuel Disturbance Attribute Data Dictionary [21]. A representation of the processes is shown in Figure 1 that describes the methods used to create the fuels estimate for this study.

2.1.1. Disturbances

Disturbances from wildfires across CONUS were incorporated by using data shared by the Monitoring Trends in Burn Severity (MTBS) [22] program, which maps the burn severity and extent of large fires across all lands in the US. At the time of analysis, the MTBS dataset included fires of an area larger than 500 acres through 2019. Therefore, for the year 2020, the MTBS dataset was augmented with data from all fires of size <500 acres from the National Interagency Fire Center (NIFC).
To ensure the consistency of fire severity characterizations between the MTBS and NIFC datasets, burn severity was informed by calculating the normalized burn ratio (NBR) [23] for one pre-fire and one ninety-day-window post-fire cloud-filtered composite image corresponding to each fire. The pre-fire NBR was then subtracted from the post-fire NBR to create the relative difference normalized burn ratio (RdNBR) index [23]. “Miller’s threshold” [23] was then applied to the RdNBR image to create a five-class burn severity classification.
For non-wildfire disturbances, including harvest, fuel mitigation treatments, and prescribed burns, there are no uniform naming or reporting conventions for forest management practices across the U.S. and the quality of data entry varies considerably from state to state. To ensure that every feature is assigned a standardized disturbance class, all unique treatment names from every dataset were compiled for review by forestry field experts who are included in the authorship of this paper. Each unique disturbance name in the document was assigned a LANDFIRE disturbance type, and assigned the appropriate three digit LANDFIRE disturbance code that captures disturbance types, severity, and time since disturbance. A distribution associated with the types and severity of disturbances is reported in Table 1.
The disturbance types most frequently found in our dataset and listed in Appendix B were fire (disturbance type 1), mechanical add (disturbance type 2, when fuels are mechanically mowed or chipped and transitioned to surface fuels), and mechanical remove (disturbance type 3, when fuels are removed via cutting, felling, burning, or harvest). We assigned a disturbance type of “other” (disturbance type 8) to chemical treatments and grazing. We excluded treatments or activities included in the datasets that would not have impacted fuels (including but not limited to seeding, habitat restoration, and invasive species removal). Treatment disturbances, such as hand thinning, piling, prescribed fire, and other treatments where canopy cover is not altered were assigned a disturbance value of 1 (low severity); mechanical thinning and harvest were assigned a disturbance value of 2 (medium severity); and clear cuts were assigned a disturbance value of 3 (high severity). For wildfire disturbances, we followed the MTBS conventions, whereby fire severity class 2 are low severity, 3 are medium severity, and 4 are high severity classifications. Classes 1 (unburned/unchanged) and 5 (increased greenness) were considered undisturbed. The code for time-since-disturbance was determined based on the year of treatment and the LANDFIRE zone. Time-since-disturbance was categorized as 1 (disturbances that occurred in 2020), 2 (disturbances that occurred in 2015–2019), and 3 (disturbances that occurred in 2011–2014). Due to differences in overall fire risk topographies, for disturbances that occurred in the LANDFIRE Southeast Super Zone (Zones 46, 55, 56, 58, and 99), the time-since-disturbance categories are 1 (disturbances that occurred in 2020), 2 (disturbances that occurred in 2017–2019), and 3 (disturbances that occurred in 2011–2016). Finally, the treatment and wildfire layers were combined into a single disturbance layer using a priority ranking ruleset informed by LANDFIRE analysts (Smail, personal comm. 2021) to ensure the most fuels-relevant disturbance value is assigned in cases of spatial overlap.
Validation with the CONUS scale is most practically accomplished with remote sensing techniques. The Hansen Global Forest Change dataset [24] provides a ‘loss year’ band that represents the year(s) when there was detectable canopy loss during the period 2000–2020 at the 30 m per pixel scale. We leveraged this band to create a forest loss bitmask for 2011–2020 and applied it to screen our final aggregate disturbance layer to remove false positives of moderate and high severity harvest [24].

2.1.2. Fuel Layers

Using LANDFIRE v2.0.0 as the base, four canopy fuel layers (canopy cover, canopy height, canopy base height, canopy bulk density) and one surface fuel layer (40 Scott and Burgan Fire Behavior Fuel Model, hereafter FM40) [25] were generated with an effective year of 2021 for use as inputs into the fire models. Fuels were only transitioned in areas that were disturbed between 2011 and 2020. Initial layers that represented lookup rulesets in the LANDFIRE Total Fuel Change Tool (LFTFCT) database were generated. First, canopy cover and height midpoint layers are derived from the LANDFIRE Fuel Vegetation Cover (FVC) and Fuel Vegetation Height (FVH) rasters based on the LFTFCT lookup table values. Next, using the new updated disturbance layer, a canopy guide layer was generated by using the LFTFCT master lookup table applied to unique combinations of the disturbance code, biophysical Settings (BPS), fuel vegetation cover (FVC), fuel vegetation height (FVH), and fuel vegetation type (FVT). The four canopy fuel layers are then generated using the following regression equation:
C a n o p y   F u e l = C X + H y + b
where C is the canopy cover midpoint, H is the canopy height midpoint, x and y are the scale factors, and b is an intercept value derived from a lookup of unique disturbance code and FVT combinations from the LFTFCT lookup table. For canopy cover and canopy height regressions, the cover and height midpoint values are derived from the initial FVC and FVH midpoint layers described above, while for canopy base height and canopy bulk density, the midpoint values are derived from the new canopy cover and height layers that were generated in the step described above. Additionally, canopy bulk density uses a ruleset to create two stand height coefficients from the canopy height midpoint value for pixels following the rules described in [26]. Each canopy fuel regression output is post-processed to ensure values are within the LFTFCT’s valid value range (CC: 0–95; CH: 0–510; CBH: 0–100; CBD: 0–45), scaled properly, and binned, if necessary, to defined midpoint values [21]. Finally, the LFTFCT canopy guide layer is applied to each layer using rulesets based on canopy cover thresholds [21].
The FM40 surface fuel estimates are generated in the same way as the canopy guide, using the LFTFCT master lookup table applied to unique combinations of the disturbance code, BPS, FVC, FVH, and FVC. Products generated include the necessary LANDFIRE fuel and vegetation datasets for the workflow described here, derived fire severity, canopy cover and canopy height midpoint, as well as disturbance estimates. Included with the 2021 fuel profile used in this study are the following five updated 2021 fuel layers: FM40, canopy cover (CC), canopy height (CH), canopy base height (CBH), and canopy bulk density (CBD). Figure 2 highlights the spatial location of the canopy and surface fuel updates across the CONUS, with Figure 3 highlighting the update of surface fuels in a more local context.

2.1.3. WUI Surface Fuel Updates

Typically, homes and other buildings in the built environment, including the WUI, are classified as non-burnable fuels within LANDFIRE. However, in order to allow the estimate of wildfire hazards within the WUI under the full range of fire weather conditions, those properties within the WUI need to be replaced by a burnable fuel estimate to permit the wildfire behavior model to estimate how wildfire could move through the WUI more accurately.
The first step of developing the WUI fuel model was to derive a current map of WUI areas. WUI areas are defined by the following two factors: building density and the distance from wildland vegetation [27]. We used the 2016 NLCD existing vegetation cover layer to identify areas of wildland vegetation, and derived our own building-density layer from MapBox building footprints (Appendix A), following evidence from Caggiani et al. [8] that such higher-resolution analyses enable more precise evaluations of wildfire risks. The WUI influence zone, WUI intermix, and WUI interface layers were defined as the following [16]:
  • Influence zone is >75% land coverage of wildland vegetation within 1 mile of a residence.
  • Intermix is >1 residence per 40 acres and groups of residences larger than 50 acres, with >50% land coverage of wildland vegetation.
  • Interface is defined as >1 residence per 40 acres and groups of residences larger than 50 acres, with <50% land coverage of wildland vegetation, and within 1 mile of wildland vegetation.
Non-burnable pixels were converted to a burnable FM40 fuel type in the WUI intermix and interface only, as much of the WUI influence zone is already estimated as burnable in LANDFIRE and does not need to have non-burnable cells converted to burnable cells to enable the fire behavior model in those areas. Any unnecessary conversions within the influence zone could potentially result in biased fire behavior by changing the FM40 fuel types in those areas.
Properties within the WUI with a non-burnable classification in the 2021 fuel profile were replaced by an effective fuel type by estimating it from a statistical analysis of 549 historical fire perimeters in the WUI from 2014–2019 (see Appendix C). These past fires were used to train a random forest machine learning algorithm to predict the appropriate fuel classification. One must note that the fuel layers do not take into account fuel estimates for the structures themselves in the WUI that could lead to increased house-to-house ignition probability; such an approach could be incorporated into a future effort. To convert non-burnable pixels in the WUI intermix and interface to allow the fire behavior model in those regions, we used a machine learning approach, as described below.
The 2021 FM40 fuel types derived in the fuel workflow (see Figure 1) described above are used as the response variable. The training and testing datasets were composed of pixels in the WUI intermix and interface that were within the fire perimeters from our disturbance dataset (2011–2020) or within a 1 km buffer around the fire perimeter, in order to capture areas that remained unburned in those incidents. Other variables included vegetation products from LANDFIRE v2.0.0 [11], Landsat data derived from 4-month composites encompassing each training fire’s ignition date (coastal, blue, green, red, NIR, NDVI, SWIR1, SWIR2, NDVI, MNDWI, BAI), GRIDMET data derived from 1-month composites encompassing each training fire’s fire ignition date (tmin, tmax, fm1000, vs, mndwi, erc, bi), topography variables from USGS (slope, elevation, aspect), building density per 1 km2, the number of structures destroyed per fire, and fire severity. A random-forest model was trained only on burnable FM40 categories. The prediction area was limited to 2021 FM40 urban/developed (FM40 class 91) and agricultural (FM40 class 93) land in the WUI intermix and interface. We ran a stratified k-fold cross validation training using 10 folds for each dataset, with the training data split 80−20% in each run. The best model had a k-fold training accuracy of 73.0%, with a mean precision (true negative rate) of 76.4% and a mean recall (true positive rate) of 73.0%. Overall, k-fold training had a mean model accuracy of 71.2% (68.7–73.7% confidence interval). The overall training accuracy was 96.9%, with a training Kappa coefficient of 96.8%. For model testing, 30% of the sampled data was withheld. The independent validation dataset showed 71.7% accuracy, with a testing Kappa coefficient of 73.0%. A framework for documenting the classification process to replace non-burnable FM40 classes with burnable classes in the newly defined WUI is presented in Figure 4, with feature importance highlighted in Figure 5.
Overall, the vegetation, topographic, and weather variables had higher importance than the fire-related or building density-related variables in the model (Figure 5). In general, non-burnable WUI intermix and interface pixels were frequently replaced with grass or grass-shrub fuel types (FM40 classes 101 and 121, and occasionally 103 in the southeast and 183 in the mid-Atlantic). The predicted pixels were replaced in the 2021 FM40 fuel layer to create the final 2021 FM40, with the surface fuel model updated for both disturbances and WUI areas (see Figure 6, for example).

2.1.4. Vegetation Changes and Impacts on Fuels

Changes in the composition and volume of vegetation due to climate change’s impacts have been discussed in depth by a number of researchers, including Westerling et al. [2], Radeloff et al. [7], Krawchuk et al. [28], and their importance to estimates of fire intensity has been discussed more recently in a review article by Bowman et al. [29]. These studies typically examine those vegetative changes over time periods of 75 to 150 years, while the current study is focused on 30 years only. To investigate the size and scope of vegetation changes on a 30-year time period at 30 m horizontal resolution, we originally planned to utilize the Land Use and Carbon Scenario Simulator [30], a Monte-Carlo based state-and-transition simulation model, to project changes to 14 carbon pools from 2021–2051. While we observed statistically significant changes in above ground modeled carbon pool volumes over 30 years across CONUS, we struggled to accurately translate this from those carbon pools to the canopy and surface vegetation classes needed to drive the fire behavior model we employ in this study. While research continues on this and several alternate ways of estimating the vegetation and fuel changes anticipated across CONUS over 30 years in a changing climate, we have elected to hold the fuel constant between the 2022 and 2052 simulations for the purposes of this study. Future wildfire exposure estimated by the model described in this study will then be independent of future vegetation changes, and will depend only on the future weather impacts on the fuel conditions and fire behavior alone.

3. Fire Weather and Climate Change

The primary inputs needed to drive the fire spread model are fuels, topography, and weather. This section details the integration of climate weather into the development of the larger FSF-WFM. The weather that can drive the growth and distribution of wildfire can be separated into the following two categories: (1) the weather before the onset of a wildfire that impacts fuel condition by making the fuels drier or wetter, and (2) the ‘fire weather’ that occurs at ignition, which can increase intensity and drive fire across the landscape. To represent a wide range of possible weather-driven fire conditions across the landscape within the simulations employed here, we used a decade of high spatial and hourly resolution weather data. Wind speed and direction, relative humidity, and temperature inputs were assembled from the Real Time Mesoscale Analysis (RTMA) dataset [5], which provides hourly estimates of sensible weather variables on a 2.5 km grid for CONUS. The RTMA surface weather data reanalysis from 2011–2020 was augmented by Oregon State PRISM (Parameter-elevation Regressions on Independent Slopes Model) [31] precipitation data to fill in gaps in the RTMA data. Ten years was chosen to represent a wide range of weather conditions, while overlapping the time period for which the fuel state is represented (i.e., LANDSFIRE 2016 augmented to 2020). While a 20- or 30-year time series would provide a more complete sampling of the possible meteorological conditions, the 10-year time series does include multiple La Niña and El Niño phases and allows for the computations to be completed in a reasonable span of time, given the available resources. Additionally, since this study does not set out to replicate or predict anomalously large or intense fires (e.g., plume fires) in a deterministic sense, the Monte Carlo approach used will deemphasize those extreme or infrequent conditions and instead emphasizes the much more frequent medium-large fires (i.e., larger than those that are easily suppressed, but smaller than the rare extreme fires).
To represent the 2052 weather, we have considered the 2048–2057 time series, created by scaling the hourly 2022 RTMA time series, to forecast 2052 conditions. To do this, we used the International Panel on Climate Change’s (IPCC) Fifth Coupled Model Intercomparison Project (CMIP5) ensemble results [32] following the Representative Concentration Pathway 4.5 (RCP 4.5), as downscaled within the daily Multivariate Adaptive Constructed Analogs (MACA) v2 product [33] to represent the expected weather conditions in 2052 across CONUS. The RCP 4.5 climate model results were chosen to be relatively conservative in outlook, and to be consistent with previous and similar work conducted for future flood risk authors [34,35].
Surface winds were held constant from the 2022 to the 2052 simulation period to preserve the realistic and high-resolution aspects of the NOAA RMTA time series in the future, to reduce uncertainties in future fire behavior and in recognition that future winds are likely to change far less significantly with climate change than other weather parameters [32,33]. The ELMFIRE fire behavior model is necessarily very sensitive to winds, and downscaled climate model results have difficulty resolving the local and orographic effects in the wind fields to a sufficient fidelity to support such fire models [36]. Even if they captured the spatial variability adequately, the high-resolution winds generated by an atmospheric model driven by boundary conditions generated from the climate model outputs would still require extensive verification and validation to be able to use them for our simulations and justify the results. Since the goal is not to recreate any particular fire event, but to use the weather time series to support a range of conditions suitable for Monte Carlo simulation, we concluded that holding the winds constant from the 2011–2020 time series to drive 2052 fire behavior would be a reasonable approach.
With winds held constant, the other 2022 weather variables underwent scaling to create a 2048–2057 hourly times series used to derive the 2052 wildfire hazards. The MACAv2 downscaled CMIP5 RCP4.5 outputs at daily resolution were used to scale the RTMA hourly time series of air temperature, relative humidity, and precipitation by computing bias adjustments between the present-day 2022 and forecast 2052 conditions (Appendix D). The biases were distributed throughout the day via gamma distribution to maintain the diurnal signal in precipitation and humidity, while allowing for the overall scaling to be representative of the climate change impacts on these variables. Extreme values in biases were adjusted inward (towards the center of the distributions) to allow for consistent statistics, while preserving the general climate variability. Air temperature adjustments at the hourly resolution were likewise adjusted with a simpler gaussian distribution that brought daily average values of the 2011–2020 RTMA hourly time series in line with the future 2048–2057 MACAv2 daily values.
The result for the 2052 weather time series is a 10-year duration, hourly resolution representation of the estimated future weather conditions at 2.5 km horizontal resolution that are characterized predominately by 1.7–2.8 deg C (3–5 deg F) average warmer temperatures across CONUS. This allows the impact of higher air temperatures from climate change on fuel conditions in 2052 to be largely isolated and evaluated, since winds and fuels are both held constant from 2022. The greatest deficiency of this approach is that it is not possible to evaluate the climate impacts of geographically coherent but temporally variable features, such as more severe or longer droughts, or greater incidences or intensities of atmospheric rivers or hurricanes. As such, these estimates are limited almost entirely to the effects caused by higher air temperatures on fuel conditions, and so must be considered an underestimate of the total possible effects of climate change on wildfire probability. Subsequent versions of this model are intended to address these deficiencies.

4. Ignition and Spatial Fire Occurrence Patterns

One of the primary indicators of where future fires will occur is informed through historical fire occurrence data. The spatial component of the fire occurrence model is built from the Fire Occurrence Database (FOD; https://www.fs.usda.gov/rds/archive/Catalog/RDS-2013-0009.5, accessed on 1 June 2022) developed by the USDA Forest Service [37,38]. The FOD includes 27 years (1992–2018) of fire occurrence data, encompassing 2.17 million georeferenced wildfire records that total 165 million acres burned. Following the best practices for annualized burn probability modeling [18], this database was filtered to remove small fires, defined as those that are less than 100 acres (Class A, B, and C fires). We acknowledge the choice of the 100 acre cutoff is somewhat arbitrary, and different thresholds (e.g., 300 acres [18], 247 acres [39]) have been used in other research and models, but was chosen as a convenient approximation of the typical scale of wildfires whose growth are often limited by human fire suppression activities.
A recognized best practice is to develop an ignition density grid using a kernel density tool [18]. The ignition density kernel formula used (see equation below from [40]) was implemented in the wildfire behavior model to generate the ignition density grid for this work, where r is the search radius (bandwidth) and di is the distance from point i to the centroid of a given cell.
D e n s i t y = 1 r 2   Σ i = 1 n 3 π   x   ( 1 ( d i r ) 2 )   f o r   d i < r

Modeling Temporal Fire Occurrence Patterns

The previous section describes how the spatial fire occurrence is modeled, but it does not address when large fires may occur. One of the strongest predictors of temporal occurrence of both the number of large fires and acres burned is the National Fire Danger Rating System (NFDRS) Energy Release Component (ERC) percentile based on fuel model G, or ERC(G)’ [41] (note ERC(G) refers to raw ERC values (Btu/ft2) and ERC(G)’ refers to ERC percentiles). ERC is 4% of the energy per unit area (Btu/ft2) that would be released during a fire. ERC depends on live and dead fuel loading by size class (as characterized by an NFDRS fuel model), as well as fuel moisture content of live and dead fuels. Although NFDRS fuel model G, which shows the best correlation with fire occurrence and burned area, contains loadings across all dead fuel size classes and live herbaceous/live woody loadings, it has a heavy loading in the 1000-hr size class. For that reason, ERC(G) is primarily a function of weather conditions over the preceding 45 days and can be thought of as a measure of intermediate to long-term dryness and as it is calculated solely from fuel moisture content, ERC is not a function of wind speed, slope, or spread rate.
Fire occurrence is normally assessed in terms of ERC percentile, as opposed to raw ERC (Btu/ft2), because ERC percentile shows better correlation with fire occurrence and size than raw ERC, since the same amount of precipitation that corresponds to wet conditions in one region may correspond to dry conditions in another region.
l o g 10   n = 0.02768   ×   E R C   ( G ) 0.2333
Figure 7 shows the number of large fires in the Western US as a function of ERC(G)’. The data in Figure 7 are demonstrably well-fit (R2 = 0.94) by the correlation in the equation above [41], which is used in the wildfire behavior model to calculate fire occurrence from ERC(G)’.

5. Wildfire Behavior Model

In the development of the FSF-WFM, we employed the open-source wildfire behavior model, ELMFIRE, which is a highly parallelized model that was used to both simulate fire spread and quantify the wildland fire hazard via Monte Carlo simulations. ELMFIRE is a Rothermal-based, level set model used to track boundaries across the landscape based on the numerical solutions of [42] and is fully described in Lautenberger [43].
The overall fire hazard and probability modeling methodology, as shown graphically in Figure 8 and described in this section, is based on the work of Finney et al. [44], best practices described by Scott et al. [18], and a relatively recent review of simulation-based burn probability modeling [45]. Consequently, the contribution of this work is not developing new techniques or approaches to fire probability and hazard modeling, but rather implementing computationally efficient and scalable modeling techniques based on existing fire probability and hazard modeling paradigms pioneered by the aforementioned authors. These scalable computing techniques make it possible to conduct CONUS scale fire probability and hazard simulations at 30 m resolution in a reasonable amount of time, using commodity-style computational resources. The CONUS domain was subdivided into 48 km by 48 km tiles, which were likewise surrounded by 8 similar tiles in a 3 × 3 grid pattern, to aid in the distributed compute workflow.
Inputs to ELMFIRE include fuels, weather time series, and ignition locations. The ignition locations were based on historical (1992–2018) fire locations described in the previous section, and limited to fire sizes of greater than 100 acres. This limitation allows the implicit inclusion of the effect of human-driven fire suppression activities in the model output to create a “real world” estimate of fire exposure—i.e., wildfires that are actively prevented from growing large. For example, the State of Rhode Island has exhibited remarkable fire suppression over the past decades and has been able to eliminate all fires over 100 acres during the 1992–2018 time period, driving the effective burn probability in Rhode Island to zero for all properties in our simulations.
For each ignition location, a weather “draw” was randomly selected for that fire that would be carried forward many hours in simulation, and could extend anywhere in the 3 × 3 (144 × 144 km) tile domain. Those simulated fires that grew to sufficient size (100 acres) were tracked and the locations, fire length, and durations were noted. This process was repeated over 100 million times, and resulted in approximately 8–10 million tracked fires of significance per simulation (2022 and 2052). The result is a statistically well-characterized set of simulated wildfires, from which the probabilistic exposure of properties and buildings to wildfire hazard based on likelihood (i.e., burn probability), flame length (i.e., intensity), and ember cast may be derived. The likelihood of a 30 m pixel burning is the number of times that the pixel had ignited over the course of all the simulations. The flame length is a measure of fire intensity, captured as binned flame lengths (see Table 2) over the distribution of all fires within the pixel, and may be expressed as the mean, median, or maximum flame length. The ember cast is a binned measure of the number of times embers, pushed ahead of a simulated fire by the fire weather time series, land in a pixel and results in an ignition of the fuels in that pixel.

5.1. Fire Spread Model

The 2D fire simulator ELMFIRE is used here to drive a stochastic fire spread analysis that is used to generate the CONUS burn probability and hazard estimates. ELMFIRE’s computational engine is similar to other two-dimensional fire simulators, such as FARSITE [46], in that it calculates surface fire spread rate using the Rothermel surface spread model [47,48], assumes that each point along the fire front behaves as an independent elliptical wavelet [49], with length to breadth ratio determined empirically [48,50], simulates transition from surface to crown fire using the Van Wagner criterion [51] (with crown fire spread rates calculated from Cruz et al. [52]), and models ember-driven ignition or “spotting” as a stochastic process with lognormal spotting distance distribution [53,54]. ELMFIRE tracks the fire front using a narrow band level set method [55], a numerical technique for tracking curved surfaces on a regular grid.
To demonstrate how ELMFIRE simulates fire spread, Figure 9 shows 24-h of fire progression from an individual ignition site. The black contour lines in Figure 9a represent the fire front position at 2-h intervals. Figure 9a also shows which parts of the burned area experienced surface fire (blue), passive crown fire (green), or active crown fire (red). Figure 9b similarly shows fire perimeter contours and flame length variation within the fire perimeter. Flame length is highest in the areas that burn as heading fires or that experience crown fire and lowest in the areas that burn as a flanking, backing, or surface fire. In this example, the fire area after 24 h of spread is approximately 560 acres.
The Monte Carlo fire spread analysis conducted here involves running millions of fire spread simulations (similar to that shown in Figure 9) sequentially over many years (2011–2021, and 2048–2057), and across all tiles in the CONUS domain. Each tile is a 144 km by 144 km tile within CONUS, consisting of a 48 km central tile surrounded by its eight neighboring tiles of the same size. For each year and tile, fuel, topography, and yearly weather, fuel moisture, and ERC percentile inputs are assembled. Starting at the beginning of the simulation year, ignition locations are determined using the spatial and temporal fire occurrence modeling techniques described earlier. Fires are ignited only in the central 48 km tile, but are allowed to spread into the adjacent eight tiles within the simulation. The progression of each fire is modeled for a randomized spread duration up to 7 days from the time of ignition, to roughly approximate the varying duration of the observed wildfires. For each pixel within the modeled fire perimeter, the burn incidence is recorded, and the binned distributions of discrete ember count and flame length are also recorded for each pixel. This ignition-burn-record process is repeated for each day in each simulation year, building up the probabilistic estimates of burn probability, flame length, and ember spread. Since fires can start in one tile and spread to adjacent tiles, each tile is post-processed concurrently with its eight neighbors.
The primary outputs after processing are conventional annualized wildfire hazard maps at 30 m resolution within CONUS, composed of the following elements:
  • Burn probability—an estimate of the likelihood that a region on the landscape burns in any single year during the simulation period.
  • Fire intensity—the distribution of conditional (i.e., upon burning) flame lengths for each pixel, within discrete flame length bins.
  • Exposure to embers—similar to fire intensity, a distribution of ember exposure per pixel to characterize the relative intensity of ember exposure from all modeled fires.

5.2. Validation

To validate the results from the fire behavior model, we compared the model fires against historical fires’ intensity and size in aggregate. Example results from a tile-by-tile comparison of modeled and historical fires are generated as each geographic tile is run, as shown below. The modeled fire sizes are larger than in the FOD because (a) there is no fire suppression element applied within ELMFIRE, and the (b) simulation end time was randomized. To partially compensate for these limitations, as stated previously, the ignition layer was limited to sources of historical fires that were a minimum of 100 acres. This assumes that suppression measures would be effective in keeping such fires small, and of short duration. The resulting comparison of the modeled fires’ sizes and intensities (Figure 10) shows that the modeled fires without explicit suppression and with randomized durations up to 7 days are systematically larger than the observed wildfires. The area of non-zero burn probabilities in the resulting hazard layers should, therefore, be considered an overestimate of the likely range of wildfire spread, which creates distributions that err on the side of caution when understanding wildfire exposure (i.e., there are likely fewer false negatives). The introduction of active fire suppression within the model is the subject of further research and may be incorporated into future versions.

6. Results

The construction of a national-scale, property-specific wildfire hazard model using an open-source fire behavior model, driven by openly available inputs, has been proven possible by our development of the FSF-WFM. The ability to extend the wildfire hazard into the WUI by replacing nonburnable LANDFIRE fuel designations with estimates derived from historical fire behavior in WUI areas was also shown to be feasible. Using the model in a Monte Carlo simulation, driven by historical ignition locations across CONUS to provide 30 m-resolution hazards, it was shown to be practical using commodity-scale computing hardware. This same scheme was shown to be applicable to both current (2022) and future (2052) scenarios, given the future estimates of climate-adjusted weather conditions.
The results of the FSF-WFM model implementation are freely and publicly available through riskfactor.com, (accessed on 8 August 2022) and show property-by-property assessments of exposure to wildfire hazard. Figure 11 shows a representative parcel from the over 143 million available, and shows the levels of resolution and discrimination among properties that are available. These results are summarized at the state level in Table 3 and Table 4, and Figure 12A,B, which will be discussed in more detail below.
The spatial variability in the distributions of the hazard at 30 m resolution, including within the WUI and the prevalence of hazard in the Eastern as well the Western U.S., highlight the importance of understanding wildfire risk at a property level across CONUS. While this paper focuses on the methodology and defers a thorough analysis of results to a later study, we present some general results to provide the reader with a sense of feasibility of the FSF-WFM to address current and future wildfire exposure. Overall, the results estimate that 71.8 million properties have a burn probability of >0 in the current environment (2022) and that probability increases by 11% over the next 30 years, and grows to 79.8 million properties in CONUS in 2052. Many of those properties have low, but not zero, burn probabilities from the model so we choose to describe two general levels of wildfire hazard based on a cumulative burn probability of 3% over the 30-year period, which we label “any exposure”, and a cumulative likelihood of 10% over the 30-year period, which we label “major exposure”. When looking at those two categories, we find about 20.2 million properties in the CONUS being subject to “any exposure” and 5.9 million properties being at “major exposure” to wildfire over the 30-year period (2022–2052). These property counts represent about 15% and 5% of all property parcels in the CONUS, which further highlights the large exposure of properties in the US to wildfire exposure. For further context, flooding, which is generally referred to as the most widespread climate peril in the US, impacts about 21.8 million properties at the “any flood” level (equivalent to 6% 30-year aggregate) and about 14.6 million properties at the “significant flood” level (equivalent to 26% 30-year aggregate) [56].

7. Any Exposure

Table 3 and Figure 12A report the results of the model when applied against individual property structures and parcel centroids (on parcels without buildings). The results indicate that the top five states in regards to “any exposure” are Texas, Florida, California, North Carolina, and Alabama. In those 5 states alone, there are nearly 30 million properties with at least a 1% cumulative probability over the next 30 years of being impacted by a wildfire. Figure 12A (upper) further illustrates that the distribution of properties “any exposure” of wildfire are disproportionately located in Texas, California, and the Southeastern US. When taking into account “any exposure” of wildfire relative to the total housing stock in Figure 12A (lower), the Mountain West states of Montana, Idaho, Wyoming, and Utah emerge as a cluster of disproportionate potential impact, along with New Mexico, Oklahoma, Mississippi, and Alabama across the southern tier of the country. The Midwest and Northeast are relatively lower in regards to “any exposure” to wildfire over the next 30 years, which is expected given the climate conditions that generally drive the peril.

7.1. Major Exposure

Table 4 and Figure 12B report the results for only those properties at “major exposure” to wildfire (3% cumulative likelihood over the 30-year period). When only looking at this subset of properties, California stands out as having the most exposure, with over 2.5 million properties in this category. Texas, Florida and Arizona, at 1.7, 1.5, and nearly 1 million properties at “major exposure”, respectively, together with California, account for over 6.5 million properties that meet the threshold of having at least 3% cumulative wildfire exposure over the next 30 years. Figure 12B (upper) highlights the fact that when shifting from “any exposure” to “major exposure”, the majority of that exposure is held in the Western US, with Florida, Mississippi, New Jersey, and North Carolina standing out as states in the eastern half of the country with higher levels of exposure than the surrounding areas. Figure 12B (lower) shifts that impact slightly when accounting for the exposure as a proportion of properties in the state. Using that metric, Arizona, Utah, and Wyoming carry the most exposure to wildfire hazard, followed by their western neighbors, California and Nevada.
The estimated geographic distribution of change in wildfire exposure due to climate change is shown in Figure 13. The percentage increase between the current year and 30 years into the future in the average burn probabilities of properties with at least 0.03% risk is at least 100% in many of the counties across the country. The annual burn probability of 0.03% corresponds to at least a 1% cumulative likelihood over a 30-year period. With higher burn probabilities, a higher incidence of losses is expected over time, as properties are exposed more often to wildfires.
Finally, a fire factor risk assessment was created on a property-level basis across CONUS. Property parcel geometries are provided by the Lightbox public-record property boundaries database. Building footprint geometries are defined by Mapbox. First Street performed a geometric intersection to match parcels to building footprints. Footprints that cross parcel boundaries were subdivided, such that no footprint geometry crosses parcel boundaries. Since some parcels intersect multiple footprint geometries, the building footprint with the largest area was designated the primary footprint.
To evaluate the exposure to wildfire flames and embers, each hazard layer was queried at the geometric centroid of each building footprint and parcel. For scoring purposes, at properties with a building footprint, the statistic at the primary footprint centroid was recorded; for parcels without a building footprint, the parcel centroid was recorded. The assignment of a 30-year, climate-adjusted aggregated wildfire risk score was then computed by calculating the likelihood and nature of exposure through burn probabilities in and belongingness to an ember zone for a building or parcel as representative of the risk for each property for 2022 and then for 2052, and then linearly interpolating this across that 30-year period.
The annual risk, as defined by burn probability in and belongingness to an ember zone for each year, was summed across the 30-year period and was used to derive the total chance of exposure over that 30-year period, which includes climate change effects. The fire factor scoring rubric is included in Table 5.

7.2. Assumptions and Limitations

The wildfire hazard estimates from the methodology described in this research paper offer insights into the current and future wildfire exposure at 30 m resolution across CONUS, using widely accepted input layers from LANDFIRE and using the Rothermal-based ELMFIRE fire behavior model that has already undergone peer-review and validation. The resulting estimates of wildfire hazard exposure provide a first view of national level, high precision, property-level exposure estimates across the US in a framework that takes into account both current and future changing exposure to wildfire. The results identify at least some level of exposure in many places that are generally not thought of as having a wildfire problem, but they also underscore the fact that there is a tremendous amount of “major exposure” in the Western US, and specifically in the WUI areas in California and the Mountain West States. These insights are intended to complement the work carried out by the WRC program by providing a property-level equivalent to the community level tool already in the public domain, using similar but independent Rothermel-based fire behavior modeling. Nevertheless, there are a number of acknowledged limitations in our methodology, many of which have already been noted, but the implications of which are discussed in the following list:
  • Lack of explicit fire suppression: since the fire behavior model ELMFIRE does not explicitly include suppression effects, the model tends to overestimate the size and intensity of wildfires, which leads to an overestimate of the extent of wildfire exposure.
  • Variable length of wildfire burn time: ELMFIRE randomizes the length of the time for each modeled wildfire, leading to overestimates in the size and intensity of wildfires. The amount of time and the number of simulated fires needed to drive the Monte Carlo simulation towards stable statistics varies geographically across the model domain.
  • Extremely large fires: the simulation method does not capture the behavior of extremely large fires, since the fire weather forcing the simulation is not coupled with the fire behavior model.
  • House to house ignition: while the replacement of the non-burnable fuels in the LANDFIRE representation of the WUI with estimates of burnable fuels allows wildfires to propagate through the WUI more accurately, the ignition and subsequent contribution to wildfire by the buildings/houses themselves to the hazard within the WUI is not yet included in FSF-WFM.
  • Vegetation changes: the vegetation between 2022 and 2052 was held constant, although it is anticipated that changes in vegetation composition and density, and thus fuels, will be driven to some degree by climate change. Keeping the 2022 fuels constant for the assessment of 2052 future exposure underestimated the total possible changes due to the climate, but focuses attention on the direct effects of the climate and future weather on the state of those fuels, which has significant implications for wildfire ignitions, intensity, and spread.
  • Future weather approximation: A comprehensive sensitivity analysis to the bias-adjustment techniques used for climate adjustment is warranted. In addition, the high quality of the winds in the 2048–2057 simulations (the same as the 2011–2022 observations) is an advantage over using modeled winds, but is nevertheless an assumption. Most importantly, since the length and severity of droughts captured in the 2011–2020 time series do not change for the 2048–2057 simulation, the possible impact of those droughts, as they increase in frequency and severity, is unresolved.
  • Incomplete fuels/disturbances for fuel updates: disturbances are not evenly reported across the US, and some areas (e.g., private lands in the Eastern US) are not well known.
  • Ignition locations: using only historical fire ignition locations limits the possible impact of climate change on plausible fire locations, and the omission of random lightning strikes leaves some areas under-sampled. Additionally, a nuance of the decision to build the ignition density surface from only >100 acre fire occurrence data is that ignition density will be zero in areas that have not experienced fires >100 acres, even if those areas have experienced fires <100 acres. Dillon et al. [39] noted that in areas where management strategies have previously been successful at limiting large fire occurrence, burn probability modeling based only on large fire occurrence may underestimate burn probability. For that reason, Dillon et al. [39] developed an ignition density surface weighted as 98% large fire occurrence and 2% small fire occurrence, and such an approach could likely be used in future work.
  • No future land use changes: to focus on the impacts of climate change on the existing parcels under future wildfire exposure, we have elected to keep the built environment constant, and to assume no changes in land use or condition. This simplifying assumption is useful for its stated purpose, but we also recognize that changes in land use will also precipitate changes in likely future ignition locations, WUI locations, fuel conditions and types.

8. Discussion and Concluding Points

The methodology presented computes the physical hazard associated with wildfire incidence for the contiguous United States at 30 m resolution, and is expressed through hazards quantifying burn probability, flame length, and ember spread for the years 2022 and 2052, based on 10-year representative Monte Carlo simulations of wildfire behavior. This methodology uses updated fuels estimates that integrate known disturbances, current and estimated future weather characteristics that are useful for understanding aggregate wildfire exposure at a high resolution, and uses a model of wildfire behavior that integrates ignition, time of burn, and spread. This work does not develop new techniques or approaches to fire probability and hazard modeling, but rather integrates several existing methods and implements a computationally efficient and scalable modeling techniques to allow for new high-resolution, CONUS-wide hazard generation—all based on existing data, fire science, and hazard modeling paradigms developed by others in the wildfire science community. We have extended these approaches to estimate not only updated, current wildfire hazards but also extending those to estimate climate change’s future impacts on these hazards.
The methodology for the augmentation of the US Forest Service’s LANDFIRE-based estimates of fuel types, densities, and conditions at a 30 m resolution is presented using an open-source, Rothermel-based wildfire behavior model, ELMFIRE, for computation. The replacement of non-burnable fuel types in LANDFIRE that represent the built environment within the wildland–urban interface (WUI), with fuel inputs from the results of machine-learning estimates trained on data from historical fires, allow the propagation of wildfire through the WUI in a way that more closely resembles the observed conditions, and often results in non-zero burn probabilities for these areas. This serves as a notable improvement and opportunity for future fire models to replicate such an approach to improve their modeling. The wildfire hazard derivation overall is heavily dependent upon the updated LANDFIRE 2016 fuel layers, and significant effort was undertaken to assemble all known disturbances throughout 2020. Combined, this provides a repeatable methodology for future research looking to incorporate current fuel estimates, when annually updated LANDFIRE data are not available.
Other inputs required for ELMFIRE include topography from the USGS National Elevation Database, and weather (winds, air temperatures, humidity, and precipitation), for which the 2011–2020 NOAA RTMA hourly time series was selected. This 10-year time series provided an adequate range of possible weather conditions for the Monte Carlo simulation, where ELMFIRE was run approximately 100 million times to produce an estimate of the 2022 wildfire hazards for CONUS. To enable an estimate of the future hazard, this same hourly time series was bias-adjusted using MACAv2 daily downscaled IPCC CMIP5 RCP4.5 climate model ensemble results. Since accurate winds are crucial to the accurate prediction of wildfire behavior, and winds have a direct and significant influence on ELMFIRE results, we elected to hold winds constant between the 2022 and 2052 simulations, and bias-adjust only air temperature, humidity, and precipitation. This choice reduced the uncertainties introduced into the hazards from the fire behavior model, and instead focuses on the impact of climate change on the condition of the fuels for the 2048–2057 Monte Carlo simulations. Vegetation was likewise held constant between the 2022 and 2052 Monte Carlo simulations, as a reasonable but conservative approximation over 30 years’ time. The differences in wildfire hazards in the 2052 estimates are then based solely upon climate’s impact on the state of the fuels and generally hotter, drier conditions are thought to influence greater burn probabilities in the 2052 estimates. Due to the vegetation being held constant, these 2052 estimates should be considered conservative estimates of future wildfire exposure.
Fire ignition locations for the simulations were kept the same for 2011–2020, as for the 2048–2057 Monte Carlo simulations, and were created from the historical origins of significant fires greater than 100 acres. This lower limit on fire size was used to implicitly account for fire suppression activities that are not currently modeled in ELMFIRE. Over 100 million fires were modeled for each simulation period, and 8–10% of those model fires grew and were tracked at 30 m resolution across the landscape for up to 7 days apiece. Outputs were aggregated to create burn probability, flame length, and ember spread hazard estimates at 30 m horizontal resolution for CONUS. These hazard estimates are conducive to the assessment of the exposure of US properties to wildfire flames and/or embers. Comparisons with historical wildfire intensities and sizes show that the lack of explicit fire suppression effects in the FSF-WFM produces overestimates of fire sizes and intensities, so the resulting wildfire hazards should be considered to be conservative overestimates. Comparisons to historical wildfire losses and the US Forest Service’s WFC products generally show consistency at the state and community levels, but additional validation using historical losses at the building level should be undertaken in the future. The FSF-WFM wildfire hazards will produce fewer false negatives of risk assessments at the property level, and when combined with specific building vulnerability, could be used to provide similarly conservative estimates of climate-adjusted wildfire losses at the building level.
Wildfire hazards are estimated to be non-zero for 71.8 million of the over 140 million properties in CONUS, and will include an additional 11% properties over the next 30 years, due to climate change impacts on fuel conditions. While most of the overall wildfire risk is associated with properties west of 100 degrees W longitude in the American West, much of the change in wildfire exposure is observed east of the Mississippi River in areas not normally associated with large wildfire exposure. Over 5.9 million properties are found to have a “major” aggregate wildfire exposure of 10% over the 30-year analysis period from 2022–2052, which invites further investigation at the hyper-local level to discover ways to mitigate that exposure. Since the fuels and winds have been held the same between 2022 and 2052 in our simulations, the implication is that any increase in wildfire exposure is due to the future weather’s increased impacts on fuel conditions. Thus the influence of climate change on fuel conditions is the primary cause of the estimated increase in wildfire exposure throughout the country.
The FSF-WFM represents the first national-scale, property-level wildfire exposure model that has been developed using a geographically-consistent approach. The ability to consistently assess wildfire exposure, and thus risk for every property across the CONUS, should give local, state, and national government decision makers another data tool to help guide the allocation of resources, allow property owners to better assess their risk and implement meaningful solutions to reduce that risk, and provide financial markets with the opportunity to price risk into the cost of property more effectively through insurance, mortgage, and other financial products.

Author Contributions

Conceptualization, E.J.K. and D.S.; methodology, D.S., M.A., J.R.P., C.L., C.R.L., K.D.W., G.W.J. and O.M.D.; software, C.L., K.D.W., B.M., G.W.J., K.M.; validation, G.W.J., C.R., K.M., A.M., H.H. and M.B.; formal analysis, C.L., C.R.L., M.A.; investigation, B.W., J.R.P., B.M., M.B., K.L., H.H., A.M., F.C.; resources, J.R.P., M.A., O.M.D.; data curation, N.F., C.R., C.R.L.; writing—original draft preparation, E.S., J.R.P., E.J.K., C.R.L. and C.L.; writing—review and editing, E.S.; visualization, C.R.L., C.L., M.A., N.F. and J.R.P.; supervision, D.S.; project administration, E.J.K., M.A. and D.S.; funding acquisition, E.J.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding, and was entirely supported by the philanthropic donors to the First Street Foundation. Cloud compute credits in support of this research were graciously provided through Amazon Web Services by the Amazon Sustainability Data Initiative (ASDI).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The resulting FSF-WFM hazard layers and associated property-specific vulnerability and economic assessments will be freely and publicly available for noncommercial use at https://riskfactor.com, (accessed on 8 August 2022). The public availability of this climate information is meant to inform the public, enable new research efforts on wildfire risk, level the playing field with private commercial interests that already have access to this kind of information, and help address the privatization of climate impact information.

Acknowledgments

The authors are grateful for the philanthropic support provided to the First Street Foundation to enable this research. We are also thankful for the helpful comments from three anonymous reviewers of this manuscript. The First Street Foundation Wildfire Model is the product of a collaborative partnership between the First Street Foundation and the Pyregence Consortium, which includes the Spatial Informatics Group (SIG), Reax Engineering, and Eagle Rock Analytics. The Pyregence Consortium is grateful for the California Energy Commission’s support [EPC-18-026] of the Consortium that paved the way for this new study. The authors thank Jeff Knickerbocker (SIG) for his compute and hardware support, and Teal Dimitrie, Jean-Pierre Wack and Alecio O’Day for their project management of the Pyrgence Consortium’s efforts. Special gratitude is extended to the following individuals who provided valuable insights related to the FSF-WFM and its inputs: Janice Cohen of the National Center for Atmospheric Research; LeRoy Westerling and John Abatzoglou of the University of California Merced; Mark Finney, Gregory Dillion, Karen Short, and Eva Karau of the Rocky Mountain Research Station; Joel Reynolds, Marybeth Keifer, and Windy Bunn of the National Park Service; Ben Sleeter and Nathan Wood from the U.S. Geological Survey; Joshua Picotte, Tobin Smail, and Inga La Puma of the LANDFIRE team. The authors are grateful to Ana Pinheiro-Privette from the Amazon Sustainability Data Initiative for her support of their compute and data services workloads on AWS.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Data Sources Used in the Development of the FSF-WFM

NameSubjectSource
LANDFIREFuelsLandFire.gov (accessed on 1 June 2022)
USGS NEDTopographyhttps://www.usgs.gov/programs/national-geospatial-program/national-map (accessed on 1 June 2022)
USDA Forest Service’s Fire Occurrence Database (FOD)Ignition Locationshttps://doi.org/10.2737/RDS-2013-0009.5 (accessed on 1 June 2022)
MTBS, NIFCHistorical fireshttps://www.mtbs.gov and https://www.nifc.gov (accessed on 1 June 2022)
NOAA RTMAWeather, 2011–2020https://mtarchive.geol.iastate.edu/ and https://www.ncei.noaa.gov/has/HAS.FileAppRouter?datasetname=9950_01&subqueryby=STATION&applname=&outdest=FILE (accessed on 1 June 2022)
MACAv2ClimateMACA data portal (accessed on 1 June 2022)
Future 2052 WeatherWeather estimates for 2048-2057; Derived from MACAv2 and NOAA RTMA time seriesderived
Future 2052 FuelsAssumed to be same as 2022 for V1Held constant
Property BoundariesLightboxcommercial
Building FootprintsMapboxcommercial
Building DensityDerived from Building footprint information derived

Appendix B. Treatment Disturbance Inputs

Dataset NameSource Data Link
Hazardous Fuels Treatments—Fire—USFShttps://apps.fs.usda.gov/arcx/rest/services/EDW/EDW_HazardousFuelsTreatments_01/MapServer/3 (accessed on 1 June 2022)
Hazardous Fuels Treatments—Other—USFShttps://apps.fs.usda.gov/arcx/rest/services/EDW/EDW_HazardousFuelsTreatments_01/MapServer/4 (accessed on 1 June 2022)
Hazardous Fuels Treatments—Mechanical—USFShttps://apps.fs.usda.gov/arcx/rest/services/EDW/EDW_HazardousFuelsTreatments_01/MapServer/5 (accessed on 1 June 2022)
Hazardous Fuels Treatments—All Other Values—USFShttps://apps.fs.usda.gov/arcx/rest/services/EDW/EDW_HazardousFuelsTreatments_01/MapServer/7 (accessed on 1 June 2022)
Timber Harvest—USFShttps://apps.fs.usda.gov/arcx/rest/services/EDW/EDW_TimberHarvest_01/MapServer/8 (accessed on 1 June 2022)
CALMAPPER Treatment Projectshttps://egis.fire.ca.gov/arcgis/rest/services/CalMapper/CalMAPPER_Public/FeatureServer/2 (accessed on 1 June 2022)
CALFIRE Priority Treatment Projectshttps://services1.arcgis.com/jUJYIo9tSA7EHvfZ/ArcGIS/rest/services/PriorityProjects2019/FeatureServer/0 (accessed on 1 June 2022)
CALFIRE Timber Harvesthttps://services1.arcgis.com/jUJYIo9tSA7EHvfZ/ArcGIS/rest/services/CAL_FIRE_Timber_Harvesting_Plans_All_WGS84/FeatureServer/0 (accessed on 1 June 2022)
CALFIRE Fire Perimeters—Prescribed Firehttps://egis.fire.ca.gov/arcgis/rest/services/FRAP/FirePerimeters_FS/FeatureServer/1 (accessed on 1 June 2022)
CALFIRE Forest Health CCI Awarded Projects 2017-18https://services1.arcgis.com/jUJYIo9tSA7EHvfZ/ArcGIS/rest/services/Forest_Health_CCI_Awards_2017_2018/FeatureServer/0 (accessed on 1 June 2022)
Treatment areas—National Park Servicehttps://mapservices.nps.gov/arcgis/rest/services/WildlandFire/WildlandFire/FeatureServer/5 (accessed on 1 June 2022)
Treatment areas—ID—BLMhttps://navigator.blm.gov/api/share/e5fb96b234d32c5a (accessed on 1 June 2022)
Treatment areas & Harvest—OR Dept Forestryhttps://gisapps.odf.oregon.gov/data/FernsNoapsPolygons.Zip (accessed on 1 June 2022)
Treatment areas—CA—BLMhttps://navigator.blm.gov/api/share/a446b6874c2a37fc (accessed on 1 June 2022)
Treatment areas—NM—BLMhttps://gis.blm.gov/nmarcgis/rest/services/Range/BLM_NM_Vegetation_Treatments/MapServer/0 (accessed on 1 June 2022)
Harvest—WA Dept Natural Resourceshttps://gis.dnr.wa.gov/site2/rest/services/Public_Forest_Practices/WADNR_PUBLIC_FP_FPA/MapServer/6/ (accessed on 1 June 2022)

Appendix C. Fires Used to Estimate Fuels in WUI Areas

Incident NameStateIgnition DateLat.Long.Acres BurnedTotal Structures DamagedTotal Structures DestroyedTotal Structures Threatened
ShockeyCA9/23/201232.618−116.33526671045125
Bastrop County ComplexTX9/4/201130.13−97.23531,838017091160
Pine CreekOR7/14/201444.808−120.27331,0330016
Highway 613 FireMS10/31/201430.507−88.5266350030
Carlton ComplexWA7/14/201448.248−119.96276,09104711103
Mills CanyonWA7/8/201447.626−120.29721,95203571
AnacondaUT7/20/201440.562−112.23711420030
High RangeID8/3/201445.743−116.49353280330
Happy Camp ComplexCA8/14/201441.707−123.196118,49126767
Knf BeaverCA7/30/201441.89−122.87134,27406235
Snag CanyonWA8/3/201447.167−120.47512,508122279
Johnson BarID8/3/201446.096−115.61415,1700057
RainID8/3/201445.583−115.1854772004
Assayii LakeNM6/13/201436.032−108.84413,1760550
SlideAZ5/20/201435.009−111.80222,69800350
FrenchCA7/28/201437.294−119.3614,53400106
WayCA8/18/201435.735−118.461394712121500
EilerCA7/31/201440.799−121.55830,967030755
Taylor Mountain RoadUT7/5/201440.531−109.57329653350
Triple GFL5/9/201526.118−81.591736000
Grand LakeFL4/19/201525.75−80.45513680011
Lime HillOR8/5/201544.37−117.3312,210054
Dry GulchOR9/12/201544.829−117.13918,36900507
MannID8/18/201544.263−116.8415270030
Mm43 Hwy 52ID6/25/201543.977−116.411,0222010
CelebrationID6/6/201543.26−116.4977281000
SodaID8/10/201543.319−116.861282,88811145
Sleepy HollowWA6/28/201547.455−120.375323827350
I-90WA7/19/201547.013−119.95913970020
Highway 8WA8/4/201545.802−120.18435,29600350
Brown RanchTX8/11/201529.993−100.42817,8810322
County Line 2OR8/12/201544.829−121.41268,189071452
Roosa GapNY5/3/201541.638−74.42127470011
Pipeline 1PA5/3/201541.123−75.677666000
NorthCA7/17/201534.372−117.4744366523700
Gilmore GulchWA7/5/201546.16−116.96480740011
TucannonWA8/29/201546.359−117.678280900140
Ridge RoadND4/14/201548.079−103.093390000
PowerlineOK1/26/201535.364−95.88411830011
HighwayCA4/19/201533.907−117.624121200252
Z Bar 7OK3/31/201536.664−96.1495908000
2230 RoadOK4/4/201536.454−96.1582650060
Wf West End 2015TX2/13/201529.588−94.3416590000
Razor FirePA4/18/201540.779−75.682728004
Boars HammockFL4/26/201526.885−81.253790000
Tallgrass EastKS4/14/201538.41−96.5251745000
Wf Texas Point NortheastTX10/4/201529.705−93.934635000
GreenwoodOK3/23/201536.054−96.3195774000
West ProngOK3/24/201536.413−96.064367601500
Trail 12FL5/5/201528.788−82.3661041000
StationWY10/11/201542.882−106.1898459446392
Big Spring BranchWV11/17/201537.701−81.8251044000
Little Horse CreekWV11/17/201538.132−81.8511145000
Little JerrellWV11/18/201537.985−81.6461193000
Trace ForkWV11/14/201537.434−81.934784000
Kearny RiverAZ6/17/201533.068−110.9215435550
WillowAZ8/8/201534.837−114.54460844031710
GoodellWA8/11/201548.683−121.22766240050
Stouts CreekOR7/30/201542.859−122.98527,57000645
Route ComplexCA7/31/201540.601−123.54135,44402475
GrenadeCA4/29/201533.404−117.5141776000
River ComplexCA7/31/201540.914−123.36478,53100506
SolimarCA12/26/201534.303−119.342108300103
CuestaCA8/17/201535.356−120.612241501339
ParkhillCA6/20/201535.367−120.4241795518100
TassajaraCA9/19/201536.391−121.58910851210
LowellCA7/25/201539.212−120.8692633131800
TeslaCA8/19/201537.636−121.5942508010
LumpkinCA9/11/201539.527−121.327113700200
WraggCA7/22/201538.481−122.069845552700
RockyCA7/29/201538.91−122.4596,1258966959
ValleyCA9/12/201538.788−122.61377,5079520199150
RoughCA7/31/201536.852−118.884146,369041536
WashingtonCA6/19/201538.642−119.69918,48502251
ButteCA9/9/201538.266−120.59272,894489016400
CorrineCA6/19/201537.179−119.5106403250
WillowCA7/25/201537.282−119.479599000455
Cape HornID7/5/201547.998−116.5211505114309
SlideID8/14/201546.096−115.38213,5090029
I-90 SpragueWA8/1/201547.314−117.9341771002
Carpenter Rd.WA8/15/201548.05−118.09162,4880431005
Lawyer 2ID8/11/201546.23−116.10811,3780025
MunicipalID8/15/201546.469−116.191969511302
WoodratID8/11/201546.167−115.77165130081
Tepee SpringsID8/12/201545.318−116.11694,878061410
EagleOR8/11/201545.028−117.37314,5020152
Canyon Creek ComplexOR8/12/201544.301−118.85109,78610054722
Black CanyonWA8/14/201547.976−120.05361,379000
First CreekWA8/14/201547.929−120.24479712219556
Chelan ComplexWA8/14/201547.912−119.84621,7741552948
West Fork Fish CreekMT8/14/201546.909−114.80414,49505372
North StarWA8/13/201548.415−118.94218,547014225
Marble ValleyWA8/14/201548.404−117.89234312241326
RennerWA8/14/201548.758−118.19313,97500120
Blue CreekWA7/20/201546.037−118.085990012250
9 MileWA8/13/201548.971−119.296505201080
LimebeltWA8/14/201548.507−119.694137,0980020
Hidden PinesTX10/13/201530.081−97.18338072141406
Tunk BlockWA8/14/201548.478−119.339180,11101453000
Liberty HillLA10/13/201532.345−92.9077110312
LakeCA6/17/201534.147−116.76230,421047390
SunlandWA5/29/201647.045−119.9919400020
16 MilePA4/20/201641.199−75.1497896911287
Bear TownPA4/20/201641.181−75.222649000
Sams Point Fire-Verkeerder FireNY4/23/201641.681−74.3431929007
Road 10WA8/2/201647.23−119.35727500887
Elmer CityWA9/11/201647.978−118.942561901140
Rocky Mtn Fire 2016VA4/16/201638.31−78.665929900337
Fifteen MileOR7/1/201645.638−121.00640440045
Range 12WA7/30/201646.495−119.869167,60400250
County Line Road FireNC3/10/201635.011−79.5131704000
Mcbee CommandWA7/15/201646.249−119.51918130020
Cellar MountainVA3/17/201637.93−79.128737008
South Ward GapWA7/31/201646.177−119.825418402100
KahlotusWA8/22/201646.645−118.63393860430
StarbuckWA7/18/201646.527−118.0872414000
Eades HollowVA11/21/201637.775−78.85115640017
RattlesnakeOR7/24/201644.835−121.11892960041
Table RockID6/30/201643.591−116.131248112100
Cottonwood CaSD10/16/201643.905−101.86241,775020
CliftonID8/23/201642.159−112.0123560020
Henrys CreekID8/21/201643.447−111.76552,98808125
SalvageID6/24/201642.8−114.66918471011
RockNV7/29/201639.871−119.896238700800
MetzCA5/22/201636.387−121.2173826005
Bug CreekAZ6/28/201634.294−112.117118400105
LongviewAZ6/6/201631.633−110.54811050040
RidgeAZ5/25/201631.529−110.33813910053
CrutchTX3/23/201635.615−101.11745,052040
OptimaOK12/16/201636.692−101.095084000
PoplarNC3/31/201536.106−82.337768003
Chestnut KnobNC11/6/201635.619−81.657641800417
HortonNC11/22/201636.146−81.568148000325
Bench BluffTN11/12/201635.594−85.2421715000
Pinnacle MountainSC11/9/201635.055−82.7217869011136
Rd 80KS3/17/201638.164−96.39163,061000
Bar-Dew LakeOK3/19/201636.828−96.03914,8060062
BearOK2/17/201636.309−96.1759660080
Pawnee CoveOK2/18/201636.228−96.4113418050200
PharoahOK2/18/201635.535−96.09713,57901450
Sand CreekOK2/18/201635.3−96.104383900170
Double HeaderOK3/6/201635.35−96.08314130015
KatieOK3/6/201635.258−96.09915732035
MustangOK3/6/201636.669−96.02510,0600082
Hall HornOK3/16/201636.557−96.2985180000
WalkerOK3/19/201636.354−96.2352340000
VarsityOK4/7/201635.728−96.43912160055
BurmacKS3/23/201638.087−97.66810,66812110
Burley HillKS4/5/201639.016−96.64516,381100
Quinton FireOK2/17/201635.14−95.3991219005
Mason FireOK2/19/201635.261−95.4731915000
Round Prairie Road FireOK2/19/201634.677−95.1732858000
CyclopsAL11/5/201633.845−87.033660004
Mount PleasantVA11/19/201637.738−79.17811,00100120
GapCA8/27/201641.856−123.03633,940014160
WillardCA9/11/201640.378−120.749282807625
Kewa FireWA8/2/201648.183−118.28419850590
Cayuse MtnWA8/22/201647.847−118.03897440231535
HartWA8/21/201647.822−118.12528192040605
WhitWY8/2/201644.409−109.36112,73108165
Cliff CreekWY7/17/201643.296−110.38236,13101135
ChimneyCA6/1/201635.855−118.02514770024
ErskineCA6/23/201635.569−118.33448,066752862500
CedarCA8/16/201635.791−118.57129,1910122599
SlateCA10/4/201636.082−118.5562121000
MeadowCA10/30/201635.975−118.5794346000
PioneerID7/18/201644.139−115.585189,59606465
Roaring LionMT7/31/201646.177−114.24880963662347
Spokane ComplexWA8/22/201647.492−117.2896839217303
SherpaCA6/15/201634.497−120.033754905271
ReyCA8/18/201634.586−119.72533,32355301
PilotCA8/7/201634.308−117.2478267005600
Blue CutCA8/16/201634.324−117.50636,8568321611
BogartCA8/30/201633.986−116.933147502426
FishCA6/20/201634.181−117.939452800869
DeerCA7/1/201635.222−118.688188500300
SageCA7/9/201634.366−118.5741002102500
Little ValleyNV10/14/201639.266−119.8392964040200
ColdCA8/3/201638.536−122.07762890252
TrailheadCA6/28/201638.963−120.835743002600
SoberanesCA7/22/201636.322−121.701132,3805682010
SandCA7/22/201634.391−118.3541,56162010,300
ClaytonCA8/13/201638.915−122.5873792293021500
LomaCA9/26/201637.116−121.8184380128325
Border 3CA6/19/201632.611−116.57279583171000
MormonAZ5/15/201634.961−111.5737897000
GooseCA7/30/201637.015−119.466248719400
TenderfootAZ6/8/201634.232−112.708436303300
CurryCA7/2/201636.087−120.4528370025
ChimneyCA8/13/201635.738−121.07546,95024701898
JuniperAZ5/20/201633.864−110.92632,29310141
ElkAZ7/21/201634.174−109.8641965000
Beaver CreekCO6/19/201640.957−106.50544,221017131
FultonAZ9/12/201634.281−110.89323700296
Beulah HillCO10/3/201638.07−104.9285769014750
JunkinsCO10/17/201638.14−105.13619,023026745
TopockAZ4/6/201634.74−114.5114221012
I40TX3/23/201635.244−100.35514,780014200
350 ComplexOK4/5/201636.647−99.26658,05500825
Anderson Creek FireOK3/23/201637.107−98.835374,52305410,000
Big CreekMO2/13/201636.634−92.834031000
Bob WhiteWV4/3/201637.957−81.78240010
Upper Conley HollowWV4/4/201637.889−82.0951438000
Jimmie Creek RdKY10/26/201637.373−82.386500000
Raven RockVA11/2/201637.172−82.61227312104
Bridge Creek RoadTN10/30/201635.241−85.5591777000
Spruce Pine Rd St Rt 7KY4/13/201637.522−82.894786000
Little Shepherd TrailKY10/26/201636.951−83.1136751000
Big BranchKY11/17/201637.062−82.956762000
Poe RoadTN11/11/201635.261−85.254758000
MowbrayTN11/9/201635.286−85.2087210050
Bolts Br.KY11/24/201637.087−83.6561069000
Lake ChinnabeeAL11/28/201633.473−85.87112542035
Caney HeadAL3/20/201633.391−85.843957000
Halls TopTN4/4/201635.875−83.14224641171
Eagles NestKY11/2/201637.53−83.39228570050
Jetts Creek FireKY11/6/201637.508−83.5633021000
Bowlings CreekKY11/21/201637.361−83.431023000
Moore PeachVA4/10/201636.625−82.9861345000
Sr116TN11/3/201636.169−84.3182222000
Charles Branch LaneTN11/8/201636.204−84.3451071000
Timber RidgeGA11/12/201634.828−83.363100200400
Neddy Mountain RoadTN11/11/201635.948−83.07278800150
Silver MineNC4/21/201635.908−82.79160820015
Party RockNC11/5/201635.472−82.2418572031050
State LineTN4/16/201635.926−82.92211110017
Sr116-DevoniaTN11/4/201636.116−84.413077000
TellicoNC11/3/201635.299−83.58914,17201336
Maple SpringsNC11/4/201635.392−83.93276961029
Dick’s CreekNC10/23/201635.399−83.2498330031
Dobson 3NC11/8/201635.506−83.2447410050
Hwy 190KY11/2/201636.74−83.724957000
Railroad Grade 2016TN4/18/201636.227−82.11317900025
Kentucky RidgeKY11/7/201636.686−83.8581133000
Old RoughyNC11/9/201635.371−83.855340044
East Miller CoveTN11/17/201635.744−83.799133100100
Quarry CreekTN11/16/201635.349−84.281643009
Cobbly NobTN11/28/201635.779−83.342732231080
Stinking CreekTN11/9/201636.453−84.19910,768000
BotelerNC10/25/201635.068−83.673862600314
KnobNC11/2/201635.114−83.5371132000
Camp BranchNC11/23/201635.179−83.558323402140
Wild GooseLA2/9/201631.406−92.89812600320
Chimney Tops 2TN11/23/201635.687−83.50314,99825720662800
Knox Bell LineKY10/29/201636.907−83.6051272000
Rock MountainGA11/9/201634.99−83.52225,22400250
North PeakNC3/23/201635.753−81.986680002
Old 50FL9/27/201628.547−80.906804001
IslandFL5/5/201629.318−81.767527000
SkiboMN5/6/201647.498−92.04476300110
Clear CreekNC11/20/201635.72−82.113349300392
TombstoneNC3/8/201635.545−81.7251747000
Chetco BarOR7/12/201742.238−124.049194,87793012,483
HelenaCA8/31/201740.775−123.06218,70981415350
CanyonCA9/25/201733.861−117.662740601910
Canyon 2CA10/9/201733.823−117.734910258265000
Minerva 5CA7/29/201739.903−120.944454500395
DetwilerCA7/16/201737.55−120.12183,297211311500
GateCA5/20/201732.654−116.829226500315
RailroadCA8/29/201737.441−119.61312,765019511
Lilac 5CA12/7/201733.299−117.2034159691931500
MissionCA9/3/201737.236−119.466100689250
EarthstoneNV7/3/201739.591−119.51735,29901131
PreacherNV7/24/201738.855−119.588533000800
PraterNV8/6/201739.551−119.66915720030
Long ValleyCA7/11/201739.992−119.9280,456310500
Winnemucca RanchNV7/4/201739.756−119.644415305300
Cold SpringsNV7/14/201739.646−119.938155700100
OperaCA4/30/201733.997−117.3011070000
CutterNV10/3/201738.829−119.604406500150
RoadrunnerCA7/30/201736.015−118.93324360010
PierCA8/29/201736.122−118.70836,626021360
WinterCA7/6/201738.526−122.05424850063
AtlasCA10/9/201738.364−122.23751,6641297905000
JonesOR8/11/201744.004−122.51210,260015
NunsCA10/9/201738.349−122.50356,883000
HatcheryNV7/4/201738.977−114.0911142000
CreekCA12/5/201734.294−118.35215,833811232500
East ForkMT8/27/201748.231−109.57621,1655580
RyeCA12/5/201734.43−118.6354895365460
Gibralter RidgeMT8/8/201748.86−114.849629901145
CaribouMT8/11/201748.979−115.35128,101040570
West ForkMT8/30/201748.519−115.60621,15400709
Canyon CreekWA7/15/201748.271−120.07212320285
East SaddleWA8/12/201746.782−119.34917,3182240
WallCA7/7/201739.463−121.406648812915400
CascadeCA10/9/201739.359−121.37516,15502001000
CherokeeCA10/9/201739.591−121.58584150353
TubbsCA10/9/201738.568−122.6836,9811457629,192
SulfurCA10/9/201738.992−122.6662591702051720
Oil WellNV7/17/201740.92−115.7257240028200
Redwood Valley IncidentCA10/9/201739.339−123.21336,545090100
PocketCA10/9/201738.77−122.88318,691000
PalmerCA9/2/201733.991−117.121414810150
Mecca FireOR6/26/201744.781−121.23125150023
Emerson 0638 RnOR7/25/201744.69−121.01210,6830125
ThomasCA12/4/201734.459−119.303281,982280106318,000
Nena SpringsOR8/9/201744.974−121.19870,074410199
Pilot ValleyNV8/13/201741.106−114.095257826157
AnaOR7/8/201743.009−120.76958012455
Eagle CreekOR9/2/201745.618−121.94248,816095526
Sheep GapMT8/29/201747.475−115.04624,7020080
Silver DollarWA7/2/201746.57−119.77930,7890030
Horn Butte 0594 RnOR7/21/201745.682−120.06993250035
Glade 3WA7/30/201746.144−120.05410,5820120
Morgan CreekOR8/3/201744.426−117.2352329000
Martin CanyonID7/23/201743.49−114.1684053000
LagoonID7/26/201742.955−114.43814840310
Mammoth CaveID8/4/201743.155−114.19350,3910350
BreezeID6/26/201743.302−115.8718630040
North DelphiaMT7/14/201746.548−108.27837670040
Sage HillsMT7/20/201745.763−108.348119700100
Lincoln BeachUT6/23/201740.071−111.84222980014
MulberryAZ5/6/201731.901−110.61118460420
LizardAZ6/7/201731.986−110.00615,79100108
EncinoAZ6/21/201731.65−110.6481357215250
CajeteNM6/15/201735.809−106.559143300233
SawmillAZ4/23/201731.822−110.68747,35700415
Alice CreekMT7/22/201747.142−112.43829,97104240
SunriseMT7/17/201747.07−114.83826,89600382
TarinaCA6/30/201735.385−118.7931257006
Lolo PeakMT7/15/201746.666−114.24262,3162101962
MendenhallMT8/26/201745.653−110.1811960230
JulyMT7/3/201747.89−108.57511,40906101
HonditoNM5/15/201736.608−106.026949000
HillCA6/26/201735.405−120.48119005430
WolfSD3/4/201743.998−102.16917970045
Hodgeman CountyKS3/6/201738.143−99.8538518980
South WenasWA6/27/201746.72−120.601295600175
SpartanWA6/26/201747.32−120.15587750091
SheepWA7/24/201746.758−120.5441564000
Monument HillWA8/17/201747.304−119.71864371923175
MeyersMT7/14/201745.989−113.55268,71101344
Jolly MountainWA8/11/201747.341−120.97838,159005624
Rattlesnake HillsWA7/6/201746.515−120.43235530030
PerrytonTX3/6/201735.99−100.36290,211011200
MonitorWA11/1/201747.503−120.4119600300
BrianheadUT6/17/201737.788−112.69374,2765261526
Thirty SevenCA10/9/201738.155−122.47417734080
303TX2/28/201733.401−102.5359601000
SlinkardCA8/29/201738.659−119.571881400510
KeystoneWY7/3/201741.174−106.28127840180
TrippTX2/10/201734.448−100.7892573004
Dumas ComplexTX3/6/201735.357−101.72226,15500150
PrisonTX2/28/201734.522−101.80424204131143
2018 North Sargent WfTX10/16/201728.813−95.6239590060
Oks—283OK3/7/201736.689−99.75568,55800300
Lefors EastTX3/7/201735.365−100.52468,701000
Beaver MountainOK1/31/201735.168−95.3394754000
Powder MillOK2/2/201734.964−95.38215510047
HighlandsKS3/7/201738.182−97.91974189131100
Jupiter HillsKS3/4/201738.109−97.848128311100
Legion LakeSD12/11/201743.669−103.39454,86803203
Sugar CoveNC1/28/201735.75−82.1436380014
Turn Table FireSC4/2/201733.312−79.8781868000
Dobson KnobNC4/9/201735.812−81.99317200045
Big Branch FireKY4/9/201737.181−83.0546510032
Ne 212th StFL3/31/201729.469−81.9566100040
Sod Farm 2FL4/16/201728.913−81.453901000
Lost CreekOK3/3/201735.403−96.14124200025
SpocogeeOK3/1/201736.071−96.3363180045
Gun RangeOK3/21/201735.646−96.0631524003
Cod DrFL7/8/201728.767−82.2586400012
ConnerFL3/26/201729.252−81.9186760010
310 West Of ComoMS1/29/201734.518−90.075745216
Cr630 EFL2/15/201727.782−81.315509601420
Oks—StarbuckOK3/7/201737.081−99.893657,299001000
BonitaNM6/3/201736.58−106.14977540065
Garfield RoadFL3/22/201730.418−82.02272114213
West MimsGA4/6/201730.651−82.294166,73704920
AppleCA6/9/201839.924−122.3492849090
CreekCA6/24/201840.486−122.5181353011610
Middle RidgeOK3/21/201735.577−94.6268501000
Persimmon RidgeOK3/21/201735.653−95.0885333000
SunCA10/7/201840.23−122.14339210070
Lost FireOK3/21/201734.708−95.7574178024
Potato HillsOK3/23/201734.694−95.2262503000
MontecitoWA6/28/201846.175−119.76218770050
Wagon WheelWA9/1/201846.35−119.51340630090
Milepost Twenty TwoWA6/20/201846.966−120.0574060016
BofferWA8/11/201846.142−119.1414645070
ConradWA7/1/201846.739−120.665461101220
Milepost 90WA8/1/201845.681−120.93710,7570070
Lee Williams RdFL3/5/201726.133−81.6377288161000
South Valley RoadOR8/1/201845.371−121.16120,471019212
Jackson RanchOK3/23/201735.797−96.21742430040
Substation 0730 RnOR7/17/201845.5−120.93969,1098521363
Mile Marker 44WA9/1/201846.149−120.5294063000
Boxcar 0410 RnOR6/21/201845.022−121.00499,8740055
Tenino FireOR8/16/201844.708−121.3718821000
Graham 0420 OdOR6/21/201844.55−121.42102011204
Angel SpringsWA8/2/201847.775−118.0294718014170
EagleCA7/13/201841.268−120.10521160018
Soap LakeWA6/11/201847.436−119.4921580035
Chelan HillsWA7/27/201847.782−119.962185048100
Rocky ReachWA7/13/201847.527−120.327334600313
BoylstonWA7/19/201846.85−120.1166,292161
KeithlyID7/25/201844.464−116.84317,5880031
Silver StateNV7/14/201840.886−115.6653766000
RockyNV6/23/201840.378−118.26216411110
OwyheeNV7/21/201841.948−116.07753470060
South SugarloafNV8/17/201841.716−116.019241,426317116
GoodwinAZ6/24/201734.381−112.29928,1923331400
La TunaCA9/1/201734.23−118.31670351101376
PowerlineID8/4/201742.699−112.60554,3780135
White CreekNC3/16/201735.837−81.8834166007
ShoestringID8/5/201742.878−114.57435,543012
Penn Swamp FireNJ7/20/201739.677−74.6383587000
WeogufkeeOK3/20/201735.233−95.9042226000
HolidayFL4/5/201726.003−80.46888580020
30th AveFL4/20/201726.177−81.605646318143884
PowerlineWY8/12/201844.461−108.9261837000
Raintree BlvdFL5/13/201727.074−82.0523319000
Flat Rock FireNY7/12/201844.87−73.637658002
Tye RiverVA5/3/201837.905−79.15417610029
Spring CreekCO6/27/201837.543−105.144107,1081192252878
BlaineCA8/13/201733.984−117.291117460441
RooseveltWY9/15/201843.06−110.38755,3301571153
8 MileID9/22/201842.595−111.5410010595
MilesOR7/16/201842.828−122.69940,343021011
Ramsey CanyonOR8/22/201842.587−122.992212701540
NatchezCA7/15/201841.895−123.56638,80000104
Taylor CreekOR7/15/201842.488−123.61957,505003292
KlondikeOR7/16/201842.418−123.873178,311001940
FergusonCA7/13/201837.635−119.80797,3070115236
KerlinCA9/4/201840.625−123.512177505100
HirzCA8/9/201840.984−122.27946,70001171
DeltaCA9/5/201841.007−122.46263,732745330
RiverCA7/27/201839.055−123.01948,92002305
CountyCA6/30/201838.683−122.15592,4504311516
NelsonCA8/10/201838.312−121.999220511260
WhalebackCA7/27/201840.627−120.82418,64000460
HatCA8/9/201840.998−121.489197100380
BoydsWA8/11/201848.632−118.1545196010529
Crescent MountainWA7/29/201848.384−120.44653,258001196
CampCA11/8/201839.748−121.565153,68775118,83817,500
StoneCA8/15/201841.425−121.01339,45502119
ChavesNV6/3/201839.3−119.41236521098
Upper ColonyNV6/17/201838.812−119.41112550092
DonnellCA8/1/201838.383−119.82236,1510136305
Rattlesnake CreekID7/23/201845.232−116.376846100618
Rabbit FootID8/2/201844.848−114.2333,787001446
PineryAZ5/12/201831.982−109.35314740025
ViewpointAZ5/11/201834.693−112.34353891400
TinderAZ4/27/201834.59−111.1116,0830961700
Hub PointAZ7/27/201834.251−110.294674000
Soldier CanyonNM6/7/201833.185−105.759138600100
PiersonNM4/17/201832.499−103.4181060002
Ute ParkNM5/31/201836.533−105.02830,1773152952
ValleyCA7/6/201834.105−116.946125000500
Harbor BayTX4/13/201835.619−101.63114280850
CranstonCA7/25/201833.715−116.70513,0966126230
HolyCA8/6/201833.704−117.46822,845182413,300
StoneCA6/4/201834.55−118.292165900150
Badger HoleCO4/17/201837.433−102.08849,1461240
CharlieCA9/22/201834.521−118.559336700100
Cr 26TX4/14/201835.256−100.0841386006
MillironTX4/13/201834.913−100.01120,43702175
34 ComplexOK4/12/201836.585−99.35257,533055150
HillCA11/8/201834.207−118.953442724437
FrontCA8/19/201835.119−120.0971126005
PerryNV7/27/201839.802−119.49653,734316418
AirlineCA6/4/201836.391−120.9621477001
Lake ChristineCO7/3/201839.419−107.03712,506961329
ChateauCO6/29/201838.815−105.3141408754
Carson MidwayCO3/16/201838.527−104.7164773020
RheaOK4/12/201836.003−99.003277,9490503500
OrganNM6/24/201832.454−106.5274880001
Harman RoadTX7/18/201831.333−97.972309401115
Owl CreekNV8/30/201840.658−115.52611650012
Lime Rock Rd (19)FL6/24/201829.765−84.8661190436400
Cougar CreekWA7/28/201847.816−120.47842,681003000
RozellMO2/15/201836.587−92.8621970060
Bald MountainUT8/24/201839.925−111.69321,016012600
Pole CreekUT9/6/201839.982−111.53102,426012626
Hill TopUT8/6/201839.737−111.441178424320
Dollar RidgeUT7/1/201840.111−110.87769,81764531023
Little Shepherds TrailKY5/1/201836.976−83.034541000
Keepers Branch FireSC3/4/201833.189−79.542819006
Range TwoNV9/30/201840.669−115.44893610860
Dog HeadNM6/14/201634.851−106.319,8160731950
RanchCA7/27/201839.269−122.775427,0480086
416CO6/1/201837.493−107.90355,123003386
Grass ValleyWA8/11/201847.94−119.16676,074120330
West 60OK3/7/201836.789−96.48718,7150010
CarrCA7/23/201840.715−122.593233,71028216085013
Badger CreekWY6/10/201841.055−106.1120,75243553
TomahawkAR4/12/201836.056−92.6715330815
WoolseyCA11/8/201834.125−118.82497,962365164357,000
Flag Pond (11)FL3/21/201826.139−81.5952562000
Buffalo CorralAZ7/14/201931.573−110.3771144000
KincadeCA10/23/201938.672−122.77677,7856037590,015
EasyCA10/30/201934.267−118.8292105112635
SandalwoodCA10/10/201933.999−117.084104816760
TickCA10/24/201934.451−118.3944932483110,425
MariaCA11/1/201934.314−119.06610,036052722
SaddleridgeCA10/10/201934.318−118.5159656933825,760
Black BridgeCO4/4/201938.074−103.1771690007
116th Ave Se (11)FL3/21/201826.034−81.5629,262030
335TX4/13/201832.676−99.7162481015100
Blue Creek #2OK3/14/201835.052−95.5942578000
DeerteOK3/15/201835.367−95.8711833008
HenryOK3/20/201835.581−96.4211600017
BrewsterOK3/16/201835.685−96.29940490020
Flying GOK3/12/201836.124−96.203188600184
WalkerOK3/24/201836.35−96.2332626009
DrumbOK3/24/201836.467−96.330,419009
Onion PrairieOK3/6/201836.641−96.01739110020
New Years WfTX1/1/201829.706−93.90856130012
Farmers RoadTX1/22/201832.738−97.596237900252
CarbonTX4/13/201835.214−100.06812,14800500
PembertonAZ8/6/201934.731−112.6761211000
Green RavineUT9/3/201940.683−112.2282260003
Goose PointUT8/21/201940.076−111.83291900023
DoveAZ5/24/201933.799−112.42910780010
White WingAZ5/30/201933.796−112.43717970030
TenajaCA9/4/201933.549−117.2581820201200
WalkerCA9/4/201940.09−120.58558,7520978
Long ValleyCA8/24/201939.882−119.99524514180
BriceburgCA10/6/201937.608−119.932555501160
WoodburyAZ6/8/201933.52−111.175130,243001537
Stuckey Rd MaMT9/2/201947.586−111.3574077090
MatsonWA10/7/201946.619−118.9387150086
SandCA6/8/201938.904−122.259247307125
MuseumAZ7/21/201935.263−111.622011000
Elmer CityWA6/23/201948.064−118.94319962245
Desert CanyonWA7/23/201947.711−120.14815050016
243 CommandWA6/4/201946.855−119.78518,8911236
Cut AcrossMT4/7/201945.63−106.7185800130
BoulderCA6/5/201935.32−119.931199000
DeckerCO9/8/201938.441−105.995987604142
BoulinAZ8/6/201935.382−112.03640940012
ColdwaterAZ5/30/201934.481−111.33516,8240014
CaveCA11/25/201934.489−119.76827611015,000
Left HandWA7/23/201946.915−120.975323400347
North HillsMT7/26/201946.765−111.944414400600
Tx Point East Christmas EveTX12/24/201929.716−93.911630043
Williams FlatsWA8/2/201947.977−118.49844,6800356
Pedro MountainWY8/24/201942.337−106.82621,9100890
Cove CreekID8/4/201945.346−114.46552730039
ChanningTX2/16/201935.663−102.27578550020
BurnsideOK3/19/201934.593−97.02911170010
East Kennedy CreekKS4/2/201938.29−95.8151003000
Clark Branch #2KY9/17/201937.603−82.6211315000
344 DFL9/11/201930.356−84.6322339000
Cr 2297 Allenton (03)FL3/30/201930.115−85.4827250052
Dry HollowWV11/28/201938.831−79.30814120016
Kennedy PeakVA11/14/201938.762−78.466769000
Spring Hill FireNJ3/30/201939.79−74.4518182000

Appendix D. 2052 Weather Time Series Estimation

We used these publicly-available online repositories of NOAA RTMA data:
  • Iowa state (Iowa State) [57] covering 2011-current
  • National Centers for Environmental Information (NCEI) National Digital Guidance Database (NDGD) (NOAA NCEI, 2022) covering 2011- current
Since approximately 3% to 4% of hourly records are missing from both the Iowa State and NCEI datasets, we merged these datasets to minimize the amount of missing data.
Mid-century fire weather inputs were generated by scaling or “nudging” 2011–2020 Real Time Mesoscale Analysis (RTMA) temperature, relative humidity, and precipitation by computing distribution changes between 2021 and 2051 ensemble regional model output. While wind is an essential fire weather variable, it was excluded from nudging because the literature does not support a statistically robust climate change signal for it (e.g., Torralba et al., 2017 [58]).
We used the Multivariate Adaptive Constructed Analogs (MACA) version 2 as our source for 2020 and 2050 calibration data because it demonstrates skillful performance for fire weather variables, particularly relative humidity and wind (Abatzoglou and Brown, 2012 [33]). MACA is statistically downscaled Coupled Model Intercomparison Project version 5 (CMIP5; Taylor et al., 2012 [32]) available in a ~4-km grid over the contiguous United States (CONUS). It is available at daily resolution and at temporal ranges for present-day conditions (1950–2005) and future experiments under Representative Concentration Pathways (RCPs) 4.5 and 8.5. RCP4.5, which represents an intermediate scenario, was used for this experiment. Temporal ranges for the “2020” and “2050” calibration data were defined as 2011–2020 and 2048–2057 ensemble means, respectively. Note that the CCSM4 MACA ensemble member was not included in order to maintain consistency across variables, because it was not available for relative humidity at the MACA Data Portal (https://climate.northwestknowledge.net/MACA/data_portal.php, accessed on 8 August 2022).
Nudging was performed independently at each geographic grid point. In order to facilitate this, the 4-km rectilinear MACA data had to be interpolated to the 2.5-km curvilinear RTMA grid. The National Center for Atmospheric Research Command Language (NCL)’s regridding package from the Earth System Modeling Framework was used (https://www.ncl.ucar.edu/Applications/ESMF.shtml, accessed on 8 August 2022). We used bilinear interpolation for temperature and relative humidity, and conservative regridding for precipitation since it does not represent a smoothly varying field.
Nudging code was obtained from the bias_correction Python3 package (https://github.com/pankajkarman/bias_correction, accessed on 8 August 2022). Bias correction techniques like those employed here are used to correct modeled projections by adjusting them with statistical scaling factors derived between historical model output and observations; put another way, they are used to downscale climate change projections to station locations (e.g., Fang et al., 2015 [59], Luo et al., 2018 [60]). Likewise, we use them here to nudge RTMA output with respect to scaling factors computed between 2020 and 2050 MACA output. We summarize the techniques here.

Appendix D.1. Precipitation

RTMA precipitation projections were obtained via gamma distribution mapping, a technique developed by Switanek, et al., 2017 [61]). Gamma distribution mapping is unlike other bias correction methods because it accounts for the frequency of rain days and the likelihood of events (Switanek, et al., 2017 [61]).
After summing hourly RTMA precipitation into daily accumulations to match the temporal resolution of the MACA data, this procedure goes as follows:
  • Define rain days (RDs) as those which feature non-zero precipitation accumulations for each distribution. Then compute the expected number of projected 2050 RTMA rain days (RDP):
    R D P = R D R T M A × R D M A C A F R D M A C A H  
    where RDRTMA, RDMACAF, and RDMACAH are the number of rain days in RTMA, 2050 MACA, and 2020 MACA data, respectively. In the event MACA 2050 has fewer RDs than MACA 2020, this results in a downward adjustment in RDs in the nudged RTMA distribution. Note that projected RTMA RDs can never increase above original RTMA RDs, an important limitation in this approach.
  • Fit gamma probability density functions (PDFs) to each of the modeled rain-day distributions. Cumulative distribution functions (CDFs) and their inverses (ICDFs) are found from these PDFs.
  • Calculate the relative scaling (SFR) between the fitted RTMA and 2020 MACA distributions at all CDF values corresponding to the precipitation events of the RTMA time series:
    S F R = I C D F R T M A ( C D F R T M A ) I C D F M A C A H ( C D F R T M A )  
    where SFR is an array of relative scaling factors, ICDFRTMA and ICDFMACAH are the ICDFs for the fitted 2020 RTMA and MACA distributions, respectively, and CDFRTMA is the estimated CDF for RTMA precipitation values. As an example, say the largest value in the RTMA time series is equal to 40 mm and corresponds to a CDF value of 0.997 (in other words, ICDFRTMA(0.997) will yield 40 mm), but ICDFMACAH(0.997) yields only 6 mm (we discuss these systematic magnitude discrepancies in Assumptions and Limitations). The most extreme nudged value will have a relative scaling factor of 6.67 (40/6).
  • Calculate recurrence intervals (RIs) from the three fitted CDFs, then find the adjusted RI for 2020 RTMA:
    R I S C = m a x ( 1 , R I M A C A F × R I R T M A R I M A C A H )  
    where RIMACAF, RIRTMA, and RIMACAH correspond to the RIs for the RTMA and 2050 and 2020 MACA data, respectively. The value is always greater than or equal to 1 to ensure that the subsequent step yields values between 0 and 1. This step adjusts the RI of 2050 MACA events by differences in the extremity of 2020 RTMA and MACA events. For example, if the RI for the most extreme 2020 RTMA value is shorter than that of the most extreme 2020 MACA value, then the RI for the most extreme 2050 MACA value will be shortened accordingly. RISC is then used to find the corresponding scaled CDF values with
    C D F S C = 1 1 R I S C  
    where CDFSC reflects the scaling of the change in MACA event likelihoods between 2050 and 2020 with respect to RTMA likelihoods.
  • The initial array of projected RTMA values is:
    R T M A I = I C D F M A C A F ( C D F S C ) × S F R  
    where ICDFMACAF is the ICDF for the 2050 MACA data. In the event that (A1) adjusted the RDs for the projected time series, RTMAI is linearly interpolated along the adjusted length of RDP.
  • The array of projected values is then placed back into the RTMA time series at the corresponding locations. For example, if the largest annual RTMA value occurred on April 21st, then the largest projected RTMA value will be placed on that same day.
  • Last, the ratio between the projected and original RTMA values for each day is applied to the hourly RTMA precipitation time series.

Appendix D.2. Temperature

RTMA temperature projections were obtained via a normal distribution mapping technique (Switanek, et al., 2017 [61]). After applying cubic spline interpolation to daily minimum and maximum MACA temperatures to obtain hourly distributions (we discuss this further in Assumptions and Limitations), the steps go as follows:
  • All three distributions are linearly detrended to get a more accurate measure of the natural variability. While trends are added back at the end, all subsequent steps use the detrended time series.
  • Gaussian PDFs are fit to each of the annual distributions, then corresponding CDFs and ICDFs are derived from these.
  • An absolute scaling factor is found:
    S F A = [ I C D F R T M A ( C D F R T M A ) I C D F M A C A H ( C D F M A C A H ) ] × ( σ M A C A F σ M A C A H )  
    where σ MACAF and σ MACAH are the standard deviations of the 2050 and 2020 MACA data and all other terms are defined as before.
  • RIs are calculated from the CDFs, scaled with (A3), then used to find the modified CDF with
    C D F S C = 0.5 + s g n ( C D F M A C A F 0.5 ) × | 0.5 1 R I S C |  
    Equation (A7) differs from (A4) because the normal distribution is two-tailed.
  • The initial array of projected RTMA values is
    R T M A I = I C D F M A C A F ( C D F S C ) + S F A  
  • As before, reinsert projected values into the original RTMA time series. Then add the trend of the RTMA time series back into the projected one.

Appendix D.3. Relative Humidity

Mid-century relative humidity projections used a comparatively simpler approach called empirical quantile mapping (Gudmundsson et al., 2012 [62]). Unlike the previous two methods, it makes no assumptions about the underlying distribution of the data. Briefly, it goes as follows:
  • An empirical CDF (ECDF) is approximated for the 2020 MACA distribution.
  • The corresponding percentiles from the ECDF are found in the RTMA data.
  • The difference between the value at each equivalent percentile is computed between the 2050 and 2020 MACA data.
  • This difference is applied to the RTMA data.
In this way, the order of observations in the 2020 record is maintained (i.e., underlying variability) while shifting towards a future climate.

References

  1. Burke, W.D.; Tague, C.; Kennedy, M.C.; Moritz, M.A. Understanding How Fuel Treatments Interact With Climate and Biophysical Setting to Affect Fire, Water, and Forest Health: A Process-Based Modeling Approach. Front. For. Glob. Chang. 2021, 3, 143. [Google Scholar] [CrossRef]
  2. Westerling, A.L.; Hidalgo, H.G.; Cayan, D.R.; Swetnam, T.W. Warming and earlier spring increase western US forest wildfire activity. Science 2006, 313, 940–943. [Google Scholar] [CrossRef]
  3. Vose, J.M.; Peterson, D.L.; Domke, G.M.; Fettig, C.J.; Joyce, L.A.; Keane, R.E.; Luce, C.H.; Prestemon, J.P.; Band, L.E.; Clark, J.S.; et al. Forests. In Impacts, Risks, and Adaptation in the United States: Fourth National Climate Assessment, Volume II; Reidmiller, D.R., Avery, C.W., Easterling, D.R., Kunkel, K.E., Lewis, K.L.M., Maycock, T.K., Stewart, B.C., Eds.; U.S. Global Change Research Program: Washington, DC, USA, 2018; pp. 232–267. [Google Scholar] [CrossRef]
  4. Western Forestry Leadership Coalition. The True Cost of Wildfire in the Western U.S. 2010. Available online: https://www.blm.gov/or/districts/roseburg/plans/collab_forestry/files/TrueCostOfWilfire.pdf (accessed on 1 January 2020).
  5. NOAA National Centers for Environmental Information. 2022. Available online: https://www.ncei.noaa.gov/has/HAS.FileAppRouter?datasetname=9950_01&subqueryby=STATION&applname=&outdest=FILE (accessed on 1 January 2020).
  6. Abatzoglou, J.T.; Williams, A.P. Impact of anthropogenic climate change on wildfire across western US forests. Proc. Natl. Acad. Sci. USA 2016, 113, 11770–11775. [Google Scholar] [CrossRef] [PubMed]
  7. Radeloff, V.C.; Helmers, D.P.; Kramer, H.A.; Mockrin, M.H.; Alexandre, P.M.; Bar-Massada, A.; Butsic, V.; Hawbaker, T.J.; Syphard, A.D.; Stewart, S.I. Rapid Growth of US Wildland-Urban Interface Raises Wildfire Risk. Proc. Natl. Acad. Sci. USA 2018, 115, 3314–3319. [Google Scholar] [CrossRef]
  8. Caggiano, M.D.; Hawbaker, T.J.; Gannon, B.M.; Hoffman, C.M. Building loss in wui disasters: Evaluating the core components of the wildland–urban interface definition. Fire 2020, 3, 73. [Google Scholar] [CrossRef]
  9. Blach, J.K.; Bradley, B.A.; Abatzoglou, J.T.; Nagy, R.C.; Fusco, E.J.; Mahood, A.L. Human-started wildfires expand the fire niche across the United States. Proc. Natl. Acad. Sci. USA 2017, 114, 2946–2951. [Google Scholar] [CrossRef] [PubMed]
  10. Scott, J.H.; Gilbertson-Day, J.W.; Moran, C.; Dillon, G.K.; Short, K.C.; Vogler, K.C. Wildfire Risk to Communities: Spatial Datasets of Landscape-Wide Wildfire Risk Components for the United States; Forest Service Research Data Archive: Fort Collins, CO, USA, 2020. [CrossRef]
  11. LANDFIRE. Existing Vegetation Type Layer, LANDFIRE 2.0.0, U.S. Department of the Interior, Geological Survey, and U.S. Department of Agriculture. 2020. Available online: http://landfire.cr.usgs.gov/viewer/ (accessed on 1 July 2021).
  12. Oliveira, S.; Rocha, J.; Sá, A. Wildfire risk modeling. Curr. Opin. Environ. Sci. Health 2021, 23, 100274. [Google Scholar] [CrossRef]
  13. Smith, A.B. 2021 US Billion Dollar Weather and Climate Disasters in Historical Context including New County-Level Exposure, Vulnerability and Projected Damage Mapping. In Proceedings of the 102nd American Meteorological Society Annual Meeting, Houston, TX, USA, 23–27 January 2022; AMS: Providence, RI, USA, 2022. [Google Scholar]
  14. Meldrum, J.R.; Barth, C.M.; Goolsby, J.B.; Olson, S.K.; Gosey, A.C.; White, J.; Brenkert-Smith, H.; Champ, P.A.; Gomez, J. Parcel-Level Risk Affects Wildfire Outcomes: Insights from Pre-Fire Rapid Assessment Data for Homes Destroyed in 2020 East Troublesome Fire. Fire 2022, 5, 24. [Google Scholar] [CrossRef]
  15. Helmbrecht, D.; Gilbertson-Day, J.; Scott, J.H.; Hollingsworth, L. Wildfire Risk to Residential Structures in the Island Park Sustainable Fire Community: Caribou-Targhee National Forest; US Department of Agriculture, Forest Service, Fire Sciences Lab.: Missoula, MT, USA, 2016; 34p.
  16. Glickman, D.; Babbitt, B. Urban Wildland Interface Communities within the Vicinity of Federal Lands That Are at High Risk from Wildfire. Fed. Regist. 2001, 66, 751–777. [Google Scholar]
  17. Iglesias, V.; Stavros, N.; Balch, J.K.; Barrett, K.; Cobian-Iñiguez, J.; Hester, C.; Kolden, C.A.; Leyk, S.; Nagy, R.C.; Reid, C.E.; et al. Fires that matter: Reconceptualizing fire risk to include interactions between humans and the natural environment. Environ. Res. Lett. 2022, 17, 045014. [Google Scholar] [CrossRef]
  18. Scott, J.H.; Short, K.C.; Finney, M.A. “FSim: The Large Fire Simulator Guide to Best Practices”, Pyrologix LLC. Available online: https://pyrologix.com/wp-content/uploads/2019/11/FSimBestPractices_0.3.1.pdf (accessed on 5 March 2018).
  19. Ott, C.W.; Adhikari, B.; Alexander, S.P.; Hodza, P.; Xu, C.; Minckley, T.A. Predicting fire propagation across heterogeneous landscapes using WyoFire: A Monte Carlo-driven wildfire model. Fire 2020, 3, 71. [Google Scholar] [CrossRef]
  20. Lautenberger, C. Mapping Areas at Elevated Risk of Large-Scale Structure Loss Using Monte Carlo Simulation and Wildland Fire Modeling. Fire Saf. J. 2017, 91, 768–775. [Google Scholar] [CrossRef]
  21. LFD: Landfire Fuel Dictionary. Available online: https://landfire.gov/DataDictionary/LF200/LF2016Remap_CC.pdf (accessed on 1 January 2020).
  22. Finco, M.; Quayle, B.; Zhang, Y.; Lecker, J.; Megown, K.A.; Brewer, C.K. Monitoring Trends and Burn Severity (MTBS): Monitoring wildfire activity for the past quarter century using landsat data. In Proceedings of the Moving from Status to Trends: Forest Inventory and Analysis (FIA) Symposium 2012, Baltimore, MD, USA, 4–6 December 2012; Gen. Tech. Rep. NRS-P-105. U.S. Department of Agriculture, Forest Service, Northern Research Station: Newtown Square, PA, USA, 2012; pp. 222–228. [Google Scholar]
  23. Miller, J.D.; Thode, A.E. Quantifying burn severity in a heterogeneous landscape with a relative version of the delta Normalized Burn Ratio (dNBR). Remote Sens. Environ. 2007, 109, 66–80. [Google Scholar] [CrossRef]
  24. Hansen, M.C.; Potapov, P.V.; Moore, R.; Hancher, M.; Turubanova, S.A.; Tyukavina, A.; Thaus, D.; Stehmans, S.V.; Goetzt, J.S.; Loveland, T.R.; et al. High-Resolution Global Maps of 21st-Century Forest Cover Change. Science 2013, 342, 850–853. [Google Scholar] [CrossRef] [PubMed]
  25. Scott, J.H.; Burgan, R.E. Standard Fire Behavior Fuel Models: A Comprehensive Set for Use with Rothermel’s Surface Fire Spread Model; Gen. Tech. Rep. RMRS-GTR-153; U.S. Department of Agriculture, Forest Service, Rocky Mountain Research Station: Fort Collins, CO, USA, 2005; 72p.
  26. Reeves, M.C.; Ryan, K.C.; Rollins, M.C.; Thompson, T.G. Spatial fuel data products of the LANDFIRE project. Int. J. Wildland Fire 2009, 18, 250–267. [Google Scholar] [CrossRef]
  27. Martinuzzi, S.; Stewart, S.I.; Helmers, D.P.; Mockrin, M.H.; Hammer, R.B.; Radeloff, V.C. The 2010 Wildland-Urban Interface of the Conterminous United States; Research Map NRS-8; U.S. Department of Agriculture, Forest Service, Northern Research Station: Newtown Square, PA, USA, 2015; 124p. [CrossRef]
  28. Krawchuk, M.A.; Moritz, M.A.; Parisien, M.A.; Dorn, J.V.; Hayhoe, K. Global Pyrogeography: The Current and Future Distribution of Wildfire. PLoS ONE 2009, 4, e5102. [Google Scholar] [CrossRef] [PubMed]
  29. Bowman, D.M.J.S.; Kolden, C.A.; Abatzoglou, J.T.; Johnston, F.H.; van der Werf, G.R.; Flannigan, M. Vegetation fires in the Anthropocene. Nat. Rev. Earth Environ. 2020, 1, 500–515. [Google Scholar] [CrossRef]
  30. LUCAS: Land Use and Carbon Scenario Simulator. Available online: https://www.usgs.gov/centers/western-geographic-science-center/science/lucas-model (accessed on 1 January 2020).
  31. PRISM: Parameter-elevation Relationships on Independent Slopes Model. Available online: https://prism.oregonstate.edu/ (accessed on 1 January 2020).
  32. Taylor, K.E.; Stouffer, R.J.; Meehl, G.A. An Overview of CMIP5 and the experiment design. Bull. Amer. Meteor. Soc. 2012, 93, 485–498. [Google Scholar] [CrossRef]
  33. Abatzoglou, J.T.; Brown, T.J. A comparison of statistical downscaling methods suited for wildfire applications. Int. J. Climatol. 2012, 32, 772–780. [Google Scholar] [CrossRef]
  34. Bates, P.D.; Niall, Q.; Chris, S.; Andrew, S.; Oliver, W.; Jeison, S.; James, S.; Gaia, O.; Jeff, N.; Guy, S.; et al. Combined Modeling of US Fluvial, Pluvial, and Coastal Flood Hazard Under Current and Future Climates. Water Resour. Res. 2021, 57, e2020WR028673. [Google Scholar] [CrossRef]
  35. Kearns, E.J.; Amodeo, M.; Chadwick, S.; Eby, M.; Porter, J.R. Making Climate Change Personal: Enabling Action Through Communication of Flood Risk Science to Individuals. In Proceedings of the AGU Fall Meeting Abstracts, Online, 1–17 December 2020; Volume 2020, p. SY040-06. [Google Scholar]
  36. Cohen, J.D. Reducing the wildland fire threat to homes: Where and how much? In Proceedings of the Symposium on Fire Economics, Planning, and Policy: Bottom Lines, San Diego, CA, USA, 5–9 April 1999; Gen. Tech. Rep. PSW-GTR-173. U.S. Department of Agriculture, Forest Service, Pacific Southwest Research Station: Albany, CA, USA, 1999; pp. 189–195. [Google Scholar]
  37. Short, K.C. A spatial database of wildfires in the United States, 1992–2011. Earth Syst. Sci. Data 2014, 6, 1–27. [Google Scholar] [CrossRef]
  38. Short, K.C. Spatial Wildfire Occurrence Data for the United States, 1992–2018. [FPA_FOD_20210617], 5th ed.; Forest Service Research Data Archive: Fort Collins, CO, USA, 2021. [CrossRef]
  39. Dillon, G.K.; Menakis, J.; Fay, F. Wildland Fire Potential: A Tool for Assessing Wildfire Risk and Fuels Management Needs. In Proceedings of the Large Wildland Fires Conference, Missoula, MT, USA, 19–23 May 2014; Proc. RMRS-P-73. Keane, R.E., Jolly, M., Parsons, R., Riley, K., Eds.; U.S. Department of Agriculture, Forest Service, Rocky Mountain Research Station: Fort Collins, CO, USA, 2015; pp. 60–76. [Google Scholar]
  40. Silverman, B.W. Density Estimation for Statistics and Data Analysis; Chapman and Hall: New York, NY, USA, 1986. [Google Scholar]
  41. Riley, K.L.; Abatzoglou, J.T.; Grenfell, I.C.; Klene, A.E.; Heinsch, F.A. The relationship of large fire occurrence with drought and fire danger indices in the western USA, 1984–2008: The role of temporal scale. Int. J. Wildland Fire 2013, 22, 894–909. [Google Scholar] [CrossRef]
  42. Rehm, R.G.; McDermott, R.J. Fire-Front Propagation Using the Level Set Method; US Department of Commerce, National Institute of Standards and Technology: Gaithersburg, MD, USA, 2009; pp. 1–12.
  43. Lautenberger, C. Wildland Fire Modeling with an Eulerian Level Set Method and Automated Calibration. Fire Saf. J. 2013, 62, 289–298. [Google Scholar] [CrossRef]
  44. Finney, M.A.; McHugh, C.W.; Grenfell, I.C.; Riley, K.L.; Short, K.C. A simulation of probabilistic wildfire risk components for the continental United States. Stoch. Environ. Res. Risk Assess. 2011, 25, 973–1000. [Google Scholar] [CrossRef]
  45. Parisien, M.A.; Dawe, D.A.; Miller, C.; Stockdale, C.A.; Armitage, O.B. Applications of simulation-based burn probability modeling: A review. Int. J. Wildland Fire 2019, 28, 913–926. [Google Scholar] [CrossRef]
  46. Finney, M.A. FARSITE, Fire Area Simulator–Model Development and Evaluation, No. 4; US Department of Agriculture Forest Service, Rocky Mountain Research Station: Fort Collins, CO, USA, 1998.
  47. Rothermel, R.C. A Mathematical Model for Predicting Fire Spread in Wildland Fuels; Research Paper Int-115; USDA Forest Service: Fort Collins, CO, USA, 1972.
  48. Anderson, H.E. Predicting Wind-Driven Wild Land Fire Size and Shape; Research Paper INT-RP-305; United States Department of Agriculture Forest Service, Intermountain Forest and Range Experiment Station: Fort Collins, CO, USA, 1983. [Google Scholar]
  49. Richards, G.D. A General Mathematical Framework for Modelling Two-Dimensional Wildland Fire Spread. Int. J. Wildland Fire 1995, 5, 63–72. [Google Scholar] [CrossRef]
  50. Finney, M.A. Landscape fire simulation and fuel treatment optimization. In Methods for Integrating Modeling of Landscape Change: Interior Northwest Landscape Analysis System; US Department of Agriculture, Forest Service, Pacific Northwest Research Station: Portland, OR, USA, 2004; pp. 117–131. [Google Scholar]
  51. Van Wagner, C.E. Conditions for the Start and Spread of Crown Fire. Can. J. For. Res. 1977, 7, 23–34. [Google Scholar] [CrossRef]
  52. Cruz, M.G.; Alexander, M.E.; Wakimoto, R.H. Development and testing of models for predicting crown fire rate of spread in conifer forest stands. Can. J. For. Res. 2005, 35, 1626–1639. [Google Scholar] [CrossRef]
  53. Sardoy, N.; Consalvi, J.L.; Kaiss, A.; Fernandez-Pello, A.C.; Porterie, B. Numerical study of ground-level distribution of firebrands generated by line fires. Combust. Flame 2008, 154, 478–488. [Google Scholar] [CrossRef]
  54. Perryman, H.A.; Dugaw, C.J.; Varner, J.M.; Johnson, D.L. A cellular automata model to link surface fires to firebrand lift-off and dispersal. Int. J. Wildland Fire 2013, 22, 428–439. [Google Scholar] [CrossRef]
  55. Sethian, J.A. A fast marching level set method for monotonically advancing fronts. Proc. Natl. Acad. Sci. USA 1996, 93, 1591–1595. [Google Scholar] [CrossRef]
  56. First Street Foundation. The First National Flood Risk Assessment: Defining America’s Growing Risk. 2020; Available online: https://assets.firststreet.org/uploads/2020/06/first_street_foundation__first_national_flood_risk_assessment.pdf (accessed on 8 August 2022).
  57. Iowa State, Iowa Environmental Mesonet. 2022. Available online: https://mtarchive.geol.iastate.edu/ (accessed on 1 June 2022).
  58. Torralba, V.; Doblas-Reyes, F.J.; Gonzalez-Reviriego, N. Uncertainty in recent near-surface wind speed trends: A global reanalysis intercomparison. Environ. Res. Lett. 2017, 12, 114019. [Google Scholar] [CrossRef]
  59. Fang, G.H.; Yang, J.; Chen, Y.N.; Zammit, C. Comparing bias correction methods in downscaling meteorological variables for a hydrologic impact study in an arid area in China. Hydrol. Earth Syst. Sci. 2015, 19, 2547–2559. [Google Scholar] [CrossRef]
  60. Luo, M.; Liu, T.; Meng, F.; Duan, Y.; Frankl, A.; Bao, A.; De Maeyer, P. Comparing bias correction methods used in downscaling precipitation and temperature from regional climate models: A case study from the Kaidu River Basin in Western China. Water 2018, 10, 1046. [Google Scholar] [CrossRef]
  61. Switanek, M.B.; Troch, P.A.; Castro, C.L.; Leuprecht, A.; Chang, H.; Mukherjee, R.; Demaria, E.M.C. Scaled distribution mapping: A bias correction method that preserves raw climate model projected changes. Hydrol. Earth Syst. Sci. 2017, 21, 2649–2666. [Google Scholar] [CrossRef]
  62. Gudmundsson, L.; Bremnes, J.B.; Haugen, J.E.; Engen-Skaugen, T. Technical Note: Downscaling RCM precipitation to the station scale using statistical transformations—A comparison of methods. Hydrol. Earth Syst. Sci. 2012, 16, 3383–3390. [Google Scholar] [CrossRef]
Figure 1. The modifications of the LANDFIRE v2.0.0 fuels (LF), including use of disturbance data (FVC—fuel vegetation cover; FVH—fuel vegetation height; FVT—fuel vegetation type; BPS—biophysical settings) to modify surface (FM40) and canopy fuel classes (canopy cover—CC; canopy height—CH; canopy base height—CGH; canopy bulk density—CBD), and the replacement of non-burnable fuel classes in the WUI with fuel classes that were found to approximate observed fire behavior in past WUI fires. Computation was performed using Google Earth Engine (EE). This figure represents steps 1 and 2 in the fuel methodology; steps 3 and 4 are included in Figure 4.
Figure 1. The modifications of the LANDFIRE v2.0.0 fuels (LF), including use of disturbance data (FVC—fuel vegetation cover; FVH—fuel vegetation height; FVT—fuel vegetation type; BPS—biophysical settings) to modify surface (FM40) and canopy fuel classes (canopy cover—CC; canopy height—CH; canopy base height—CGH; canopy bulk density—CBD), and the replacement of non-burnable fuel classes in the WUI with fuel classes that were found to approximate observed fire behavior in past WUI fires. Computation was performed using Google Earth Engine (EE). This figure represents steps 1 and 2 in the fuel methodology; steps 3 and 4 are included in Figure 4.
Fire 05 00117 g001
Figure 2. Disturbances: locations of updated canopy and surface fuels.
Figure 2. Disturbances: locations of updated canopy and surface fuels.
Fire 05 00117 g002
Figure 3. Modified surface fuel estimates: (left) LANDFIRE 2016 canopy cover (green) and identified disturbances 2017–2020 (gray) are shown. (Right) new canopy cover estimates (green) are reduced in disturbed areas.
Figure 3. Modified surface fuel estimates: (left) LANDFIRE 2016 canopy cover (green) and identified disturbances 2017–2020 (gray) are shown. (Right) new canopy cover estimates (green) are reduced in disturbed areas.
Fire 05 00117 g003
Figure 4. Random forest model training and classification to replace non-burnable FM40 classes with burnable classes in the newly defined WUI.
Figure 4. Random forest model training and classification to replace non-burnable FM40 classes with burnable classes in the newly defined WUI.
Fire 05 00117 g004aFire 05 00117 g004b
Figure 5. WUI FM40 random forest importance values.
Figure 5. WUI FM40 random forest importance values.
Fire 05 00117 g005
Figure 6. Classification of building structures from unburnable to burnable fuel types.
Figure 6. Classification of building structures from unburnable to burnable fuel types.
Fire 05 00117 g006
Figure 7. Number of fires in the Western US as a function of ERC(G) percentile (reproduced from Riley et al., 2013 [41]).
Figure 7. Number of fires in the Western US as a function of ERC(G) percentile (reproduced from Riley et al., 2013 [41]).
Fire 05 00117 g007
Figure 8. Wildfire behavior model, process diagram. The “fire factor” result on the right hand side is an expression of wildfire hazard as a function of burn probability, flame length, and ember presence and is available for viewing for every property in the U.S. at riskfactor.com (accessed on 8 August 2022).
Figure 8. Wildfire behavior model, process diagram. The “fire factor” result on the right hand side is an expression of wildfire hazard as a function of burn probability, flame length, and ember presence and is available for viewing for every property in the U.S. at riskfactor.com (accessed on 8 August 2022).
Fire 05 00117 g008
Figure 9. Example ELMFIRE fire spread simulation for individual fire ignition. (a) Fire type (surface fire, passive crown fire, or active crown fire). (b) Flame length in feet.
Figure 9. Example ELMFIRE fire spread simulation for individual fire ignition. (a) Fire type (surface fire, passive crown fire, or active crown fire). (b) Flame length in feet.
Fire 05 00117 g009
Figure 10. Historical and modeled fire sizes (in acres) versus intensity (flame length in ft).
Figure 10. Historical and modeled fire sizes (in acres) versus intensity (flame length in ft).
Fire 05 00117 g010
Figure 11. From riskfactor.com, (accessed on 8 August 2022): (bottom) a representative description of property-level exposure (in Naples, FL), showing the high-resolution, property-specific nature of the estimates that are produced by the FSF-WFM. Colors depict the probability (%) at 30 m resolution of being impacted by wildfire hazard during the year 2022. The likelihood numbers at right are estimates of annual likelihood of wildfire for the property at the red pin (center), and which is predicted to more than double by 2052. (Top) CONUS burn probability at 30 m resolution for 2022; gray areas show areas with negligible exposure.
Figure 11. From riskfactor.com, (accessed on 8 August 2022): (bottom) a representative description of property-level exposure (in Naples, FL), showing the high-resolution, property-specific nature of the estimates that are produced by the FSF-WFM. Colors depict the probability (%) at 30 m resolution of being impacted by wildfire hazard during the year 2022. The likelihood numbers at right are estimates of annual likelihood of wildfire for the property at the red pin (center), and which is predicted to more than double by 2052. (Top) CONUS burn probability at 30 m resolution for 2022; gray areas show areas with negligible exposure.
Fire 05 00117 g011
Figure 12. (A) Geographic distribution of “any exposure” to wildfire (>1% cumulative exposure): those individual properties with >1% cumulative exposure over 2022–2052 were counted on a state by state basis, and compared to the total number of all properties. (B). Geographic distribution of “major exposure” to wildfire (>3% cumulative exposure): those individual properties with >3% cumulative exposure over 2022–2052 were counted on a state by state basis, and compared to the total number of all properties.
Figure 12. (A) Geographic distribution of “any exposure” to wildfire (>1% cumulative exposure): those individual properties with >1% cumulative exposure over 2022–2052 were counted on a state by state basis, and compared to the total number of all properties. (B). Geographic distribution of “major exposure” to wildfire (>3% cumulative exposure): those individual properties with >3% cumulative exposure over 2022–2052 were counted on a state by state basis, and compared to the total number of all properties.
Fire 05 00117 g012aFire 05 00117 g012b
Figure 13. Percent Increase in annual likelihood of wildfire exposure among properties with at least 0.03% annual likelihood of wildfire exposure.
Figure 13. Percent Increase in annual likelihood of wildfire exposure among properties with at least 0.03% annual likelihood of wildfire exposure.
Fire 05 00117 g013
Table 1. Distribution of disturbance types and severity.
Table 1. Distribution of disturbance types and severity.
Disturbance Type%
Fire92.8%
Mechanical add0.2%
Mechanical reduce5.8%
Other1.2%
Disturbance Severity%
Low65.0%
Medium20.2%
High14.8%
Table 2. Hazards quantified by the FSF-WFM: following the ELMFIRE simulations for 2022 and 2052, the wildfire hazard for any 30 m pixel within the CONUS domain is represented by the combination of the burn probability, mean or maximum flame length (intensity), and exposure to embers.
Table 2. Hazards quantified by the FSF-WFM: following the ELMFIRE simulations for 2022 and 2052, the wildfire hazard for any 30 m pixel within the CONUS domain is represented by the combination of the burn probability, mean or maximum flame length (intensity), and exposure to embers.
Measure of ExposureDescriptionUnits
Burn probabilityLikelihood that a pixel catches fire out of all the simulations normalized by likelihood of ignitions%
Max flame lengthMaximum flame length experienced at a pixel across all simulationsft
Sum of flame lengthSum of flame length for all simulations that experience fireft
Mean flame lengthSum of flame length divided by times burnedft
Binned counts of flame lengths(0,2,4,6,8,12,20,+)ft
Ember luxDimensionless number that is a proxy for the count of embers landing in a pixel. Does not reflect mass of embers, whether they are still burning, or distance traveled<none>
Ember likelihoodLikelihood that an ember falls into a pixel across all simulations, similar to times burned%
Max embersMax number of emberscount
Table 3. Top 25 state ranking by “any risk” (cumulative burn probability of >1%).
Table 3. Top 25 state ranking by “any risk” (cumulative burn probability of >1%).
StateTotal PropertiesAny RiskPct Any Risk
Texas11,957,7079,450,09179.03
Florida8,975,2807,197,68580.19
California11,341,3837,131,84962.88
North Carolina5,451,2783,126,13057.35
Alabama3,019,3002,727,45590.33
Georgia4,413,8392,482,09156.23
Arizona3,225,7632,463,01976.35
Virginia3,795,4182,265,92759.70
South Carolina2,616,0912,068,04879.05
Colorado2,491,6102,000,32180.28
Oklahoma2,215,7551,901,85085.83
Tennessee3,278,7391,879,31657.32
New Jersey3,449,5411,859,39553.90
Mississippi1,904,4941,695,46289.02
Arkansas1,923,5561,558,00581.00
Missouri3,191,5021,503,14347.10
Minnesota2,964,7081,472,20649.66
New Mexico1,495,3921,380,73692.33
Louisiana2,365,2071,254,93653.06
Utah1,363,4631,153,35684.59
Kansas1,633,5211,087,98866.60
New York5,376,613999,21718.58
Washington3,031,769996,96032.88
Oregon1,807,336911,74550.45
Idaho1,036,925878,06884.68
Table 4. Top 25 state ranking by “major risk” (cumulative burn probability of >3%).
Table 4. Top 25 state ranking by “major risk” (cumulative burn probability of >3%).
StateTotal PropertiesMajor RiskPct Major Risk
California11,341,3832,554,77722.53
Texas11,957,7071,686,57114.10
Florida8,975,2801,540,41317.16
Arizona3,225,763998,24130.95
Oklahoma2,215,755451,92820.40
Utah1,363,463425,16331.18
New Mexico1,495,392409,53827.39
Nevada1,209,308314,20325.98
Idaho1,036,925196,01418.90
Washington3,031,769187,2756.18
Colorado2,491,610177,0817.11
New Jersey3,449,541171,5684.97
Montana894,052167,04018.68
South Dakota666,388164,70224.72
Mississippi1,904,494121,3676.37
Wyoming339,209113,57033.48
North Carolina5,451,27896,7741.78
Kansas1,633,52181,3094.98
Oregon1,807,33670,6803.91
Alabama3,019,30053,7261.78
Nebraska1,138,19150,1504.41
South Carolina2,616,09146,2921.77
North Dakota679,02333,7374.97
Louisiana2,365,20723,7831.01
Minnesota2,964,70819,9290.67
Table 5. Fire factor assignment.
Table 5. Fire factor assignment.
Fire FactorCriteria
1No modeled exposure of being in a wildfire (burn probability) and no modeled exposure to embers (including not being in an ember zone) are considered to have minimal risk.
2Located in areas exposed to embers through the created ember zone or in an area with less than a 1% cumulative chance of burning over 30 years.
3With 1–3% chance of burning over 30 years.
4With a 4–6% chance of burning over 30 years.
5With a 7–9% chance of burning over 30 years.
6With a 10–14% chance of burning over 30 years.
7With a 14–20% chance of burning over 30 years.
8With a 21–26% chance of burning over 30 years.
9With a 27–36% chance of burning over 30 years.
10With more than a 36% chance of burning over 30 years.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kearns, E.J.; Saah, D.; Levine, C.R.; Lautenberger, C.; Doherty, O.M.; Porter, J.R.; Amodeo, M.; Rudeen, C.; Woodward, K.D.; Johnson, G.W.; et al. The Construction of Probabilistic Wildfire Risk Estimates for Individual Real Estate Parcels for the Contiguous United States. Fire 2022, 5, 117. https://doi.org/10.3390/fire5040117

AMA Style

Kearns EJ, Saah D, Levine CR, Lautenberger C, Doherty OM, Porter JR, Amodeo M, Rudeen C, Woodward KD, Johnson GW, et al. The Construction of Probabilistic Wildfire Risk Estimates for Individual Real Estate Parcels for the Contiguous United States. Fire. 2022; 5(4):117. https://doi.org/10.3390/fire5040117

Chicago/Turabian Style

Kearns, Edward J., David Saah, Carrie R. Levine, Chris Lautenberger, Owen M. Doherty, Jeremy R. Porter, Michael Amodeo, Carl Rudeen, Kyle D. Woodward, Gary W. Johnson, and et al. 2022. "The Construction of Probabilistic Wildfire Risk Estimates for Individual Real Estate Parcels for the Contiguous United States" Fire 5, no. 4: 117. https://doi.org/10.3390/fire5040117

APA Style

Kearns, E. J., Saah, D., Levine, C. R., Lautenberger, C., Doherty, O. M., Porter, J. R., Amodeo, M., Rudeen, C., Woodward, K. D., Johnson, G. W., Markert, K., Shu, E., Freeman, N., Bauer, M., Lai, K., Hsieh, H., Wilson, B., McClenny, B., McMahon, A., & Chishtie, F. (2022). The Construction of Probabilistic Wildfire Risk Estimates for Individual Real Estate Parcels for the Contiguous United States. Fire, 5(4), 117. https://doi.org/10.3390/fire5040117

Article Metrics

Back to TopTop