TCMGIS-II based prediction of medicinal plant distribution for conservation planning: a case study of Rheum tanguticum
© Yu et al; licensee BioMed Central Ltd. 2010
Received: 6 November 2009
Accepted: 25 August 2010
Published: 25 August 2010
Many medicinal plants are increasingly endangered due to overexploitation and habitat destruction. To provide reliable references for conservation planning and regional management, this study focuses on large-scale distribution prediction of Rheum tanguticum Maxim. ex Balf (Dahuang).
Native habitats were determined by specimen examination. An improved version of GIS-based program for the distribution prediction of traditional Chinese medicine (TCMGIS-II) was employed to integrate national geographic, climate and soil type databases of China. Grid-based distance analysis of climate factors was based on the Mikowski distance and the analysis of soil types was based on grade division. The database of resource survey was employed to assess the reliability of prediction result.
A total of 660 counties of 17 provinces in China, covering a land area of 3.63 × 106 km2, shared similar ecological factors with those of native habitats appropriate for R. tanguticum growth.
TCMGIS-II modeling found the potential habitats of target medicinal plants for their conservation planning. This technology is useful in conservation planning and regional management of medicinal plant resources.
More than one-tenth of plant species are used in drugs and health products . The demand for herbal drugs and health products is steadily growing . Thus, many medicinal herbs are threatened by overexploitation, habitat destruction and lack of proper cultivation practices. Some wild species are disappearing at alarming rates [3, 4]. Rheum tanguticum Maxim. ex Balf (Dahuang) is one of those species. R. tanguticum belongs to the family Polygonaceae and is a high-altitude perennial herb sensitive to high temperature, mainly found in the alpine regions of temperate and subtropical Asia, especially in Southwest and Northwest China (e.g. Sichuan, Gansu and Qinghai) [5, 6]. As a source for rhubarb according to the Chinese Pharmacopoeia and a purgative and anti-inflammatory agent , R. tanguticum has been overexploited, suffering from replant diseases, inadequate seed dispersal, low reproductive efficiency and narrow distribution and habitat fragmentation, leading to its declines in the wild resources [6, 8].
In-situ conservation, which considered as the method of conserving endangered species in their wild habitats, is promising in protecting indigenous species and maintaining natural communities along with their intricate network of relationships . As habitat degradation and destruction is increasing, ex-situ conservation regarded as the process of cultivating and naturalizing endangered species outside of their original habitats, has become a practical alternative [10–12], especially for those overexploited and endangered medicinal plants with slow growth, small abundance and replant diseases [10, 13], e.g. Paris species in family Trilliaceae and Panax species in family Araliaceae . Ex-situ cultivation becomes an immediate action to sustain medicinal plant resources [11, 12].
Understanding the geographical distribution of plant species is essential for their ex-situ conservation activities [1, 15]. Although many plant species can be successfully introduced, cultivated and naturalized in a wide range of habitats across countries and continents , their growth and distribution in different habitats are based on local indicators , e.g. soil properties, climate conditions and environmental features . Aguilar-Stoen and Moe (2007) found that many medicinal plants thriving in harsh habitats and disturbed areas are of high medicinal efficacy because rocky and dry habitats stimulate their secondary metabolites . Many plants are only found in places where the habitat is congruent with their growth , e.g. the propagation and quality of Banksia serrata varied among habitats . Variations in growth and metabolites of medicinal plants among niches make ex-situ conservation habitat-specific.
Geographical prediction of plant distribution is important to resource conservation planning and regional management decisions . Geographic Information System (GIS) is useful in predicting the spatial distribution of target species . GIS assesses multiple interdependent abiotic factors, e.g. solar radiation, air temperature, precipitation and soil properties , affecting plant distribution, models the environmental niches of target plants  and refines their distribution maps for conservation planning .
A GIS-based computer program (TCMGIS-I) was developed specially for the distribution prediction of Chinese medicine (CM) [25, 26]. Integrating national geographic, climate and soil type databases of China, TCMGIS-I was able to determine the impacts of environmental gradients and predict the large-scale distribution of target medicinal plants . Tests with some common medicinal plants (e.g. Panax ginseng, Panax quinquefolium, Glycyrrhiza uralensis and Artemisia annua) demonstrated that TCMGIS-I prediction was consistent with the actual plants' distribution patterns [27–30].
While TCMGIS-I captures data from literature, TCMGIS-II can perform more precise variable extraction from the native habitats of target medicinal plants. Factors such as elevation, air temperature, solar radiation, precipitation and soil properties are considered by TCMGIS-II. Moreover, TCMGIS-II defines the native habitats of a target plant through specimen examination and extracts the target variables of native habitats from its databases.
The present study aims to determine (1) the most important ecological factor(s) on the distribution of R. tanguticum, (2) whether the prediction results are consistent with survey data and (3) the implications of the prediction results for the conservation planning of R. tanguticum.
Based on a spatially referenced GIS model, TCMGIS-II integrated four databases, including the national geographic, climate and soil type databases of China which were used to generate distribution models and the database of resource survey which was used to assess the quality of a model.
The geographic database of China was a digital chart (scale 1:1,000,000) at national, provincial, regional and county levels, including a series of vector maps of layers, i.e. manuals on roads, contours, geology and administrative boundaries, with all points covered with a geographic coordinate system (e.g. latitude, longitude and elevation).
The climate database of China was derived from the national climate data coving from the period of 1971 to 2000 extracted from the climate records of the state meteorological administration of China. The database included climate attributes related to plant growth, e.g. sunshine duration, relative humidity, annual precipitation, accumulated temperature, mean annual temperature, mean March temperature, annual maximum/minimum temperature and annual mean maximum/minimum temperature. The climate data were available in GIS along with data of latitude, longitude and elevation.
The soil type database of China covered a total of 2,444 counties, containing a series of vector soil maps (scale 1:1,000,000) and soil attributes and mapping unit boundaries. The soil data were classified into 12 orders, 29 suborders, 61 groups, 235 subgroups and 909 families as the basic elements of the map layers .
The database of resource survey was generated with the third national resource survey of CM in China, covering a total of 11,118 plant species in 2312 genera of 385 families, including 298 fungi, 114 algae, 43 mosses, 55 lichens, 455 ferns, 126 gymnosperms and 10,027 angiosperms , as well as descriptions on the abundance and distribution patterns of 138 rare and endangered medicinal plants, 126 of which were converted into digital charts (scale 1:1,000,000).
TCMGIS-II identified, analyzed and displayed geographically referenced information, using two major data models (i.e. raster and vector). Raster model in 1.0 × 1.0 km2 grids detected the grids sharing similar ecological factors with those of the native habitats of a target medicinal plant. Vector model stacked the layers of those factors to determine the distribution areas and ranges.
Extraction of ecological factors from native habitats
Variables extracted from the native habitats of Rheum tanguticum Maxim. ex Balf based on TCMGIS-II combined geographic, climate and soil type databases
Mean ± SE
C v (%)
3630 ± 44
63.7 ± 2.2
2450 ± 13
574 ± 7
9517 ± 951
Mean annual temperature
8.6 ± 0.1
Mean March temperature
-4.5 ± 0.2
-19.1 ± 0.2
17.2 ± 0.2
Mean minimum temperature
-11.2 ± 0.2
Mean maximum temperature
10.4 ± 0.2
6.8 ± 0.1
Data normalization and distance analysis
Where x ij is the grid value and y ij is a target variable.
When q = 1, it is Manhattan distance.
When q = 2, it is Euclidean distance.
Long distance indicates low similarity rates while short distance indicates high similarity rates.
Spatial distribution division and model quality assessment
Division on spatial distribution of R. tanguticum was established according to the grid-based clustering. The areas sharing similar ecological factors with those of native habitats were favorable for R. tanguticum distribution. The spatially predicted areas were divided into three types, namely the favorable (with similarity rate ≥95%), suitable (with similarity rate 90-95%), and slightly appropriate (with similarity rate < 90%) for R. tanguticum distribution.
To assess the reliability of the spatial prediction on R. tanguticum distribution, we employed the database of resource survey as a measure. The overlapping part between distribution range predicted by TCMGIS-II and that recorded by resource survey indicates the congruency, the part with prediction result without survey data suggests the potential distribution of R. tanguticum, and the rest part with survey data beyond prediction result indicates the contradiction between prediction result and survey data.
Where σ is the standard deviation and μ is the mean.
Target variables extracted from native habitats
TCMGIS-II extracted the target variables from 206 plots in the native habitats of R. tanguticum (Figure 1, Table 1). The results showed that the target variables varied significantly among different native habitats (Table 1, P < 0.001), with coefficient of variation ranging from 7.6% in sunshine duration to 143.4% in accumulated temperature, and the native habitats exhibited high elevation and abundant sunshine with moderate cool and dry climate in mild acid and basic soils (Table 1). Using PCA, we extracted two principal components (PCs) which accounted for 93.8% of the contribution of target variables in terms of R. tanguticum distribution (Figure 2). The PC1 (PC1 = 60.3%) was mainly related to temperatures (e.g. annual maximum, annual mean maximum, mean annual and acuminated temperatures) and the PC2 (PC2 = 33.5%) was mainly contributed by annual precipitation and relative humidity. However, elevation and annual precipitation were negatively correlated to PC1, and sunshine duration was negatively contributed to PC2 (Figure 2).
Prediction result of potential distributions
Comparison between prediction results and survey data
The ecological factors from native habitats suggest that R. tanguticum grows at high plateau (e.g. alpine meadow, grassland and shrub) with cool climate, abundant sunshine, moderate precipitation and basic soils (e.g. humus-rich loam and sandy loam) and that its distribution is mainly influenced by temperature (e.g. annual maximum, mean annual and acuminated temperatures), annual precipitation and relative humidity. The prediction results by TCMGIS-II confirmed the distribution data.
Many plant species have evolved to be habitat-specific and sensitive to environmental conditions , and those growing at the sites congruent with their native habitats are the most potent . For example, R. tanguticum from Gansu and Qinghai is recorded as a source of rhubarb in the Chinese Pharmacopoeia due to its high potency [7, 32]. The present study found that a large portion of predictive distributions were beyond what survey data covered (e.g. Xinjiang, Inner-Mongolia and Shanxi provinces), agreeing with the notion that prediction of distribution may help locate habitats for conservation [24, 36], giving insights into the discovery of potential habitats for R. tanguticum cultivation.
Interestingly, a small portion of survey data does not overlap with prediction result, e.g. Muli in Sichuan and Zhongdian in Yunnan. According to the Chinese Pharmacopoeia, there are three prescribed sources (i.e. R. tanguticum, R. palmatum and R. officinale) for rhubarb . The survey data cover the three Rheum species. On the other hand, the databases of TCMGIS-II include many abiotic factors (e.g. topographic features, climate conditions and soil properties) but not the effects of dynamic biotic interactions and species-specific features on a large scale. Many plant species are sensitive to both abiotic and biotic factors, such as competitor plants and symbiotic species [37, 38].
In the present study, the distribution of R. tanguticum predicted by TCMGIS-II program was confirmed by the resource survey data. We expect that the TCMGIS-II modeling is useful in conservation planning and regional management for the threatened medicinal plants . Both conservation and sustainable utilization of medicinal plants require robust large-scale assessment of their distribution and regionalization . Lack of data and limit of model validity are barriers for the studies on distribution of medicinal plants on a large scale . Thus, more data and model verification are necessary for further studies and GIS developments.
TCMGIS-II program was confirmed to be useful in the discovery of potential habitats congruent with the native habitats of target medicinal plants. This technology provides reliable references for the conservation planning and regional management of endangered and threatened medicinal plant resources.
geographic information system
a GIS-based program for the distribution prediction of traditional Chinese medicine
the improved version of TCMGIS-I program
- C v :
coefficient of variation
- One-way ANOVA:
one-way analysis of variance
principal components analysis
a principal component
We thank Prof Yulin Lin for his help with specimen identification, Prof Chengzhong Sun and Prof Runhuai Zhao for their constructive comments on TCMGIS-II modeling, Prof Christine Leon and Mr Chun Un for their assistance in polishing the manuscript. This study was supported by the National Key Technology R&D Programs in the 11th Five-year Plan of China (2006BAI09B02, 2006BAI21B07) and China Post-doctoral Foundation (20090450329).
- Russell-Smith J, Karunaratne NS, Mahindapala R: Rapid inventory of wild medicinal plant populations in Sri Lanka. Biol Conserv. 2006, 132: 22-32. 10.1016/j.biocon.2006.03.009.View ArticleGoogle Scholar
- Hamilton AC: Medicinal plants, conservation and livelihoods. Biodivers Conserv. 2004, 13: 1477-1517. 10.1023/B:BIOC.0000021333.23413.42.View ArticleGoogle Scholar
- Mentreddy SR: Review - medicinal plant species with potential antidiabetic properties. J Sci Food Agric. 2007, 87: 743-750. 10.1002/jsfa.2811.View ArticleGoogle Scholar
- Reddy KN, Reddy CS: First red list of medicinal plants of Andhra Pradesh, India - conservation assessment and management planning. Ethnobot Lealf. 2008, 12: 103-107.Google Scholar
- Wu Z, Raven PH, Garden MB: Flora of China. 1994, Beijing: Science PressGoogle Scholar
- Wang XM, Ren Y: Rheum tanguticum, an endangered medicinal plant endemic to China. J Med Plants Res. 2009, 3: 1195-1203.Google Scholar
- Chinese Pharmacopoeia Committee: Pharmacopoeia of the People's Republic of China. 2010, Beijing: Chinese Medical Technological PressGoogle Scholar
- Li M, Li LX, Liu Y: Study survey on rhubarb in recent years. World Sci Tech/Mod Trad Chin Med. 2006, 8: 34-39.Google Scholar
- Meilleur BA, Hodgkin T: In situ conservation of crop wild relatives: status and trends. Biodivers Conserv. 2004, 13: 663-684. 10.1023/B:BIOC.0000011719.03230.17.View ArticleGoogle Scholar
- Shinwari ZK, Gilani SS: Sustainable harvest of medicinal plants at Bulashbar Nullah, Astore (Northern Pakistan). J Ethnopharmacol. 2003, 84: 289-298. 10.1016/S0378-8741(02)00333-1.View ArticlePubMedGoogle Scholar
- Barazani O, Perevolotsky A, Hadas R: A problem of the rich: prioritizing local plant genetic resources for ex situ conservation in Israel. Biol Conserv. 2008, 141: 596-600. 10.1016/j.biocon.2007.10.014.View ArticleGoogle Scholar
- Rita A, Silvano M: Ex situ conservation and rare plants propagation in the Lecce botanical garden: reproductive biology problems. Caryologia. 2006, 59: 345-349.Google Scholar
- Chen SL, Wei JH, Huang LF, Guo BL, Xiao PG: Probing into the theory and practice of wild medicinal materials tending. China J Chin Mater Med. 2004, 29: 1123-1126.Google Scholar
- Long CL, Li H, Ouyang ZQ, Yang XY, Li Q, Trangmar B: Strategies for agrobiodiversity conservation and promotion: a case from Yunnan, China. Biodivers Conserv. 2003, 12: 1145-1156. 10.1023/A:1023085922265.View ArticleGoogle Scholar
- Williams JN, Seo CW, Thorne J, Nelson JK, Erwin S, O'Brien JM, Schwartz MW: Using species distribution models to predict new occurrences for rare plants. Divers Distrib. 2009, 15: 565-576. 10.1111/j.1472-4642.2009.00567.x.View ArticleGoogle Scholar
- Ginnett TF, Demment MW: Sexual segregation by Masai giraffes at two spatial scales. Afr J Ecol. 1999, 37: 93-106. 10.1046/j.1365-2028.1999.00163.x.View ArticleGoogle Scholar
- Ghimire SK, McKey D, Aumeeruddy-Thomas Y: Heterogeneity in ethnoecological knowledge and management of medicinal plants in the Himalayas of Nepal: implications for conservation. Ecol Soc. 2004, 9: 9-24.Google Scholar
- Pulliam HR: On the relationship between niche and distribution. Ecol Lett. 2000, 3: 349-361. 10.1046/j.1461-0248.2000.00143.x.View ArticleGoogle Scholar
- Aguilar-Stoen M, Moe SR: Medicinal plant conservation and management: distribution of wild and cultivated species in eight countries. Biodivers Conserv. 2007, 16: 1973-1981. 10.1007/s10531-006-9125-7.View ArticleGoogle Scholar
- Whelan RJ, De Jong NH, Von der Burg S: Variation in bradyspory and seedling recruitment without fire among populations of Banksia serrata (Proteaceae). Aust J Ecol. 1998, 23: 121-128. 10.1111/j.1442-9993.1998.tb00710.x.View ArticleGoogle Scholar
- Rodriguez JP, Brotons L, Bustamante J, Seoane J: The application of predictive modelling of species distribution to biodiversity conservation. Divers Distrib. 2007, 13: 243-251. 10.1111/j.1472-4642.2007.00356.x.View ArticleGoogle Scholar
- Guisan A, Thuiller W: Predicting species distribution: offering more than simple habitat models. Ecol Lett. 2005, 8: 993-1009. 10.1111/j.1461-0248.2005.00792.x.View ArticleGoogle Scholar
- Arundel ST: Using spatial models to establish climatic limiters of plant species' distributions. Ecol Model. 2005, 182: 159-181. 10.1016/j.ecolmodel.2004.07.016.View ArticleGoogle Scholar
- Chefaoui RM, Hortal J, Lobo JM: Potential distribution modelling, niche characterization and conservation status assessment using GIS tools: a case study of Iberian Copris species. Biol Conserv. 2005, 122: 327-338. 10.1016/j.biocon.2004.08.005.View ArticleGoogle Scholar
- 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 Trad Herb Drugs. 2007, 38: 481-487.Google Scholar
- Sun CZ, Liu ZQ, Chen SL, Wei JH, Zhao RH, Wang JY, Zhou YQ: Design and realization of traditional Chinese medicine adaptability analyzing system based on GIS. World Sci Tech/Mod Trad Chin Med. 2006, 8: 112-117.Google Scholar
- Wang Y, Wei JH, Chen SL, Sun CZ, Zhao RH, Liu ZQ, Xiao XH, Wang JY, Zhou YQ: Analysis of Panax Ginseng's producing area based on TCMGIS-I. Asia Pac Trad Med. 2006, 6: 73-78.Google Scholar
- Chen SL, Zhou YQ, Xie CX, Zhao RH, Sun CZ, Wei JH, Liu ZQ, Gao WW: Suitability evaluation of Panax quinquefolium's producing area based on TCMGIS-I. Zhongguo Zhong Yao Za Zhi. 2008, 33: 741-745.PubMedGoogle Scholar
- Wang JY, Zhao RH, Sun CZ, Chen SL, Liu ZQ, Wei JH, Xiao XH, Zhou YQ: Suitability evaluation of Glycyrrhiza uralensis Fisch.'s distributive area based on TCMGIS-I. Mod Chin Med. 2006, 8: 4-8.Google Scholar
- Huang LF, Xie CX, Duan BZ, Chen SL: Mapping the potential distribution of high artemisinin-yielding Artemisia annua L. (Qinghao) in China with a geographic information system. Chin Med. 2010, 5: 18-25. 10.1186/1749-8546-5-18.PubMed CentralView ArticlePubMedGoogle Scholar
- Shi XZ, Yu DS, Warner ED, Pan XZ, Petersen GW, Gong ZG, Weindorf DC: Soil database of 1: 1,000,000 digital soil survey and reference system of the Chinese genetic soil classification system. Soil Sur Horiz. 2004, 45: 129-136.View ArticleGoogle Scholar
- China National Group Corporation of Traditional and Herbal Medicine: Regionalization of Traditional Chinese Medicine. 1995, Beijing: Science PressGoogle Scholar
- Warren LT: Clustering of time series data - a survey. Pattern Recognit. 2005, 38: 1857-1874. 10.1016/j.patcog.2005.01.025.View ArticleGoogle Scholar
- Bedeian AW, Mossholder KW: On the use of the coefficient of variation as a measure of diversity. Organ Res Methods. 2000, 3: 285-297. 10.1177/109442810033005.View ArticleGoogle Scholar
- Heinanen S, von Numers M: Modelling species distribution in complex environments: an evaluation of predictive ability and reliability in five shorebird species. Divers Distrib. 2009, 15: 266-279. 10.1111/j.1472-4642.2008.00532.x.View ArticleGoogle Scholar
- Chen SL, Xie CX, Yao H, Song JY, Sun C: Methodological innovations in traditional Chinese medicine resources. World Sci Tech/Mod Trad Chin Med. 2008, 10: 1-9.Google Scholar
- Kala CP: Status and conservation of rare and endangered medicinal plants in the Indian trans-Himalaya. Biol Conserv. 2000, 93: 371-379. 10.1016/S0006-3207(99)00128-7.View ArticleGoogle Scholar
- Dhyani PP, Kala CP: Current research on medicinal plants: five lesser known but valuable aspects. Curr Sci. 2005, 88: 335-Google Scholar
- Thuiller W, Albert C, Araujo MB, Berry PM, Cabeza M, Guisan A, Hickler T, Midgely GF, Paterson J, Schurr FM, Sykes MT, Zimmermann NE: Predicting global change impacts on plant species' distributions: future challenges. Perspect Plant Ecol Evol Syst. 2008, 9: 137-152. 10.1016/j.ppees.2007.09.004.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.