Beyond Global Models: Mapping the Spatially Contingent Relationship Between Soil Sand Content and Woody Invasion

Riparian ecosystems are being increasingly threatened by hydrological alteration and biological invasions, yet the role of local environmental heterogeneity in shaping invasion dynamics remains poorly understood. To address this, we tested the hypothesis that invasion patterns are spatially structured and therefore cannot be fully captured by global statistical models. We evaluated this hypothesis by analysing the relationship between soil sand content and the abundance of Gleditsia triacanthos in a riparian forest of the Esteros de Farrapos and Islands of the Uruguay River National Park, Uruguay. Generalized Linear Mixed Model revealed no significant relationship between soil sand content and G. triacanthos abundance (χ2 = 1.93, p = 0.17). In contrast, spatially explicit analyses showed that relationships between sand content and abundance were spatially contingent. Positive linear relationships predominated in areas with low sand content (mean 24.5%, n = 12), while negative relationships were restricted to the highest sand levels (mean 87.6%, n = 3). Intermediate sand-content zones (mean 47%, n = 16) showed no consistent patterns. These results suggest that invasion patterns vary across spatial contexts and may reflect the influence of different processes operating locally, indicating that relying solely on global analyses risks misinterpreting drivers and overlooking fine-scale variation. Our findings emphasize that understanding invasive species in heterogeneous systems requires considering whether mechanisms operate at local or broad scales, and that explicitly analyzing spatial structure can guide both hypothesis formulation and field study design.