- 1Department of Mathematics and Statistics, University of Strathclyde, Glasgow, United Kingdom
- 2Biomathematics and Statistics Scotland, Edinburgh, United Kingdom
- 3Plant Health Centre Scotland, James Hutton Institute, Dundee, United Kingdom
- 4Department of Environment and Geography, University of York, York, United Kingdom
Introduction: Wood borer pests represent a growing threat to forestry worldwide, with potential for an increase in climate change induced outbreaks that could significantly impact ecosystem functioning.
Methods: This study reviews the current state-of-the-art in bark beetle modeling as presented in peer-reviewed English language scientific papers and reviews cited in Clarivate Analytics Web of Science (WoS) Core Collection in the period from 2006 to 2023, but earlier influential papers cited in reviews from this period are also in scope. We categorize studies by modeling methodology, by host and by modeled processes.
Results: Ash, pine, and spruce hosts account for 88% of studies, the majority of which focus on the continental USA, British Columbia, Canada and Central Europe. In terms of methodology, statistical methods are the most commonly employed technique, and the majority of articles model just one or two key processes with pest demography and environmental factors being the most studied across different hosts. Commonalities in methods used across pest-host systems include: use of phenology-like models; application of static species distribution models (SDMs) to understand climate impacts; modeling of spread via local and long-range kernels; and use of economic cost-benefit analysis as a tool to guide management.
Discussion: We identify several gaps in current research including quantifying economic consequences of wood borer pests and the need for greater understanding of ecological impacts and resilience. Agent- and individual-based models may also provide useful tools for understanding the complexity of socio-ecological system dynamics. However, such developments should be in tandem with wider use of techniques for parameter estimation and uncertainty quantification, including Bayesian inference in particular of the dynamics of spatial spread. A related challenge is better quantification of pathways - including trade - for entry by invasive pests, coupled with a greater understanding of potential vulnerabilities of forest systems to environmental drivers like climate change, coupled with potentially multiple endemic, emerging and novel pests. Addressing these challenges would enable both better mitigation of risks associated with wood borer beetle infestations and better management of outbreaks when they do occur.
1 Introduction
Biological pest species pose threats to the stability and functioning of ecosystems worldwide. Non-native pests can have complex and significant impacts on invaded ecosystems by affecting a variety of interacting mechanisms, e.g carbon sequestration, nutrient cycling (Ehrenfeld, 2010). Stresses induced by climate change can result in large non-normative outbreaks of native species that devastate ecosystem functioning (Logan et al., 2003). A systematic review found the distribution of evidence for direct (29.9%), indirect (29.8%), and interacting effects (40.8%) of climate change related to forest pests (Seidl et al., 2017). Development in global trade and human activity are increasing rates of arrival and opening new routes for potentially damaging pests (Dalmazzone and Giaccaria, 2014).
Wood borer pests exemplify such issues, representing a growing threat to forestry worldwide (Linnakoski and Forbes, 2019). For example, since 2002, Emerald ash borer (EAB), Agrilus planipennis has emerged as a major invasive pest of Fraxinus species in the USA (Cappaert et al., 2005), where it has caused an estimated US$ 10.7–12.5 billion (bn) of economic damage in the period 2010–2020 (Kovacs et al., 2011). EAB is also present in Western Russia and is steadily spreading west toward continental Europe (Straw et al., 2013). European ash is suffering an ongoing outbreak of fungal pathogen Hymenoscyphus fraxineus—the causal agent of ash dieback. In the UK alone, damages from ash dieback are estimated at GBP£15bn (Hill et al., 2019), which illustrates the potential impact of EAB introduction even without considering the increased impact of introduction into H. fraxinus infested populations (Musolin et al., 2017). EAB shows how invasive species not considered pests within the native range can cause major damage in invaded environments.
In contrast, an endemic forest pest which experiences large outbreaks under favorable climatic conditions is the European spruce beetle (ESB), Ips typographus. Endemic to central Europe, ESB outbreaks are responsible for the mortality of an estimated 150 million m3 of Norway spruce from 1950 to 2000 (Schelhaas et al., 2003; Økland et al., 2016, 2012). Outbreaks occur after disturbance events such as storms or drought (Seidl and Rammer, 2017), and beetle induced mortality can alter forest hydrology (Beudert et al., 2018), carbon dynamics (Seidl et al., 2014), and food webs (Seidl et al., 2017). Climate change is predicted to increase the likelihood of disturbance events, and ESB outbreaks are likely to become increasingly devastating, e.g. through greater voltinism, the number of generations produced in a year (Wermelinger et al., 2012).
Mathematical and statistical modeling can provide critical insights and predictions that will aid the management of such risks. Modeling has provided insights into our understanding of outbreak dynamics and biological invasions (Logan et al., 2003; Sakai et al., 2001; Murray, 2001), been used to predict outbreak cycles (Kunegel-Lion and Lewis, 2020; Wildemeersch et al., 2019), determine spread rates and how to slow them (Wildemeersch et al., 2019; Sharov and Liebhold, 1998), and to detect tipping-point behavior (Vindstad et al., 2019). Such tools are essential for understanding and predicting destructive outbreaks in a climatically changing, increasingly international world (Ramsfield et al., 2016). However, despite the potential for lessons to be learned across pest-host systems, the literature on modeling of wood borer pests is fractured with a widely dispersed literature and a diverse group of disciplines. As far as we know, no attempts have been made to review modeling within wood borer pests and their societal impacts to date.
Therefore, here we conduct a systematic review that synthesizes relevant advances and techniques used to model key processes for wood borer pests. Studies reviewed are categorized attending to the identified themes of: pest spread; pest demography; environmental factors affecting infestation; and the economics and management of infestation. We synthesized studies based on the modeling techniques used to investigate the dynamics of wood borer bettle infestations, which are grouped into four high-level categories: mathematical models; statistical models; simulation studies; and Geographic Information Systems (GIS) studies (including use of imagery). We also examine the commonalities and differences in the use of different techniques between hosts, and the data requirements. Finally, we identify gaps for future research.
2 Methods
We conduct a systematic search of the scientific literature for articles related to the modeling of wood-boring species infestations and invasions. We followed established guidelines (Pullin et al., 2022) to conduct the systematic review, and our methods are outlined below.
2.1 Search terms
In order to construct a “search-string” to be used when searching for the most relevant studies in the bibliographic databases, search terms were developed based on a previous review of modeling approaches within animal health for which there is significant overlap (MacCalman et al., 2016). We developed a search string with multiple parts to capture the fact that we wanted to search for modeling studies of wood-boring species infestation risk and impacts. These search terms are grouped as follows: (1) terms that describe the pest of interest; (2) terms that capture the ecosystem of interest, forests; (3) terms that capture the focus on models; (4) and terms that cover the risk and economic impact of infestations. We also included parts in the search string related to the invasion process: introduction, establishment, and spread. Other terms used cover climatic effects. The full list of search terms and final search strings used is shown in the Supplementary material.
2.2 Search strategy and screening processes
Search terms are applied to the Clarivate Analytics Web of Science (WoS) Core Collection. We exclude other databases such as Scopus, ScienceDirect (which are operated by journal publishing houses) and Google Scholar because of the absence of features required for a rigorous, replicable search (Gusenbauer and Haddaway, 2020). The WoS database uses Boolean logic to group terms during searches, which we used to connect our search themes. As described in detail in the Supplementary material, the final search was given as “TS=(bark beetles AND forest) AND (models ORrisk) AND (introduction OR establishment OR spread) AND (climatic OR economic)”. In this search, we exclude pre-2006 articles and include only primary research articles, conference proceedings and review papers. The latter were only extracted to be able to include in our database relevant articles pre-2006. Therefore, the exclusion of pre-2006 papers was only partial in that pre-2006 references in post-2006 reviews were included in the database for their screening. This process has the effect of focusing on the more influential pre-2006 papers. When the original review was conducted, AI translation tools were less reliable than they have become, and therefore we excluded non-English language papers. The search for literature was conducted in two stages with the first search conducted on 28th March 2019 with papers from 2006 to 2019 which returned 1,655 articles for inspection, and the second search conducted on 19th September 2024 with papers from 2019 to 2023 which returned 1,068 articles for inspection. This second search was conducted to include articles published after 2019, ensuring the systematic review is up to date, after a change of lead analyst partway through the project. Any duplicates were then removed, and the resulting reference list was placed into a prepared spreadsheet for screening.
The literature sources were screened by title/abstract using a set of eligibility criteria (Table 1). Papers were excluded at this stage if the screening showed that: (1) an abstract was unavailable; (2) the focal study subjects was not wood borer beetles; (3) modeling aspects were not related to wood borer beetles; and (4) the study was solely empirical or observational [see Supplementary material]. The literature that met this eligibility criteria were then filtered into a list for inclusion in the full-text screening phase. Papers were selected if including modeling of wood borer processes: pest spread (including dispersal), pest demography, environmental factors, and economic considerations and management. Papers can use one or more of the modeling techniques to investigate these processes, typically categorized as: statistical, simulation and agent-based model (ABM, also known as individual-based model, IBM), mathematical, GIS-type (including imagery). This categorization reveals patterns across techniques, species, and processes (see Section 3).
The total articles accepted for review was 468 (Figure 1). A full list of the article titles and metadata, including abstracts, process categorisations, specific techniques used, article authors and purpose of the studies, can be found in the Supplementary material.
We conducted a quality assurance exercise to check the robustness of the screening process using a random sample of 160 articles from the post-2006 search. This corroboration exercise tested the screening of Abstracts (at both levels 1 and 2, as explained in Table 1). This allowed calculation of a κ-coefficient of agreement between the screening process (Cohen, 1960) of the full dataset and an independent screening of the sub-sample. The κ-coefficient is commonly used during the study selection process to demonstrate the consistency of decisions between two independent reviewers when screening titles, abstracts, or full texts for relevance, and compares the expected and observed level of agreement. These corroboration efforts produced a κ-coefficient of 0.687 and 0.637 for Abstract levels 1 and 2, respectively. These agreement coefficients lie within the range 0.64–0.69 deemed acceptable in the literature for a systematic review (McHugh, 2012).
3 Modeling of wood borer dynamics from a host perspective
Figures 2–4 show how species, processes, and techniques are represented within the literature. These findings indicate that most of the articles reviewed focus on pine and spruce as host species (39% & 34% of all articles, respectively, see Figure 1. We found that the majority of articles (52%) modeled two key processes, while 26% modeled one key process and 20% modeled 3 processes. Only 3% of all articles modeled four processes (see Figure 3A). The majority of articles used one or two modeling techniques, accounting for 92% of the total; whereas 7% used 3 techniques and <1% used four techniques (see Figure 3B). We found within each of the key processes modeled, demographic and environmental factors were generally the most numerous (see Figure 4A). We also found that statistical methods were generally the most commonly used method, both for modeling alongside other methods, but also as the sole method used, comprising over 35% of analyses (see Figure 4B).
Figure 2. Proportion of total articles used in the review with respect to pest host. Ash hosts are given in red, Pine in blue, Spruce in green, and any host not covered by the former in purple. The total number of articles for each host grouping is given within the respective bars.
Figure 3. Proportion of total articles used in the review with respect to (A) Numbers of Processes Modelled (Pest spread; Pest demography; Environmental factors; and Economic considerations and management) and (B) Number of Methods Used (mathematical models; statistical models; simulation studies; and Geographic Information Systems (GIS) studies (including use of imagery) used within the article.
Figure 4. Proportion of articles within major modeling techniques identified in the review. (A) shows the proportion of articles within major key processes (pest spread, pest demography, environmental factors, and economic considerations and management) in each panel with respect to the number of other modeled processes and by pest host (shown by colors). (B) shows the proportion of articles within major categories of modeling methods (statistical, simulation, mathematical, GIS) used in each panel with respect to the number of articles of other methods used and by pest host. Total number of articles given within each bar, respectively, of categorization.
We present our findings attending to the three main host species: ash (Fraxinus), spruce (Picea), pine (Pinus) given their dominance in the reviewed articles (see Figure 1). For each host, we break down selected highlights from the reviewed literature according to the key processes (Pest spread, Pest demography, Environmental factors, Economic considerations and management). Note that for our data regarding the “Other” host species, no other host was the subject of more than ten articles. For example, our data contained only seven studies of fir (Abies), and oak (Quercus), from which it is not possible to draw general conclusions.
3.1 Ash—Fraxinus
For Ash, our review procedure found only studies for a single pest species, EAB, providing a relatively simple place to start summarizing the literature.
3.1.1 Pest spread
In our review of the literature we found that satellite imagery and GIS tools have been relatively infrequently used to study EAB compared with the use of such methods in other pest-species interactions (see Figure 4B). Spectral analysis of ash tree crowns from different vegetation indices were compared to identify whether there is an EAB infestation (Zhang et al., 2014). The original site of introduction in the USA and initial spread has been reconstructed from data and statistical interpolation using GIS tools and dendrochronological analysis (Siegert et al., 2014). Lyttek et al. (2019) combined a partial differential equations model with GIS data to predict EAB spread. GIS has also been integrated with network automata to produce a novel framework termed Geographic Network Automata to model dispersal patterns (Anderson and Dragićević, 2020). However, most studies have focused on mathematical modeling of dispersal.
A variety of methods are used for such studies in relation to dispersal modeling, but a common approach is to use a negative exponential function. Mercader et al. (2009) incrementally measure distance and larvae within infested trees from a known initially infested location for which a negative exponential is a statistically good fit. This approach is largely adopted within spatial simulation studies accounting for distance between cells (Kovacs et al., 2011, 2014) and cellular automata/lattice models (Mercader et al., 2011, 2016; McCullough and Mercader, 2012) representing geographical US sub-counties (which are irregular cells) (Muirhead et al., 2006). This approach is also combined with attraction dynamics between EAB females and stressed ash trees (Mercader et al., 2012). Non-exponential kernels can also be used to represent such natural dispersal. Local dispersal has been modeled using a Gaussian kernel (analogous to normal distribution/random walk) (Withrow et al., 2015; Epanchin-Niell and Liebhold, 2015), a Laplace (double exponential distribution) (Lutscher and Musgrave, 2017), and an inverse power law distribution (Prasad et al., 2010). Fixed rate dispersal has been used within a wider modeling framework (Anderson and Dragićević, 2016a; Barlow et al., 2014).
Long-distance dispersal events can be incorporated with local spread in a single kernel function using distributions with non-exponentially bounded tails (Muirhead et al., 2006). However, often these events are associated with human activity so modeling is often conducted separately from local dispersal. Simple methods involve the addition of a fixed proportion traveling farther (Mercader et al., 2011; McCullough and Mercader, 2012). Gravity models predict long distance dispersal where EAB is ‘attracted” to focal sites, such as camp sites (Muirhead et al., 2006; Prasad et al., 2010) as a way to represent human mediated transport. Dispersal of EAB within infested wood through camp sites has been identified as a source of long distance dispersal and captured in models in a variety of ways. One study examined infested camp site locations and used bootstrapping methods to estimate an average distance from an initial location (Tobin et al., 2010). Human transport networks have also been investigated through network analysis (pathway analysis) as a means of spreading EAB (Yemshanov et al., 2012). Stochastic models have used random number generators to draw probabilities for long distance dispersal determined by proximity by road (Anderson and Dragićević, 2015). Hope et al. (2021) modeled both the short- and long-distance spread of EAB to evaluate the effectiveness of different levels of regulation. The short-distance spread used a beta distribution to simulate insect flight, typically within 20 km, whilst the long-distance spread was modeled using a Cauchy distribution to simulate human-assisted spread over much longer distances.
3.1.2 Pest demography
The simplest approach to representing demography is presence-absence within a gridded lattice, or municipality (Jones, 2019; Anderson and Dragićević, 2016a; Yemshanov et al., 2015) and is often used when the focus is on pest dispersal, or the effects of EAB presence. The presence-absence approach can be extended to include the probability of establishment in a cell after arrival, measuring local suitability suggestive of potential population size (Kovacs et al., 2011, 2014).
Empirical methods were developed to make use of data obtained by felling trees and counting larvae beneath the bark, thus enabling models that represent numbers of larvae per tree (Fahrner et al., 2017). These observations are conducted within a known core population and repeated over increasing distances toward the invasion front to produce data that can yield maximum likelihood estimation of density-dependent dispersal (Mercader et al., 2012). Another population modeling approach uses a step function to represent growth rates above and below a carrying capacity based on consumption of phloem (Mercader et al., 2011). This EAB growth function was adopted in further simulation approaches (McCullough and Mercader, 2012; Mercader et al., 2012).
Differential equation models—continuous in space and time—were used to study population dynamics, where logistic growth was considered with statistically estimated parameters for growth and carrying capacity (Barlow et al., 2014; McDermott and Finnoff, 2016; Jones et al., 2016). Integro-difference equations—discrete in time, but continuous in space—were also used to model population dynamics in a similar way (Lutscher and Musgrave, 2017), and a purely discrete in time difference equation for growth is used to iterate growth dynamics within a wider simulation (Fahrner et al., 2017). Some studies used EAB as a variable in relation to ash dynamics and management effort (Jones et al., 2016; Jones and McDermott, 2015).
A key theme in modeling EAB-ash interactions is the availability of phloem as a resource. One population modeling approach uses a step function to represent growth rates above and below a carrying capacity based on consumption of phloem (Mercader et al., 2011). This EAB growth function was adopted in further simulation approaches (McCullough and Mercader, 2012; Mercader et al., 2012). Simulation studies used statistically estimated parameters to account for phloem area (Fahrner et al., 2017; Mercader et al., 2011; McCullough and Mercader, 2012; Mercader et al., 2016; Kovacs et al., 2014). Joint EAB-ash dynamics were considered using coupled systems of differential equations (Jones et al., 2016; McDermott and Finnoff, 2016; Jones and McDermott, 2015), and in one case in a disease context as susceptible and infected (Barlow et al., 2014).
A comprehensive ABM constructed using Spatial Modeling Environment (SME, spatially similar to cellular automata) which can account for GIS layers and system dynamics using STELLA® (Costanza and Gottlieb, 1998). This modeling framework has been used on a small scale within DuPage County, Illinois, USA (BenDor and Metcalf, 2006). The framework accounts for the spread and population dynamics of EAB, as well as the distribution and population dynamics of ash. Information from GIS raster layers (scale 60 m × 60 m) about land use has been used to determine the suitability of ash growth. Parameters can also be adjusted to take into account quarantine zones for EAB management scenarios. Another ABM which examined EAB infestation in Ontario, Canada has been developed (Anderson and Dragićević, 2015). The model uses GIS layers for ash distributions, local features, and to initialize infested trees. Both trees and EAB are modeled as agents with their own dynamic processes. Dispersal is split into local radial dispersal and long-distance events generated probabilistically. Beetle demographics, such as mating, emergence, and death, are determined probabilistically. This model was the basis for other additions to the processes which included climate (Anderson and Dragićević, 2016a), the inclusion of EAB parasitoid predators (Anderson and Dragićević, 2016b), and the use network connections to model spread dynamics, i.e., infestation between trees (Anderson and Dragićević, 2018).
Finally, another key theme in modeling these EAB-ash interactions was the seasonality of activity. Degree-days are used in order to predict the emergence of beetles. MacDonald et al. (2022) used growing degree-day models, including double sine, single sine and standard degree-day models, to predict first adult emergence of the year. The sine wave method was similarly used by Webb et al. (2021) to determine emergence thresholds. Moreover, degree days were also adopted by Barker et al. (2023) as part of a Degree-Days, Risk, and Phenological event mapping (DDRP) platform to model the EAB life cycle. To compare seasonal activity between various ranges, Maxent used growing degree days combined with other climatic parameters to predict the distribution (Meshkova et al., 2023).
3.1.3 Environmental factors
Models were developed to predict cold-induced mortality for EAB using under-bark temperature. These consist of non-linear statistical models and Newtonian cooling models to predict how EAB would respond to the cold (Cuddington et al., 2018), and are used to develop a species distribution model for the North American continent. A variant of examining cold temperatures involved modeling snow cover effects on under bark temperature with thermal conduction models and iterative temperature equations (DeSantis et al., 2013). Species distribution models have been developed to predict EAB invasive potential using incidence data and climatic factors using Maxent and GARP methods (Sobek-Swant et al., 2012; Liang and Fei, 2014). Both studies also looked at future conditions under IPCC representative concentration pathways of future emissions. Species distribution models were similarly developed with Maxent methods by Dang et al. (2021) using future climate variables from the WorldClim dataset.
3.1.4 Economic considerations and management
The effectiveness of EAB management was examined in detail, examples include the SLAM strategy (SLow Ash Mortality) implemented in 2008 (Mercader et al., 2012), or studies that examined the effects of tree removal on beetle production (Fahrner et al., 2017). Other papers used a cost-benefit analysis to assess the efficiency of quarantined locations (Withrow et al., 2015; Vannatta et al., 2012).
Due to the EAB outbreak and subsequent invasion in the USA, economic analysis of management decisions has been at the forefront of modeling research. Bioeconomic models have been used to determine the optimal allocation of a given budget to the prevention and control EAB infestations, i.e. surveillance, treatment and removal (Bushaj et al., 2021), to improve efficiency on the allocation of surveillance efforts using a portfolio framework (Yemshanov et al., 2014, 2015), or to explore the decision of when to act, with the objective of minimized total cost (Kovacs et al., 2014). A more recent paper optimized the selection of daily routes and the number of host trees to inspect in order to minimmise the impact of EAB infestations using linear programming (Yemshanov et al., 2020).
In these bioeconomic studies, societal costs are typically those resulting from the removal or pesticide intervention and the damages associated with ecosystem services impacts (Kovacs et al., 2011; McDermott and Finnoff, 2016; Kovacs et al., 2014; McCullough and Siegert, 2007; Withrow et al., 2015). (Epanchin-Niell and Liebhold 2015) modeled long-term damages using a lag spread model to show how these are critically affected by invasion lags, damage persistence, and spread rates. Monetary estimates of damages have been provided with multiple techniques. McConnell et al. (2019) used an Impact Analysis for PLANning (IMPLAN) model, that takes an input-output approach to estimate the economic effects of an EAB infestation to the hardwood timber industry of Louisiana. This used a PERT-Beta distribution to predict Ash mortality and a double declining balance method for the depreciation of timber value over time. Environmental valuation techniques, based on econometric modeling have also been used. For example, hedonic pricing was used to isolate the effect of dying/dead trees on house prices from other structural or location features (house characteristics) (Li et al., 2019). Jones and McDermott (2018) used a dose-response relationship between EAB outbreaks and ambient air pollution concentrations to estimate the human health impacts. In a similar manner, Jones (2019) estimated EAB outbreak damages using increased healthcare costs from hospital visits due to loss of shade. Arnberger et al. (2020) used a discrete choice experiment, in order to better understand societal preferences toward EAB-impacted forest scenarios in Vienna, Austria, and Minneapolis, USA. A latent-class analysis was used to investigate visitors' preferences heterogeneity toward viewscape, social and biophysical elements of forest landscapes. More generally, the National Tree Benefits Calculator was used to estimate damages based on differences in ecosystem services provided under different management scenarios (Arbab et al., 2022).
3.2 Pine—Pinus
Studies where pine is the host make up the majority of studies (see Figure 2). Major wood borer pests of pine consist of Mountain pine beetle (MPB), Dendroctonus ponderosae; Southern pine beetle (SPB), Dendroctonus frontalis; and various Ips species, e.g. Pine engraver beetle, Ips pini, and Eastern 5-spined engraver Ips grandicollis. Studies of MPB are the most numerous, followed by SPB and a few related to various Ips species. Demographic modeling has been studied in more detail for MPB than for EAB. The reader should assume the study species is MPB unless otherwise stated.
3.2.1 Pest spread
A variety of methods are used to capture the dynamics of spread through pine; however, in general, estimation of spreading speed is not the aim of these studies. Both diffusion, accounting for local dispersal, and chemotaxis, accounting for the effects of attack pheromones, are found in numerous parameterised partial differential equation models, which also include growth terms (Strohm et al., 2016; Haran et al., 2015; White and Powell, 1997; Powell et al., 1996; Logan et al., 1998). Density-dependent dispersal is modeled using the product of diffusion and a host susceptibility function (Powell and Bentz, 2014). In such scenarios, weakened trees release attractive volatiles, and beetles move up the resulting concentration gradient. A similar dynamic for attack pheromones has also been included in some models; however, there is a threshold value at which the pheromones become repulsive. MPB population sizes have been estimated using dispersal models based on redistribution kernels and covariance structure (Koch et al., 2021). Integro-difference models use dispersal kernels to capture beetle movement with Gaussian distributions (Goodsman et al., 2016, 2018a). SPB dispersal has been similarly quantified by analyzing mark recapture data (Turchin and Thoeny, 1993). These normal-type distributions are best used to study local dynamics as they don't allow significant long range dispersal. However, a historical study examining clustering of MPB occupancy at the landscape scale reveals fractal distributions (Gamarra and He, 2008). This suggests MPB outbreaks can be modeled using fractal dispersal kernels in meta-populations. Simulation studies have made direct use of GIS/imagery data to identify and track “red attack” which can be traced in time to discern MPB spread (Wulder et al., 2009), and similarly using composite vegetation indices (Kaiser et al., 2013).
Statistical methods to quantify spread focus on spatial clustering or the extent of outbreaks using GIS/imagery data. Point process regression models are used in conjunction with environmental variables to determine pathways of spread (Giroday et al., 2012). GLMs were used to investigate how landscape and stand-level variables influence damage using a Gaussian distribution and also the spread and development of beetles using Akaike's information criterion (Zhan et al., 2023; Yu et al., 2022a). Cluster analysis can be used to estimate the spatial pattern of outbreak development (Aukema et al., 2006), using K- means clustering and weather data (Chapman et al., 2012), and have been used to analyse SPB and Ips sexdentatus outbreaks (Egan et al., 2016). Spatial autocorrelation is also used to this effect (Weed et al., 2017). Regression models can determine the likelihood of outbreak occurrence, particularly in relation to climatic factors (Sambaraju et al., 2012; Giroday et al., 2012; Horn et al., 2012), in Ips confusus (Santos and Whitham, 2010), and for SPB (Ungerer et al., 1999). ROC curve analysis was used to compare the Ips sexdentatus susceptibility maps produced by frequency ratio, analytical hierarchy process and logistic regression models (Sivrikaya et al., 2022).
3.2.2 Pest demography
Deterministic mathematical models are applied extensively to capture population dynamics. Discrete time Nicholson-Bailey dynamics using stochastic processes for beetle and susceptible host presence allow interactions between hosts and pest, which can be used to investigate attack threshold dynamics for a tree using parameter estimates from sampled data (Goodsman et al., 2016), and different strengths of density-dependence (Goodsman et al., 2017). Since tree mortality is related to host development, models which account for the age-structure of host interacting with beetle populations have been developed using Leslie matrices. These have been parameterised using maximum likelihood (Heavilin and Powell, 2008) and discrete annual phenology to measure outbreak intensity and recovery (Duncan et al., 2015). Differential equations are studied widely to understand pest dynamics, including: investigating aggregation dynamics using systems of reaction-diffusion equations (Strohm et al., 2016); the effectiveness of time-delay in disease control (Dong et al., 2023); SIR infection dynamics of hosts by pests through neighboring stands (Křivan et al., 2016); Allee effects and temperature-dependent emergence using age-structure (Powell and Bentz, 2014); and to develop models of phenotypic phenology using maximum likelihood parameter estimations for beetle development time (Yurk and Powell, 2010). Integral projection models aid in examining age structure and phenology of beetles (Goodsman et al., 2018a). Deterministic methods are employed to model SPB population dynamics using difference equations to determine the effects of temperature on population dynamics (Friedenberg et al., 2008); and to investigate the population dynamics of temperature modulated density-dependence (Turchin et al., 1991). Dynamic state variables can be applied to investigate how beetle fitness was impacted from energetic costs due to host selection (Chubaty et al., 2009). Two studies examining the case of fungi associated with MPB were found in this review. Respectively a deterministic and probabilistic model are applied to study Grosmannia clavigera and Ophiostoma montium within beetle outbreaks (Addison et al., 2015), and how temperature affects the mutualism (Addison et al., 2013). Euclidean distances were used in a path-finding algorithm to analyse dispersal patterns and investigate the role of human-vectored long-distance dispersal (Lee et al., 2022).
More complex models are less tractable to mathematical analysis, and researchers therefore tend to make use of simulation instead. Mathematical models are often simulated to corroborate results or investigate different parameter selections (Powell et al., 1996; White and Powell, 1997; Logan et al., 1998; Logan and Bentz, 1999; Powell et al., 2000; Powell and Logan, 2005; Powell and Bentz, 2014). A simulator called MPB-R has been developed, which simulates temperature-dependent phenology and beetle emergence and attack synchrony to predict subsequent generational growth rates (Bentz et al., 2016, 1991; Hicke et al., 2006). Bone et al. (2006) have developed a detailed cellular automata model that accounts for MPB phenology and dispersal, over-wintering temperature changes, and stand susceptibility to give an output of stand mortality. The model applies fuzzy set theory to account for uncertainty within the definitions of stand susceptibility, and from the uncertainty of GIS data input, e.g. interpretations of pixels. The model is then used to predict the spread of MPB infestation from GIS data. Brush and Lewis (2023) created a new model using survival, full and Hill functions to capture the key elements of MPB population dynamics. They used the model to simulate different scenarios to investigate forest resilience and the Allee threshold required to overcome them.
ABM models have also been developed for MPB dynamics. ForestSimMPB is an ABM which allows the user to input environmental data via GIS and is able to represent MPB populations at two spatial scales. Since it is a bottom-up model, outbreaks begin on the microscale, but evolve into larger infestations at the macroscale (Pérez and Dragicevic, 2011; Perez and Dragicevic, 2010). Lombardo et al. (2018) constructed an ABM of a SPB development rate model to predict larval development using seasonal temperature data, comparing the southern edge of the range with the northernmost range. Furthermore, Weiss et al. (2019) contructed a model to simulate the dynamics of the Monochamus alternatus beetle. Its various components describe the host availability and suitability of surrounding trees; take into account wind speed and direction; and simulate the entire beetle life cycle, including dispersal to nearby trees.
Statistical methods are primarily used to determine the effects of external factors, such as climate, on beetle demographics. Regression studies are most common when investigating pest demography. Studies which examine aspects of beetle life history include using regression to determine: the maximum rate of population growth and strength of density dependence (Trzcinski and Reid, 2009); changes in cyclical patterns (An and Gan, 2022); overwinter beetle survival rates (Bone and Nelson, 2019); sex-ratio skew within adult and larval populations (James et al., 2016); to determine the relationship between flight activity and climatic conditions (Brockerhoff et al., 2017); parameters for beetle infestation of a single tree in order to scale up to the stand level (Björklund et al., 2009); to compare the abundance of beetles between preferred, less preferred and mixed host trees (van Halder et al., 2022); and density of beetles emerging per m2 using experimental data (Negrón, 2019). Effects of temperature on beetle development are well studied and include: temperature dependent population growth rates (Weed et al., 2015); effects of cooling temperature on the mortality of eggs and larvae (Bleiker et al., 2017); and temperature dependent phenology and development using autologistic regression (Aukema et al., 2008). Predictions for attack densities have been calculated for Ips sexdentatus using regression models to determine temperature effects on development rate (Pineau et al., 2017b), and how offspring fitness can be affected by beetle density (Pineau et al., 2017a).
3.2.3 Environmental factors
Modeling that accounts for environmental factors focuses on temperature dependence for development, but also extends to temperature-related climatic factors. Effects of temperature dependence on beetle population dynamics can be captured using a variety of methods. Temperature dependent phenology is captured via parameter estimations for stochastic differential equations for individual beetle development (Yurk and Powell, 2010), and through integral projection models with parameterised air temperature models investigating cold mortality (Goodsman et al., 2018a). Statistical methods which determine the effects of temperature on phenology include: using regression models to determine the effects of temperature on outbreaks using GIS information and long term datasets (Weed et al., 2015); assessing how temperature affects non-linear population dynamics using regression models comparable to a Ricker model (Goodsman et al., 2018b); to determine the cooling mortality point in beetle eggs (Bleiker et al., 2017); and, in a spatio-temporal context with GIS temperature data, to determine phenology during an outbreak (Aukema et al., 2008). Bayesian statistics is used to analyse temperature and development times and uncover thresholds from MPB populations in the southern range of the USA (McManis et al., 2018). Seasonality can be studied using parameterised differential equation models where growth rate is dependent on temperature data linked to seasonality (Powell and Bentz, 2014; Powell and Logan, 2005; Hicke et al., 2006). Such studies use development rate curves developed from temperature data (Logan and Bentz, 1999).
Climate implicitly includes temperature, but also accounts for factors such as rainfall and/or humidity. Climatic effects of cooling on the development of SPB can be captured using difference models where growth is a statistical function of climate variables (Friedenberg et al., 2008), where the proportion that survive is determined by the cooling mortality temperature (Trân et al., 2007), and using a parameterised regressed Ricker's equation with temperature and rainfall estimates (Turchin et al., 1991).
Climatic factors can be investigated alongside other environmental factors, often characterized by GIS data. Spatial mortality measures of MPB have been determined using GLM models with binary response variables to predict new areas of MPB mortality (Liang et al., 2014). Logistic regression models are used to explain the spatio-temporal patterns of MPB outbreaks using annual temperature and drought variables (Preisler et al., 2012), to determine how temperature and elevation impact outbreaks (Sambaraju et al., 2012), predict current and future outbreaks from climate-beetle relationships and stand conditions (Buotte et al., 2016). Cooke and Carroll (2017) have used recruitment curves with correlated thermal response functions to make predictions of MPB spread. More recently, machine learning tools in the form of classification regression trees (CART) have been used with aerial imagery, climate, and elevation data to determine the mortality of Whitebark pine as a result of outbreaks (Jewett et al., 2011). This method has also been used to predict the extent of infestation within the Colorado Range Front of ponderosa pine (Negrón and Popp, 2004). Regression models and aerial imagery can determine which environmental variables are associated with outbreaks in MPB (Giroday et al., 2012; Wulder et al., 2006), and in Jeffrey pine beetle Dendroctonus jeffreyi (Egan et al., 2016). Various Ips species have applied regression methods to predict the mortality of Jeffrey pine in drought stricken Arizona (Negrón et al., 2009).
Climatic factors affecting SPB outbreaks can be included using logistic regression models to determine the probability of outbreaks and how severity is impacted by climatic and topographic conditions (Hernandez et al., 2012). Classification And Regression Trees (CART) methods can assess the predictive power of environmental covariates (Duehl et al., 2011). Other methods for determining outbreak success include PCA analysis of stand characteristics (Aoki et al., 2018).
Species distribution models (SDM) are a popular method to model distributions of MPB and Dendroctonus rhizophagus. Regression models can estimate lower temperature thresholds for SPB mortality to spatially determine northern range limits (Ungerer et al., 1999). Canonical discriminant functions have been used to develop maps for both SPB and MPB (Williams and Liebhold, 2002). This method uses discrete state variables as a linear function of covariates, in this case, infestation and climatic conditions. MaxEnt is widely used to develop distribution maps for several species, including D. rhizophagus using WorldClim climate variables (Fick and Hijmans, 2017) to predict distributions under current and future climate warming scenarios (Smith et al., 2013), in Ips mannsfeldi (Sarikaya et al., 2018), the Pine Wood Nematode (Bursaphelenchus xylophilus) (Wang et al., 2023; Tuomola et al., 2021; Lee et al., 2021) and comparing MPB, Dendroctonus breviconnis, and Ips pini (Evangelista et al., 2011). Comparisons between the SDM methods of MaxEnt BIOCLIM and ENFA using D. rhizophagus suggests subtle differences between the methods. In the context of modeling MPB a comparison of Maxent, GLM, and boosted regression trees (BRT) showed Maxent and BRT to have similar performance (Sidder et al., 2016).
Indirect studies of outbreaks which examine the effects on forests are widely studied using GIS and statistical methods. Efforts to predict the mortality of pine stands have been formed using vegetation indices from aerial imagery, regression, and GAM (Meddens and Hicke, 2014), and using vegetation indices from a time series of MODIS aerial imagery to predict host mortality (Meddens et al., 2012; Franklin et al., 2003). This has also been accomplished using GIS data and CART algorithms to predict stand mortality by examining forest stand characteristics (Vorster et al., 2017). Using MODIS imagery data, the spatial association of scalable hexagons (SASH) method (Potter et al., 2016), which identifies non-random temporal patterns in ecological phenomena, has been used to characterize stand mortality related to MPB outbreaks. One study used the statistical method, least absolute shrinkage and selection operator (LASSO), to determine the effects of climatic factors on outbreaks (Reyes et al., 2012). CART and Random Forest methods, alongside GIS techniques, can estimate the effects of MPB infestation and climatic factors which give rise to wildfires (Mietkiewicz and Kulakowski, 2016), quantifying the links between MPB outbreaks and wild fire (Liang et al., 2016). Partial least squares regression, using Landsat imagery, has been used to determine the importance of key climate variables that drive Tomicus spp. outbreaks (Yu et al., 2022b).
Simulation studies which do not examine beetle infestation directly include: Terrestrial Regional Ecosystem Exchange Simulator (TREES) to assess the effects on tree mortality from attack, specifically, blue stain fungal attack by Grosmannia clavigera (Millar et al., 2017); and to predict which stands are susceptible to attack under different fire management regimes using Spatially Explicit Model for LANdscape Dynamics (SEM-LAND) (Li et al., 2005). Integrated Biosphere Simulator (IBIS) is an ecosystem simulator which allows the user to estimate climate modeling and vegetation phenology in a single framework. Landry et al. (2016) use an extension called Marauding Insect Module (MIM), which accounts for MPB outbreaks, to estimate the effects of merchantable timber biomass, ecosystem carbon, surface albedo, and radiative forcing. Several studies used CLIMEX to simulate the impacts of climate on the development and distribution of beetle populations (Zhou et al., 2019; Yoon et al., 2023; Watt et al., 2011).
3.2.4 Economic consideration and management
Studies modeling management practices for forestry are more dominant than economic studies that take a broader societal perspective. Sims et al. (2010) used a deterministic bioeconomic model that accounts for timber and non-timber benefits to evaluate the optimal management choices of a central forest manager under MPB outbreaks driven by climatic events. A sensitivity analyses allowed to evaluate the effects of changes in the relative value of salvage timber and the public's valuation of timber vs. non-timber benefits. Waring et al. (2009) simulated outbreaks of SPB and the Mexican pine beetle Dendroctonus mexicanus using Forest Vegetation Simulator (FVS) with climatic data from weather stations to predict the dynamics of beetle infestation and economic costs in terms of revenues forgone. Partial differential equations have been simulated to model MPB within a managed stand where parameters were adjusted between different levels of effort to determine the effectiveness of different control strategies (Strohm et al., 2016). The effectiveness of thinning vs. non-thinning strategies has also been carry out by comparing three population modeling scenarios with hypothetical expansion factors (Coggins et al., 2011), and determining the change in beetle density through non-linear (exponential) regression models (MacQuarrie and Cooke, 2011). An individual-based model was developed by Xia et al. (2022) to simulate the spread of pine wilt disease and evaluate the effectiveness of infected tree removal among other control practices. Linear regression and aerial imagery were applied to evaluate the determinants of MPB population growth and compare the effectiveness of multiple management strategies (Trzcinski and Reid, 2009). Logarithmic regression can assess beetle attack densities on individual trees using field data, and findings can inform stand prioritization for MPB control (Björklund et al., 2009). McNichol et al. (2019) used Kaplan-Meier survival analysis and regression to examine the effects of prescribed burning as a management strategy.
Models were also used to inform strategies when the risk of wildfire outbreaks may be associated to beetle infestation. The interaction between beetles, climate, and management is complex and many such studies are simulation based. LANDIS-II has the capability of running simulations with extensions for beetle dynamics (MPB, Jeffrey pine beetle, Fir engraver beetle Scolytus ventralis), carbon dynamics, wildfire, and differing management strategies (Scheller et al., 2018); FVS, an MPB extension, and a fire and fuels extension can be used to develop a risk rating for wildfire occurrence under different thinning practices (Ager et al., 2007; Donato et al., 2013). BioSIM is another simulator which is a decision support tool for pest management planning. This simulator simulates weather suitability for MPB outbreaks to determine how inter-annual changes in weather affect stand mortality from outbreaks (Creeden et al., 2014).
Similarly to EAB, monetary estimates of societal damages have been estimated using a hedonic pricing approach that computes the change in house prices as a result of damage from MPB infestation on the surrounding environment (Price et al., 2010; Cohen et al., 2016). Early work used contingent valuation to estimate the willingness to pay of forest owners to participate in a cost sharing thinning programme for SPB control, using regression analysis to evaluate its socioeconomic determinants (Watson et al., 2013).
3.3 Spruce—Picea
Spruce pests found in our review are European spruce beetle (ESB) Ips typographus, Spruce beetle (SB) Dendroctonus rufinpennis, and two studies on Great spruce bark beetle Dendroctonus micans. Many of the studies originate in Central Europe, e.g., the Czech Republic, Slovakia, Germany, and Poland, where there are large coniferous forests. The majority of the studies focus on ESB, therefore, the reader should assume the subject is ESB unless otherwise stated.
3.3.1 Pest spread
Few models explicitly represent spread, but rather examine aspects of spread within a wider context. Diffusion has been used to model radial spread from timber storage as point sources using mark-recapture data to parameterise the diffusion coefficient (Skarpaas and Økland, 2009), and a parameterised diffusion coefficient for flying beetles that land on an infested tree (Byers, 1996). Dispersal kernels have also been calculated from mark-recapture methods (Dolezal et al., 2016). A parameterised gamma distribution can capture dispersal and aggregation within simulations (Økland et al., 2016).
Statistical methods focus on spatial or spatio-temporal aggregation to give an indication of pest spread. Infestation is often examined in the context of the forest and not primarily about pest dynamics. Non-linear regression is used to predict infestation gradients using distances from outbreak points (Potterf et al., 2019). Hierachial clustering was used to identify high risk pests based on their ecological similarities and geographic distributions to predict future establishment areas (Duffy et al., 2021). Spatial correlations of D. micans have been calculated to determine if outbreak patterns are related to stand characteristics, and if the densities of an obligate predator Rhizophagus grandis can estimate D. micans density (Gilbert and Grégoire, 2003). Concentric spread from point locations has been analyzed over consecutive years using the statistical process known as Kriging (Hlásny and Turčáni, 2013). These authors use sanitary felling data to establish a time series of spreading, which can be used to infer a rate of spread.
A few studies incorporate statistical methods and GIS/imagery data to model spread. This spread was captured by identifying the health of the host pine trees. Linear and partial least squares discriminant analysis was used to determine if spectral vegetation indices can distinguish between healthy, infested and susceptible trees (Trubin et al., 2023; Abdullah et al., 2019b). Principal component analysis and a random forest model were also used to assess the capability of spectral vegetation indices to determine the health or attack phase of a tree (Abdullah et al., 2019a; Mandl and Lang, 2023). Ali et al. (2021) investigated if canopy chlorophyll content could be used to estimate the health of a tree using an Invertible Forest Reflectance radiative transfer model.
3.3.2 Pest demography and environmental factors
We treat pest demography and environmental factors together since very few studies solely examined population processes without an environmental or climatic component using GIS data. Population growth rates have been estimated statistically using Generalized Linear Models (GLMs) to assess the impacts of biological factors using temperature and precipitation data, and biotic factors e.g. voltinism within the population (Marini et al., 2013). Other efforts to estimate growth rates include using oviposition data for different temperature regimes with a statistical exponential model (Wermelinger and Seifert, 1999, 1998). Similarly, Generalized Linear Mixed Models (GLMMs) are applied to predict the proportion of univoltine to multivoltine spruce beetles within stands by including stand characteristics e.g. host height, and air temperature (Hansen et al., 2001). Williams et al. (2021) also used GLMMs to investigate how beetle population dynamics are affected by air temperature and elevation (Williams et al., 2021). They've also been applied to investigate the effects of thermal variability between summer and winter on the flight behavior of SB and associated species (Dell and Davis, 2019). A brief review of parameters which are important when modeling life-cycle dynamics of ESB, such as thermal sum days, and voltinism can be found in (Netherer and Pennerstorfer, 2001).
Pest densities have been estimated by counting larvae within individual wind-felled trees and using linear regression to estimate the density across stand levels from survey data (Borkowski and Podlaski, 2011). Some examples of interactions between beetle species have been examined. Økland et al. (2009), 2011) used experimental data to parameterise simulations of beetle populations using iterative equations to determine if ESB could become invasive in the USA with spruce beetle (SB) facilitating through positive interaction. Competitive interactions between ESB and the six toothed spruce bark beetle Pityogenes chalcographus are studied using GLMs to predict tree deaths from infestation during a prescribed burning programme (Eriksson et al., 2006). Predictions for temperature-dependent development time of each larval instar for D. micans and it's obligate predator R. grandis are estimated from experimental data using non-linear regression (Gent et al., 2017). Jakoby et al. (2019) used a time-varying distributed delay model to predict the phenology of beetle populations and potential climatic effects.
Spatio-temporal predictions of infestations using environmental variables were examined using a variety of methods. Regression models used a combination of point pattern occurrence estimates with stand characteristics and storm damage (Stadelmann et al., 2014) and using stand characteristics and weather data (Stadelmann et al., 2013). Stepwise multiple regression can be applied alongside spatial autocorrelation and kriging to predict population densities of D. micans in relation to stand characteristics (Gilbert and Grégoire, 2003). Logistic regressions, using meteorological data, were used to identify forest stands susceptible to infestation (Mezei et al., 2019). GLM methods are used to predict population densities from stand characteristics and storm damage data (Kärvemo et al., 2014a; Eriksson et al., 2005). Predicting the effects of damage from infestation due to drought is studied using GLMMs to examine forest management data, climate and storm data, and water content of the soil (Pasztor et al., 2014; Overbeck and Schmidt, 2012).
Next, a number of studies that assess the impact of abiotic factors, such as storms, on the risk of outbreaks are discussed. Discrete time survival analysis methods can quantify risk from beetle population structure, variation in precipitation, spatial connectivity, and outbreak history (Seidl et al., 2016). CART methods are used to develop hazard models for SB outbreaks within individual stands using their characteristics (Reynolds and Holsten, 1994), and to predict historic outbreaks from climatic conditions (Hebertson and Jenkins, 2008). Probability models have been developed to estimate the chance a stand is infested using cumulative link models (proportional odds models) using stand characteristics and soil sample data as explanatory variables (Blomqvist et al., 2018).
Spatial correlation and GIS techniques are harnessed to identify or predict factors which may drive infestation. Causal relationships between environmental factors and SB outbreaks can be estimated using Random Forest methods, such as between temperature and drought variables, while simultaneously quantifying spatio-temporal synchrony (Hart et al., 2017, 2014), alongside LANDSAT imagery data to classify the extent of infestation over a defined time period (Latifi et al., 2014), and to map host mortality from SB attack (Woodward et al., 2018). A similar method, known as Boosted Regression Trees, can be combined with GLM methods to create risk maps of infestation using similar data (Kärvemo et al., 2014b). Similarly, evaluation of hyperspectral data from HyMAP can be achieved within ArcGIS to identify and map infestation (Lausch et al., 2013b). Methods have been developed using algorithms to analyse aerial imagery and correct the time-lag between forest disturbance and beetle outbreaks, which is verified on a historic data set (Kautz, 2014). LANDSAT data can be used to predict infestation and outbreaks by identifying the healthy (green) vs attacked trees from color gradients using Bayesian GLM and logistic regression (Hais et al., 2016). A number of techniques generate predicted future distributions of beetle infestation. Probability models and GIS data are used to create outbreak maps from stand infestation distance and temperature data (Kautz et al., 2013). Ortiz et al. (2013) illustrate this by comparing three different types of predictive distribution models (Maxent, GLM, Random Forests) from RapidEye and TerraSAR-X remote sensing tools.
Several ABM/IBM models have been developed to investigate various outbreaks. Although such models do not represent individual beetles, they do typically represent the individual characteristics of a location and can be linked to detailed mechanistic models of beetle development and/or behavior. A widely used example is a comprehensive phenological model for ESB called PHENIPS (Baier et al., 2007), which has been independently validated using climate data (Berec et al., 2013) and extensively linked to various ABMs such as iLand (Dobor et al., 2020). The model combines topographical (e.g. slope and elevation), and thermoclimatic (e.g. air temperature and solar radiation) variables with a parameterised beetle developmental model. PHENIPS can be applied using GIS data to create spatio-temporal models and is used within many studies to account for the phenology of ESB during outbreaks. Authors examine the effects of drought stress on tree stands (Matthews et al., 2018), with GAMMs to determine the effects of sanitary felling on outbreaks (Mezei et al., 2017a), using non-linear regression to determine which PHENIPS variables were most important at predicting timber loss from outbreaks (Mezei et al., 2017b), to predict future forest productivity under future IPCC RCP climate scenarios and to determine which modules are important for forest managers (Thiele et al., 2017), to compare population density before and after disturbance events (Fleischer et al., 2016). PHENIPS can also be combined with other models, such as LandClim, to further bolster the beetle dynamics module included within LandClim (Temperli et al., 2013).
An individual-based gap model for individual trees and forests FAREAST has a module for SB dynamics called UVAFME (University of Virginia Forest Model Enhanced) (Xiaodong and Shugart, 2005; Foster et al., 2018). UVAFME incorporates the probability of SB infestation from climatic inputs, stand and individual tree characteristics. Other ABMs include SAMBIA, Infestation Pattern Simulation (IPS), IPS-SPREADS and PICUS. SAMBIA is a bottom-up spatially explicit model at the stand scale that models beetles, trees, natural enemies, and forest management (Fahse and Heurich, 2011). SAMBIA accounts for beetle life cycle dynamics e.g. dispersal, breeding, natural and tree or management induced mortality. The model has been validated against historical outbreaks using spatio-temporal data collected in the Bavarian National Forest in Germany with results suggesting more research is required due to the complexity of interacting factors (Lausch et al., 2013a). PICUS is a spatially explicit forest patch model which accounts for different forest processes, such as ESB outbreaks in the form of additional modules based on PHENIPS (Seidl et al., 2007). PICUS has been used to examined the effects of ESB outbreaks on carbon sequestration and timber production under different climate and management scenarios (Seidl et al., 2008). It can also be used to integrate with another forest model EFISCEN, which operates at a larger scale than PICUS, to assess management strategies and climate change impacts (Seidl et al., 2009). Pietzsch et al. (2023) used IPS-SPREADS to simulate behavior and spread of ESB populations. The IPS-SPREADS model accounts for the interactions between individual spruce trees and beetles. Markov chains then used the simulation results to predict infestations further into the future. IPS is a bottom-up spatially explicit model which was developed to assess dispersal and infestation based on the interaction between tree and beetle (Kautz et al., 2014). Other modules have been developed to include windthrow dynamics, in particular, to investigate how cluster size and the spatial extent of windthrow affect ESB population dynamics (Potterf and Bone, 2017).
3.3.3 Economic considerations and management
Modeling studies focus on the impacts of management strategies on outbreaks. Dobor et al. (2020) and Zimova et al. (2020) used the landscape simulation model iLand to investigate the management strategies of sanitation logging and shortened rotation lengths, respectively. One study produced a bark beetle susceptibility index as input to a decision support system called Heureka so that users can simulate different management strategies to see how forest stands are affected (Nordkvist et al., 2023). More complex ecosystem simulators are used to investigate bark beetle attacks such as LandClim, and LPJ-GUESS. LandClim is a spatially explicit forest simulator used to investigate forest dynamics with fire, beetle, and wind disturbance. The beetle module has been parameterised using experimental data for SB outbreaks and climatic variables to predict outbreak dynamics in Colorado (Temperli et al., 2015). LPJ-GUESS is an assessment model for making management decisions regarding biodiversity, biogeochemistry, and vegetation while taking into account climatic changes. Jönsson et al. (2007), 2012, 2015) have adapted this model for use in beetle-attacked forests by including a phenological model for ESB to assess the importance of climate within outbreaks.
A combination of Generalized Additive Mixed Models (GAMMs) and PHENIPS can model the spatio-temporal effects of outbreaks from different sanitary felling zones, taking account of stand characteristics and climatic factors (Mezei et al., 2017a). Marini et al. (2012) used an iterative statistical model with historical forest health data on timber loss attributed to beetle attacks to inform management on the role of both biotic and abiotic controls on population dynamics. Management has been modeled in the SAMBIA framework where sanitation felling of infested trees occurs before beetle emergence by adjusting the beetle mortality component with a probability of detection and subsequent removal (Fahse and Heurich, 2011). PICUS and EFISCEN have been used to simulate 110 years of forest management by replacing spruce with non-coniferous alternatives to reduce outbreaks and investigate carbon sequestration dynamics (Seidl et al., 2008). Sanitary felling of infested trees is a common method used to curtail ongoing outbreaks. Ogris and Jurc (2010) use a machine learning algorithm called M5' to predict where sanitary felling should be conducted within a forest. The method uses 45 covariates to develop the prediction, including climatic variables, forest characteristics, topography, and soil, which they validated against existing spatial data.
Another important area of research is the study of pathways for pests to enter an ecosystem via global trade routes. (Skarpaas and Økland 2009) used a Gompertz model of bark beetle dynamics and diffusive dispersal to assess the risk of introductions and examine different pre-emptive measures available (e.g. storage enclosure) to mitigate this risk. A hedonic pricing framework has been used to assess the impact on property value from the joint effects of SB and wildfire risk using spatial econometric techniques (Hansen and Naughton, 2013). The effects of bark beetle attacks on timber revenues were simulated by Subramanian et al. (2016) using a Heureka-Standwise model.
4 A synthesis of wood borer pest modeling research
Our review identifies the key modeling techniques used to investigate the dynamics wood borer beetle infestations and the key pest-host systems and system processes considered in the literature. We systematically categorize studies by modeling methodology, by host and by modeled processes (see Figure 4). A growing number of hosts are being studied; before 2019, ash, pine, and spruce accounted for >99% of studies, whereas currently they only account for 88% of studies. The host species in the remaining studies are quite varied; the next most common hosts, oak and fir, only account for seven articles each.
The global distribution of studies was very uneven, with most papers addressing the insect infestations in the Northern Hemisphere, mostly the USA, Canada and Europe. Focal location and pest species were dominated by studies of pine within the continental USA and British Columbia, Canada, where vast pine forests are located, and significant impacts of wood borer beetles are observed. However, it is also likely that this focus also results from disparities in research funding within but especially between countries. Spruce studies focus on forest disturbance and impacts of infestation on such disturbances, with the effects of silvicultural practices and wind throw damage of particular interest. Studies of European spruce beetle (ESB) were located in Central Europe, while spruce beetle (SB) studies were conducted in the continental USA, and Emerald Ash Borer studies both in the USA and Russia. There were also several pest-host species combinations for which no models were found.
Our review focuses on the current state-of-the-art of bark beetle modeling approaches, rather than on historical developments. We carried out the temporal analysis using our classification of models, but we could not find any clear temporal trends, although there are possible shifts toward using more computationally demanding and data-demanding models.
4.1 Commonality of modeling challenges and approaches
The complexity of real-world wood borer pest-host systems requires modelers to make simplifying assumptions and these are typically based on common approaches. Such simplifications are driven by the questions addressed and thus vary between pest-host systems, as discussed. First, reducing the number of interacting variables of host and pest, e.g. investigating forest dynamics with pest presence/absence as a factor e.g. Blomqvist et al. (2018); Vorster et al. (2017); Jones (2019), or investigating beetle dynamics with simplified environmental (including host) characteristics, e.g. assuming the host has spatially uniform presence (Goodsman et al., 2018a; McDermott and Finnoff, 2016). Second, limiting scope according to study scale e.g. detailed population dynamics within a stand (Fahse and Heurich, 2011), but cruder demographics at larger scales e.g. probability of infestation across a national park (Potterf et al., 2019).
However, many modeling techniques are used across all host systems, reflecting that host-parasite systems share common life cycle elements; laying eggs beneath tree bark to feed on phloem until emergence. Much of the life cycle can be modeled using similar methods with minor system specific adjustments. For example, modeling dispersal via kernels which are used within all the pest species described (see Sections 3.1.1, 3.2.1, 3.3.1). This suggests that existing modeling studies can be used as a starting point for modeling emerging wood borer pests problems.
At large scales, static species distribution model (SDM) methods, such as MaxEnt and Genetic Algorithm for Rule Set Production (GARP), are applied to pests for all three main host species. These models used environmental covariates to establish species distributions, and crucially, can be used to estimate the suitability of other regions. For instance, the suitability of North America for EAB has been predicted from the native range using MaxEnt (Liang and Fei, 2014; Dang et al., 2021), or to predict future suitability for a pest (Smith et al., 2013; Ning et al., 2021).
Study of interventions to control infestation is common across systems. From a management perspective, cost-benefit analysis is a useful tool in guiding the application of limited resources (Yemshanov et al., 2015; Arbab et al., 2022). However, goals of management were markedly different across host species, for example, in ash the control of EAB spread is prioritized with plans such as SLow Ash Mortality (SLAM) (McCullough and Mercader, 2012), whereas in pine and spruce, the focus is use of silvicultural practices such as stand thinning to reduce outbreaks and fire hazards (Donato et al., 2013; Mezei et al., 2017a; Subramanian et al., 2016).
4.2 Modeling approaches driven by specific wood borer pest-host system properties
Differences in study scale vary by host with local (e.g. stand level) population dynamics models of ash and spruce pests less well researched as those in MPB. This is likely due to lack of data availability for associated pest species, which in turn, could be due to study focus, e.g. predicting ESB infestation from previous storm occurrences only requires infestation point data (Stadelmann et al., 2014). Use of remote sensing data was less prominent within ash studies, as these perform better in large coniferous forests, such as pine and spruce, that do not seasonally defoliate. In contrast, ash-EAB studies were often conducted on the smaller state-wide, or county-wide scale (Zhang et al., 2014; Siegert et al., 2014; Lyttek et al., 2019) compared to large-scale spruce or pine forests, for instance, in national parks (Potter et al., 2016; Liang et al., 2016; Pietzsch et al., 2023).
We identified a difference in research questions across systems. In the case of ash and EAB, the focus is primarily mortality of infested stands and the desire to control economic impacts (Kovacs et al., 2014; McConnell et al., 2019). This reflects the fact that infested ash is present within urban areas where municipal considerations, such as property valuations, removal and treatment of infestation, are a concern. Studies in pine and spruce focused on environmental impacts rather than economics, reflecting the scales over which environmental effects can be estimated effectively, and the relative absence of humans dwelling within vast forests. For spruce, the key issue is the broader impact of tree mortality on forest ecosystem function and productivity, e.g. carbon cycling (Seidl et al., 2008) and biodiversity (Thrippleton et al., 2023). For pine, the broader context of forest ecosystem science (Buotte et al., 2016; An and Gan, 2022) has been important, as is understanding pest-host dynamics e.g. Logan et al. (1998); Powell et al. (2000); Ozcan et al. (2022).
4.3 Model parameterisation and data gaps
Robust parameterisation of the bark beetle outbreak models is an essential and yet often overlooked part of the research linking to policy recommendations. Many outbreak models are highly nonlinear and sensitive to thresholds in beetle demography, host condition, and climate. Thus, without empirically constrained parameters, they cannot reliably support management decisions. In our review, we found different approaches from purely strategic, non-calibrated or literature-based mechanistic models and heuristically tuned landscape simulators to formally estimated frequentist risk models and fully Bayesian-calibrated mechanistic or hierarchical models, cf. (Table 2). More recently, ML-based predictors trained on large remote-sensing and monitoring datasets are increasingly playing a role. However, we found that many older models did not fit parameters, and most models that did make use of data were based on standard statistical frameworks rather than models of spread and establishment. Most long-term landscape and Earth-system models still depend predominantly on literature values or heuristic calibration. However, fine-scale and regional models now exploit diverse data (tree inventories, aerial surveys, remote sensing, climate, management records).
Table 2. Review of parameter estimation techniques with examples of papers deploying a particular approach.
Data used within modeling studies, whether directly to estimate model parameters or indirectly via statistical inference of parameters, comes from four different sources—laboratory experiments, field studies, remote sensing and operational/administrative data. Laboratory gathered data often investigates the effects of temperature on the development of pest species, e.g. Gent et al. (2017); Yurk and Powell (2010), which can be used directly to inform model parameters. Ideally, modeling and experimental methods should be considered in tandem with experimental studies designed explicitly to inform modeling, but it is often possible to make use of data collected for other reasons. Generally, the aim of field studies and remote sensing is more varied, often comparing effects across populations from different sites, and is usually not focused on the needs of models. Remote sensing is widely used within pine and spruce studies, but not for ash, where only five studies were found to use such techniques. Operational/administrative data can be used to inform on the scale and costs of impacts and interventions, including their distribution in time and space. They can also inform on key risk factors, e.g. trade or movement of people, to understand introduction or dispersal risk. When using field, remote sensing or operational data to inform statistical inference, care must be taken to adequately account for variability and bias induced by the data collection process. To aid future modeling, we now describe data and parameterisation in terms of: host and pest distributions; demographic processes including dispersal, environmental drivers, ecological impacts; and economic costs. The sources which estimated parameters were drawn from a small pool of empirical studies and adding to this limited data would also be a valuable goal of future work.
4.3.1 Distribution of host and pest
Understanding the distribution of hosts—and ideally the extent of pest problems—is critical to modeling impacts in the real world and requires ongoing research and surveillance. Remote sensing from satellite data still has important limitations, but can nonetheless be a useful tool to help in determining the spatial distribution of ground cover, e.g. susceptible forests and used to determine the occupancy of wood borer beetles and track spread. Databases such as MODIS/VIIRS, Sentinel-1, Sentinel-2, Landsat and Google satellite imagery can provide recent information about canopy cover. Land-use data (that classify land cover as e.g. agricultural, urban, forest - for example, the Copernicus CORINE program for Europe) can be used as environmental covariates in models to provide information about how pests are affected by land use (Bossard et al., 2000).
4.3.2 Demographic processes including births, deaths and dispersal
Reproduction rates have been estimated from laboratory or field studies, for example, temperature-dependent growth and development time (Powell and Bentz, 2014; Dell and Davis, 2019), and fecundity from egg counts (Lutscher and Musgrave, 2017). Wood borer species have a “rue carrying capacity”, where only a finite number of larvae can live within an occupied tree (Raffa et al., 2015; Anderbrant, 1990; Coulson, 1979) that has been estimated using diameter at breast height (DBH) of trees (Fahrner et al., 2017; Blomqvist et al., 2018; McCullough and Siegert, 2007). Quantifying competition/mutualism parameters remains a fundamental ongoing challenge and methods of estimation are not universally applied. Doing so would allow more effective assessment of the benefits of pest control from a competing species, or the effect of mutualism between MPB and blue stain fungi e.g. Grosmannia clavigera and Ophiostoma montium (Addison et al., 2013, 2015). Dispersal parameters are particularly difficult to obtain for flying insects such as wood borer beetles. Mark-recapture methods have been used widely for insect species, such as ESB (Dolezal et al., 2016), but logistical challenges remain, particularly with long distance dispersal events (Hagler and Jackson, 2001; Nathan et al., 2003). Flight mill experiments have been popular for estimating flight distance, but more of these would provide a greater pool to draw on for dispersal estimates (Taylor et al., 2010; Evenden et al., 2014).
4.3.3 Environmental drivers
Data from weather stations is helpful in determining environmental conditions at particular locations. Information databases can be used, such as WorldClim (Fick and Hijmans, 2017) or IPCC CMIP databases with predictions for future climatic conditions for data on a global scale (Meehl et al., 2000). Using such global datasets allows comparisons between different regions within the same framework of environmental covariates. This type of data was common among studies within this review, but challenges remain related to the scaling in landscape contexts. Often data must be aggregated from a fine scale to a coarse scale, which can lead to statistical bias, e.g. the “middle number problem” of ecology, where modeled relationships or parameter choices may not be valid in other locations or systems (Newman et al., 2019). CLIMEX is a useful software that also makes use of historical data and future climate estimates to simulate future distributions of pests (Zhou et al., 2019).
4.3.4 Economic costs
Economic information is required for cost-benefit analyses and portfolio optimisation. Most studies use the costing of a single tree, which is then scaled to the desired forest patch. These studies are primarily in ash where the cost is direct application of control measures (McCullough and Mercader, 2012; Hope et al., 2021). Information on house prices were used to value depreciation as a result of reducing the aesthetic appeal for potential buyers (Cohen et al., 2016; Kovacs et al., 2011, 2014). This type of information will have considerable uncertainty from effects not covered by the analysis, for example, the strength of demand and availability within the housing market. Fuchs et al. considered historical fluctuating timber prices as market risk along with bark beetle calamity prices in order to assess management options using portfolio optimisation (Fuchs et al., 2022).
5 Discussion
5.1 Research gaps
Despite the many advances in applying modeling to wood borer infestations, there remain significant gaps. Although the dynamics of host-parasite interactions have been well studied via simulation modeling (forward modeling), techniques to parameterise models (reverse engineer model parameters from data, also known as inverse modeling) in specific host-pest systems are required. Furthermore, ecological effects beyond host-parasite interactions e.g. competition and mutualism, were noticeably absent. We only found evidence of modeling the interaction between MPB and blue stain fungi e.g. Grosmannia clavigera and Ophiostoma montium (Addison et al., 2013, 2015). Examples of host parasitoid systems were only found in the context of biological control for ash-EAB systems (Anderson and Dragićević, 2016a).
Current modeling focuses on what happens after pest arrival, not on quantifying pathways for entry, although it is possible that our search criteria excluded such studies. Species distribution models (SDMs) have been used to quantify climatic suitability for wood borer pests such studies do not identify arrival sites and we found only a single study that examined trade as a point of entry for potential pests (Skarpaas and Økland, 2009). Modeling trade in general, and for entry points that could establish an invasive species, has been studied elsewhere (Dalmazzone and Giaccaria, 2014; Levine and D'Antonio, 2003). Applying these tools to international trade in timber products would help toward quantifying establishment risk from potential pests (Koch et al., 2011; Krishnankutty et al., 2020).
Economic models covered implementation of management strategies (McCullough and Mercader, 2012; Seidl et al., 2009; Strohm et al., 2016), healthcare costs (Jones and McDermott, 2018), and house price depreciation (Kovacs et al., 2014; Cohen et al., 2016), but modeling including the impact of pests on the economics of the timber industry were not well represented. Timber loss was included for pine and spruce hosts, with few studies including the impact of beetle outbreaks on salvage logging operations. We believe that inclusion of timber economics in studies modeling wood borer infestations is currently attainable and would be desirable for the timber industry such as from Fuchs et al. (2022).
5.2 Opportunities for further research
Sustained integration between researchers within different systems is needed. For example, spread models within EAB-ash systems is well developed, but to a much lesser extent when the host is spruce. Using techniques applied to EAB-ash systems would therefore contribute to a greater understanding of spruce-pest systems. The same can be said for pine-pest systems, where a deep understanding of many methods and problems gained via years of applications to MPB should be invaluable in the study of other wood borer-host systems. We nevertheless acknowledge that some relevant models may be missed given the search terms used in our search strategy, for example, the RITY models by Ogris et al. (2019).
A greater degree of generality within models would allow methods to be transferred between species. This is not always possible due to the particular nature of each host-pest system. The PHENIPS model for ESB is an example of a detailed framework which could be used within another system (Baier et al., 2007). However, even when such general frameworks are available, the key problem is how to parameterise them for specific systems. For example, Bayesian methods are not widely used in modeling wood borer pest, but could allow a better use of available data to estimate model parameters and quantify uncertainty in these estimates. A key advantage of Bayesian approaches is that they provide a toolkit to combine widely different sources of information e.g. various sources of hard data and expert opinion. Many of the techniques used in disease transmission modeling (see e.g. Swallow et al., 2022) could inform woodborer modeling, allowing the development of hierarchical models that account for variation at different scales, e.g. through fixed and random effects. However, restricted access to data sets hinders model parameterisation and access via repositories, such as the GBIF, would be highly desirable.
Despite some promising results e.g. classification of ash (62% successful detection) among other broadleaf species in typically urban settings (Zhang et al., 2014), use of remote sensing to detect infestation in highly heterogeneous contexts is a challenge. Key aspects of this challenge for remote sensing technologies include: host distribution, where for mixed stands it is also necessary to distinguish between host species; variation in foliage in deciduous host species; and variation in foliage due to pest damage.
Predicting potential invasion sites outside the native range using SDMs is an emerging area of research. However, static SDM methods should be used with caution since (i) they require focal species distributions to be at equilibrium, which is clearly not true for ongoing invasions (Václavík and Meentemeyer, 2012), and (ii) data from native range distributions, likely to be close to equilibrium, can be poor. Application of dynamic species distribution models to data on observed invasions is a promising area of research. Catterall et al. (2012) have developed a spatio-temporal model which uses the same types of data as SDMs, but allows for temporal variation and can provide similar outputs.
Use of complex ABMs allows a large number of factors within systems to be studied and facilitates the creation of more realistic scenarios. One of the criticisms of ABMs has been that an understanding of the effects any individual process are difficult to attain. However, developments in mathematical analysis are emerging, which allow analysis to be conducted (Cornell et al., 2019), opening further opportunities to analyse complex models and gain a deeper understanding of the processes involved. Furthermore, uncertainty quantification techniques, including history matching, can now be applied to assess complex ABMs against data (Vernon et al., 2022).
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
VK: Conceptualization, Data curation, Formal analysis, Writing – original draft. MK: Data curation, Formal analysis, Writing – original draft. GM: Conceptualization, Writing – review & editing. JT: Writing – review & editing. AK: Conceptualization, Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. MK, JT, GM, and AK contribution was funded by two UKRI projects, NE/V019988/1, Learning to adapt to an uncertain future: linking genes, trees, people and processes for more resilient treescapes (newLEAF) and NE/V020099/1, Connected Treescapes. GM was also funded as part of UKRI newLEAF project NE/V020005/1. VK was partly funded by the Scottish Plant Health Center project, PHC2018/14. GM and VK were supported by the Scottish Government's Rural and Environment Science and Analytical Services Division (RESAS). GM was also funded by the UKRI project NE/V020005/1 (newLEAF).
Acknowledgments
Yewen Chen and her contribution is the corroboratory analysis.
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors JT, AK declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/ffgc.2026.1706950/full#supplementary-material
References
Abdullah, H., Darvishzadeh, R., Skidmore, A. K., and Heurich, M. (2019a). Sensitivity of landsat-8 oli and tirs data to foliar properties of early stage bark beetle (Ips typographus, L.) infestation. Remote Sens. 11:398. doi: 10.3390/rs11040398
Abdullah, H., Skidmore, A. K., Darvishzadeh, R., and Heurich, M. (2019b). Sentinel-2 accurately maps green-attack stage of European spruce bark beetle (Ips typographus, L.) compared with landsat-8. Remote Sens. Ecol. Conservat. 5, 87–106. doi: 10.1002/rse2.93
Addison, A., Powell, J. A., Bentz, B. J., and Six, D. L. (2015). Integrating models to investigate critical phenological overlaps in complex ecological interactions: the mountain pine beetle-fungus symbiosis. J. Theor. Biol. 368, 55–66. doi: 10.1016/j.jtbi.2014.12.011
Addison, A. L., Powell, J. A., Six, D. L., Moore, M., and Bentz, B. J. (2013). The role of temperature variability in stabilizing the mountain pine beetle–fungus mutualism. J. Theor. Biol. 335, 40–50. doi: 10.1016/j.jtbi.2013.06.012
Ager, A. A., McMahan, A., Hayes, J. L., and Smith, E. L. (2007). Modeling the effects of thinning on bark beetle impacts and wildfire potential in the Blue Mountains of eastern Oregon. Landsc. Urban Plan. 80, 301–311. doi: 10.1016/j.landurbplan.2006.10.010
Ali, A. M., Abdullah, H., Darvishzadeh, R., Skidmore, A. K., Heurich, M., Roeoesli, C., et al. (2021). Canopy chlorophyll content retrieved from time series remote sensing data as a proxy for detecting bark beetle infestation. Remote Sens. Appl.-Soc. Environm. 22:100524. doi: 10.1016/j.rsase.2021.100524
An, H., and Gan, J. (2022). Periods and amplitudes of southern pine beetle infestations under climate change. Climate 10:126. doi: 10.3390/cli10090126
Anderbrant, O. (1990). Gallery construction and oviposition of the bark beetle Ips typographus (Coleoptera: Scolytidae) at different breeding densities. Ecol. Entomol. 15, 1–8. doi: 10.1111/j.1365-2311.1990.tb00777.x
Anderson, T., and Dragićević, S. (2015). An agent-based modeling approach to represent infestation dynamics of the emerald ash borer beetle. Ecol. Inform. 30, 97–109. doi: 10.1016/j.ecoinf.2015.09.003
Anderson, T., and Dragićević, S. (2016a). A geosimulation approach for data scarce environments: modeling dynamics of forest insect infestation across different landscapes. ISPRS Int. J. Geo-Inform. 5:9. doi: 10.3390/ijgi5020009
Anderson, T., and Dragićević, S. (2020). A geographic network automata approach for modeling dynamic ecological systems. Geogr. Anal. 52, 3–27. doi: 10.1111/gean.12183
Anderson, T. M., and Dragićević, S. (2016b). Geospatial pest-parasitoid agent based model for optimizing biological control of forest insect infestation. Ecol. Modell. 337, 310–329. doi: 10.1016/j.ecolmodel.2016.07.017
Anderson, T. M., and Dragićević, S. (2018). Network-agent based model for simulating the dynamic spatial network structure of complex ecological systems. Ecol. Modell. 389, 19–32. doi: 10.1016/j.ecolmodel.2018.10.008
Aoki, C. F., Cook, M., Dunn, J., Finley, D., Fleming, L., Yoo, R., et al. (2018). Old pests in new places: effects of stand structure and forest type on susceptibility to a bark beetle on the edge of its native range. For. Ecol. Manage. 419–420, 206–219. doi: 10.1016/j.foreco.2018.03.009
Arbab, N., Grabosky, J., and Leopold, R. (2022). Economic assessment of urban ash tree management options in new jersey. Sustainability 14:2172. doi: 10.3390/su14042172
Arnberger, A., Schneider, I. E., Eder, R., and Choi, A. (2020). Differences in urban forest visitor preferences for emerald ash borer-impacted areas. Forestry 93, 225–238. doi: 10.1093/forestry/cpz072
Aukema, B. H., Carroll, A. L., Zheng, Y., Zhu, J., Raffa, K. F., Moore, R. D., et al. (2008). Movement of outbreak populations of mountain pine beetle: influences of spatiotemporal patterns and climate. Ecography 31, 348–358. doi: 10.1111/j.0906-7590.2007.05453.x
Aukema, B. H., Carroll, A. L., Zhu, J., Raffa, K. F., Sickley, T. A., and Taylor, S. W. (2006). Landscape level analysis of mountain pine beetle in British Columbia, Canada: spatiotemporal development and spatial synchrony within the present outbreak. Ecography 29, 427–441. doi: 10.1111/j.2006.0906-7590.04445.x
Baier, P., Pennerstorfer, J., and Schopf, A. (2007). PHENIPS—A comprehensive phenology model of Ips typographus (L.) (Col., Scolytinae) as a tool for hazard rating of bark beetle infestation. For. Ecol. Manage. 249, 171–186. doi: 10.1016/j.foreco.2007.05.020
Barker, B. S., Coop, L., Duan, J. J., and Petrice, T. R. (2023). An integrative phenology and climatic suitability model for emerald ash borer. Front. Insect Sci. 3:1239173. doi: 10.3389/finsc.2023.1239173
Barlow, L.-A., Cecile, J., Bauch, C. T., and Anand, M. (2014). Modelling interactions between forest pest invasions and human decisions regarding firewood transport restrictions. PLoS ONE 9:e0090511. doi: 10.1371/journal.pone.0090511
BenDor, T. K., and Metcalf, S. S. (2006). The spatial dynamics of invasive species spread. Syst. Dynam. Rev. 22, 27–50. doi: 10.1002/sdr.328
Bentz, B. J., Duncan, J. P., and Powell, J. A. (2016). Elevational shifts in thermal suitability for mountain pine beetle population growth in a changing climate. Forestry: Int. J. Forest Res. 89, 271–283. doi: 10.1093/forestry/cpv054
Bentz, B. J., Logan, J. A., and Amman, G. D. (1991). Temperature-dependent development of the mountain pine beetle (Coleoptera: Scolytidae) and simulation of its phenology. Can. Entomol. 123, 1083–1094. doi: 10.4039/Ent1231083-5
Berec, L., Doležal, P., and Hais, M. (2013). Population dynamics of Ips typographus in the Bohemian Forest (Czech Republic): Validation of the phenology model PHENIPS and impacts of climate change. For. Ecol. Manage. 292, 1–9. doi: 10.1016/j.foreco.2012.12.018
Berryman, A. (1974). Dynamics of bark beetle populations: towards a general productivity model. Environm. Entomol. 3, 579–585. doi: 10.1093/EE/3.4.579
Beudert, B., Bernsteinová, J., Premier, J., and Bässler, C. (2018). Natural disturbance by bark beetle offsets climate change effects on streamflow in headwater catchments of the bohemian forest. Silva Gabreta 24, 21–45.
Björklund, N., Lindgren, B. S., Shore, T. L., and Cudmore, T. (2009). Can predicted mountain pine beetle net production be used to improve stand prioritization for management? For. Ecol. Manage. 257, 233–237. doi: 10.1016/j.foreco.2008.08.034
Bleiker, K. P., Smith, G. D., and Humble, L. M. (2017). Cold tolerance of mountain pine beetle (Coleoptera: Curculionidae) eggs from the historic and expanded ranges. Environm. Entomol. 46, 1165–1170. doi: 10.1093/ee/nvx127
Blomqvist, M., Kosunen, M., Starr, M., Kantola, T., Holopainen, M., and Lyytikäinen-Saarenmaa, P. (2018). Modelling the predisposition of Norway spruce to Ips typographus L. infestation by means of environmental factors in southern Finland. Eur. J. Forest Res. 137, 675–691. doi: 10.1007/s10342-018-1133-0
Bone, C., Dragicevic, S., and Roberts, A. (2006). A fuzzy-constrained cellular automata model of forest insect infestations. Ecol. Modell. 192, 107–125. doi: 10.1016/j.ecolmodel.2005.09.013
Bone, C., and Nelson, M. F. (2019). Improving mountain pine beetle survival predictions using multi-year temperatures across the western usa. Forests 10:866. doi: 10.3390/f10100866
Borkowski, A., and Podlaski, R. (2011). Statistical evaluation of Ips typographus population density: a useful tool in protected areas and conservation-oriented forestry. Biodivers. Conserv. 20, 2933–2951. doi: 10.1007/s10531-011-0121-1
Bossard, M., Feranec, J., and Otahel, J. (2000). Corine Land Cover Technical Guide: Addendum 2000. Technical Report No 40. Copenhagen: European Environment Agency.
Bretfeld, M., Speckman, H., Beverly, D., and Ewers, B. (2021). Bayesian predictions of bark beetle attack and mortality of three conifer species during epidemic and endemic population stages. Front. For. Glob. Change 4:679104. doi: 10.3389/ffgc.2021.679104
Brockerhoff, E. G., Chinellato, F., Faccoli, M., Kimberley, M., and Pawson, S. M. (2017). Effects of elevation and aspect on the flight activity of two alien pine bark beetles (Coleoptera: Curculionidae, Scolytinae) in recently-harvested pine forests. For. Ecol. Manage. 384:132–136. doi: 10.1016/j.foreco.2016.10.046
Brush, M., and Lewis, M. A. (2023). Coupling mountain pine beetle and forest population dynamics predicts transient outbreaks that are likely to increase in number with climate change. Bull. Math. Biol. 85:11. doi: 10.1007/s11538-023-01215-7
Buotte, P. C., Hicke, J. A., Preisler, H. K., Abatzoglou, J. T., Raffa, K. F., and Logan, J. A. (2016). Climate influences on whitebark pine mortality from mountain pine beetle in the Greater Yellowstone Ecosystem. Ecol. Appl. 26, 2507–2524. doi: 10.1002/eap.1396
Bushaj, S., Buyuktahtakin, I. E., Yemshanov, D., and Haight, R. G. (2021). Optimizing surveillance and management of emerald ash borer in urban environments. Nat. Resour. Model. 34:12267. doi: 10.1111/nrm.12267
Byers, J. A. (1996). An encounter rate model of bark beetle populations searching at random for susceptible host trees. Ecol. Modell. 91, 57–66. doi: 10.1016/0304-3800(95)00166-2
Cappaert, D., McCullough, D. G., Poland, T. M., and Siegert, N. W. (2005). Emerald ash borer in North America: a research and regulatory challenge. Am. Entomol. 51, 152–165. doi: 10.1093/ae/51.3.152
Catterall, S., Cook, A. R., Marion, G., Butler, A., and Hulme, P. E. (2012). Accounting for uncertainty in colonisation times: a novel approach to modelling the spatio-temporal dynamics of alien invasions using distribution data. Ecography 35, 901–911. doi: 10.1111/j.1600-0587.2011.07190.x
Chapman, T. B., Veblen, T. T., and Schoennagel, T. (2012). Spatiotemporal patterns of mountain pine beetle activity in the southern Rocky Mountains. Ecology 93, 2175–2185. doi: 10.1890/11-1055.1
Chubaty, A. M., Roitberg, B. D., and Li, C. (2009). A dynamic host selection model for mountain pine beetle, Dendroctonus ponderosae Hopkins. Ecol. Modell. 220, 1241–1250. doi: 10.1016/j.ecolmodel.2009.01.039
Coggins, S. B., Coops, N. C., Wulder, M. A., Bater, C. W., and Ortlepp, S. M. (2011). Comparing the impacts of mitigation and non-mitigation on mountain pine beetle populations. J. Environ. Manage. 92, 112–120. doi: 10.1016/j.jenvman.2010.08.016
Cohen, J., Blinn, C. E., Boyle, K. J., Holmes, T. P., and Moeltner, K. (2016). Hedonic valuation with translating amenities: mountain pine beetles and host trees in the colorado front range. Environm. Resource Econ. 63, 613–642. doi: 10.1007/s10640-014-9856-y
Cooke, B. J., and Carroll, A. L. (2017). Predicting the risk of mountain pine beetle spread to eastern pine forests: considering uncertainty in uncertain times. For. Ecol. Manage. 396, 11–25. doi: 10.1016/j.foreco.2017.04.008
Cornell, S. J., Suprunenko, Y. F., Finkelshtein, D., Somervuo, P., and Ovaskainen, O. (2019). A unified framework for analysis of individual-based models in ecology and beyond. Nat. Commun. 10:4716. doi: 10.1038/s41467-019-12172-y
Costanza, R., and Gottlieb, S. (1998). Modelling ecological and economic systems with stella: part II. Ecol. Modell. 112, 81–84. doi: 10.1016/S0304-3800(98)00073-8
Coulson, R. N. (1979). Population dynamics of bark beetles. Annual Rev. Entomol. 24, 417–447. doi: 10.1146/annurev.en.24.010179.002221
Creeden, E. P., Hicke, J. A., and Buotte, P. C. (2014). Climate, weather, and recent mountain pine beetle outbreaks in the western United States. For. Ecol. Manage. 312, 239–251. doi: 10.1016/j.foreco.2013.09.051
Cuddington, K., Sobek-Swant, S., Crosthwaite, J. C., Lyons, D. B., and Sinclair, B. J. (2018). Probability of emerald ash borer impact for Canadian cities and North America: a mechanistic model. Biol. Invasions 20, 2661–2677. doi: 10.1007/s10530-018-1725-0
Dalmazzone, S., and Giaccaria, S. (2014). Economic drivers of biological invasions: a worldwide, bio-geographic analysis. Ecol. Econ. 105, 154–165. doi: 10.1016/j.ecolecon.2014.05.008
Dang, Y.-Q., Zhang, Y.-L., Wang, X.-Y., Xin, B., Quinn, N. F., and Duan, J. J. (2021). Retrospective analysis of factors affecting the distribution of an invasive wood-boring insect using native range data: the importance of host plants. J. Pest Sci. 94, 981–990. doi: 10.1007/s10340-020-01308-5
Das, A. K., Baldo, M., Dobor, L., Seidl, R., Rammer, W., Modlinger, R., et al. (2025). The increasing role of drought as an inciting factor of bark beetle outbreaks can cause large-scale transformation of central European forests. Landscape Ecol. 40:108. doi: 10.1007/s10980-025-02125-w
Dell, I. H., and Davis, T. S. (2019). Effects of site thermal variation and physiography on flight synchrony and phenology of the north american spruce beetle (coleoptera: Curculionidae, scolytinae) and associated species in colorado. Environ. Entomol. 48, 998–1011. doi: 10.1093/ee/nvz067
DeSantis, R. D., Moser, W. K., Gormanson, D. D., Bartlett, M. G., and Vermunt, B. (2013). Effects of climate on emerald ash borer mortality and the potential for ash survival in North America. Agricult. Forest Meteorol. 178–179, 120–128. doi: 10.1016/j.agrformet.2013.04.015
Dobor, L., Hlasny, T., Rammer, W., Zimova, S., Barka, I., and Seidl, R. (2020). Spatial configuration matters when removing windfelled trees to manage bark beetle disturbances in central european forest landscapes. J. Environ. Manage. 254:109792. doi: 10.1016/j.jenvman.2019.109792
Dolezal, P., Okrouhlik, J., and Davidkova, M. (2016). Fine fluorescent powder marking study of dispersal in the spruce bark beetle, Ips typographus (Coleoptera: Scolytidae). Eur. J. Entomol. 113:1. doi: 10.14411/eje.2016.001
Donato, D. C., Simard, M., Romme, W. H., Harvey, B. J., and Turner, M. G. (2013). Evaluating post-outbreak management effects on future fuel profiles and stand structure in bark beetle-impacted forests of Greater Yellowstone. For. Ecol. Manage. 303, 160–174. doi: 10.1016/j.foreco.2013.04.022
Dong, R., Sui, H., and Ding, Y. (2023). Mathematical modeling and stability analysis of the delayed pine wilt disease model related to prevention and control. Mathematics 11:3705. doi: 10.3390/math11173705
Duehl, A. J., Koch, F. H., and Hain, F. P. (2011). Southern pine beetle regional outbreaks modeled on landscape, climate and infestation history. For. Ecol. Manage. 261, 473–479. doi: 10.1016/j.foreco.2010.10.032
Duffy, C., Tuffen, M. G., Fealy, R., and Griffin, C. T. (2021). Application of hierarchical clustering to identify high risk pests to sitka spruce: Ireland as a case study. Forestry 94, 86–101. doi: 10.1093/forestry/cpaa014
Duncan, J. P., Powell, J. A., Gordillo, L. F., and Eason, J. (2015). A model for mountain pine beetle outbreaks in an age-structured forest: predicting severity and outbreak-recovery cycle period. Bull. Math. Biol. 77, 1256–1284. doi: 10.1007/s11538-015-0085-5
Egan, J. M., Sloughter, J. M., Cardoso, T., Trainor, P., Wu, K., Safford, H., et al. (2016). Multi-temporal ecological analysis of Jeffrey pine beetle outbreak dynamics within the Lake Tahoe Basin. Popul. Ecol. 58, 441–462. doi: 10.1007/s10144-016-0545-2
Ehrenfeld, J. G. (2010). Ecosystem consequences of biological invasions. Annual Rev. Ecol. Evol. System. 41, 59–80. doi: 10.1146/annurev-ecolsys-102209-144650
Epanchin-Niell, R. S., and Liebhold, A. M. (2015). Benefits of invasion prevention: Effect of time lags, spread rates, and damage persistence. Ecol. Econ. 116, 146–153. doi: 10.1016/j.ecolecon.2015.04.014
Eriksson, M., Lilja, S., and Roininen, H. (2006). Dead wood creation and restoration burning: Implications for bark beetles and beetle induced tree deaths. For. Ecol. Manage. 231, 205–213. doi: 10.1016/j.foreco.2006.05.050
Eriksson, M., Pouttu, A., and Roininen, H. (2005). The influence of windthrow area and timber characteristics on colonization of wind-felled spruces by Ips typographus (L.). For. Ecol. Manage. 216, 105–116. doi: 10.1016/j.foreco.2005.05.044
Evangelista, P. H., Kumar, S., Stohlgren, T. J., and Young, N. E. (2011). Assessing forest vulnerability and the potential distribution of pine beetles under current and future climate scenarios in the Interior West of the US. For. Ecol. Manage. 262, 307–316. doi: 10.1016/j.foreco.2011.03.036
Evenden, M. L., Whitehouse, C. M., and Sykes, J. (2014). Factors influencing flight capacity of the mountain pine beetle (Coleoptera: Curculionidae: Scolytinae). Environm. Entomol. 43, 187–196. doi: 10.1603/EN13244
Fahrner, S. J., Abrahamson, M., Venette, R. C., and Aukema, B. H. (2017). Strategic removal of host trees in isolated, satellite infestations of emerald ash borer can reduce population growth. Urban Forest. Urban Green. 24, 184–194. doi: 10.1016/j.ufug.2017.03.017
Fahse, L., and Heurich, M. (2011). Simulation and analysis of outbreaks of bark beetle infestations and their management at the stand level. Ecol. Modell. 222, 1833–1846. doi: 10.1016/j.ecolmodel.2011.03.014
Fernández-Carrillo, A., Franco-Nieto, A., Yagüe-Ballester, M. J., and Gómez-Giménez, M. (2024). Predictive model for bark beetle outbreaks in European forests. Forests 15:1114. doi: 10.3390/f15071114
Fick, S. E., and Hijmans, R. J. (2017). Worldclim 2: new 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 37, 4302–4315. doi: 10.1002/joc.5086
Fleischer, P., Fleischer, P., Ferenčík, J., Hlaváč, P., and Kozánek, M. (2016). Elevated bark temperature in unremoved stumps after disturbances facilitates multi-voltinism in Ips typographus population in a mountainous forest. Forest. J. 62, 15–22. doi: 10.1515/forj-2016-0002
Foster, A. C., Shuman, J. K., Shugart, H. H., and Negron, J. (2018). Modeling the interactive effects of spruce beetle infestation and climate on subalpine vegetation. Ecosphere 9:e02437. doi: 10.1002/ecs2.2437
Franklin, S. E., Wulder, M. A., Skakun, R. S., and Carroll, A. L. (2003). Mountain Pine Beetle Red-Attack Forest Damage Classification Using Stratified Landsat TM Data in British Columbia, Canada. Bethesda: American Society for Photogrammetry and Remote Sensing.
Friedenberg, N. A., Sarkar, S., Kouchoukos, N., Billings, R. F., and Ayres, M. P. (2008). Temperature extremes, density dependence, and Southern Pine Beetle (Coleoptera: Curculionidae) population dynamics in East Texas. Environ. Entomol. 37, 650–659. doi: 10.1093/ee/37.3.650
Fuchs, J. M., Hittenbeck, A., Brandl, S., Schmidt, M., and Paul, C. (2022). Adaptation strategies for spruce forests-economic potential of bark beetle management and douglas fir cultivation in future tree species portfolios. Forestry 95, 229–246. doi: 10.1093/forestry/cpab040
Gamarra, J. G. P., and He, F. (2008). Spatial scaling of mountain pine beetle infestations. J. Animal Ecol. 77, 796–801. doi: 10.1111/j.1365-2656.2008.01389.x
Gent, C. A., Wainhouse, D., Day, K. R., Peace, A. J., and Inward, D. J. G. (2017). Temperature-dependent development of the great European spruce bark beetle Dendroctonus micans (Kug.) (Coleoptera: Curculionidae: Scolytinae) and its predator Rhizophagus grandis Gyll. (Coleoptera: Monotomidae: Rhizophaginae). Agricult. Forest Entomol. 19, 321–331. doi: 10.1111/afe.12212
Gilbert, M., and Grégoire, J.-C. (2003). Site condition and predation influence a bark beetle's success: a spatially realistic approach. Agricult. Forest Entomol. 5:87-96. doi: 10.1046/j.1461-9563.2003.00167.x
Giroday, H.-M. C., Carroll, A. L., and Aukema, B. H. (2012). Breach of the northern Rocky Mountain geoclimatic barrier: initiation of range expansion by the mountain pine beetle. J. Biogeog. 39, 1112–1123. doi: 10.1111/j.1365-2699.2011.02673.x
Goodsman, D. W., Aukema, B. H., McDowell, N. G., Middleton, R. S., and Xu, C. (2018a). Incorporating variability in simulations of seasonally forced phenology using integral projection models. Ecol. Evol. 8, 162–175. doi: 10.1002/ece3.3590
Goodsman, D. W., Cooke, B. J., and Lewis, M. A. (2017). Positive and negative density-dependence and boom-bust dynamics in enemy-victim populations: a mountain pine beetle case study. Theoret. Ecol. 10, 255–267. doi: 10.1007/s12080-017-0327-2
Goodsman, D. W., Grosklos, G., Aukema, B. H., Whitehouse, C., Bleiker, K. P., McDowell, N. G., et al. (2018b). The effect of warmer winters on the demography of an outbreak insect is hidden by intraspecific competition. Global Change Biol. 24, 3620–3628. doi: 10.1111/gcb.14284
Goodsman, D. W., Koch, D., Whitehouse, C., Evenden, M. L., Cooke, B. J., and Lewis, M. A. (2016). Aggregation and a strong Allee effect in a cooperative outbreak insect. Ecol. Appl. 26:2623-2636. doi: 10.1002/eap.1404
Gusenbauer, M., and Haddaway, N. R. (2020). Which academic search systems are suitable for systematic reviews or meta-analyses? Evaluating retrieval qualities of google scholar, pubmed, and 26 other resources. Res. Synthesis Methods 11, 181–217. doi: 10.1002/jrsm.1378
Hagler, J. R., and Jackson, C. G. (2001). Methods for marking insects: current techniques and future prospects. Annual Rev. Entomol. 46, 511–543. doi: 10.1146/annurev.ento.46.1.511
Hais, M., Wild, J., Berec, L., Br°una, J., Kennedy, R., Braaten, J., et al. (2016). Landsat imagery spectral trajectories—important variables for spatially predicting the risks of bark beetle disturbance. Remote Sens. 8:687. doi: 10.3390/rs8080687
Hansen, E. M., Bentz, B. J., and Turner, D. L. (2001). Temperature-based model for predicting univoltine brood proportions in spruce beetle (Coleoptera: Scolytidae). Can. Entomol. 133, 827–841. doi: 10.4039/Ent133827-6
Hansen, W. D., and Naughton, H. T. (2013). The effects of a spruce bark beetle outbreak and wildfires on property values in the wildland–urban interface of south-central Alaska, USA. Ecol. Econ. 96, 141–154. doi: 10.1016/j.ecolecon.2013.10.009
Haran, J., Roques, A., Bernard, A., Robinet, C., and Roux, G. (2015). Altitudinal barrier to the spread of an invasive species: could the pyrenean chain slow the natural spread of the pinewood nematode? PLoS ONE 10:e0134126. doi: 10.1371/journal.pone.0134126
Harati, S., Perez, L., and Molowny-Horas, R. (2020). Integrating neighborhood effect and supervised machine learning techniques to model and simulate forest insect outbreaks in British Columbia, Canada. Forests 11:1215. doi: 10.3390/f11111215
Hart, S. J., Veblen, T. T., Eisenhart, K. S., Jarvis, D., and Kulakowski, D. (2014). Drought induces spruce beetle (Dendroctonus rufipennis) outbreaks across northwestern Colorado. Ecology 95, 930–939. doi: 10.1890/13-0230.1
Hart, S. J., Veblen, T. T., Schneider, D., and Molotch, N. P. (2017). Summer and winter drought drive the initiation and spread of spruce beetle outbreak. Ecology 98, 2698–2707. doi: 10.1002/ecy.1963
Heavilin, J., and Powell, J. (2008). A novel method of fitting spatio-temporal models to data, with applications to the dynamics of mountain pine beetles. Nat. Res. Model. 21, 489–524. doi: 10.1111/j.1939-7445.2008.00021.x
Hebertson, E. G., and Jenkins, M. J. (2008). Climate factors associated with historic spruce beetle (Coleoptera: Curculionidae) outbreaks in Utah and Colorado. Environ. Entomol. 37, 281–292. doi: 10.1603/0046-225X(2008)37[281:CFAWHS]2.0.CO;2
Hernandez, A. J., Saborio, J., Ramsey, R. D., and Rivera, S. (2012). Likelihood of occurrence of bark beetle attacks on conifer forests in Honduras under normal and climate change scenarios. Geocarto Int. 27, 581–592. doi: 10.1080/10106049.2011.650652
Hicke, J. A., Logan, J. A., Powell, J., and Ojima, D. S. (2006). Changing temperatures influence suitability for modeled mountain pine beetle (Dendroctonus ponderosae) outbreaks in the western United States. J. Geophys. Res.: Biogeosci. 111:0101. doi: 10.1029/2005JG000101
Hill, L., Jones, G., Atkinson, N., Hector, A., Hemery, G., and Brown, N. (2019). The £15 billion cost of ash dieback in Britain. Curr. Biol. 29, R315–R316. doi: 10.1016/j.cub.2019.03.033
Hlássny, T., and Turčásni, M. (2013). Persisting bark beetle outbreak indicates the unsustainability of secondary Norway spruce forests: case study from Central Europe. Ann. For. Sci. 70, 481–491. doi: 10.1007/s13595-013-0279-7
Hope, E. S., McKenney, D. W., Pedlar, J. H., Lawrence, K., and MacDonald, H. (2021). Canadian efforts to slow the spread of emerald ash borer (Agrilus planipennis Fairmaire) are economically efficient. Ecol. Econ. 188:107126. doi: 10.1016/j.ecolecon.2021.107126
Horn, A., Kerdelhue, C., Lieutier, F., and Rossi, J.-P. (2012). Predicting the distribution of the two bark beetles tomicus destruens and tomicus piniperda in europe and the mediterranean region. Agric. For. Entomol. 14, 358–366. doi: 10.1111/j.1461-9563.2012.00576.x
Jakoby, O., Lischke, H., and Wermelinger, B. (2019). Climate change alters elevational phenology patterns of the European spruce bark beetle (Ips typographus). Glob. Chang. Biol. 25, 4048–4063. doi: 10.1111/gcb.14766
James, P. M. A., Janes, J. K., Roe, A. D., and Cooke, B. J. (2016). Modeling landscape-level spatial variation in sex ratio skew in the Mountain Pine Beetle (Coleoptera: Curculionidae). Environ. Entomol. 45, 790–801. doi: 10.1093/ee/nvw048
Jewett, J. T., Lawrence, R. L., Marshall, L. A., Gessler, P. E., Powell, S. L., and Savage, S. L. (2011). Spatiotemporal relationships between climate and whitebark pine mortality in the greater yellowstone ecosystem. Forest Sci. 57, 320–335. doi: 10.1093/forestscience/57.4.320
Jones, B. A. (2019). Tree shade, temperature, and human health: evidence from invasive species-induced deforestation. Ecol. Econ. 156, 12–23. doi: 10.1016/j.ecolecon.2018.09.006
Jones, B. A., and McDermott, S. M. (2015). Linking environmental management to health outcomes: a case study of the emerald ash borer. Appl. Econ. Letters 22, 1409–1414. doi: 10.1080/13504851.2015.1034836
Jones, B. A., and McDermott, S. M. (2018). Health impacts of invasive species through an altered natural environment: assessing air pollution sinks as a causal pathway. Environm. Resource Econ. 71, 23–43. doi: 10.1007/s10640-017-0135-6
Jones, B. A., McDermott, S. M., and Chermak, J. M. (2016). PLAN or get SLAM'ed: Optimal management of invasive species in the presence of indirect health externalities. J. Environ. Manage. 180, 538–550. doi: 10.1016/j.jenvman.2016.05.026
Jönsson, A. M., Harding, S., Bärring, L., and Ravn, H. P. (2007). Impact of climate change on the population dynamics of Ips typographus in southern Sweden. Agricult. Forest Meteorology 146, 70–81. doi: 10.1016/j.agrformet.2007.05.006
Jönsson, A. M., Lagergren, F., and Smith, B. (2015). Forest management facing climate change - an ecosystem model analysis of adaptation strategies. Mitigat. Adaptat. Strateg. Global Change 20, 201–220. doi: 10.1007/s11027-013-9487-6
Jönsson, A. M., Schroeder, L. M., Lagergren, F., Anderbrant, O., and Smith, B. (2012). Guess the impact of Ips typographus–An ecosystem modelling approach for simulating spruce bark beetle outbreaks. Agricult. Forest Meteorol. 166–167:188–200. doi: 10.1016/j.agrformet.2012.07.012
Kaiser, K. E., McGlynn, B. L., and Emanuel, R. E. (2013). Ecohydrology of an outbreak: mountain pine beetle impacts trees in drier landscape positions first. Ecohydrology 6, 444–454. doi: 10.1002/eco.1286
Kärvemo, S., Rogell, B., and Schroeder, M. (2014a). Dynamics of spruce bark beetle infestation spots: importance of local population size and landscape characteristics after a storm disturbance. For. Ecol. Manage. 334, 232–240. doi: 10.1016/j.foreco.2014.09.011
Kärvemo, S., Van Boeckel, T. P., Gilbert, M., Grégoire, J.-C., and Schroeder, M. (2014b). Large-scale risk mapping of an eruptive bark beetle – Importance of forest susceptibility and beetle pressure. For. Ecol. Manage. 318, 158–166. doi: 10.1016/j.foreco.2014.01.025
Kautz, M. (2014). On correcting the time-lag bias in aerial-surveyed bark beetle infestation data. For. Ecol. Manage. 326, 157–162. doi: 10.1016/j.foreco.2014.04.010
Kautz, M., Schopf, R., and Imron, M. A. (2014). Individual traits as drivers of spatial dispersal and infestation patterns in a host–bark beetle system. Ecol. Modell. 273, 264–276. doi: 10.1016/j.ecolmodel.2013.11.022
Kautz, M., Schopf, R., and Ohser, J. (2013). The “sun-effect”: microclimatic alterations predispose forest edges to bark beetle infestations. Eur. J. For. Res. 132, 453–465. doi: 10.1007/s10342-013-0685-2
Koch, D., Lewis, M. A., and Lele, S. (2021). The signature of endemic populations in the spread of mountain pine beetle outbreaks. Bull. Math. Biol. 83:65. doi: 10.1007/s11538-021-00899-z
Koch, F. H., Yemshanov, D., Colunga-Garcia, M., Magarey, R. D., and Smith, W. D. (2011). Potential establishment of alien-invasive forest insect species in the United States: where and how many? Biol. Invasions 13, 969–985. doi: 10.1007/s10530-010-9883-8
Kovacs, K. F., Haight, R. G., Mercader, R. J., and McCullough, D. G. (2014). A bioeconomic analysis of an emerald ash borer invasion of an urban forest with multiple jurisdictions. Resour. Energy Econ. 36, 270–289. doi: 10.1016/j.reseneeco.2013.04.008
Kovacs, K. F., Mercader, R. J., Haight, R. G., Siegert, N. W., McCullough, D. G., and Liebhold, A. M. (2011). The influence of satellite populations of emerald ash borer on projected economic costs in U.S. Communities, 2010–2020. J. Environm. Managem. 92, 2170–2181. doi: 10.1016/j.jenvman.2011.03.043
Krishnankutty, S. M., Bigsby, K., Hastings, J., Takeuchi, Y., Wu, Y., Lingafelter, S. W., et al. (2020). Predicting establishment potential of an invasive wood-boring beetle, Trichoferus campestris (Coleoptera: Cerambycidae) in the United States. Ann. Entomol. Soc. Am. 113, 88–99. doi: 10.1093/aesa/saz051
Křivan, V., Lewis, M., Bentz, B. J., Bewick, S., Lenhart, S. M., and Liebhold, A. (2016). A dynamical model for bark beetle outbreaks. J. Theor. Biol. 407, 25–37. doi: 10.1016/j.jtbi.2016.07.009
Kunegel-Lion, M., and Lewis, M. A. (2020). Mountain pine beetle outbreak duration and pine mortality depend on direct control effort. J. Environ. Manage. 260:110167. doi: 10.1016/j.jenvman.2020.110167
Landry, J.-S., Parrott, L., Price, D. T., Ramankutty, N., and Matthews, H. D. (2016). Modelling long-term impacts of mountain pine beetle outbreaks onmerchantable biomass, ecosystem carbon, albedo, and radiative forcing. Biogeosciences 13:5277–5295. doi: 10.5194/bg-13-5277-2016
Latifi, H., Schumann, B., Kautz, M., and Dech, S. (2014). Spatial characterization of bark beetle infestations by a multidate synergy of SPOT and Landsat imagery. Environ. Monit. Assess. 186, 441–456. doi: 10.1007/s10661-013-3389-7
Lausch, A., Heurich, M., and Fahse, L. (2013a). Spatio-temporal infestation patterns of Ips typographus (L.) in the Bavarian Forest National Park, Germany. Ecol. Indic. 31, 73–81. doi: 10.1016/j.ecolind.2012.07.026
Lausch, A., Heurich, M., Gordalla, D., Dobner, H. J., Gwillym-Margianto, S., and Salbach, C. (2013b). Forecasting potential bark beetle outbreaks based on spruce forest vitality using hyperspectral remote-sensing techniques at different scales. For. Ecol. Manage. 308, 76–89. doi: 10.1016/j.foreco.2013.07.043
Lee, D.-S., Choi, W. I., Nam, Y., and Park, Y.-S. (2021). Predicting potential occurrence of pine wilt disease based on environmental factors in south korea using machine learning algorithms. Ecol. Inform. 64:101378. doi: 10.1016/j.ecoinf.2021.101378
Lee, S.-J., Cho, H., Choi, Y., Choi, W. I., Chung, H. I., Lim, N. O., et al. (2022). Path-finding algorithm as a dispersal assessment method for invasive species with human-vectored long-distance dispersal event. Divers. Distrib. 28:1214–1226. doi: 10.1111/ddi.13524
Levine, J. M., and D'Antonio, C. M. (2003). Forecasting Biological Invasions with Increasing International Trade. Conserv. Biol. 17, 322–326. doi: 10.1046/j.1523-1739.2003.02038.x
Lewis, M., Nelson, W., and Xu, C. (2010). A structured threshold model for mountain pine beetle outbreak. Bullet. Mathem. Biol. 72, 565–589. doi: 10.1007/s11538-009-9461-3
Li, C., Barclay, H. J., Hawkes, B. C., and Taylor, S. W. (2005). Lodgepole pine forest age class dynamics and susceptibility to mountain pine beetle attack. Ecol. Complexity 2, 232–239. doi: 10.1016/j.ecocom.2005.03.001
Li, X., Holmes, T. P., Boyle, K. J., Crocker, E. V., and Nelson, C. D. (2019). Hedonic analysis of forest pest invasion: the case of emerald ash borer. Forests 10:820. doi: 10.3390/f10090820
Liang, L., and Fei, S. (2014). Divergence of the potential invasion range of emerald ash borer and its host distribution in North America under climate change. Clim. Change 122, 735–746. doi: 10.1007/s10584-013-1024-9
Liang, L., Hawbaker, T. J., Chen, Y., Zhu, Z., and Gong, P. (2014). Characterizing recent and projecting future potential patterns of mountain pine beetle outbreaks in the Southern Rocky Mountains. Appl. Geograp. 55, 165–175. doi: 10.1016/j.apgeog.2014.09.012
Liang, L., Hawbaker, T. J., Zhu, Z., Li, X., and Gong, P. (2016). Forest disturbance interactions and successional pathways in the Southern Rocky Mountains. For. Ecol. Manage. 375, 35–45. doi: 10.1016/j.foreco.2016.05.010
Linnakoski, R., and Forbes, K. M. (2019). Pathogens–the hidden face of forest invasions by wood-boring insect pests. Front. Plant Sci. 10:90. doi: 10.3389/fpls.2019.00090
Liu Y. and Estep, D.. (2025). A multiscale stochastic cellular automata model for dispersion of mountain pine beetles. R. Soc. Open Sci. 12:241934. doi: 10.1098/rsos.241934
Logan, J. A., and Bentz, B. J. (1999). Model analysis of mountain pine beetle (Coleoptera: Scolytidae) seasonality. Environ. Entomol. 28, 924–934. doi: 10.1093/ee/28.6.924
Logan, J. A., Régnière, J., and Powell, J. A. (2003). Assessing the impacts of global warming on forest pest dynamics. Front. Ecol. Environm. 1, 130–137. doi: 10.1890/1540-9295(2003)001[0130:ATIOGW]2.0.CO;2
Logan, J. A., White, P., Bentz, B. J., and Powell, J. A. (1998). Model analysis of spatial patterns in mountain pine beetle outbreaks. Theor. Popul. Biol. 53, 236–255. doi: 10.1006/tpbi.1997.1350
Lombardo, J. A., Weed, A. S., Aoki, C. F., Sullivan, B. T., and Ayres, M. P. (2018). Temperature affects phenological synchrony in a tree-killing bark beetle. Oecologia 188, 117–127. doi: 10.1007/s00442-018-4164-9
Lutscher, F., and Musgrave, J. A. (2017). Behavioral responses to resource heterogeneity can accelerate biological invasions. Ecology 98, 1229–1238. doi: 10.1002/ecy.1773
Lyttek, E., Lal, P., Nieddu, G., Forgoston, E., and Wieczerak, T. (2019). Modeling Agrilus planipennis F. (Coleoptera: Buprestidae) spread in New Jersey. J. Econ. Entomol. 112, 2482–2488. doi: 10.1093/jee/toz122
MacCalman, L., McKendrick, I. J., Denwood, M., Gibson, G., Catterall, S., Innocent, G., et al. (2016). MAPRA: modelling animal pathogens: review and adaptation. EFSA Support. Publicat. 13:1112E. doi: 10.2903/sp.efsa.2016.EN-1112
MacDonald, B., Baydack, R., Westwood, A. R., and Walker, D. (2022). Predicting emerald ash borer adult emergence and peak flight activity in Winnipeg, Manitoba, Canada. Front. Ecol. Evol. 10:846144. doi: 10.3389/fevo.2022.846144
MacQuarrie, C. J., and Cooke, B. J. (2011). Density-dependent population dynamics of mountain pine beetle in thinned and unthinned stands. Can. J. Forest Res. 41, 1031–1046. doi: 10.1139/x11-007
Mandl, L., and Lang, S. (2023). Uncovering early traces of bark beetle induced forest stress via semantically enriched sentinel-2 data and spectral indices. PFG-J. Photogrammet. Remote sens. Geoinformat. Sci. 91, 211–231. doi: 10.1007/s41064-023-00240-4
Marini, L., Ayres, M. P., Battisti, A., and Faccoli, M. (2012). Climate affects severity and altitudinal distribution of outbreaks in an eruptive bark beetle. Clim. Change 115, 327–341. doi: 10.1007/s10584-012-0463-z
Marini, L., Lindelöw, Å., Jönsson, A. M., Wulff, S., and Schroeder, L. M. (2013). Population dynamics of the spruce bark beetle: a long-term study. Oikos 122, 1768–1776. doi: 10.1111/j.1600-0706.2013.00431.x
Matthews, B., Netherer, S., Katzensteiner, K., Pennerstorfer, J., Blackwell, E., Henschke, P., et al. (2018). Transpiration deficits increase host susceptibility to bark beetle attack: Experimental observations and practical outcomes for Ips typographus hazard assessment. Agricultural and Forest Meteorol. 263:69–89. doi: 10.1016/j.agrformet.2018.08.004
McConnell, T. E., VanderSchaaf, C. L., and Tanger, S. M. (2019). Potential changes to louisiana hardwood timber industry economic contributions following emerald ash borer (eab) invasion: An input-output approach. J. Econ. Entomol. 112, 2751–2760. doi: 10.1093/jee/toz212
McCullough, D. G., and Mercader, R. J. (2012). Evaluation of potential strategies to SLow Ash Mortality (SLAM) caused by emerald ash borer (Agrilus planipennis): SLAM in an urban forest. Int. J. Pest Managem. 58, 9–23. doi: 10.1080/09670874.2011.637138
McCullough, D. G., and Siegert, N. W. (2007). Estimating potential emerald ash borer (Coleoptera: Buprestidae) populations using ash inventory data. J. Econ. Entomol. 100, 1577–1586. doi: 10.1603/0022-0493(2007)100[1577:EPEABC]2.0.CO;2
McDermott, S. M., and Finnoff, D. C. (2016). Impact of repeated human introductions andthe Allee effect on invasive species spread. Ecol. Modell. 329, 100–111. doi: 10.1016/j.ecolmodel.2016.03.001
McHugh, M. L. (2012). Interrater reliability: the kappa statistic. Biochemia Medica 22, 276–282. doi: 10.11613/BM.2012.031
McManis, A. E., Powell, J. A., and Bentz, B. J. (2018). Developmental parameters of a southern mountain pine beetle (Coleoptera: Curculionidae) population reveal potential source of latitudinal differences in generation time. Can. Entomol. 151, 1–15. doi: 10.4039/tce.2018.51
McNichol, B. H., Montes, C. R., Barnes, B. F., Nowak, J. T., Villari, C., and Gandhi, K. J. K. (2019). Interactions between southern Ips bark beetle outbreaks, prescribed fire, and loblolly pine (Pinus taeda L.) mortality. For. Ecol. Manage. 446, 164–174. doi: 10.1016/j.foreco.2019.05.036
Meddens, A. J. H., and Hicke, J. A. (2014). Spatial and temporal patterns of Landsat-based detection of tree mortality caused by a mountain pine beetle outbreak in Colorado, USA. For. Ecol. Manage. 322, 78–88. doi: 10.1016/j.foreco.2014.02.037
Meddens, A. J. H., Hicke, J. A., and Ferguson, C. A. (2012). Spatiotemporal patterns of observed bark beetle-caused tree mortality in British Columbia and the western United States. Ecol. Appl. 22, 1876–1891. doi: 10.1890/11-1785.1
Meehl, G. A., Boer, G. J., Covey, C., Latif, M., and Stouffer, R. J. (2000). The coupled model intercomparison project (cmip). Bull. Am. Meteorol. Soc. 81, 313–318. doi: 10.1175/1520-0477(2000)081<0313:TCMIPC>2.3.CO;2
Mercader, R. J., McCullough, D. G., Storer, A. J., Bedford, J. M., Heyd, R., Siegert, N. W., et al. (2016). Estimating local spread of recently established emerald ash borer, Agrilus planipennis, infestations and the potential to influence it with a systemic insecticide and girdled ash trees. For. Ecol. Manage. 366, 87–97. doi: 10.1016/j.foreco.2016.02.005
Mercader, R. J., Siegert, N. W., Liebhold, A. M., and McCullough, D. G. (2009). Dispersal of the emerald ash borer, Agrilus planipennis, in newly-colonized sites. Agricult. Forest Entomol. 11, 421–424. doi: 10.1111/j.1461-9563.2009.00451.x
Mercader, R. J., Siegert, N. W., Liebhold, A. M., and McCullough, D. G. (2011). Influence of foraging behavior and host spatial distribution on the localized spread of the emerald ash borer, Agrilus planipennis. Popul. Ecol. 53, 271–285. doi: 10.1007/s10144-010-0233-6
Mercader, R. J., Siegert, N. W., and McCullough, D. G. (2012). Estimating the influence of population density and dispersal behavior on the ability to detect and monitor agrilus planipennis (Coleoptera: Buprestidae) populations. J. Econ. Entomol. 105, 272–281. doi: 10.1603/EC11172
Meshkova, V., Borysenko, O., Kucheryavenko, T., Skrylnyk, Y., Davydenko, K., and Holusa, J. (2023). Potential westward spread of emerald ash borer, Agrilus planipennis Fairmaire, 1888 (Coleoptera: Buprestidae) from Eastern Ukraine. Forests 14:736. doi: 10.3390/f14040736
Mezei, P., Blaženec, M., Grodzki, W., Škvarenina, J., and Jakuš, R. (2017a). Influence of different forest protection strategies on spruce tree mortality during a bark beetle outbreak. Ann. For. Sci. 74:65. doi: 10.1007/s13595-017-0663-9
Mezei, P., Jakuš, R., Pennerstorfer, J., Havašová, M., Škvarenina, J., Ferenčík, J., et al. (2017b). Storms, temperature maxima and the Eurasian spruce bark beetle Ips typographus–An infernal trio in Norway spruce forests of the Central European High Tatra Mountains. Agricult. Forest Meteorol. 242, 85–95. doi: 10.1016/j.agrformet.2017.04.004
Mezei, P., Potterf, M., Skvarenina, J., Rasmussen, J. G., and Jakus, R. (2019). Potential solar radiation as a driver for bark beetle infestation on a landscape scale. Forests 10:604. doi: 10.3390/f10070604
Mietkiewicz, N., and Kulakowski, D. (2016). Relative importance of climate and mountain pine beetle outbreaks on the occurrence of large wildfires in the western USA. Ecol. Appl. 26, 2525–2537. doi: 10.1002/eap.1400
Millar, D. J., Ewers, B. E., Mackay, D. S., Peckham, S., Reed, D. E., and Sekoni, A. (2017). Improving ecosystem-scale modeling of evapotranspiration using ecological mechanisms that account for compensatory responses following disturbance. Water Resour. Res. 53, 7853–7868. doi: 10.1002/2017WR020823
Muirhead, J. R., Leung, B., Overdijk, C., Kelly, D. W., Nandakumar, K., et al. (2006). Modelling local and long-distance dispersal of invasive emerald ash borer Agrilus planipennis (Coleoptera) in North America. Divers. Distribut. 12, 71–79. doi: 10.1111/j.1366-9516.2006.00218.x
Müller, M., Olsson, P.-O., Eklundh, L., Jamali, S., and Ardö, J. (2022). Features predisposing forest to bark beetle outbreaks and their dynamics during drought. Forest Ecol. Managem. 523:120480. doi: 10.1016/j.foreco.2022.120480
Murray, J. D. (2001). Mathematical Biology: I. An Introduction, volume 3. Cham: Springer Science & Business Media.
Musolin, D. L., Selikhovkin, A. V., Shabunin, D. A., Zviagintsev, V. B., and Baranchikov, Y. N. (2017). Between Ash dieback and Emerald Ash Borer: Two Asian Invaders in Russia and the Future of Ash in Europe.
Nathan, R., Perry, G., Cronin, J. T., Strand, A. E., and Cain, M. L. (2003). Methods for estimating long-distance dispersal. Oikos 103, 261–273.
Negrón, J. F. (2019). Biological Aspects of Mountain Pine Beetle in Lodgepole Pine Stands of Different Densities in Colorado, USA. Forests 10:18. doi: 10.3390/f10010018
Negrón, J. F., McMillin, J. D., Anhold, J. A., and Coulson, D. (2009). Bark beetle-caused mortality in a drought-affected ponderosa pine landscape in Arizona, USA. For. Ecol. Manage. 257, 1353–1362. doi: 10.1016/j.foreco.2008.12.002
Negrón, J. F., and Popp, J. B. (2004). Probability of ponderosa pine infestation by mountain pine beetle in the Colorado Front Range. For. Ecol. Manage. 191, 17–27.
Netherer, S., and Pennerstorfer, J. (2001). Parameters relevant for modelling the potential development of Ips typographus (Coleoptera: Scolytidae). Integrat. Pest Managem. Rev. 6, 177–184. doi: 10.1023/A:1025719301446
Newman, E. A., Kennedy, M. C., Falk, D. A., and McKenzie, D. (2019). Scaling and complexity in landscape ecology. Front. Ecol. Evol. 7:293. doi: 10.3389/fevo.2019.00293
Ning, H., Tang, M., and Chen, H. (2021). Impact of climate change on potential distribution of chinese white pine beetle Dendroctonus armandi in China. Forests 12:544. doi: 10.3390/f12050544
Nordkvist, M., Eggers, J., Fustel, T. L.-A., and Klapwijk, M. J. (2023). Development and implementation of a spruce bark beetle susceptibility index: a framework to compare bark beetle susceptibility on stand level. Trees Forests People 11:100364. doi: 10.1016/j.tfp.2022.100364
Ogris, N., Ferlan, M., Hauptman, T., Pavlin, R., Kavčič, A., Jurc, M., et al. (2019). Rity – a phenology model of Ips typographus as a tool for optimization of its monitoring. Ecol. Modell. 410:108775. doi: 10.1016/j.ecolmodel.2019.108775
Ogris, N., and Jurc, M. (2010). Sanitary felling of Norway spruce due to spruce bark beetles in Slovenia: a model and projections for various climate change scenarios. Ecol. Modell. 221, 290–302. doi: 10.1016/j.ecolmodel.2009.05.015
Økland, B., Erbilgin, N., Skarpaas, O., Christiansen, E., and Løangström, B. (2011). Inter-species interactions and ecosystem effects of non-indigenous invasive and native tree-killing bark beetles. Biol. Invasions 13, 1151–1164. doi: 10.1007/s10530-011-9957-2
Økland, B., Krokene, P., and Lange, H. (2012). The effect of climate change on the spruce bark beetle. Science Nordic 27, 1–5.
Økland, B., Nikolov, C., Krokene, P., and Vakula, J. (2016). Transition from windfall- to patch-driven outbreak dynamics of the spruce bark beetle Ips typographus. For. Ecol. Manage. 363, 63–73. doi: 10.1016/j.foreco.2015.12.007
Økland, B., Skarpaas, O., and Kausrud, K. (2009). Threshold facilitations of interacting species. Popul. Ecol. 51:513. doi: 10.1007/s10144-009-0141-9
Ortiz, S. M., Breidenbach, J., and Kändler, G. (2013). Early detection of bark beetle green attack using TerraSAR-X and RapidEye data. Remote Sens. 5, 1912–1931. doi: 10.3390/rs5041912
Overbeck, M., and Schmidt, M. (2012). Modelling infestation risk of Norway spruce by Ips typographus (L.) in the Lower Saxon Harz Mountains (Germany). For. Ecol. Manage. 266:115–125. doi: 10.1016/j.foreco.2011.11.011
Ozcan, G. E., Sivrikaya, F., Sakici, O. E., and Enez, K. (2022). Determination of some factors leading to the infestation of Ips sexdentatus in crimean pine stands. For. Ecol. Manage. 519:120316. doi: 10.1016/j.foreco.2022.120316
Pasztor, F., Matulla, C., Rammer, W., and Lexer, M. J. (2014). Drivers of the bark beetle disturbance regime in Alpine forests in Austria. For. Ecol. Manage. 318, 349–358. doi: 10.1016/j.foreco.2014.01.044
Perez, L., and Dragicevic, S. (2010). Modeling mountain pine beetle infestation with an agent-based approach at two spatial scales. Environm. Modell. Softw. 25, 223–236. doi: 10.1016/j.envsoft.2009.08.004
Pérez, L., and Dragicevic, S. (2011). ForestSimMPB: a swarming intelligence and agent-based modeling approach for mountain pine beetle outbreaks. Ecol. Inform. 6, 62–72. doi: 10.1016/j.ecoinf.2010.09.003
Pietzsch, B. W., Wudel, C., and Berger, U. (2023). Nonparametric upscaling of bark beetle infestations and management from plot to landscape level by combining individual-based with markov chain models. Eur. J. For. Res. 142, 129–144. doi: 10.1007/s10342-022-01512-1
Pineau, X., Bourguignon, M., Jactel, H., Lieutier, F., and Sallé, A. (2017a). Pyrrhic victory for bark beetles: successful standing tree colonization triggers strong intraspecific competition for offspring of Ips sexdentatus. For. Ecol. Manage. 399, 188–196. doi: 10.1016/j.foreco.2017.05.044
Pineau, X., David, G., Peter, Z., Sallé, A., Baude, M., Lieutier, F., et al. (2017b). Effect of temperature on the reproductive success, developmental rate and brood characteristics of Ips sexdentatus (Boern.). Agricult. Forest Entomol. 19, 23–33. doi: 10.1111/afe.12177
Potter, K. M., Koch, F. H., Oswalt, C. M., and Iannone, B. V. (2016). Data, data everywhere: detecting spatial patterns in fine-scale ecological information collected across a continent. Landsc. Ecol. 31, 67–84. doi: 10.1007/s10980-015-0295-0
Potterf, M., and Bone, C. (2017). Simulating bark beetle population dynamics in response to windthrow events. Ecol. Complex. 32, 21–30. doi: 10.1016/j.ecocom.2017.08.003
Potterf, M., Nikolov, C., Kočická, E., Ferenčík, J., Mezei, P., and Jakuš, R. (2019). Landscape-level spread of beetle infestations from windthrown- and beetle-killed trees in the non-intervention zone of the Tatra National Park, Slovakia (Central Europe). For. Ecol. Manage. 432, 489–500. doi: 10.1016/j.foreco.2018.09.050
Powell, J., Kennedy, B., White, P., Bentz, B., Logan, J., and Roberts, D. (2000). Mathematical elements of attack risk analysis for mountain pine beetles. J. Theor. Biol. 204, 601–620. doi: 10.1006/jtbi.2000.1080
Powell, J. A., and Bentz, B. J. (2014). Phenology and density-dependent dispersal predict patterns of mountain pine beetle (Dendroctonus ponderosae) impact. Ecol. Modell. 273, 173–185. doi: 10.1016/j.ecolmodel.2013.10.034
Powell, J. A., and Logan, J. A. (2005). Insect seasonality: circle map analysis of temperature-driven life cycles. Theor. Popul. Biol. 67, 161–179. doi: 10.1016/j.tpb.2004.10.001
Powell, J. A., Logan, J. A., and Bentz, B. J. (1996). Local projections for a global model of mountain Pine Beetle Attacks. J. Theor. Biol. 179, 243–260. doi: 10.1006/jtbi.1996.0064
Prasad, A. M., Iverson, L. R., Peters, M. P., Bossenbroek, J. M., Matthews, S. N., Davis Sydnor, T., et al. (2010). Modeling the invasive emerald ash borer risk of spread using a spatially explicit cellular model. Landsc. Ecol. 25, 353–369. doi: 10.1007/s10980-009-9434-9
Preisler, H. K., Hicke, J. A., Ager, A. A., and Hayes, J. L. (2012). Climate and weather influences on spatial temporal patterns of mountain pine beetle populations in Washington and Oregon. Ecology 93, 2421–2434. doi: 10.1890/11-1412.1
Price, J. I., McCollum, D. W., and Berrens, R. P. (2010). Insect infestation and residential property values: a hedonic analysis of the mountain pine beetle epidemic. Forest Policy Econ. 12, 415–422. doi: 10.1016/j.forpol.2010.05.004
Pullin, A. S., Frampton, G. K., Livoreil, B., and Petrokofsky, G. (2022). Guidelines and Standards for Evidence Synthesis in Environmental Management. Available online at: https://www.environmentalevidence.org/information-for-authors (Accessed September 1, 2024).
Raffa, K. F., Gregoire, J.-C., and Lindgren, B. S. (2015). Natural History and Ecology of Bark Beetles. London: Elsevier, 1–40.
Ramazi, P., Kunegel-Lion, M., Greiner, R., and Lewis, M. A. (2021). Predicting insect outbreaks using machine learning: a mountain pine beetle case study. Ecol. Evolut. 11, 13014–13028. doi: 10.1002/ece3.7921
Ramsfield, T. D., Bentz, B. J., Faccoli, M., Jactel, H., and Brockerhoff, E. G. (2016). Forest health in a changing world: effects of globalization and climate change on forest insect and pathogen impacts. Forestry: Int. J. Forest Res. 89, 245–252. doi: 10.1093/forestry/cpw018
Reyes, P. E., Zhu, J., and Aukema, B. H. (2012). Selection of spatial-temporal lattice models: assessing the impact of climate conditions on a mountain pine beetle outbreak. J. Agric. Biol. Environ. Stat. 17, 508–525. doi: 10.1007/s13253-012-0103-0
Reynolds, K., and Holsten, E. (1994). Classification of spruce beetle hazard in Lutz spruce (Picea x lutzii) stands on the Kenai Peninsula, Alaska. Can. J. Forest Res. 24, 1015–1021. doi: 10.1139/x94-132
Sakai, A. K., Allendorf, F. W., Holt, J. S., Lodge, D. M., Molofsky, J., With, K. A., et al. (2001). The population biology of invasive species. Annual Rev. Ecol. Systemat. 32, 305–332. doi: 10.1146/annurev.ecolsys.32.081501.114037
Sambaraju, K. R., Carroll, A. L., Zhu, J., Stahl, K., Moore, R. D., and Aukema, B. H. (2012). Climate change could alter the distribution of mountain pine beetle outbreaks in Western Canada. Ecography 35, 211–223. doi: 10.1111/j.1600-0587.2011.06847.x
Santos, M. J., and Whitham, T. G. (2010). Predictors of Ips confusus outbreaks during a record drought in Southwestern USA: implications for monitoring and management. Environ. Manage. 45, 239–249. doi: 10.1007/s00267-009-9413-6
Sarikaya, O., Karaceylan, I. B., and Sen, I. (2018). maximum entropy modeling (Maxent) of current and future distributions of Ips Mannsfeldi (Wachtl, 1879) (Curculionidae: Scolytinae) in Turkey. Appl. Ecol. Environm. Res. 16, 2527–2535. doi: 10.15666/aeer/1603_25272535
Schelhaas, M.-J., Nabuurs, G.-J., and Schuck, A. (2003). Natural disturbances in the European forests in the 19th and 20th centuries. Global Change Biol. 9, 1620–1633. doi: 10.1046/j.1365-2486.2003.00684.x
Scheller, R. M., Kretchun, A. M., Loudermilk, E. L., Hurteau, M. D., Weisberg, P. J., and Skinner, C. (2018). Interactions among fuel management, species composition, bark beetles, and climate change and the potential effects on forests of the Lake Tahoe Basin. Ecosystems 21, 643–656. doi: 10.1007/s10021-017-0175-3
Seidl, R., Baier, P., Rammer, W., Schopf, A., and Lexer, M. J. (2007). Modelling tree mortality by bark beetle infestation in Norway spruce forests. Ecol. Modell. 206, 383–399. doi: 10.1016/j.ecolmodel.2007.04.002
Seidl, R., Müller, J., Hothorn, T., Bässler, C., Heurich, M., and Kautz, M. (2016). Small beetle, large-scale drivers: how regional and landscape factors affect outbreaks of the European spruce bark beetle. J. Appl. Ecol. 53, 530–540. doi: 10.1111/1365-2664.12540
Seidl, R., and Rammer, W. (2017). Climate change amplifies the interactions between wind and bark beetle disturbances in forest landscapes. Landsc. Ecol. 32, 1485–1498. doi: 10.1007/s10980-016-0396-4
Seidl, R., Rammer, W., Jäger, D., and Lexer, M. J. (2008). Impact of bark beetle (Ips typographus L.) disturbance on timber production and carbon sequestration in different management strategies under climate change. For. Ecol. Manage. 256, 209–220.
Seidl, R., Schelhaas, M.-J., Lindner, M., and Lexer, M. J. (2009). Modelling bark beetle disturbances in a large scale forest scenario model to assess climate change impacts and evaluate adaptive management strategies. Regional Environm. Change 9, 101–119. doi: 10.1007/s10113-008-0068-2
Seidl, R., Schelhaas, M.-J., Rammer, W., and Verkerk, P. J. (2014). Increasing forest disturbances in Europe and their impact on carbon storage. Nat. Climate Change 4, 806–810. doi: 10.1038/nclimate2318
Seidl, R., Thom, D., Kautz, M., Martin-Benito, D., Peltoniemi, M., Vacchiano, G., et al. (2017). Forest disturbances under climate change. Nat. Climate Change 7, 395–402. doi: 10.1038/nclimate3303
Sharov, A. A., and Liebhold, A. M. (1998). Model of slowing the spread of gypsy moth (lepidoptera: Lymantriidae) with a barrier zone. Ecol. Appl. 8, 1170–1179. doi: 10.1890/1051-0761(1998)008[1170:MOSTSO]2.0.CO;2
Sidder, A. M., Kumar, S., Laituri, M., and Sibold, J. S. (2016). Using spatiotemporal correlative niche models for evaluating the effects of climate change on mountain pine beetle. Ecosphere 7:e01396. doi: 10.1002/ecs2.1396
Siegert, N. W., McCullough, D. G., Liebhold, A. M., and Telewski, F. W. (2014). Dendrochronological reconstruction of the epicentre and early spread of emerald ash borer in North America. Divers. Distribut. 20, 847–858. doi: 10.1111/ddi.12212
Sims, C., Aadland, D., and Finnoff, D. (2010). A dynamic bioeconomic analysis of mountain pine beetle epidemics. J. Econ. Dynam. Cont. 34, 2407–2419. doi: 10.1016/j.jedc.2010.06.010
Sivrikaya, F., Ozcan, G. E., Enez, K., and Sakici, O. E. (2022). Comparative study of the analytical hierarchy process, frequency ratio, and logistic regression models for predicting the susceptibility to Ips sexdentatus in crimean pine forests. Ecol. Inform. 71:101811. doi: 10.1016/j.ecoinf.2022.101811
Skarpaas, O., and Økland, B. (2009). Timber import and the risk of forest pest introductions. J. Appl. Ecol. 46, 55–63. doi: 10.1111/j.1365-2664.2008.01561.x
Smith, S. E., Mendoza, M. G., Zuniga, G., Kalbrook, K., Hayes, J. L., and Byrne, D. N. (2013). Predicting the distribution of a novel bark beetle and its pine hosts under future climate conditions. Agricult. Forest Entomol. 15, 212–226. doi: 10.1111/afe.12007
Sobek-Swant, S., Kluza, D. A., Cuddington, K., and Lyons, D. B. (2012). Potential distribution of emerald ash borer: what can we learn from ecological niche models using Maxent and GARP? For. Ecol. Manage. 281, 23–31. doi: 10.1016/j.foreco.2012.06.017
Sommerfeld, A., Rammer, W., Heurich, M., Hilmers, T., Müller, J., and Seidl, R. (2021). Do bark beetle outbreaks amplify or dampen future bark beetle disturbances in central europe? J. Ecol. 109, 737–749. doi: 10.1111/1365-2745.13502
Soukhovolsky, V., Kovalev, A., Tarasova, O., Modlinger, R., Křenová, Z., Mezei, P., et al. (2022). Wind damage and temperature effect on tree mortality caused by Ips typographus L.: phase transition model. Forests 13:180. doi: 10.3390/f13020180
Stadelmann, G., Bugmann, H., Wermelinger, B., and Bigler, C. (2014). Spatial interactions between storm damage and subsequent infestations by the European spruce bark beetle. For. Ecol. Manage. 318, 167–174. doi: 10.1016/j.foreco.2014.01.022
Stadelmann, G., Bugmann, H., Wermelinger, B., Meier, F., and Bigler, C. (2013). A predictive framework to assess spatio-temporal variability of infestations by the European spruce bark beetle. Ecography 36, 1208–1217. doi: 10.1111/j.1600-0587.2013.00177.x
Straw, N. A., Williams, D. T., Kulinich, O., and Gninenko, Y. I. (2013). Distribution, impact and rate of spread of emerald ash borer Agrilus planipennis (Coleoptera: Buprestidae) in the Moscow region of Russia. Forestry: Int. J. Forest Res. 86, 515–522. doi: 10.1093/forestry/cpt031
Strohm, S., Reid, M. L., and Tyson, R. C. (2016). Impacts of management on Mountain Pine Beetle spread and damage: a process-rich model. Ecol. Modell. 337, 241–252. doi: 10.1016/j.ecolmodel.2016.07.010
Subramanian, N., Bergh, J., Johansson, U., Nilsson, U., and Sallnaes, O. (2016). Adaptation of forest management regimes in southern sweden to increased risks associated with climate change. Forests 7:8. doi: 10.3390/f7010008
Swallow, B., Birrell, P., Blake, J., Burgman, M., Challenor, P., Coffeng, L. E., et al. (2022). Challenges in estimation, uncertainty quantification and elicitation for pandemic modelling. Epidemics 38:100547. doi: 10.1016/j.epidem.2022.100547
Taylor, R. A. J., Bauer, L. S., Poland, T. M., and Windell, K. N. (2010). Flight performance of Agrilus planipennis (Coleoptera: Buprestidae) on a flight mill and in free flight. J. Insect Behav. 23, 128–148. doi: 10.1007/s10905-010-9202-3
Temperli, C., Bugmann, H., and Elkin, C. (2013). Cross-scale interactions among bark beetles, climate change, and wind disturbances: a landscape modeling approach. Ecol. Monographs 83, 383–402. doi: 10.1890/12-1503.1
Temperli, C., Veblen, T. T., Hart, S. J., Kulakowski, D., and Tepley, A. J. (2015). Interactions among spruce beetle disturbance, climate change and forest dynamics captured by a forest landscape model. Ecosphere 6:art231. doi: 10.1890/ES15-00394.1
Thiele, J. C., Nuske, R. S., Ahrends, B., Panferov, O., Albert, M., Staupendahl, K., et al. (2017). Climate change impact assessment–A simulation experiment with Norway spruce for a forest district in Central Europe. Ecol. Modell. 346, 30–47. doi: 10.1016/j.ecolmodel.2016.11.013
Thrippleton, T., Temperli, C., Krumm, F., Mey, R., Zell, J., Stroheker, S., et al. (2023). Balancing disturbance risk and ecosystem service provisioning in swiss mountain forests: an increasing challenge under climate change. Region. Environm. Change 23:29. doi: 10.1007/s10113-022-02015-w
Tobin, P. C., Diss-Torrance, A., Blackburn, L. M., and Brown, B. D. (2010). What does “local” firewood buy you? Managing the risk of invasive species introduction. J. Econ. Entomol. 103, 1569–1576. doi: 10.1603/EC10140
Trân, J. K., Ylioja, T., Billings, R. F., Régnièère, J., and Ayres, M. P. (2007). Impact of minimum winter temperatures on the population dynamics of dendroctonus frontalis. Ecol. Appl. 17, 882–899. doi: 10.1890/06-0512
Trubin, A., Kozhoridze, G., Zabihi, K., Modlinger, R., Singh, V. V., Surovy, P., et al. (2023). Detection of susceptible norway spruce to bark beetle attack using planetscope multispectral imagery. Front. Forests Global Change 6:1130721. doi: 10.3389/ffgc.2023.1130721
Trzcinski, M. K., and Reid, M. L. (2009). Intrinsic and extrinsic determinants of mountain pine beetle population growth. Agricult. Forest Entomol. 11, 185–196. doi: 10.1111/j.1461-9563.2008.00408.x
Tuomola, J., Gruffudd, H., Ruosteenoja, K., and Hannunen, S. (2021). Could pine wood nematode (Bursaphelenchus xylophilus) cause pine wilt disease or even establish inside healthy trees in finland now-or ever? Forests 12:1679. doi: 10.3390/f12121679
Turchin, P., Lorio, P. L., Taylor, A. D., and Billings, R. F. (1991). Why do populations of southern pine beetles (Coleoptera: Scolytidae) fluctuate? Environ. Entomol. 20, 401–409. doi: 10.1093/ee/20.2.401
Turchin, P., and Thoeny, W. T. (1993). Quantifying dispersal of southern pine beetles with mark-recapture experiments and a diffusion model. Ecol. Appl. 3, 187–198. doi: 10.2307/1941801
Ungerer, M. J., Ayres, M. P., and Lombardero, M. J. (1999). Climate and the northern distribution limits of Dendroctonus frontalis Zimmermann (Coleoptera: Scolytidae). J. Biogeography 26, 1133–1145. doi: 10.1046/j.1365-2699.1999.00363.x
Václavík, T., and Meentemeyer, R. K. (2012). Equilibrium or not? Modelling potential distribution of invasive species in different stages of invasion. Divers. Distribut. 18, 73–83. doi: 10.1111/j.1472-4642.2011.00854.x
van Halder, I., Sacristan, A., Martin-Garcia, J., Pajares, J. A., and Jactel, H. (2022). Pinus pinea: a natural barrier for the insect vector of the pine wood nematode? Ann. For. Sci. 79:1. doi: 10.1186/s13595-022-01159-3
Vannatta, A. R., Hauer, R. H., and Schuettpelz, N. M. (2012). Economic analysis of emerald ash borer (Coleoptera: Buprestidae) management options. J. Econ. Entomol. 105, 196–206. doi: 10.1603/EC11130
Vernon, I., Owen, J., Aylett-Bullock, J., Cuesta-Lazaro, C., Frawley, J., Quera-Bofarull, A., et al. (2022). Bayesian emulation and history matching of june. Phil. Trans. R. Soc. A 380:20220039. doi: 10.1098/rsta.2022.0039
Vindstad, O. P. L., Jepsen, J. U., Ek, M., Pepi, A., and Ims, R. A. (2019). Can novel pest outbreaks drive ecosystem transitions in northern-boreal birch forest? J. Ecol. 107, 1141–1153. doi: 10.1111/1365-2745.13093
Vorster, A. G., Evangelista, P. H., Stohlgren, T. J., Kumar, S., Rhoades, C. C., Hubbard, R. M., et al. (2017). Severity of a mountain pine beetle outbreak across a range of stand conditions in Fraser Experimental Forest, Colorado, United States. For. Ecol. Manage. 389:116–126. doi: 10.1016/j.foreco.2016.12.021
Wang, J., Deng, J., Yan, W., and Zheng, Y. (2023). Habitat suitability of pine wilt disease in northeast china under climate change scenario. Forests 14:1687. doi: 10.3390/f14081687
Waring, K. M., Reboletti, D. M., Mork, L. A., Huang, C.-H., Hofstetter, R. W., Garcia, A. M., et al. (2009). Modeling the impacts of two bark beetle species under a warming climate in the Southwestern USA: ecological and economic consequences. Environ. Manage. 44, 824–835. doi: 10.1007/s00267-009-9342-4
Watson, A. C., Sullivan, J., Amacher, G. S., and Asaro, C. (2013). Cost sharing for pre-commercial thinning in southern pine plantations: willingness to participate in Virginia's pine bark beetle prevention program. Forest Policy Econ. 34, 65–72. doi: 10.1016/j.forpol.2013.05.004
Watt, M. S., Ganley, R. J., Kriticos, D. J., and Manning, L. K. (2011). Dothistroma needle blight and pitch canker: the current and future potential distribution of two important diseases of Pinus species. Can. J. Forest Res. 41, 412–424. doi: 10.1139/X10-204
Webb, C. R., Mona, T., and Gilligan, C. A. (2021). Predicting the potential for spread of emerald ash borer (Agrilus planipennis) in Great Britain: what can we learn from other affected areas? Plants People Planet 3, 402–413. doi: 10.1002/ppp3.10195
Weed, A. S., Ayres, M. P., Liebhold, A. M., and Billings, R. F. (2017). Spatio-temporal dynamics of a tree-killing beetle and its predator. Ecography 40, 221–234. doi: 10.1111/ecog.02046
Weed, A. S., Bentz, B. J., Ayres, M. P., and Holmes, T. P. (2015). Geographically variable response of Dendroctonus ponderosae to winter warming in the western United States. Landsc. Ecol. 30, 1075–1093. doi: 10.1007/s10980-015-0170-z
Weiss, J., Sheffield, K., Weeks, A., and Smith, D. (2019). Modelling the incursion and spread of a forestry pest: Case study of monochamus alternatus hope (coleoptera: Cerambycidae) in victoria. Forests 10:198. doi: 10.3390/f10020198
Wermelinger, B. (2004). Ecology and management of the spruce bark beetle Ips typographus–a review of recent research. Forest Ecol. Managem. 202, 67–82. doi: 10.1016/j.foreco.2004.07.018
Wermelinger, B., Epper, C., Kenis, M., Ghosh, S., and Holdenrieder, O. (2012). Emergence patterns of univoltine and bivoltine Ips typographus (L.) populations and associated natural enemies. J. Appl. Entomol. 136, 212–224. doi: 10.1111/j.1439-0418.2011.01629.x
Wermelinger, B., and Seifert, M. (1998). Analysis of the temperature dependent development of the spruce bark beetle Ips typographus (L) (Col., Scolytidae). J. Appl. Entomol. 122, 185–191. doi: 10.1111/j.1439-0418.1998.tb01482.x
Wermelinger, B., and Seifert, M. (1999). Temperature-dependent reproduction of the spruce bark beetle Ips typographus, and analysis of the potential population growth. Ecol. Entomol. 24, 103–110. doi: 10.1046/j.1365-2311.1999.00175.x
White, P., and Powell, J. (1997). Phase transition from environmental to dynamic determinism in mountain pine beetle attack. Bull. Math. Biol. 59, 609–643. doi: 10.1016/S0092-8240(97)00025-6
Wildemeersch, M., Franklin, O., Seidl, R., Rogelj, J., Moorthy, I., and Thurner, S. (2019). Modelling the multi-scaled nature of pest outbreaks. Ecol. Modell. 409:108745. doi: 10.1016/j.ecolmodel.2019.108745
Williams, D. T., Cull, T., and Forster, J. (2021). Investigating the abundance and flight period of bark beetles (coleoptera: Curculionidae: Scolytinae) over elevational gradients in sitka spruce forests. Agric. For. Entomol. 23, 134–145. doi: 10.1111/afe.12412
Williams, D. W., and Liebhold, A. M. (2002). Climate change and the outbreak ranges of two North American bark beetles. Agricult. Forest Entomol. 4, 87–99. doi: 10.1046/j.1461-9563.2002.00124.x
Withrow, J. R., Smith, E. L., Koch, F. H., and Yemshanov, D. (2015). Managing outbreaks of invasive species – A new method to prioritize preemptive quarantine efforts across large geographic regions. J. Environ. Manage. 150, 367–377. doi: 10.1016/j.jenvman.2014.11.001
Woodward, B. D., Evangelista, P. H., and Vorster, A. G. (2018). Mapping progression and severity of a southern colorado spruce beetle outbreak using calibrated image composites. Forests 9:336. doi: 10.3390/f9060336
Wulder, M. A., Ortlepp, S. M., White, J. C., Coops, N. C., and Coggins, S. B. (2009). Monitoring the impacts of mountain pine beetle mitigation. For. Ecol. Manage. 258, 1181–1187. doi: 10.1016/j.foreco.2009.06.008
Wulder, M. A., White, J. C., Bentz, B., Alvarez, M. F., and Coops, N. C. (2006). Estimating the probability of mountain pine beetle red-attack damage. Remote Sens. Environ. 101, 150–166. doi: 10.1016/j.rse.2005.12.010
Xia, C., Chon, T.-S., Takasu, F., Choi, W. I., and Park, Y.-S. (2022). Simulating pine wilt disease dispersal with an individual-based model incorporating individual movement patterns of vector beetles. Front. Plant Sci. 13:886867. doi: 10.3389/fpls.2022.886867
Xiaodong, Y., and Shugart, H. H. (2005). FAREAST: a forest gap model to simulate dynamics and patterns of eastern Eurasian forests. J. Biogeography 32, 1641–1658. doi: 10.1111/j.1365-2699.2005.01293.x
Yemshanov, D., Haight, R. G., Koch, F. H., Lu, B., Venette, R., Lyons, D. B., et al. (2015). Optimal allocation of invasive species surveillance with the maximum expected coverage concept. Divers. Distribut. 21, 1349–1359. doi: 10.1111/ddi.12358
Yemshanov, D., Haight, R. G., MacQuarrie, C. J. K., Koch, F. H., Liu, N., Venette, R., et al. (2020). Optimal planning of multi-day invasive species surveillance campaigns. Ecol. Solut. Evid. 1:12029. doi: 10.1002/2688-8319.12029
Yemshanov, D., Koch, F. H., Lu, B., Lyons, D. B., Prestemon, J. P., Scarr, T., et al. (2014). There is no silver bullet: the value of diversification in planning invasive species surveillance. Ecol. Econ. 104, 61–72. doi: 10.1016/j.ecolecon.2014.04.024
Yemshanov, D., Koch, F. H., Lyons, D. B., Ducey, M., and Koehler, K. (2012). A dominance-based approach to map risks of ecological invasions in the presence of severe uncertainty. Divers. Distribut. 18, 33–46. doi: 10.1111/j.1472-4642.2011.00848.x
Yoon, S., Jung, J.-M., Hwang, J., Park, Y., and Lee, W.-H. (2023). Ensemble evaluation of the spatial distribution of pine wilt disease mediated by insect vectors in South Korea. For. Ecol. Manage. 529:120677. doi: 10.1016/j.foreco.2022.120677
Yu, L., Zhan, Z., Ren, L., Liu, Y., Huang, H., and Luo, Y. (2022a). Effects of stand and landscape level variables on shoot damage ratios caused by shoot beetles in Southwest China. For. Ecol. Manage. 507:120030. doi: 10.1016/j.foreco.2022.120030
Yu, L., Zhan, Z., Zhou, Q., Gao, B., Ren, L., Huang, H., et al. (2022b). Climate drivers of pine shoot beetle outbreak dynamics in Southwest China. Remote Sens. 14:2728. doi: 10.3390/rs14122728
Yurk, B. P., and Powell, J. A. (2010). Modeling the effects of developmental variation on insect phenology. Bull. Math. Biol. 72, 1334–1360. doi: 10.1007/s11538-009-9494-7
Zhan, Z., Yu, L., Ren, L., Gao, B., Li, H., Wang, L., et al. (2023). The association between stand and landscape level factors and red turpentine beetle damage in different infestation stages. For. Ecol. Manage. 531:120790. doi: 10.1016/j.foreco.2023.120790
Zhang, K., Hu, B., and Robinson, J. (2014). Early detection of emerald ash borer infestation using multisourced data: a case study in the town of Oakville, Ontario, Canada. J. Appl. Remote Sens. 8:083602. doi: 10.1117/1.JRS.8.083602
Zhou, Y., Ge, X., Zou, Y., Guo, S., Wang, T., and Zong, S. (2019). Climate change impacts on the potential distribution and range shift of Dendroctonus ponderosae (Coleoptera: Scolytidae). Forests 10:860. doi: 10.3390/f10100860
Zhou, Y., Ge, X., Zou, Y., Guo, S., Wang, T., and Zong, S. (2021). Prediction of the potential global distribution of the Asian Long horned beetle Anoplophora glabripennis (Coleoptera: Cerambycidae) under climate change. Agricult. Forest Entomol. 23, 557–568. doi: 10.1111/afe.12461
Keywords: ash, economic impacts, Emerald ash borer, environmental drivers, European spruce beetle, forest pest modeling, pine, spruce
Citation: Keenan VA, Kanyamibwa M, Marion G, Touza J and Kleczkowski A (2026) A review of process modeling for wood borer pests. Front. For. Glob. Change 9:1706950. doi: 10.3389/ffgc.2026.1706950
Received: 16 September 2025; Revised: 17 December 2025;
Accepted: 08 January 2026; Published: 04 February 2026.
Edited by:
Milica Zlatkovic, University of Novi Sad, SerbiaReviewed by:
Stephen Joseph Mayor, Ontario Forest Research Institute, CanadaCerian Webb, University of Cambridge, United Kingdom
Copyright © 2026 Keenan, Kanyamibwa, Marion, Touza and Kleczkowski. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Vincent A. Keenan, dmluY2VudGtlZW5hbkBob3RtYWlsLmNvLnVr; Julia Touza, anVsaWEudG91emFAeW9yay5hYy51aw==
Vincent A. Keenan1,2,3*