基于两种轨迹模型的褐飞虱迁飞轨迹比较研究

异地预测是迁飞性害虫发生预测的重要内容,迁飞轨迹模拟和预测是能较好地反映害虫迁飞时空动态的一种异地预测方法。褐飞虱作为我国水稻生产上的一种重要迁飞性害虫,其迁飞轨迹的准确预报,可为其灾变预警和有效防控提供科学依据。为了比较选择一些准确性好、分辨率高、易于推广应用的害虫迁飞轨迹模型,选取2006年7月初发生在湖南省洪江市的一次褐飞虱重大北迁过程作为典型个例,运用中尺度天气研究和预报模式WRF,结合NCEP气象再分析数据,利用HYSPLIT和FLEXPART两种轨迹计算模式对褐飞虱迁飞轨迹进行模拟,并验证模型模拟和计算的准确度和精确度。研究结果表明:(1)WRF-HYSPLIT和WRF-FLEXPART两种轨迹计算模式在虫源地、迁飞路径(迁飞方位角和走向)、迁飞高度、迁飞速率和迁飞距离计算上总体趋势一致,但存在一定的差异,后者的起伏变化大于前者。(2)尽管两种耦合模式在调用WRF模式输出的预报场物理变量方面有很多相同之处,但WRF-FLEXPART耦合模式在运行计算过程中比WRF-HYSPLIT耦合模式多考虑了对流参数、地表胁迫和各种地形参数,因而能更全面地反映中尺度天气过程(特别是对流性天气过程)对昆虫起飞、空中飞行和降落的动力作用,也能更真实地反映地表物理过程、大气湍流结构和地形起伏对褐飞虱种群迁飞的影响。(3)从褐飞虱种群对生境和取食条件选择上看,两种模式模拟的各高度迁入种群的虫源区、迁飞路径和降落地都是合理的、准确的。但从褐飞虱迁出、空中飞行和降落所处的三维流场来看,WRF-FLEXPART模式轨迹走向与盛行气流方向的吻合度要明显高于WRF-HYSPLIT模式。(4)两种模式均可作为业务工具在迁飞性害虫测报中推广应用。     

Live pine pollen in rainwater: reconstructing its long-range transport

Raindrops brim with pollen even when there is no ambient local pollen. How does this nonlocal pollen get inside rain? The likely answer is long-range transport beneath or inside clouds. To test this hypothesis, we captured rain-delivered pollen on Ocracoke Island, NC, USA over a 12-day interval before local pine pollen release then reconstructed its trajectory and its atmospheric exposure conditions. Findings were as follows: four rain episodes yielded a total of 632 pollen grains of which 6.7% germinated. To find pollen sources, we first identified pollen-releasing forested areas using a predictive heat sum equation for each rain episode. Next, we constructed the backward trajectory for air parcels carrying rain-delivered pollen from those forests using the MLDP atmospheric transport and dispersion model. Nonlocal sources were located at distances up to 300 km from Ocracoke Island and distances lessened with each successive episode. Below-cloud transport time was 8 and 17 h for Episodes A and B, respectively. Pollen grains were exposed to harsh atmospheric conditions during below-cloud transport, yet some grains still germinated. Atmospheric turbulence patterns changed for each episode, so distance from pollen source was poorly correlated with pollen transport time. Pollen germination was not closely correlated with either distances or transport time. In-cloud transport was more likely for pollen sampled during Episodes C and D. Pine pollen, although rarely allergenic, brings fresh insights into how atmospheric events can trigger human respiratory distress.

     

Population Structure of Mountain Pine Beetle Symbiont Leptographium longiclavatum and the Implication on the Multipartite Beetle-Fungi Relationships

Over 18 million ha of forests have been destroyed in the past decade in Canada by the mountain pine beetle (MPB) and its fungal symbionts. Understanding their population dynamics is critical to improving modeling of beetle epidemics and providing potential clues to predict population expansion. Leptographium longiclavatum and Grosmannia clavigera are fungal symbionts of MPB that aid the beetle to colonize and kill their pine hosts. We investigated the genetic structure and demographic expansion of L. longiclavatum in populations established within the historic distribution range and in the newly colonized regions. We identified three genetic clusters/populations that coincide with independent geographic locations. The genetic profiles of the recently established populations in northern British Columbia (BC) and Alberta suggest that they originated from central and southern BC. Approximate Bayesian Computation supports the scenario that this recent expansion represents an admixture of individuals originating from BC and the Rocky Mountains. Highly significant correlations were found among genetic distance matrices of L. longiclavatum, G. clavigera, and MPB. This highlights the concordance of demographic processes in these interacting organisms sharing a highly specialized niche and supports the hypothesis of long-term multipartite beetle-fungus co-evolutionary history and mutualistic relationships.     

Dispersal and influences on movement for Anoplophora glabripennis calculated from individual mark-recapture

We conducted an individual mark‐release‐recapture experiment on the beetle, Anoplophora glabripennis Motchulsky (Coleoptera: Cerambycidae). This invasive beetle has been introduced from Asia to Europe and North America and poses a serious threat to several important species of tree. Eradication efforts may benefit from knowledge of dispersal behaviour. Trees were cut and held to determine emergence rate of A. glabripennis. Unique marks were painted onto 912 beetles released into a group of 165 trees in Gansu, China. Data on subsequent sightings of beetles were used in a truncated diffusion model to calculate flight distances. Characteristics of the trees and climatic information were used in statistical tests for influence on movement. A total of 2245 sightings of beetles were observed and 29% of marked beetles were resighted. The scanning technique using binoculars was 90% effective in finding beetles and provided 81% accuracy for determining the sex of the beetles. Experimental manipulation of density quantified how A. glabripennis congregated on unoccupied trees and were repulsed from crowded hosts. The seasonal emergence rate of adults declined exponentially from July 20 to August 5. The results suggested A. glabripennis fly to nearby host trees at a rate of 34% per day. Median flight distance was estimated at 20 m per day. Statistical analysis with a generalized linear model tested the beetle's propensity to leave a tree and distance of flight. Generally, beetle movement showed a significant response to beetle density, weather conditions, beetle size, and tree size, in that order. The techniques developed here improve on previous recapture techniques to quantify dispersal and can be useful for analysing populations of other organisms.     

The role of atmospheric dispersion models and ecosystem sensitivity in the determination of characterisation factors for acidifying and eutrophying emissions in LCIA

Background, aim and scope  The methodological choices and framework to assess environmental impacts in life cycle assessment are still under discussion. Despite intensive developments worldwide, few attempts have been made hitherto to systematically present the role of different factors of characterisation models in life cycle impact assessment (LCIA). The aim of this study is to show how European average and country-dependent characterisation factors for acidifying and eutrophying emissions differ when using (a) acidifying and eutrophying potentials alone, (b) depositions from an atmospheric dispersion model or (c) critical loads in conjunction with those depositions. Furthermore, in the latter case, the contributions of emissions, an atmospheric transport model and critical loads to changes in characterisation factors of NO₂ are studied. In addition, the new characterisation factors based on the accumulated exceedance (AE) method are presented using updated emissions, a new atmospheric transport model and the latest critical loads. Materials and methods  In this study, characterisation factors for acidifying and eutrophying emissions are calculated by three different methods. In the ‘no fate’ (NF) methods, acidifying and eutrophying potentials alone are considered as characterisation factors. In the ‘only above terrestrial environment’ (OT) approach, characterisation factors are based on the deposition of the acidifying or eutrophying substances to terrestrial land surfaces. The third method is the so-called AE method in which critical loads are used in conjunction with depositions. The results of the methods are compared both at the European and the country level using weighted mean, weighted standard deviation, minimum and maximum values. To illustrate the sensitivity of the AE method, changes in European emissions, employed atmospheric dispersion model and the critical loads database are conducted step-by-step, and the differences between the results are analysed. Results and discussion  For European average characterisation factors, the three characterisation methods of acidification produce results in which the contributions of NH₃, NO₂ and SO₂ to the acidification indicator do not differ much within each method when 1 kg of each acidifying substance is emitted. However, the NF methods cannot describe any spatial aspects of environmental problems. Both OT and AE methods show that the spatial aspects play an important role in the characterisation factors. The AE method results in greater differentiations between country-dependent characterisation factors than does the OT method. In addition, the results of the AE and OT methods differ from each other for individual countries. A major shortcoming of the OT approach is that it does not consider the sensitivity of the ecosystems onto which the pollutants are deposited, whereas the AE approach does. In the case of the AE method, a new atmospheric dispersion model, new information on emissions and critical loads have a different influence on the characterisation factors, depending on the country. The results of statistics show that the change in the atmospheric dispersion model has a greatest influence on the results, since ecosystem-specific depositions are taken into account for the first time. Conclusions and recommendations  The simple NF methods can be used in a first approximation to assess the impacts of acidification and terrestrial eutrophication in cases where we do not know where the emissions occur. The OT approach is a more advanced method compared with the NF method, but its capability to describe spatial aspects is limited. The AE factors are truly impact-oriented characterisation factors and the information used here represents the current best knowledge about the assessment practice of acidification and terrestrial eutrophication in Europe. The key message of this study is that there is no shortcut to achieving advanced characterisation of acidification and terrestrial eutrophication: an advanced methodology cannot develop without atmospheric dispersion models and information on ecosystem sensitivity.     

Mountain pine beetle selectivity in old‐growth ponderosa pine forests,Montana, USA

A historically unprecedented mountain pine beetle (MPB) outbreak affected western Montana during the past decade. We examined radial growth rates (AD 1860–2007/8) of co‐occurring mature healthy and MPB‐infected ponderosa pine trees collected at two sites (Cabin Gulch and Kitchen Gulch) in western Montana and: (1) compared basal area increment (BAI) values within populations and between sites; (2) used carbon isotope analysis to calculate intrinsic water‐use efficiency (iWUE) at Cabin Gulch; and (3) compared climate‐growth responses using a suite of monthly climatic variables. BAI values within populations and between sites were similar until the last 20–30 years, at which point the visually healthy populations had consistently higher BAI values (22–34%) than the MPB‐infected trees. These results suggest that growth rates two–three decades prior to the current outbreak diverged between our selected populations, with the slower‐growing trees being more vulnerable to beetle infestation. Both samples from Cabin Gulch experienced upward trends in iWUE, with significant regime shifts toward higher iWUE beginning in 1955–59 for the visually healthy trees and 1960–64 for the MPB‐infected trees. Drought tolerance also varied between the two populations with the visually healthy trees having higher growth rates than MPB‐infected trees prior to infection during a multi‐decadal period of drying summertime conditions. Intrinsic water‐use efficiency significantly increased for both populations during the past 150 years, but there were no significant differences between the visually healthy and MPB‐infected chronologies.     

Diurnal flight pattern of Platypus koryoensis (Coleoptera: Platypodinae) in relation to abiotic factors in Korea

The diurnal flight pattern of Platypus koryoensis (Murayama) was examined using sticky traps attached to the trunks of oak trees in central Korea in 2010 and 2011. The flight activities of the beetle were estimated on the basis of 3-h intervals for trap catches from 11:00 to 14:00 on the next day, on June 25–26 and July 1–2, 2010, and on the basis of 2-h intervals for trap catches from 5:00 to 17:00, between June 9 and July 21, 2011 (the peak flight period of the beetle). Over 77% of the beetles were caught from 9:00 to 13:00, with the daily variations in the facing slope. The beetles began to be caught when the air temperature reached around 16 °C, and were the most active when the temperature ranged from 20 to 27 °C. No beetles were caught during rainfall, suggesting that rainfall is one of the factors that hinder beetle flight. The beetles were caught by traps in the east-facing plot earlier than those in the south- and west-facing plots suggesting that the flight behavior of the beetle can be affected by the light. Direction of the beetle flight during the peak of daily flight (from 09:00 to 13:00) was downward along the slope.     

Atmospheric trajectories and spring bird migration across the Gulf of Mexico

We examined the relationship between the longitude of peak arrival of trans-Gulf migrants on the northern coast of the Gulf of Mexico in spring and wind trajectories over the Gulf at three different altitudes (500, 1,500, and 2,500 m above ground level). We used data from 10 WSR-88D radars (weather surveillance radar-1988-Doppler) from Brownsville, Texas, to Key West, Florida, to record the time and longitude of peak arrival on the northern Gulf coast for four spring migrations (2001–2004). We used the National Oceanic Atmospheric Administration Air Resources Laboratory HYSPLIT transport and dispersion model at the READY Web site to generate backward, 24-h atmospheric trajectories based on archived atmospheric data for each trans-Gulf flight. The trajectories began at the geographic location where radar indicated the greatest concentrations of arriving migrants. Although the longitude of peak arrival varied, peak densities of most trans-Gulf migrants arrived on the northern coast near longitude 95°W. Regression analyses showed that the relationship between the longitude of peak trans-Gulf arrival and the direction of atmospheric trajectory was significant but weak at the 500-m level, where few migrants occurred, and was insignificant for the 1,500- and 2,500-m altitudes, where migrant densities were greater. We conclude that winds aloft over the Gulf have little influence on the longitude of peak trans-Gulf arrival on the northern coast of the Gulf of Mexico, and we speculate that the arrival pattern may reflect the trans-Gulf migration pathways that evolved during the Last Glacial Maximum.     

Mountain pine beetle develops an unprecedented summer generation in response to climate warming

The mountain pine beetle (MPB; Dendroctonus ponderosae) is native to western North America, attacks most trees of the genus Pinus, and periodically erupts in epidemics. The current epidemic of the MPB is an order of magnitude larger than any previously recorded, reaching trees at higher elevation and latitude than ever before. Here we show that after 2 decades of air-temperature increases in the Colorado Front Range, the MPB flight season begins more than 1 month earlier than and is approximately twice as long as the historically reported season. We also report, for the first time, that the life cycle in some broods has increased from one to two generations per year. Because MPBs do not diapause and their development is controlled by temperature, they are responding to climate change through faster development. The expansion of the MPB into previously inhospitable environments, combined with the measured ability to increase reproductive output in such locations, indicates that the MPB is tracking climate change, exacerbating the current epidemic.     

Modelling analysis of source regions of long-range transported birch pollen that influences allergenic seasons in Lithuania

This study analyses the spatial and temporal distribution of regional and long-range transported birch (Betula L.) pollen in Lithuania and the neighbouring countries. The potential long-range transport cases of birch pollen in Lithuania were analysed for the whole period of available observations, 2004–2007. The birch pollen was recorded at three measurement stations in Lithuania by using Hirst-type volumetric spore traps. The phenological observations in Lithuania were also used for the detection of potential long-range transport-induced episodes. Two variants of the regional and continental scale atmospheric dispersion model SILAM (Lagrangian and Eulerian) in an adjoint mode (used for inverse dispersion modelling and data assimilation), and the trajectory model HYSPLIT were employed to evaluate the source origins of the observed pollen. During four seasons in 2004–2007, we found in total 24 cases, during which remarkable pollen concentrations were recorded before the local flowering season. According to modelling, most of these were originated from the sources outside Lithuania: Latvia, southern Sweden, Denmark, Belarus, Ukraine and Moldova, possibly, also coastal regions of Germany and Poland. Two episodes were attributed to local early-flowering birch trees. The spatial and temporal patterns of the long-range transport of early pollen to Lithuania were found out to be highly variable; the predicted source regions for the cases considered were similar only for some dates in 2004 and 2006. During the analysed period, we found both cases, in which the predictions of the SILAM model variants and those of the HYSPLIT model were similar, and cases, in which there were substantial differences. In general, for complicated atmospheric circulation patterns the model predictions can be drastically different, with a tendency of trajectory model to fail reproducing the key episode features.     

Long-term black carbon dynamics in cultivated soil

Black carbon (BC) is a quantitatively important C pool in the global C cycle due to its relative recalcitrance compared with other C pools. However, mechanisms of BC oxidation and accompanying molecular changes are largely unknown. In this study, the long-term dynamics in quality and quantity of BC were investigated in cultivated soil using X-ray photoelectron spectroscopy (XPS), Fourier-transform infrared (FTIR) and nuclear magnetic resonance (NMR) techniques. BC particles and changes in BC stocks were obtained from soil collected in fields that were cleared from forest by fire at 8 different times in the past (2, 3, 5, 20, 30, 50, 80 and 100 years before sampling) in western Kenya. BC contents rapidly decreased from 12.7 to 3.8 mg C g⁻¹ soil during the first 30 years following deposition, after which they slowly decreased to a steady state at 3.5 mg C g⁻¹ soil. BC-derived C losses from the top 0.1 m over 100 years were estimated at 6,000 kg C ha⁻¹. The initial rapid changes in BC stocks resulted in a mean residence time of only around 8.3 years, which was likely a function of both decomposition as well as transport processes. The molecular properties of BC changed more rapidly on surfaces than in the interior of BC particles and more rapidly during the first 30 years than during the following 70 years. The Oc/C ratios (Oc is O bound to C) and carbonyl groups (C=O) increased over the first 10 and 30 years by 133 and 192%, respectively, indicating oxidation was an important process controlling BC quality. Al, Si, polysaccharides, and to a lesser extent Fe were found on BC particle surfaces within the first few years after BC deposition to soil. The protection by physical and chemical stabilization was apparently sufficient to not only minimize decomposition below detection between 30 and 100 years after deposition, but also physical export by erosion and vertical transport below 0.1 m.     

A field experiment on dispersal of newly emerged adults ofMonochamus alternatus (Coleoptera: Cerambycidae)

Newly-emerged adults of Monochamus alternatus aged 1 to 5 days were code-numbered with lacquer paint and released by placing them on the trunks of one or two trees in a Pinus thunbergii stand at weekly intervals during the beetle emergence period from 1980 to 1983. Beetles were captured at weekly intervals from one week after the first day of release. Determinations were made on the distance and direction of beetle dispersal during a week after release and analysed by a method of Inoue (1978). When the stand canopy was closed, the rate of beetle's stay on trees was 0.56 per week. The beetles dispersed at random by walk and flight. When the pine stand was sparse, the rate of beetle's stay on trees was 0.02–0.30 per week. They dispersed at random by flight. The average distances traversed were estimated to be 7.1–37.8 m for the first week after emergence. Using other method, the average distance traversed was estimated to be 10–20 m for each week through the first 3 weeks after release. The results of stepwise multiple regression analysis and a simple field experiment suggested that the dispersal of newly-emerged beetles was affected by stand density, number of beetles emerging from individual dead trees, maximum air temperature, and precipitation.     

Phase transition from environmental to dynamic determinism in mountain pine beetle attack

Mountain pine beetle outbreaks are responsible for widespread tree mortality in pine forests throughout western North America. Intensive outbreaks result in significant economic loss to the timber industry and massive changes to the forest habitat. Because of the time and space scales involved in a beetle outbreak, mathematical models are needed to study the evolution of an outbreak. In this paper we present a partial differential equation model of the flight phase of the mountain pine beetle which includes chemotactic responses and tree defense. We present a numerical method for integrating this model and use this method to investigate the relationship between emergence rate, forest demographic and patterns of beetle attack. In particular we look at how emergence rate affects the beetles' ability to successfully attack strong trees, which may be an indicator of an epidemic outbreak.     

A dispersion modelling approach to determine the odour impact of intensive poultry production units in Ireland

The use of atmospheric dispersion modelling has become more common for the determination of odour impacts from existing poultry production facilities and the assessment of setback distances for new facilities. Setback distances for broiler, layer and turkey units were determined using the atmospheric dispersion model ISCST3 and the Environmental Protection Agency (EPA, Ireland) recommended criterion (C(98,1-h)6.0 ou(E) m(-3)) and a new odour annoyance criterion (C(98,1-h) 9.7 ou(E) m(-3)) developed in this study. For a typical size unit in Ireland, maximum setback distances of 660, 665 and 1035 m were calculated for 40,000 broilers, 40,000 layers and 10,000 turkeys respectively at the current limit (C(98,1-h) 6.0 ou(E) m(-3)). However, if the suggested odour impact criterion (C(98,1-h) 9.7 ou(E) m(-3)) is implemented, the maximum setback distances decrease to 460, 500 and 785 m for broilers, layers and turkeys, respectively.     

Analysis of atmospheric dispersion of olive pollen in southern Spain using SILAM and HYSPLIT models

SILAM atmospheric dispersion model and the HYSPLIT trajectory model were used to detect the source areas and calculate transport dynamics for airborne olive pollen observed in the city of Córdoba, southwest of Iberian Peninsula. The ECMWF weather data with 3-h time interval and spatial resolution of 25 × 25 km² and 75 hybrid vertical levels were used as meteorological inputs in both models to produce a coherent set of results in order to compare these two different approaches. Seven episodes recorded before and after the local flowering season in 2006 were analyzed using both models. The results provided an indication of the origins of olive pollen recorded in the city of Córdoba, revealing the influence of three main source areas at specific periods. One area was located nearby, to the southwest of the city (early May), another in the south of the province (mid-May) and the third to the east (late May/early June). The SILAM model yielded more detailed and quantitative results when identifying olive pollen sources and charting transport dynamics. The results from the HYSPLIT trajectory approach and SILAM footprints were qualitatively similar. However, a weak point of back trajectories was their lower sensitivity to details of the transport, as well as the necessity of subjective analysis of the trajectory plots, which were subject for possible misinterpretations. Information on both pollen source locations and local tree flowering phenology was required in order to ensure consistent analysis of the influence of olive sources for both models. Further than this, due to the fact that both models are widely used in other research areas, the results of this work could have a widespread range of application, such as to simulate the transport of radionuclides, e.g., in emergency preparedness exercises.     

Damaged forests provide an opportunity to mitigate climate change

British Columbia (BC) forests are estimated to have become a net carbon source in recent years due to tree death and decay caused primarily by mountain pine beetle (MPB) and related post‐harvest slash burning practices. BC forest biomass has also become a major source of wood pellets, exported primarily for bioenergy to Europe, although the sustainability and net carbon emissions of forest bioenergy in general are the subject of current debate. We simulated the temporal carbon balance of BC wood pellets against different reference scenarios for forests affected by MPB in the interior BC timber harvesting area using the Carbon Budget Model of the Canadian Forest Sector (CBM‐CFS3). We evaluated the carbon dynamics for different insect‐mortality levels, at the stand‐ and landscape level, taking into account carbon storage in the ecosystem, wood products and fossil fuel displacement. Our results indicate that current harvesting practices, in which slash is burnt and only sawdust used for pellet production, require between 20–25 years for beetle‐impacted pine and 37–39 years for spruce‐dominated systems to reach pre‐harvest carbon levels (i.e. break‐even) at the stand‐level. Using pellets made from logging slash to replace coal creates immediate net carbon benefits to the atmosphere of 17–21 tonnes C ha^?1, shortening these break‐even times by 9–20 years and resulting in an instant carbon break‐even level on stands most severely impacted by the beetle. Harvesting pine dominated sites for timber while using slash for bioenergy was also found to be more carbon beneficial than a protection reference scenario on both stand‐ and landscape level. However, harvesting stands exclusively for bioenergy resulted in a net carbon source unless the system contained a high proportion of dead trees (>85%). Systems with higher proportions of living trees provide a greater climate change mitigation if used for long lived wood products.     

环境因素对黄曲条跳甲种群扩散的影响

应用标记回收技术,研究了寄主及气象因素对黄曲条跳甲种群扩散的影响.结果表明,黄曲条跳甲的自然扩散是由环境引起(τ≈1),温度、风速与之有一定的相关性(P<0.05),高温不利于黄曲条跳甲的扩散,其扩散主要沿逆风或与风向垂直方向进行,而湿度与黄曲条跳甲成虫扩散无显著相关(P>0.05).寄主是影响黄曲条跳甲扩散的重要因子,在几种寄主植物中,菜心上收集的标记成虫数最多,其次为芥菜,再次为白菜,芥蓝上基本没有.     

Stage‐specific probabilistic phenology model of Cydia pomonella (Lepidoptera: Tortricidae) using laboratory maximum likelihood parameter estimates

In this study, a probabilistic degree‐day phenology model has been developed for the codling moth, Cydia pomonella, and calibrated using data from laboratory growth studies. The model is further used to predict the succession and overlapping of certain biological events of C. pomonella in probabilistic‐physiological time scale in northern Greece fruit orchards. The model satisfactorily predicts the stage‐specific pest population dynamics, including egg laying and hatching, the occurrence of larvae and pupae stages and the emergence of adults. According to the model projections for the adult flights, there is a very high probability, p = 0.999, of observing adults of the first flight generation until 333 degree‐days (DD), but a very low probability of finding adults of the second flight generation. Moreover, at 575 DD, the probability of finding an individual to lay eggs is p = 0.15. However, there is nearly the same probability of egg hatch, p = 0.36, and larval completion p = 0.313, while at the same time, the probability of pupal completion is very low, p = 0.001. The above model predictions were validated using field data for the adult stage emergence as well as for the percentage of larval damage providing satisfactory results considering that larval emergence prediction was close to actual fruit damage observed in field. This information is very important considering that IPM programs rely on the use of biorational compounds, such as IGRs and bio‐toxins which are stage selective and often have a shorter residual activity than the preceding broad‐spectrum insecticides.     

Broadscale Specificity in a Bark Beetle–Fungal Symbiosis: a Spatio-temporal Analysis of the Mycangial Fungi of the Western Pine Beetle

Whether and how mutualisms are maintained through ecological and evolutionary time is a seldom studied aspect of bark beetle–fungal symbioses. All bark beetles are associated with fungi and some species have evolved structures for transporting their symbiotic partners. However, the fungal assemblages and specificity in these symbioses are not well known. To determine the distribution of fungi associated with the mycangia of the western pine beetle (Dendroctonus brevicomis), we collected beetles from across the insect’s geographic range including multiple genetically distinct populations. Two fungi, Entomocorticium sp. B and Ceratocystiopsis brevicomi, were isolated from the mycangia of beetles from all locations. Repeated sampling at two sites in Montana found that Entomocorticium sp. B was the most prevalent fungus throughout the beetle’s flight season, and that females carrying that fungus were on average larger than females carrying C. brevicomi. We present evidence that throughout the flight season, over broad geographic distances, and among genetically distinct populations of beetle, the western pine beetle is associated with the same two species of fungi. In addition, we provide evidence that one fungal species is associated with larger adult beetles and therefore might provide greater benefit during beetle development. The importance and maintenance of this bark beetle–fungus interaction is discussed.