Ef fects of meteorological factors on the defoliation dynamics of the larch caterpillar ( Dendrolimus superans Butler)in the Great Xing’an boreal forests

Lei Fang · Yue Yu · Guofei Fang · Xu Zhang ·Zhijun Yu · Xudong Zhang · Ellen Crocker ·Jian Yang

Abstract Larch caterpillars are widely distributed in the Great Xing’an boreal forests; however, the relationship between caterpillar defoliation dynamics and climatic factors is poorly understood. The aims of this study are to investigate the primary weather conditions that might inf luence forest defoliation and to identify the most important life stage of the larch caterpillar at which forest defoliation might be mitigated by incorporating more inhibitory inf luences from climatic factors. The life cycle of the larch caterpillar was partitioned into four stages and multiple linear regression and mixed ef fect models were combined with a relative weight analysis approach to evaluate the importance and inf luence of meteorological variables on defoliation dynamics. The results show that warmer temperatures in growing seasons and overwintering periods can increase the defoliation area, while rainy and humid growing seasons decrease the defoliation area. Total precipitation during the early instar larval period had the greatest power to explain the variance in defoliation dynamics and had a very strong inhibitory ef fect, followed by the accumulative temperatures of the late instar larval period which had a positive impact,and precipitation during the middle instar larval period which had a negative impact. Weather conditions during the early instar larval period had the greatest inf luence on the area defoliated and accounted for 40% of the explained variance. This study demonstrates that climatic warming and drying will increase the risk of larch caterpillar outbreaks in the Great Xing’an Mountains.

Keywords Larch caterpillar · Defoliation area · Climatic warming · Life cycle · Great xing’an mountains · Relative importance

Introduction

Forest insects are critical destructive agents that can directly exert stress ef fects on forest growth or even lead to substantial tree mortality, thereby decreasing timber production and ecosystem services related to water and soil retention and carbon sinks (Dale et al. 2001; Kurz et al. 2008; Medvigy et al. 2012; Zhang et al. 2014; Flower and Gonzalez-Meler 2015). As a natural disturbance, forest insects may indirectly interact with other disturbances (e.g., wildf ire, drought and pathogens), and amplify their negative ef fects on ecosystem resilience, which may trigger irreparable ecosystem shifts (Kane et al. 2017; Seidl et al. 2017). Monitoring and predicting outbreaks of forest insects are topics that have received increasing attention due to the benef its produced by the implementation of practical governance measures in advance to improve forest health and maintain valuable ecological functions.

Climatic changes, especially warming and drought, are important driving factors that regulate insect outbreaks through direct impacts on insect physiology and behavior and indirect impacts on host plants and predators (Liebhold et al. 2000; Bale et al. 2002; Bentz et al. 2010; Kolb et al. 2016). For instance, minimum winter temperatures are one of the predominant environmental obstacles that forest insects must overcome in order to survive (Han et al.2005; Bale and Hayward 2010). Many studies have shown that high-latitude regions will experience stronger warming trends in winter (Dale et al. 2001; Bale and Hayward 2010;Venette 2013), which may signif icantly shorten the duration of winter and decrease the inhibitory ef fects of cold temperatures on forest insects. Under such trends, it is predicted that more insects will successfully survive winter,negatively af fecting forest health. Warming is expected to dramatically inf luence the population dynamics of forest insects; therefore, temperature is frequently regarded as an important predictor of forest insect outbreaks (Yamamura et al. 2005; Bentz et al. 2016; Zidon et al. 2016).

Similarly, changes in precipitation patterns and associated water def icits are also identif ied as critical factors controlling the breeding and growth of insects. Heavy rainfall or hailstorms during the breeding period can mechanically injure or kill eggs that adhere to foliage or twigs. Cold and rainy weather can inhibit the dispersal and feeding of newly hatched larvae and slow their growth (McManus et al. 1980;Rouault et al. 2006; Jaworski and Hilszczański 2013). Moist environments may induce the dispersal of diseases caused by bacteria, fungi, or viruses, which may inf luence the population and growth of forest insects (McManus et al. 1980; Han et al. 2004). Moreover, physiological stress caused by warming and/or water def icits that increase host plant susceptibility to forest insects are also well-documented globally (Bale et al. 2002; Rouault et al. 2006; van Mantgem et al. 2009;Allen et al. 2010; Kolb et al. 2016).

Defoliating insects such as the pine caterpillar, gypsy moth and spruce budworm are the major insects causing disturbance across Eurasian and North American boreal forests, and these groups have received increasing public attention as a result of their pervasive inf luence and the weak resistance of boreal ecosystems to insect disturbance (Gray 2008; Bouchard and Pothier 2010; Weed et al. 2013; Zhang et al. 2014). Previous studies have shown that the population dynamics of forest insects are largely driven by the spatial synchrony with extreme climatic events (Han et al. 2005;Wang et al. 2007; Hódar et al. 2012) and by phenological asynchrony with host plants or predators (Liebhold et al.2000; Bjornstad et al. 2010; Pureswaran et al. 2015). These studies have provided notable conclusions that have helped to disentangle the mechanisms underlying how environmental factors regulate insect population dynamics. Most of those discussions have focused on forest insects that have relatively short life histories, while the inf luences of climatic factors on the population dynamics of forest insects with longer life histories, i.e., insects that require several years to complete a whole lifecycle, have received relatively little attention.

Forest insects in boreal regions usually require multiple years to develop body growth and physiological functions(Bale et al. 2002). This multiyear life history increases the complexity of research endeavors attempting to clarify the climate-driven outbreak mechanisms of these forest insects. Previous studies have emphasized the linkages between inter- or intra-annual temperature variations and population dynamics of forest insects, while much less is known on other climatic factors that may also exert similar or consistent inf luences on the insect population at dif ferent developmental stages. Since forest insects may develop various adaption strategies to overcome severe exogenous environments (Bale et al. 2002; Pimentel et al. 2010), clarifying the linkages between weather conditions and population dynamics may be helpful for forming a comprehensive understanding about regulatory mechanisms that occur at the regional scale. Specif ically, identifying an insect’s sensitive life-history stage in terms of local weather conditions and quantifying the relative importance of those variables on population dynamics are of notable practical signif icance for forest insect management.

Compared with boreal ecosystems in other regions, the forests in the Great Xing’an Mountains, which serve as an important ecological barrier to extreme climate (e.g., the Siberian cold current) and natural disasters (e.g., f loods) in northeastern China, may be particularly vulnerable to climate change ef fects because they are located within a border zone between boreal and temperate climates (Ni 2011). The primary type of forest insect disturbance that has caused serious forest losses is the larch caterpillar,Dendrolimus superansButler, which requires two years to complete its lifecycle in this region (Yuan et al. 2008; Chen et al. 2010).An extreme outbreak of larch caterpillar occurred from 1989 to 1991, which caused over one million ha of forest damage in the Great Xing’an Mountains (Liu et al. 1994). The mean population density of insects was approximately 283 individuals per tree, and the maximum might have reached 2000 per tree (Liu et al. 2002). Thus, the objectives of this study were: (1) to identify the primary weather conditions that govern f luctuations in larch caterpillar populations and are linked to its special life-history traits; and, (2) to determine the critical life stage in which the insect population would be sensitive to the inf luence of climatic factors. To achieve these goals, historical larch caterpillar records and daily meteorological observations from six forestry bureaus in the Great Xing’an boreal forests were used to examine these relationships. We evaluated the relative inf luence of variables related to temperature, relative humidity and precipitation, which were calculated from daily observations of four life stages of the larch caterpillar, on episodic outbreaks of this defoliator based on a linear regression model and linear mixed ef fect model.

Materials and methods

Study site

The study area is located within six forestry bureaus (local administrative agencies) of the Great Xing’an Mountains in Northeast China (121° 10′ 41″–127° 01′ 06″ E, 50° 05′53″–53° 33′ 38″ N), Jiagedaqi, Songling, Amuer, Xinlin,Hanjiayuan and Shibazhan, which occupy more than half of the total area and are administered by Heilongjiang Province (Fig. 1). The topography in the Great Xing’an Mountains is gently undulating with an average slope of 10°. The elevation ranges from 130 m in the northeast to 1500 m in the southwest. The climate is mid-latitude continental cold-temperate, the mean annual temperature is − 1.8 °C, and the monthly mean January and July temperatures are − 20 °C and 30 °C, respectively. The mean annual precipitation is approximately 550 mm and presents great inter-annual variations, with more than 60% falling between May and October. The frost-free period is less than 130 days in the south and less than 100 days in the north. The forests in the study area are classif ied as cooltemperate coniferous forests and form the southern extension of the eastern Siberian boreal forests. Dahurian larch(Larix gmeliniiRupr.) and white birch (Betula platyphyllaSukaczev) are the dominant canopy species, mixed with Scotch pine (Pinus sylvestrisvar. mongolicaLitv.), Korean spruce (Picea koraiensisNakai), Mongolian oak (Quercus mongolicaFisch. ex Ledeb.), and two species of aspen(Populus davidianaDode andPopulus suaveolensFisch.)(Fang et al. 2015). As a result of heavy forest harvesting since the middle of the last century and frequent wildf ires,the forest landscape of the study area is primarily composed of middle and young age stands.

Life cycle of the larch caterpillar

The larch caterpillar,D. superans, is an endemic species in China that can seriously endanger coniferous species, such as Dahurian larch (its preferred host plant), Scotch pine and Korean spruce (Liu et al. 1994; Yu et al. 2016). Generally the insect produces one generation per year in the southern part of northeastern China, although it has one generation every two years in our study area due to the relatively short period of appropriate temperatures for development occur(Liu et al. 1994; Tian et al. 1998). Previous studies have found that the larch caterpillar has ninth instar larvae and begins to overwinter during the second to fourth instar stages at the end of September (Fig. 2). The caterpillars then leave the host tree and remain in ground litter or grasses from October to the following April. Low temperatures during winter can freeze the larvae, which then thaw when temperatures rebound. At the end of April of the following year,the caterpillars will ascend the tree and consume the foliage until the second overwintering period, which begins during the f ifth to seventh instar stages (Yue et al. 1996; Chen et al. 2010). After entering its second diapause period, the surviving larvae will continue their development and pupate at the end of May. After a 20-day pupal stage, adults lay eggs on the needles during the second half of June. The second generation of larch caterpillar eggs will hatch during June and July and eventually change to the larval stage in approximately 15 days. By reviewing its life cycle, we can identify key temporal nodes at which climatic factors may regulate the population by directly impacting their physiology or by modifying the phenology of their host plants, thus indirectly impacting the larch caterpillars.

Data source

Larch caterpillars can impact substantially forest growth in a short time, although these are usually nonlethal to host tree species; however, when the degree of defoliation is over 75%, tree mortality may occur (Yu and Yin 1980). With regards to the low mortality induced by the larch caterpillar,the defoliation area is considered a more reliable indicator to ref lect the adverse consequences of larch caterpillars from a forest health perspective than tree mortality. In this study, the area of defoliated forests was used as a proxy to ref lect the dynamics of the larch caterpillar population density. Historical defoliation data between 2003 and 2017 were obtained from annual reports of the six local forest pest control and quarantine stations, which are completed and maintained by the these forestry bureaus (Fig. 3). Thirty-four valid records of annual defoliation areas of larch caterpillar were obtained from these six forestry bureaus. These records were double-checked by a higher administrative department as the f igures of defoliation area related to f inancial inputs for pest control. Years without valid records were considered those without pronounced outbreak events.

Fig. 2 Life cycle of the larch caterpillar (two years per generation) in the Great Xing’an boreal forests

Fig. 3 Annual areas of defoliation for 2003–2017 for six forestry bureaus in the Great Xing’an boreal forests

The estimation of defoliation area was based on f ield surveys carried out along the predesigned routes that covered the stands where larch caterpillar outbreaks occurred as well as new 100 m × 100 m sampling plots. Populations of larch caterpillars were checked using the shaking-of f collection method for 20 larch trees within each plot. A forest sub-compartment was classif ied as a defoliated area when the population density of the insect was greater than 20 caterpillars per tree according to a national guideline of defoliation severity classif ication. Twenty caterpillars per tree is an empirical population threshold for reporting pine caterpillar disaster,and a lower population density does not warrant a report to higher levels of the forest insect pest control department(see http://www.fores tpest.org/forec astin g/forec astme thods/). The defoliation area was estimated based on the proportion of defoliated sub compartments. Although this estimation was somewhat coarse, such defoliation data are the most reliable of ficial records of larch caterpillar outbreaks in the study area.

Historical weather data for the f ive forestry bureaus were obtained from the Climatic Data Center, National Meteorological Information Center, and the China Meteorological Administration ( http://data.cma.cn/). These weather stations are part of a network of meteorological stations across China. Daily statistics of weather variables, such as average temperature, precipitation, and relative humidity were obtained from 1975 to 2017 (Fig. 4). The historical weather observations show increasing average annual temperatures and decreasing relative humidities and precipitation, which suggest a warmer and drier climate in the Great Xiang’an Mountains.

By reviewing the biological characteristics of the larch caterpillar, previous studies have revealed that larvae at the seventh to ninth instar stages consume more than 97% of the total leaf amount during the entire life cycle (Yu and Yin 1980; Liu et al. 1994; Yue et al. 1996). Because larch caterpillar f ield surveys were mainly carried out between May and June, middle and late instar larvae coexisted in the forests. It was hypothesized that the defoliation area during a specif ic year was primarily caused by late instar larvae. The life cycle of the late instar larvae was followed, and based on a review of previous studies, 15 weather variables were selected that were considered to impact survival rates during key life stages (Table 1). These included the breeding period,the early instar larval period, the middle instar larval period and the late instar larval period. The 15 weather variables were grouped into four categories based on their temporal information that closely relates each variable to the four key life stages. At the same time, these variables were also separated into three categories, as they ref lect the corresponding environmental conditions of temperature, precipitation and relative humidity.

Statistical analyses

A multiple linear regression (MLR) model was used to analyze the relative importance of weather conditions on the forest area defoliated by larch caterpillars. To avoid overf itting as a result of the limited dataset, we used a maximum of three explanatory variables to establish the models. Because explanatory variables selected from the same weather category often represent high collinearity, one meteorological variable was chosen from each of the three weather categories (i.e., temperature, precipitation, relative humidity) to establish a candidate model. To reduce redundant candidate models, we f irst used the Pearson correlation coef ficient(r) generated by the bivariate correlation test to exclude explanatory variables that have very low linear correlations (r< 0.30) with the response variable. By visualizing the relationships between variables, exponential transformation was carried out for accumulated daily temperatures below the supercooling point (SCP, − 22 °C) during the f irst overwintering period (EATO) and accumulated daily temperatures below the SCP during the second overwintering period (MATO) to ensure that these two variables were linearly correlated with the response variable. Cook’s distance was calculated for each correlation test to eliminate outliers that have Cook’s distance greater than three times the means.A stepwise regression approach was also used as a supplement to choose a candidate model because some predictor variables may not pass the bivariate correlation test but can still explain a high portion of the remaining variation when interacting with the most dominant explanatory variable.

Fig. 4 Trends of a annual mean temperature, b relative humidity,and c precipitation represent a warmer and drier climate in the Great Xing’an Mountains since 1975, based on f ive weather stations

Although such modeling strategy could reduce collinearity among explanatory variables, a variance inf lation factor (VIF) was still used to detect the collinearity between explanatory variables following the protocol of Zuur et al.( 2010). A general rule of thumb is if the VIF is greater than 10, then collinearity is considered high, and if greater than 4,then collinearity occurs. To ensure that signif icant collinearity problems did not occur, we simply excluded the models if their explanatory variables had VIF values higher than 4. The VIF values were calculated using the “car” package in the R statistical software program (R Core Development Team).

The coef ficient of determination (R 2 ) is the most popular quantif ication of model performance and thus is frequently used as a basis to evaluate the relative importance of explanatory variables. However, this measure may be problematic when independently quantifying the unique contribution of an explanatory variable to the model f itting because of the collinearity among explanatory variables. Although a stringent process was carried out to mitigate collinearity,there was still a need to discriminate the direct ef fects from the joint ef fects with other variables when partitioning the explained variance. The proportional contribution of each explanatory variable to the coef ficient of determination was calculated based on the relative weight analysis approach proposed by Johnson ( 2000). The “relaimpo” package was used to calculate the relative weights for each variable contained in the regression model (Grömping 2006). To evaluate the individual strengths of three dif ferent categories of meteorological variables on modeling the defoliation area of larch caterpillars, the accumulated relative weights of the explanatory variables was calculated. It was assumed that meteorological variables that had higher accumulated values of relative weight would provide more explanatory power. Similarly, the accumulated relative weights for each stage of the life cycle were also compared to demonstrate the ef fect of weather conditions on specif ic life stages to identify the factor with the greatest importance for determining the defoliation area of the larch caterpillar.

The inf luence of the meteorological variables on the defoliation area of larch caterpillars was evaluated based on the values and directions (i.e., positive or negative) of the estimated slope coef ficients. Since the defoliation area may be related to other biotic or abiotic conditions belonging to a given forestry bureau, such as the inherent nature of local forests, the topography and human activities associated with pest control, these underlying variations are not well ref lected by the meteorological variables. The inf luence(slopes) of weather conditions on the defoliation area may be biased if the forestry bureau-related variations (intercepts) are assumed to be uniform in the models. The nested structure of our sampling data allows us to add a f lexible random ef fect in the regression model to address individual dif ferences among forestry bureaus. Such mixed ef fects modeling consists of f ixed and random components. Here,the explanatory variables derived from weather observations are def ined as f ixed ef fects and specif ied the forestry bureaus as the random ef fect. We used the “lem4” package to run the linear mixed ef fects (LME) model with random intercepts(Bates et al. 2014).

Table 1 Review of the potential eff ects of weather conditions on larch caterpillars (and some other forest insects) at diff erent periods of the life cycle

Table 1 (continued)

Many studies have noted the importance of standardizing variables into comparable scales for multivariate analysis because the explanatory variables have dif ferent ranges and represent dif ferent meteorological meanings. Thus, the raw data was standardized based on a Z-score scaling method to ensure that their inf luence was measured at comparable scales. A negative Z-score indicates that the raw value is below the mean, and a positive Z-score indicates that the raw value is above the mean. However, current LME functions cannot quantify the statistical signif icance of the f ixed ef fects from explanatory variables in terms of thep-value;therefore, we performed a likelihood ratio (LR) test to verify the signif icance (atalpha= 5% level) of an explanatory variable when introducing this variable into the LME model.Signif icant variables were not selected based solely on the quantif ication of the LR test outputs because the reliability of testing f ixed ef fects is still debatable (Bolker et al. 2009).Akaike information criterion (AIC) was also used to evaluate the signif icance of a variable to the LME model. In general,lower AIC values represent less information loss. Similar to the LR test, we calculated the dif ference in the AIC when applied to the LME model with and without a specif ic variable. A positive AIC dif ference indicates that the variable can reduce information loss. The combined subset of variables selected from these two test methods were identif ied as key regulators to provide a more comprehensive interpretation.

Results

Evaluation of model performance

Nine of f ifteen meteorological variables were well correlated(|r|> 0.30) with the defoliation area of the larch caterpillar(Appendix Fig. S1); they included four temperature variables, three precipitation variables and two relative humidity variables. Together with one model selected by the stepwise regression approach, 25 (i.e., 4 × 3 × 2 + 1) MLR models and corresponding LME models were generated (Table 2). The maximum VIF values of the MLR models were all less than 2, which indicated weak collinearity among the explanatory variables. According to the model outputs, 12 MLR models out of 25 explained 28.3–42.1% of the variance in the annual defoliation area among dif ferent forestry bureaus at thep< 0.05 level. In addition, four MLR models could explain approximately 25% variance in defoliation area at thep< 0.10 level. The MLR model, whose explanatory variables consisted of the accumulated temperature (daily temperature above 8 °C) of the growing season (LATG), total precipitation during the growing season (LTP) and mean daily relative humidity during the growing season (EMRH),explained the maximum variance in the defoliation area of the larch caterpillar, while variable combinations such as the accumulated temperature (daily temperature above 8 °C) of the growing season (MATG)-total precipitation during the growing season (ETP)-EMRH, LATG-total precipitation during the growing season (MTP)-EMRH, LATG-MTPmean daily relative humidity during the growing season(MMRH) and MATG-MTP-EMRH explained over 38.0%of the variance. In addition, the results showed that MLR models with temperature variables associated with the growth period had more explanatory power than models with temperature variables associated with the overwintering period. By incorporating a random ef fect in terms of spatial location, the results showed residual standard errors(RSE) improved for seven of these 25 models, including four of the signif icant models.

Importance analysis

According to the relative importance of explanatory variables in terms of their proportional contribution to the coeff icient of determination, the relative importance of selected variables varied considerably between dif ferent models. The results clearly show that the ETP, which represented the total precipitation of the early instar larvae period, exerted the most relative importance on improving the model explanatory power because it independently contributed approximately 15% (13.5–25.7%) to the explained variance for most of the models (Fig. 5). In addition to the ETP, the LATG(11.1–17.8%), MTP (8.3–16.3%) and EMRH (7.2–15.9%)also independently contributed over 10% explained variance for most of the models when they were included. The two winter temperature variables contributed the least explained variance (EATO: 1.3–7.5%; MATO: 0.5–5.3%) when they were combined with two other kinds of meteorological variables.

The relative importance of explanatory variables were grouped based on their temporal attributes (Fig. 6 a) and physical meanings (Fig. 6 b) and the relative proportion of the total explained variance for each group was calculated.The results indicate that nearly 40% of the total explained variance was contributed by variables related to the early instar larval period, especially the ETP and EMRH(Fig. 6 a). Four meteorological variables belonging to the middle instar larval period together contributed over 31%of the total explained variance, while the LATG, whichwas the only variable of the late instar larval period, contributed over 24% of the total explained variance. In contrast, the only variable selected from the breeding period contributed less than 5% to the total explained variance,indicating the relatively weak importance of weather conditions during this period for determining the defoliation area two years later. These results emphasize the key role of weather conditions during the early and middle instar larval periods in determining the population dynamics of the late instar larvae in future springs. When the relative weights of the explanatory variables were grouped into the three categories of relative humidity, temperature and precipitation (Fig. 6 b), variables related to precipitation contributed 40% of the total explained variance, while variables related to temperature and relative humidity explained approximately 33.6% and 26%, respectively. The precipitation variables were more strongly correlated with the larch caterpillar population density than the other two weather parameters.

Table 2 Parameters evaluating the performance of the multiple linear regression (MLR) model and linear mixed ef fect (LME)model

Fig. 5 Relative importance of nine selected variables based on the proportion of explained variance. See Table 1 for def initions of the variable abbreviations

Fig. 6 Relative importance of the meteorological variables grouped into a four categories based on life stages and b three categories based on meteorological attributes. See Fig. 2 and Table 1 for more information about the variable description at dif ferent life stages

Ef fects of weather conditions on defoliation area

The results of the tests of signif icance based on the LR test and AIC showed that 35 out of 72 slope coef ficients estimated from the LME models were statistically signif icant(Table 2). The directions of ef fects from selected meteorological variables were consistent with the correlation tests.The results show that temperature-related variables had signif icant positive ef fects on the defoliated area (Fig. 7).Higher EATO and MATO indicated relatively warmer winters, and higher MATG and LATG hotter growing seasons.Although the standardized slope coef ficients varied among the models, two cumulative temperature variables of growing seasons, LATG and MATG, clearly exhibited a stronger inf luence on the area of defoliation than the two cumulative low winter temperature variables, EATO and MATO, which can ref lect the degree of coldness in two overwintering periods. Our results also clearly show that variables related to precipitation and relative humidity exhibited strong inhibitory ef fects on the defoliation area of the larch caterpillar,especially the total amount of precipitation during the early instar larvae period (i.e., ETP), which produced the steepest slope in the negative direction. In general, precipitation variables have a stronger negative inf luence than relative humidity variables at the same life stage. Overall, the results demonstrate that a growing season with frequent hot and dry periods will increase the defoliation area.

Fig. 7 Inf luence of the nine selected variables on the defoliation area of larch caterpillar. Only coef ficients passing signif icance tests of the linear mixed ef fects models are shown(variables with the signif icance symbol in Table 2). See Table 1 for def initions of the variable abbreviations

Discussion

The results indicate that weather conditions are closely associated with the area defoliated in the study area. Nine variables exerted a notable inf luence on caterpillar outbreak,and the variables ETP, LATG, MTP and EMRH exerted a considerable inf luence on improving model performance.Weather variables during the early instar larval and breeding stages played key roles in regulating the populations of larch caterpillars. Among these variables, total precipitation at the early instar larval period was the most important variable because it contributed the greatest proportion of the explained variance in the linear models. The coef ficients generated by the mixed ef fects models indicate that the defoliation area could decrease signif icantly with increasing total precipitation during the early and middle instar larval periods, as well as during the breeding period. This f inding suggests that increased precipitation had a very strong inhibitory ef fect on larch caterpillar outbreaks. In contrast,drought (low precipitation) in the growing season can exacerbate defoliation, especially under high temperatures and low relative humidity. Hot droughts during the growing season can increase the suitability of host plants for feeding and growth of insect larvae. Water def icits induced by hot,dry conditions can lead to stomatal closure of leaves, reducing transpiration and maintaining water potential; however,the trees will become warmer and nutrients will be concentrated at the same time, which is favorable to defoliating insects (Fleming and Candau 1998). Previous studies revealed that heavy rainfall or continuous drizzle may inhibit the development of forest insects. Heavy rainfall may kill the eggs or young larvae due to strong mechanical damage or water immersion (Chen et al. 2010; Khaliq et al. 2014).This hypothesis is strongly supported by the negative correlations of precipitation variables observed in this study.

The ef fects of temperature on forest insects have attracted considerable attention in recent decades. Insect disturbances by both defoliators and wood-boring insects are projected to intensify under warmer climates during both the growing season and overwintering period. Our results clearly show that larch caterpillar outbreaks in the study area follow such dynamics. Although a cool climate is usually observed during the growing season in boreal forests, such as in our study area, under a warmer growing season, the feeding activities of insects will increase which will cause more severe forest defoliation. Warmer temperatures will enable larch caterpillars to develop faster and attain higher survival rates than under normal conditions. The notably positive inf luence and explanatory power of temperature in the growing season of the late and middle instar larval periods supports this notion.In addition, the modulation ef fect of winter temperature is another important limitation for larch caterpillar outbreaks.The results clearly show positive but nonlinear relationships between the accumulated daily temperatures below the supercooling point of − 22 °C and the area defoliated.The duration of the supercooling period could dramatically determine the area defoliated, suggesting that warming in winter could facilitate outbreaks of larch caterpillars the following year because more insects survived overwintering. This f inding is in accordance with numerous studies,highlighting the ecological importance of the supercooling period as a predominant agent regulating population dynamics of forest insects (Renault et al. 2002). Among the temperature variables, our results did not identify signif icant relationships between extreme low winter temperatures and area of defoliation, suggesting that the larvae of larch caterpillars during diapause are insensitive or tolerant to ephemeral coldness in winter.

Relative humidity ref lects the saturation status of water vapor in the atmosphere at a given temperature, which is dif ferent from the amount of moisture available as measured by precipitation. At a lower relative humidity, moisture evaporates at a faster rate. The ef fects of relative humidity on forest insect outbreaks have received less attention than the ef fects of temperature and precipitation, primarily because of the high variability of relative humidity and the dif ficulty interpreting its ef fects (Tauber et al. 1998; Jaworski and Hilszczański 2013). Low relative humidity can lead to water def icits in host plants under the combined impact of temperature and precipitation; therefore, relative humidity can indirectly inf luence insect development and survivability(Rouault et al. 2006). Although relative humidity had less overall relative importance than precipitation or temperature for explaining variance in the area of defoliation, relative humidity of the early and middle instar larval periods had consistently negative ef fects on defoliation area relative to precipitation. The EMRH had a stronger ef fect and explanatory power than most of the temperature and precipitation variables. Humid environments, especially those that have almost continuous drizzle, may exacerbate spore dispersal ofBeauveria bassiana,recognized as an important pathogenic bacterium that can trigger a sharp decrease in larch caterpillar populations because the early instar larvae are vulnerable to disease (Xu et al. 2008a; b). Nevertheless,exposure to high or insuf ficient humidity for lengthy periods can generate adverse ef fects on development or survivability of insects because such conditions can directly modulate the interior water balance of the insect body (Tauber et al.1998; Han et al. 2004; Chang and Gao 2008 ) Persistent low humidity can lead to high mortality in insects via dehydration (Rouault et al. 2006). From this perspective, relationships between relative humidity and defoliation area were expected to be nonlinear as a result of adverse ef fects at both extremes. However, we did not think that such nonlinearity needed to be ref lected in our data because the relative humidity of the growing season seldom reaches the critical point at which larch caterpillars are killed by persistent dehydration. Moreover, defoliators can obtain moisture by feeding to maintain water balance (Tauber et al. 1998; Rouault et al. 2006).

Our study provides information for the development of scientif ic pest control strategies to maintain forest health and resilience. The results highlight the important role of weather conditions in modifying the intensity of defoliation caused by larch caterpillars. Furthermore, this research indicates that weather conditions can be linked to specif ic life-history traits of the caterpillar, which may be useful in developing pest control plans. For instance, the results show that weather conditions at the early instar larval stage had the strongest inf luence on the extent of the defoliation.Extreme weather conditions may disrupt the activities of insects and the biological habits upon which insects depend to successfully survive, develop, or mate. Current pest control ef forts can be improved by incorporating the inhibitory ef fects of such weather conditions. For example, the introduction of bacteria that function as pathogens, important biological control agents of the larch caterpillar, can be carried out on days with higher relative humidity during the early instar larval period. Additionally, our results demonstrate the ecological importance of supercooling periods in regulating defoliation. Such natural regulation mechanisms can be applied to improve the ef fectiveness of pest control techniques by destroying the insect’s overwintering shelters (e.g., surface litter) in geographical hotspots of high insect populations during the early instar larval period with extreme cold temperatures.

There are several limitations regarding the use of defoliation data to predict the outbreak of larch caterpillars using empirical models. The f irst concerns the quality of the defoliation data. Although uniform standards of f ield pest investigation are strictly enforced by the national forest service and local pest control stations, these datasets may be subject to bias due to changes in the forestry bureaus, dif ferent f ield investigators, and/or the legacy ef fects of past insect outbreaks and pest control intensity. The other limitation involves the temporal and spatial integrality of the defoliation datasets. Since we only collected some temporally fragmented datasets based on of ficial forest insect reports, our results may represent only a general overview of the relationships and relative importance of meteorological factors involved in regulating larch caterpillar defoliation dynamics.Although the results are reasonable and were mostly crossvalidated by previous studies, it is possible that extrapolating insect outbreaks based on the empirical models generated in our study may have missing defoliation information.However, several recent studies claim that remote sensing can be applied to retrieve continuous forest disturbance history (e.g., insect, f ire and harvesting) at a decadal time scale(Neigh et al. 2014; McDowell et al. 2015; Meigs et al. 2015).Monitoring forest changes related to larch caterpillars over a greater spatial area would be benef icial for obtaining a more comprehensive understanding of the relationships between weather conditions and larch caterpillar outbreaks.

Furthermore, the impacts of forest insects can be inf luenced by biotic factors associated with stand structure and composition (Charbonneau et al. 2012; Robert et al. 2012).Changes in stand structure and composition as a result of forest management operations such as thinning, harvesting and planting may alter the stand’s susceptibility to insects and ultimately cause spatial variability in tree mortality and defoliation severity. Although such impacts work at f iner spatial scales than the one used in this study, we believe that they may induce uncertainty that our models cannot explain.

Conclusions

Meteorological factors are closely associated with the defoliation dynamics of larch caterpillars in the Great Xing’an boreal forests. By considering the special life-history traits of the caterpillar in the study area, the analyses revealed that the meteorological variables calculated from the precipitation at the early instar larval period explained most of the variance in the defoliation area. Precipitation and relative humidity had signif icant negative ef fects on the regulation of defoliation by larch caterpillars. In contrast, temperature variables had signif icant positive ef fects, especially the accumulative temperatures during the late and middle instar larval periods. The duration of supercooling periods during the early and middle instar larval periods can substantially limit the area defoliated. Overall, the study demonstrates that the early instar larval period is the key period during which the population of larch caterpillars may be controlled by exploiting the regulatory ef fects of weather. In this context,the results show that larch caterpillar outbreaks may occur if the early and middle instar larvae experience severe drought during their growing seasons or warm winters during overwintering. These results may provide valuable information to help forestry agencies prevent and control forest insect pests and adapt to novel insect dynamics against the background of global climate change.

Compliance with ethical standardsConf licts of interest The authors declare that they have no conf licts of interest.Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source,provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons.org/licen ses/by/4.0/.