Gonzalez-Garcia, A.; Neyret, M.; Lopez-Tejedor, A.; Prima, M. C.; Si-Moussi, S.; Renaud, J.; Gueguen, M.; Lavorel, S.

Protected areas cannot halt biodiversity loss in isolation; integrating them with surrounding human-dominated landscapes is critical. However, this integration is challenged by substantial landscape heterogeneity at their borders, hindering our understanding of cross-border changes in ecosystem service provision. We introduce a novel framework for characterizing these dynamics by analyzing ecosystem service gradients along protected area borders. For 16 protected areas in the French Alps, we assessed 12 ecosystem services using a mix of established biophysical models and novel connectivity-based models for mobile species. These were aggregated into three stakeholder-driven domains reflecting respectively rural, cultural, and urban management priorities. Automated polynomial regression analysis classified borders into five gradient types. The most common were 'Decreasing Gradients', representing a decline in ecosystem services outside the protected area, and 'Increasing Gradients', with the opposite pattern. Our analysis reveals these patterns are driven by specific landscape configurations, uncovering frequent trade-offs between the three management priorities, where, for instance, landscapes supporting rural priorities often degrade cultural and urban ones. We also identify key opportunities for synergies, by identifying areas where ecosystem services for all three priority domains increase simultaneously outside the protected area. This spatially explicit typology provides a powerful diagnostic tool for designing targeted interventions, such as prioritizing habitat restoration where ecosystem services decline or managing agricultural landscapes to mitigate conflicts across management priorities, supporting a more effective integration of protected areas into the wider landscape.