Strong indirect interactions of Tarsonemus mites (Acarina: Tarsonemidae) and Dendroctonus frontalis (Coleoptera: Scolytidae)

Phoretic mites of bark beetles are classic examples of commensal ectosymbionts. However, many such mites appear to have mutualisms with fungi that could themselves interact with beetles. We tested for indirect effects of phoretic mites on Dendroctonus frontalis, which attacks and kills pine trees in North America. Tarsonemus mites are known to carry ascospores of Ophiostoma minus, which tends to outcompete the mutualistic fungi carried by D. frontalis. Experimental additions and removals of mites from beetles demonstrated that Tarsonemus propagate O. minus in beetle oviposition galleries. Furthermore, the abundance of Tarsonemus and O. minus tended to covary in nature. These results verified a strong mutualism between Tarsonemus and O. minus. Results also indicated that O. minus is an antagonist of D. frontalis: beetle larvae seldom survived in the presence of O. minus (compared to 83% survival elsewhere). Apparently, this is an indirect result of O. minus outcompeting the two species of mycangial fungi that are critical to beetle nutrition. Thus, Tarsonemus mites close a loop of species interactions that includes a commensalism (mites and beetles), a mutualism (mites and O. minus), asymmetric competition (O. minus and mycangial fungi), and another mutualism (mycangial fungi and beetles). This interaction system produces negative feedback that could contribute to the endogenous population dynamics of D. frontalis. Reproductive rate of Tarsonemus was more temperature‐sensitive than beetle generation time (which constrains the time for mite reproduction within a tree). This differential temperature sensitivity produces a narrow range of temperatures (centred at 27°C) in which mite reproduction per D. frontalis generation can attain its maximum of 100 mites/beetle. Consequently, seasonal oscillations in temperature are predicted to produce oscillations in the D. frontalis community, and climatic differences between regions could influence the community to dampen or exacerbate the cyclical outbreak dynamics of D. frontalis.     

Cold tolerance and invasive potential of the redbay ambrosia beetle (<Emphasis Type="Italic">Xyleborus glabratus</Emphasis>) in the eastern United States

Native Lauraceae (e.g. sassafras, redbay) in the southeastern USA are being severely impacted by laurel wilt disease, which is caused by the pathogen Raffaelea lauricola T. C. Harr., Fraedrich and Aghayeva, and its symbiotic vector, the redbay ambrosia beetle (Xyleborus glabratus Eichhoff). Cold temperatures are currently the only viable limitation to the establishment of X. glabratus in northern populations of sassafras. The observed lower lethal temperature of X. glabratus (? 10.0 °C) is warmer than its supercooling point (? 22.0 °C), indicating the beetle is a freeze intolerant and chill susceptible species. Empirically derived X. glabratus lower lethal temperature thresholds were combined with host distribution and microhabitat-corrected climate data to produce species distribution models for X. glabratus in the eastern USA. Macroclimate data (30-year mean annual minimum temperature) were corrected (? 1.2 °C) to account for thermal buffering afforded to X. glabratus while living inside sassafras trees. Only 0.1% of the current US sassafras spatial extent experiences sufficiently harsh winters (locales where mean annual minimum winter temperatures ≤ ? 6.2 °C for ≥ 12 h) to exclude X. glabratus establishment in our species distribution model. Minimum winter temperatures will likely cause some X. glabratus mortality in ~ 52% of the current spatial extent of sassafras, although current data do not allow a quantification of X. glabratus mortality in this zone. Conversely, ~ 48% of the current spatial extent of sassafras is unlikely to experience sufficiently cold winter temperatures to cause any significant impediment to X. glabratus spread or establishment. A modest climate change scenario (RCP4.5) of + 1.4 °C would result in 91% of the current spatial extent of sassafras in the eastern USA occurring where winter minimum temperatures are unlikely to cause any mortality to X. glabratus.     

Predation and bark beetle dynamics

Bark beetle populations may undergo dramatic fluctuations and are often important pests in coniferous forests. Their dynamics are thought to be primarily driven by factors affecting the resistance of the host tree to attack, i.e., bottom-up forces, while natural enemies are usually assigned a minor role in these systems. I present behavioral experiments that suggest that the clerid beetle Thanasimus dubius may be an important source of mortality for the bark beetle Dendroctonus frontalis during attack of the host tree, and determine the nature of the functional response of T. dubius under conditions close to natural. I also examine the numerical response of T. dubius to large-scale fluctuations in D. frontalis density, and the relationship between bark beetle population trends and predator density, and find that beetle populations tend to decline when predator densities are high. Combined with the effects of clerid larvae on bark beetle broods, these results suggest that top-down forces generated by natural enemies could also be an important component of bark beetle dynamics. The implications of these results for bark beetle dynamics are discussed in relation to the prolonged life-cycle of clerid beetles. Received: 23 January 1997 / Accepted: 5 April 1997     

Climate change and the outbreak ranges of two North American bark beetles

Abstract

• 1 One expected effect of global climate change on insect populations is a shift in geographical distributions toward higher latitudes and higher elevations. Southern pine beetle Dendroctonus frontalis and mountain pine beetle Dendroctonus ponderosae undergo regional outbreaks that result in large‐scale disturbances to pine forests in the south‐eastern and western United States, respectively.
• 2 Our objective was to investigate potential range shifts under climate change of outbreak areas for both bark beetle species and the areas of occurrence of the forest types susceptible to them.
• 3 To project range changes, we used discriminant function models that incorporated climatic variables. Models to project bark beetle ranges employed changed forest distributions as well as changes in climatic variables.
• 4 Projected outbreak areas for southern pine beetle increased with higher temperatures and generally shifted northward, as did the distributions of the southern pine forests.
• 5 Projected outbreak areas for mountain pine beetle decreased with increasing temperature and shifted toward higher elevation. That trend was mirrored in the projected distributions of pine forests in the region of the western U.S. encompassed by the study.
• 6 Projected outbreak areas for the two bark beetle species and the area of occurrence of western pine forests increased with more precipitation and decreased with less precipitation, whereas the area of occurrence of southern pine forests decreased slightly with increasing precipitation.
• 7 Predicted shifts of outbreak ranges for both bark beetle species followed general expectations for the effects of global climate change and reflected the underlying long‐term distributional shifts of their host forests.

     

Influence of temperature on the northern distribution limits of Scirpophaga incertulas Walker (Lepidoptera: Pyralidae) in China

We explored the influence of temperature on the northern distribution limits of Scirpophaga incertulas Walker, an important agricultural pest of rice in Asia. We analyzed ≥48 years of records from 186 climate stations of Mainland China to estimate the annual probabilities of reaching the lower lethal temperature for S. incertulas. The relevant climatic metric, minimum annual temperature, approximated a normal distribution. Consequently, the probability density function for any site could be characterized with the mean and standard deviation of minimum annual temperatures. We used the local regression method to map the mean and standard deviation of minimum annual temperatures throughout Mainland China and then calculated isolines representing annual probabilities for reaching or exceeding the lower lethal temperature of S. incertulas. In addition, we calculated and mapped the number of generations per year based on the annual accumulative degree days and the sum of effective temperatures required to complete one generation. The empirical northern distribution limits of S. incertulas were generally congruent with the theoretical limits based on winter survival, with exceptions within the Shandong and Sichuan provinces, which are apparently thermally suitable but where the host plant is not cultivated. The expected number of generations per year was 3–5 within most of the range of S. incertulas in China. In central China, the expected number of generations per year was about 3. A climate warming scenario of 4 °C in minimum and maximum daily temperatures predicted an increase in the expected number of generations per year in central China from about 3 to 4.     

Influence of temperature on the reproductive success,brood development and brood fitness of the eastern larch beetle Dendroctonus simplex LeConte

相似文献   

Climatic controls on distribution of Fagus crenata forests in Japan

Abstract. We used classification tree analysis to develop a climate‐based distribution model for Fagus crenata forests in Japan. Four climatic variables judged likely to affect the distribution of the species (summer and winter precipitation, minimum temperature of the coldest month and Kira's warmth index) were chosen as independent variables for the model. Latitudinal and longitudinal information was also used to examine effects of spatial autocorrelation on the model. The climatic factors associated with the distribution of the forests were analysed using a classification tree to devise prediction rules. Predicted areas of high probability for forest occurrence lay mainly on the Sea of Japan side of northern Honshu and southern Hokkaido. This is consistent with actual forest distribution. Some areas with high predicted probabilities of F. crenata forest occurrence were beyond the current natural northern range limits of these forests. Since these areas were widely scattered, it was assumed that the species has been hindered from colonizing them due to dispersal limitations. Deviance‐weighted scores, used to compare magnitudes of the contributions of predictor variables, revealed winter precipitation as the most influential factor, followed by the warmth index, the minimum temperature of the coldest month and summer precipitation. Attempts were made to generate ecological explanations for the effects of the four climatic factors on the distribution of F. crenata forests.     

Impact of climatic variation on populations of pine processionary moth Thaumetopoea pityocampa in a core area of its distribution

• 1 There is growing appreciation of climatic effects on insect population dynamics at the margins of distribution limits. Climatic effects might be less important and/or involve different drivers and processes near the centre of distributions.
• 2 We evaluated the effects of interannual variation in temperatures, radiation and precipitation on populations of pine processionary moth Thaumetopoea pityocampa in Central–South Portugal, a low altitude area with a relatively mild Mediterranean climate near the centre of the north–south range of the species.
• 3 We tested for effects of climate on mortality of young larvae, growth rates, final larval mass and fecundity.
• 4 Results indicated high mortality of early instars associated with low minimum and maximum daily temperatures and low precipitation. Low minimum temperatures were further associated with high parasitism by the larval parasitoid Phryxe caudata (Rondani) (Diptera, Tachinidae). Furthermore, larval growth rates were higher with high solar radiation during December and January, which was itself negatively related to precipitation and air temperature. Slow larval growth rates led to lower final mass at pupation in the spring, and smaller egg masses and smaller initial colony sizes during the next autumn.
• 5 Thus climatic factors, and temperature in particular, apparently contributed to population dynamics of T. pityocampa in the core of its distribution, as well as at its northern limits. The most specific climatic parameters of importance, however, and the connections between climate, physiology and insect demographics in the core area were clearly different from northern areas.

     

Geographical distributions,relative abundance and coexistence of Drosophila aldrichi and Drosophila buzzatii in Australia

Abstract Climatic data and collection records for the cactophilic Drosophila aldrichi and Drosophila buzzatii for 97 localities were used to examine the effects of geographical location, season, host plant species and climatic factors on their range and relative abundance. Temporal variation in relative abundance was assessed from monthly collections over 4 years at one locality. Effects of weather variables over the 28 days before each collection were examined. A generalized linear model of the spatial data showed significant geographical variation in relative abundance, and significant climatic effects, with the proportion of D. aldrichi higher in the warm season, and increasing as temperature variation decreased and moisture indices increased. The temporal data gave generally concordant results, as D. aldrichi proportion was higher in summer and autumn, and increased as maximum and minimum temperatures increased, and as variation in maximum temperature decreased. In a laboratory competition experiment, D. aldrichi eliminated D. buzzatii at 31°C, but was itself eliminated at 18°C and 25°C. The range of D. buzzatii is constrained only by availability of its host plant, Opuntia species, although its relative abundance is reduced in the northern part of its distribution. The range of D. aldrichi, from central Queensland to northern NSW, Australia, is entirely within that of D. buzzatii, and its relative abundance decreases from north to south. Both climate and weather, particularly temperature variability, have direct effects on the relative abundances of the two species, and both likely act indirectly by influencing the outcome of interspecific competition.     

Complex emergence patterns in a bark beetle predator

1 The emergence pattern of Thanasimus dubius (F.) (Coleoptera: Cleridae), a common predator of the southern pine beetle, Dendroctonus frontalis Zimmermann (Coleoptera: Scolytidae), was studied under field conditions across different seasons. A simple statistical model was then developed to characterize the emergence data, using the truncated geometric distribution. Data are also presented on the mortality of T. dubius eggs at various temperatures and humidities in an effort to explain certain aspects of emergence behaviour. 2 Emergence of T. dubius from a given tree usually occurred in several discrete episodes across a two‐year period, with most individuals emerging in spring or autumn. Almost no emergence occurred in July and August, which may be an adaptation to avoid high temperature mortality. Emergence patterns appeared similar across seasons, with the time of year serving mainly to shift the pattern through time. 3 Cycles in D. frontalis abundance may be the result of delayed density dependence generated by its natural enemy complex. The predator T. dubius is likely to be an important component of this delayed density dependence, because of its lengthy development time and apparent impact on D. frontalis.     

Testing temperature effects on woodboring beetles associated with oak dieback

A warmer climate may potentially have a strong effect on the health status of European oak forests by weakening oak trees and facilitating mass reproduction of wood boring insects. We did a laboratory experiment in Slovakia to study the response of major pest beetles of oak and their parasitoids to different temperature regimes as background for predicting climatic effects and improving management tools of European oak forests. With higher temperatures the most important oak pest Scolytus intricatus emerged much earlier, which indicate that completion of a second generation and increased damage further north in European oak forests may be possible. Lower temperatures gave longer larval galleries and more offspring per parents but still lower beetle production due to semivoltine life cycle. For buprestids and longhorn beetles warmer temperatures resulted in more emerging offspring and a shift towards earlier emergence in the same season, but no emergence in the first season indicated that a change to univoltine populations is not likely. Reduced development success of parasitoids at the highest temperatures (25/30 °C) indicates a loss of population regulation for pest beetle populations. A warmer climate may lead to invasion of other population-regulating parasitoids, but also new serious pest may invade. With expected temperature increases it is recommended to use trap trees both in April and in June, and trap trees should be removed within 2 months instead 1 year as described in the current standard.     

Northward expansion and climatic factors affecting the distribution limits of Cheilomenes sexmaculata (Coleoptera: Coccinellidae) in Japan

The distributional expansion of the ladybird beetle Cheilomenes sexmaculata (Fabricius) in Japan was analyzed based on literature surveys and observations of 1,312 individual specimens collected from 1918 to 2005. An expansion toward higher latitudes, from Kyushu (33°N) to central Honshu (36°N), was observed between 1910 and the early 1990s. Multiple regression analysis showed that the presence of this beetle was associated with an annual mean temperature but not with the mean minimum temperature in the coldest month, and 95 % occurrence probability of this beetle was obtained in a locality at 15.1 °C annual mean temperature. The northward expansion of this species synchronizes with increasing areas where an annual mean temperature is more than 15 °C. Therefore, climatic warming over the past 100 years is considered as a major factor influencing the northward expansion of this species.     

Synergistic effects of α-pinene and exo-brevicomin on pine bark beetles and associated insects in Arizona

The southern pine beetle (Dendroctonus frontalis) and western pine beetle (Dendroctonus brevicomis) cause significant mortality to pines in the southern and western United States. The effectiveness of commercial lures at capturing these bark beetles in Arizona has not been tested and may vary from other regions of their distribution. We conducted experiments using baited Lindgren funnel traps to investigate (i) if D. frontalis is more attracted to the standard commercial lure for D. brevicomis (frontalin + exo‐brevicomin + myrcene) than the D. frontalis lure (frontalin + terpene blend), (ii) whether replacement of myrcene with α‐pinene changes trap catches of Dendroctonus and associated insects, and (iii) whether the attraction to these lures varies across the geographical range of ponderosa pine forests throughout Arizona. In 2005, we tested various combinations of frontalin, exo‐brevicomin, myrcene and α‐pinene to D. frontalis, D. brevicomis and associated species. Dendroctonus frontalis, D. brevicomis and the predator Temnochila chlorodia were most attracted to lures with exo‐brevicomin. The replacement of the myrcene component with α‐pinene in the D. brevicomis lure resulted in the capture of twice as many bark beetles and Elacatis beetles. However, T. chlorodia did not differentiate between monoterpenes. In 2006, traps were set up in 11 locations around Arizona to test the relative attraction of lure combinations. In 9 out 11 locations, the D. brevicomis lure with α‐pinene was more attractive than the lure with myrcene or a terpene blend. These results suggest that the D. brevicomis lure with α‐pinene rather than myrcene is more effective lure to capture D. brevicomis and D. frontalis in Arizona. However, geographical variation in attractiveness to lures is evident even within this region of the beetles’ distributions. Differential attraction of Dendroctonus and their predators to these lures suggests potential use in field trapping and control programmes.     

Continuous population estimates forDendroctonus frontalis Zimm. (Coleoptera: Scolytidae) occurring in infestations

Infestations of Dendroctonus frontalisZimm. are often observed to enlarge continuously by the colonization of new hosts in a pattern similar to a forest fire. This pattern of infestation growth presents unique problems in quantitatively estimating populations of D. frontalis. Beetle populations on each infested tree in an infestation go through five processes: attack, oviposition, reemergence, survivorship, and emergence. These processes, which have been described mathematically in the literature, each take several days for completion. In order to follow the distribution and abundance of D. frontalis throughout the course of development of a spot, we need a daily estimate of the number of beetles involved in each process on every tree. Since it is not practical to sample each tree daily, we developed a procedure whereby quantitative estimation procedures for within-tree populations were used in combination with the mathematical models for the life processes to produce a daily record of the number of adults successfully attacking trees, the number of eggs oviposited, the number of beetles reemerging, number of beetles surviving within the trees, and the number of beetles emerging. These daily estimates were then summarized for all trees in the spot for the duration of the infestation. The daily record of populations of D. frontalis, used with information on infestation geometry, were suggested to be of value in describing and elucidating several important facets of population dynamics including dispersal patterns within infestations, between tree beetle loss (mortality), and time lags among the various population processes. The information reported can be used to develop simulation models of population dynamics or to validate existing models.     

The suppression of Dendroctonus frontalis and subsequent wildfire have an impact on forest stand dynamics

Question: Interacting disturbance effects from Dendroctonus frontalis outbreaks and wildfire are thought to maintain Pinus spp. composition in the southeastern U. S. Our objective was to assess forest composition, structure, and succession following the interaction of two frequently occurring disturbance events in southern Pinus spp. forests: cut‐and‐leave suppression, a commonly used means for managing D. frontalis outbreaks, and wildfire. Location: Western Gulf Coastal Plain, Louisiana, USA. Method: Pinus taeda stands with cut‐and‐leave suppression and subsequent wildfire were compared to stands undisturbed by D. frontalis but with the same wildfire events twenty years after Pinus spp. mortality. The woody plant community was assessed in three different size classes and used to predict future forest types with the Forest Vegetation Simulator (50 years). Results: P. taeda is the most abundant (> 50%) species of saw‐ and poletimber‐sizes following cut‐and‐leave suppression with wildfire and in stands only with fire. Using canonical correspondence analysis, vegetation assemblages were primarily explained by slope position and elevation (7.6% variation explained). Fire intensity and stand age also accounted for variance in the ordination (4.4% and 3.1%, respectively). Dominant and co‐dominant P. taeda forest types were predicted by the model to be the most abundant forest types in each disturbance regime. In addition, new regeneration represents high hazard for future mortality from D. frontalis. Conclusion: Our study demonstrates that cut‐and‐leave suppression with additional wildfire disturbance maintains P. taeda composition, and does not alter forest composition differently from stands receiving only wildfire. As a result, predicted Pinus spp. basal area under both disturbances is great enough to facilitate future bark beetle disturbance.     

Partial age-specific life tables for Monochamus titillator in Dendroctonus frontalis infested loblolly pines

Partial age-specific life tables were constructed for Monochamus titillator(Fabricius) (Coleoptera: Cerambycidae) in Dendroctonus frontalisZimmermann (Coleoptera: Scolytidae) infested loblolly pine (Pinus taedaL.) trees. Successive life stages including egg, early larvae, mid-stage larvae, late larvae, and adult emergence were sampled within six sample trees. Generation mortality ranged from 60.94% to 98.61% in sample trees. Highest mortality typically occurred to eggs and mid-stage larvae. Possible mortality factors included resinosis, predaceous beetles, parasitoids, and woodpeckers. Dendroctonus frontalisbrood stages were determined for consecutive M. titillatorsampling. Monochamus titillatorand D. frontaliscoexisted and likely interacted in the phloem of host trees for at least 20 days.     

Importance of resin ducts in reducing ponderosa pine mortality from bark beetle attack

The relative importance of growth and defense to tree mortality during drought and bark beetle attacks is poorly understood. We addressed this issue by comparing growth and defense characteristics between 25 pairs of ponderosa pine (Pinus ponderosa) trees that survived and trees that died from drought-associated bark beetle attacks in forests of northern Arizona, USA. The three major findings of our research were: (1) xylem resin ducts in live trees were >10% larger (diameter), >25% denser (no. of resin ducts mm⁻²), and composed >50% more area per unit ring growth than dead trees; (2) measures of defense, such as resin duct production (no. of resin ducts year⁻¹) and the proportion of xylem ring area to resin ducts, not growth, were the best model parameters of ponderosa pine mortality; and (3) most correlations between annual variation in growth and resin duct characteristics were positive suggesting that conditions conducive to growth also increase resin duct production. Our results suggest that trees that survive drought and subsequent bark beetle attacks invest more carbon in resin defense than trees that die, and that carbon allocation to resin ducts is a more important determinant of tree mortality than allocation to radial growth.     

Cold hardiness of Lymantria monacha and L. dispar (Lepidoptera: Erebidae) eggs to extreme winter temperatures: implications for predicting climate change impacts

1. It has been predicted that temperature increases of 3.6–5.8 °C would shift the northern distribution limit of Lymantria monacha (Linnaeus) and Lymantria dispar (Linnaeus) by 500–700 km, but these predictions ignore the effects of minimum winter temperatures. It was hypothesised that winter cold can limit range expansion due to high egg mortality in cold temperatures. 2. The present study determined the supercooling points of overwintering eggs of these forest pests, and compared these with recent minimum winter temperatures in the areas of origin of three populations. Eggs from one L. monacha and one L. dispar population from the species' core distribution area in Germany were included, as well as L. monacha eggs from Finland, near the northern border of the species' distribution. 3. The median supercooling points of both species were more than 10 °C lower than the median minimum winter temperatures of their areas of origin, and the median supercooling points of Finnish and German L. monacha eggs did not differ significantly. The median supercooling point of German L. monacha eggs differed from that of German L. dispar eggs. 4. Previous literature on the topic is referenced, and translations of the old German and Russian sources are given. Based on these results, it is argued that the frequent claim that L. monacha eggs can survive cold down to ?40 °C is unsupported, with a value near ?30 °C being a more likely limit. 5. Winter cold alone can limit the predicted range shifts of these species to 200–300 km under 3.6–5.8 °C increase scenarios, which is less than half the value of earlier estimates.     

Thermal plasticity of photosynthesis: the role of acclimation in forest responses to a warming climate

The increasing air temperatures central to climate change predictions have the potential to alter forest ecosystem function and structure by exceeding temperatures optimal for carbon gain. Such changes are projected to threaten survival of sensitive species, leading to local extinctions, range migrations, and altered forest composition. This study investigated photosynthetic sensitivity to temperature and the potential for acclimation in relation to the climatic provenance of five species of deciduous trees, Liquidambar styraciflua, Quercus rubra, Quercus falcata, Betula alleghaniensis, and Populus grandidentata. Open‐top chambers supplied three levels of warming (+0, +2, and +4 °C above ambient) over 3 years, tracking natural temperature variability. Optimal temperature for CO₂ assimilation was strongly correlated with daytime temperature in all treatments, but assimilation rates at those optima were comparable. Adjustment of thermal optima was confirmed in all species, whether temperatures varied with season or treatment, and regardless of climate in the species' range or provenance of the plant material. Temperature optima from 17° to 34° were observed. Across species, acclimation potentials varied from 0.55 °C to 1.07 °C per degree change in daytime temperature. Responses to the temperature manipulation were not different from the seasonal acclimation observed in mature indigenous trees, suggesting that photosynthetic responses should not be modeled using static temperature functions, but should incorporate an adjustment to account for acclimation. The high degree of homeostasis observed indicates that direct impacts of climatic warming on forest productivity, species survival, and range limits may be less than predicted by existing models.     

Some like it hot! Rapid climate change promotes changes in distribution ranges of Nezara viridula and Nezara antennata in Japan

We compared past and current limits of the distribution range of the southern green stink bug, Nezara viridula (L.) (Heteroptera: Pentatomidae), in central Japan. In the early 1960s, the northern limit of the range was in Wakayama Prefecture and was limited by a +5 °C isothermal line for the mean January temperature. In 2006–2007, a new survey demonstrated that this northern limit had shifted northwards by 85 km (i.e., at a mean rate of 19.0 km per decade). The shift was most likely promoted by milder winter conditions. The mean January to February temperature in the region was 1.03–1.91 °C higher in 1998–2007 than in 1960–1969. The number of cold days (with the mean temperature below +5 °C) also significantly decreased, while the annual lowest temperature significantly increased. Nezara viridula was found mostly close to those locations where (i) the mean January temperature exceeded +5 °C, (ii) the mean number of cold days did not exceed 26 in January to February, and (iii) where the mean annual lowest temperature did not drop below –3.0 °C. The general linear model shows that the mean January temperature and number of cold days are the most important factors controlling the northern limit of distribution of N. viridula. All the climatic data suggest that over the last 45 years, environmental conditions have become more favourable for overwintering of N. viridula at many locations in central Japan. This has probably promoted the northward spread of the species, representing a direct response to climate warming. A sympatrically distributed congeneric, Nezara antennata Scott, seems to respond to the warming by a retreat from the ocean coast towards cooler elevated habitats, which might be a complex response to elevated temperature and interspecific mating with N. viridula. It is suggested that the range changes in both species will continue in response to further climate change.