Sampling bias and the robustness of ecological metrics for plant-damage-type association networks

Plants and their insect herbivores have been a dominant component of the terrestrial ecological landscape for the past 410 million years and feature intricate evolutionary patterns and co-dependencies. A complex systems perspective allows for both detailed resolution of these evolutionary relationships as well as comparison and synthesis across systems. Using proxy data of insect herbivore damage (denoted by the damage type or DT) preserved on fossil leaves, functional bipartite network representations provide insights into how plant-insect associations depend on geological time, paleogeographical space, and environmental variables such as temperature and precipitation. However, the metrics measured from such networks are prone to sampling bias. Such sensitivity is of special concern for plant-DT association networks in paleontological settings where sampling effort is often severely limited. Here, we explore the sensitivity of functional bipartite network metrics to sampling intensity and identify sampling thresholds above which metrics appear robust to sampling effort. Across a broad range of sampling efforts, we find network metrics to be less affected by sampling bias and/or sample size than richness metrics, which are routinely used in studies of fossil plant-DT interactions. These results provide reassurance that cross-comparisons of plant-DT networks offer insights into network structure and function and support their widespread use in paleoecology. Moreover, these findings suggest novel opportunities for using plant-DT networks in neontological terrestrial ecology to understand functional aspects of insect herbivory across geological time, environmental perturbations, and geographic space.