Open Access
How to translate text using browser tools
1 March 2013 Arthropod Biodiversity and Community Structures of Organic Rice Ecosystems in Guangdong Province, China
Jie Zhang, Xue Zheng , Hu Jian, Xiaowa Qin, Fenghui Yuan, Runjie Zhang
Author Affiliations +
Abstract

The diversity and community structure of arthropods in an organic double-cropped rice ecosystem in Guangdong Province, China was studied. We compared the arthropod communities in the early season (Apr–Jul) crop to those in the late season (Aug–Nov) crop in 2009. The comparisons were undertaken using a combination of community assessment approaches, including morphospecies richness, the Shannon-Weaver diversity index, H', the Pielou-evenness index, J, the Simpson dominance index C, the Jaccard similarity index q and the compositions of the sub communities. We collected 114 species of arthropods, which consisted of including 58 species of spiders, 16 species of predatory insects, 25 species of phytophagous insects, 15 species of neutral/other insects, in early season crop. Subsequently we collected 109 species of arthropods, which consisted of 50 species of spiders, 19 species of predatory insects, 24 species of phytophagous insects, and 16 species of neutral/other insects, in the late season crop. There were no significant differences (P < 0.05) between the arthropod communities of the early and late season rice crops with respect to the Shannon-Weaver diversity index, the Pielou evenness index and the Simpson dominance index. Moreover the Jaccard similarity index in early and late season rice was as high, i.e., 0.70. The spider sub community had the greatest number of species in both rice crops, but the phytophagous insect sub community had the largest number of individuals in both rice crops. The dominance of predatory insects in the early season rice crop was significantly lower (P < 0.05) than in late season crop, but there was no significant difference in the composition of the neutral/other subcommunity between the early and late season rice crops.

Rice (Oryza sativa L.; Poales: Poaceae) is a major food crop of the world and its cultivation has been carried out in all regions with warm and abundant moisture weather conditions, mainly subtropical regions. Conventional rice cultivation has often accomplished high yields and stable crop production, but has been heavily dependent on continuous and excessive inputs of chemical pesticides, which lead to pest resistance, resurgence, pesticide residue, ground water contamination, and risks to human health and animal habitats (Nagata 1982; Hirai 1993). Organic cultivation of rice has been regarded as a sustainable system because it avoids the problems such as “3Rs” (pest resistance, resurgence, pesticide residue) and other problems of culture heavily dependent on various chemical inputs (Regannold et al. 1990). For example, Kajimura et al. (1993) reported that the population densities of the rice brown planthopper, Nilaparvata lugens (Stål) (Hemiptera: Delphacidae), and the white-backed planthopper, Sogatella furcifera (Horvath) (Hemiptera: Delphacidae), were much lower in an organically farmed than in chemically fertilized rice fields. The maintenance of biodiversity within agricultural environments is widely recognized as being essential for their agronomic sustainability (Swift & Anderson 1994; Maison et al. 1997). An important principle of integrated pest management is to maximize natural control, and, therefore, the temporal changes in arthropod abundance, diversity, species richness and community structures are important considerations in designing pest management strategies. After rice establishment, arthropod species colonize and over time progressively increase in diversity. In rice fields, predators, pollinators and soil microorganisms are all key components of biodiversity (Altieri & Nicholls 1999). Their communities may vary with the environment, varieties, cropping patterns, and cultivation practices. Rice fields often support high levels of biodiversity, which play an important role in the agricultural productivity of these systems (Cohen et al. 1994; Schoenly et al. 1998; O'Malley 1999). Organically grown rice fields reportedly have significantly higher morphospecies richness and diversity than the conventional rice fields (Wilson et al. 2008) and other organic agroecosystems also have a greater richness of arthropod species than their conventional counterparts (Dritschilo & Wanner 1980; Brown & Adler 1989; Goh & Lange 1989; Kromp 1989, 1990). Extensive studies in Asian rice fields have demonstrated that when predator communities are conserved through minimizing pesticide use, the impact of pests such as the N. lugens is often reduced to negligible levels (Way & Heong 1994; Schoenly et al. 1996; Settle et al. 1996).

So far, there has been no relevant study of double-cropped organic rice in Guangdong, China. The objective of this study was to identify arthropod species, describe the structure of arthropod community, examine arthropod abundances of early and late season organic rice, in order to provide theoretical basis for the sustainable control of organic rice pests.

MATERIALS AND METHODS

Study Sites

Three separate rice plots were selected each 30 × 45 m and at least 500m distant from each other, and each plot was bordered on all sides by an unplanted walkway 40 cm wide. These experimental sites were located at Huizhou city (N 23°09′50″ E 114°29′10″), Guangdong province of China, which receives ∼ 1,630 mm of annual rainfall, with 20– 22 °C annual average temperature, red soil and in 2009 the dominant natural vegetation was barnyard grass (Echinochloa sp.). The first early-season rice crop was grown during Apr through Jul followed in the same fields by the late season crop during Aug through Nov.

At all sites the rice cultivar, ‘Haina’ was planted and cultured with organic methods. No synthetic agrochemicals had been used on the land for at least 5 years, and none were applied before or during the production of the organic rice crop at any stage. Agronomic practices such as irrigation for growing rice were the same as followed by local farmers.

Sampling

AT 15 d after transplanting (DAT), the arthropod community was sampled in both the early and late season rice crops at two weeks intervals. Five samples were taken at random in each rice plot. All samples were collected near the center of the plot, at least 5m from the edge in order to reduce edge effects. Arthropods inhabiting the rice field and those on the water surface were sampled using a portable vacuum-suction machine, modified according to Carino et al. (1979). This apparatus collects arthropods through a pyramid-shaped, mylar-covered enclosure (0.5 × 0.5 × 0.9 m high, and 0.25 m2 planar area) fitted with a collecting bag. The enclosure usually covered 4–5 stalks of rice plant after transplanting but only 2–3 after the rice plant had reached maximum tillering. Arthropods inside the enclosure were drawn through a rubber collection hose into a plastic reservoir with a nylon mesh retainer. Sampling duration was fixed at 1 to 2 min depending on the growth stage of the rice. All arthropods (with the exception of parasitoids) inside the enclosure were collected, and then transferred to sample jars containing 70% ethanol, which were returned to the laboratory for identification.

Sampling at each sampling time was repeated 3 times on each of the 3 separate rice plots, and 5 samples were taken each repetition in each area for a total of 15 samples per plot per each sampling time. From Apr 2009 to Nov 2009, the sampling dates were regularly spaced every 2 wk and a total of 180 samples were obtained. All sampled arthropods were examined at low magnification (6.5X) and identified to species.

Data Analysis

The arthropods were separated into 4 functional groups: spiders, predatory insects, phytophagous insects, and neutral and other insects. Neutral insects in this report are a category that consists of insects which do not harm to rice either directly or indirectly.

Alpha species diversity was calculated using the Shannon-Weaver diversity index H', the Pielou-eveness index J, the Simpson dominance index C, and the Jaccard similarity index q (Gianni et al. 2011; Ricardo & Francisco 2011).

Differences between the diversity indices and composition of the arthropod guilds or functional groups in the 2 seasons in the organic rice plots were evaluated with the Tukey test. Species abundance and composition were analyzed using analysis of variance (ANOVA of arcsine, logarithmic and, square root transformed percentages), and means were separated with the Tukey-test as calculated by SPSS 16.0 software. Throughout the text, results are shown as means ±SE.

RESULTS

Abundance of Arthropods

A total of 16,902 individuals were identified as belonging to 135 species of arthropods in 2 classes, 10 orders and 47 families. The morphospecies collected and corresponding taxa of the voucher specimens are shown in Table 1.

The species richness and total number of arthropods in early season rice crop were a little higher than those in the late season crop. There were 114 species of arthropods with 3,177 individuals (the mean number) in the early season crop, and 109 species of arthropods with 2,457 individuals (the mean number) in the late season crop.

Species Richness and Diversity

Some common community indices, specifically the Shannon-Weaver H', the Simpson Dominance C, the Pielou evenness J and the Jaccard Similarity q indices, were calculated for the early and late season rice crops, and are shown in Table 2. There were no significant differences (H': F(2,6) = 2.62, P > 0.05; C: F(2,6) = 4.99, P > 0.05; J: F(2,6) = 0.78, P > 0.05) of these indices between the early and late season crops.

The temporal dynamics of the main indices of arthropod community diversity in the early rice crop and the late rice crop are shown in Fig. 1. The Shannon-Weaver diversity index and the Pielou evenness index of the arthropod community did not differ significantly between the early and late crops at 15, 29, 57 and 85 DAT. However at 43 DAT, the 2 indices of arthropod community in the early crop were significantly lower (P < 0.05) than the community in the late crop, but at 71 DAT, the reverse was true. The Simpson dominance index of the arthropod community in the early crop was higher than in the late crop at 15, 29, 43, 57 and 85 DAT, however, it was significantly lower (P < 0.05) for this community in the late crop than in the early crop at 71 DAT. The Jaccard similarity index q was high between the early and late season rice.

Dominance Distribution

Overall, 114 species of arthropods (58 species of spiders, 16 species of predatory insects, 25 species of phytophagous insects, 15 species of neutral or other insects) and 109 species of arthropods (50 species of spiders, 19 species of predatory insects, 24 species of phytophagous insects, 16 species of neutral or other insects) were observed in the early and late crops, respectively (Table 3).

There was almost significant dominance of the phytophagous insect functional groups in both early and late season rice (Fig. 2). Spiders and predatory insects displayed the second highest level of dominance among the 4 functional groups. The dominance of predatory insects in the early crop was significantly lower than in the late crop (F(2,6) = 1.25, P < 0.05), but there was no significant difference in the other arthropod functional groups between the early and late crops (F(2,6) = 2.47, P > 0.05) (Table 4; Fig. 2).

DISCUSSION

A total of 135 species of arthropods was collected in the present study. The number of arthropod species found in this study was higher than those found by Li et al. (2007) in Hangzhou, Fuyang, but obviously lower than those from tropical areas (Settle et al. 1996; Schoenly et al. 1998). These differences was probably are a function of the different study locations, management regimes and sampling strategies. Our study area at Huizhou was located in the southern part of subtropical China (N 23° 09′50″ E 114° 29′0″), which receives about 2,000 mm of rainfall annually and has an average annual temperature of 22 °C.

With respect to the arthropod community in the present study, we found no significant difference in diversity of between the early and late season rice crops, dominance distribution or evenness of arthropods. Clearly the similarity of the arthropod communities in the early and late rice crops was high. Samples taken from both early and late season organic rice yielded faunal assemblages in our study, showed greater morphospecies richness and higher Shannon diversity indices were found compared with those collected from the conventional rice fields (Tao et al. 1996; Li et al. 2007). This is consistent with the researches in Australian (Wilson et al. 2008) and Jiangmen organic rice fields (Zhong et al. 2005).

TABLE 1.

LIST OF THE ARTHROPOD MORPHOSPECIES AND THE HIGHER TAXONOMIC RANKS OF THE VOUCHER SPECIMENS COLLECTED IN PLOTS OF ORGANICALLY GROWN RICE AT HUIZHOU, GUANGDONG PROVINCE, CHINA.

t01a_01.gif

(CONTINUED)

t01b_01.gif

(CONTINUED)

t01c_01.gif

Diversity of arthropod community varies widely in different seasons and rice growth stages, and the general trend toward arthropod diversity and evenness in early season rice are less than in late season rice (You & Wu 1989). Results in the present study are consistent with this point. The diversity index and evenness index in early season rice were 1.01 and 0.63, respectively, and 1.14 and 0.71, respectively in late season rice.

In the present study, the majority of arthropods were phytophages in both early and late season rice. Taxonomically, the phytophagous insects in the two season rice were dominated by hemipterans. Planthoppers and leafhoppers were the most important components of the phytophagous fauna, and they were high percentages of the phytophages. Some important hemipteran pests such as the N. lugens have been known to cause huge losses to rice production (Heong et al. 1992; Qiu et al. 2004; Backus et al. 2005; Wang & Wang 2007) at 43 DAT, and the dominance index of the arthropod community in early rice was significantly higher than in late rice. This result might have been caused by the high N. lugens population densities in the above mentioned studies, but in our study, the levels of N. lugens remained below the treatment threshold.

Spider populations (including Linyphiidae, Tetragnathidae and Lycosidae) showed positive responses in both early and late season rice. The dominant species of spiders in early season rice were Tetragnatha shikokiana Yaginuma and Hylyphantes graminicola,; whereas the dominant species in the late season rice were Ummeliata insecticeps Boes. et Str. and Pirata subpiraticus Boes. et Str. These spider species were important factors in controlling planthopper and leafhopper populations. Similar responses were also observed in other studies. Such positive response of spider populations have direct impacts on hopper survival and spiders are the key factor in population control of hoppers (Kenmore et al. 1984). Since spiders constitute more than 20% of arthropods in the present study, their impacts would be underestimated by merely comparing numerical relationships.

In this study, predatory insects were mainly composed of the Hemiptera and Coleoptera. Cyrtorrhinus livdipennis Reuter (Hemiptera: Miridae) and Paederus fuscipes Curtis (Coleoptera: Staphylinidae) were the dominant species in both early and late season rice. Although predatory insects represented a lower percentage than spiders in early season rice, they still played an important role in the control of the pests.

TABLE 2.

DIVERSITY INDICES OF THE ARTHROPOD COMMUNITY IN THE EARLY SEASON (APR–JUL) AND LATE SEASON (AUG–NOV) CROPS OF DOUBLE-CROPPED ORGANICALLY GROWN RICE AT HUIZHOU, GUANGDONG PROVINCE, CHINA.

t02_01.gif

Fig. 1.

Temporal dynamics (± SE) of main indices [diversity index (A), dominant index (B) and evenness index (C)] of arthropod community diversity in the early season (Apr–Jul) and late season (Aug–Nov) crops of double-cropped organically grown rice at Huizhou, Guangdong Province, China.

f01_01.jpg

The regulation effects of the neutral insects on pest abundance were mainly realized by natural enemies (Wu et al. 1994). The functional group of natural enemies grew faster than that of pest insects at the earlier rice stages mainly by using neutral insects as prey. The results showed that extremely large populations of Poduridae species and Chironomidae species dominated faunal assemblages in both early and late season rice fields. Although the percentages of neutral and other insects represented a low level in the study, they played an important role in the community food web in paddy fields.

Based on analysis of the arthropod biodiversity and community structure of early and late season organic rice ecosystems, we thus drew the following conclusions: (1) The Shannon-Weaver diversity and Pielou evenness indices in late season rice were a little higher than those of the early season rice; the Simpson dominance index in late season rice was a little lower than that of the early season rice; the Jaccard similarity index between early and late season rice was high up to 0.70; (2) the preponderance of spiders (as having the largest guild membership) was found in both early and late season rice, followed by phytophagous insects, predatory insects, and neutral and other insects; (3) the numerical dominance of phytophagous insect individuals was found in both early and late season rice, followed by numerical dominance of spider individuals.

Fig. 2.

Dominance distribution (± SE) of functional groups of arthropods in the early season (Apr–Jul) and late season (Aug–Nov) crops of double-cropped organically grown rice at Huizhou, Guangdong Province, China.

f02_01.jpg

TABLE 3.

NUMBER OF SPECIES AND INDIVIDUALS FOR EACH SUB-COMMUNITY SAMPLED IN THE EARLY SEASON (APR–JUL) AND LATE SEASON (AUG–NOV) CROPS OF DOUBLE-CROPPED ORGANICALLY GROWN RICE AT HUIZHOU, GUANGDONG PROVINCE, CHINA.

t03_01.gif

TABLE 4.

COMPOSITION OF THE ARTHROPOD SUB-COMMUNITIES IN THE EARLY SEASON (APR–JUL) AND LATE SEASON (AUG–NOV) CROPS OF DOUBLE-CROPPED ORGANICALLY GROWN RICE AT HUIZHOU, GUANGDONG PROVINCE, CHINA.

t04_01.gif

Our results clearly revealed the early and late rice crop arthropod community structures and the dynamics of phytophagous insects, spiders, predatory insects, neutral and other insects in a Guangdong organically grown double-cropped rice ecosystem. These results may provide useful foundation for exploring integrated pest management strategies appropriate for organically grown rice.

ACKNOWLEDGMENTS

We thank the following experts for their assistance in identification of the arthropods: Prof. Gu Dexiang for the identification of the spiders; Prof. Chen Zhenyao for help in identifying the Heteroptera; Prof. Pang Hong and Jia Fenglong for identification of some Coleoptera specimens; Prof. Zhang Dandan for help in identification of the lepidopteran moths. The study was supported by National Science and Technology Support Project (2008BADA5B05).

REFERENCES CITED

1.

M. A. Altieri , and C. I. Nicholls 1999. Biodiversity, ecosystem function, and insect pest management in agricultural systems, pp. 69–84 In W. W. Collins and C. O. Qualset [eds.], Biodiversity in agroecosystems. CRC Press, Boca Raton, FL. Google Scholar

2.

E. A. Backus , M. S. Serrano , and C. M. Ranger 2005. Mechanisms of hopperburn: an overview of insect taxonomy, behavior, and physiology. Annu. Rev. Entomol. 50: 125–151. Google Scholar

3.

M. W. Brown , and C. R. L. Adler 1989. Community structure of phytophagous arthropods on apple. Environ. Entomol. 18: 600–607. Google Scholar

4.

F. O. Carino , P. E. Kenmore , and V. A. Dyck 1979. The farmcop suction sampler for hoppers and predators in flooded rice fields. Int. Rice Res. Newsl. 4: 21–22. Google Scholar

5.

J. E. Cohen , K. Schoenly , K. L. Heong , H. Justo , G. Arida , A. T. Barrion , and J. A. Litsinger 1994. A food web approach to evaluation of the effect of insecticide spraying on insect pest population dynamics in a Philippine irrigated rice ecosystem. J. Appl. Ecol. 31: 747–763. Google Scholar

6.

W. Dritschilo , and D. Wanner 1980. Ground beetle abundance in organic and conventional corn fields. Environ. Entomol. 9: 629–631. Google Scholar

7.

Q. H. Gianni , J. C. Francisco , and C. F. M. Ivan 2011. Species diversity of myrmecofauna and araneofauna associated with agroecosystem and forest fragments and their interaction with carabidae and staphylinidae. Florida Entomol. 94: 500–509. Google Scholar

8.

K. S. Goh , and W. H. Lange 1989. Arthropods associated with insecticide-treated and untreated artichoke fields in California. J. Econ. Entomol. 82: 621–625. Google Scholar

9.

K. Hirai 1993. Recent trends of insecticide susceptibility in the brown planthopper Nilaparvata lugens (Stål) (Homoptera: Delphacidae). Appl. Entomol. Zool. 28: 339–346. Google Scholar

10.

K. L. Heong , G. B. Aquino , and A. T. Barrion 1992. Population dynamics of plant- and leafhoppers and their natural enemies in rice ecosystems in the Philippines. Crop Prot. 11: 371–379. Google Scholar

11.

T. Kajimura , Y. Maeoka , I.N. Widiarta , T. Sudo , K. Hidaka , F. Nakasuji , and K. Nagai 1993. Effect of organic farming of rice plants on population density of leafhoppers and planthoppers. I. Population density and reproductive rate. Japanese J. Appl. Entomol. Zool. 37: 137–144. Google Scholar

12.

P. E. Kenmore , F. O. Carino , G. A. Perez , V. A. Dyck , and A. P. Gutierrez 1984. Population regulation of the rice brown planthopper (Nilaparvata lugens Stål) with rice fields in the Philippines. J. Plant Prot. Trop. 1: 19–38. Google Scholar

13.

B. Kromp 1989. Carabid beetle communities (Carabidae, Coleoptera) in biologically and conventionally farmed agroeco-systems. Agri. Eco. Environ. 27: 241–251. Google Scholar

14.

B. Kromp 1990. Carabid beetles (Carabidae, Coleoptera) as bioindicators in biological and conventional farming in Austrian potato fields. Biol. Fert. Soils. 9: 182–187. Google Scholar

15.

F. F. Li , G. Y. Ye , Q. Wu , Y. F. Peng , and X. X. Chen 2007. Arthropod abundance and diversity in Bt and non-Bi rice fields. Environ. Entomol. 36: 646–654. Google Scholar

16.

P. A. Matson , W. J. Parton , A. G. Power , and M. J. Swift 1997. Agricultural intensification and ecosystem properties. Science. 277: 504–509. Google Scholar

17.

T. Nagata 1982. Insecticide resistance and chemical control of the rice planthopper Nilaparvata lugens Stål. Bull. Kyushu Natl. Agri. Experi. Stat. 22: 49–164. Google Scholar

18.

R. E. O'Malley 1999. Agricultural wetland management for conservation goals, pp. 857–885 In D. Batzer , R. Rader , S. Wissinger [eds.], Invertebrates in freshwater wetlands of North America. Wiley, New York. Google Scholar

19.

J. P. Ricardo , and J. P. F. Francisco 2011. Diversity and community structure of opiinae (Hymenoptera: Braconidae) in the forest estate of Artikutza (Spain). Florida Entomol. 94: 472–479. Google Scholar

20.

H. M. Qim , J. C. Wu , G. Q. Yang , B. Dong , and D. H. Li 2004. Changes in the uptake function of the rice root to nitrogen, phosphorus and potassium under brown planthopper, Nilaparvata lugens (Stål) (Homoptera: Delphacidae) and pesticide stresses, and effect of pesticides on rice-grain filling in field. Crop Prot. 23: 1041–1048. Google Scholar

21.

J. P. Regannold , R. I. Papendick , and J. F. Parr 1990. Sustainable agriculture. Sci. Amer. 262: 112–120. Google Scholar

22.

M. J. Swift , and J. M. Anderson 1994. Biodiversity and ecosystem function in agricultural systems, pp. 15–41 In E. D. Schulze and H. A. Mooney [eds.], Biodiversity and ecosystem function. Springer, Berlin. Google Scholar

23.

K. Schoenly , J. E. Cohen , K. L. Heong , G. Arida , A. T. Barrion , and J. A. Litsinger 1996. Quantifying the impact of insecticides on food web structure of rice-arthropod populations in Philippine farmers' irrigated fields: a case study, pp. 343–351 In G. Polis and K. Winemiller [eds.], Food Webs: Integration of Patterns and Dynamics. Chapman and Hall, New York. Google Scholar

24.

K. G. Schoenly , H. D. Justo , A. T. Barrion , M. K. Harris , and D. G. Bottrell 1998. Analysis of invertebrate biodiversity in a Philippine farmer's irrigated rice field. Environ. Entomol. 27: 1125–1136. Google Scholar

25.

W. H. Settle , H. Ariawan , E. T. Astuti , W. Cahyna , A. L. Hakim , D. Hindayana , and A. Srilestari 1996. Managing tropical rice pests through conservation of generalist natural enemies and alternative prey. Ecol. 77: 1975–1988. Google Scholar

26.

F. L. Tao , G. W. Liang , and X. F. Pang 1996. Studies on the seasonal dynamics of arthropod communities in rice field indifferent habitats. J. South China Agric. Univ. 17: 25–30. Google Scholar

27.

Y. H. Wang , and M. H. Wang 2007. Factors affecting the outbreak and management tactics of brown planthopper in China recent years. Pestic. Sci. Admin. 25: 49–54. Google Scholar

28.

M. J. Way , and K. L. Heong 1994. The role of biodiversity in the dynamics and management of insect pests of tropical irrigated rice-a review. Bull. Entomol. Res. 84: 567–587. Google Scholar

29.

A. L. Wilson , R. J. Watts , and M. M. Stevens 2008. Effects of different management regimes on aquatic macroinvertebrate diversity in Australian rice fields. Ecol. Res. 23: 565–572. Google Scholar

30.

J. C. Wu , G. W. Hu , J. Tang , Z. L. Shu , J. S. Yang , Z. N. Wan , and Z. C. Ren 1994. Studies on the regulation effect of neutral insect on the community food web in paddy field. Acta Ecol. Sinica. 14: 381–386. Google Scholar

31.

M. S. You , and Z. F. Wu 1989. The diversity of the arthropod communities in paddy fields. J. Fujian Agric. Forest Univ. 18: 532–538. Google Scholar

32.

P. S. Zhong , G. W. Liang , and L. Zeng 2005. Biodiversity of Major Natural Enemies in Organic Farming Rice Fields. Chinese J. Biol. Contr. 21: 155–158. Google Scholar
Jie Zhang, Xue Zheng , Hu Jian, Xiaowa Qin, Fenghui Yuan, and Runjie Zhang "Arthropod Biodiversity and Community Structures of Organic Rice Ecosystems in Guangdong Province, China," Florida Entomologist 96(1), 1-9, (1 March 2013). https://doi.org/10.1653/024.096.0101
Published: 1 March 2013
KEYWORDS
abundance
abundancia
arthropod biodiversity
biodiversidad de artrópodos
community dynamics
dinámica de la comunidad
ecosistemas de arroz orgánico
Back to Top