Biodiversity Patterns of Macroinvertebrate Assemblages in Natural and Artificial Lentic Waters on an Oceanic Island

The Azorean islands have been historically affected by human activities, mainly due to the combined effects of habitat degradation and fragmentation, and the introduction of exotic species. We here aim to analyze the role of environmental characteristics and spatial descriptors in supporting regional biodiversity of macroinvertebrates by considering natural ponds and artificial tanks. After the monthly variation of macroinvertebrate assemblages was assessed in three temporary and two permanent ponds in the Azorean island of Terceira during a complete inundation-desiccation annual cycle, the assemblage differences of 12 ponds (three temporary and nine permanent ponds) and 8 closely-located artificial tanks were analyzed across a range of landscape disturbances. Macroinvertebrate assemblages were found to differ according to hydroperiod and sampled months. Although the former explained the highest variance, macroinvertebrate differentiation by hydroperiod was also dependent on the study month. Our results also revealed a consistent monthly pattern of species replacement. However, the contribution of nestedness to the macroinvertebrate β-diversity was notable when temporary ponds were close to desiccation, probably indicating a deterministic loss of species due to the impoverished water conditions of the ponds facing desiccation. When the macroinvertebrate assemblages were analyzed in relation to physico-chemical variations and spatial descriptors, the artificial tanks were not clearly segregated from the natural ponds, and only differentiated by pH differences. In contrast, those natural ponds exhibiting high concentrations of total phosphorous (likely signs of anthropization) also discriminated the ordination of ponds in a distance-based redundancy analysis, and showed impoverished assemblages in comparison with well-preserved ponds. The macroinvertebrate assemblages of the natural ponds showed a significant spatial pattern, but this spatial influence was not significant when tanks and ponds were considered together. Our results suggest that tanks may act as possible reservoirs of biodiversity during the desiccation period of temporary ponds, but are unable to establish successful populations. These fishless permanent tanks can complement the conservation of a biodiversity that is largely maintained by the pristine high-altitude natural ponds. The establishment of a guideline for conservation management that also considers the artificial tanks is necessary to benefit the local and regional Azorean macroinvertebrate diversity.

The Azorean islands have been historically affected by human activities, mainly due to the combined effects of habitat degradation and fragmentation, and the introduction of exotic species. We here aim to analyze the role of environmental characteristics and spatial descriptors in supporting regional biodiversity of macroinvertebrates by considering natural ponds and artificial tanks. After the monthly variation of macroinvertebrate assemblages was assessed in three temporary and two permanent ponds in the Azorean island of Terceira during a complete inundation-desiccation annual cycle, the assemblage differences of 12 ponds (three temporary and nine permanent ponds) and 8 closely-located artificial tanks were analyzed across a range of landscape disturbances. Macroinvertebrate assemblages were found to differ according to hydroperiod and sampled months. Although the former explained the highest variance, macroinvertebrate differentiation by hydroperiod was also dependent on the study month. Our results also revealed a consistent monthly pattern of species replacement. However, the contribution of nestedness to the macroinvertebrate β-diversity was notable when temporary ponds were close to desiccation, probably indicating a deterministic loss of species due to the impoverished water conditions of the ponds facing desiccation. When the macroinvertebrate assemblages were analyzed in relation to physico-chemical variations and spatial descriptors, the artificial tanks were not clearly segregated from the natural ponds, and only differentiated by pH differences. In contrast, those natural ponds exhibiting high concentrations of total phosphorous (likely signs of anthropization) also discriminated the ordination of ponds in a distance-based redundancy analysis, and showed impoverished assemblages in comparison with well-preserved ponds. The macroinvertebrate assemblages of the natural ponds showed a significant spatial pattern, but this spatial influence was not significant when tanks and ponds were considered together. Our results suggest that tanks may act as possible reservoirs of biodiversity during the desiccation period of temporary ponds, but are unable

INTRODUCTION
Freshwater ecosystems are currently highly vulnerable to external perturbations associated with anthropogenic changes, mainly due to the introduction of exotic species, landscape transformations and climate change (Sala et al., 2000;Pyšek et al., 2010). Indeed, landscape transformation (e.g., to urban or agricultural areas) has largely contributed to the global loss and degradation of freshwater habitats and is a major threat to freshwater organisms (McKinney, 2002). In contrast to the general negative role of human disturbance in shaping biodiversity, artificial ponds or tanks associated with agricultural systems have sometimes demonstrated that they are able to sustain biodiversity (Abellán et al., 2006;Declerck et al., 2006;Céréghino et al., 2008;Thiere et al., 2009). These humanmade ponds usually have different environmental characteristics in comparison with natural ponds (e.g., Hill et al., 2016, but see Deacon et al., 2018Deacon et al., , 2019 for similar chemical properties), given that these artificial systems usually have concrete sides and reduced vegetation cover, and possibly higher contaminant inputs compared to natural ponds (Hassall, 2014). Artificial ponds usually show reduced freshwater biodiversity, can support the occurrence of exotic species and may act as possible ecological traps (Oertli and Parris, 2019). However, artificial ponds that maintain good water quality (e.g., low concentrations of nutrient and other pollutants, vegetated bed/margins, low electrical conductivity) may provide habitat diversification, refuge and supplies for aquatic macroinvertebrates, amphibians, and terrestrial species (Mitsch and Gosselink, 2000;Oertli and Parris, 2019). Hence, artificial permanent waterbodies for cattle drinking-water have received increasing attention as regards conservation due to their secondary function as refuge habitat for macroinvertebrates (Ruggiero et al., 2008). In the current scenario of a changeable world, decisions regarding humanmade ponds are essential for conservation programs. However, the particular environmental conditions that favor or deplete biodiversity in artificial ponds across different land-uses, and their positive or negative contribution to biodiversity in a natural pond network remain quite unexplored.
Temporary ponds are characterized by recurrent inundation and desiccation (Williams, 1997;Florencio et al., 2011;Céréghino et al., 2012), but they can usually persist for centuries in the same region, in addition to possessing a worldwide distribution (Williams, 1997;Williams et al., 2001). However, the high and valuable biodiversity (high species richness in relation with the size of the pond, rare taxa, and uniqueness) of temporary ponds contrasts with their sensitivity and vulnerability to external perturbation (Williams, 2006). Consequently, temporary ponds are considered priority habitats for conservation by the European Union (code 3170 of the Habitat Directive). The inundation of temporary ponds usually starts at the onset of the rainy season, whereas the duration of the recurrent dry period is less predictable (Williams, 1997). These temporary ponds support singular macroinvertebrate taxa that often cannot survive in other types of aquatic ecosystems (Collinson et al., 1995;Bilton et al., 2001;Williams, 2006). The macroinvertebrate species that inhabit these ponds must cope with pond desiccation by adjusting their life cycles to the water permanence (hydroperiod) and employing particular strategies to survive pond desiccation (Williams, 2006). Moreover, these macroinvertebrates usually cannot cope with the presence of fish, having evolved in the absence of these top predators, which in most occasions are restricted to permanent ponds (Wellborn et al., 1996). In addition, the heterogeneity of a pond network is reflected in its hydroperiod gradient (Florencio, 2010). Therefore, widely variable hydroperiods, together with good pond connectivity, are key factors that contribute to the conservation of a high diversity of macroinvertebrates (Urban, 2004;Jeffries, 2005).
The proper management and conservation of biodiversity against possible threats require the acquisition of a good understanding of the regional species diversity, and the measurement of β-diversity is a useful concept to assess that (Socolar et al., 2016) as it allows the proper assessment of seasonal and spatial changes in aquatic invertebrate assemblage compositions (Florencio et al., 2016a). In metacommunity ecology, β-diversity may be defined as the variation in species composition among sites in a geographical region (Legendre et al., 2005; but see e.g., Koleff et al., 2003;Tuomisto, 2010;Anderson et al., 2011). The observed differences in β-diversity values may be the result of the combination of species between sites (replacement), and of the loss (or gain) of species between sites (nestedness), which it is known as a form of β-diversity partitioning (Baselga, 2010;Baselga and Orme, 2012). Alternatively, Carvalho et al. (2012) also provide this partition of β-diversity, but incorporate a richness difference component instead of nestedness. In ponds, those systems with high species replacement would demand the conservation of several ponds with variable richness and environmental conditions, while systems with low species replacement and high nestedness contribution to the β-diversity would prioritize the conservation of those ponds with the highest richness (see Baselga, 2010).
The Azorean archipelago is an ideal model system to assess the impact of landscape alterations on biodiversity, given that the Azores has suffered intensive landscape transformations, as well as the introduction of exotic species (approximately 58% of the arthropods are exotic species; Borges et al., 2005). Such transformations from native to agricultural areas have drastically restricted the extent of the original native forests to high altitude areas, much less accessible to humans (Triantis et al., 2010;Borges et al., 2020). Natural freshwater ecosystems in Macaronesian archipelagos (Azores, Madeira, Selvagens, Canary Islands, and Cabo Verde, sensu Engler, 1914) are considered priority conservation areas by the EU Water Framework Directive (2000/60/EC) (WFD), and are mainly located in the high-altitude native forests. Although these ecosystems are essential for the conservation of a unique freshwater community, limited knowledge of their temporal and spatial dynamics hinders the establishment of proper conservation guidelines (Hughes and Malmqvist, 2005). The present study focuses on the landscape matrix of Terceira Island, Azores, where natural and artificial ponds are located, to analyze the role of environmental characteristics and humanmade ponds in supporting regional biodiversity. The spatiotemporal variations in the natural temporary and permanent ponds in Azores were first determined by using monthly data of three well-preserved temporary and two permanent ponds that were sampled monthly, during a complete year of inundationdesiccation cycling. Secondly, changes in macroinvertebrate assemblages were addressed across a range of anthropogenic disturbance, using all accessible ponds within Terceira Island as a study case, encompassing 12 natural ponds (three temporary and nine permanent ponds), and eight closely-located artificial tanks constructed to supply water for cattle. We hypothesized that (i) pond differences and seasonal variations in macroinvertebrate assemblages would indicate different pre-dominant taxonomical groups during the hydrological year, and that such differences would be explained by the pond characteristics; (ii) artificial tanks could act as reservoirs of aquatic biodiversity, and be inhabited by similar macroinvertebrate assemblages as natural ponds. To explore the latter, we have disentangled the environmental characteristics that favored macroinvertebrate diversity in these artificial tanks.

Study Site
The Azores is an archipelago located in the North Atlantic, about 1,600 km from the European coast and 3,900 km from North America, between latitudes 36 • 55 ′ -39 • 43 ′ N and longitudes 24 • 45 ′ -31 • 17 ′ W (Figure 1). It comprises nine main islands and some islets, all of volcanic origin, which are organized in three groups: the western group (Corvo and Flores islands); the central group (Faial, Pico, São Jorge, Graciosa, and Terceira islands); and the Eastern group (São Miguel and Santa Maria islands). The climate is temperate oceanic, characterized by stable temperatures, substantial precipitation (mean annual precipitation of 740-2,400 mm; Bettencourt, 1979) and high relative atmospheric humidity, which can reach more than 95% in high-altitude native forests.
This study was carried out in Terceira Island which has an area of 402 km 2 and a maximum elevation of 1,023 meters above sea level. The Azorean archipelago, and specifically Terceira Island, was colonized by the Portuguese in the 15th century. Since the first settlements, the landscape was gradually transformed to accommodate agricultural activities and intensive pasturelands for cattle and crops after 1,950 (Triantis et al., 2010). Native forests were gradually destroyed as a consequence of these major land-use changes and currently the few remnants of the original forest or even secondary patches are located in the most inaccessible areas of the island (Figure 1), representing scarcely 6% of the Terceira Island surface area (Gaspar et al., 2008). Eight of the 12 studied ponds were located in natural and seminatural areas of Terceira Island, at between 387 and 912 m.a.s.l., including two protected areas with still pristine native forest, considered two of the most pristine areas in the Azores (Gaspar et al., 2011): "Caldeira de Santa Bárbara e Mistérios Negros" and "Terra Brava." These ponds are small (<7,000 m 2 ) and shallow. The dominant aquatic vegetation is formed by the genera Juncus (Juncaceae), Eleocharis (Cyperaceae), and Sphagnum (Sphagnaceae), and by the species Littorella uniflora (Plantaginaceae), Hydrocotile vulgaris (Apiaceae), and Polytrichum commune (Polytrichaceae). The surrounding area is dominated mainly by native and endemic plants (e.g., Juniperus brevifolia, Laurus azorica, Ilex perado subsp. azorica, and Erica azorica), with a dominance of two types of forests: "Juniperus-Ilex Montane Forests" and "Juniperus Montane Woodlands" (Elias et al., 2016) that also include bryophyte communities on all substrates (Gabriel and Bates, 2005). The other four ponds and the eight tanks were located in human-modified landscapes (Figure 1, see details below).
Our study period was from November 2013 to August 2014. First, to understand the variations in the macroinvertebrate assemblages of the natural ponds on Terceira Island, monthly samples were collected in three temporary ponds (VF1, VF3, VF5) and two permanent ponds (NEG, SER) between November 2013 and August 2014 (Figure 1, see also Supplementary Figure 1, ESM1 in the Supplementary Material) (hereafter referred to as "monthly sampling"). These ponds were selected based on their hydroperiod gradient and spatial location within the protected area "Caldeira de Santa Bárbara e Mistérios Negros." Secondly, we sampled almost all natural ponds (12 out of the 16 known ponds as four of them were inaccessible; local names included in Supplementary Table 1, ESM1 in the Supplementary Material), and eight artificial tanks in Terceira Island in May 2014 (hereafter referred to as "extensive sampling") to assess the environmental characteristics that would determine differences in the macroinvertebrate assemblages across a gradient of landscape disturbance (Figure 1 The artificial tanks were located in agricultural areas, and were selected as the closest sites to a natural pond in order to consider their role as possible reservoirs of the aquatic fauna. The tanks have an artificial substrate of cement, with a surface area of 4-6 m 2 , and 0.4-1 m of maximum depth. (2) agricultural areas, including arable lands, fruit orchards and permanent crops; (3) exotic forests dominated mainly by Pittosporum undulatum; (4) forestry plantations of Cryptomeria japonica and Eucalyptus spp.; (5) intensive managed pasturelands exposed to intensive cattle grazing throughout the whole year; (6) seminatural pasturelands, which were only managed with low cattle grazing activity and usually located in high elevation areas; (7)

Sampling Methods
Macroinvertebrates were sampled using a dip net of 40 × 25 cm with 1 mm mesh size. Approximately 1.5 m stretches of water were netted in each sampling unit, which encompassed three successive sweep nettings, collecting benthic, and open-water macroinvertebrates. Sampling was performed from the littoral zone to the deepest point of the pond when accessible, and considering different aquatic plant covers to include all the possible microhabitats. Because the dip netting efficiency seems to be better in small ponds (Heyer et al., 1994), the number of sampling units was proportional to the pond size (see Florencio et al., 2009). Hence, the number of samples per pond ranged from 5 to 19 sampling units. In total, 557 samples were obtained during the monthly sampling and 169 during the extensive sampling. Most individuals were sorted in the field, counted, and released, to minimize the researcher's impact in the field. Individuals of unidentified species were preserved in 70% ethanol for subsequent identification in the laboratory. Individuals were identified to species level, except for the genus Dryops which was identified to genus level, Culicidae to family level and Oligochaeta to Class level; adult and larval stages were considered separately given their different ecological requirements (hereafter referred to as "species" for simplicity). The Dipteran family Chironomidae was not considered because of sampling and mesh size limitations.

Physico-Chemical Characteristics and Land-Use Variables of Ponds
The following environmental variables were measured in situ to characterize the study ponds during monthly and extensive samplings: electrical conductivity (EC) standardized at 20 • C (Multi-range Conductivity meter HANNA HI 98127), pH (pH meter HANNA HI 98311), dissolved oxygen concentration, and temperature (Oxi 315i WTW) in the water column. Two replicates were obtained for each measurement per pond and month, but given the homogeneity in the measurements, the two pond values per month were averaged before statistical analyses. We also collected 1,500 ml of the water column that was later filtered in the laboratory (Whatman GF/C filters, 47 mm diameter) to measure the concentrations of chlorophylla (spectrophotometric methods, Gonçalves, 2008), dissolved inorganic nitrogen (DIN) (following Golterman, 1991) and dissolved inorganic phosphate reactive to Mo (i-P) according to Murphy and Riley (1962). The concentration of total phosphorous in the water (TP) was also analyzed as i-P after acid digestion of the unfiltered water sample treated with K 2 S 2 O 8 (Golterman, 2004). DIN and i-P concentrations were only measured in the extensive sampling.
The land-use map of Terceira Island was constructed based on land-use cover classes of CORINE 2006 Land Cover (see Bossard et al., 2000) and DROTH (2008), at 30 × 30 meters resolution (Figure 1). The land-use map was divided into eight categories: (1) urban and industrial areas; (2) agricultural areas, including arable land, fruit orchards and permanent crops; (3) exotic forests dominated mainly by Pittosporum undulatum; (4) forestry plantations of Cryptomeria japonica and Eucalyptus spp.; (5) intensive managed pasturelands exposed to intensive cattle grazing throughout the whole year, characterized by common exotic grasses and leguminous forbs; (6) seminatural pasturelands, which were only managed with low cattle grazing activity and usually located in high altitude areas; (7) uncovered areas including rocky and open spaces without vegetation, and (8) native vegetation including native forests dominated by endemic plant species, but also including high-elevation anthropized pasturelands that have been abandoned and colonized by naturalized vegetation (native and exotic plants).
The percentage area occupied by the different land-uses in Terceira Island was calculated to quantify the impact of anthropogenic disturbance on the analyzed ponds and tanks. For this, each pond/tank center was established as the centroid of a 250-meter-radius buffer in order to estimate the local influence of the surrounding land-uses on the macroinvertebrate assemblages. Only intensively managed pasturelands, seminatural pasturelands and native vegetation surrounded the study ponds, and thus their percentage areas were considered as explanatory variables. The land-use map and the occupied percentage area calculation were performed using the GIS-based software ArcGIS 10 ArcMap v. 10.1.

Spatial Variables Based on Interpond Distances
Regarding the extensive sampling, we retained 11 orthonormal spatial descriptors of the 12 natural ponds and seven tanks (N = 19, one tank was excluded because no individuals were collected) using Principal Coordinates of Neighborhood Matrix (PCNMs) map distances. We used as truncation distance the longest distance among ponds to maintain pond connectivity (Borcard and Legendre, 2002). The extraction of these spatial descriptors was performed in R software 2.14.2, using the "pcnm" command in the "vegan" package (Oksanen et al., 2019).

Temporal Biodiversity Patterns
For the monthly sampling data, abundance-based rarefaction, and sample-based rarefaction were used to determine the completeness of the natural pond inventories (see ESM2 in the Supplementary Material). Given that differences in sampling efforts were not observed among the 5 monthly-sampled ponds (ESM2 in the Supplementary Material), raw data were used to perform all statistical analyses. To analyze the differences in the macroinvertebrate assemblage compositions among the sampled ponds and months, we constructed a matrix of species (columns) and ponds (rows) including all the sampling months by averaging the number of individuals of each species per each pond and month. Then, triangular matrices were obtained using the Bray-Curtis index of similarity (the inverse of the Bray-Curtis index of dissimilarity). Assemblage similarity differences between ponds and months, as well as the possible interaction between these two factors (ponds × months), were analyzed using a two-way crossed permutational multivariate analysis of variance (PERMANOVA; . This non-parametric procedure allowed us to analyze differences between ponds and months based on the Bray-Curtis similarity distance. Statistically significant differences were tested using permutations of group membership; 9,999 permutations were performed using the software PRIMER v.6 (Anderson et al., 2008).
Differences in macroinvertebrate β-diversity patterns were investigated across the sampled months and ponds using presence-absence matrices and pooling adult and larval stages for each independent taxon. Then, triangular matrices were calculated using the Sørensen index of dissimilarity on the incidence data. The β-diversity partitioning approach proposed by Baselga (2010) was used to calculate the β-diversity patterns. Because the temporary ponds were drying out during the summer months, the per-month matrices varied in matrix size. Therefore, to determine β-diversity partitioning between ponds and to obtain comparable values of β-diversity for different matrix sizes, a multiple-site dissimilarity procedure was used (Baselga, 2012). To calculate the multiple-site β-diversity partitioning, the minimum number of ponds sampled in a month (three) was used to resample the total number of ponds sampled per month (1,000 random samples). To determine multiple-site β-diversity partitioning, the command "beta.sample" (betapart package, Baselga and Orme, 2012) was used, implemented in R software 2.14.2. Multiple-site β-diversity (β SOR ) was partitioned into two additive components that accounted for dissimilarity due to species replacement (β SIM ) and dissimilarity due to nestedness (β SNE ), respectively, in accordance with the formula β SOR = β SIM + β SNE (Baselga, 2010).

Environmental and Land-Use Effects on Macroinvertebrate Assemblages
Using monthly data, an environmental matrix per pond and month was constructed. With this matrix, triangular matrices of similarity were calculated using the Euclidean distance between each pair of ponds for each environmental variable. In order to analyze the environmental variables that influenced the monthly pattern of β-diversity, the multiple-site β-diversity partitioning, β SIM and β SNE , were regressed against the average Euclidean distances of each environmental variables per month. A forward stepwise regression was performed following Blanchet et al. (2008), using the β SIM and β SNE per month as response-variables and the average Euclidean distances of the environmental variables per month (N = 10) as predictor variables. This complete procedure is considered to be an effective procedure for controlling Type I error (Peres-Neto and Legendre, 2010). These analyses were performed using the Statistica V.8 software.
Regarding the extensive sampling, environmental data and the land-use variables were used to construct a matrix including the 12 natural ponds and seven artificial tanks as cases. This matrix was used first to analyze the explanatory variables that determined the differences in macroinvertebrate assemblages of natural and artificial ponds. Secondly, a matrix was constructed that excluded the artificial tanks to analyze the explanatory variables that determined the differences in macroinvertebrate assemblages only among the natural ponds. The studied environmental variables were the concentration of TP, i-P, and DIN, electrical conductivity at 20 • C, dissolved oxygen concentration, pH, and chlorophyll-a concentration in the water column. The studied land-use variables were the percentage area of native vegetation, seminatural pasturelands, and intensive pasturelands. All variables, with exception of pH and the land-use variables, were log (x+1) transformed to satisfy normality assumptions. Triangular matrices were constructed using Euclidean distances of the environmental variables to determine if the explanatory variables (environmental and landuse characteristics) of ponds and tanks differed. Then, an ANOSIM analysis was performed using the aquatic systems (pond or tank) as factor, and 9,999 permutations for the significance level. Moreover, each individual variable was analyzed using a Student t-test, and alternatively a Mann-Whitney test when the normality assumption was not satisfied, to determine whether it differed significantly between natural ponds and artificial tanks. Moreover, in order to detect if the macroinvertebrate assemblages differed between ponds and tanks, a matrix was constructed including the average number of individuals of each species (columns) per pond (rows). Then, a triangular similarity matrix was built using the Bray-Curtis index on the abundance matrix. Subsequently, an ANOSIM test was performed using pond and tank as grouping factor, and 9,999 permutations to assess the significance level. Finally, an exploratory SIMPER analysis was performed to detect the species mainly contributing to the differences between ponds and tanks.
Furthermore, to investigate the impact of the explanatory variables and the spatial descriptors driving the macroinvertebrate assemblages across the extensive sampling, the relationship between the biological Bray-Curtis similarity matrix, the aforementioned matrix of transformed explanatory variables and the spatial descriptors (PCNMs) was analyzed. To do that, a resemblance matrix was also constructed based on Bray-Curtis similarity that excluded the tanks to analyze the variables driving the differences among the assemblages of ponds. We first visualized the relationship between all biological data and the explanatory variables as principal component ordinations using a distance-based redundancy analysis (dbRDA) (McArdle and , which constrains to linear combinations of the predictor variables (Anderson et al., 2008). A Student t-test, or alternatively a Mann-Whitney test, was performed, using the scores of the axes 1 and 2 of the dbRDA (dbRDA1 and dbRDA2) to analyze if they discriminated natural ponds and artificial tanks. Secondly, a distance-based linear model (DistLM) was performed for data partitioning (analogous to linear multiple regression) using forward stepwise and Adjusted R 2 (AdjR 2 ) criteria for variable selection. This procedure was performed separately for the aforementioned explanatory variables and for the spatial descriptors, retaining the significant variables that contributed to the dissimilarities of the macroinvertebrate assemblages. This DistLM model only included predictor variables that improved the explained sum of squares that would be expected by adding some random variable, taking into account the number of variables in the model (Anderson et al., 2008). We thus estimated the sequential partial increase in explained variability using the AdjR 2 . The DistLM analysis was performed, first to analyze the differences between the tanks and natural ponds, and secondly, among the natural ponds alone. All these multivariate analyses were performed using the Primer V.6. software, the DistLM and dbRDA analyses were performed using the add-on package PERMANOVA+ (Anderson et al., 2008), while t-test and Mann-Whitney' test were performed in R software 2.14.2. In order to explore the partial contribution of the significant explanatory variables and the significant spatial descriptors (PCNMs) retained in the DistLM analyses on the macroinvertebrate assemblages, multiple regression on distance matrices (MRM) was performed, an extension of Mantel test (Legendre et al., 1994). We considered Spearman correlations (r s ) and performed a forward-selection procedure to identify the significant variables (Legendre et al., 1994). The significance of the MRM models was assessed using 1,000 permutations using the "MRM" command ("ecodist" package, Goslee and Urban, 2007) in R software 2.14.2. Two successive models were constructed: (i) the environmental model, only using the significant explanatory variables to measure environmental influence on the macroinvertebrate assemblages, and (ii) the environmental + spatial model, using the significant PCNM spatial descriptors in addition to the explanatory variables to obtain partial effects. This procedure was performed for the considered tanks and natural ponds (N = 19), and among the natural ponds alone (N = 12).

General Diversity Patterns
A total of 4,247 individuals (3,295 individuals in the monthly study and 952 in the extensive study) were collected, belonging to 17 species; Basommatophora (1), Ephemeroptera (1), Odonata (3), Heteroptera (2), Coleoptera (9), and Trichoptera (1) (see ESM3 in the Supplementary Material). In the monthly sampling, temporary ponds had the highest species richness (temporary ponds = 13 species, permanent ponds = 9 species). In contrast, permanent ponds had the highest proportion of native individuals. Overall, temporary ponds had the highest occurrence of endemic species (represented by Hydroporus guernei Régimbart, 1891 -Dytiscidae and Limnephilus atlanticus Nybom, 1948 -Limnephilidae), and exhibited the highest proportion of endemic individuals. Temporary and permanent ponds had a similar small proportion of exotic individuals (1 and 0.8%, respectively, Figure 2A). The seasonal pattern of macroinvertebrates revealed high abundance of the Fam. Dytiscidae, mainly consisting of the endemic H. guernei (61.3%), with both adults and larvae detected across three different periods (November, February and May, see Figure 3, and ESM5 in the Supplementary Material). Moreover, the Families Libellulidae and Corixidae were also abundant when most temporary ponds were close to the desiccation (April and May, Figure 3).
Notably, natural ponds harbored the highest proportion of individuals of native (84%) and endemic species (13%), in comparison with the artificial tanks (79% and 8%, respectively). Exotic species were more frequent in tanks (13%) in comparison with the natural ponds (3%, Figure 2B). The two types of systems shared three exotic species, i.e., Helisoma trivolvis (Say, 1817), Helochares lividus (Forster, 1771), and Cercyon haemorrhoidalis (Fabricius, 1775). Additionally, individuals of the gastropods Planorbidae and the coleopterans Hydrophilidae were mostly detected in tanks, where only Culicidae and Baetidae individuals were found. In contrast, the nymphs of the odonates Coenagrionidae, Libellulidae and the endemic L. atlanticus were only detected in natural ponds (Figure 4).

Temporal Biodiversity Patterns
Significant effects of both pond and month factors were detected in the macroinvertebrate assemblages, as well as a significant interaction between these two factors (PERMANOVA, Pond Factor: Pseudo-F = 3.801, P-value < 0.001; Month Factor: Pseudo-F = 12.623, P-value < 0.001; Month × Pond Factor: Pseudo-F = 2.548 P-value < 0.001). While pond explained a higher variance (14.47%) than the month factor (6.56%), the high explanatory variance of the interaction between the two factors, month and pond (13.38%), indicated that the macroinvertebrate differences between ponds were dependent on the study month. Analyzing the monthly β-diversity partitioning, the contribution of species replacement to the β-diversity was higher than the contribution of the nestedness component across almost the whole seasonal pattern. However, in August, when most temporary ponds had dried out and permanent ponds had notably reduced their surface area, the value of the nestedness contribution to β-diversity increased, while the value of species replacement notably decreased (Figure 5).

Environmental and Land-Use Effects on Macroinvertebrate Assemblages
In the monthly sampling, the dissolved oxygen concentration, the TP concentrations, the electrical conductivity and the pH significantly explained the monthly pattern of species   FIGURE 5 | Average values of the components of the macroinvertebrate β-diversity, the species replacement (β SIM ) and the contribution of the nestedness (β SNE ) to the β-diversity after β-diversity partitioning. Bars indicate the standard deviation (SD); July and August do not show SD because only a pair-wise comparison between the three ponds that maintained water by that time was performed.
replacement, while only chlorophyll-a concentration significantly explained the monthly pattern of the contribution of nestedness to the observed β-diversity (Table 1).
In most of the natural ponds, more than 98% of the surrounding area comprised natural vegetation (FUN, NGR, PHN, and SER). However, GIN pond was located in a disturbed area only surrounded by intensive pasturelands ( Table 2). Tanks were located in pasturelands, most of them contiguous to natural ponds located in areas of native vegetation (Figure 1). The surrounding area of TPI, TA4, and TA5 included ∼50% of native vegetation, while some tanks (e.g., TA3, TA6, and TA7) were located in areas dominated by intensive pasturelands (Table 2). Hence, significant differences in the environmental explanatory variables were observed between these two groups of ponds and tanks (ANOSIM, Global R = 0.979, P-value < 0.001). In contrast, statistically significant but weak differences were observed between the macroinvertebrate assemblages of the natural ponds and artificial tanks (ANOSIM, Global R = 0.186, P-value < 0.05). The SIMPER analysis revealed that the only taxa with a contribution >10% to this dissimilarity were the larvae of Corixa affinis and the Oligochaeta, mainly occurring in the artificial tanks (see ESM5 in the Supplementary Material).
In the dbRDA ordination, we observed that pH was the main difference between the assemblage compositions of ponds and tanks (Figure 6). The scores of axes 1 and 2 of the dbRDA significantly differed between natural ponds and tanks (dbRDA1 t = 2.879, P-value < 0.01; dbRDA2 W = 38, Pvalue < 0.0001). Hence, the DistLM analysis revealed that pH was the only significant variable explaining the differences among the macroinvertebrate assemblages of ponds and tanks (5% of explained variability, Table 3). However, another four variables were also retained in the DistLM analysis, which explained 12% of the total variability (Table 3). Only using natural ponds, the dbRDA1 axis could also discriminate them into two groups: the ponds with relatively high concentrations of DIN and dissolved oxygen, and mostly surrounded by native vegetation (group 1: PHN, SER, FUN, NGR, NEG, VF1, VF3, and VF5 ponds), and those with high concentrations of TP and chlorophyll-a (group 2: GIN, PAT, BOI, and ARE ponds). When the DistLM analysis was performed only using these natural ponds, the concentration of TP in the water column was the single significant explanatory DO is dissolved oxygen concentration, EC is electrical conductivity, TP is total phosphorous concentration, Chl. is chlorophyll-a concentration, DIN is dissolved inorganic nitrogen concentration and i-P indicates dissolved inorganic phosphate concentration. The percentage areas of native vegetation (Nat), seminatural pasturelands (Sem), and intensive pasturelands (Int) are also indicated as land-use variables. A Student t-test was performed to analyze if the environmental variables were significantly different between natural ponds and artificial tanks, and a non-parametric Mann-Whitney test for Sem and Int, which did not satisfy normality assumption. Diff indicates the statistic values of these two parametric (t) and non-parametric (W) tests, respectively; the P-values are also indicated as * < 0.05,** < 0.01,*** < 0.001; all the non-significant differences were P-values > 0.06 (indicated as n.s.).
variable (10% of explained variability, Table 3). However, another three variables were retained in the DistLM analysis, which explained 21% of the total variability (Table 3).

Environmental Variables vs. Spatial Descriptors
In the extensive sampling, seven out of the 11 spatial descriptors were retained in the DistLM analysis. Only the PCNM7 had a significant influence on the macroinvertebrate assemblages (Table 3). However, when the partial effect of the PCNM7 was analyzed taking into account the pH, the model revealed that only the pH was significant in explaining the dissimilarity of the macroinvertebrate assemblages (Table 4). When the DistLM analysis was repeated with only the natural ponds, two significant spatial descriptors (PCNM2 and PCNM5) were detected out of the five PCNMs retained. The partial effect of these two spatial descriptors, taking into account the concentration of TP in the water column, revealed the important effect of these two PCNMs on the dissimilarities in the macroinvertebrate assemblages, as they explained 34% of the total variability in the environmental + spatial model (Table 4).

DISCUSSION
In order to understand the extent of the anthropogenic influence on the unique freshwater ecosystems of Macaronesia, we have first investigated the ecological dynamic of these singular ecosystems and their macroinvertebrate assemblages. In the Azorean island of Terceira, a consistent seasonal pattern of macroinvertebrate species replacement has been detected across a hydroperiod gradient that disappeared when most temporary ponds were desiccated in summer. This desiccation period was thus observed to be an important factor for the community structure of the Azorean macroinvertebrates, as occurs in other sensitive regions of temporary pond networks, such as in Mediterranean-climate areas (Florencio et al., 2011). The role of artificial tanks as possible reservoirs for macroinvertebrates might be especially relevant as refuges during aestivation for those organisms inhabiting temporary ponds, which need to cope with pond desiccation in summer (see Deacon et al., 2019;Samways et al., 2020). Our results suggest that the studied artificial tanks maintain good water quality and are suitable for macroinvertebrate species, thus increasing the heterogeneity of habitats as permanent waterbodies in the peak of the hydroperiod gradient. This is even more important if we consider that fishless permanent waterbodies are scarce in the Azorean islands (Florencio and Lamelas-López, 2016;Raposeiro et al., 2017).

Characteristics of The Azorean Natural Ponds
The natural ponds in Terceira Island were characterized by well-oxygenated and slightly acidic waters, with low values of conductivity, which indicates that the ponds are primarily flooded with rainwater (Florencio and Lamelas-López, 2016). Differences in macroinvertebrate composition among ponds were observed to be higher than seasonal differences (monthly), which can be associated with the limited fluctuations in temperatures throughout the year under a temperate oceanic climate. This result contrasts with the typical seasonal variations of European ponds, where changes in the environmental characteristics are reflected in seasonal variations of assemblage compositions (e.g., Jeffries, 2003). Seasonal differences among ponds in the macroinvertebrate assemblages reflected the desiccation of the temporary ponds and the reduction of the water level in the permanent ponds. Interestingly, the larvae and adults of the endemic beetle H. guernei peaked three times during the hydrological cycle of temporary ponds, i.e., November (after pond inundation), February, and May (close to pond desiccation). This result seems to indicate at least three successive periods of reproduction for this species, which is important for conservation regarding: (1) the Endangered category of this species in the IUCN Red List of Threatened Species , and (2) the protected character of these high-altitude natural ponds in the Azores. We acknowledge that a strong signal of seasonal variation in assemblage composition was detected though the number of ponds was low. It is likely, then, that the sampled ponds provided a good representation of the study area. Seasonal β-diversity patterns revealed that the differences in assemblage compositions were mainly associated with species replacement, and that pond desiccation favored the contribution of nestedness to the observed β-diversity, a pattern that has also been found in other temporary pond networks (Florencio et al., 2011(Florencio et al., , 2016a. The seasonal variation of the dissolved oxygen and TP concentrations, electrical conductivity and pH determined the seasonal changes in species replacement, probably associated with the variation in the timing and amount of rainfall. However, the proximity to the dry-period increased the contribution of nestedness to the β-diversity and also reduced species replacement. Close to pond desiccation, high concentrations of TP and chlorophyll-a were reached in concordance with a reduction in the inundated surface area, which also seemed to influence the pH and electrical conductivity of the water. Pond desiccation is an asynchronous process within a temporary pond network, as it depends on the hydroperiod of each pond. In Mediterranean ponds, water depth and electrical conductivity have also been detected as key variables explaining changes in assemblage composition and β-diversity associated with species replacement (Florencio et al., 2014). In our results, most temporary ponds were desiccated in August, when the values of nestedness contribution to β-diversity were maximal, mainly influenced by a high chlorophyll-a concentration. Fish were detected in the only two permanent ponds that held water during August during the study period (Florencio and Lamelas-López, 2016), and these were the only ponds that were able to act as refuge for macroinvertebrates in summer.

Artificial Tanks as Reservoirs for Biodiversity
Fishless artificial tanks associated with agricultural systems acquire an important role in preserving the macroinvertebrate assemblages of Azorean ponds. This possibility is supported by the similar composition that has been detected in the macroinvertebrate assemblages of ponds and tanks, with only small variations explained by the notably higher pH of artificial tanks. This could be explained by the longer persistence of water, smaller rain catchment (i.e., smaller basin) and the construction materials of the tanks (basic materials; Hassall, 2014). Hence, the high pH values in tanks could have favored the gastropods of the Fam. Planorbidae (e.g., Şahin and Zeybek, 2016). The conservation of macroinvertebrates in artificial tanks was also supported by their low concentration of nutrients and chlorophyll-a, reflecting a low probability of eutrophication (see Campbell et al., 2009). Consequently, tanks can be considered potential suitable habitats for at least 10 species of macroinvertebrates (see ESM3 in the Supplementary Material). Notably, some native species were only captured in tanks, such as the beetles Hygrotus confluens (Fabricius, 1787) or the single species of the Ephemeroptera order known in Azores, Cloeon 3 | Partial effects of the explanatory variables (Environmental) and the spatial descriptors based on PCNMs (Spatial) using distance-based linear models (DistLM), performing a forward stepwise process based on the Adjusted R 2 (AdjR 2 ).

Environmental
AdjR 2  Environmental and Spatial DistLM analyses have been performed separately. The Pseudo-F and the P-values are also indicated. The AdjR 2 indicates the explained variability that added the inclusion of each explanatory variable in a sequential order. These analyses were performed using the biological resemblance matrix (Bray-Curtis similarity index) and the environmental and land-use variables to analyze the explanatory variables that mainly contributed to the differences in the assemblage compositions between ponds and tanks (Ponds vs. tanks), but also among the ponds alone (Only ponds). DO is dissolved oxygen concentration, TP is total phosphorous concentration, Chl is chlorophyll-a concentration, DIN is dissolved inorganic nitrogen concentration, i-P is dissolved inorganic phosphate concentration. For the land-use variables, Nat is the percentage area of native vegetation and Sem is the percentage area of seminatural pasturelands. We only used those retained variables that were significant in the DistLM analyses (see Table 3), including only the explanatory variables (Environmental) in the MRM model, and including the spatial descriptors in this model (Environmental + Spatial). TP is total phosphorous concentration in the water column, r 2 is the model explanatory capacity (ranged 0-1), rs indicates the coefficients of Spearman correlations.
dipterum (Linnaeus, 1761). This is possibly associated with the preference of these species for higher water temperature, as has been registered in tanks (see Stauder, 1991;McKee and Atkinson, 2000). In Azores, Ephemeroptera mainly occurs in human-made systems (Brinck and Scherer, 1957), contributing to increase the total biodiversity. We also found adults and larvae of the endemic species Hydroporus guernei and larvae of the Limnephilus atlanticus  in tanks, although these species were mainly recorded in natural temporary ponds.
Other species can benefit from the long hydroperiod of artificial tanks, such as the dragonfly Anax imperator Leach, 1815 (Fam. Aeshnidae), which usually prefers permanent waters to develop its instars and complete metamorphosis (see Corbet, 1957). However, natural ponds harbored the highest proportion of individuals belonging to native non-endemic and endemic species, most of which (8 out of 13 species) were also detected in artificial tanks. Additionally, although exotic species were more frequently detected in tanks (13% of all the collected individuals in tanks belong to exotic species), these exotic species do not seem to have so far successfully colonized natural ponds, which are only harboring 3% of individuals belonging to exotic species. This is probably due to the pristine conditions of habitats surrounding the natural ponds. However, the artificial tanks had a higher proportion of exotic species than natural ponds despite the short distance between them (180-2,130 meters). It is possible that, due to their location in the pristine native forests, the natural ponds could be able to resist the invasion of exotic species (i) as well-preserved communities (biotic resistance to invasion, sensu Elton, 1958), (ii) because high-altitude native forests could act as physical barriers to the colonization of exotic species (Florencio et al., 2016b), (iii) and possibly also due to the temporary character of some natural ponds, which harbor a singular fauna adapted to pond desiccation (Williams, 2006). Hence, when comparing two distinct types of systems (natural vs. artificial), the macroinvertebrate dissimilarities between them seemed to be lower than those found between ponds located in well-preserved areas and those with signs of anthropogenic disturbance. The former group exhibited slightly higher concentrations of DIN and dissolved oxygen, which could be related to a different exposure to the wind, given that these ponds are generally located at higher altitude in the native forest. The latter group of ponds were associated with human activities, located at lower altitude where agricultural activities and cattle ranching are more common, which was reflected in the slightly higher concentrations of TP and chlorophyll-a. Nonetheless, both natural ponds and tanks exhibited a relatively low concentration of nutrients (DIN and TP) that fell within the range of other wetlands and ponds located in protected areas of comparable rainfall regimes (Plenzler and Michaels, 2015). Natural ponds displayed a high degree of endemicity that was mainly associated with temporary ponds, and a low species richness of macroinvertebrates. These two phenomena of high endemicity and "disharmony" (impoverished biota in comparison with the adjacent mainland, even lacking entire taxonomical groups) are typical of the insular biota, and have been described for the freshwater invertebrates of the Azores (Raposeiro et al., 2012;Florencio and Lamelas-López, 2016). The macroinvertebrate assemblages of these natural ponds were spatially structured, but this spatial pattern was not detected when tanks and ponds were considered together. This result suggests that the aquatic organisms benefit from tanks as refuges (e.g., to cope with the desiccation of temporary ponds), though they are not suitable habitats to establish successful populations (but see Svensson, 1977).

CONCLUSIONS
Despite storing water for livestock farming, it has been shown that artificial tanks can play a beneficial role in preserving macroinvertebrate assemblages. These freshwater systems of anthropogenic origin can harbor a wide range of taxa (see Hassall, 2014;Hill et al., 2015), enhancing regional diversity without replacing the essential conservation value of natural ponds (Reyne et al., 2020). In the present study, the artificial tanks in Terceira Island provided permanent waters that were able to harbor endemic and native species, increasing the heterogeneity of habitats along the hydroperiod gradient. The good water quality of the study tanks has also been shown, as they exhibited low concentrations of nutrients and chlorophyll-a, and low conductivity. These tanks increased the regional biodiversity despite their concrete walls without aquatic vegetation, both considered possible disruptive elements for macroinvertebrate diversity in urban/rural ponds (Oertli and Parris, 2019). However, we suggest that these concrete walls could have prevented grazing cattle from entering the water, which can significantly affect water quality (Campbell et al., 2009). Nevertheless, the occurrence of submerged and littoral vegetation might provide food, refuge and structure (e.g. for reproduction) for some macroinvertebrates (Fuentes-Rodríguez et al., 2013) and amphibians (Swartz and Miller, 2019). In our study, the lack of aquatic vegetation in the tanks could have prevented the proper establishment of odonate species (see Foote and Hornung, 2005;Hykel et al., 2020), as well as the occurrence of other aquatic species (e.g., aquatic beetle and bug species, see Deacon et al., 2018). Additionally, it is not known if tanks located at greater distances from the protected areas would also act as refuge habitats of macroinvertebrate assemblages.
Our study suggests some recommendations to preserve the macroinvertebrate diversity of the Azorean ponds. (1) Maintain optimal environmental conditions in tanks like those found in natural and non-disturbed ponds (except for the pH of the water, due to the nature and physical structure of the tanks). (2) Maintain a high level of conservation of the pristine areas that contain the natural ponds, which allow most species to develop their life cycles and sustain high biodiversity, and protect these ponds against the arrival of potential invaders (see Samways et al., 2020). (3) The identification of anthropogenic pressures on the most disturbed ponds is essential to establish management strategies for conservation, mainly in the agricultural areas.
(4) Further studies should evaluate possible strategies for the management of the artificial tanks to establish suitable conservation guidelines at local and landscape scales, that integrates both artificial tanks and natural ponds. Some examples of these strategies may be the inclusion of aquatic vegetation in the artificial tanks, as well as building facilities in their concrete walls (e.g., ramps) to favor the entrance/exit of some aquatic organisms. (5) The design of corridors from highaltitude native forests to the coast (see Aparício et al., 2018), including some of the study ponds and close-located tanks, could promote the conservation of the aquatic macroinvertebrates of the Terceira Island. We also highlight the importance of conducting seasonal studies that cover complete inundationdesiccation periods of ponds to properly represent the entire macroinvertebrate community, and to establish appropriate guidelines for conservation regarding the whole community. The singularity and vulnerability of the Azorean ponds must be highlighted, mainly regarding the scarce number of temporary ponds in Terceira Island, and we recommend including them in conservation programs, such as the Ramsar Convention.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
MF and LL-L conceived the idea. LL-L performed the field sampling with the help of MF. VG and LS carried out some chemical determinations in the laboratory. LL-L and MF wrote the manuscript with significant contributions by PB, and all the authors revised the different versions of the manuscript.

FUNDING
We are grateful to the Fundación Barrié for the academic fellowship of LL-L. MF's grant has been funded by the Conselho Nacional de Desenvolvimento Científico e Tecnológico -CNPq (401045/2014-5), program Ciência sem Fronteiras. MF's current contract was supported by the Universidad Autónoma de Madrid. MF is also grateful to the project ClimaRiskinPond (PID2019-104580GA-I00/AEI/10.13039/501100011033) funded by the Spanish Ministerio de Ciencia e Innovación, and PB was supported by the project FCT-UIDB/00329/2020.

ACKNOWLEDGMENTS
We are grateful to Rui Elias and Rosalina Gabriel who provided the technical equipment for pH, dissolved oxygen concentration, and electrical conductivity measurements in the water column. We also appreciate the contributions of Adolfo Cordero for his support and comments, and to Herlander Lima for performing some chemical analyses in the laboratory. We are also grateful to Luis Ansias, Jose Sarangollo, Herlander Lima, Raquel Daza, Lucía Lamelas, David Rodilla, Fernando Pereira, and Paco Gutiérrez for their assistance in fieldwork.

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fevo.

2020.605176/full#supplementary-material
Supplementary Material 1 | Photographs of some of the study natural ponds (Supplementary Figure 1), some of the artificial tanks ( Supplementary Figure 2), and a list of the local names of the study natural ponds (Supplementary Table 1).
Supplementary Material 2 | Completeness of the pond inventories through abundance-and sample-based rarefaction curves to assess if the pond inventories are comparable in sampling efforts.
Supplementary Material 3 | Check-list including the relative abundance of the taxa captured during the extensive sampling. The relative abundance was calculated as the average number of individuals captured, considering the number of sampling units per pond and tank.
Supplementary Material 4 | Check-list of the captured taxa, including the number of individuals per pond and month, during the monthly sampling ("monthly sampling" spreadsheet), and the environmental variables of these ponds per sampling month ("variables monthly sampling" spreadsheet). A indicates Adult, L Larva, -Not captured. DO is dissolved oxygen concentration, EC is electrical conductivity, TP is total phosphorous concentration, and Chl is chlorophyll-a concentration.
Supplementary Material 5 | Independent percentage of contribution of each species to the macroinvertebrate dissimilarities detected between natural ponds and artificial tanks after a SIMPER analysis.