Direct and indirect climate controls predict heterogeneous early-mid 21st century wildfire burned area across western and boreal North America

Predicting wildfire under future conditions is complicated by complex interrelated drivers operating across large spatial scales. Annual area burned (AAB) is a useful index of global wildfire activity. Current and antecedent seasonal climatic conditions, and the timing of snowpack melt, have been suggested as important drivers of AAB. As climate warms, seasonal climate and snowpack co-vary in intricate ways, influencing fire at continental and sub-continental scales. We used independent records of seasonal climate and snow cover duration (last date of permanent snowpack, LDPS) and cell-based Structural Equation Models (SEM) to separate direct (climatic) and indirect (snow cover) effects on relative changes in AAB under future climatic scenarios across western and boreal North America. To isolate seasonal climate variables with the greatest effect on AAB, we ran multiple regression models of log-transformed AAB on seasonal climate variables and LDPS. We used the results of multiple regressions to project future AAB using GCM ensemble climate variables and LDPS, and validated model predictions with recent AAB trends. Direct influences of spring and winter temperatures on AAB are larger and more widespread than the indirect effect mediated by changes in LDPS in most areas. Despite significant warming trends and reductions in snow cover duration, projected responses of AAB to early-mid 21^st century are heterogeneous across the continent. Changes in AAB range from strongly increasing (one order of magnitude increases in AAB) to moderately decreasing (more than halving of baseline AAB). Annual wildfire area burned in coming decades is likely to be highly geographically heterogeneous, reflecting interacting regional and seasonal climate drivers of fire occurrence and spread