Integrating theory and experiments to link local mechanisms and ecosystem-level consequences of vegetation patterns in drylands Ricardo Martinez-Garcia1;∗, Ciro Cabal2, Justin M. Calabrese3;4;5, Emilio Hernández-García6, Corina E. Tarnita2, Cristóbal López6;y, Juan A. Bonachela7;y 1 ICTP South American Institute for Fundamental Research & Instituto de Física Teórica - Universidade Estadual Paulista, São Paulo, SP, Brazil. 2 Department of Ecology and Evolutionary Biology, Princeton University, Princeton, NJ, USA. 3 Center for Advanced Systems Understanding (CASUS), Görlitz, Germany. 4 Helmholtz-Zentrum Dresden Rossendorf (HZDR), Dresden, Germany. 5 Department of Ecological Modelling, Helmholtz Centre for Environmental Research – UFZ, Leipzig, Germany. 6 IFISC, Instituto de Física Interdisciplinar y Sistemas Complejos (CSIC-UIB), Palma de Mallorca, Spain. 7 Department of Ecology, Evolution, and Natural Resources, Rutgers University, New Brunswick, NJ, USA ∗ Correspondence: [email protected] y equal contribution. Abstract Self-organized spatial patterns of vegetation are frequent in water-limited regions and have been suggested as important indicators of ecosystem health. However, the mechanisms underlying their emergence remain unclear. Some theories hypothesize that patterns could result from a water- mediated scale-dependent feedback (SDF), whereby interactions favoring plant growth dominate at short distances and growth-inhibitory interactions dominate in the long range. However, we know little about how net plant-to-plant interactions may shift from positive to negative as a function of inter-individual distance, and in the absence of strong empirical support, the relevance of this SDF for vegetation pattern formation remains disputed. These theories predict a sequential change in pattern shape from gapped to labyrinthine to spotted spatial patterns as precipitation declines. Nonetheless, alternative theories show that the same sequence of patterns could emerge even if net interactions between plants were always inhibitory (purely competitive feedbacks, PCF). Importantly, although these alternative hypotheses lead to visually indistinguishable patterns they predict very different desertification dynamics following the spotted pattern, limiting their potential use as an ecosystem-state indicator. Moreover, vegetation interaction with other ecosystem components can introduce additional spatio-temporal scales that reshape both the patterns and the desertification dynamics. Therefore, to make reliable ecological predictions for a focal ecosystem, it is crucial that models accurately capture the mechanisms at play in the system of interest. Here, we review existing theories for vegetation self-organization and their conflicting predictions about desertification dynamics. We further discuss possible ways for reconciling these predictions and potential empirical tests via manipulative experiments to improve our understanding of how vegetation self-organizes and better predict the fate of the ecosystems where they form. arXiv:2101.07049v3 [q-bio.PE] 15 May 2021 Keywords: ecological patterns, competition, scale-dependent feedback, ecological transitions, spatial self-organization, mathematical models. 1 Introduction From microbial colonies to ecosystems extending over continental scales, complex biological systems often feature self-organized patterns, regular structures that cover large portions of the system and emerge from nonlinear interactions among its components (Meinhardt, 1982; Camazine et al., 2003; Sole and Bascompte, 2006; Pringle and Tarnita, 2017). Importantly, because harsh environmental conditions provide a context in which self-organization becomes important for survival, emergent patterns contain 1 crucial information about the physical and biological processes that occur in the systems in which they form (Sole and Bascompte, 2006; Meron, 2018; Zhao et al., 2019). A well-known example of ecological self-organization is vegetation pattern formation in water-limited 4/27/2021 Bing Mapas: indicaciones, planificación de viajes, cámaras de tráfico y mucho más ecosystems (Deblauwe et al., 2008; Rietkerk and van de Koppel, 2008). Flat landscapes can show veg- etationNotas spots regularly distributed on a matrix of bare soil, soil-vegetation labyrinths, and gaps of bare soil regularly interspersed on a homogeneous layer of vegetation (Fig. 1). Importantly, water availability strongly influences the specific shape of the pattern. In agreement with model predictions (von Hard- enberg et al., 2001; Meron et al., 2004),4/27/2021 a Fourier-basedBing analysis Mapas: indicaciones, planificación of satellite de viajes, cámaras de tráfico imagery y mucho más covering extensive areas of Sudan revealed that more humid regions are dominated by gapped patterns, whereas spotted patterns dominate in more arid conditions12.27, (Deblauwe 3.18 et al., 2011). However, imagery time series are not long enough to observe whether vegetation cover in a specific region undergoes these transitions between patterns when aridity increases over time. a b 4/27/2021 Bing Mapas: indicaciones, planificación de viajes, cámaras de tráfico y mucho más4/27/2021 Bing Mapas: indicaciones, planificación de viajes, cámaras de tráfico y mucho más 11.02, 28.18 11.02, 28.18 250 pies 100 m © 2021 Maxar c d 250 pies 100 m © 2021 Maxar https://www.bing.com/maps 1/1 200m https://www.bing.com/maps 1/1 250 pies 100 m 500 pies 200 m © 2021 Maxar © 2021 Maxar Figure 1: Aerial images of self-organized vegetation patterns. In all panels, vegetated regions are darker and bare-soil regions, lighter. a) Spot pattern in Sudan (11◦34’55.2” N; 27◦54’47.52”E). b) Labyrinthine pattern in Mali (12◦27’50”N; 3◦18’30”E). c) Gap pattern in Niger (11◦1’12”N; 28◦10’48”E). d) Band pattern in Sudan (11◦3’0”N; 28◦17’24”E). Microsoft product screen shot(s) reprinted with permission from Microsoft Corporation. Image Copyrights ©2021 Maxar. Following the spotted pattern, if precipitation continues to decrease, models predict that patterned ecosystems undergo a transition to a desert state. This observed correlation between pattern shape and water availability suggests that the spotted pattern could serve as a reliable and easy-to-identify early- https://wwwwarning.bing.com/maps indicator of this ecosystem shifthttps://www (Scheffer.bing.com/maps 1/1 and Carpenter, 2003; Rietkerk et al., 2004;1/1 Scheffer et al., 2009; Dakos et al., 2011, 2015). This has reinforced the motivation to develop several models aiming to explain both the formation of spatial patterns of vegetation and their dependence on environmental variables (von Hardenberg et al., 2001; Rietkerk et al., 2002; Meron et al., 2004; Borgogno et al., 2009; Martínez-García et al., 2014; Gowda et al., 2014). Although some attempts to test model predictions with satellite imagery exist (Weissmann et al., 2017; Bastiaansen et al., 2018), theoretical studies using models remain the dominant approach to study this hypothesized transition. All these frameworks successfully reproduce the sequence of gapped, labyrinthine, and spotted pat- terns found in satellite imagery (Fig. 2a). However, they disagree in their predictions regarding the nature of the desertification transition that follows the spotted pattern. For example, Rietkerk et al. (2002) and von Hardenberg et al. (2001) predicted that ecosystems may undergo abrupt desertification, including a hysteresis loop, following the spotted pattern. Martínez-García et al. (2013a) and Yizhaq and Bel (2016) predicted that desertification could also occur gradually with progressive loss of vegetation biomass. Using alternative modeling approaches, other studies have supported the idea that whether an ecosystem will collapse gradually or abruptly is system-dependent and determined by the intensity of 2 stochasticity (Weissmann et al., 2017), vegetation and soil type (Kéfi et al., 2007b), colonization rates (Corrado et al., 2015), and intensity of external stresses, such as grazing (Kéfi et al., 2007a). Because drylands cover ∼ 40% of Earth’s land surface and are home to ∼ 35% of the world population (Mortimore et al., 2009), determining whether they will respond abruptly or gradually to aridification is critical both from an ecosystem-management and socio-economic point of view. Active lines of theoretical research aiming to address this question have focused on understanding how different components of the ecosystem may interact with each other to determine an ecosystem’s response to aridification (Bonachela et al., 2015; Yizhaq and Bel, 2016), as well as on designing synthetic feedbacks (in the form of artificial microbiomes) that could prevent ecosystems collapses or make such transitions smoother (Villa Martín et al., 2015; Conde-Pueyo et al., 2020; Vidiella et al., 2020). The question has also attracted considerable attention from empirical researchers (Maestre et al., 2016), and recent evidence suggests that certain structural and functional ecosystem attributes respond abruptly to aridity (Berdugo et al., 2020). Despite current efforts, whether desertification is more likely to occur gradually or abruptly remains largely unknown. Here, we outline and rank strategies that will help answer this question. In section 2, we discuss the ecological rationale behind existing models for vegetation self-organization. We review such models in section 3 and summarize their opposing predictions about the ecosystem collapse in section 4. In section 5, we describe possible manipulative experiments and empirical measures to