Landscape level analysis of mountain pine beetle in British Columbia, Canada: spatiotemporal development and spatial synchrony within the present outbreak

Eruptive herbivores can exert profound landscape level influences. For example, the ongoing mountain pine beetle outbreak in British Columbia, Canada, has resulted in mortality of mature lodgepole pine over > 7 million ha. Analysis of the spatio-temporal pattern of spread can lend insights into the processes initiating and/or sustaining such phenomena. We present a landscape level analysis of the development of the current outbreak. Aerial survey assessments of tree mortality, projected onto discrete 12x12 km cells, were used as a proxy for insect population density. We examined whether the outbreak potentially originated from an epicenter and spread, or whether multiple localized populations erupted simultaneously at spatially disjunct locations. An aspatial cluster analysis of time series from 1990 to 2003 revealed four distinct time series patterns. Each time series demonstrated a general progression of increasing mountain pine beetle populations. Plotting the geographical locations of each temporal pattern revealed that the outbreak occurred first in an area of west-central British Columbia, and then in an area to the east. The plot further revealed many localized infestations erupted in geographically disjunct areas, especially in the southern portion of the province. Autologistic regression analyses indicated a significant, positive association between areas where the outbreak first occurred and conservation lands. For example, the delineated area of west-central British Columbia is comprised of three conservation parks and adjacent working forest. We further examined how population synchrony declines with distance at different population levels. Examination of the spatial dependence of temporal synchrony in population fluctuations during early, incipient years (i.e. 1990-1996) suggested that outbreaking mountain pine beetle populations are largely independent at scales > 200 km during non-epidemic periods. However, during epidemic years (i.e. 1999-2003), populations were clearly synchronous across the entire province, even at distances of up to 900 km. The epicentral pattern of population development can be used to identify and prioritize adjacent landscape units for both reactive and proactive management strategies intended to minimize mountain pine beetle impacts.