Spatial Distribution Pattern of the Populations of Carex siderosticta at Mt. Geumjeong and Mt. Ahop
Spatial Distribution Pattern of the Populations of Carex siderosticta at Mt. Geumjeong and Mt. Ahop
Journal of Life Science. 2015. Apr, 25(4): 369-375
Copyright © 2015, Korean Society of Life Science
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : November 07, 2014
  • Accepted : March 04, 2015
  • Published : April 30, 2015
Export by style
Cited by
About the Authors
만규 허

식물 보존의 권고지역에서 식물 집단의 공간적 데이터는 여러 목적에서 중요하다. 부산광역시 금정산과 아홉산 의 대사초( Carex siderosticta )의 지리적 거리에 따른 공간적 분포를 분석하였다. 공간적 양상의 분석 방법은 여러 패치 척도, 분산 척도에 의거한 plot(플롯)의 크기에 따라 집단의 균질성 또는 운집을 분석하였다. 대사초의 많은 자연 플롯은 산림군락에서 균질하지 않았다. 예를 들면 균질한 플롯은 6.0 m × 6.0 m 이내였다. 플롯의 크기가 6.0 m × 12.0 m 이상이면 운집되었다. 대사초의 이웃 패치는 평균 7.5 m에서 9.0 m사이였다. 그런데 자연집단이 인간의 활동에 의해 교란되면 운집은 9.0 m보다 짧은 거리에서 일어났다. 결론적으로 대사초의 지리적 분포는 플롯의 밀도에 균질하지 않았고 금정산과 아홉산 집단에서 집단 크기에 따라 다르며 인간의 간섭은 플롯에서 밀도 효과를 일으킨다.
Knowledge of specific habitats is also critical in the study of the distribution of plants. Habitats include a combination of physical factors that represent the environmental conditions in which organisms live [10] . Each plant community was considered as a discrete patch occupied by individuals of different species from a limited regional pool because all species were assumed to be in competition with each other. The spatial prediction of species distributions from survey data has recently been recognized as a significant component of conservation planning [8 , 18] .
Quantitative examination of spatially explicit data in ecology is broadly categorized as “spatial analysis” [13 , 15] . The analysis of the spatial pattern of individuals of a particular species has long been a concern of ecologists [21] .
Many ecologists have adopted several different major schools of spatial analysis from other disciplines [15] . The first of these comes from geography, and its methods include the use of statistics (e.g., Moran’s I ) to measure spatial auto-correlation [6 , 19 , 20] . It measures the degree to which the occurrence of an event in a real unit constrains, or makes more probable, the occurrence of an event in a neighboring areal unit.
Over the last decade, statistical ecologists have developed a number of elegant tools to incorporate spatial (or temporal) variables in analyses of multivariate ecological data sets. The presence of spatial (or temporal) autocorrelation in community composition data can be tested using multivariate Mantel correlograms. Variation partitioning [3 , 17] provides a method for distinguishing the separate (and combined) influences of environmental, spatial, and temporal variables on the variability of multi-species distribution and abundance data.
Carex siderosticta Hance is an herbaceous and belongs to the family, Cyperaceae. The fruits are achene and yellow-flower blooms July–August in forests. Flowering culms and vegetative culms spaced; flowering culms clothed by bladeless sheaths at base, pale brown, without leaves. Leaf blades of vegetative culms are oblong-lanceolate, sometimes with white stripes. Flowering culms grow up to 30 cm tall. C. siderosticta is a species of sedge native to East Asia [23] .
In this report, the several statistical tools of percentage distribution and population structure of the geographical areas are used to study the spatial distribution of C. siderosticta in Busan.
Mt. Geumjeong and Mt. Ahop locate in south of the Korean. A sample of a large (more than 500 individuals) natural population of C. siderosticta collected at both mountains and was used in this study. It is expected to provide useful experimental conditions because of the large undisturbed and isolated site. Most temperate Carex species including C. siderosticta have shoots formed during the previous year, some emerging in autumn, others remaining below ground until spring [2] . The maximum shoot life span for temperate species appears to be approximately 24 months but mortality is very high; sometimes 90% of shoots do not live for the whole 2-year life span. Mortality is caused by differences in time of emergence, flowering, animal grazing, the age of the genet, and internal competition through the rhizome system. Thus, C. siderosticta is ideal species to study the spatial distribution pattern of the population levels.
The purpose of this paper was to describe a statistical analysis for detecting a species association, which is valid even when the assumption of within- species spatial randomness is violated. The purpose of this study is addressed: is there a spatial structure within two populations of C. siderosticta ? and 2) if so, what is the spatial pattern and is it the same for all populations?
Materials and Methods
- Study area
I conducted the spatial analysis in the communities of Carex siderosticta at Mt. Geumjeong in Busan-si. The mountain (801.5 m) is highest in Busan. Mt. Ahop (350.0 m) locates in Gijang-gun, Busan-si, Korea. Mt. Geumjeong and Mt. Ahop is about 20 km away. It has a temperate climate with a little hot and long summer. In this region the mean annual temperature is 14.7℃ with the maximum temperature being 29.4℃ in August and the minimum −0.6℃ in January. Mean annual precipitation is about 1519.1 mm with most rain falling period between June and August.
- Sampling procedure
Spatial ecologists use artificial sampling units (so-called quadrats) to determine abundance or density of species. The number of events per unit area are counted and divided by area of each square to get a measure of the intensity of each quadrat. I established many 1.5 m × 1.5 m quadrats with an area of 12 m × 12 m each around one area at Mt. Ahop and Mt. Geumjeong. I randomly located quadrates in each plot which I established populations. The quadrat sizes were 1.5 m × 1.5 m, 1.5 m × 3 m, 3 m × 3 m, 3 m × 6 m, 6 m × 6 m, 6 m × 12 m, and 12 m × 12 m. I mapped all plants to estimate C. siderosticta density per plot.
- Index calculation and data analysis
The spatial pattern of C. siderosticta was analyzed according to the Neatest Neighbor Rule [5 , 15] with Microsoft Excel 2010.
Average viewing distance ( rA ) was calculated as follows:
PPT Slide
Lager Image
Where ri is the distance from the individual to its nearest neighbor. N is the total number of individuals within the quadrat.
The expectation value of mean distance of individuals within a quadrat ( rB ) was calculated as follows:
PPT Slide
Lager Image
Where D is population density and D is the number of individuals per plot size.
R = rA / rB
When R > 1, it is a uniform distribution, R = 1, it is a random distribution, R < 1, it is an aggregated distribution.
The significance index of the deviation of R that departs from the number of “1” is calculated from the following formula [14] .
PPT Slide
Lager Image
PPT Slide
Lager Image
When CR > 1.96, the level of the significance index of the deviation of R is 5%, and When CR > 2.58, the level is 1%.
I calculated the degree of population aggregation under different sizes of plots by dispersion indices: index of clumping or the index of dispersion (C), aggregation index ( CI ), mean crowding (M*), patchiness index (PAI), negative binominal distribution index K , Ca indicators ( Ca is the name of one index) [16] and Morisita index (IM) were calculated with Microsoft Excel 2010. The formulae are as follows:
  • Index of dispersion:C=S2/m
  • Aggregation index
  • Mean crowding
  • Patchiness index
  • Aggregation intensity
  • Ca indicatorsCa= 1/k
Where S 2 is variance and m is mean density of C. siderosticta .
When C . M *. PAI > 1, it means aggregately distributed, when C . M *. PAI < 1, it means uniformly distributed, when CI , PA , Ca > 0, it means aggregately distributed, and when CI , PA , Ca < 0 it means uniformly distributed.
The mean aggregation number to find the reason for the aggregation of C. siderosticta was calculated [1] .
δ = mr / 2 k
Where r is the value of chi-square when 2k is the degree of freedom and k is the aggregation intensity
- Spatial structure
Numerical simulations of previous analyses were performed to investigate the significant differences at various distance scales, i.e., 1.0 m, 1.5 m, 2.0 m, and so on. However, no significant population structure was found within the 1.5 m distance classes by means of Moran’s I , and a significant population structure was revealed beyond 1.5-m. Thus, the distance classes are 0-1.5 m (class I), 1.5-3.0 m (class II), 3.0-4.5 m (class III), 4.5-6.0 m (class IV), 6.0-7.5 m (class V), 7.5-9.0 m (class VI), 9.0-10.5 m (class VII), 10.5-12.0 m (class VIII), 12.0-13.5 m, 13.5-15.0 m (class IX), and 15.0-16.5 m (class X). The codes of classes are the same as in the distance classes and are listed Table 1 .
Spatial patterns ofCarex siderostictaindividuals at different sampling quadrat sizes in Mt. Geumjeong and Mt. Ahop
PPT Slide
Lager Image
Spatial patterns of Carex siderosticta individuals at different sampling quadrat sizes in Mt. Geumjeong and Mt. Ahop
The spatial structure was quantified by Moran’s I , a coefficient of spatial autocorrelation (SA) [19] . As applied in this study, Moran’s I quantifies the similarity of pairs of spatially adjacent individuals relative to the population sample as a whole. The value of I ranges between +1 (completely positive autocorrelation, i.e., paired individuals have identical values) and −1 (completely negative autocorrelation). Each plant was assigned a value depending on the presence or absence of a specific individual. If the i th plant was a homozygote for the individual of interest, the assigned pi value was 1, while if the individual was absent, the value 0 was assigned [20] .
Pairs of sampled individuals were classified according to the Euclidian distance, dij, so that class k included dij satisfying k − 1 < dij < k + 1, where k ranges from 1 to 10. The interval for each distance class was 1.5 m. Moran’s I statistic for class k was calculated as follows:
I ( k ) = n∑i∑j ( i≠j ) WijZiZj/S∑Zi 2
where Zi is pi − p (p is the average of pi ); Wij is 1 if the distance between the i th and j th plants is classified into class k; otherwise, Wij is 0; n is the number of all samples and S is the sum of Wij {∑i∑j(i ≠j)Wij} in class k. Under the randomization hypothesis, I (k) has the expected value u1 = −1/(n − 1) for all k. Its variance, u2, has been given, for example, in Sokal and Oden [19] . Thus, if an individual is randomly distributed for class k, the normalized I (k) for the standard normal deviation (SND) for the plant genotype, g (k) = { I (k) − u1}/u2 1/2 , asymptotically has a standard normal distribution [6] . Hence, SND g(k) values exceeding 1.96, 2.58, and 3.27 are significant at the probability levels of 0.05, 0.01, and 0.001, respectively.
- The spatial pattern of individuals
Population densities (D) varied from 0.547 to 5.778, with a mean of 2.644 ( Table 1 ). The D value of Mt. Geumjeong area (3.108) is higher than Mt. Ahop area (2.180). There was shown significant difference between both mountain areas. The values ( R ) of spatial distance (the rate of observed distance-to-expected distance) among the nearest individuals were higher than 1 and the significant index of R (Cr) was > 2.58. When this parameter was applied to two areas, the small plots (1.5 m × 1.5 m, 1.5 m × 3 m, 3 m × 3 m, and 3 m × 6 m) of C. siderosticta were uniformly distributed in the forest community ( Table 1 ). However, C. siderosticta were aggregately distributed in large plots (6 m × 12 m, 12 m × 12 m, and 12 m × 24 m) ( Table 1 )
- The degree of population aggregation
Dispersion index ( C ) were higher than 1 except two quadrats (1.5 m × 1.5 m and 1.5 m × 3 m for Mt. Geumjeong area) ( Table 2 ). Thus aggregation indices ( CI ) were positive except two plots at Mt. Geumjeong which indicate a clumped distribution. The mean crowding (M*) and patchiness index (PAI) showed positive values. In Mt. Ahop, the three indices, C , M *, PAI were >1 and their values of PI and Ca except two plots were also shown greater than zero, thus it means aggregately distributed. The most individuals of C. siderosticta were clustered and the distribution pattern of the C. siderosticta was quadrat-sampling dependent. As the sizes of quadrat were greater, the PI values of C. siderosticta showed high.
Changes in gathering strength ofCarex siderostictaat different sampling quadrat sizes
PPT Slide
Lager Image
Changes in gathering strength of Carex siderosticta at different sampling quadrat sizes
Morisita index ( IM ) is related to the patchiness index ( PAI ) and showed an overly steep slope at the plot 3 m x 3 m in Mt. Geumjeong and at the plot 1.5 m x 3 m in Mt. Ahop. When the area was smaller than 1.5 m x 3 m, the degree of aggregation increased significantly with increasing quadrat sizes, while the patchiness indices did not change from the plot 6 m × 6 m to 12 m × 24 m.
The mean aggregation number (δ) analysis showed that the reasons for aggregation of C. siderosticta differed in quadrats with different plot sizes. The most clusters at 25 quadrat was determined by environmental factors. When the size was one 2.5 m × 5 m plot at west, the cluster was determined by both species characteristics and environmental factors.
- Analysis of spatial autocorrelation
The spatial autocoefficient, Moran’s I is presented in Table 3 . Separate counts for each type of joined individuals and for each distance class of separation were tested for significant deviation from random expectations by calculating the SND. Moran’s I of C. siderosticta significantly differed from the expected value in only 16 of 40 cases (40%). Five of these values (31.3%) were negative, indicating a partial dissimilarity among pairs of individuals in the 10 distance classes. Eleven of the significant values (68.7%) were positive, indicating similarity among individuals in the first 4 distance classes, i.e., pairs of individuals can separate by more than 10 m. Namely, significant aggregations were partially observed within IV classes. As a matter of course, the negative SND values at classes VI, VII, VIII, and X. Thus, dissimilarity among pairs of individuals could found by more than 15.5 m.
Spatial autocorrelation coefficients (Moran’sI) among two populations ofCarex siderostictafor ten distance classes
PPT Slide
Lager Image
*: p<0.05, **: p<0.01, ***: p<0.001.
The comparison of Moran’s I values to a logistic regression indicated that a highly significant percentage of individual dispersion in C. siderosticta populations at Mt. Geumjeong could be explained by isolation by distance.
When R = 1, it is a random distribution; R < 1, it is an aggregation; R > 1, it is a uniform distribution [14] . According to this rule, all plots and areas of C. siderosticta at Mt. Geumjeong are uniform distribution ( Table 1 ). However, According to dispersion indices of Llord [16] , many plots are not uniform distribution ( Table 2 ) and not consistent with the rule. R = rA / rB [5] . N for rA is total numbers of within the plot and rB is concerned with plot size. Although, a large plot has large N, the plot size is not N. As D is the number of individuals per plot size, the nearest neighbor rule by Clark and Evans [5] is good for spatial pattern. 16 plots (66.7%) showed were uniform distributed. In only 8 plots (33.3%), the three indices, C , M *, PAI were >1, and PI and Ca > 0, thus it means aggregately distributed. Aggregation is mainly caused by the environmental factors [14] . When δ > 2, the aggregation was mainly caused by both species characteristics and environmental factors [14] . Most 23 plots except one had low δ < 2. I recognized that the important environmental factors might be considered competition, growth rate, little decomposition, light, and below-ground resources. The characteristics of the C. siderosticta concerned included primarily their life history, artificial disturbance, and population density. Life history theory seeks to understand the variation in traits such as growth rate, number and size of offspring and life span observed in nature, and to explain them as evolutionary adaptations to environmental conditions [22] . Artificial disturbances are important environmental factors affecting C. siderosticta such as constitutional roads at east area, temple construction at west area, and farming at south area. At the plots which had fewer C. siderosticta , the cluster was mainly determined by C. siderosticta themselves. Although the small proportion of seeds from mother removed by ants, 99% seeds in C. siderosticta fell within 20 cm of the scape and 91% within 10 cm of the scape [12] . In addition, the mean value of the aggregation index changed irregularly with population growth rate [4] .
PPT Slide
Lager Image
The curves of patchiness according to geographic distances for two communities of Carex siderosticta using values of Green index. X axis is quadrat sizes (m × m).
PPT Slide
Lager Image
The changes of the mean aggregation numbers for communities of Carex siderosticta. X axis is quadrat sizes (m × m).
A significant positive value of Moran’s I indicated that pairs of individuals separated by distances that fell within distance class V had similar individuals, whereas a significant negative value indicated that they had dissimilar individuals ( Table 3 ). The overall significance of individual correlograms was tested using Bonferroni’s criteria. The results revealed that patchiness similarity was shared among individuals within up to a scale of a 7.5 m~10 m distance. Thus it was looked for the presence of dispersion correlations between neighbors at this scale.
The results from this study are consistent with the supposition that a plant population is subdivided into local demes, or neighborhoods of related individuals [7 , 11] . Previous reports on the local distribution of genetic variability suggested that microenvironmental selection and limited gene flow are the main factors causing substructuring of alleles within a population [9] .
In conclusion, C. siderosticta populations within Mt. Geumjeong was observed a strong spatial structure. Neighboring patches of C. siderosticta are predominantly 7.5 m to 10 m apart on average. The present study demonstrates that a spatial structure of C. siderosticta in the Mt. Geumjeong populations could be explained by isolation by distance, limited gene flow, and topography. However, if the natural populations were disturbed by human activities, the aggregation was occurred in more short distance than a scale of a 7.5 m~10 m distance. The results of this study were used as systematic conservation planning which is an effective way to seek and identify efficient and effective types of reserve design to capture or sustain the highest priority biodiversity values and to work with communities in support of local ecosystems. Conservation biology is an objective science when biologists advocate for an inherent value in nature.
Arbous A. G. , Kerrich J. E. 1951 Accident statistics and the concept of accident proneness Biometrics 7 340 - 342    DOI : 10.2307/3001656
Bernard J. M. 1990 Life history and vegetative reproduction in Carex Can. J. Bot. 68 1441 - 1448    DOI : 10.1139/b90-182
Borcard D. , Legendre P. 2012 Is the Mantel correlogram powerful enough to be useful in ecological analysis? A simulation study Ecology 93 1473 - 1481    DOI : 10.1890/11-1737.1
Caswell1 H. , Takenori Takada T. , Hunte C. M. 2004 Sensitivity analysis of equilibrium in density-dependent matrix population models Letters 7 380 - 387
Clark P. J. , Evans F. C. 1954 Distance to nearest neighbor as a measure of spatial relationships in populations Ecology 35 445 - 453    DOI : 10.2307/1931034
Cliff A. D. , Ord J. K. 1971 Spatial autocorrelation Pion London
Ehrlich P. R. , Raven P. H. 1969 Differentiation of populations Science 165 1228 - 1232    DOI : 10.1126/science.165.3899.1228
Elith J. , Burgman M. A. , Brigham C. A. , Schwanz M. W 2002 Population Viability in Plants Springer Habitat models for PVA
Epperson B. K. , Chung M. G. , Telewski F. W. 2003 Spatial pattern variation in a contact zone of Pinus ponderosa and P. arizonica (Pinaceae) Am. J. Bot. 90 25 - 31    DOI : 10.3732/ajb.90.1.25
Evans K. L. , Duncan R. P. , Blackburn T. M. , Crick H. Q. P. 2005 Investigating geographic variation in clutch size using a natural experiment Funct. Ecol. 19 616 - 624    DOI : 10.1111/j.1365-2435.2005.01016.x
Garnier L. K. M. , Durand J. , Dajoz I. 2002 Limited seed dispersal and microspatial population structure of an agamospermous grass of West African savannahs. Hyparrhenia diplandra (Poaceae) Am. J. Bot. 89 1785 - 1791    DOI : 10.3732/ajb.89.11.1785
Guitián P. , Medrano M. M. , Guitián J. 2003 Seed dispersal in Erythronium dens-canis L. (Liliaceae): variation among habitats in a myrmecochorous plant Plant Ecology 169 171 - 177    DOI : 10.1023/A:1026043411357
Legendre P. , Fortin M. J. 1989 Spatial pattern and ecological analysis Vegetatio 80 107 - 138    DOI : 10.1007/BF00048036
Lian X. , Jiang Z. , Ping X. , Tang S. , Bi J. , Li C. 2012 Spatial distribution pattern of the steppe toad-headed lizard (Phrynocephalus frontalis) and its influencing factors Asian Herpetol. Res. 3 46 - 51    DOI : 10.3724/SP.J.1245.2012.00046
Liebhold A. M. , Gurevitch J. 2002 Integrating the statistical analysis of spatial data in ecology Ecography 25 553 - 557    DOI : 10.1034/j.1600-0587.2002.250505.x
Lloyd M. 1967 Mean crowding J. Ani. Ecol. 36 1 - 30    DOI : 10.2307/3012
Peres-Neto P. R. , Legendre P. , Dray S. , Borcard D. 2006 Variation partitioning of species data matrices: estimation and comparison of fractions Ecology 87 2614 - 2625    DOI : 10.1890/0012-9658(2006)87[2614:VPOSDM]2.0.CO;2
Scott J. M. , Heglund P. J. , Samson F. , Haufler J. , Morrison M. , Raphael M. , Wall B. 2002 Predicting Species Occurrences: Issues of Accuracy and Scale Island Press Covelo, CA 868 -
Sokal R. R. , Oden N. L. 1978 Spatial autocorrelation in biology 1 Methodology. Biol. J. Linn. Soc. 10 199 - 228    DOI : 10.1111/j.1095-8312.1978.tb00013.x
Sokal R. R. , Oden N. L. 1978 Spatial autocorrelation in biology 2. Some biological implications and four applications of evolutionary and ecological interest Biol. J. Linn. Soc. 10 229 - 249    DOI : 10.1111/j.1095-8312.1978.tb00014.x
Sone K. , Ohkubo K. , Matsuo T. , Hata K. 2013 Spatial distribution pattern of pine trees killed by pine wilt disease in a sparsely growing, young pine stand J. Plant Stud. 2 36 - 41
Souza A. F. , Martins F. R. 2004 Microsite specialization and spatial distribution of Geonoma brevispatha, a clonal palm in south-eastern Brazil Ecol. Res. 19 521 - 532    DOI : 10.1111/j.1440-1703.2004.00670.x
Yan-Cheng T. , Qiu-Yun X. 1989 Cytological studies of Carex siderosticta Hance (Cyperaceae) and its importance in phytogeography Cathaya 1 49 - 60