1. Introduction
In all mite families, the Tetranychidae family, which encompasses > 1300 species of spider mites (Acari: Tetranychidae), holds a significant position regarding its impact on agricultural production [
1]. Specifically, more than 100 pest spider mite species could harm agricultural and forestry crops [
2]. The two-spotted mite (
Tetranychus urticae), a generalist herbivore that could feed on > 1150 host plant species (including a multitude of important crops such as corn, soybean, tomato, cucumber, and so on) [
1,
3] is the most widely distributed and destructive spider mite species worldwide. Furthermore, multiple
T. urticae populations could strongly resist major acaricides and insecticides in various regions [
4]. Since 1983,
T. urticae has been widespread in most regions in China and is currently the dominant spider mite species [
5,
6]. Specifically, in the Tibetian Autonomous Region of China,
T. urticae was reported to harm numerous vegetable crops and was also difficult to control [
7]. However, few studies on the biological management of spider mites have been conducted in Tibet, the world’s highest plateau.
Based on the drawbacks of insecticide or pesticide resistance and pollution, biological control is considered one of the most environmentally friendly and economical methods for managing insect pests [
8]. Predatory mites are one of the most important natural biological control agents. Specifically, phytoseiid mites are vital biological control agents, with
Amblyseius swirskii Athias-Henriot (Acari: Phytoseiidae) as the most successful predatory mite, having been employed in > 50 countries globally [
9,
10]. Other biological factors used for controlling the disastrous harm caused by two-spotted mites include predatory mite species such as
Neoseiulus californicus and many other mites of the genus
Amblyseius [
11,
12,
13].
With an average elevation of >4000 m above sea level, the Tibetan Plateau is the highest of China’s three major landform ladders. It is characterized by low temperatures and hypoxia, the most substantial threats to the lives of arthropods and other animals in the area [
14]. Invertebrates often develop various useful strategies to combat adverse highland conditions. For example, to combat hypoxia conditions (2% O
2), the bean weevil increases the synthesis of metabolites, including carbohydrates, amino acids, and organic acids [
15]. Furthermore, since the cytochrome oxidase activity is higher in Tibetan locusts, they show a higher hypoxia tolerance than lowland locusts [
16]. Additionally, as a means of surviving the high-latitude environment, the enemy insect
Bombus pyrosoma would upregulate energy metabolism genes [
17]. Despite these useful insights, there are no recent explorations on the effects of highland conditions on the fitness and prey capacity of predatory mites.
The predator
N. californicus was first reported by McGregor in 1954 in California, United States [
18,
19]. As a polyphagous predator,
N. californicus mainly fed on Tetranychid spider mites (especially
T. urticae, and other spider mites such as
T. evansi) [
20,
21,
22]. It could also consume various small insects, including other mite species or thrips; hence, McMurtry classified it as a type II specialist predatory mite (selective predators of Tetranychid mites) [
10].
On the other hand, predatory mites of the genus
Amblyseius, popular for feeding on various prey including a multitude of mites from many families (such as Tetranychidae, Eriophyidae, Pyroglyphidae, and so on) and many other arthropods including enormous pest insects (such as whiteflies, thrips, nematodes, and so on) [
10,
23], were always considered as generalist predators, hence were listed as subtype III-b in type III lifestyle predators by McMurtry [
23]. However, although the species
N. californicus and other predators also belong to the genus
Amblyseius, the impact of highland conditions on them has received little research attention, and their predatory capacity remains largely unknown.
Herein, four predatory mite species, including species N. californicus and three species from the genus Amblyseius (A. swirskii, A. barkeri, and A. cucumeris) were examined to determine the most suitable predatory mite species for controlling T. urticae in the Tibetan Plateau, the world’s highest plateau. The fitness and predatory capacities of the four species under highland conditions were assessed. The results showed that in the Tibetan Plateau, A. swirskii and A. barkeri were the most effective and least effective predatory mites in preying on T. urticae, respectively.
2. Materials and Methods
2.1. Establishment of the T. urticae and predatory mite populations
The
T. urticae population was collected from soybean leaves in the solar greenhouse of the Lhasa Agricultural Science and Technology Park (roughly E91°07′, N29°38′). Spider mite species identification was confirmed by COI sequencing [
24]. After introducing the spider mites into the laboratory, the isofemale line of
T. urticae mites was established and grown on the leaves of kidney beans (
Phaseolus vulgaris L.) in controlled conditions. The
T. urticae population was used as the target or objective pest mite to determine the prey capacity of the predatory mites.
Herein, four predatory mite species, including
A. swirskii,
Amblyseius barkeri (Hughes),
Amblyseius cucumeris (Oudemans), and
N. californicus were examined. All predatory mite species were purchased from Fujian Yanxuan Biological Control Technology Company (Fuzhou, China) and were categorized based on three characteristics: spermatheca, abdominal shield, and dorsal shield [
25].
2.2. Fitness experiment for the four predatory mite species
The fitness of all four predatory species was assessed through a life-table experiment after rearing in highland conditions in Lhasa, Tibet Autonomous Region, China (roughly E91°07′, N29°38′) for more than ten generations (approximately six months). The rearing conditions were as follows: 1. Temperatures of 25.5 ± 1°C; 2. Relative Humidity (RH) of 80 ± 5%; 3. A photoperiod of 16:8 h (L: D); and 4. A mean altitude of 3642 m above sea level. First, thirty new 1-day-old emerged female adults and thirty new 1-day-old emerged male adults were selected for each mite species to oviposit for 4 h in one dish (9 cm). Following that, 30 eggs of each species were randomly selected to conduct the subsequent test, and one egg was placed in one dish with a kidney bean leaf (about 7 cm) as a replicate, with each mite species having 30 replicates. In each dish, five T. urticae eggs and five T. urticae nymphs (mixed age stages) were added every five days as a diet. To calculate the juvenile survival and developmental rate of different predatory mite species, the stages of all predatory mites were monitored and recorded every 12 h. Ten predatory mites were regarded as a replicate for juvenile survival calculation. Each mite species experiment was performed in triplicate.
Once a female predatory mite emerged, one newly emerged male adult was introduced to mate with one female adult. Ten female adults of each mite species were randomly selected each day to record fecundity. Given that the pre-adult Survival Rate (SR) of A. barkeri was extremely low, ten 1-day-old female A. barkeri adults were selected from other leaf dishes. Fecundity and hatchability data were first tested for normality (Kolmogorov–Smirnov test) and then homogeneity of group variances (Levene’s test). Subsequently, one-way ANOVA analysis (SPSS 21.0) was employed to determine the variance of fecundity and the developmental rate among all four predatory mites. Given that the SRs did not follow a normal distribution, they were analyzed by a Kruskal–Wallis test, and multiple comparisons were performed using Dunn’s test with Bonferroni correction (SPSS v21.0). The Log-rank (Mantel-Cox) test was used to compare the survival curves of each predatory mite’s female adults (GraphPad Prism v9.0.2).
2.3. Functional response experiment
The functional response experiments on the ability of the four predatory mite species to prey on
T. urticae (different stages) were conducted under 25℃ conditions. First, an equivalent amount of
T. urticae at three variant stages (egg, nymph, and adult) was selected at six density gradients (1, 3, 5, 7, 9, and 11), and the same number of
T. urticae eggs, nymphs, and adults were all placed in a dish, while one 1-day-old female adult predatory mite was placed in the same dish after 24 h of starvation. After 24 h, the predation activities of the predatory mites on different stages of
T. urticae were recorded. The preying of predatory mites on different stages of
T. urticae was measured using a predator-prey model, specifically the Holling type II functional response [
26]:
Where Na = the number of T. urticae mites preyed on, ɑ = the attack rate (proportion of prey captured by each predator per unit of searching time), Th = the handling time (the time it takes for predatory mites to identify, kill, and consume T. urticae mites), N=the prey density, and T = the time it takes for predators to find the prey (in this test, T is 1 d).
2.4. Predator interference by different densities of predatory mites
The predatory abilities of different mites at various densities were measured. Herein, six predatory mite density gradients were used (1, 3, 5, 7, 9, and 11). Each mite species was prepared per each density, and these combinations were replicated three times. For every replicate, forty
T. urticae female adults were first put in a leaf dish. One selected density of 1-day-old female adults from a single predatory species was then placed in the same dish, and the predation situation was recorded after 24 h. Predator interference analysis was performed using Watt’s model of the effect of the densities of attacked and attacking species on the number of the attacked (Watt, 1959):
Where A = the number of T. urticae female adults preyed on, ɑ = the attack rate without competition, X = the density of predators per leaf dish, and b = the intraspecific competition parameter.
Data on the predatory capacity, handling time, attack coefficient, and max predatory amount were tested for normality (Kolmogorov–Smirnov test) and homogeneity of group variances (Levene’s test) if they followed a normal distribution and One-way ANOVA analysis (SPSS 21.0) was performed to determine the variance among all four predatory mites. On the other hand, the data were subjected to a Kruskal–Wallis test, and multiple comparisons were performed using Dunn’s test with Bonferroni correction if they followed a non-normal distribution (SPSS 21.0).
4. Discussion
Herein, the fitness and prey capacity of four predatory mite species (A. swirskii, A. barkeri, A. cucumeris, and N. californicus) were assessed using the life-table and functional response experiments. According to the results, A. swirskii had the highest fecundity and pre-adult SR, N. californicus took the least time in transitioning from the egg stage to emergence, and A. barkeri could not survive in highland conditions, as most of its nymphs died before emergence. Furthermore, based on the functional response results, A. swirskii had the highest prey capacity on T. urticae mites compared to the other three predatory species on many parameters. Our findings collectively implied that among the four tested groups, A. swirskii could be the best predatory mite species to biologically control two-spotted mites in the Tibetan Plateau.
The two-spotted mite (
T. urticae) is one of the world’s most disastrous agricultural arthropod pests, damaging a wide range of host plants, including many vital crops such as vegetables, cotton, maize, and so on [
27,
28]. In China, chemical acaricides and insecticides were the primary methods for controlling
T. urticae. However, over the past few decades,
T. urticae rapidly developed a global resistance to these chemical agents, especially in China [
29,
30,
31,
32]. In Tibet, the world’s largest highland, chemical pesticides were widely applied to combat pests, and the resulting significant residues became an issue of concern [
33]. Interestingly, no treatments other than chemical control were used in Tibet to effectively prevent the harm caused by two-spotted mites. Despite being important biological control agents against many pest species, including spider mites, predatory mites were rarely used in Tibet, and pertinent research was also extremely limited even in our studies; hence, the fitness and predatory capacity of many predatory mites on
T. urticae in this region remained unclear [
33,
34,
35]. Adding to the literature, this study comprehensively explored the fitness and predatory capacity of the four widely used predatory mite species in the Tibetan area.
Multiple studies have reported that high-attitude conditions could significantly influence arthropods’ development, morphological characteristics, and body size, among many other aspects [
36,
37,
38,
39,
40]. For example, a significant link was found between insects’ body size and altitude in a study on wing-reduced stonefly in New Zealand [
40]. Furthermore, insects’ fitness was found to decrease with increasing elevation, as demonstrated by the reduced fecundity of the willow leaf beetle (
Chrysomela aeneicollis) under high altitude conditions [
41]. Herein, we confirmed that a high attitude significantly adversely affected
A. barkeri, resulting in an extremely low SR before emergence (13.33% ± 5.77%), whereas the other three predator species had a better SR as they all had a relatively better fitness, especially the
A. swirskii female adult mites, which had the highest fecundity and pre-adult SR (
Figure 2 and
Figure 3).
It has been reported that
A. swirskii is one of the most successful biocontrol agents globally, as it can prey on many severe pests and is easy to rear [
9]. Furthermore, another study in China indicated that
A. swirskii had an advantage in the predation of whiteflies over other two predatory mite species;
Amblyseius orientalis and
Neoseiulus californicus [
42]. These results are consistent with our findings, as
A. swirskii showed a significantly higher predatory capacity on
T. urticae than
N. californicus at 15 d post-release (
Figure 9A), and its fitness in Tibet was comparable to that of similar strains reared under normal conditions [
43]. As for
N. californicus, our results indicated that although most fitness parameters including developmental time, pre-adult SR, and longevity in Tibet were comparable to those of mites reared in low altitude conditions [
44,
45,
46], its fecundity in Tibet (Mean ± SEM = 10.10 ± 0.43) was significantly lower compared to that under relatively lower elevation laboratory conditions (Mean ± SEM =38.31±2.52) [
44]. This finding was unsurprising as, just like other arthropods, mites’ fitness could be influenced by many abiotic and biotic factors [
44,
47,
48].
Phytoseiid mites, including
N. californicus are widely distributed biological agents for controlling spider mites. They have been used for decades and have attracted substantial attention given their satisfactory capacity to control the harm caused by
T. urticae [
44,
49,
50,
51]. Notably, the prey capacity of Phytoseiid mites on spider mites varies under various conditions, including different temperature ranges [
47,
51], biotic factors such as host plants and prey quality and prey densities [
23]. Here, we elucidated the influence of high altitudes on the prey capacities of four predatory mite species and deduced that compared to the other three predators,
A. swirskii has the highest predation capacity on adult
T. urticae mites at 15 d post-release (Mean ± SEM = 14.28 ± 2.24), indicating an applicability advantage over the other three species in controlling two-spotted mites in the Tibetan Plateau.
Overall, four predatory species, including one species N. californicus, and three species from the genus Amblyseius (A. swirskii, A. barkeri, and A. cucumeris) were investigated herein to compare their fitness and prey capacities on T. urticae in the Tibetan Plateau. We found that A. swirskii had a relatively higher fitness and better predatory capacity than the other predatory mites for controlling T. urticae in a highland area. Furthermore, we discovered that the density of nine predatory mites per plant was the most effective in controlling spider mites and that A. barkeri was an unsuitable biological control agent in highland conditions. Our findings could inform the development of feasible biological control methods for managing T. urticae in the Tibetan Plateau using suitable predatory mite species. They could also be useful in delaying the emergence or increase of resistance to chemical pesticides in two-spotted mites.
Figure 1.
Developmental time of four predatory mite species in Tibetan Plateau. The same letters above columns represent homogeneous groups in post-hoc tests (P > 0.05) following an ANOVA.
Figure 1.
Developmental time of four predatory mite species in Tibetan Plateau. The same letters above columns represent homogeneous groups in post-hoc tests (P > 0.05) following an ANOVA.
Figure 2.
Fecundity of four predatory mite species in Tibetan Plateau. The same letters above columns represent homogeneous groups in post-hoc tests (P > 0.05) following an ANOVA.
Figure 2.
Fecundity of four predatory mite species in Tibetan Plateau. The same letters above columns represent homogeneous groups in post-hoc tests (P > 0.05) following an ANOVA.
Figure 3.
Pre-adult survival rate of four predatory mite species in Tibetan Plateau. The same letters above columns represent homogeneous groups in post-hoc tests (P > 0.05) following by Kruskal–Wallis tests and Dunn’s tests with Bonferroni correction for multiple comparisons.
Figure 3.
Pre-adult survival rate of four predatory mite species in Tibetan Plateau. The same letters above columns represent homogeneous groups in post-hoc tests (P > 0.05) following by Kruskal–Wallis tests and Dunn’s tests with Bonferroni correction for multiple comparisons.
Figure 4.
Longevity of four predatory mite species in Tibetan Plateau. Survival curves for individual hosts were compared using the Kaplan–Meier method and log-rank test.
Figure 4.
Longevity of four predatory mite species in Tibetan Plateau. Survival curves for individual hosts were compared using the Kaplan–Meier method and log-rank test.
Figure 5.
Predatory capacity (preyed spider mites per predator) after release for 1 d of four predatory mite species in Tibetan Plateau. The same letters above columns represent homogeneous groups in post-hoc tests (P > 0.05) following an ANOVA.
Figure 5.
Predatory capacity (preyed spider mites per predator) after release for 1 d of four predatory mite species in Tibetan Plateau. The same letters above columns represent homogeneous groups in post-hoc tests (P > 0.05) following an ANOVA.
Figure 6.
Predatory capacity (preyed spider mites per predator) after release for 15 d of four predatory mite species in Tibetan Plateau. The same letters above columns represent homogeneous groups in post-hoc tests (P > 0.05) following an ANOVA.
Figure 6.
Predatory capacity (preyed spider mites per predator) after release for 15 d of four predatory mite species in Tibetan Plateau. The same letters above columns represent homogeneous groups in post-hoc tests (P > 0.05) following an ANOVA.
Figure 7.
Predatory capacity (preyed spider mites per predator) between release for 1 d and 15 d of four predatory mite species in Tibetan Plateau. Asterisks above points indicate statistically significant differences (Mann–Whitney U-test, *P < 0.05, **P < 0.01, ***P < 0.001, ns, not significant).
Figure 7.
Predatory capacity (preyed spider mites per predator) between release for 1 d and 15 d of four predatory mite species in Tibetan Plateau. Asterisks above points indicate statistically significant differences (Mann–Whitney U-test, *P < 0.05, **P < 0.01, ***P < 0.001, ns, not significant).
Figure 8.
Predatory capacity (A), handling time (B), attack coefficient (C) and max predatory amount (D) after release for 1 d of four predatory mite species in Tibetan Plateau. The max predatory amount data are analyzed by one-way ANOVA, SPSS 21.0, while the data of predatory capacity, handling time and attack coefficient are analyzed by Kruskal–Wallis tests, SPSS 21.0 (*P < 0.05, **P < 0.01, ns, not significant).
Figure 8.
Predatory capacity (A), handling time (B), attack coefficient (C) and max predatory amount (D) after release for 1 d of four predatory mite species in Tibetan Plateau. The max predatory amount data are analyzed by one-way ANOVA, SPSS 21.0, while the data of predatory capacity, handling time and attack coefficient are analyzed by Kruskal–Wallis tests, SPSS 21.0 (*P < 0.05, **P < 0.01, ns, not significant).
Figure 9.
Predatory capacity (A), handling time (B), attack coefficient (C) and max predatory amount (D) after release for 15 d of four predatory mite species in Tibetan Plateau. Results are shown as box and whiskers; the same letters above columns represent homogeneous groups in post-hoc tests (P > 0.05) following an ANOVA, except for attack coefficient, which is analyzed by Kruskal–Wallis tests, SPSS 21.0 (*P < 0.05).
Figure 9.
Predatory capacity (A), handling time (B), attack coefficient (C) and max predatory amount (D) after release for 15 d of four predatory mite species in Tibetan Plateau. Results are shown as box and whiskers; the same letters above columns represent homogeneous groups in post-hoc tests (P > 0.05) following an ANOVA, except for attack coefficient, which is analyzed by Kruskal–Wallis tests, SPSS 21.0 (*P < 0.05).
Table 1.
Parameter estimates of four predatory species’ functional response when feeding on different stages of T. urticae (green) by Holling Ⅱ functional response model (after release for 1d), the 1st repeat.
Table 1.
Parameter estimates of four predatory species’ functional response when feeding on different stages of T. urticae (green) by Holling Ⅱ functional response model (after release for 1d), the 1st repeat.
Predatory mite species |
Tetranychus urticae (green) stage |
Functional response model (Na=) |
R2 |
Chi-sqaure value |
Amblyseius swirskii |
adult |
0.238N0/(1+0.206 N0) |
0.992 |
28.547 |
nymph |
1.346 N0/(1+0.751 N0) |
0.751 |
16.336 |
egg |
1.178 N0/(1+0.506 N0) |
0.767 |
17.043 |
Amblyseius cucumeris |
adult |
0.238 N0/(1+0.329 N0) |
0.783 |
17.794 |
nymph |
0.759 N0/(1+0.289 N0) |
0.925 |
24.793 |
egg |
0.928 N0/(1+0.186 N0) |
0.922 |
24.672 |
Amblyseius barkeri |
adult |
0.246 N0/(1+0.241 N0) |
0.988 |
28.335 |
nymph |
0.925 N0/(1+0.561 N0) |
0.873 |
22.091 |
egg |
0.808 N0/(1+0.343 N0) |
0.879 |
22.382 |
Neoseiulus californicus |
adult |
0.246 N0/(1+0.241 N0) |
0.985 |
28.145 |
nymph |
0.925 N0/(1+0.561 N0) |
0.965 |
27.032 |
egg |
0.808 N0/(1+0.343 N0) |
0.826 |
19.776 |
Table 2.
Parameter estimates of four predatory species’ functional response when feeding on different stages of T. urticae (green) by Holling Ⅱ functional response model (after release for 1d), the 2nd repeat.
Table 2.
Parameter estimates of four predatory species’ functional response when feeding on different stages of T. urticae (green) by Holling Ⅱ functional response model (after release for 1d), the 2nd repeat.
Predatory mite species |
Tetranychus urticae (green) stage |
Functional response model (Na=) |
R2 |
Chi-sqaure value |
Amblyseius swirskii |
adult |
0.204N/(1+0.048N) |
0.977 |
27.677 |
nymph |
0.894N/(1+0.132N) |
0.893 |
23.111 |
egg |
0.605N/(1+0.027N) |
0.965 |
27.033 |
Amblyseius cucumeris |
adult |
0.430N/(1+0.079N) |
0.963 |
26.894 |
nymph |
0.629N/(1+0.063N) |
0.963 |
26.887 |
egg |
0.835N/(1+0.043N) |
0.986 |
28.169 |
Amblyseius barkeri |
adult |
0.235N/(1+0.160N) |
0.978 |
27.745 |
nymph |
0.453N/(1+0.117N) |
0.949 |
26.092 |
egg |
0.973N/(1+0.244N) |
0.875 |
22.189 |
Neoseiulus californicus |
adult |
0.199N/(1+0.026N) |
0.973 |
27.444 |
nymph |
1.020N/(1+0.352N) |
0.760 |
16.744 |
egg |
1.098N/(1+0.146N) |
0.912 |
24.111 |
Table 3.
Parameter estimates of four predatory species’ functional response when feeding on different stages of T. urticae (green) by Holling Ⅱ functional response model (after release for 1d), the 3rd repeat.
Table 3.
Parameter estimates of four predatory species’ functional response when feeding on different stages of T. urticae (green) by Holling Ⅱ functional response model (after release for 1d), the 3rd repeat.
Predatory mite species |
Tetranychus urticae (green) stage |
Functional response model (Na=) |
R2 |
Chi-sqaure value |
Amblyseius swirskii |
adult |
0.431N/(1+0.064N) |
0.961 |
26.774 |
nymph |
0.935N/(1+0.180N) |
0.941 |
25.690 |
egg |
0.885N/(1+0.126N) |
0.946 |
25.925 |
Amblyseius cucumeris |
adult |
0.426N/(1+0.117N) |
0.916 |
24.335 |
nymph |
0.662N/(1+0.118N) |
0.907 |
23.869 |
egg |
1.028N/(1+0.112N) |
0.920 |
24.533 |
Amblyseius barkeri |
adult |
0.529N/(1+0.385N) |
0.886 |
22.573 |
nymph |
0.703N/(1+0.172N) |
0.914 |
24.212 |
egg |
0.882N/(1+0.122N) |
0.940 |
25.616 |
Neoseiulus californicus |
adult |
0.445N/(1+0.196N) |
0.919 |
24.501 |
nymph |
0.657N/(1+0.083N) |
0.934 |
25.321 |
egg |
1.059N/(1+0.133N) |
0.917 |
24.369 |
Table 4.
Parameter estimates of four predatory species’ functional response when feeding on different stages of T. urticae (green) by Holling Ⅱ functional response model (after release for 15d), the 1st repeat.
Table 4.
Parameter estimates of four predatory species’ functional response when feeding on different stages of T. urticae (green) by Holling Ⅱ functional response model (after release for 15d), the 1st repeat.
Predatory mite species |
Tetranychus urticae (green) |
Functional response equation |
R2 |
Chi-sqaure value |
Amblyseius swirskii |
adult |
0.209N/(1+0.075N) |
0.982 |
27.993 |
nymph |
0.631N/(1+0.091N) |
0.939 |
25.573 |
egg |
1.064N/(1+0.097N) |
0.929 |
25.036 |
Amblyseius cucumeris |
adult |
0.214N/(1+0.063N) |
0.995 |
28.684 |
nymph |
0.645N/(1+0.096N) |
0.915 |
24.256 |
egg |
0.867N/(1+0.075N) |
0.979 |
27.810 |
Amblyseius barkeri |
adult |
0.240N/(1+0.340N) |
0.784 |
17.827 |
nymph |
0.466N/(1+0.234N) |
0.899 |
23.455 |
adult |
0.701N/(1+0.225N) |
0.914 |
24.222 |
Neoseiulus californicus |
adult |
0.853N/(1+0.491N) |
0.844 |
20.699 |
nymph |
1.093N/(1+0.131N) |
0.921 |
24.613 |
adult |
1.170N/(1+0.503N) |
0.820 |
19.477 |
Table 5.
Parameter estimates of four predatory species’ functional response when feeding on different stages of T. urticae (green) by Holling Ⅱ functional response model (after release for 15d), the 2nd repeat.
Table 5.
Parameter estimates of four predatory species’ functional response when feeding on different stages of T. urticae (green) by Holling Ⅱ functional response model (after release for 15d), the 2nd repeat.
Predatory mite species |
Tetranychus urticae (green) |
Functional response equation |
R2 |
Chi-sqaure value |
Amblyseius swirskii |
adult |
0.459N/(1+0.018N) |
0.913 |
24.517 |
nymph |
0.919N/(1+0.154N) |
0.934 |
25.315 |
egg |
0.731N/(1+0.029N) |
0.930 |
25.067 |
Amblyseius cucumeris |
adult |
0.196N/(1+0.010N) |
0.968 |
27.169 |
nymph |
1.171N/(1+0.308N) |
0.821 |
19.563 |
egg |
1.099N/(1+0.100N) |
0.977 |
27.654 |
Amblyseius barkeri |
adult |
0.428N/(1+0.092N) |
0.956 |
26.504 |
nymph |
0.663N/(1+0.105N) |
0.930 |
25.101 |
adult |
0.644N/(1+0.062N) |
0.939 |
25.545 |
Neoseiulus californicus |
adult |
0.197N/(1+0.010N) |
0.963 |
26.912 |
nymph |
0.861N/(1+0.076N) |
0.988 |
28.329 |
adult |
1.104N/(1+0.123N) |
0.898 |
23.367 |
Table 6.
Parameter estimates of four predatory species’ functional response when feeding on different stages of T. urticae (green) by Holling Ⅱ functional response model (after release for 15d), the 3rd repeat.
Table 6.
Parameter estimates of four predatory species’ functional response when feeding on different stages of T. urticae (green) by Holling Ⅱ functional response model (after release for 15d), the 3rd repeat.
Predatory mite species |
Tetranychus urticae (green) |
Functional response equation |
R2 |
Chi-sqaure value |
Amblyseius swirskii |
adult |
0.940N/(1+0.195N) |
0.918 |
24.449 |
nymph |
0.903N/(1+0.153N) |
0.917 |
24.385 |
egg |
0.608N/(1+0.030N) |
0.947 |
26.031 |
Amblyseius cucumeris |
adult |
0.455N/(1+0.195N) |
0.910 |
24.001 |
nymph |
1.046N/(1+0.228N) |
0.852 |
21.062 |
egg |
0.833N/(1+0.054N) |
0.950 |
26.173 |
Amblyseius barkeri |
adult |
0.432N/(1+0.126N) |
0.927 |
24.914 |
nymph |
0.657N/(1+0.112N) |
0.959 |
26.676 |
adult |
0.653N/(1+0.076N) |
0.907 |
23.858 |
Neoseiulus californicus |
adult |
0.438N/(1+0.164N) |
0.913 |
24.186 |
nymph |
0.830N/(1+0.078N) |
0.903 |
23.628 |
adult |
0.605N/(1+0.051N) |
0.864 |
21.650 |