Skip to main content

Impact of global warming scenarios on life-history traits of Tetranychus evansi (Acari: Tetranychidae)



The tomato red spider mite, Tetranychus evansi Baker & Pritchard (Acari: Tetranychidae), is an agricultural pest of solanaceous crops. Although T. evansi is of South American subtropical origin, it has recently expanded its distribution range to many tropical and temperate areas around the world. Its potential distribution range in response to scenarios of global warming was recently modeled, confirming its current and possible future distributions. Here, we experimentally investigated the biological traits of T. evansi in the context of the current and future global warming (2100) scenarios. Using an environmental simulation system, we tested the life-history traits of T. evansi under current summer temperatures (as of June, July, and August 2016) and under expected temperature increases based on two IPCC scenarios: RCP2.6 (+ 1 °C) and RCP8.5 (+ 3.7 °C). The mites were introduced into each scenario on 1 June and their sequential progeny were used for testing in each following month.


The mite could develop and reproduce under all scenarios. There was a decrease in the duration of lifespan and female fecundity at RCP8.5 during June and August, but this may be compensated for by the high intrinsic rate of increase, which implies faster population growth and shorter generation time.


Our study and other reports reveal the high adaptability of T. evansi to a wide range of summer temperatures; this may explain its current distribution. We anticipate that global warming will favor the spread of T. evansi and may further expand its distribution to a large area of the globe. These findings should be of ecological and practical relevance for designing prevention and control strategies.


Global warming is expected to have measurable impacts on terrestrial organisms [1,2,3,4]. According to the Intergovernmental Panel on Climate Change (IPCC) and Representative Concentration Pathways (RCPs) for greenhouse gas concentration scenarios, the global mean surface temperature is projected to rise by 0.3 to 4.8 °C by the end of the 21st century [5]. Biological consequences of global warming may include invasion and extinction rate acceleration; interruption of biological timing (reproduction, migration, growing season length, outbreaks, and distributions); disruption of ecological interactions (phenology, food sources, and predators); and alteration of ecosystem composition and functions [1, 6,7,8]. However, the magnitude of these effects will depend largely on a species’ ability to adapt and cope with global warming or to migrate into more suitable environments [6, 9]. Therefore, understanding the way organisms may adjust to global warming is a key challenge to research into climate change and biological invasion. To visualize the impacts of, and species adjustments to, global warming, accumulation of detailed ecological and biological data are important for predicting potential problems and formulating appropriate countermeasures to conserve rare species and protect against invasive ones.

The tomato red spider mite, Tetranychus evansi Baker & Pritchard (Acari: Tetranychidae), was reported for the first time in Brazil by Silva [10]. At the time, the species was misidentified as Tetranychus marianae McGregor, but it was later described as T. evansi by Baker and Pritchard [11]. For a summary of the taxonomic status of T. evansi, see Navajas et al. [12]. Tetranychus evansi is a pest that is destructive to several economically important members of the family Solanaceae, including tomato, eggplant, potato, tobacco, and other nightshade species [13,14,15]. Recently, T. evansi emerged as an invasive pest distributed in the tropical and temperate areas of nearly 43 countries; it is associated with 136 host plants of 36 plant families [16,17,18,19]. Although T. evansi is not a serious pest in its native habitat, it causes severe losses—sometimes reaching 100%—in tomato and other solanaceous plants, and it also disrupts the community composition of other Tetranychus species in invaded areas [18, 20, 21]. Tetranychus evansi does not enter diapause and is able to produce throughout the year if environmental conditions are favorable [13, 14]. Characteristics such as a higher intrinsic rate of increase (particularly at high temperatures: the thermal optimum is ~ 35 °C) than other Tetranychus species, absence of effective biological control agents, resistance to pesticides, and the ability to manipulate plant defenses account for the high invasive potential of T. evansi [15, 20, 22,23,24,25,26].

In Japan, T. evansi was first reported in Osaka in 2001. Its distribution has now expanded throughout the country, particularly in areas with temperate or tropical climates [18, 23, 27]. According to Boubou et al. [28, 29], two distinct genetic lineages of T. evansi have colonized areas outside their region of origin: Lineage 1 has the most invasive potential and has been recorded in several countries around the world, whereas lineage 2 has been found only in some parts of southern Europe. Accordingly, the population of T. evansi in East Asia (including Japan) belongs to lineage 1. This lineage has broader adaptability to the climatic gradient in temperate ecosystems and can exploit a wider range of host plants than lineage 2. It might have been accidentally introduced into Japan from neighboring countries in Asia (e.g., Taiwan) or from African or European countries [12, 17, 18, 30,31,32].

Predicting the geographical distribution range of a species such as T. evansi, with its strong invasive abilities, is essential [18]. Using available data, Meynard et al. [31] statistically predicted the current distributions of T. evansi in its native range and its invaded range, as well as its expected future expansion range in response to climate change. However, to date, there has been no conclusive, solid, and empirical evidence demonstrating the likely impacts of global warming on the biological activities of T. evansi. Assessing population growth, survival, reproduction, and the rate of increase, as influenced by scenarios of global warming, on an individual species basis could be a valid approach [2]. Here, we experimentally examined the life-history traits of T. evansi under three temperature scenarios: the current temperature, plus two global warming scenarios according to the IPCC [5], namely the mean values for “a stringent mitigation scenario” (RCP2.6) and “high greenhouse gas emissions” (RCP8.5) in 2100. For the current temperature scenario, Tokyo summer (June, July, and August) temperatures of 2016 were used. For the future scenarios, current temperatures were increased by 1 °C to represent RCP2.6 or by 3.7 °C to represent RCP8.5. We selected summer because this is the season of greatest population expansion of several mite species in the field [33,34,35,36,37].


Accuracy of temperature simulation

The set temperatures values under the current and 2100 IPCC scenarios are shown in Fig. 1. More than 98% of the variance in temperature set values could be explained by changes in the measured values (Fig. 2).

Fig. 1
figure 1

Simulated global warming scenarios. “Current,” natural temperature measured at 10-min interval (Japan Meteorological Agency); RCP2.6, current temperature increased by 1 °C; RCP8.5, current temperature increased by 3.7 °C. Mites were introduced to these conditions on 1 June and their sequential progeny were used in each subsequent month

Fig. 2
figure 2

Pearson correlation coefficient (r) of the relationships between set and measured values of temperature in the simulation system under the three scenarios

Relative humidity

Relative humidity was not controlled (Additional file 1: Figure S1). The measured relative humidity during the experiment was ~ 80%, ~ 70%, and ~ 65% under the current, RCP2.6, and RCP8.5 scenarios, respectively.

Development and lifespan

Two-way ANOVA for the impact of month and RCP scenario on immature development time, lifespan, and reproduction (eggs/female) showed significant main effects and interactions of both factors (Table 1). Different temperature scenarios and months did not have significant effects on hatchability percentage (% hatch) or survival to adulthood (% survival) of immature T. evansi2, P > 0.05, Table 2).

Table 1 Two-way ANOVA of factors (months and Representative Concentration Pathway [RCP] scenarios) affecting immature development time (egg-to-adult), lifespan, and reproduction of Tetranychus evansi
Table 2 Effects of simulated global warming scenarios on number of days to development at each stage, survival rate to adult, and adult longevity and lifespan of Tetranychus evansi

There were significant differences in the durations of the immature stages and in adult longevity and lifespan among different temperature scenarios during June (eggs: F2, 175 = 621.7; larvae: F2, 167 = 17.93; protonymphs: F2, 159 = 43.66; deutonymphs: F2, 155 = 64.03; egg-to-adult: F2, 155 = 497.4; longevity: F2, 129 = 14.95; lifespan: F2, 149 = 17.91; Pall < 0.0001) (Table 2). Post-hoc comparisons revealed no significant differences between the current and RCP2.6 scenarios in the lengths of the egg and larval periods and in adult longevity and lifespan (P > 0.05), but significant differences were found between these two scenarios for the duration of the protonymph (P < 0.0001) and deutonymph (P = 0.0055) stages and for total egg-to-adult period duration (P < 0.0001) (Table 2). The durations of the immature stages, as well as adult longevity and lifespan, of mites reared under the current scenario were all significantly different from those of mites reared under the RCP8.5 scenario (P < 0.0001). The immature development and adult longevity and lifespan were generally longer in mites reared under RCP2.6 than in those reared under RCP8.5 (protonymphs, P = 0.0343; others, P < 0.0001).

During July there were significant differences among different temperature scenarios in the durations of the egg incubation period (F2, 180 = 140.4, P < 0.0001), the larval period (F2, 180 = 13.26, P < 0.0001), the protonymph stage (F2, 177 = 34.5, P < 0.0001), and the egg-to-adult period (F2, 173 = 229.8, P < 0.0001), but not in the length of the deutonymph stage (F2, 173 = 1.466, P = 0.234). Adult longevity differed significantly among scenarios (F2, 162 = 3.316, P = 0.0388) but lifespan did not (F2, 171 = 1.656, P = 0.194). Post-hoc comparison showed that, with the exception of deutonymph period duration and lifespan (P > 0.05), all other developmental periods differed significantly among different temperature scenarios in July (P < 0.05).

In August there were significant differences among temperature scenarios in the durations of the immature stages and the adult longevity and lifespan of T. evansi (eggs: F2, 187 = 37.98; larvae: F2, 180 = 32.47; protonymphs: F2, 162 = 19.24; deutonymphs: F2, 156 = 32.06; egg-to-adult: F2, 156 = 265.6; longevity: F2, 132 = 27.36; lifespan: F2, 152 = 31.01; Pall < 0.0001). The larval, deutonymph, and egg-to-adult periods, but not the egg incubation period (P = 0.0571), were significantly longer in mites reared at the current temperature than in those reared at RCP2.6 (Pall < 0.0001), whereas the protonymph stage was significantly shorter in the former. The immature stages (except protonymphs) of mites reared under the current temperature were significantly longer than those of mites reared at RCP8.5. The egg incubation, protonymph, and egg-to-adult periods of mites reared at RCP8.5 were significantly shorter than those of mites reared at RCP2.6 (P < 0.0001). Adult longevity and lifespan did not differ significantly between the current and RCP2.6 scenarios (P > 0.05), but both were significantly greater under these two scenarios than under RCP8.5 (P < 0.0001). Male egg-to-adult duration was slightly shorter than that of females under all test conditions (data not shown).

Reproductive phases, fecundity, and sex ratio

The durations of the reproductive phases (pre-oviposition period [PrOP], oviposition period [OP], and post-oviposition period [PsOP]), as well as fecundity (total eggs/female [TEF]), varied among different temperature scenarios (Table 3). In June there were significant differences among temperature scenarios in the durations of the reproductive phases and in fecundity (PrOP: F2, 140 = 7.148, P = 0.0011; OP: F2, 129 = 15.4, P < 0.0001; PsOP: F2, 129 = 5.931, P = 0.0034; TEF: F2, 128 = 13.31, P < 0.0001). Duration of the oviposition and post-oviposition periods, as well as total eggs/female, did not differ significantly between mites under the current and RCP2.6 scenarios during June (P > 0.05), whereas the pre-oviposition period was significantly longer under current conditions (P = 0.0182). During June, mites reared under RCP8.5 had a significantly shorter oviposition period and fewer total eggs/female than those under the current and RCP2.6 scenarios (P < 0.0001). There were no significant differences in sex ratio among different temperature scenarios in June (χ2, P > 0.05).

Table 3 Effects of simulated global warming scenarios on length of reproductive phases, fecundity, and sex ratio of Tetranychus evansi

In July, there were significant differences among temperature scenarios in the durations of the pre-oviposition and oviposition periods, as well as in total eggs/female, but not in the post-oviposition period (PrOP: F2, 170 = 12.85, P < 0.0001; OP: F2, 162 = 5.052, P = 0.0074; PsOP: F2, 162 = 1.101, P = 0.335; TEF: F2, 162 = 7.182, P = 0.0010). In July, mites under RCP2.6 or RCP8.5 had pre-oviposition periods significantly shorter than those under the current temperature; the oviposition period under RCP2.6 was significantly longer than that under the current scenario but similar to that under RCP8.5, and total eggs/female under RCP2.6 and RCP8.5 was significantly greater than under the current scenario. There were no significant differences in sex ratio among different temperature scenarios in July (χ2, P > 0.05).

In a pattern similar to that in July, there were significant differences among temperature scenarios in the durations of the pre-oviposition and oviposition periods and in total eggs/female, but not in the duration of the post-oviposition period (PrOP: F2, 152 = 3.785, P = 0.0249; OP: F2, 132 = 28.11, P < 0.0001; PsOP: F2, 132 = 1.89, P = 0.155; TEF: F2, 132 = 37.13, P < 0.0001). The oviposition period was significantly shorter, and total eggs/female significantly fewer, under RCP8.5 than under the other two temperature scenarios (P < 0.0001). The ratio of females to males under RCP8.5 was slightly smaller than those under the current and RCP2.6 scenarios, but the differences were not significant (χ2, P > 0.05).

Life tables

A steeper decline in the age-specific survival rate (lx) curve was observed in mites under RCP8.5 in June and August (Fig. 3) owing to higher mortalities among adults in both months. In June and August the declines under the current and RCP2.6 scenarios were similar. In July, the lx curves all declined in a similar manner. The age-specific fecundity (mx, females/female/day) curves revealed that the onset of reproduction and the reproduction peak approached earlier in mites reared under RCP8.5 than in mites reared under RCP2.6 or the current scenario, regardless of the month. The highest peak in June occurred under RCP2.6 (7.8 females/female/day), in July under RCP2.6 (11.2 females/female/day), and in August under the current scenario (10 females/female/day).

Fig. 3
figure 3

Age-specific survival rate (lx) and age-specific fecundity (mx) of Tetranychus evansi under different global warming scenarios

In June, there were significant differences among temperature scenarios of current and RCP8.5 in all life table parameters (Table 4). The mean values of R0 (41.5), and T (17.8) were lower under RCP8.5 and differed significantly from those under RCP2.6. The intrinsic rate of increase and finite rate of increase under RCP8.5 were significantly higher than those under current scenario. To a similar extent, in July there were significant differences among temperature scenarios in the life table values; the means of the life-table values under RCP2.6 and RCP8.5 differed significantly from those under the current temperature. In August, there were significant differences among temperature scenarios of current and RCP2.6 in T (P < 0.0001), rm (P = 0.0085), λ (P = 0.0085) and Dt (P = 0.0087) but not R0 (P > 0.05).

Table 4 Effects of simulated global warming scenarios on life-table parameters of Tetranychus evansi


The tomato red spider mite, T. evansi, is a serious agricultural pest with high invasion potential and is currently spread over many countries worldwide. Unlike other Tetranychus species, T. evansi develops at a range of temperature between 12 and 45 °C; it has a high thermal optimum (~ 35 to 43 °C) and a high rate of population increase (rm = ~ 0.4 at 35 °C) [15, 22, 23]. These features enhance the ability of T. evansi to invade and adapt to a variety of geographical areas with varying climates and to become an emerging agricultural pest [12, 31]. It is therefore not surprising that we found here that T. evansi developed and reproduced even under the 2100 RCP8.5 scenario, in which the temperature was 3.7 °C above the current temperature and the maximum temperature experienced by the mites was 41.1 °C, during August. In general, development was favored by high temperatures (i.e., faster development with increased temperature), but reproduction was still high within moderate temperature ranges (i.e., under the current and RCP2.6 scenarios).

Development and reproduction

The egg hatchability of T. evansi decreased slightly (but not significantly across scenarios within months) with increasing temperature and was as low as 91.3% at RCP8.5 in August. Different strains of T. evansi reared at a constant range of temperatures (15 to 35 °C) have shown 95% to 99% hatchability [23]. Our results suggested that the hatchability, as well as several life-history traits, of T. evansi exposed to natural fluctuations in temperature to high levels in August (peaking at more than 40 °C under RCP8.5) might be hindered. In a manner similar to other mite species, T. evansi generally developed faster as the temperature increased. Under the current temperature scenario, the egg-to-adult period was 16.1 days in June, 11.0 days in July, and 9.2 days in August, when the monthly average temperatures were 22.4, 25.4, and 27.1 °C, respectively. The egg-to-adult period was ~ 1 day shorter under RCP2.6 (+ 1 °C) than in the current scenario. Under RCP8.5 (+ 3.7 °C), the egg-to-adult period was 2 to 4 days shorter than in the current scenario and 1 to 3 days shorter than under RCP2.6. Under semi-field conditions in Mauritius, Moutia [13] reported that the egg-to-adult period of T. evansi was about 6.5 days in summer (mean temperature 22.8 °C) and 18.5 days in winter (mean temperature 19.4 °C). The value of 6.5 days [13] is much shorter than the values we obtained here—for example, under the current scenario in June, when the average temperature was ~ 22 °C. This difference undoubtedly is a result of differences in the diurnal temperature fluctuation patterns between the tropical climate of Mauritius (20.3484° S, 57.5522° E) and the temperate to humid subtropical climate of Tokyo (35.6895° N, 139.6917° E). Our T. evansi males developed slightly faster than the females (data not shown), and this trend was common among all conditions tested. This phenology is common among Tetranychus species and has been reported before in T. evansi [15, 23]. The immature survival rate (i.e., rate of survival to adulthood) ranged from 80% to 97%, with no significant effects of different temperature scenarios within the same month. To a similar extent, different strains of T. evansi reared at temperatures between 15 and 35 °C have survival rates between 88 and 98%, with no significant differences [23].

The lifespan of mites reared at RCP8.5 in June was about two-thirds those of mites reared under the other scenarios, and the total eggs/female was about half that under the current and RCP2.6 scenarios. This trend was reversed during July, in which month the lifespan was similar among the different temperature scenarios. Interestingly, total eggs/female during July was significantly greater for females at RCP2.6 (~ 106 eggs/female) and RCP8.5 (~ 109 eggs/female) than under the current scenario (~ 78 eggs/female). In accordance with our experimental design, whereby the mites were introduced to each climate change scenario at the beginning of June and their offspring were used in the experiment the following month, the mites under RCP2.6 and RCP8.5 may have had the advantage of being adapted to temperature increases. Nevertheless, with the further increase in temperature during August under RCP8.5 the mites suffered a huge reduction in their lifespans and reproduction as compared with under other scenarios (see Tables 2, 3, 4). This can be explained by the maximum thermal limit for T. evansi, which might have been reached under RCP8.5 in August, resulting in large mortality rates among adults (Fig. 3). The percentage of females among the offspring ranged from 79% to 92% (see Table 3). A female-biased sex ratio is common among different strains of T. evansi, and to some extent it is higher than in other Tetranychus species [14, 15, 22, 23].

Life tables

Life tables are very convenient tools for assessing environmental effects on population development, survival, and reproduction [38] The net reproductive rate (females/female/generation) was highest under current and RCP2.6 in June and August, but it differed significantly from that in the current scenario in July and from that under RCP8.5 only in June and August (see Table 4). Generation time (days) was shortest under RCP8.5 and ranged from 17.8 days in June to 12.9 days in August. The intrinsic rate of population increase was, in general, higher at increased temperature owing to the short development time and early peak reproduction. This might also explain the wide temperature range at which T. evansi is able to develop and cause economic crop losses. The finite rate of increase (females/female/day) displayed a trend similar to rm. Population doubling times were similar between RCP2.6 and RCP8.5 temperature scenarios in June, significantly shorter under the two scenarios in July, and similar between current and RCP8.5 and between RCP2.6 and RCP8.5 in August. Because of the fluctuating nature of temperature regimens that mimic natural conditions, the life-history and life-table trait measurements that we obtained here are, to some extent, different from those in previous reports that have relied solely on laboratory studies at constant temperatures [15, 22, 23]. Environmental factors other than temperature might be involved (e.g., host plant, see Murungi et al. [39]). However, these differences among results further highlight the importance of and the need for gathering biological measurements under conditions that mimic the natural diurnal fluctuations in environmental factors such as temperature.


Invasive alien pest species pose great challenges to world agriculture, biodiversity, and ecosystems. Empirical evidence of the likely impact of global warming on agricultural pests is essential. Our findings indicated that T. evansi is able to rapidly adapt to the increases in temperature; therefore, with long-term adaptation, the mite will be able to spread more widely and to a broad range of environmental temperatures. As well as possessing the advantages of a high thermal optimum, short generation time, and high rate of population growth, T. evansi—as is becoming evident in other organisms—might have already undergone what is known as “micro-evolutionary change in situ”; through its journey from its native range to almost everywhere around the world this species may already have encountered large environmental variations [1, 6, 7, 40, 41]. This speculation is supported by the fact that T. evansi has already spread to a geographical area that has climatic gradients different from those in its native habit [12]. Our findings and previous modeling results [31] are in agreement that the greatest T. evansi threat will be shifted to higher latitudes. However, it might be difficult to generalize our findings to other pest mites, because development and reproduction temperatures are species-specific and have to be assessed on an individual species basis.

Materials and methods

Mite culture

The founder population of T. evansi was originally collected from black nightshade plant (Solanum nigrum L.) in Tokyo (Japan, 35°35′ N, 139°36′ E) in November 2006 and was maintained on black nightshade at 25 ± 1 °C under a 16:8-h (light:dark) photoperiod at the Laboratory of Applied Entomology and Zoology, Faculty of Agriculture, Ibaraki University, Japan [17, 24]. A colony of T. evansi was brought to our laboratory at Tokyo University of Agriculture and Technology in May 2017 and thereafter maintained on detached leaves of eggplant (Solanum melongena L. cv. Senryo #2) at 25 ± 1 °C with a 16:8-h photoperiod. The mites were maintained on eggplant for about 3 months before they were used in the experiments.

Climate data

Temperature data for Tokyo in 2016 (June, July, and August) in 10-min measurement intervals were obtained from the Japan Meteorological Agency ( To examine the effects of predicted global warming, the average temperature increase values in two scenarios defined by the IPCC [5] were used. Under RCP2.6 the temperature was projected to increase by 1 °C, and under RCP8.5 by 3.7 °C, on average. An environmental simulation system (ESS) was used to simulate temperature data and create the three scenarios. ESS is a computer-based closed system for simulating natural climate conditions [37, 42, 43]. Briefly, the system software created a schedule of 10-min intervals with corresponding set values of air temperature (set temperature, TSVs) and continuously measured the process values of air temperature (measured temperature, TPVs). TPV was adjusted to TSV by switching on or off an air heater (TSR210-A; Tescom, Tokyo, Japan) and a refrigerator (JF-NU40B; Haier Japan Sales, Osaka, Japan) every 10 min. The daily natural photoperiod was also created by the system by turning on or off light-emitting diodes. Relative humidity was not simulated but was recorded throughout the experiments (Additional file 1: Figure S1).

Development and reproduction

Three units of the ESS were used simultaneously. Each unit was assigned to one of the temperature scenarios. A random-aged egg-laying adult females of T. evansi were introduced onto a detached eggplant leaf (ca. 40 × 40 mm) placed upside down on a water-saturated cotton pad in a polystyrene square Petri dish (100 × 15 mm). Then, on 1 June, each dish was transferred to one of the ESS units; the females were left to lay eggs for 24 h and then removed. The eggs were categorized into two groups: one group was kept on the same detached leaf until hatching and used to calculate hatchability percentages. The other group (approximately 80 eggs) was transferred separately onto eggplant leaf discs (1 cm in diameter) and observed daily throughout the mites’ whole lifespan. Leaf discs were replaced every 3 days or when necessary. Upon adult female emergence, one adult male was introduced for copulation. To determine the sex ratio, 2 or 3 days after the onset of oviposition, the eggs laid by females in each treatment during 2 consecutive days were collected and transferred to a new detached eggplant leaf. The collected eggs were kept to develop to adulthood and then sexed. These adults were maintained under the same conditions and used to produce eggs for the consecutive-month experiment. For example, the June progeny under the current scenario were transferred to a new detached leaf on 1 July under the same scenario and allowed to lay eggs for 24 h; the adults were then discarded and the eggs treated as previously. The same observations and experimental procedures were applied to other months and scenarios.

Statistical analysis

Pearson’s correlation coefficient was used to examine the relationship between TSVs and TPVs as evidence of simulation accuracy. Percentages (hatchability, survival to adulthood, and sex ratio) were compared by using a Chi-squared test with Bonferroni correction. Two-way analysis of variance (ANOVA) was used to test the effects of month and temperature scenario on development time (egg-to-adult), lifespan, and reproduction (eggs/female). One-way ANOVA with Tukey’s post hoc test was used for means separation with α = 0.05. The development and reproduction data were square-root transformed to meet assumptions of normality when necessary and retransformed for representation purposes. The data for development, survival rate, longevity and female daily fecundity were analyzed according to the age-stage, two-sex life table approach using the computer program TWOSEX-MSChart [44,45,46]. The following parameters were calculated as described by Chi & Liu [44]: age-stage-specific survival rate (sxj, the probability that a newly laid egg will survive to age x and stage j), age-specific survival rate (lx, the percentage of females alive at age x), age-specific fecundity (mx, the number of female offspring produced by a female in a unit of time), age-stage specific fecundity (fxj, the mean fecundity of females at age x), intrinsic rate of natural increase (r), finite rate of increase (λ, the number of times the population multiplies in a unit of time), net reproductive rate (R0, number of female offspring/female/generation), mean generation time (T, mean age of mothers at time of birth of female offspring) and population doubling time (Dt, the time needed for the population to double). The age-specific survival rate (lx) and age-specific fecundity (mx) are calculated as:

$$l_{x} = \mathop \sum \limits_{j = 1}^{k} s_{xj}$$
$$m_{x} = \frac{{\mathop \sum \nolimits_{j = 1}^{k} s_{xj} f_{xj} }}{{\mathop \sum \nolimits_{j = 1}^{k} s_{xj} }}$$

where k is the number of life stages. The reproductive rate (R0) is calculated as follows:

$$R_{0} = \mathop \sum \limits_{x = 0}^{\infty } l_{x} m_{x}$$

The intrinsic rate of population increase (rm) was calculated by using the Euler–Lotka equation as follows:

$$\mathop \sum \limits_{x = 0}^{\infty } e^{{ - r_{m} \left( {x + 1} \right)}} l_{x} m_{x} = 1$$

where x is a female age in days. The finite rate of increase (λ), mean generation time (T) and population doubling time (Dt) are calculated as follows:

$$T = \ln \frac{{R_{0} }}{{r_{m} }},$$
$$\lambda = e^{{r_{m} }} ,$$
$$D_{t} = \frac{\ln 2}{{r_{m} }}.$$

The bootstrap procedure (B = 100,000) was used to estimate means and standard errors for the life tables parameters and the paired bootstrap test was used to compare different treatment based on the confidence interval of the differences using TWOSEX-MSChart [47,48,49]. ANOVAs were performed in R v. 3.4.0 software [50]. Pearson’s correlation coefficient was determined by using scipy.stats.pearsonr in Python [51]. Data were visualized with Matplotlib in Python [52].

Availability of data and materials

The data used and/or analyzed during the current study are available from the corresponding authors on reasonable request.



environmental simulation system




intergovernmental panel on climate change

l x :

age-specific survival rate

m x :

age-specific fecundity


oviposition period


pre-oviposition period


post-oviposition period

R 0 :

net reproductive rate

T :

generation time

r m :

intrinsic rate of population increase


representative concentration pathway


total eggs/female


measured temperature


set temperature

λ :

finite rate of population increase

D t :

population doubling time


  1. Bale JS, Masters GJ, Hodkinson ID, Awmack C, Bezemer TM, Brown VK, Butterfield J, Buse A, Coulson JC, Farrar J, Good JE. Herbivory in global climate change research: direct effects of rising temperature on insect herbivores. Glob Change Biol. 2002;8:1–16.

    Article  Google Scholar 

  2. Logan JA, Regniere J, Powell JA. Assessing the impacts of global warming on forest pest dynamics. Front Ecol Environ. 2003;1:130–7.

    Article  Google Scholar 

  3. Deutsch CA, Tewksbury JJ, Huey RB, Sheldon KS, Ghalambor CK, Haak DC, Martin PR. Impacts of climate warming on terrestrial ectotherms across latitude. PNAS. 2008;105:6668–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Macfadyen S, McDonald G, Hill MP. From species distributions to climate change adaptation: knowledge gaps in managing invertebrate pests in broad-acre grain crops. Agric Ecosyst Environ. 2018;253:208–19.

    Article  Google Scholar 

  5. Pachauri RK, Allen MR, Barros VR, Broome J, Cramer W, Christ R, et al. Climate change 2014: synthesis report. Contribution of Working Groups I, II and III to the fifth assessment report of the Intergovernmental Panel on Climate Change. Geneva: IPCC; 2014.

    Google Scholar 

  6. Hughes L. Biological consequences of global warming: is the signal already apparent? Trends Ecol Evol. 2000;15:56–61.

    Article  CAS  PubMed  Google Scholar 

  7. Ladányi M, Horváth L. A review of the potential climate change impact on insect populations-general and agricultural aspects. Appl Ecol Environ Res. 2010;8:143–52.

    Article  Google Scholar 

  8. Buckley LB, Arakaki AJ, Cannistra AF, Kharouba HM, Kingsolver JG. Insect development, thermal plasticity and fitness implications in changing, seasonal environments. Integr Comp Biol. 2017;57:988–98.

    Article  PubMed  Google Scholar 

  9. Arribas P, Abellán P, Velasco J, Millán A, Sánchez-Fernández D. Conservation of insects in the face of global climate change. In: Johnson N, Jones TH, editors. Global climate change and terrestrial invertebrates. Chichester: Wiley; 2017. p. 349–67.

    Google Scholar 

  10. Silva P. Um novo a´caro nocivo ao tomateiro na Bahia. Bol Inst Biol Bahia. 1954;1:1–20.

    Google Scholar 

  11. Baker EW, Pritchard AE. The tetranychoid mites of Africa. Hilgardia. 1960;29:455–574.

    Article  Google Scholar 

  12. Navajas M, De Moraes GJ, Auger P, Migeon A. Review of the invasion of Tetranychus evansi: biology, colonization pathways, potential expansion and prospects for biological control. Exp Appl Acarol. 2013;59:43–65.

    Article  PubMed  Google Scholar 

  13. Moutia LA. Contribution to the study of some phytophagous Acarina and their predators in Mauritius. Bull Entomol Res. 1958;49:59–75.

    Article  Google Scholar 

  14. Qureshi AH, Oatman ER, Fleschner CA. Biology of the spider mite, Tetranychus evansi. Ann Entomol Soc Am. 1969;62:898–903.

    Article  Google Scholar 

  15. De Moraes GJ, McMurtry JA. Effect of temperature and sperm supply on the reproductive potential of Tetranychus evansi (Acari: Tetranychidae). Exp Appl Acarol. 1987;3:95–107.

    Article  Google Scholar 

  16. Migeon A, Dorkeld F. Spider Mites Web: a comprehensive database for the Tetranychidae. 2006–2017. Accessed 5 Mar 2018.

  17. Gotoh T, Araki R, Boubou A, Migeon A, Ferragut F, Navajas M. Evidence of co-specificity between Tetranychus evansi and Tetranychus takafujii (Acari: Prostigmata, Tetranychidae): comments on taxonomic and agricultural aspects. Int J Acarol. 2009;35:485–501.

    Article  Google Scholar 

  18. Migeon A, Ferragut F, Escudero-Colomar LA, Fiaboe K, Knapp M, de Moraes GJ, Ueckermann E, Navajas M. Modelling the potential distribution of the invasive tomato red spider mite, Tetranychus evansi (Acari: Tetranychidae). Exp Appl Acarol. 2009;48:199–212.

    Article  PubMed  Google Scholar 

  19. Kazak C, Döker I, Karut K. First record of invasive tomato spider mite Tetranychus evansi (Acari: Tetranychidae) in Turkey. Int J Acarol. 2017;43:325–8.

    Article  Google Scholar 

  20. Ferragut F, Garzón-Luque E, Pekas A. The invasive spider mite Tetranychus evansi (Acari: Tetranychidae) alters community composition and host-plant use of native relatives. Exp Appl Acarol. 2013;60:321–41.

    Article  PubMed  Google Scholar 

  21. Furtado IP, De Moraes GJ, Kreiter S, Flechtmann CH, Tixier M, Knapp M. Plant inhabiting phytoseiid predators of midwestern Brazil, with emphasis on those associated with the tomato red spider mite, Tetranychus evansi (Acari: Phytoseiidae, Tetranychidae). Acarologia. 2014;54:425–31.

    Article  Google Scholar 

  22. Bonato O. The effect of temperature on life history parameters of Tetranychus evansi (Acari: Tetranychidae). Exp Appl Acarol. 1999;23:11–9.

    Article  Google Scholar 

  23. Gotoh T, Sugimoto N, Pallini A, Knapp M, Hernandez-Suarez E, Ferragut F, Ho CC, Migeon A, Navajas M, Nachman G. Reproductive performance of seven strains of the tomato red spider mite Tetranychus evansi (Acari: Tetranychidae) at five temperatures. Exp Appl Acarol. 2010;52:239–59.

    Article  CAS  PubMed  Google Scholar 

  24. Gotoh T, Fujiwara S, Kitashima Y. Susceptibility to acaricides in nine strains of the tomato red spider mite Tetranychus evansi (Acari: Tetranychidae). Int J Acarol. 2011;37:93–102.

    Article  Google Scholar 

  25. Nyoni BN, Gorman K, Mzilahowa T, Williamson MS, Navajas M, Field LM, Bass C. Pyrethroid resistance in the tomato red spider mite, Tetranychus evansi, is associated with mutation of the para-type sodium channel. Pest Manag Sci. 2011;67:891–7.

    Article  CAS  PubMed  Google Scholar 

  26. Schimmel BC, Ataide L, Chafi R, Villarroel CA, Alba JM, Schuurink RC, Kant MR. Overcompensation of herbivore reproduction through hyper-suppression of plant defenses in response to competition. New Phytol. 2017;214:1688–701.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Ehara S, Ohashi K. A new species of Tetranychus (Acari: Tetranychidae) from the Kinki district, Japan. Acta Arachnol. 2002;51:19–22.

    Article  Google Scholar 

  28. Boubou A, Migeon A, Roderick GK, Navajas M. Recent emergence and worldwide spread of the red tomato spider mite, Tetranychus evansi: genetic variation and multiple cryptic invasions. Biol Invasions. 2011;13:81–92.

    Article  Google Scholar 

  29. Boubou A, Migeon A, Roderick GK, Auger P, Cornuet JM, Magalhães S, Navajas M. Test of colonisation scenarios reveals complex invasion history of the red tomato spider mite Tetranychus evansi. PLoS ONE. 2012;7:e35601.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Migeon A, Auger P, Hufbauer R, Navajas M. Genetic traits leading to invasion: plasticity in cold hardiness explains current distribution of an invasive agricultural pest, Tetranychus evansi (Acari: Tetranychidae). Biol Invasions. 2015;17:2275–85.

    Article  Google Scholar 

  31. Meynard CN, Migeon A, Navajas M. Uncertainties in predicting species distributions under climate change: a case study using Tetranychus evansi (Acari: Tetranychidae), a widespread agricultural pest. PLoS ONE. 2013;8:e66445.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Santamaría ME, Auger P, Martínez M, Migeon A, Castañera P, Díaz I, Navajas M, Ortego F. Host plant use by two distinct lineages of the tomato red spider mite, Tetranychus evansi, differing in their distribution range. J Pest Sci. 2018;91:169–79.

    Article  Google Scholar 

  33. Takafuji A, Kamibayashi M. Life cycle of a non-diapausing population of the two-spotted spider mite, Tetranychus urticae Koch in a pear orchard. Res Popul Ecol. 1984;26:113–23.

    Article  Google Scholar 

  34. Montserrat M, Sahún RM, Guzmán C. Can climate change jeopardize predator control of invasive herbivore species? A case study in avocado agro-ecosystems in Spain. Exp Appl Acarol. 2013;59:27–42.

    Article  PubMed  Google Scholar 

  35. Montserrat M, Guzmán C, Sahún RM, Belda JE, Hormaza JI. Pollen supply promotes, but high temperatures demote, predatory mite abundance in avocado orchards. Agric Ecosyst Environ. 2013;164:155–61.

    Article  Google Scholar 

  36. Ghazy NA, Osakabe M, Negm MW, Schausberger P, Gotoh T, Amano H. Phytoseiid mites under environmental stress. Biol Control. 2016;96:120–34.

    Article  Google Scholar 

  37. Ghazy NA, Suzuki T, Amano H. Development and reproduction of Neoseiulus californicus (Acari: Phytoseiidae) and Tetranychus urticae (Acari: Tetranychidae) under simulated natural temperature. Environ Entomol. 2018;47:1005–12.

    Article  PubMed  Google Scholar 

  38. Maia ADH, Luiz AJ, Campanhola C. Statistical inference on associated fertility life parameters using jackknife technique: computational aspects. J Econ Entomol. 2000;93:511–8.

    Article  Google Scholar 

  39. Murungi LK, Nyende A, Wesonga J, Masinde P, Knapp M. Effect of African nightshade species (Solanaceae) on developmental time and life table parameters of Tetranychus evansi (Acari: Tetranychidae). Exp Appl Acarol. 2010;52(1):19–27.

    Article  PubMed  Google Scholar 

  40. Hill JK, Griffiths HM, Thomas CD. Climate change and evolutionary adaptations at species’ range margins. Annu Rev Entomol. 2011;56:143–59.

    Article  CAS  PubMed  Google Scholar 

  41. Pauls SU, Nowak C, Bálint M, Pfenninger M. The impact of global climate change on genetic diversity within populations and species. Mol Ecol. 2013;22:925–46.

    Article  PubMed  Google Scholar 

  42. Nishide Y, Suzuki T, Tanaka S. Synchrony in the hatching of eggs in the desert locust Schistocerca gregaria (Orthoptera: Acrididae): egg condition influences hatching time in the laboratory and under simulated field temperatures. Appl Entomol Zool. 2017;52:599–604.

    Article  Google Scholar 

  43. Nishide Y, Suzuki T, Tanaka S. The hatching time of Locusta migratoria under outdoor conditions: role of temperature and adaptive significance. Physiol Entomol. 2017;42:146–55.

    Article  CAS  Google Scholar 

  44. Chi H, Liu H. Two new methods for the study of insect population. Bull Inst Zool Acad Sin. 1985;24:225–40.

    Google Scholar 

  45. Chi H. Life table analysis incorporating both sexes and variable development rates among individuals. Environ Entomol. 1988;17:26–34.

    Article  Google Scholar 

  46. Chi H. TWOSEX-MSChart: a computer program for the age-stage, two-sex life table analysis. Taichung: National Chung Hsing University; 2019. Accessed 20 Feb 2019.

  47. Chi H. Timing of control based on the stage structure of pest populations—a simulation approach. J Econ Entomol. 1990;83:1143–50.

    Article  Google Scholar 

  48. Efron B, Tibshirani R. An introduction to the bootstrap. New York: Chapman & Hall; 1985.

    Google Scholar 

  49. Akca I, Ayvaz T, Yazıcı E, Smith CL, Chi H. Demography and population projection of Aphis fabae (Hemiptera: Aphididae): with additional comments on life table research criteria. J Econ Entomol. 2015;108:1466–78.

    Article  PubMed  Google Scholar 

  50. R Core Team. R: a language and environment for statistical computing. R Found Stat Comput. 2016.

  51. Seabold S, Perktold J. Statsmodels: econometric and statistical modeling with python. In: Proceedings of the 9th python in science conference 2010 Jun 28, vol. 57. Scipy; p. 61.

  52. Hunter JD. Matplotlib: a 2D graphics environment. Comput Sci Eng. 2007;9:90–5.

    Article  Google Scholar 

Download references


We thank Masao Ohyama for his technical assistance in developing the simulation system used in this study.


This study is supported by the Japan Society for the Promotion of Science JSPS (Grants Numbers JP25450069, L17551). The funding sources had no roles in the design of the study, data collection and data analysis, interpretation of data or in the writing of the manuscript.

Author information

Authors and Affiliations



NAG and TS designed the study and performed the research. NAG, TS and TG analyzed the data. NAG and TS wrote the first draft of the manuscript. All authors contributed to revisions. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Noureldin Abuelfadl Ghazy or Takeshi Suzuki.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent to publish

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Figure S1.

Relative humidity (%) measurements made during the experiments.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ghazy, N.A., Gotoh, T. & Suzuki, T. Impact of global warming scenarios on life-history traits of Tetranychus evansi (Acari: Tetranychidae). BMC Ecol 19, 48 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: