Quantitative assessment of plant-arthropod interactions in forest canopies: A plot-based approach
Martin Volf aff001; Petr Klimeš aff001; Greg P. A. Lamarre aff001; Conor M. Redmond aff001; Carlo L. Seifert aff001; Tomokazu Abe aff004; John Auga aff005; Kristina Anderson-Teixeira aff006; Yves Basset aff001; Saul Beckett aff007; Philip T. Butterill aff001; Pavel Drozd aff009; Erika Gonzalez-Akre aff006; Ondřej Kaman aff001; Naoto Kamata aff010; Benita Laird-Hopkins aff001; Martin Libra aff001; Markus Manumbor aff005; Scott E. Miller aff012; Kenneth Molem aff005; Ondřej Mottl aff001; Masashi Murakami aff004; Tatsuro Nakaji aff013; Nichola S. Plowman aff001; Petr Pyszko aff009; Martin Šigut aff009; Jan Šipoš aff014; Robert Tropek aff001; George D. Weiblen aff017; Vojtech Novotny aff001
Authors place of work:
Biology Centre of the Czech Academy of Sciences, Ceske Budejovice, Czech Republic
aff001; German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig, Leipzig, Germany
aff002; Faculty of Science, University of South Bohemia, Ceske Budejovice, Czech Republic
aff003; Faculty of Science, Chiba University, Chiba, Japan
aff004; New Guinea Binatang Research Center, Madang, Papua New Guinea
aff005; Conservation Ecology Center, Smithsonian Conservation Biology Institute; Front Royal, VA, United States of America
aff006; ForestGEO, Smithsonian Tropical Research Institute, Panama City, Panama
aff007; Maestria de Entomologia, Universidad de Panama, Panama City, Panama
aff008; Faculty of Science, University of Ostrava, Ostrava, Czech Republic
aff009; Graduate School of Agricultural and Life Sciences, The University of Tokyo, Furano, Japan
aff010; School of Biological Sciences, University of Bristol, Bristol, United Kingdom
aff011; National Museum of Natural History, Smithsonian Institution, Washington, DC, United States of America
aff012; Tomakomai Experimental Forest, Hokkaido University, Tomakomai, Japan
aff013; Institute of Botany, Czech Academy of Sciences, Brno, Czech Republic
aff014; Department of Zoology, Fisheries, Hydrobiology and Apiculture, Mendel University in Brno, Brno, Czech Republic
aff015; Department of Ecology, Faculty of Science, Charles University, Prague, Czech Republic
aff016; Bell Museum and Department of Plant & Microbial Biology, University of Minnesota, Saint Paul, MN, United States of America
Published in the journal:
PLoS ONE 14(10)
Research on canopy arthropods has progressed from species inventories to the study of their interactions and networks, enhancing our understanding of how hyper-diverse communities are maintained. Previous studies often focused on sampling individual tree species, individual trees or their parts. We argue that such selective sampling is not ideal when analyzing interaction network structure, and may lead to erroneous conclusions. We developed practical and reproducible sampling guidelines for the plot-based analysis of arthropod interaction networks in forest canopies. Our sampling protocol focused on insect herbivores (leaf-chewing insect larvae, miners and gallers) and non-flying invertebrate predators (spiders and ants). We quantitatively sampled the focal arthropods from felled trees, or from trees accessed by canopy cranes or cherry pickers in 53 0.1 ha forest plots in five biogeographic regions, comprising 6,280 trees in total. All three methods required a similar sampling effort and provided good foliage accessibility. Furthermore, we compared interaction networks derived from plot-based data to interaction networks derived from simulated non-plot-based data focusing either on common tree species or a representative selection of tree families. All types of non-plot-based data showed highly biased network structure towards higher connectance, higher web asymmetry, and higher nestedness temperature when compared with plot-based data. Furthermore, some types of non-plot-based data showed biased diversity of the associated herbivore species and specificity of their interactions. Plot-based sampling thus appears to be the most rigorous approach for reconstructing realistic, quantitative plant-arthropod interaction networks that are comparable across sites and regions. Studies of plant interactions have greatly benefited from a plot-based approach and we argue that studies of arthropod interactions would benefit in the same way. We conclude that plot-based studies on canopy arthropods would yield important insights into the processes of interaction network assembly and dynamics, which could be maximised via a coordinated network of plot-based study sites.
Ants – Arthropoda – Forests – Cherries – Interaction networks – Leaves – Trees – Herbivory
Forest canopies represent one of the most diverse environments on the planet , harbouring a large proportion of terrestrial arthropod diversity estimated at 6.8 million species . At the same time, canopies are among the least explored habitats due to the logistical challenges of accessibility . This combination of high diversity and inaccessibility has fascinated biologists for more than 150 years .
The development of single-rope climbing and fogging has provided ecologists with efficient tools for researching canopy arthropod communities, generating several influential studies e.g. [3, 4, 5]. Such studies spurred the development of new methods of access that nowadays include canopy walkways, canopy rafts, balloons, cherry pickers, or canopy cranes . Canopy studies have contributed to our understanding of species global diversity and biotic interactions [2, 6], but, as pointed out by Lowman et al. , “…the real challenge is ahead. Canopy organisms, both mobile and sessile, must be surveyed and their roles measured.”
Research into canopy arthropods has progressed from species inventories to the study of their interactions, allowing us to understand how hyper-diverse communities of canopy arthropods are maintained . Particular sampling methods are suitable for different systems and questions concerning the various roles arthropods play in forest canopies  (Table 1). Methods that allow access to individual branches or certain parts of the canopy are suitable for exploratory studies on arthropod diversity, detailed surveys focused on specific taxa, or manipulative experiments e.g. [8, 9]. They also allow for comparative studies across various canopy microhabitats and their arthropod communities . However, to fully census interactions between arthropods and plants on the level of the canopy as a whole, sampling methods must provide access to the entire canopy, from the terminal branches, through the inner canopy, to the lower branches. This is because arthropod species composition may differ considerably among various parts of the canopy , reflecting variation in resource availability and leaf traits . Neglecting some parts of the canopy, therefore, has the potential to influence the results of the census. In addition, methods suitable for censusing canopy arthropod interactions must facilitate the sampling of arthropods in such a way that enables the reliable reconstruction of the interaction network. In the tropics, transient herbivorous arthropod species (i.e. species with no lasting association to the sampled plant) can comprise up to 20% of species found on a particular tree . Thus, dead arthropods sampled from a plant do not constitute reliable interactions. To reliably reconstruct interaction networks, one needs to either sample live arthropods for feeding trials  or a use molecular detection of trophic interactions . Similarly, it is necessary to map ant nests rather than simply sample individual ants, as up to half of the ants foraging in a tree are tourists from surrounding vegetation .
Most importantly, for a quantitative analysis of arthropod interaction networks, the methods should allow structured sampling across large parts of the canopy, thus including all species in proportion to their abundance [7, 21, 28]. Previous studies often focused on sampling individual tree species, individual trees or selected constituent parts. Selective sampling is particularly beneficial for exploring insect-plant interactions in a phylogenetical or evolutionary framework as it allows the researcher to focus on particular lineages of interest [29, 30]. Methods that employ selective sampling are also valuable when assessing herbivore specialization or the effects of host-plant traits on insect community structure. This is because all focal species can be sampled with equal effort, thus allowing for direct comparisons between herbivore or host species [14, 31]. However, a drawback of selective sampling is that it does not facilitate quantitative network structure analyses, because it tends to skew interaction frequencies, over- or underestimate specialization and diversity, and biases network structure . In particular, it typically omits a high proportion of the arthropod and plant taxa co-existing at the sites, hence not reflecting species diversity and network structure at the whole forest level. We argue that for interaction network analyses, a plot-based approach, where entire plots are censused for plants and arthropods, is preferable, as it more accurately reflects the diversity and abundance of the available resources [21, 23].
Plot-based approaches applied to forest vegetation have greatly benefitted plant ecology research . We anticipate the study of arthropod interaction networks would benefit in equal measure . We accessed canopies using tree felling, canopy crane, and cherry picker techniques (Fig 1) across biogeographic regions (Palearctic, Nearctic, Neotropical, and Australian) and forest types (tropical vs. temperate, lowland vs. montane, primary vs. secondary). We compare our plot-based methods with non-plot-based sampling and highlight the strengths and limitations of the methods for sampling mobile flightless exophytic herbivores (leaf-chewing insect larvae), endophytic herbivores (miners and gallers), and flightless invertebrate predators (spiders and ants). Our aim is to stimulate plot-based research by providing practical and reproducible sampling guidelines for the analysis of arthropod interaction networks in forest canopies. We expect i) plot-based data and non-plot-based data to provide largely different estimates of interaction network structure as non-plot-based sampling skews frequencies between rare and abundant species, ii) felling to be the most efficient method in terms of sampling effort as it allows employing large teams of field workers who can simultaneously access large parts of the canopy, iii) all three methods to provide similar acessibility to the canopy with access to over 75% of foliage.
Materials and methods
During our operations, we took advantage of ongoing logging operations (Mikulcice, Toms Brook) and shifting agriculture (PNG sites); no plot was cleared solely for sampling. All projects were conducted in close collaboration with the local community and land owners. We obtained all research and export permits where required. Arthropods and plants from Papua New Guinea were sampled and exported under the permits nr. 070382, 070384, 080275, 010075, 011209, 011324, 012134, 014282, 0133004, 133005, and 018060 issued by Department of Environment and Conservation, Papua New Guinea, and 0139/2008, 0162/2010, and 0203/2013 issued by Forest Research Institute and Department of Forests, Papua New Guinea. Arthropods and plants from Panama were obtained and exported under the permits nr. SE/A-49-16, SE/AP-28-16, SC/AP-2-16, SEX/P-30-17, SEX/A-67-17, SEX/A-76-17 issued by Ministerio de Ambiente, Panama. The individuals whose faces are fully or partially visible in Fig 1 have given written informed consent (as outlined in PLOS consent form) to publish these photos.
Following a standardized protocol (Appendix 1) and workflow (Fig 2), we sampled i) lowland temperate forests in the Czech Republic (Mikulcice, Lanzhot), Japan (Tomakomai), and USA (Toms Brook); ii) lowland tropical forests in Panama (San Lorenzo) and Papua New Guinea (hereafter PNG; Wanang); and iii) highland tropical forests in PNG (Numba, Yawan) (Table 2, S1 Table).
Setting up the plot
At each location, we selected 0.1 ha plots with a vegetation structure and species composition typical for local broadleaf forests (Table 2, S1 and S2 Tables). In Wanang, Numba, and Yawan, plots were larger and subdivided into 0.1 ha sections (Table 2). Forest edges, plantations, stands with non-native vegetation, and large gaps were all avoided, as were steep slopes and swampy areas (for technical and safety reasons). We took GPS coordinates of all plot corners and used measuring tape or laser range finders to set up the plot and map all plants with DBH ≥5 cm. Each stem was tagged and identified to species level. It took 2–12 hours for three people to set up a 0.1 ha plot and map 24–251 trees within. In Tomakomai (4 trees), Toms Brook (8), and San Lorenzo (17) some trees proved to be hazardous to sample or were damaged by factors beyond our control, such as a hurricane, during the sampling. These trees were replaced by conspecifics or other broadleaf trees with a similar DBH adjacent to the plot. One non-native and one coniferous tree in Toms Brook were treated in the same way.
Timing the sampling
Arthropod abundances and species composition can vary dramatically throughout the year in seasonal forests. For example, temperate leaf-chewing insects exhibit one major peak during spring leaf-flush, and a smaller peak in late summer . Furthermore, peaks in abundance may differ among arthropod guilds, for instance leaf miners, where the major peak appears to occur later than for leaf-chewers (S1 Fig, S3 Table). A single, short sampling campaign can fail to capture all arthropod groups. Therefore, we generally sampled temperate plots at a slower pace throughout the season to mitigate this undesirable effect, returning periodically to the sites in order to sample trees (i.e. typically one to several tree individuals were sampled per day depending on their canopy sizes). Sampling effort was increased during abundance peaks if they materialised. During this period, sampling was conducted whenever the weather permitted. In this way, the variation in sampling effort mirrors the variability in insect abundance, and the probability that an insect will be sampled remains constant throughout the season. We spread the sampling seasonally within each sampled tree species to avoid a bias due to an unbalanced seasonal sampling (Appendix 1). In wet tropical forests, sampling was carried out with constant effort throughout the seasons as the effects of seasonality are much less pronounced and individual species appear throughout the year . However, a variable sampling strategy would be advisable in dry tropical and subtropical forests, where seasonality asserts greater influence . Such intense and relatively long-term sampling required careful logistical planning. This involved negotiating the research plan well in advance with land owners, crane drivers, chainsaw operators, and local managers so as to avoid clashing with other projects at the given sites. For example, for the sampling from felled trees, we specifically sought plots of forest that were scheduled for logging and paid the loggers to cut the trees on our schedule.
The requirements for accessing the forest canopy and obtaining live arthropods dramatically limit the range of suitable methods for the study of quantitative arthropod interaction networks (Table 1). We sampled arthropods from felled trees, and from standing trees using canopy cranes or cherry pickers (Fig 1). Arthropods were, as far as possible, completely sampled from all trees with DBH ≥5 cm. The percentage of the canopy accessed was visually estimated for each tree (Appendix 1). We sampled on days without strong rain or wind to mitigate safety risks and lowered arthropod activity due to harsh weather. The focal arthropod groups included all live leaf-chewing insect larvae (free feeding and semi-concealed), leaf mines, galls (insects and mites), spiders, and ants (foraging and nesting; S1 Table). Some species of galls were extremely abundant, making their complete sampling impractical. In such cases, we selected 3–5 branches each with 100–500 leaves, calculated the mean number of galls per leaf per branch, and used the resulting values to estimate the total abundance on the respective tree (Appendix 1).
Felling trees as a standardized destructive method is only suitable when it does not contribute to net deforestation. During our operations, we took advantage of ongoing logging operations (Mikulcice, Toms Brook) and shifting agriculture (PNG sites); no plot was cleared solely for sampling. All projects were conducted in close collaboration with the local community and land owners.
Sampling began with the clearing of the understory, followed by the felling of trees with DBH ≥5cm. One tree was felled at a time, starting with the shortest. Lianas on trees were cut prior to felling in order to free up the focal tree from its neighbours. Felled individuals were directed into gaps created by previous felling. Once felled, the entire tree (trunk included) was searched and all focal arthropods hand collected, a process taking anywhere from minutes to several hours, depending on the crown size. Prompt work minimized the loss of arthropods through dispersal or predation. It also prevented the contamination by foraging ants and spiders from the ground. Using division of labour, each team member focused primarily on one arthropod group, but would also contribute to the collection of secondary groups. Trees were always fully sampled on the day of felling, and necessitated teams of 7–21 members, dependent on study site and season (Table 2).
Unlike sampling from cranes and cherry pickers, felling allows the sampling of arthropods dwelling in large branches and trunks, such as nesting ants (Table 1). At felling sites, we intensively searched every tree for ant nests and foraging ants with a team of two to three collectors, as described in Klimes et al. . Foraging ants were collected first, before searching for ant nests by cutting branches, inspecting live and dead twigs, by dissecting parts of the trunk and bark, and by inspection of epiphytic aerial soil (Appendix 1).
Conversely, felling is not suitable for mobile, flying herbivores . Even non-flying herbivores may become dislodged when the crown forcefully impacts the ground. If this were a serious concern, the ratio between endophytic herbivores and leaf-chewing larvae would depend on the method. However, the ratio of leaf-chewing larvae to active miners sampled in individual 0.1 ha plots did not differ among the methods (χ2 (2) = 2.57, p = 0.2764) when compared by linear mixed-effect models using the ‘lmer4’ R package , with site as a random effect.
We sampled arthropods from canopy cranes in Tomakomai and San Lorenzo. In Tomakomai, the crane is 25 m high, covers ca 0.5 ha of forest, and is operated by researchers from the gondola. In San Lorenzo, the crane is operated by a driver. The maximum accessible height from the gondola is 40.5 m. The crane covers almost 1.0 ha of tropical forest .
There were 4–7 team members working in the field, typically including 2 members sampling from the crane (canopy team), 1–2 members sorting samples on the ground (ground team), and possibly 1–2 members accessing larger mid-story trees by climbing (climbing team). The canopy team sampled branches starting at the tip and working towards the base, in order to minimize arthropod loss during sampling. Arthropods were sampled by beating onto a beating tray, followed by a visual search and hand collection of any remaining arthropods. The canopy team was assisted by an additional member during periods of peak arthropod abundance. The samples were regularly delivered to the ground team for sorting.
Sampling from the crane was augmented with other methods. The canopy team accessed understory trees from ladders. Step ladders were ideal for sampling 3–5 m tall trees. For sampling at heights up to 8 m, or on sloped terrain, modular ladder poles were more efficient and stable. In addition, more complex forest architecture, as in San Lorenzo, required the climbing team. Using a single rope technique, they accessed those mid-story trees inaccessible from the gondola or ladders (Fig 1).
A cherry picker (an elevated truck-mounted work platform) was employed in Lanzhot. The 20 ton vehicle was transported by truck to the site, thus necessitating a forest access road. We used a Platform GENIE Z-135/70 JRT (Genie Industries, Redmond, WA, USA), which was equipped with a retractable arm enabling canopy access up to 43 m. The arm was operated by researchers directly from the basket. This four-wheel drive model can operate on gravel or clay forest roads, but not on off-road terrain. Plots were set up along a forest road with a firm dirt surface (~4 m wide, and completely covered by forest canopy) in order to provide good access to the plot from a single straight trajectory and to avoid having to manoeuvre the cherry picker between trees. Two team members sampled trees starting from the base and working towards the treetop. Arthropods were sampled using a beating tray combined with hand collection of any remaining individuals, before a final manual search by both workers. Samples were delivered to the ground team for processing before transportation to the laboratory. There were 2–6 people processing samples in the ground team, depending on insect abundance.
Leaf area estimates
We calculated the leaf area of sampled trees in order to standardize arthropod abundance and allow cross-site comparisons (Appendix 1).
At the felled sites, we quantified leaf biomass directly by defoliating each tree and weighing the fresh foliage. Mature and young leaves were sampled and weighed separately immediately following herbivore sampling. Care was taken that only leaves, with no other plant parts such as twigs and flowers, were sampled. At Mikulcice and Toms Brook sites, where team size was limited, only 50% or 25% of the canopy was defoliated on the largest trees and the results extrapolated to 100%. This measure was taken to ensure the complete sampling of large trees on the day of felling.
At the crane and cherry picker sites, defoliating trees and weighting the biomass was not possible. Instead, we visually estimated the number of young and mature leaves on standing trees. These estimates were conducted separately for every branch sampled for arthropods. The estimates were carried out for branches with ca 500 leaves each by two persons from the canopy team. The mean value of the two estimates was taken. The branch level estimates for the given tree were then summed to give an estimate for the entire tree. This method yielded more reliable results than if estimating leaves on larger branches or whole trees.
At all sites, a random sample of leaves from each tree was then arranged on a 50 x 50 cm board with white background (the “leaf frame”) and photographed. One frame each of young and mature leaves was processed for small trees (DBH <15 cm), while at least two frames were processed for larger trees. The leaf area of each sample was then calculated using ImageJ 1.48 . For felled trees, we included the weight of the sample to obtain the area to weight ratio. For the trees sampled from cranes and cherry pickers, we divided the leaf area of the sample by the number of leaves in the frame to obtain the mean area per leaf.
Finally, we calculated the total sampled leaf area for each tree using (i) the total leaf biomass and the area to weight ratio from the photographed sample for the felled trees, or (ii) the estimated total number of leaves on the tree multiplied by the mean leaf size of the photographed sample for the crane and cherry picker trees.
In Tomakomai, Mikulcice, and Lanzhot, pre-sorting, photographing, and labelling of samples was done in the field by a team consisting of 1–6 members, depending on arthropod abundance (Appendix 1). This made subsequent sorting in the lab much faster. Smaller trees in Toms Brook were treated the same way. Otherwise, samples were processed entirely in the laboratory.
We assigned all leaf-chewing insect larvae, galls, and mines to morphospecies according to their morphology . Each morphospecies was given a unique code name and was photographed. We preferred to assign initial morphotypes de novo per each individual tree sampled instead of using a complex system of morphospecies across all trees within the plot or even across multiple plots (Appendix 1). This approach is rapid and resistant to errors as even incorrect morphotyping does not generate false host plant records. It requires a second step where individual morphospecies are cross-referenced across all trees on completion of sampling. It is suitable for taxonomically poorly known and species diverse samples, where per-guild richness for an entire plot could reach hundreds of morphospecies.
We reared larval insect herbivores to adults or parasitoids (Appendix 1). Only in Toms Brook, where insect taxonomy and host associations are well known, were leaf-chewing larvae immediately stored in ethanol due to the overwhelming logistics of rearing all. We preserved larvae that died during rearing, the larvae from Toms Brook, spiders, and representative samples of all ant castes from each nest or foraging event in vials with 95% ethanol for subsequent DNA barcoding. The results of DNA barcoding along with reared adults are being used to refine morphospecies concepts and assign final identifications [21, 23, 39–41]. See Data Accessibility section for details on the publicly available sequences.
Non-plot-based studies typically focus on i) abundant tree species or ii) a taxonomically/phylogenetically representative selection of species e.g. [42, 43, 44]. In order to compare plot-based and non-plot-based methods, we derived both types of data from plot-based data on plant-caterpillar interactions in 0.8 ha of PNG highland primary rainforest Yawan . This dataset was chosen because it is species-rich, the caterpillar species exhibit various levels of host specificity, and were already identified to an acceptable level. Only living trees identified to species and caterpillar species/morphospecies with confirmed host associations were included. The pruned dataset, representing 0.8 ha of primary forest (eight x 0.1 ha plots), included 113 tree species and 186 caterpillar morphospecies.
We computed network statistics and structure from the Yawan primary plots, and compared them with networks comprised of i) the most abundant (in terms of amount of foliage based on leaf area calculations) tree species and ii) a taxonomically representative selection of primary tree species including all tree families that had at least 200 m2 of foliage sampled for arthropods.
For type (i) networks, we combined all 0.1 ha plots to represent a larger patch (0.8 ha) of rainforest and ranked tree species in order of decreasing amount of total foliage. We then selected the species whose cumulative total foliage represented 20%, 40%, 60%, and 80% of the total foliage of the 0.8 ha patch (3, 7, 15, and 31 species, respectively). In each threshold category, we rarefied the foliage amount of each species (Fsp) to equal the average total foliage of a 0.1 ha plot divided by the respective number of tree species (Fthresh). This was achieved by randomly selecting individual trees until Fsp > = Fthresh. The final trees (Tf) are only partially sampled of their caterpillars (Cf) so that Fsp = Fthresh (if we use 0.25Tf then we take 0.25Cf randomly selected caterpillars, rounded to the nearest integer). The individual trees and the caterpillars found on them made up the networks from which statistics were computed. The process was repeated 100 times for each category to account for the random tree (and caterpillar) selection.
For type (ii) networks (taxonomical selection), we limited the dataset to tree species and families that had at least 200 m2 of foliage sampled (200 m2 equals ca 0.7% of the total foliage in the 0.8 ha primary forest patch). We selected the most abundant species per family in terms of leaf area (19 species and families selected). With this tree selection, rarefaction then proceeded as per type (i) networks.
We focused our comparison between plot-based and non-plot-based networks on: i) connectance (realised proportion of possible links), ii) web asymmetry (balance between numbers of species in the two levels; positive values indicate higher proportion of higher trophic level species), iii) nestedness (temperature of the matrix; 0 means high nestedness, 100 means chaos), iv) species richness of caterpillars, v) weighted generality (mean effective number of host species per caterpillar species), and vi) weighted vulnerability (mean effective number of caterpillar species per host species) as defined and computed in ‘bipartite’ package . The network parameters were compared using 95% confidence intervals.
Furthermore, we compared the efficiency of each sampling method we used across the plots. We expressed the method efficiency as i) Foliage accessibility per plot (the average percentage of accessible foliage), ii) Area-based sampling effort (ASE) required to sample each 0.1 ha plot (in person-hours), and iii) Resource-based sampling effort (RSE) required to sample 1 m2 of foliage (in person-hours). Only time spent on sample collection and sorting in the field was counted towards the sampling effort. Workers helping with logistics (chainsaw operators or the crane driver in Panama) were excluded. We modelled the relationship between these components of sampling efficiency and the sampling Method (felling, crane, cherry picker), Forest type (temperate, tropical lowland primary, tropical lowland secondary, tropical highland primary, tropical highland secondary), the Number of stems (DBH ≥ 5 cm) per plot, and Sampled leaf area using linear mixed-effect models as implemented in the R package ‘lmer4’ . Foliage accessibility was arcsine-transformed and sampling effort log-transformed. We used Site as a random factor in all mixed-effect models (S2 Table). Model simplification by forward selection was employed to produce the most parsimonious model based on Akaike’s Information Criterion (AIC). All analyses were performed in R software version 3.4.0 .
In total, we sampled focal arthropod groups from 5.3 ha of forest, representing 6,280 trees and 167,744 m2 of foliage (Table 2). We sampled 89,243 leaf-chewing larvae, 14,547 active mines, 135,446 abandoned mines, 28,698 spiders, 35,343 ant individuals, 3,487 ant nests, and sampled or estimated abundance of 2,963,942 insect and mite galls (S2 Table).
All non-plot-based data types showed highly biased network structure towards higher connectance, higher web asymmetry, and higher nestedness temperature when compared to the plot-based data (Fig 3). Non-plot-based data using the most abundant tree species representing 20% of the foliage in the local forest had the highest caterpillar richness while those using a taxonomically representative selection of tree families had the lowest caterpillar richness. Non-plot-based data using the most abundant tree species representing 20% and 40% of the foliage in the local forest showed lower generality than the plot-based data. Differences in vulnerability were less pronounced mainly because of the high variability in plot-based data. However, vulnerability was highest in the data using the most abundant tree species representing 20% and 40% of the foliage in the local forest.
On average, Foliage accessibility was 82.5% ±3.9% (mean ±SD) foliage in felled plots, 82.7% ±3.3% foliage in plots sampled by canopy crane, and 89.3% ±6.3 foliage in plots sampled by cherry picker (S2 Fig). Foliage accessibility correlated with Method (χ2 (2) = 6.91, p = 0.0254). The optimum model, after simplification, included the fixed effects Forest type (highest in lowland and highland secondary tropical forests), Method (highest from the cherry picker), Number of stems (positive correlation), and Sampled leaf area (negative correlation) (χ2 (8) = 64.02, p < 0.0001) (S4 Table).
The average ASE required to sample a 0.1ha plot was 1583 ±579 person-hours (mean± SD) for felled trees, 1867 ±673 for sampling by canopy crane, and 1128 ±305 for sampling by cherry picker. Method did not have a significant effect on ASE (χ2 (2) = 1.49, p = 0.4740). The optimum model that explained differences in ASE included the fixed effects Number of stems (positive correlation) and Forest type (highest in lowland primary tropical forests) (χ2 (5) = 95.24, p < 0.0001; S4 Table).
The average RSE to sample 1 m2 of foliage was 0.51 ± 0.24 (mean± SD) person-hours for sampling felled trees, 1.14 ±0.15 for sampling by canopy crane, and 0.92 ±0.10 for sampling by cherry picker. Method did not have a significant effect on RSE (χ2 (2) = 3.52, p = 0.1722). The optimum model explaining differences in RSE included the fixed effects Number of stems (positive correlation), Sampled leaf area (negative correlation), and Forest type (highest in temperate forests) (χ2 (6) = 80.75, p < 0.0001; S4 Table).
We propose a plot-based approach to studying arthropod interaction networks, using three methods for sampling a continuous area of forest canopy. Plot-based standardisation means that frequent associations can be distinguished from those that are casual or rare [48, 49]. Focusing on a selection of abundant tree species or representative families sampled with a standardized sampling effort skews the proportions between rare and common interactions. As expected, this resulted in higher connectance, high web asymmetry and higher nestedness in the simulated non-plot-based data. This is not surprising as all these parameters are linked to the network size, which has been reduced under the selective sampling scenario. In addition to reducing network size, non-plot-based sampling focused on abundant or phylogenetically distinct hosts can affect the patterns recovered in host specificity and diversity. This is because such hosts typically harbor distinct arthropod communities. Locally abundant tree species tend to harbor higher diversity of herbivores than rarer hosts . Focusing on such hosts can lead to over-estimations of diversity. On the other hand, hosts from isolated or chemically distinct families can have relatively species poor herbivore communities . Emphasizing such hosts in the data can lead to under-estimations of diversity. Also, many herbivores are shared between congeneric or confamilial hosts while the amount of shared herbivores decreases with the host phylogenetic distance . Specificity of interactions can thus be biased in datasets that include skewed proportions of such hosts although this trend was not particularly pronounced in our simulated data.
Plot-based sampling provides a robust description of the community structure as one can assume that the interactions are completely censused for the proportion of the canopy successfully sampled (Fig 4). One can then test and improve the performance of models that predict trophic interactions in real communities by decomposing the effects of abundance, plant characteristics and arthropod community composition [28, 41]. Derived food-web metrics are comparable on a common area basis, and may identify processes shaping communities of canopy arthropods across various habitats, ecosystems, or geographic regions .
Plot-based analyses thus represent an ideal counterpart to those based on a stratified selection of focal species sampled with an equal sampling effort. Methods that employ equal sampling effort are advantageous for studying the host specialization of herbivores, the effects of host traits on herbivore communities, or insect-plant interactions in a phylogenetic framework e.g. [14, 30]. However, modern methods enable the measurement of host specialization with respect to host phylogeny or chemical similarity in plot-based data also . Furthermore, a plot-based approach can be used to investigate spatial distribution of arthropods across the forest canopy and their impact on competitors and other trophic levels. This is important, for instance, when considering competition among ants where canopy connectivity and structure play important roles in forming ant communities . Furthermore, herbivores may have density-dependent effects on plant survival that need to be studied in a spatially explicit framework .
One limitation of plot-based sampling methods is that the logistical challenges necessitate relatively large teams and overall effort. Despite our expectations, however, all methods demanded comparably high sampling effort, with none being significantly more efficient. Such prerequisites stem from the need to census all parts of the canopy, including those difficult to access, in order to reconstruct truly quantitative interaction networks [7, 21]. Foliage accessibility positively correlated with the number of stems in the plot, probably because many of the trees in densely vegetated plots were small and easier to access. On the other hand, the number of stems within a plot increased both types of sampling effort that we quantified. ASE (total effort per 0.1ha plot) was highest in lowland primary tropical forests characterized by relatively high stem density and large trees difficult to sample. RSE (effort per 1 m2 of foliage) was highest in temperate forests. This may be because arthropod density is generally higher in temperate forests , especially during the spring abundance peak.
High effort per site prevented a rigorous methodological comparison where the same forest is sampled by all three methods. All methods enabled access to over 80% of the foliage. But the unbalanced distribution of methods may be one reason why the cherry picker appeared to provide better access to the canopy than felling or cranes. Similarly to Corff and Marquis , we operated the cherry picker in almost optimal conditions in temperate forest where plots were close to an access road and the trees could be accessed from a straight trajectory. Operating in less favourable conditions would dramatically decrease foliage accessibility or require employing additional methods. Sampling from cranes also had to be supplemented by other techniques at both our crane sites. While sampling by other techniques represented a small proportion of sampling effort in the temperate Tomakomai forest, it considerably increased the sampling effort in San Lorenzo tropical rain forest. In San Lorenzo, only 49% of the trees (representing 58% of the foliage sampled) were accessed solely by crane.
Each method also has its own specifics unrelated to its overall efficiency. Felling generally requires larger teams [22, 23] as felled trees need to be sampled immediately. Cranes and cherry pickers allow proceeding at a slower pace with a smaller team e.g. [12, 20]. The three methods are also not completely comparable in terms of the sampled arthropod groups. All were suitable for sampling endophytic and exophytic non-flying arthropods. Less mobile flying herbivores, such as aphids or psyllids, were also well represented in our samples, although they were not the focus of our study. Felling was the only method which enabled sampling of nesting ants, which can represent an important proportion of the canopy arthropods . Quantitative sampling of highly mobile macroscopic arthropods (adult beetles, flies or true bugs) was not possible by these methods, although they were better represented in crane and cherry picker samples.
Other methods, such as fogging, may be more suitable for surveying highly mobile arthropods . Such methods can also strongly reduce the required team size and effort. To assess trophic interactions, however, they would need to be combined with a massive barcoding effort so the sampled arthropods could be reliably assigned to their host-plants. A molecular approach to assessing trophic interactions is becoming increasingly popular  and can be especially useful in well studied or less diverse communities. However, the approach may face identification limitations in diverse communities that include a high proportion of closely related and/or hybridizing hosts. Indeed, standard barcode markers may fail to provide a sufficient resolution for such hosts unless combined with specifically selected ones . The implementation of such methods for large scale plot-based sampling should, therefore, be carefully considered. Furthermore, the sampling of endophytic or semi-concealed herbivores and ant nests would require the employment of additional methods.
We suggest that a global network using the methods described for area-based sampling would provide important insights into the processes of food web assembly and dynamics . To that end, we propose a network of permanent plots where canopy arthropods and their interactions would be censused by non-destructive sampling. The network of permanent plots could benefit from collaboration with the global network of ForestGEO plots  which generates major insights into forest community ecology. We suggest that plots of 0.1 ha are an appropriate size to be sampled from cranes or cherry pickers, and which allow for repeat surveys, while keeping the required effort manageable. A single 0.1 ha plot census can yield information on more than 100,000 canopy arthropods and their interactions, thus the potential to make significant contributions to arthropod ecology research is huge.
The network of permanent plots should ideally be augmented by a larger network of temporal plots to be sampled by felling. Despite a slight revival in canopy crane construction , such platforms are still missing from vast regions, including Africa and North America. Similarly, opportunities for the use of cherry pickers remain limited in many forests. The sampling of 0.1 ha plots by felling thus seems to be the only widely applicable option in many regions. These plots could be highly replicated and ideally adjacent to the ForestGeo plots.
Sampling canopy arthropods by felling can become a salvage sampling strategy to obtain data on arthropod communities being lost due to ongoing deforestation. There has been considerable activity in the past decade focused on constructing large-scale experiments, such as planting forest stands of a given richness , or manipulation of landscape fragmentation , which deepen our understanding of how ongoing changes in forest structure affect ecological interactions. However, ecologists have been slow to take advantage of ongoing logging operations, urban development, or shifting agriculture for destructive arthropod and plant sampling to salvage the data. Yet, such data in combination with data from permanent plots would enable the exploration of trends in arthropod networks along major environmental gradients . Furthermore, the detailed data obtained by our methods could be used for modelling forest composition and arthropod interactions. Combining such models with high-throughput methods, such as remote sensing, that allow forest composition to be assessed may enable us to predict the basic characteristics of plant-arthropod interactions over large spatial scales . Ultimately, the application of the outlined methods could lead to high-impact results with far-reaching consequences, such as the prediction of the effects of forest degradation on forest arthropod communities, and the identification and preservation of arthropod diversity hotspots in the world’s forests.
S1 Fig [docx]
Seasonal trends in abundance of leaf chewing larvae and active miners.
S4 Table [docx]
Variables with a significant effects on Foliage accessibility, Area-based sampling effort, and Resource-based sampling effort as selected by forward selection in linear mixed effect models.
S5 Table [docx]
List of staff, interns, students, volunteers, and local assistants who helped with the sampling.
1. Lowman MD, Schowalter TD, Franklin J. Methods in Forest Canopy Research: Univ of California Press; 2012.
2. Hamilton AJ, Novotny V, Waters EK, Basset Y, Benke KK, Grimbacher PS, et al. Estimating global arthropod species richness: refining probabilistic models using probability bounds analysis. Oecologia. 2013;171(2):357–65. doi: 10.1007/s00442-012-2434-5 22968292
3. Erwin TL. Tropical forests: their richness in Coleoptera and other arthropod species. The Coleopterists Bulletin. 1982;36(1):74–5.
4. Lowman M. An assessment of techniques for measuring herbivory: is rainforest defoliation more intense than we thought? Biotropica. 1984;16(4):264–8.
5. Basset Y, Aberlenc HP, Delvare G. Abundance and stratification of foliage arthropods in a lowland rain forest of Cameroon. Ecological Entomology. 1992;17(4):310–8.
6. Novotny V, Miller SE, Hrcek J, Baje L, Basset Y, Lewis OT, et al. Insects on plants: Explaining the paradox of low diversity within specialist herbivore guilds. American Naturalist. 2012;179(3):351–62. doi: 10.1086/664082 22322223
7. Godfray HCJ, Lewis OT, Memmott J. Studying insect diversity in the tropics. Philosophical Transactions of the Royal Society of London Series B-Biological Sciences. 1999;354(1391):1811–24.
8. Lowman M, Foster R, Wittman P, Rinker H. Herbivory and insect loads on epiphytes, vines and host trees in the rain forest canopy of French Guiana. In: Hallé F, editor. Biologie d’une canopée de forêt equatoriale–III Rapport de la mission d’exploration scientifique de la Canopée de Guyane. Paris, Lyon: Pro-Natura International & Operation Canopee; 1996. p. 116–128.
9. Díaz I, Sieving K, Peña-Foxon M, Armesto J. A field experiment links forest structure and biodiversity: epiphytes enhance canopy invertebrates in Chilean forests. Ecosphere. 2012;3(1):1–17.
10. Wardhaugh CW, Stork NE, Edwards W. Canopy invertebrate community composition on rainforest trees: Different microhabitats support very different invertebrate communities. Austral Ecology. 2014;39(4):367–77.
11. Basset Y, Hammond PM, Barrios H, Holloway JD, Miller SE. Vertical stratification of arthropod assemblages. In: Basset Y, Novotny V, Miller SE, Kitching R, editors. Arthropods of Tropical Forests: Spatio-Temporal Dynamics and Resource Use in the Canopy. Cambridge: Cambridge University Press; 2003. p. 17–27.
12. Murakami M, Yoshida K, Hara H, Toda MJ. Spatio-temporal variation in Lepidopteran larval assemblages associated with oak, Quercus crispula: the importance of leaf quality. Ecological Entomology. 2005;30(5):521–31.
13. Basset Y. Species abundance and body size relationships in insect herbivores associated with New Guinea forest trees, with particular reference to insect host-specificity. In: Stork NE, Adis J, Didham RK, editors. Canopy Arthropods. London, New York: Chapman & Hall; 1997. p. 237–64.
14. Novotny V, Basset Y. Host specificity of insect herbivores in tropical forests. Proceedings of the Royal Society of London B: Biological Sciences. 2005;272(1568):1083–90.
15. Kress WJ, García-Robledo C, Uriarte M, Erickson DL. DNA barcodes for ecology, evolution, and conservation. Trends in Ecology & Evolution. 2015;30(1):25–35.
16. Klimes P, Fibich P, Idigel C, Rimandai M. Disentangling the diversity of arboreal ant communities in tropical forest trees. PLOS One. 2015;10(2):e0117853. doi: 10.1371/journal.pone.0117853 25714831
17. Basset Y, Cizek L, Cuenoud P, Didham RK, Guilhaumon F, Missa O, et al. Arthropod Diversity in a Tropical Forest. Science. 2012;338(6113):1481–4. doi: 10.1126/science.1226727 23239740
18. Ødegaard F. Species richness of phytophagous beetles in the tropical tree Brosimum utile (Moraceae): the effects of sampling strategy and the problem of tourists. Ecological Entomology. 2004;29(1):76–88.
19. Wardhaugh CW, Stork NE, Edwards W, Grimbacher PS. The overlooked biodiversity of flower-visiting invertebrates. PLOS One. 2012;7(9):e45796. doi: 10.1371/journal.pone.0045796 23029246
20. Corff JL, Marquis RJ. Differences between understorey and canopy in herbivore community composition and leaf quality for two oak species in Missouri. Ecological Entomology. 1999;24(1):46–58.
21. Volf M, Pyszko P, Abe T, Libra M, Kotásková N, Šigut M, et al. Phylogenetic composition of host plant communities drives plant‐herbivore food web structure. Journal of Animal Ecology. 2017;86(3):556–65. doi: 10.1111/1365-2656.12646 28146344
22. Whitfeld TJS, Novotny V, Miller SE, Hrcek J, Klimes P, Weiblen GD. Predicting tropical insect herbivore abundance from host plant traits and phylogeny. Ecology. 2012;93(8):S211–S22.
23. Redmond CM, Auga J, Gewa B, Segar ST, Miller SE, Molem K, et al. High specialization and limited structural change in plant‐herbivore networks along a successional chronosequence in tropical montane forest. Ecography. 2019;42(1):162–72.
24. Reynolds BC, Crossley D Jr. Spatial variation in herbivory by forest canopy arthropods along an elevation gradient. Environmental Entomology. 1997;26(6):1232–9.
25. Kitching R, Bergelson J, Lowman M, McIntyre S, Carruthers G. The biodiversity of arthropods from Australian rainforest canopies: general introduction, methods, sites and ordinal results. Austral Ecology. 1993;18(2):181–91.
26. Lowman MD. Leaf growth dynamics and herbivory in five species of Australian rain-forest canopy trees. Journal of Ecology. 1992;80(3):433–47.
27. Schowalter TD, Zhang Y. Canopy arthropod assemblages in four overstory and three understory plant species in a mixed-conifer old-growth forest in California. Forest Science. 2005;51(3):233–42.
28. Klimes P, Idigel C, Rimandai M, Fayle TM, Janda M, Weiblen GD, et al. Why are there more arboreal ant species in primary than in secondary tropical forests? Journal of Animal Ecology. 2012;81(5):1103–12. doi: 10.1111/j.1365-2656.2012.02002.x 22642689
29. Volf M, Julkunen‐Tiitto R, Hrcek J, Novotny V. Insect herbivores drive the loss of unique chemical defense in willows. Entomologia Experimentalis et Applicata. 2015;156(1):88–98.
30. Volf M, Segar ST, Miller SE, Isua B, Sisol M, Aubona G, et al. Community structure of insect herbivores is driven by conservatism, escalation and divergence of defensive traits in Ficus. Ecology Letters. 2018;21:83–92. doi: 10.1111/ele.12875 29143434
31. Volf M, Hrcek J, Julkunen‐Tiitto R, Novotny V. To each its own: differential response of specialist and generalist herbivores to plant defence in willows. Journal of Animal Ecology. 2015;84(4):1123–32. doi: 10.1111/1365-2656.12349 25649252
32. Anderson‐Teixeira KJ, Davies SJ, Bennett AC, Gonzalez‐Akre EB, Muller‐Landau HC, Joseph Wright S, et al. CTFS‐ForestGEO: a worldwide network monitoring forests in an era of global change. Global Change Biology. 2015;21(2):528–49. doi: 10.1111/gcb.12712 25258024
33. Novotny V, Miller SE. Mapping and understanding the diversity of insects in the tropics: past achievements and future directions. Austral Entomology. 2014;53(3):259–67.
34. Novotny V, Basset Y. Seasonality of sap-sucking insects (Auchenorrhyncha, Hemiptera) feeding on Ficus (Moraceae) in a lowland rain forest in New Guinea. Oecologia. 1998;115(4):514–22. doi: 10.1007/s004420050549 28308272
35. Ribeiro DB, Freitas AV. Large-sized insects show stronger seasonality than small-sized ones: a case study of fruit-feeding butterflies. Biological Journal of the Linnean Society. 2011;104(4):820–7.
36. Basset Y, Charles E, Novotny V. Insect herbivores on parent trees and conspecific seedlings in a Guyana rain forest. Selbyana. 1999;20(1):146–58.
37. Basset Y, Horlyck V, Wright SJ. Studying forest canopies from above: the International Canopy Crane Network. Panama City: Smithsonian Tropical Research Institute and UNEP; 2003.
39. Miller SE, Hrcek J, Novotny V, Weiblen GD, Hebert PD. DNA barcodes of caterpillars (Lepidoptera) from Papua New Guinea. Proceedings of the Entomological Society of Washington. 2013;115(1):107–10.
40. Miller SE, Rosati ME, Gewa B, Novotny V, Weiblen GD, Hebert PD. DNA Barcodes of Lepidoptera Reared from Yawan, Papua New Guinea. Proceedings of the Entomological Society of Washington. 2015;117(2):247–51.
41. Segar ST, Volf M, Isua B, Sisol M, Redmond CM, Rosati ME, et al. Variably hungry caterpillars: predictive models and foliar chemistry suggest how to eat a rainforest. Proceedings of the Royal Society of London B: Biological Sciences. 2017;284(1866):20171803.
42. Novotny V, Miller SE, Baje L, Balagawi S, Basset Y, Cizek L, et al. Guild-specific patterns of species richness and host specialization in plant-herbivore food webs from a tropical forest. Journal of Animal Ecology. 2010;79(6):1193–203. doi: 10.1111/j.1365-2656.2010.01728.x 20673235
43. Novotny V, Drozd P, Miller SE, Kulfan M, Janda M, Basset Y, et al. Why are there so many species of herbivorous insects in tropical rainforests? Science. 2006;313(5790):1115–8. doi: 10.1126/science.1129237 16840659
44. Novotny V, Miller SE, Leps J, Basset Y, Bito D, Janda M, et al. No tree an island: the plant-caterpillar food web of a secondary rain forest in New Guinea. Ecology Letters. 2004;7(11):1090–100.
45. Dormann CF, Gruber B, Fründ J. Introducing the bipartite package: analysing ecological networks. Interaction. 2008;1:0.2413793.
46. Bates D, Maechler M, Bolker B, Walker S. lme4: Linear mixed-effects models using Eigen and S4. R package version. 2014;1(7):1–23.
47. R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing. https://www.R-project.org; 2017.
48. Lewis OT, Memmott J, Lasalle J, Lyal CH, Whitefoord C, Godfray HCJ. Structure of a diverse tropical forest insect–parasitoid community. Journal of Animal Ecology. 2002;71(5):855–73.
49. Tylianakis JM, Tscharntke T, Lewis OT. Habitat modification alters the structure of tropical host–parasitoid food webs. Nature. 2007;445(7124):202–5. doi: 10.1038/nature05429 17215842
50. Straw N, Ludlow A. Small-scale dynamics and insect diversity on plants. Oikos. 1994;71(1):188–92.
51. Grandez-Rios JM, Bergamini LL, de Araújo WS, Villalobos F, Almeida-Neto M. The effect of host-plant phylogenetic isolation on species richness, composition and specialization of insect herbivores: A comparison between native and exotic hosts. PLOS One. 2015;10(9):e0138031. doi: 10.1371/journal.pone.0138031 26379159
52. Jorge LR, Prado PI, Almeida‐Neto M, Lewinsohn TM. An integrated framework to improve the concept of resource specialisation. Ecology Letters. 2014;17(11):1341–50. doi: 10.1111/ele.12347 25168335
53. Bagchi R, Gallery RE, Gripenberg S, Gurr SJ, Narayan L, Addis CE, et al. Pathogens and insect herbivores drive rainforest plant diversity and composition. Nature. 2014;506(7486):85–8. doi: 10.1038/nature12911 24463522
54. Li HQ, Chen JY, Wang S, Xiong SZ. Evaluation of six candidate DNA barcoding loci in Ficus (Moraceae) of China. Molecular ecology resources. 2012;12(5):783–90. doi: 10.1111/j.1755-0998.2012.03147.x 22537273
55. Nakamura A, Kitching RL, Cao M, Creedy TJ, Fayle TM, Freiberg M, et al. Forests and their canopies: achievements and horizons in canopy science. Trends in Ecology & Evolution. 2017;32(6):438–51.
56. Grossman JJ, Vanhellemont M, Barsoum N, Bauhus J, Bruelheide H, Castagneyrol B, et al. Synthesis and future research directions linking tree diversity to growth, survival, and damage in a global network of tree diversity experiments. Environmental and Experimental Botany. 2018;152:68–89.
57. Ewers RM, Didham RK, Fahrig L, Ferraz G, Hector A, Holt RD, et al. A large-scale forest fragmentation experiment: the Stability of Altered Forest Ecosystems Project. Philosophical Transactions of the Royal Society of London B: Biological Sciences. 2011;366(1582):3292–302. doi: 10.1098/rstb.2011.0049 22006969
58. Müller J, Brandl R. Assessing biodiversity by remote sensing in mountainous terrain: the potential of LiDAR to predict forest beetle assemblages. Journal of Applied Ecology. 2009;46(4):897–905.