Skip to main content

Mapping the potential distribution of high artemisinin-yielding Artemisia annua L. (Qinghao) in China with a geographic information system



Artemisia annua L. is an important source for artemisinin, a potent drug for treating malaria. This study aims to map and predict the potential geographic distribution of A. annua L. in China.


The Geographic Information System for traditional Chinese medicine (TCM-GIS) was developed and used to map the potential geographic distribution of A. annua L.


Climatic, edaphic and topographic characteristics of A. annua L. microhabitats in Youyang County were mapped to find distribution patterns. The maps identified that certain habitats in the Chongqing region and some potential regions, especially in Guizhou Province, possess similarity indices of ≥98%. In particular, high quality microhabitats A. annua L. were found in the Wuling mountains region.


The present study demonstrates a GIS approach to predict potential habitats for A. annua L. TCM-GIS is a powerful tool for assessing bioclimatic suitability for medicinal plants.


Artemisia annua L. (Qinghao, Annual Wormwood) is a strongly fragrant, annual herbaceous plant used in Chinese medicine [1]. A. annua L. is the only natural botanical source for artemisinin (Qinghaosu) [2, 3] and a potential source for essential oils for the perfume industry [4]. A. annua L. is now cultivated in China, Vietnam, India, Romania, Kenya and Tanzania [5]. Artemisinin, an endoperoxide sesquiterpene lactone in the aerial parts of A. annua L., is more efficacious, faster and less toxic than chloroquine in treating malaria. In addition, artemisinin is a potent anti-cancer agent, a possible antibacterial agent as well as a natural pesticide [6, 7]. Chemical and biological synthesis of artemisinin is still under development due to poor yields [811]. Therefore, wild or cultivated A. annua L. is a major source for artemisinin [2, 3, 12].

The artemisinin content is highly dependent on plant ecotypes, ecological interactions, seasonal and geographic variations [1318]. In fact, artemisinin is absent in some A. annua L. Artemisinin was first isolated in China and some Chinese germplasm has relatively higher artemisinin levels than those of Europe, North America, East Africa and Australia [2, 13, 16, 17, 19, 20]. In Youyang County, Chongqing, China, the hometown of A. annua L., the plants have high (0.9%) levels of artemisinin. In 2006 the county became a national protected geographic area recognized by the General Administration of Quality Supervision, Inspection and Quarantine of China [21]. As the demand for artemisinin remains high around the world, finding suitable geographic regions for A. annua L. is a critical research area for the World Health Organization [22].

The geographic information system (GIS) technology manages geographic information with applications for various fields such as natural resources, transportation planning, environmental studies and vegetation distribution studies [2326]. Recently updated, the geographic information system for traditional Chinese medicine (TCM-GIS) captures, stores, analyzes and displays geographically referenced information to analyze genetic, ecological and geographic patterns of the spatial distribution of a target species. Using the TCM-GIS, our previous studies analyzed the potential habitats and distributions of Chinese medicinal plants such as Glycyrrhiza uralensis Fisch., Panax quinquefolium and Panax ginseng[2729]. The present study aims to characterize the eco-environmental conditions in the A. annua L. production areas in Youyang County and predict the potential distributions of A. annua L. with a high artemisinin-yielding potentials.


Data collection

The spatial distribution of A. annua L. was based on the following four sources: (1) the flora of China [30]; (2) scientific literature concerning the geographic distribution of A. annua L. in China [31]; (3) the Chinese Virtual Herbarium (CVH) [32], (4) germplasm accessions from the Sharing Information System for Chinese Medicinal Plant Germplasm Resources [33]; (5) field data of wild A. annua L. and interviews in Youyang County in 2008. Due to the excellent quality of A. annua L. from the habitats in Youyang County [31, 3436], a total of 180 accessions of A. annua L. germplasm were collected and used in the present study.

The potential distribution mapping program TCM-GIS and geo-referenced datasets were used to develop eco-adaptation models. The TCM-GIS package included three databases, namely (1) a basic geographic information database including digital line graphics and a digital elevation model (scale: 1:1,000,000), (2) a soil database (scale: 1:4,000,000), (3) and a climate database (mean values between 1971 and 2000). All three databases were used for spatial analysis and model calibration.

Raster and vector are two main data models in the TCM-GIS. Raster layers (1 × 1 km2 resolution) were used for the eco-environmental analysis and cluster analysis. Vector layers were used to derive and identify the spatial extent and location of suitable habitats through overlay analysis. Moreover, global positioning system data on the locations of the 180 accessions were obtained for villages such as Banqiao, Zhongduo, Mawang and Nanmu and used in the TCM-GIS analysis (Figure 1).

Figure 1

Spatial distribution of A. annua L. germplasm collection sites as mapped with the TCM-GIS.

In the present study, 14 eco-environmental variables were chosen for the predication of spatial distribution in Youyang County. These variables, namely (1) average temperature in January (ATJA), (2) average temperature in February (ATF), (3) average temperature in March (ATM), (4) average temperature in April (ATAP), (5) average temperature in May (ATMA), (6) average temperature in June (ATJ), (7) average temperature in July (ATJU), (8) average temperature in August (ATA), (9) average annual temperature (AAT), (10) annual sunshine time (AST), (11) total annual precipitation (TAP), (12) relative humidity (RH), (13) altitude (AL), (14) and soil properties (SP), were classified into three categories: topography, climate and edaphology (Table 1).

Table 1 Environmental factors used in this study.

Data analysis

An optimal range was established by identifying minima and maxima for eco-environmental variables (e.g. elevation and temperature) at sample collection sites. The A. annua L. macro-habitats were characterized by examining the mean, minimal and maximal values, standard deviation (SD), standard error (SE), and coefficient of variation (CV) of these variables (Table 2). Prior to distance analysis, we normalized the raster grid data representing each variable. We derived the mean absolute deviation using the following equation:

Table 2 Summary of eco-environmental characteristics from known A. annua L. habitats (n = 180).

where xkf was the measured values of the variable f and mf is the mean for the variable f. For the determination of similarity between grid data and eco-factor ranges, the statistical distance was calculated with the Minkowski distance equation [37]:

which is a generalization of the Euclidean distance and Manhattan distance; in general the shorter the distance, the greater the similarity. The comprehensive similarity index (SI) of each factor layer was calculated with an overlay analysis with various weighting values. Finally, maps with two ranks of predictive distributions were generated, followed by a grid-based spatial cluster analysis, vector-based overlaying, intersection analysis and an area calculation (Figures 2, 3, 4, Table 3).

Figure 2

Distribution of suitable A. annua L. production areas in China with a similarity index (SI) of 99-100%

Figure 3

Distribution of suitable A. annua L. production areas in China with a similarity index (SI) of 98-99%.

Figure 4

Suitable regions for A. annua L. production with a similarity index (SI) of ≥98%.

Table 3 Major A. annua L. regions with similarity index (SI) of 99%-100%.

The most favorable region for A. annua L. growth is one that has an SI range of 99%-100%, while the second-most favorable region is one that has an SI range of 98%-99%.

Results and Discussion

Eco-environmental preferences

The climatic, edaphic and topographic characteristics of known A. annua L. habitats are listed in Table 2. While low CV values for RH (CV: 0.33), TAP (1.28), AST (3.33), ATJU (4.60), AAT (4.69), ATA (6.23), ATJ (6.77) and ATMA (6.81) suggested that these could be the major limiting factors affecting the distribution of high quality A. annua L., high CV values for AL (29.79), ATJA (21.46) and ATF (21.43) suggested otherwise. According to the CV values, weighting value for each parameter was divided into levels I (0.15), II (0.08), III (0.06) and IV (0.03) and weighting values should add up to one. In addition, datasets of eco-factors from known habitats in Youyang County were as follows: ATJA = 1.2-5.6°C, ATF = 2.0-6.0°C, ATM = 4.0-10.0°C, ATAP = 10.0-16.0°C, ATMA = 14.0-20.0°C, ATJ = 18.0-24.0°C, ATJU = 21.6-27.3°C, ATA = 20.0-26.0°C, AAT = 15.9-21.0°C, AST = 1048-1200 h, TAP = 1169-1267 mm, RH = 79.2-80.6%, AL = 498-1010 mm. Soil types were mainly yellow soil, yellow sandy soil, limestone soil, paddy soil and brown soil with pH value at 6-7 and organic matter content ≥1.3%. Thus, we assumed that these conditions were optimal for the growth of high artemisinin-yielding A. annua L.

A. annua L. is a short-day plant. Non-juvenile plants are very responsive to short photoperiodic stimuli and flower about two weeks after induction. They require about 1000 hours of sunlight per year. Our results suggest that annual sunlight time is a critical factor for the growth of A. annua L., which is consistent with previous studies [5, 38]. Previous findings that A. annua L. requires a strict watering regime during the preliminary growth stages [5, 39] are also consistent with our results.

Predictive maps

Figures 2 and 3 are the maps derived from the TCM-GIS analyses. The predicted areas were primarily located in the Wuling Mountain region in central China, covering Guizhou, Chongqing, Hunan, Hubei and Sichuan (25°14'-31°38' N to 104°31'-111°51'E). The predicted habitat density was high in northeastern Guizhou, southeastern Chongqing, northwestern Hunan, southwestern Hubei and parts of southern Sichuan.

The total favorable regions (SI 98%-99%) made up 1.60% of China's total land area covering 162 counties and cities (a total of 60,292 km2), among which Guizhou took the lead with 31,150 km2 including 68 counties and cities. The most favorable region for A. annua L. (SI 99%-100%) was in the 58 counties and cities in Guizhou Province with a predicted area of 54,350 km2. The second largest predicted area (14,330 km2) was in the 12 counties and cities in Chongqing, followed by Hunan, Hubei and Sichuan (Figure 4). The counties and cities with significant areas of potential habitat are listed in Table 3. The data indicated that Youyang County contained the largest favorable area with more than 4000 km2. Unexpectedly, the total predicted areas in Wuchuan and Zunyi Counties in Guizhou exceeded 2000 km2.

One of the world's largest artemisinin manufacturers and its affiliates operate A. annua L. farms in the Chongqing Wulingshan Mountain Range [40, 22]. Apart from this, Guizhou may be another important region for A. annua L. cultivation, particularly in the northeastern part of the province. Our model predicted that 13% of this area is potential A. annua L. habitats [41, 42]. Our model did not predict Guangxi Province, known for its habitats of A. annua L. of relatively low quality, as a region for A. annua L. cultivation possibly due to the subtropical climate, low altitude and red soil in Guangxi which are very different from those in other A. annua L. regions in China [9].

Interviews with the locals suggest that the Guizhou region and Youyang County have comparative advantages for A. annua L. growth with a high-yield variety and minimal pests. Furthermore, the northeastern Guizhou is home to wild populations of A. annua L. which may be an alternative source for artemisinin.

Using the TCM-GIS, we aimed to determine the optimal ecological factors from known habitats and the results showed that RH, TAP, AST, STJU, AAT and SP were important limiting factors. We also aimed to map the distribution of potential regions for the development of A. annua L. in China based on selected climatic, soil and topographical values. Using bioclimatic similarity theory and the TCM-GIS, we predicted the potential growing areas at the county level, particularly in northeastern Guizhou Province. The TCM-GIS is adequate for predicting and identifying potential areas for A. annua L. cultivation.

Using a higher resolution raster and vector spatial databases, we improved the resolution of species distribution considerably on the national surveys conducted in the 1960s, 1970s and 1980s. While most of the survey data were based largely on personal experiences and rough estimates, the model used in the present study is relatively objective.


The present study demonstrates a GIS approach to predict the potential habitats for A. annua L. TCM-GIS is a powerful tool for assessing bioclimatic suitability for medicinal plants.



traditional Chinese medicine geographic information system


geographic information system


similarity index


standard deviation


standard error


coefficient of variation


average temperature in January


average temperature in February


average temperature in March


average temperature in April


average temperature in May


average temperature in June


average temperature in July


average temperature in August


average annual temperature


annual sunshine time


total annual precipitation


relative humidity




soil properties


Chinese Virtual Herbarium.


  1. 1.

    Wang CW: The forests of China: with a survey of grassland and desert vegetation. Maria Moors Cabot Foundation Publication. 1961, Cambridge Massachusetts: Harvard University, 5: 155-164.

    Google Scholar 

  2. 2.

    Jain DC, Mathur AK, Gupta MM, Singh AK, Verma RK, Gupta AP, Kumar S: Isolation of high artemisinin-yielding clones of Artemisia annua. Phytochemistry. 1996, 5: 993-1001. 10.1016/S0031-9422(96)00369-X.

    Article  Google Scholar 

  3. 3.

    Klayman DL: Qinghaosu (artemisinin): an antimalarial drug from China. Science. 1985, 228: 1049-1055. 10.1126/science.3887571.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Simon JE, Charles E, Cebert L, Grant J, Janick J, Whipkey A: Artemisia annua L: a promising aromatic and medicinal. Advances in New Crops: Proceeding of the First National Symposium New Crops: Research, Development, Economics. 1990, Portland: Timber Press, Incorporated, 522-526.

    Google Scholar 

  5. 5.

    World Health Organization: WHO monograph on good agricultural and collection practices (GACP) for Artemisia annua L. Geneva. 2006

    Google Scholar 

  6. 6.

    Efferth T, Dunstan H, Sauerbrey A, Miyachi H, Chitambar CR: The anti-malarial artesunate is also active against cancer. Int J Oncol. 2001, 18: 767-773.

    CAS  PubMed  Google Scholar 

  7. 7.

    Kumar S, Khanuja SPS, Kumar TRS, Jain DC, Srivastava S, Bhattacharya AK, Saikia D, Shasany AK, Darokar MP, Sharma RP: Method for the use of alpha arteether as an anti-bacterial and anti-fungal agent. US Patent 6127405. 2000

    Google Scholar 

  8. 8.

    Avery MA, Chong WKM, Jennings-White C: The total synthesis of (+)-Artemisinin and (+)-9-Desmethylartemisinin. Tetrahedron Lett. 1987, 28: 4629-4632. 10.1016/S0040-4039(00)96582-1.

    CAS  Article  Google Scholar 

  9. 9.

    Ravindranathan T, Kumar MA, Menon RB, Hiremath SV: Stereoselective synthesis of artemisinin. Tetrahedron Lett. 1990, 31: 755-758. 10.1016/S0040-4039(00)94621-5.

    CAS  Article  Google Scholar 

  10. 10.

    Van Geldre E, Vergauwe A, Eeckhout Van den E: State of art of the production of the antimalarial compound artemisinin in plants. Plant Mol Biol. 1997, 33: 199-209. 10.1023/A:1005716600612.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Chen DH, Li GF, Ye HC: Expression of a chimeric farnesyl diphosphate synthase gene in Artemisia annua L. transgenic plants via Agrobacterium tumefaciens-mediated transformation. Plant Sci. 2000, 155: 179-185. 10.1016/S0168-9452(00)00217-X.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Kumar S, Gupta SK, Singh P, Bajpai P, Gupta MM, Singh D, Gupta AK, Ram G, Shasany AK, Sharma S: High yields of artemisinin by multi-harvest of Artemisia annua crops. Ind Crops Prod. 2004, 1: 77-90. 10.1016/j.indcrop.2003.07.003.

    Article  Google Scholar 

  13. 13.

    Delabays N, Benakis A, Collet G: Selection and breeding for high artemisinin (Qinghaosu) yielding strains of Artemisia annua. Acta Hort (ISHS). 1993, 330: 203-208.

    Article  Google Scholar 

  14. 14.

    Liersch R, Soicke H, Stehr C, Tullner HU: Formation of artemisinin in Artemisia annua during one vegetation period. Planta Med. 1986, 52: 387-390. 10.1055/s-2007-969193.

    Article  Google Scholar 

  15. 15.

    Singh A, Vishwakarma RA, Husain A: Evaluation of Artemisia annua strains for higher artemisinin production. Planta Med. 1988, 64: 475-476. 10.1055/s-2006-962515.

    Article  Google Scholar 

  16. 16.

    Charles DJ, Simon JE, Wood KV, Heinstein P: Germplasm variation in artemisinin content of Artemisia annua using an alternative method of artemisinin analysis from crude plant extracts. J Nat Prod. 1990, 53: 157-159. 10.1021/np50067a021.

    CAS  Article  Google Scholar 

  17. 17.

    Woerdenbag HJ, Pras N, Chan NG, Bang BT, Bos R, van Uden W, Van YP, van Boi NV, Batterman S, Lugt CB: Artemisinin related sesquiterpenes and essential oil in Artemisia annua during a vegetation period in Vietnam. Planta Med. 1994, 60: 272-275. 10.1055/s-2006-959474.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Gupta SK, Singh P, Bajpai P, Ram G, Singh D, Gupta MM, Jain DC, Khanuja SPS, Kumar S: Morphogenetic variation for artemisinin and volatile oil in Artemisia annua. Ind Crops Prod. 2002, 16: 217-224. 10.1016/S0926-6690(02)00049-3.

    CAS  Article  Google Scholar 

  19. 19.

    Trigg EI: Qinghaosu (artemisinin) as an antimalarial drug. Econ Med Plant Res. 1989, 3: 19-55.

    Article  Google Scholar 

  20. 20.

    Klayman DL, Lin AJ, Acton N, Scovill JP, Hoch JM, Milhous WK, Theoharides AD, Dobek AS: Isolation of artemisinin (qinghaosu) from Artemisia annua growing in the United States. J Nat Prod. 1984, 47: 715-717. 10.1021/np50034a027.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Genearal Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China. []

  22. 22.

    World Health Organization Global Malaria Programme: Proceedings of the Meeting on the production of artemisinin and artemisinin-based combination therapies: 6-7 June 2005; Arusha. 2006, Global Malaria Programme

    Google Scholar 

  23. 23.

    Goodchild MF: Geographic information systems and science: today and tomorrow. Procedia: Earth and Planetary Science. 2009, 1: 1037-1043. 10.1016/j.proeps.2009.09.160.

    Google Scholar 

  24. 24.

    Segura S, Coppens d'Eeckenbrugge G, Lòpez L, Grum M, Guarino L: Mapping the potential distribution of five species of Passiflora in Andean countries. Genet Resour Crop Evol. 2003, 12: 555-566. 10.1023/A:1024424012134.

    Article  Google Scholar 

  25. 25.

    Chen SL, Suo FM, Han JP, Xie CX, Yao H, Li XW, Li Y, Wei JH: Analysis on ecological suitability and regionalization of traditional Chinese medicinal materials. Chin Tradit Herb Drugs. 2007, 38: 481-487.

    CAS  Google Scholar 

  26. 26.

    Smith-Ramirez C, Diaz I, Pliscoff P, Valdovinos C, Mendez MA, Larrain J, Samaniego H: Distribution patterns of flora and fauna in southern Chilean coastal rain forests: Integrating natural history and GIS. Biodivers and Conserv. 2007, 16: 2627-2648. 10.1007/s10531-006-9073-2.

    Article  Google Scholar 

  27. 27.

    Wang JY, Zhao RH, Sun CZ, Chen SL Liu SQ, Wei JH, Xiao XX, Zhou YQ: Suitability evaluation of Glycyrrhiza Uralensis Fisch.'s distributive area based on TCMGIS. Mod Chin Med. 2006, 8: 4-8.

    Google Scholar 

  28. 28.

    Chen SL, Zhou YQ, Xie CX, Zhao RH, Sun CZ, Wei JH, Liu SQ, Gao WW: Suitability evaluation of Panax quinquefolium producing area based on TCMGIS. Zhongguo Zhong Yao Za Zhi. 2008, 33: 741-745.

    PubMed  Google Scholar 

  29. 29.

    Wang Y, Wei JH, Chen SL, Sun CZ, Zhao RH, Liu SQ, Xiao XX, Wang JY, Zhou YQ: Analysis of Panax Ginseng' producing area based on TCMGIS. Asia-Pacific Trad Med. 2006, 6: 73-78.

    Google Scholar 

  30. 30.

    Academiae Sinicae Edita: Flora Reipublicae Popularis Sinicae. 1989, Beijing: Science Press, 76 (2): 62-65.

    Google Scholar 

  31. 31.

    Zhong GY, Zhou HR, Ling Y, Hu M, Zhao PP: Investigation on ecological environment and quantitative analysis of Artemisinin of Sweet Wormwood (Artemisia annua). Chin Med Mat. 1998, 29: 264-267.

    Google Scholar 

  32. 32.

    The Chinese Virtual Herbarium (CVH). []

  33. 33.

    Germplasm Accessions From the Sharing Information System for Chinese Medicinal Plant Germplasm Resources. []

  34. 34.

    Yang SP, Yang X, Huang JG, Ding DR: Advances in researches on Artemisinin production. J Trop Subtrop Bot. 2004, 12: 189-194.

    CAS  Google Scholar 

  35. 35.

    Xu XY, Zheng YM, Fu SQ, Han YM, Yang YH: Studies on HPLC Fingerprint of A rtemisia annua L. from Chongqing. Lishizhen Med Mater Med Res. 2009, 20: 1188-1189.

    CAS  Google Scholar 

  36. 36.

    Xia ZZ: Studies on the methods for assay of artemisinin in Artemisia annual L. and its semi-synthetic derivatives. Master thesis. 2006, Chongqing Medical University, School of Pharmaceutical Sinence

    Google Scholar 

  37. 37.

    He YQ: Multivariate Statistical Analysis. 2004, Beijing: China Renmin University Press

    Google Scholar 

  38. 38.

    Purdue University Annual Wormwood (Artemisia annua L.). []

  39. 39.

    Wei X, Li F, Xu CQ, Fu XH, Huang ZF: Study on biological characteristics of Artemisia annua L. Guihaia. 1997, 17: 166-168.

    Google Scholar 

  40. 40.

    Holley Pharmaceuticals Company. []

  41. 41.

    Zhou AA, Zhen WX, Ge FH: Determination of artemisinin in Artemisia annuna by HPLC-ELSD. Zhong Yao Cai. 2006, 29: 242-245.

    CAS  PubMed  Google Scholar 

  42. 42.

    Wang YL, Zhang L, Zhang HM: Determination of artemisinin in Artemisia annua from different habitats in Changde region Hunan province. Pharmaceut care Res. 2007, 5: 381-382.

    Google Scholar 

Download references


The authors would like to thank the National Natural Science Foundation of China for its support through project No. 3050081.

Author information



Corresponding author

Correspondence to Shilin Chen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

LFH, SLC and CXX designed the study and performed the analyses. BZD helped with data analysis. All authors wrote the manuscript. All authors read and approved the final version of the manuscript.

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Huang, L., Xie, C., Duan, B. et al. Mapping the potential distribution of high artemisinin-yielding Artemisia annua L. (Qinghao) in China with a geographic information system. Chin Med 5, 18 (2010).

Download citation


  • Total Annual Precipitation
  • Potential Habitat
  • Artemisia Annua
  • Artemisinin Content
  • Chinese Medicinal Plant