Being Well-Connected Pays in a Disturbed World: Enhanced Herbivory in Better-Linked Habitats

Seascapes are typically comprised of multiple components that are functionally linked by the movement of organisms and fluxes of matter. Changes to the number and spatial arrangement of these linkages affect biological connectivity that, in turn, can alter ecological functions. Herbivory is one such function, pivotal in controlling excessive algal growth when systems become disturbed. Here, we used microcosm experiments to test how the change to connectivity affects herbivory under different levels of disturbance. We applied network theory to measure types of connectivity at different scales (patch and whole system) and quantified herbivory by a crustacean mesograzer exposed to excess algae, mimicking pulse and press disturbances. We demonstrate that greater connectivity significantly enhances herbivory in Clibanarius virescens: Both the number of linkages and their spatial arrangement interact to shape the response of herbivory in systems to disturbance. Our findings highlight the value of controlled experiments for advancing theories about the potential effects of connectivity on important ecological functions, such as herbivory, and justify further investigation to measure how connectivity might affect the resilience of ecosystems. We posit that the variation in the type, and scale, of spatial linkages might have profound consequences for managing the capacity of ecosystems to respond to disturbance.


Introduction
The movements of organisms, nutrients and energy connect populations and ecosystems [1,2]. Connectivity can influence the composition and diversity of assemblages, especially those responding to disturbance [3][4][5]. As a result, connectivity can also affect the distribution and efficacy of important ecological functions, such as herbivory, predation, scavenging and the dispersal of seeds, larvae and pollen [6][7][8][9][10][11][12]. Spatial linkages can also alter how ecosystems respond to disturbance by modifying key functions that underpin the capacity for resilience [13][14][15]. For example, the spatial configuration of wildlife corridors on land, and marine reserve networks in the sea, can influence dispersal patterns, modify the diversity of responses by individuals contributing to the same ecosystem function and determine the potential for ecosystems to recover from disturbance [6,16,17].
Herbivory is an important process in aquatic systems that confers resilience in the face of a disturbance by helping ecosystems to resist shifts to alternative states [3,18,19]. The relationship between herbivory and coral reef resilience is particularly well-understood; examples include: (i) erbivores in lagoons prevent the overgrowth of corals by macro-algae [3]; (ii) roving herbivores influence the rate of grazing over depth gradients [20]; and (iii) coral reef resilience can be at risk even under low levels of fishing pressure on herbivorous fish [21]. Examples from habitats such as rocky reefs and seagrass meadows also demonstrate the critical role that herbivores play in controlling the growth of algae [22][23][24], a key process necessary for resilient ecosystems. Ensuring that important ecological functions, such as herbivory, are maintained in the face of disturbances is, therefore, important for managers of ecosystems that rely on such functions.
There is evidence to suggest that connectivity influences ecological functions [25]. Field studies have documented this effect with empirical data. For example, the connectivity of larval populations affects the resilience of algal-dominated coral reefs [14], and grazing by fish that use reefs can prevent disturbed reefs from shifting to an algal-dominated state [3]. Directly measuring and testing connectivity in the field can, however, be very difficult because ecological systems are comprised both of numerous different habitat types and a diversity of distinct habitat patches. These components are potentially linked by biological connections that operate at different scales [26,27]. Controlled experiments are an effective alternative to field experiments and can help advance theories about how ecosystem functions are linked to connectivity [28][29][30][31].
Here, we test the effect of changes in connectivity on an important ecosystem function in two microcosm experiments. The experiments were designed to measure how different levels of connectivity shape the response of a key ecological process, grazing of algae (herbivory), to a disturbance event (simulated algal bloom). The microcosms comprised multiple distinct patches and we varied the level of connectivity among patches using different spatial configurations of connections between patches (by varying both the number and location (configuration) of connections). We applied a spatially distinct disturbance event to part of the system (excess algal growth mimicking blooms that occur in eutrophic waters) and measured herbivory responses to press-and pulse-style disturbances in two separate experiments. We hypothesised that (i) patches with greater connectivity would have greater rates of herbivory, because they are more accessible to grazers, (ii) increases in connectivity across a whole system would increase herbivory, (iii) how the connections were configured would influence herbivory, and (iv) the influence of configuration on herbivory would be more visible during a press disturbance.

Materials and Methods
We tested the effect of connectivity on herbivory by placing mesograzers (yellow-footed hermit crabs, Clibanarius virescens) into microcosms with different spatial configurations of patches. We then measured herbivory rates under two disturbance scenarios and under seven different types of connectivity between patches. We selected the yellow-footed hermit crab as our model species because it is commonly found on intertidal and subtidal reefs in tropical to warm-temperate regions in the Indo-Pacific region, with demonstrated evidence that disturbances affect their grazing performance [32]. Disturbance in the current study was in the form of excess algal pellets added to specific patches within the system. This disturbance represents a scenario where, for example, parts of an embayment close to major nutrient inputs (e.g., estuarine plumes, sewage or storm-water effluent) are affected by algal blooms [33,34], whereas adjacent areas are not substantially altered by excess algae (for example, Yabsley et al. [35]).

System Design
We created microcosms with different features to test how connectivity affects mesograzer consumption under conditions that mimic systems partially affected by algal blooms. Microcosms were constructed with nine equally sized patches (patch size: 12 × 12 cm; system size: 36 × 36 cm), arranged in a square, and divided by white acrylic sheets (3 mm thick, 5 cm high; see Figure S1 for an example system). System connectivity was manipulated by varying the number and location of connections between patches; connections were openings in the wall (3 × 3 cm) that allowed the animals to move directly between adjacent patches. We added excess algal pellets to certain patches within the microcosm, mimicking algal bloom conditions in a disturbed system but leaving part of the system in an undisturbed state (Figure 1). The system was spatially configured so that all disturbed patches were connected to each other, whereas the spatial arrangement of connections between the undisturbed and disturbed parts varied among systems ( Figure 1). at which crabs consumed algal pellets (mean = 2.7, SE 0.1) in pilot studies (n = 10 trials), to ensure that consumers had the capacity to measurably influence the number of algal pellets introduced. .

Figure 1.
Experimental microcosms with different levels of connectivity that did not mirror each other. Connectivity increases with respect to an increasing number of links between the undisturbed and disturbed part of the system, and with respect to alternative spatial configurations for a given number of links. Connectivity increases with respect to an increasing number of links between the undisturbed and disturbed part of the system, and with respect to alternative spatial configurations for a given number of links.
By changing the number of connections and the way they were configured between undisturbed and disturbed patches (Figure 1), we were able to use two connectivity metrics at both the patch and whole-of-system scales (see Appendix A), and we used network theory to quantify connectivity. Network theory, a form of graph theory, is used to quantify connectivity in a wide variety of systems, including social systems [36], communication and transport systems [37], and ecosystems [38,39]. Metrics have been used to index multiple aspects of connectivity, ranging from simple measures such as the number of connections directly linking a single patch, to integrated measures of connections across the entire spatial configuration of a system [27,40]. The metric of closeness centrality is the inverse of the sum of the shortest path lengths between a node of interest and all other nodes in the system [36]. Closeness centrality provides a useful description of the patterns of connectivity likely to be found in the field [41].
Using multiple metrics to analyse the connectivity of a system or habitat patch can provide detailed information to assist in the management of anthropogenic impacts on connectivity [27]. The first metric of connectivity used in our experiments at the patch scale was the number of connections that linked disturbed patches to neighbouring patches (Table 1; Appendix A). At the system scale, the first metric of connectivity was the number of connections linking the undisturbed to the disturbed part, yielding four possible configurations (Figure 1, Figure A1). Next, we created alternative configurations of connections between the undisturbed and disturbed part of the systems, excluding arrangements that mirrored each other, while keeping to a practical number of systems amenable for testing, and this resulted in a total of seven systems (Figure 1, Appendix A). We quantified the configuration using the network theory metric, closeness centrality [36] for each disturbed patch and at the whole-of-system scale (see Appendix A for details on how these metrics were calculated for our experimental systems). The configuration metric resulted in values ranging between 3.4 and 8.3 for patch-scale and between 18.1 and 58.3 for system-scale, where patches or systems with lower values are less connected than those with higher values (Table 1; Appendix A). Closeness centrality accounts for the topographical distance between patches of interest; thus, our metric for configuration accounted for the placement of the patches and the different distances that the animals needed to move to reach the patches with excess algae. For example, patch B (Table 1) does not have any direct connections from undisturbed patches and the consistently lower values for configuration reflect this; however, in one situation (model system one), patch D is more isolated than patch B and, thus, patch D has a lower value for configuration.
Crabs of approximately the same size were collected from Burleigh Heads, Australia, in May 2016. Crabs were acclimatised in a controlled-temperature laboratory at Griffith University on the Gold Coast, Australia, to the following average conditions: Temperature 24.8 • C, salinity 35.5 ppt, pH 7.8, 12 h night-day cycle from 6 am to 6 pm. Crabs were housed together to allow the opportunity for interactions, limiting the likelihood of behaviour changes in response to unfamiliar individuals being a factor during experiments. Crabs were fed every four days on the same algal pellets used during experiments (New Life Spectrum Marine Fish Formula, 1 mm sinking pellets). Experiments were conducted on scheduled feeding days during June 2016 between 08:00 and 16:00 (hours where crabs were observed to be active) inside the same controlled-temperature laboratory where they were housed. Individuals for experiments were randomly selected from aquaria prior to feeding on scheduled feeding days (individuals observed to be moulting, carrying eggs or have moved to a larger shell were excluded from experiments to account for any behaviour or size change that might affect their consumption). The model systems were placed into separate containers and filled with filtered sea water collected from the Gold Coast seaway to 4 cm (more than enough for all crabs to be fully submerged) under the same conditions (temperature/salinity/pH) they were acclimatised to. One crab was placed into each of the nine patches. We measured the rates (number of pellets h −1 ind −1 ) at which crabs consumed algal pellets (mean = 2.7, SE 0.1) in pilot studies (n = 10 trials), to ensure that consumers had the capacity to measurably influence the number of algal pellets introduced. Table 1. Experimental system setups with patch-and system-scale connectivity values for the 'number of connections' and 'configuration' of the connections. Number of connections at the system scale is the number of links between disturbed patches (A-D) and undisturbed patches (grey), and number of connections at the patch scale is the number of neighbouring patches to the patch of interest (A-D). Values for configuration at the patch scale and system scale are calculated using an adaptation of the network theory measure of closeness centrality from Freeman [36]. Red arrows indicate open connections.    Table 1. Experimental system setups with patch-and system-scale connectivity values for the 'number of connections' and 'configuration' of the connections. Number of connections at the system scale is the number of links between disturbed patches (A-D) and undisturbed patches (grey), and number of connections at the patch scale is the number of neighbouring patches to the patch of interest (A-D). Values for configuration at the patch scale and system scale are calculated using an adaptation of the network theory measure of closeness centrality from Freeman [36].

Model
Red arrows indicate open connections.

Model System
System Scale Patch  Table 1. Experimental system setups with patch-and system-scale connectivity values for the 'number of connections' and 'configuration' of the connections. Number of connections at the system scale is the number of links between disturbed patches (A-D) and undisturbed patches (grey), and number of connections at the patch scale is the number of neighbouring patches to the patch of interest (A-D). Values for configuration at the patch scale and system scale are calculated using an adaptation of the network theory measure of closeness centrality from Freeman [36].
Red arrows indicate open connections.

Model System
System Scale Patch

Press Versus Pulse Herbivory Experiments
We applied two different styles of disturbance (separately) to the microcosms to test whether responses to press-style disturbances were different to pulse disturbances. In both experiments, each of the seven systems were replicated six times, yielding 42 measurements for system-level analyses and 168 observations for patch-level analyses. A summary of the length of the experiment and distribution of algal pellets can be found in Table 2. Table 2. Summary of experimental parameters for microcosms testing the effect of connectivity on herbivory. Microcosms were designed so that connectivity varied between disturbed and undisturbed patches in various ways.

Experiment Pulse Press
Length of experiment 1 h 7 h No. of crabs per patch 1 1 No. of patches per system 9 9 No. of pellets added per hour to undisturbed patches (×5) 1 1 No. of pellets added per hour to disturbed patches (×4) 13 9 In both the press and pulse disturbance experiments, one algal pellet was added to undisturbed patches because it is less than the normal consumption of a single crab (3 pellets h −1 expected under optimal conditions [32]). In an analogous field situation, this represents a patch where competition for food might be high, or resident grazers are already keeping algae levels low, and thus, some animals need to move elsewhere to find food.
In the pulse disturbance experiment, each of the disturbed patches received 13 pellets, ensuring there was too much for only the original crab to consume within the experimental period. Therefore, grazer mediation of the pulse disturbance relied on the movement of at least five other crabs into the disturbed patches relatively quickly for the system to return to its original state.
In the press disturbance experiment, each of the disturbed patches received nine pellets every hour for a seven-hour period, again ensuring there was too much algae for only the original crab to consume within the experimental period. This represented a slower onset and less intense disturbance than the pulse experiment. Therefore, grazer mediation of the press disturbance relied on the movement of fewer (three) other crabs into the disturbed patches, and for this level of herbivory to continue in the disturbed patches every hour for seven hours.

Statistical Analyses
We aimed to assess whether the total number of connections and/or the network configuration (as our key measures of system connectivity) were affecting the herbivory rates of our model animals.
To do this, we used generalised linear mixed-effect models (GLMMs) to test how the response metric (number of pellets grazed per patch out of the number added) varied with connectivity at the patch and system scale. We considered four explanatory variables for herbivory: i) The number of connections at the patch level, and ii) the value for configuration at the patch level, and then (i) and (ii) at the system level. We included 'replicates' as a random effect, to account for the effect of individual crab behaviour. We aimed to identify a best model but could not test all combinations of explanatory variables because the patch-and system-level variables were collinear. Therefore, we used a forward stepwise model simplification procedure: We tested the four metrics as lone main effects, then identified which patchand system-level metrics had the lowest Akaike Information Criterion (AIC). We then took those 'best' main effects and tested two more models with additive and interactive effects of the best patch-and system-level metrics. All models were run with the binomial distribution for observations in the lme4 package [42] in the software R ver. 3.3.1. [43]. The design of statistical analyses was the same for the pulse and press experiments, but these were analysed separately.

Results
Higher connectivity significantly (p < 0.05) increased herbivory. This was true regardless of whether we measured herbivory at the patch or system scale or whether we applied a pulse or press disturbance. The best (based on AIC) models at the patch and system scale were used to identify additive and interactive effects; these models differed between press and pulse disturbances (Table 3). Table 3. The number of connections and configuration at the patch and system scale significantly increased herbivory (p < 0.05 for all models). Best models for herbivory in response to pulse and press disturbances are in bold, determined using the difference in Akaike Information Criterion (AIC) values from generalised linear mixed models; ∆i. Akaike weights (AIC W ) are used to show the probability that the model is better than those being compared. No. of connections = number of connections between undisturbed and disturbed patches at the system level or number of direct connections to individual patches.

Pulse Disturbance Experiment
At the patch and system scale, both the number of connections and their configuration were significant (p < 0.05) predictors of herbivory in response to a pulse disturbance event (Table 3, Figure 2). Using the strongest predictors (determined using lowest AIC: Table 3) at the patch (number of connections) and system (configuration) scale, the additive and interactive effects were significant (p < 0.05). Therefore, herbivory, in response to a pulse disturbance, was best-predicted by the interaction between the configuration of connections at the system scale and the number of individual connections at the patch scale (lowest AIC value, Table 3, Figure 2). . Figure 2. Herbivory during a pulse disturbance and the interaction between patch-and system-scale responses to an increase in connectivity. Herbivory is the mean (SE) number of algal pellets consumed after a 1 h period (colours indicate no. of connections at the patch scale). No. of connections (patch) is the number of connections that linked disturbed patches to neighbouring patches (disturbed patches had a minimum of two connections). Configuration (system) is a modified value for closeness centrality at the system scale, which quantifies a commonly used dimension of connectivity.

Press Disturbance Experiment
At the patch and system scale, both configuration and number of connections were significant (p < 0.05) predictors of herbivory in response to a press disturbance event (Table 3, Figure 3). Using the strongest predictors (determined using lowest AIC: Table 3) at the patch (configuration) and system (number of connections) scale, the additive and interactive effects were significant (p < 0.05). Therefore, herbivory, in response to a press disturbance, was best-predicted by the interaction between the configuration of connections at the patch scale and the number of connections linking undisturbed to disturbed parts at the system scale (lowest AIC value, Table 3, Figure 3).

Figure 2.
Herbivory during a pulse disturbance and the interaction between patch-and system-scale responses to an increase in connectivity. Herbivory is the mean (SE) number of algal pellets consumed after a 1 h period (colours indicate no. of connections at the patch scale). No. of connections (patch) is the number of connections that linked disturbed patches to neighbouring patches (disturbed patches had a minimum of two connections). Configuration (system) is a modified value for closeness centrality at the system scale, which quantifies a commonly used dimension of connectivity.

Press Disturbance Experiment
At the patch and system scale, both configuration and number of connections were significant (p < 0.05) predictors of herbivory in response to a press disturbance event (Table 3, Figure 3). Using the strongest predictors (determined using lowest AIC: Table 3) at the patch (configuration) and system (number of connections) scale, the additive and interactive effects were significant (p < 0.05). Therefore, herbivory, in response to a press disturbance, was best-predicted by the interaction between the configuration of connections at the patch scale and the number of connections linking undisturbed to disturbed parts at the system scale (lowest AIC value, Table 3, Figure 3). Diversity 2020, 12, x FOR PEER REVIEW 9 of 18 Figure 3. Herbivory during a press disturbance and the interaction between patch-and system-scale responses to an increase in connectivity. Herbivory is the mean (SE) number of algal pellets consumed after a 7 h period of a consistent supply of algae to disturbed patches in each system (colours indicate no. of connections at the system scale). No. of connections (system) is the number of connections linking the undisturbed to the disturbed part of the system. Configuration (patch) is a modified value for closeness centrality at the patch scale, which quantifies a commonly used dimension of connectivity.

Variability in Grazing Rates
Grazing rates, at both the patch and system scale, were variable, possibly as a result of behavioural interactions (Figures 2 and 3). We accounted for differences in the behaviour of crabs by including the replicate as a random effect in the analysis. It is, however, possible that some interactions between crabs were not explicitly addressed in the analysis. For example, competition may have impeded consumption rates by competitively inferior individuals in some patches; this is likely to occur in the field, and hence, our experiments mimic this scenario. Crowding near food patches may have constrained new arrivals to quickly access algal pellets, possibly lowering overall consumption rates for patches with many individuals. By contrast, increasing the number of connections provided alternative pathways to reach food patches, thereby making it more probable that more crabs would be the 'first-on-scene' at food patches. Herbivory during a press disturbance and the interaction between patch-and system-scale responses to an increase in connectivity. Herbivory is the mean (SE) number of algal pellets consumed after a 7 h period of a consistent supply of algae to disturbed patches in each system (colours indicate no. of connections at the system scale). No. of connections (system) is the number of connections linking the undisturbed to the disturbed part of the system. Configuration (patch) is a modified value for closeness centrality at the patch scale, which quantifies a commonly used dimension of connectivity.

Variability in Grazing Rates
Grazing rates, at both the patch and system scale, were variable, possibly as a result of behavioural interactions (Figures 2 and 3). We accounted for differences in the behaviour of crabs by including the replicate as a random effect in the analysis. It is, however, possible that some interactions between crabs were not explicitly addressed in the analysis. For example, competition may have impeded consumption rates by competitively inferior individuals in some patches; this is likely to occur in the field, and hence, our experiments mimic this scenario. Crowding near food patches may have constrained new arrivals to quickly access algal pellets, possibly lowering overall consumption rates for patches with many individuals. By contrast, increasing the number of connections provided alternative pathways to reach food patches, thereby making it more probable that more crabs would be the 'first-on-scene' at food patches.

Discussion
We show that herbivory increases with higher connectivity under both press and pulse disturbances in a controlled environment using a single model species. We tested two complementary metrics of connectivity (the number of connections and configuration); both were important predictors of herbivory at two spatial scales, the system and the patch. Herbivory by the hermit crab, Clibanarius virescens, was greater in better-connected systems and our results show that herbivory is best-predicted using both patch-and system-scale metrics of connectivity simultaneously, providing support for our first three hypotheses. We also found that configuration had a greater influence on herbivory during the press disturbance, supporting our fourth hypothesis. During the pulse disturbance, once a patch had three or four connections, the configuration of the system connections had little influence on herbivory in C. virescens (Figure 2). On the other hand, during the press disturbance, regardless of the number of system connections, the configuration of the patches continued to influence herbivory rates in C. virescens (Figure 3).
The finding that both the number and configuration of connections were important predictors for herbivory in disturbed experimental microcosms is supported by previous work. For example, the analysis of wildlife crossings planned for mitigation of the Trans-Canada highway upgrade resulted in two recommendations to increase potential movement of key species: (1) Increase the number (i.e., number of connections); and (2) change the location (i.e., configuration of connections) of crossings [44]. Analysing the effects of spatial configuration at the level of individual habitats, and at the landscape scale, has proven to be important in both terrestrial [44,45] and marine ecosystems [46][47][48]. The spatial configuration of habitat patches within marine networks is a key consideration for managers interested in improving the success of larval dispersal across marine networks [40]. The performance of restored habitats is also predicted to depend on the placement of restoration initiatives at sites with configurations that enhance the exchange of energy or individuals with other habitat patches [47,49]. Landscape configuration can also be used to predict movement patterns of important waterbird species, which help to maintain ecosystem functioning in fragmented landscapes [50]. The configuration of connections is an important consideration because it can influence how resources are accessed, used or even isolated within landscapes, thereby potentially contributing to the resilience of individual ecosystems and habitat patches [51].
The interaction between metrics for connectivity best-predicted the herbivory in response to both press and pulse disturbances. Determining which metric is the best predictor for a function, and considering connectivity at both the system and patch scale, might therefore be an important consideration for the management of any ecosystem [52,53]. For example, during the pulse disturbance, when a patch had a low number of connections (two) and low system configuration (20), grazing was low (<8 pellets h −1 ). However, in a patch with the same low number of connections (two), but high system configuration (60), more grazing occurred (>11 pellets h −1 ) (Figure 2). Knowledge of the spatial configuration of an entire system can, therefore, provide more accurate information about how much grazing might occur in a patch of interest, compared with simply measuring connections to neighbouring patches. This is supported by previous findings, which show that metrics indexing both system and patch measures of connectivity led to different conclusions about the viability of a population subject to habitat fragmentation, when compared with results from patch measures alone [54]. During the press disturbance, although the system with three connections had a high rate for grazing overall, patches with a lower value for configuration in that system had lower grazing rates than patches in a system with only one connection (Figure 3). These findings provide further evidence for management to consider whether they are interested in the functional response of a patch of interest, or the functional response of an overall system. Finally, when connectivity was high at both the patch and system scale, regardless of whether it was a press or pulse disturbance, herbivory remained high.
The variability in grazing rates, during both the press and pulse disturbance experiments, primarily reflects the complex nature of the relationship between connectivity and herbivory. It may, however, be further influenced by behavioural responses. Competition for a common resource by different species of herbivores in natural ecosystems is an important process that can control the growth and cover of algae in lake [55] and marine ecosystems [56]. The species of herbivore can also impact the type of algal consumption needed to maintain healthy reefs [57]. We acknowledge that the relationship between herbivory and the behaviour of herbivores is a complex process and can differ depending on the herbivore of interest, e.g., insect versus mammalian herbivory [58]. Therefore, it would be important to take into account the species of herbivores present in a natural ecosystem and their known behaviours when using this study as a baseline to test different levels and configurations of connectivity on the ecosystem function of herbivory.
Our results demonstrate that herbivory responds positively to increased connectivity and we have used herbivory to represent a positive effect that counteracts disturbance (i.e., excess algal growth) (Figure 4a). However, responses depend on whether the function responds positively or negatively to a disturbance in the system and the number of disturbances present in the system (summarised in Figure 4). In other situations, grazing can be detrimental to ecosystem integrity, such as excessive activity of sea urchins on coral reefs that can lower coral recruitment [59,60], which would result in a negative relationship with increasing connectivity (Figure 4b). Testing two (or more) opposing stressors together has the potential to determine whether intermediate levels of connectivity maximise resilience. For example, a combination of the positive effect of herbivore grazing for reef resilience [61] with the negative impacts from the spread of an invasive species [62] might result in a unimodal response (e.g., Figure 4c). A unimodal relationship has also been modelled in food web studies to highlight the opposing effects of prey protection with low connectivity and prey exploitation at higher levels of connectivity [63]. Further investigation is required to determine the impact that two opposing stressors might have in response to increasing connectivity.
A widely propagated design principle for experiments in ecology emphasises the pivotal role of robust inferences based on experimental conditions. Typically, ecologists seek to identify general relationships and distinguish such 'generalities' from inferences based on an idiosyncratic and narrow range of conditions [64]. It must be acknowledged that whilst no single experimental design will always be sufficient to define the broader generalities that ecologists seek [65], designs that achieve transferability of results to a broad range of settings and conditions will offer more robust and widely applicable insights. We used a single ecological function, herbivory, to test the effects of connectivity on ecological processes. Increasing connectivity in our experiment increased the rate of herbivory, and it is not implausible that other functions (e.g., predation, scavenging) may respond to connectivity in a similar manner under similar experimental conditions. For example, the availability of open corridors (higher connectivity) had a positive outcome for predators in a metacommunity experiment [12], further highlighting that extending our experimental design to other consumers and resources represents fertile ecological ground for future work.
Our conclusions are drawn from an experimental test performed in a controlled environment on one species, so further investigation is warranted to expand the findings beyond the limitations of a microcosm approach. The use of agal pellets as a proxy for herbivory in our model systems is a simplified way of measuring grazing. However, we feel its use is justified for several reasons. The algal pellets we used have a composition that reflects the natural algal diet of our model species and are of standardised weight and size, making it possible to accurately quantify and compare small changes in consumption rates between repeat measurements in a way that is difficult to achieve with live algae. Given our species responded well to this food source during acclimatisation, we believe it serves as a practical way of measuring consumption for the purposes used in this study. Future studies could consider building upon these findings by adding more complexity to food sources and model systems or by adding stressors. Despite the limitations and simple nature of our study, the techniques we have applied here could be replicated to test the way that other ecosystem functions and different species will respond to an increasing number and configuration of connectivity under other, single and multiple stressors/disturbances, and at different magnitudes. Different ecosystem processes, such as dispersal, habitat preference, competition, scavenging and predation, could affect these conclusions because ecosystems exhibit a wide range of characteristics [66]. Multiple stressors occurring simultaneously in a system may also result in antagonistic, synergistic or additive effects [67]. It would therefore be useful to test the connectivity relationship using multiple stressors that are likely to work antagonistically, synergistically and/or additively to determine how the functional response changes.
Diversity 2020, 12, x FOR PEER REVIEW 12 of 18 likely to work antagonistically, synergistically and/or additively to determine how the functional response changes.

Conclusions
Our study suggests that connectivity at the system and patch scale are both important considerations for understanding grazer-mediated responses to press and pulse disturbances. We successfully designed and tested microcosm systems that allow testing of the ways that simple ecological systems respond to changes in linkages between spatially distinct units (mimicking habitat patches, islands or other parts of mosaics). Our experimental design opens opportunities to test multiple concepts of network theory for ecological systems in a way that is rarely possible to perform in the field (for logistical and ethical reasons). Furthermore, our experiments have practical applications for real-world scenarios. For example, in situations where effort has been made to restore

Conclusions
Our study suggests that connectivity at the system and patch scale are both important considerations for understanding grazer-mediated responses to press and pulse disturbances. We successfully designed and tested microcosm systems that allow testing of the ways that simple ecological systems respond to changes in linkages between spatially distinct units (mimicking habitat patches, islands or other parts of mosaics). Our experimental design opens opportunities to test multiple concepts of network theory for ecological systems in a way that is rarely possible to perform in the field (for logistical and ethical reasons). Furthermore, our experiments have practical applications for real-world scenarios. For example, in situations where effort has been made to restore important habitats that have failed (with little understanding as to why), our experiments show that there is potential value in analysing connectivity and configuration together. This type of analysis may help to uncover reasons for conservation failure where connectivity to other habitats is important for resistance or recovery from a disturbance. In a world where effective and efficient conservation designs will become increasingly critical in the face of mounting human threats, enhancing such designs by applying more robust connectivity criteria will be a useful input that can benefit from conceptual advances made in controlled experiments.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A.
Appendix A.1. Network Theory Metrics Used to Measure Configuration at the Patch and System Scale Closeness centrality [1] is a network theory metric that we adapted to calculate a value for the configuration of the patches and systems for analysis of our experimental microcosms. Closeness centrality is the inverse of the sum of the shortest path between the patch of interest and all other patches in the system, measured as the number of links (or topographical distance) in the microcosms used here.

Appendix A.2. Patch-Scale Configuration Metric
At the patch scale, we adapted closeness centrality to the following measure for configuration: where: p k = disturbed patch of interest (A, B, C or D ( Figure A1)) p i = all undisturbed patches (U1, U2, U3, U4 and U5 ( Figure A1)) d p i , p k = the number of links (topographical distance) between p k and p i We modified the original equation for closeness centrality developed by Freeman [36] in two ways: 1.
In the original equation, p i refers to all other patches within the system; however, because we were only interested in the grazing power of the crabs coming from the undisturbed patches to each of the disturbed patches, our measure for patch configuration considers only the connectivity between the disturbed patch of interest and all undisturbed patches ( Figure A1);

2.
Due to the limited number of patches and overall size of our system, closeness values were in the order of 1 × 10 −3 ; therefore, we multiplied the original equation by 100 to obtain an easily interpretable value.

Appendix A.3. System-Scale Configuration Metric
At the whole-of-system scale, we varied the number and spatial arrangement of the connections linking disturbed and undisturbed patches. We combined two measures to represent 'configuration' at the system scale, to account for the way these connections increased in connectivity, and to account for the direction of the migration of animals through the system: (1) the average configuration of each undisturbed patch to each of the disturbed patches; and (2) the probability that a crab would walk in the right direction toward a disturbed patch to graze.
For the first measure, the average configuration, we were interested in how closely connected all of the undisturbed patches were to the disturbed patches; thus, we adapted closeness centrality as follows: x ̅ closeness (U1: U5) = [∑ d(p u , p d ) where: p u = undisturbed patch of interest (U1, U2, U3, U4 or U5 ( Figure A2)) p d = all disturbed patches (white patches in Figure A2) d(p u , p d ) = the number of links (topographical distance) between p u and p d (a) (b) Figure A2. (a): Example showing the patches (indicated by arrows) used to calculate configuration for U1 (circled) (e.g., distance between U1 and all disturbed patches (in white)), with the average distance from U1 to U5 used to provide a configuration value; (b): Example system showing the probability in each undisturbed patch that the crab will move in the right direction toward a disturbed patch.

Appendix A.3. System-Scale Configuration Metric
At the whole-of-system scale, we varied the number and spatial arrangement of the connections linking disturbed and undisturbed patches. We combined two measures to represent 'configuration' at the system scale, to account for the way these connections increased in connectivity, and to account for the direction of the migration of animals through the system: (1) the average configuration of each undisturbed patch to each of the disturbed patches; and (2) the probability that a crab would walk in the right direction toward a disturbed patch to graze.
For the first measure, the average configuration, we were interested in how closely connected all of the undisturbed patches were to the disturbed patches; thus, we adapted closeness centrality as follows: x closeness (U1 : U5) = where: p u = undisturbed patch of interest (U1, U2, U3, U4 or U5 ( Figure A2)) p d = all disturbed patches (white patches in Figure A2) d p u , p d = the number of links (topographical distance) between p u and p d i=1 where: p u = undisturbed patch of interest (U1, U2, U3, U4 or U5 ( Figure A2)) p d = all disturbed patches (white patches in Figure A2) d(p u , p d ) = the number of links (topographical distance) between p u and p d (a) (b) Figure A2. (a): Example showing the patches (indicated by arrows) used to calculate configuration for U1 (circled) (e.g., distance between U1 and all disturbed patches (in white)), with the average distance from U1 to U5 used to provide a configuration value; (b): Example system showing the probability in each undisturbed patch that the crab will move in the right direction toward a disturbed patch. Figure A2. (a): Example showing the patches (indicated by arrows) used to calculate configuration for U1 (circled) (e.g., distance between U1 and all disturbed patches (in white)), with the average distance from U1 to U5 used to provide a configuration value; (b): Example system showing the probability in each undisturbed patch that the crab will move in the right direction toward a disturbed patch.
For the second measure, we then summed the probability that each of the five crabs in undisturbed patches (U1:U5) would take the correct path toward a disturbed patch to graze ( Figure A2). For example, if a patch has two connections and only one of those connections takes the crab in the right direction (towards a disturbed patch with excess algae) to graze, then that patch has a probability of 1 2 ; therefore, our equation is as follows: No. o f correct paths Total number o f paths (A3) Finally, these two measures were multiplied to provide our overall metric for configuration at the system scale: Configuration (system) = x closeness (U1 : U5) × P U1:U5 (A4) Connectivity metrics for each of the systems at the patch and system scale are shown in Table 1.