Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Traits influence detection of exotic plant species in tropical forests

  • Decky I. Junaedi ,

    Contributed equally to this work with: Decky I. Junaedi, Michael A. McCarthy, Gurutzeta Guillera-Arroita

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Software, Visualization, Writing – original draft, Writing – review & editing

    deck001@lipi.go.id, deqee82@gmail.com

    Affiliations Centre of Excellence for Biosecurity Risk Analysis (CEBRA), School of Biosciences, The University of Melbourne, Parkville, Victoria, Australia, School of Biosciences, The University of Melbourne, Parkville, Victoria, Australia, Cibodas Botanic Gardens, Indonesian Institute of Sciences (LIPI), Cipanas, Cianjur, West Java, Indonesia

  • Michael A. McCarthy ,

    Contributed equally to this work with: Decky I. Junaedi, Michael A. McCarthy, Gurutzeta Guillera-Arroita

    Roles Conceptualization, Formal analysis, Methodology, Supervision, Validation, Writing – original draft, Writing – review & editing

    Affiliation School of Biosciences, The University of Melbourne, Parkville, Victoria, Australia

  • Gurutzeta Guillera-Arroita ,

    Contributed equally to this work with: Decky I. Junaedi, Michael A. McCarthy, Gurutzeta Guillera-Arroita

    Roles Conceptualization, Formal analysis, Software, Validation, Writing – original draft, Writing – review & editing

    Affiliation School of Biosciences, The University of Melbourne, Parkville, Victoria, Australia

  • Jane A. Catford ,

    Roles Conceptualization, Funding acquisition, Methodology, Supervision, Writing – original draft, Writing – review & editing

    ‡ These authors also contributed equally to this work.

    Affiliations School of Biosciences, The University of Melbourne, Parkville, Victoria, Australia, Biological Sciences, University of Southampton, Southampton, United Kingdom, Fenner School of Environment and Society, The Australian National University, Canberra, ACT, Australia

  • Mark A. Burgman

    Roles Conceptualization, Formal analysis, Funding acquisition, Methodology, Project administration, Resources, Supervision, Validation, Writing – review & editing

    ‡ These authors also contributed equally to this work.

    Affiliations Centre of Excellence for Biosecurity Risk Analysis (CEBRA), School of Biosciences, The University of Melbourne, Parkville, Victoria, Australia, School of Biosciences, The University of Melbourne, Parkville, Victoria, Australia, Centre for Environmental Policy, Faculty of Natural Sciences, Imperial College London, South Kensington Campus, London, United Kingdom

Abstract

Detecting exotic plant species is essential for invasive species management. By accounting for factors likely to affect species’ detection rates (e.g. survey conditions, observer experience), detectability models can help choose search methods and allocate search effort. Integrating information on species’ traits can refine detectability models, and might be particularly valuable if these traits can help improve estimates of detectability where data on particular species are rare. Analysing data collected during line transect distance sampling surveys in Indonesia, we used a multi-species hierarchical distance sampling model to evaluate how plant height, leaf size, leaf shape, and survey location influenced plant species detectability in secondary tropical rainforests. Detectability of the exotic plant species increased with plant height and leaf size. Detectability varied among the different survey locations. We failed to detect a clear effect of leaf shape on detectability. This study indicates that information on traits might improve predictions about exotic species detection, which can then be used to optimise the allocation of search effort for efficient species management. The innovation of the study lies in the multi-species distance sampling model, where the distance-detection function depends on leaf traits and height. The method can be applied elsewhere, including for different traits that may be relevant in other contexts. Trait-based multispecies distance sampling can be a practical approach for sampling exotic shrubs, herbs, or grasses species in the understorey of tropical forests.

Introduction

Invasive exotic plant species are a priority for environmental management because they threaten native biodiversity [1] and alter ecosystem processes [2]. Once exotic species are established and naturalised, their management can be very costly [3], and the effort needed for eradication increases markedly with species abundance [4]. For these reasons, early and strategic invasive species management is important, and effective detection of exotic species is key to this endeavour.

Information about exotic invasive species detectability supports decisions at different management stages, from early detection for potential control or eradication, to evaluation of eradication or containment. Early detection of exotic species significantly reduces the cost of management [5], given that preventative management is often more efficient than curative management [6]. Given limited available resources for exotic species management [7, 8], the ratio of detection cost to the benefit of eradicating / managing these species is crucial when aiming for efficient invasive species monitoring [9].

Detectability refers to “the quality of being detectable” [10] or the extent to which something is observable. At the species level, detectability can be measured by the probability that the species is detected at a site under a prescribed search strategy (e.g. half an hour of unstructured search), given its presence. Detectability can also be defined at the level of individuals, e.g. the probability than an individual is detected along a line transect based on its (perpendicular) distance from the transect. Accounting for detectability provides various benefits to invasive exotic species management. First, imperfect detection can be accounted for in models that estimate the occurrence and abundance of species [1113] and communities [1416], leading to more reliable information to guide decisions about suitable management actions [17]. Second, estimates of detectability can guide optimal survey effort allocation and strategies for species detection given limited resources [9, 1820].

Vegetation surveys use various methods, including unstructured searches [21], fixed plots, plotless (point-based) methods and nearest neighbour techniques [22]. Distance sampling is another potential approach for surveying plant species. Distance sampling is a method for estimating species abundance that uses data on distance to detected targets collected along line transects, and accounts for imperfect detection by modelling how the probability of detecting an individual decays with its perpendicular distance from the transect [23]. The method has been primarily used to survey animals so far, but can be practical for plant surveys too [24]. Being plotless, it can be simpler to implement than plot-based approaches [25], and yet can be as efficient [26]. Apart from its value for abundance estimation, the detectability information from distance sampling can guide survey effort for rare or exotic plants. Lastly, this method is also suitable for individual level detectability, while several recent studies in tropical ecosystem were mostly conducted for or related to species or community level detectability [2729].

A number of factors can affect plant detectability. For instance, Chen et al. [11] show how site characteristics (e.g. elevation) and survey methods affect the detectability of vascular plant species in Switzerland. Garrard et al. [12] modeled the detectability of plant species in an Australian temperate grassland and showed that the time it takes to detect the first individual of a target species is explained by site characteristics, species abundance, and surveyor experience. In addition, Garrard et al. [12] showed that traits, such as color and flowering time, affect detectability. Indeed, flowers and fruits can be important cues for plant detection, but reproductive organs are not always present when field surveys are conducted. Further study is needed to assess the potential influence of other traits on plant detectability, as well as to expand the scope of detectability studies to other ecosystem types such as tropical forests.

Plant height, leaf size and leaf shape are plausible candidates for explaining variation in detectability among species. These characteristics may be important for detecting invasive plant species, especially in tropical rainforest with dense and complex vegetation structure. Leaf shape, leaf area and plant height are possibly important when distinguishing exotic from native plant species. For example, assuming surveyor can reliably distinguish exotics from natives, it is logical that taller individuals will be easier to detect than shorter ones, especially in dense stands of vegetation like those found in the lower layers of tropical rainforest. Similarly, leaves with distinct shapes or sizes may be easier to detect.

This study tests whether common leaf visual characteristics and plant height correlate with the detectability of individuals of multiple exotic species when surveyed using the distance sampling method in tropical secondary rainforests. Identifying traits that affect plant detectability in general is useful because borrowing explanatory power across species can: (1) increase the accuracy of abundances estimates for species that are detected rarely; and (2) help estimate the effort required to detect a species, even if data on detectability for that species is sparse. The innovation of this study lies in the multi-species detection model where the distance-detection function depends on leaf traits and height, which could serve such applications elsewhere.

Materials and methods

Locations

Line transect surveys were conducted in native Indonesian rainforest ecosystems located next to Cibodas Botanic Gardens (Cibodas), Kuningan Botanic Gardens (Kuningan), Baturraden Botanic Gardens (Baturraden) and Eka Karya Bali Botanic Gardens (Eka Karya) (S1 Fig). These four botanic gardens have different ages: Cibodas was established in 1852 and Eka Karya in 1959, whereas Baturraden and Kuningan are considerably younger having been established in 2004 and 2005, respectively.

This field study was approved and permitted by authorities that managed all the accessed lands in this study. These authorities are: Mount Gede-Pangrango National Park (TNGGP) and Cibodas Botanic Gardens—LIPI (Cibodas); Nature Conservation Agency (BKSDA) of Bali and Eka Karya Bali Botanic Gardens—LIPI (Bali); Mount Ciremai National Park (TNGC) and Kuningan Botanic Gardens (Kuningan); and Baturraden Botanic Gardens (Baturraden). This study did not collect or sample any protected species or its propagules. The participant that appears in this manuscript has given written informed consent (as outlined in the PLOS consent form) to publish his image.

These botanic gardens (Table 1) contain native and exotic plant species in their living collections, and are all located next to native tropical montane rainforest ecosystems. The surrounding native forest contains various exotic species, some of which appear to have escaped from the gardens. These ecosystems are dominated by families’ characteristic of tropical montane rainforests, including Fagaceae, Euphorbiaceae, Meliaceae and Lauraceae [3032]. The climate at the four sampling locations is similar, except for Eka Karya which is relatively drier [32, 33].

thumbnail
Table 1. Information about the four tropical rainforest study sites in Indonesia, each located next to a botanic garden.

https://doi.org/10.1371/journal.pone.0202254.t001

Distance sampling

Line transect distance sampling surveys were conducted at all four sampling locations, from the border of each botanic garden towards the interior of the adjacent native rainforest, following the method of Buckland [23] (lay out in S2 Fig). Total transect line length across all sites was 3350 m. Different transect lengths were used at each location because of variations in rainforest sizes (constraining transect length), land contour (affecting survey feasibility), and detection rate of the surveyed exotics. We decided to stop the sampling on a particular transect line when detection rates of exotic species on that transect line were low. We considered the detection rate low when we did not detect the same species within 100 m of the last detection along that transect. Records were limited to detections within a 20 m wide strip (10 m perpendicular to either side of the transect lines). This corresponds to a total area of 6.7 ha (3350 m x 20 m = 67000 m2): 2.7 ha in Cibodas, 1.8 ha in Eka Karya, 0.9 ha in Kuningan, and 1.3 ha in Baturraden. Perpendicular distances were measured using a Bosch GLM50 laser rangefinder.

Leaf shape, leaf size and plant height

Our candidate species traits were leaf size, leaf shape and plant height, with methodology for measuring traits based on Pérez-Harguindeguy et al. [34]. With a limited number of exotic species recorded during the surveys, we limited the number of species traits to reduce the risk of over-fitting in the detection model. While other traits undoubtedly influence detection, we focused on these traits because we expected them to be particularly relevant to detection of exotics species in the understorey. We conducted multiple measurements for each trait (on each individual) and used average values following the suggestion of Pérez-Harguindeguy, Díaz [34]. One leaf sample set was collected from each of the detected plants. Each set consisted of between 1 and 10 leaves (mean ≈ 3–4), depending on leaf availability during sampling. Plant height was defined as “the shortest distance between the upper boundary of the main photosynthetic tissues (excluding inflorescences) on a plant and the ground level” [34]. Traits values may be subject to phenotypic plasticity [35]. However, we did not take account for this and assumed that the sampled trait values reflected a single value for each species.

Along with the exotic species data samples, we also collected digital photo images and data on physical leaf samples from dominant native species. Native species dominance was determined based on written descriptions of the vegetation communities (S1 File). The native species considered dominant were 4 tree species and several understorey species (3 vines, 7 shrubs, and 8 herbs).

We collected healthy, young but fully expanded leaves with no physical, predator, or disease caused defects [34]. We photographed samples next to a scale in a laboratory using a Panasonic Lumix DMC-TZ10 digital camera for leaf area, leaf size and leaf shape calculation (examples for every species in S4 Fig).

Leaf area and complexity data for exotic and native samples were measured by analysing the leaf image in the IMAGEJ image processing program [36]. The distribution of leaf area among species was quite skewed, so we quantified leaf size as the square-root of leaf area; this helped make the distribution of the data more even so the potential for leverage was reduced. The use of the square-root of leaf area, which approximates the average linear dimension of the leaves, is also a natural way to consider the size of a leaf. We quantified leaf complexity as the ‘leaf dissection index’ of McLellan and Endler [37], which is computed as leaf perimeter divided by the square-root of leaf area. High values of leaf complexity indicate more complex leaf shape or extreme width-length leaf ratio. If leaf complexity of exotics is considerably different from natives, then leaf complexity may improve exotic species detectability.

Detected exotic leaf sizes were generally similar to natives (Fig 1A); a Welch two-sample t-test of the difference between native and exotic leaf size was not significantly different from 0 (95% confidence interval: [-3.29 cm, 0.93 cm]; p-value = 0.266; sample sizes: exotics 25 species, native 22 species). Similarly, there were no significant differences between exotics and natives leaf shape (95% confidence interval of difference: [-2.41 cm, 0.31 cm]; p-value = 0.125, sample sizes: exotics 25, native 22) (Fig 1).

thumbnail
Fig 1. Leaf traits data distribution of exotic and dominant native species collected from four locations: Cibodas, Eka Karya Bali, Baturraden and Kuningan.

(A) Distribution of square-root leaf area of native and exotic species. (B) Leaf shape of native and exotic species. Total samples: 1890 samples (1338 samples from 25 exotic species + 552 samples from 22 native species).

https://doi.org/10.1371/journal.pone.0202254.g001

Detection model

We used a half-normal detection model to describe how the probability of detection of species i at study site j declines with d, the perpendicular distance from the transect: Eq 1

This model assumes that detection of the species on the transect line is perfect, and involves a scale parameter σ, which controls the width of the detection function. Large values of σ imply that the probability of detection remains high even if the plant is far away from the transect; small values imply that detectability decays quickly with distance from the transect.

In our analysis, we truncated the detection function (zeroed beyond 10 metres), to match the field sampling protocol (only detections within 10 meters from the transect were recorded).

We analysed our multi-species multi-site dataset in a single model, with the scale parameter described as a function of study site and species characteristics using a log regression as follows: Eq 2 where b0j controls the average detectability across all species at each site (each of the 4 gardens), each of the bx parameters is a regression coefficient that represents how the corresponding explanatory variable influences detectability, and ɛi is a random effect term that allows for un-modelled variation in detectability among species (i.e., variation not accommodated by the explanatory variables). Random effect values were modelled as being drawn from a normal distribution with mean zero and a standard deviation that was estimated from the data. The explanatory variables included were: mean leaf shape (S), mean height (H) and mean leaf size (square-root leaf area) (L) at the species level. Explanatory variables were centered and, to account for possible non-linear effects, we included quadratic terms of S (S2) and L (L2). In tropical rainforest ecosystems, native species have particular leaf sizes or shapes ranges due to their response to high rainfall [38]. Therefore, an exotic species with leaves that differ substantially in size or shape from those of natives (either smaller or larger) might be easier to detect, so a quadratic relationship might be expected.

We conducted our analysis in a Bayesian framework using JAGS [39], called from R with package jagsUI [40]. We provide our JAGS code in S2 File. We ran 3 MCMC chains and drew 10,000 samples per chain, omitting the first 5,000 iterations as a burn-in. We assessed convergence of the MCMC chains using the Gelman-Rubin R-hat diagnostic [41]; we assumed convergence when R-hat values were smaller than 1.1. We used uninformative priors for all the parameters. Specifically, for the regression coefficients we used a normal distribution with a mean of 0 and standard deviation equal to 1000 (i.e. precision 1 x 10−6). The prior for the standard deviation of the random effect was a uniform distribution on the range 0 to 10. The likelihood function of the distance sampling model [42] is not readily available as a sampling distribution in JAGS, so in our code we implemented it using the “zeroes trick” which allows arbitrary sampling distributions to be used [43]. We limited the analysis to species that were detected at least five times.

Results

We recorded 1440 individuals in total, detecting 55 exotic species of which 25 had at least five detections. The 30 species omitted for analyses made up only 102 detections. Cestrum aurantiacum Lindl. at Cibodas was the species most frequently detected with 399 detections (Fig 2).

thumbnail
Fig 2. The number of times each exotic species was detected at each of the four locations for the species included in the detectability model.

BBG: Bali Botanic Gardens, BRBG: Baturraden Botanic Gardens, CBG: Cibodas Botanic Gardens, KBG: Kuningan Botanic Gardens.

https://doi.org/10.1371/journal.pone.0202254.g002

We did not detect large exotic tree species in this study, but we detected small exotic trees (Calliandra callothyrsus and Zapoteca tetragona). In higher elevation tropical forest, trees tend to be smaller than at lower elevations [32]. We are also certain that there are no large exotic trees in the study site because previous studies also noted their absence [4448].

The probability of detection of individuals increased significantly with height; the regression coefficient bH was estimated with mean = 0.110 and 95% credible interval [0.028, 0.200] (Table 2), with the detection probability changing by up to 0.4 over the range of the data.

thumbnail
Table 2. Summary statistics for the posterior distribution of parameters in the detection model (Eq 1 and Eq 2).

https://doi.org/10.1371/journal.pone.0202254.t002

Detectability varied substantially among the four locations. Detectability was lowest at Kuningan and highest at Cibodas. The standard deviation of the random effect was not estimated to be particularly large, suggesting that there are modest differences in the detectability of species that remain unexplained by species height, leaf area and leaf shape (Fig 3).

thumbnail
Fig 3. Predicted detection (sigma, in log scale) of 25 detected exotic species.

Predicted detection based on fixed model only labelled with ‘f’ (e.g. sp1f). Predicted detection includes random effects in the detectability model labelled with ‘r’ (e.g. sp1r).

https://doi.org/10.1371/journal.pone.0202254.g003

We failed to detect an effect of leaf shape on detectability; the 95% credible intervals for the coefficient for leaf shape and squared leaf shape include zero (Table 2). However, we did find evidence of a linear effect of leaf size (bA), with detectability generally increasing with leaf size, particularly up to 11–13 cm; the effects were relatively small, with the estimated change in the probability of detection generally being < 0.2 over the range of the data (Fig 4). The credible interval of the quadratic term for leaf size (bA2) included zero, suggesting that the expected relationship between leaf size and detectability (increased detectability of species with extreme leaf sizes; see Detection model in the methods) might not exist. Indeed, the estimated effect suggests, if anything, the opposite relationship, with slightly greater detectability at intermediate leaf sizes, although any non-monotonic relationship that might exist appears to be small (Fig 4).

thumbnail
Fig 4. The predicted detectability of exotic species at a particular survey location.

(a) Predicted detectability as a function of plant height. (b) Predicted detectability as a function of leaf size. Predicted detectability at different distances (2 m, 5 m, and 10 m) from the transect represented in one curve. The shaded areas show 95% credible intervals.

https://doi.org/10.1371/journal.pone.0202254.g004

Discussion

As expected, taller species were detected more easily than shorter species (Fig 4A). A previous study that assessed the effect of plant life form on detectability did not find such an effect [25]. This can be explained because it was a study of detectability at the species level and detection probability of species varies with local abundance [49]. Larger species will tend to be less abundant than smaller species. Therefore, the limited variation in detection probability among life forms in Chen et al. [25] might be due to differences in abundance. In addition, we can explain the important effect of height in our study given that it was conducted in the lower layer of the tropical rainforest, which is dense and filled with numerous different plant species (Fig 5), particularly in secondary tropical rainforest [50]. Tropical rainforest is the most diverse and complex ecosystem type among all tropical land ecosystems [51]. Thus, considering the dense and crowded community in the lower layers of tropical secondary rainforest, it is not surprising that taller exotic species were easier to detect than shorter ones.

thumbnail
Fig 5. An example of lower layer rainforest (understorey).

This native rainforest is located next to Cibodas Botanic Gardens and managed by Gede-Pangrango National Park. Picture was taken in May 2015 by DIJ.

https://doi.org/10.1371/journal.pone.0202254.g005

The probability of detecting exotic plant species increased with leaf size. This is reasonable given that larger leaves might be more prominent, similar to the effect of plant height. Species with smaller leaves might be more easily obscured by other leaves and simply harder to recognize during surveys. The anticipated effect of species being more easily detected if they had extreme leaf sizes might not exist, perhaps because the range of leaf sizes of exotic species was quite similar to that of native species (Fig 1), so exotic species with extreme leaf sizes might not stand out as being distinctively exotic.

Initially, we expected that exotics with higher leaf shape scores would be easier to detect. However, we did not find evidence for this. There are several potential reasons why leaf shape might not influence detectability. The similar leaf shapes of natives and exotics mean that exotic species might not stand out from natives during surveys. More generally, it may be hard to detect a correlation between leaf shape and detectability because leaf shape varies among and within species [52]. However, similar results were obtained when leaf shape was introduced in the model as an individual level characteristic (unpublished data). Alternatively, other leaf shape measures from morphometric analysis may lead to different results than when using the leaf dissection index. Morphometric analysis might reveal leaf shapes that are more detectable during surveys. Leaf complexity may also be important when combined with leaf colour (i.e. both of them collectively affect detectability). Colour is an important plant attribute for detection purposes [12, 53], but we did not consider it in this study because insufficient data were available.

This study did not consider the spatial distribution of detected exotic species (e.g. randomly distributed vs clumpy distributed; [54]). Detected exotic species in this study might not be distributed randomly, a major assumption in the standard distance sampling method [23]. We also did not include environmental characteristics or surveyor experience, unlike Garrard et al. [12]. We controlled surveyor experience by having the same surveyor for all samplings. In our study, environmental characteristics were similar across sites and the same person conducted all the surveys. We did not consider topographic factors in this study. All four botanic gardens are situated in similar contexts (at the foot of a mountain and within relatively similar elevation ranges and topographic situations) and we assumed that topography did not strongly affect detectability across different study sites. Lastly, height and leaf size themselves might not be the actual reasons for improved detectability, but rather other important traits correlate with them (e.g. showy flowers or fruit color). We focused on leaf traits and height, ignoring flower or fruit as detection cues, because during the observation periods, most of the detected exotics were without flowers or fruits.

Distance sampling and plant abundance

Distance sampling is often used for animal sampling to estimate abundance [15, 16, 55], but has surprisingly been rarely used for either exotic or tropical rainforest species level vegetation analysis. Distance sampling has similar accuracy to standard survey methods [26], but can be simpler, easier to use and less labour-intensive than standard methods like quadrats or Whittaker’s plots [22]. These standard sampling techniques were primarily developed for use in temperate systems and mostly for multi-species detection [22]. However, standard plot techniques are less suitable and likely to be more laborious and time consuming for exotic plant detection in dense tropical rainforests, which have complex structures and profiles [50]. Our multi-species detection model might be usefully incorporated into distance sampling to help estimate abundance, particularly of rare (and potentially sessile) species.

Value of traits for exotic plant detectability model

A key value of this study is that it uses information on traits to help predict detectability. The contribution of leaf area and plant height to exotic species detectability are useful for invasive species management. Given trait-based detectability predictions, this study demonstrates that general visual traits (height and leaf size) can be important factors to consider for exotic species survey effort estimation and detection optimization, aside from sampling / survey designs [5658] and appropriate sampling time [9, 18].

Since most of the detected exotics in this study were shrubs, herbs, and small trees, leaf traits and height significantly contributed to the detectability of these exotic species. We only tested these traits because the detected exotic species were limited to particular life forms, and we aimed to keep the number of traits manageable to avoid over-fitting given the moderate number of species (three traits, 25 species). Other traits may influence detectability, particularly of those species in the higher level of tropical forests. Bark characteristics or branching patterns may be more useful when detecting taller trees. In general, traits that increase the visibility of the species, and traits that help distinguish the exotic species from the natives will be important. Further studies to separate these two factors in the detection model framework would be useful.

Trait-based studies have greatly increased the understanding of invasion processes [59] and for risk assessments [60] in tropical ecosystems, yet have rarely been used for informing exotic species detection in this bioregion. Trait-based detection studies might provide additive solutions for detection optimization. The methodology can be applied to different traits that may be relevant.

Conclusions

We conclude that vegetative traits such as plant height and leaf size affect exotic plant detection in tropical forests. Even though the effect of leaf size is relatively minor, this trait remains a relevant predictor of exotic species detection in tropical forests. Further study is needed to assess how other potentially relevant vegetative traits, such as branching patterns and wood-bark texture, influence detectability. Our study demonstrates the utility of traits when predicting exotic species detectability, crucial information when determining the required detection effort.

Supporting information

S1 Fig. Map of sampling locations: Cibodas, Kuningan, Baturraden and Eka Karya Bali.

https://doi.org/10.1371/journal.pone.0202254.s001

(TIF)

S2 Fig. Line transect distance sampling lay out.

Black circles represent individuals of exotic species and dashed arrows show their perpendicular distance to the transect. Surveys were conducted from the border of the botanic gardens towards the native rainforest interior. Only detections within 10 m of the transect line were recorded. The number of transects, distance between transects and their length varied among locations depending on the field conditions.

https://doi.org/10.1371/journal.pone.0202254.s002

(TIF)

S3 Fig. History plot of Markov chain Monte Carlo (MCMC) simulation based on the constructed detectability model.

The plots demonstrate the converged chains in the model simulation for all variables (nodes) included (b0, bS, bS2, bH, bA, and bA2) (left part of the graph). The credible interval for all involved variables are not containing zero values except for bA, bS and bS2, suggesting weak contribution detected from leaf size (bA), shape (bS) and its quadratic term (bS2) to detectability in this study (right part of the graph).

https://doi.org/10.1371/journal.pone.0202254.s003

(TIF)

S4 Fig. Examples of used leaf photos for leaf area and leaf shape calculation of 25 detected exotic species.

The photo ordered from smallest to the largest leaf shape (complexity) value.

https://doi.org/10.1371/journal.pone.0202254.s004

(TIF)

S1 File. List of reference publications used for native species composition information data.

https://doi.org/10.1371/journal.pone.0202254.s005

(PDF)

S2 File. R code script for Bayesian analysis of the tested detection model in this study.

https://doi.org/10.1371/journal.pone.0202254.s006

(PDF)

S3 File. Total data set used for the analysis.

https://doi.org/10.1371/journal.pone.0202254.s007

(CSV)

Acknowledgments

DIJ thanks the managements and numerous staff members of Cibodas Botanic Gardens—LIPI, Eka Karya Bali Botanic Gardens—LIPI, Gede-Pangrango National Park (TNGGP), BKSDA Bali, Ciremai National Park (TNGC), Baturraden Botanic Gardens (Kebun Raya Baturraden), and Kuningan Botanic Gardens (Kebun Raya Kuningan) for their supports during fieldworks. The authors also thank Stuart Jones for dedicated peer review and comments on the draft manuscript. DIJ thank two anonymous reviewers (thesis external reviewers) for their constructive comments and suggestions for the manuscript improvement. DIJ is supported by a Melbourne Research Scholarship (MRS) and Centre of Excellence for Biosecurity Risk Analysis (CEBRA) during the completion of this work. We acknowledge support from the Australian Research Council (DE120102221 to JAC) and ARC Centre of Excellence for Environmental Decisions (CEED).

References

  1. 1. Gurevitch J, Padilla DK. Are invasive species a major cause of extinctions? Trends Ecol Evol. 2004;19(9):470–4. pmid:16701309
  2. 2. Ehrenfeld JG. Ecosystem consequences of biological invasions. Annu Rev Ecol Evol Syst. 2010;41:59–80.
  3. 3. Olson LJ. The economics of terrestrial invasive species: a review of the literature. Agricultural and Resource Economics Review. 2006;35(1):178.
  4. 4. Dodd AJ, Ainsworth N, Burgman MA, McCarthy MA. Plant extirpation at the site scale: implications for eradication programmes. Diversity and Distributions. 2015;21(2):151–62.
  5. 5. McNeely JA. An introduction to human dimensions of invasive alien species. The Great Reshuffling: Human Dimensions of Invasive Alien Species IUCN Publishers, Gland, Switzerland. 2001:5–22.
  6. 6. Leung B, Springborn MR, Turner JA, Brockerhoff EG. Pathway-level risk analysis: the net present value of an invasive species policy in the US. Front Ecol Environ. 2014;12(5):273–9.
  7. 7. Pimentel D, Zuniga R, Morrison D. Update on the environmental and economic costs associated with alien-invasive species in the United States. Ecol Econ. 2005;52(3):273–88.
  8. 8. Peh KSH. Invasive species in Southeast Asia: the knowledge so far. Biodivers Conserv. 2010;19(4):1083–99.
  9. 9. Bogich TL, Liebhold AM, Shea K. To sample or eradicate? A cost minimization model for monitoring and managing an invasive species. J Appl Ecol. 2008;45(4):1134–42.
  10. 10. Collins English Dictionary–Complete and Unabridged, 12th Edition 2014. 2014. Detectability.
  11. 11. Chen G, Kéry M, Plattner M, Ma K, Gardner B. Imperfect detection is the rule rather than the exception in plant distribution studies. J Ecol. 2013;101(1):183–91.
  12. 12. Garrard GE, McCarthy MA, Williams NS, Bekessy SA, Wintle BA. A general model of detectability using species traits. Methods Ecol Evol. 2013;4(1):45–52.
  13. 13. Lahoz‐Monfort JJ, Guillera‐Arroita G, Wintle BA. Imperfect detection impacts the performance of species distribution models. Global ecology and biogeography. 2014;23(4):504–15.
  14. 14. Archaux F, Henry P-Y, Gimenez O. When can we ignore the problem of imperfect detection in comparative studies? Methods Ecol Evol. 2012;3(1):188–94.
  15. 15. Sollmann R, Gardner B, Williams KA, Gilbert AT, Veit RR. A hierarchical distance sampling model to estimate abundance and covariate associations of species and communities. Methods in Ecology and Evolution. 2015.
  16. 16. Sollmann R, Gardner B, Williams KA, Gilbert AT, Veit RR. A hierarchical distance sampling model to estimate abundance and covariate associations of species and communities. Methods Ecol Evol. 2016.
  17. 17. Olson CA, Beard , Karen H, Koons , David N, Pitt , William C. Detection probabilities of two introduced frogs in Hawaii: implications for assessing non-native species distributions. In: Center UNW, editor. Nebraska: Staff Publications; 2012.
  18. 18. Mehta SV, Haight RG, Homans FR, Polasky S, Venette RC. Optimal detection and control strategies for invasive species management. Ecol Econ. 2007;61(2):237–45.
  19. 19. Christy MT, Yackel Adams AA, Rodda GH, Savidge JA, Tyrrell CL. Modelling detection probabilities to evaluate management and control tools for an invasive species. J Appl Ecol. 2010;47(1):106–13.
  20. 20. Hauser CE, McCarthy MA. Streamlining ‘search and destroy’: cost‐effective surveillance for invasive species management. Ecol Lett. 2009;12(7):683–92. pmid:19453617
  21. 21. Bouxin G, Le Boulengé E. A phytosociological system based on multi-scaled pattern analysis: a first example. Vegetatio. 1983;54(1):3–16.
  22. 22. Ellenberg D, Mueller-Dombois D. Aims and methods of vegetation ecology. New York: John Wiley & Sons; 1974.
  23. 23. Buckland ST, Anderson D.R., Burnham K. P., Laake J. L., Borchers D. L., and Thomas L. Introduction to distance sampling estimating abundance of biological populations. New York: Oxford University Press; 2001.
  24. 24. Buckland ST, Borchers DL, Johnston A, Henrys PA, Marques TA. Line transect methods for plant surveys. Biometrics. 2007;63(4):989–98. pmid:18078477
  25. 25. Chen G, Kery M, Zhang J, Ma K. Factors affecting detection probability in plant distribution studies. J Ecol. 2009;97(6):1383–9.
  26. 26. Kissa D, Sheil D. Visual detection based distance sampling offers efficient density estimation for distinctive low abundance tropical forest tree species in complex terrain. For Ecol Manage. 2012;263:114–21.
  27. 27. Oo WP, Koike F. Dry forest community types and their predicted distribution based on a habitat model for the central dry zone of Myanmar. For Ecol Manage. 2015;358:108–21.
  28. 28. Amaral CH, Roberts DA, Almeida TIR, Souza Filho CR. Mapping invasive species and spectral mixture relationships with neotropical woody formations in southeastern Brazil. ISPRS Journal of Photogrammetry and Remote Sensing. 2015;108:80–93.
  29. 29. Pouteau R, Meyer J-Y, Stoll B. A SVM-based model for predicting distribution of the invasive tree Miconia calvescens in tropical rainforests. Ecol Modell. 2011;222(15):2631–41.
  30. 30. Werner WL. A comparison between two tropical montane ecosystems in Asia: Pidurutalagala (Ceylon/Sri Lanka) and Pangrango-Gede (Java). Mountain Research and Development. 1986:335–44.
  31. 31. Yamada I. Forest ecological studies of the montane forest of Mt Pangrango, West Java. 1. Stratification and floristic composition of the montane rain forest near Cibodas South East Asian Studies. 1975;13:402–26.
  32. 32. Hamzah A, Toha M, Van Steenis C. Mountain flora of Java. 1972.
  33. 33. Whitten T, Soeriaatmadja RE, Afiff SA. Ecology of Java & Bali: Oxford University Press; 1996.
  34. 34. Pérez-Harguindeguy N, Díaz S, Garnier E, Lavorel S, Poorter H, Jaureguiberry P, et al. New handbook for standardised measurement of plant functional traits worldwide. Australian Journal of Botany. 2013;61(3):167–234.
  35. 35. Turner KG, Fréville H, Rieseberg LH. Adaptive plasticity and niche expansion in an invasive thistle. Ecol Evol. 2015;5(15):3183–97. pmid:26357544
  36. 36. Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat Meth. 2012;9(7):671–5.
  37. 37. McLellan T, Endler JA. The relative success of some methods for measuring and describing the shape of complex objects. Syst Biol. 1998;47(2):264–81.
  38. 38. Givnish T. Leaf and canopy adaptations in tropical forests. Physiological ecology of plants of the wet tropics: Springer; 1984. p. 51–84.
  39. 39. Plummer M. JAGS: A program for analysis of Bayesian graphical models using Gibbs sampling. 2003.
  40. 40. Kellner K. A Wrapper Around 'rjags' to Streamline 'JAGS' Analyses. 2015.
  41. 41. Gelman A, Rubin DB. Inference from iterative simulation using multiple sequences. Statistical science. 1992:457–72.
  42. 42. Buckland ST, Rexstad EA, Marques TA, Oedekoven C. Distance Sampling: Methods and Applications: Springer; 2015.
  43. 43. Spiegelhalter D, Thomas A, Best N, Lunn D. WinBUGS user manual. version; 2003.
  44. 44. Junaedi DI. Exotic plants in the Cibodas Botanic Gardens remnant forest: inventory and cluster analysis of several environmental factors. Buletin Kebun Raya. 2014;17(1):1–8.
  45. 45. Mutaqien Z, Zuhri M. Establishing a long-term permanent plot in remnant forest of Cibodas Botanic Garden, West Java. Biodiversitas. 2011;12(4):218–24.
  46. 46. Rozak AH, Gunawan H. ALTITUDINAL GRADIENT AFFECTS ON TREES AND STAND ATTRIBUTES IN MOUNT CIREMAI NATIONAL PARK, WEST JAVA, INDONESIA. Jurnal Penelitian Kehutanan Wallacea. 2015;4(2):93–9.
  47. 47. Sutomo S, Undaharta N, Bangun T, Lugrayasa I. STUDI AWAL KOMPOSISI DAN DINAMIKA VEGETASI POHON HUTAN GUNUNG POHEN CAGAR ALAM BATUKAHU BALI. Bumi Lestari. 2012;12(2).
  48. 48. Zuhri M, Mutaqien Z. The Spread of Non-native Plant Species Collection of Cibodas Botanical Garden into Mt. Gede Pangrango National Park. Journal of Tropical Life Science. 2013;3(2):74–82.
  49. 49. McCarthy MA, Moore JL, Morris WK, Parris KM, Garrard GE, Vesk PA, et al. The influence of abundance on detectability. Oikos. 2013;122(5):717–26.
  50. 50. Whitmore TC. An Introduction to Tropical Rain Forests. Oxford: Oxford University Press; 1998. 251 p.
  51. 51. May RM. How many species are there on earth? Science(Washington). 1988;241(4872):1441–9.
  52. 52. Kidner CA, Umbreen S. Why is leaf shape so variable. International Journal of Plant Development. 2010;3:64–75.
  53. 53. Kendal D, Hauser CE, Garrard GE, Jellinek S, Giljohann KM, Moore JL. Quantifying plant colour and colour difference as perceived by humans using digital images. PLoS One. 2013;8(8):e72296. pmid:23977275
  54. 54. Condit R, Ashton PS, Baker P, Bunyavejchewin S, Gunatilleke S, Gunatilleke N, et al. Spatial patterns in the distribution of tropical tree species. Science. 2000;288(5470):1414–8. pmid:10827950
  55. 55. Isaac NJ, Cruickshanks KL, Weddle AM, Marcus Rowcliffe J, Brereton TM, Dennis RL, et al. Distance sampling and the challenge of monitoring butterfly populations. Methods Ecol Evol. 2011;2(6):585–94.
  56. 56. Huebner CD. Detection and monitoring of invasive exotic plants: a comparison of four sampling methods. Northeastern Naturalist. 2007;14(2):183–206.
  57. 57. Hui C, Foxcroft LC, Richardson DM, MacFadyen S. Defining optimal sampling effort for large‐scale monitoring of invasive alien plants: a Bayesian method for estimating abundance and distribution. J Appl Ecol. 2011;48(3):768–76.
  58. 58. Moore AL, McCarthy MA. Optimizing ecological survey effort over space and time. Methods Ecol Evol. 2016.
  59. 59. Dawson W, Burslem DFRP, Hulme PE. The comparative importance of species traits and introduction characteristics in tropical plant invasions. Diversity and Distributions. 2011;17(6):1111–21.
  60. 60. Chong KY, Corlett RT, Yeo DCJ, Tan HTW. Towards a global database of weed risk assessments: a test of transferability for the tropics. Biol Invasions. 2011;13(7):1571–7.