R_JAGS code for estimation and analysis of species-area-relationship (SAR) parameters from NEON (National Ecological Observatory Network) data on plant surveys

<p><span>Invasive species science is heavily geared toward the invasive agent. </span>However, management to protect native species also requires a proactive approach focused on understanding the features affecting community vulnerability to invasion impacts<span>. </span><span>Vulnerability </span><span>is likely the result of </span><span>factors acting across spatial scales, from </span><span>local to regional, and it is the combined effects of these factors that will determine the magnitude of vulnerability.</span><span> We introduce an analytical framework that quantifies the scale-dependent impact of biological invasions from the shape of the native species-area-relationship (SAR). We leverage newly available, biogeographically extensive vegetation data from the US National Ecological Observatory Network to assess plant community vulnerability to invasion impact as a function of factors acting across scales. We analyzed more than 1000 SARs widely distributed across the USA along environmental gradients and under different levels of invasion. </span>Results show that a decrease in native richness is consistently associated with invasive species cover<span>, but it is only at relatively high levels of invasion that native richness is compromised. After accounting for variation in baseline ecosystem diversity, net primary productivity, and human modification, ecoregions that are colder and wetter seem to be most vulnerable to losses of native plant species at the local level, while warmer and wetter areas seem most susceptible at the landscape level. We also document how the combined effects of cross-scale factors result in a heterogenous spatial pattern of vulnerability. </span><span>This pattern </span><span>cannot be predicted by analyses at any single scale, underscoring the importance of accounting for factors acting across scales. Simultaneously assessing differences in vulnerability between distinct plant communities at local, landscape and regional scales provided outputs that can be used to inform policy and management aimed at reducing vulnerability to the impact of plant invasions.</span></p>